Repressive KRAB domain-containing zinc-finger proteins (KRAB-ZFPs) are abundant in mammalian genomes and contribute both to the silencing of transposable elements (TEs) and to the regulation of developmental stage- and cell type-specific gene expression. Here we describe studies of zinc finger protein 92 (Zfp92), an X-linked KRAB-ZFP that is highly expressed in pancreatic islets of adult mice, by analyzing global Zfp92 knockout (KO) mice. Physiological, transcriptomic and genome-wide chromatin binding studies indicate that the principal function of ZFP92 in mice is to bind to and suppress the activity of B1/Alu type of SINE elements and modulate the activity of surrounding genomic entities. Deletion of Zfp92 leads to changes in expression of select LINE and LTR retroelements and genes located in the vicinity of ZFP92-bound chromatin. The absence of Zfp92 leads to altered expression of specific genes in islets, adipose and muscle that result in modest sex-specific alterations in blood glucose homeostasis, body mass and fat accumulation. In islets, Zfp92 influences blood glucose concentration in postnatal mice via transcriptional effects on Mafb, whereas in adipose and muscle, it regulates Acacb, a rate-limiting enzyme in fatty acid metabolism. In the absence of Zfp92, a novel TE-Capn11 fusion transcript is overexpressed in islets and several other tissues due to de-repression of an IAPez TE adjacent to ZFP92-bound SINE elements in intron 3 of the Capn11 gene. Together, these studies show that ZFP92 functions both to repress specific TEs and to regulate the transcription of specific genes in discrete tissues.

We use ATAC-seq to examine chromatin accessibility for four different tissues in Drosophila melanogaster: adult female brain, ovaries, and both wing and eye-antennal imaginal discs from males. Each tissue is assayed in eight different inbred strain genetic backgrounds, seven associated with a reference quality genome assembly. We develop a method for the quantile normalization of ATAC-seq fragments and test for differences in coverage among genotypes, tissues, and their interaction at 44099 peaks throughout the euchromatic genome. For the strains with reference quality genome assemblies, we correct ATAC-seq profiles for read mis-mapping due to nearby polymorphic structural variants (SVs). Comparing coverage among genotypes without accounting for SVs results in a highly elevated rate (55%) of identifying false positive differences in chromatin state between genotypes. After SV correction, we identify 1050, 30383, and 4508 regions whose peak heights are polymorphic among genotypes, among tissues, or exhibit genotype-by-tissue interactions, respectively. Finally, we identify 3988 candidate causative variants that explain at least 80% of the variance in chromatin state at nearby ATAC-seq peaks.

MicroRNAs (miRNAs) are a well-characterized class of small RNAs (sRNAs) that regulate gene expression post-transcriptionally. miRNAs function within a complex milieu of other sRNAs of similar size and abundance, with the best characterized being tRNA fragments or tRFs. The mechanism by which the RNA-induced silencing complex (RISC) selects for specific sRNAs over others is not entirely understood in human cells. Several highly expressed tRNA trailers (tRF-1s) are strikingly similar to microRNAs in length but are generally excluded from the microRNA effector pathway. This exclusion provides a paradigm for identifying mechanisms of RISC selectivity. Here, we show that 5' to 3' exoribonuclease XRN2 contributes to human RISC selectivity. Although highly abundant, tRF-1s are highly unstable and degraded by XRN2 which blocks tRF-1 accumulation in RISC. We also find that XRN mediated degradation of tRF-1s and subsequent exclusion from RISC is conserved in plants. Our findings reveal a conserved mechanism that prevents aberrant entry of a class of highly produced sRNAs into Ago2.

The current model for Okazaki fragment maturation in bacteria invokes RNA cleavage by RNase H, followed by strand displacement synthesis and 5' RNA flap removal by DNA polymerase I (Pol I). RNA removal by Pol I is thought to occur through the 5'-3' flap endo/exonuclease (FEN) domain, located in the N-terminus of the protein. In addition to Pol I, many bacteria encode a second, Pol I-independent FEN. The contribution of Pol I and Pol I-independent FENs to DNA replication and genome stability remains unclear. In this work we purified Bacillus subtilis Pol I and FEN, then assayed these proteins on a variety of RNA-DNA hybrid and DNA-only substrates. We found that FEN is far more active than Pol I on nicked double-flap, 5' single flap, and nicked RNA-DNA hybrid substrates. We show that the 5' nuclease activity of B. subtilis Pol I is feeble, even during DNA synthesis when a 5' flapped substrate is formed modeling an Okazaki fragment intermediate. Examination of Pol I and FEN on DNA-only substrates shows that FEN is more active than Pol I on most substrates tested. Further experiments show that ΔpolA phenotypes are completely rescued by expressing the C-terminal polymerase domain while expression of the N-terminal 5' nuclease domain fails to complement ΔpolA. Cells lacking FEN (ΔfenA) show a phenotype in conjunction with an RNase HIII defect, providing genetic evidence for the involvement of FEN in Okazaki fragment processing. With these results, we propose a model where cells remove RNA primers using FEN while upstream Okazaki fragments are extended through synthesis by Pol I. Our model resembles Okazaki fragment processing in eukaryotes, where Pol δ catalyzes strand displacement synthesis followed by 5' flap cleavage using FEN-1. Together our work highlights the conservation of ordered steps for Okazaki fragment processing in cells ranging from bacteria to human.

Gene-body methylation (gbM) refers to sparse CG methylation of coding regions, which is especially prominent in evolutionarily conserved house-keeping genes. It is found in both plants and animals, but is directly and stably (epigenetically) inherited over multiple generations in the former. Studies in Arabidopsis thaliana have demonstrated that plants originating from different parts of the world exhibit genome-wide differences in gbM, which could reflect direct selection on gbM, but which could also reflect an epigenetic memory of ancestral genetic and/or environmental factors. Here we look for evidence of such factors in F2 plants resulting from a cross between a southern Swedish line with low gbM and a northern Swedish line with high gbM, grown at two different temperatures. Using bisulfite-sequencing data with nucleotide-level resolution on hundreds of individuals, we confirm that CG sites are either methylated (nearly 100% methylation across sampled cells) or unmethylated (approximately 0% methylation across sampled cells), and show that the higher level of gbM in the northern line is due to more sites being methylated. Furthermore, methylation variants almost always show Mendelian segregation, consistent with their being directly and stably inherited through meiosis. To explore how the differences between the parental lines could have arisen, we focused on somatic deviations from the inherited state, distinguishing between gains (relative to the inherited 0% methylation) and losses (relative to the inherited 100% methylation) at each site in the F2 generation. We demonstrate that deviations predominantly affect sites that differ between the parental lines, consistent with these sites being more mutable. Gains and losses behave very differently in terms of the genomic distribution, and are influenced by the local chromatin state. We find clear evidence for different trans-acting genetic polymorphism affecting gains and losses, with those affecting gains showing strong environmental interactions (G×E). Direct effects of the environment were minimal. In conclusion, we show that genetic and environmental factors can change gbM at a cellular level, and hypothesize that these factors can also lead to transgenerational differences between individuals via the inclusion of such changes in the zygote. If true, this could explain genographic pattern of gbM with selection, and would cast doubt on estimates of epimutation rates from inbred lines in constant environments.

Methyl salicylate is an important inter- and intra-plant signaling molecule, but is deemed undesirable by humans when it accumulates to high levels in ripe fruits. Balancing the tradeoff between consumer satisfaction and overall plant health is challenging as the mechanisms regulating volatile levels have not yet been fully elucidated. In this study, we investigated the accumulation of methyl salicylate in ripe fruits of tomatoes that belong to the red-fruited clade. We determine the genetic diversity and the interaction of four known loci controlling methyl salicylate levels in ripe fruits. In addition to Non-Smoky Glucosyl Transferase 1 (NSGT1), we uncovered extensive genome structural variation (SV) at the Methylesterase (MES) locus. This locus contains four tandemly duplicated Methylesterase genes and genome sequence investigations at the locus identified nine distinct haplotypes. Based on gene expression and results from biparental crosses, functional and non-functional haplotypes for MES were identified. The combination of the non-functional MES haplotype 2 and the non-functional NSGT1 haplotype IV or V in a GWAS panel showed high methyl salicylate levels in ripe fruits, particularly in accessions from Ecuador, demonstrating a strong interaction between these two loci and suggesting an ecological advantage. The genetic variation at the other two known loci, Salicylic Acid Methyl Transferase 1 (SAMT1) and tomato UDP Glycosyl Transferase 5 (SlUGT5), did not explain volatile variation in the red-fruited tomato germplasm, suggesting a minor role in methyl salicylate production in red-fruited tomato. Lastly, we found that most heirloom and modern tomato accessions carried a functional MES and a non-functional NSGT1 haplotype, ensuring acceptable levels of methyl salicylate in fruits. Yet, future selection of the functional NSGT1 allele could potentially improve flavor in the modern germplasm.

The prototype enzymes of the ubiquitous type IA topoisomerases (topos) family are Escherichia coli topo I (topA) and topo III (topB). Topo I shows preference for relaxation of negative supercoiling and topo III for decatenation. However, as they could act as backups for each other or even share functions, strains lacking both enzymes must be used to reveal the roles of type IA enzymes in genome maintenance. Recently, marker frequency analysis (MFA) of genomic DNA from topA topB null mutants revealed a major RNase HI-sensitive DNA peak bordered by Ter/Tus barriers, sites of replication fork fusion and termination in the chromosome terminus region (Ter). Here, flow cytometry for R-loop-dependent replication (RLDR), MFA, R-loop detection with S9.6 antibodies, and microscopy were used to further characterize the mechanism and consequences of over-replication in Ter. It is shown that the Ter peak is not due to the presence of a strong origin for RLDR in Ter region; instead RLDR, which is partly inhibited by the backtracking-resistant rpoB*35 mutation, appears to contribute indirectly to Ter over-replication. The data suggest that RLDR from multiple sites on the chromosome increases the number of replication forks trapped at Ter/Tus barriers which leads to RecA-dependent DNA amplification in Ter and to a chromosome segregation defect. Overproducing topo IV, the main cellular decatenase, does not inhibit RLDR or Ter over-replication but corrects the chromosome segregation defect. Furthermore, our data suggest that the inhibition of RLDR by topo I does not require its C-terminal-mediated interaction with RNA polymerase. Overall, our data reveal a pathway of genomic instability triggered by R-loops and its regulation by various topos activities at different steps.

Most non-spherical bacteria rely on the actin-like MreB cytoskeleton to control synthesis of a cell-shaping and primarily rod-like cell wall. Diverging from simple rod shape generally requires accessory cytoskeletal elements, which locally interfere with the MreB-guided cell wall synthesis. Conserved and widespread representatives of this accessory cytoskeleton are bactofilins that polymerize into static, non-polar bundles of filaments. Intriguingly, many species of the Actinobacteria and Rhizobiales manage to grow rod-like without MreB by tip extension, yet some of them still possess bactofilin genes, whose function in cell morphogenesis is unknown. An intricate representative of these tip-growing bacteria is Rhodomicrobium vannielii; a member of the hitherto genetically not tractable and poorly studied Hyphomicrobiaceae within the MreB-less Rhizobiales order. R. vannielii displays complex asymmetric cell shapes and differentiation patterns including filamentous hyphae to produce offspring and to build dendritic multicellular arrays. Here, we introduce techniques to genetically access R. vannielii, and we elucidate the role of bactofilins in its sophisticated morphogenesis. By targeted mutagenesis and fluorescence microscopy, protein interaction studies and peptidoglycan incorporation analysis we show that the R. vannielii bactofilins are associated with the hyphal growth zones and that one of them is essential to form proper hyphae. Another paralog is suggested to represent a novel hybrid and co-polymerizing bactofilin. Notably, we present R. vannielii as a powerful new model to understand prokaryotic cell development and control of multipolar cell growth in the absence of the conserved cytoskeletal element, MreB.

Females with polycystic ovary syndrome (PCOS), the most common endocrine disorder in women, have an increased risk of developing cardiometabolic disorders such as insulin resistance, obesity, and type 2 diabetes (T2D). While only diagnosable in females, males with a family history of PCOS can also exhibit a poor cardiometabolic profile. Therefore, we aimed to elucidate the role of sex in the cardiometabolic comorbidities observed in PCOS by conducting bidirectional genetic risk score analyses in both sexes. We first conducted a phenome-wide association study (PheWAS) using PCOS polygenic risk scores (PCOSPRS) to identify potential pleiotropic effects of PCOSPRS across 1,380 medical conditions recorded in the Vanderbilt University Medical Center electronic health record (EHR) database, in females and males. After adjusting for age and genetic ancestry, we found that European (EUR)-ancestry males with higher PCOSPRS were significantly more likely to develop hypertensive diseases than females at the same level of genetic risk. We performed the same analysis in an African (AFR)-ancestry population, but observed no significant associations, likely due to poor trans-ancestry performance of the PRS. Based on observed significant associations in the EUR-ancestry population, we then tested whether the PRS for comorbid conditions (e.g., T2D, body mass index (BMI), hypertension, etc.) also increased the odds of a PCOS diagnosis. Only BMIPRS and T2DPRS were significantly associated with a PCOS diagnosis in EUR-ancestry females. We then further adjusted the T2DPRS for measured BMI and BMIresidual (regressed on the BMIPRS and enriched for the environmental contribution to BMI). Results demonstrated that genetically regulated BMI primarily accounted for the relationship between T2DPRS and PCOS. Overall, our findings show that the genetic architecture of PCOS has distinct sex differences in associations with genetically correlated cardiometabolic traits. It is possible that the cardiometabolic comorbidities observed in PCOS are primarily explained by their shared genetic risk factors, which can be further influenced by biological variables including sex and BMI.

Changes in gene regulation represent an important path to generate developmental differences affecting anatomical traits. Interspecific divergence in gene expression often results from changes in transcription-stimulating enhancer elements. While gene repression is crucial for precise spatiotemporal expression patterns, the relative contribution of repressive transcriptional silencers to regulatory evolution remains to be addressed. Here, we show that the Drosophila pigmentation gene ebony has mainly evolved through changes in the spatial domains of silencers patterning its abdominal expression. By precisely editing the endogenous ebony locus of D. melanogaster, we demonstrate the requirement of two redundant abdominal enhancers and three silencers that repress the redundant enhancers in a patterned manner. We observe a role for changes in these silencers in every case of ebony evolution observed to date. Our findings suggest that negative regulation by silencers likely has an under-appreciated role in gene regulatory evolution.

How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.

Molecular chaperones play a key role in maintaining proteostasis and cellular health. The abundant, essential, cytosolic Hsp90 (Heat shock protein, 90 kDa) facilitates the folding and activation of hundreds of newly synthesized or misfolded client proteins in an ATP-dependent folding pathway. In a simplified model, Hsp70 first helps load client onto Hsp90, ATP binding results in conformational changes in Hsp90 that result in the closed complex, and then less defined events result in nucleotide hydrolysis, client release and return to the open state. Cochaperones bind and assist Hsp90 during this process. We previously identified a series of yeast Hsp90 mutants that appear to disrupt either the 'loading', 'closing' or 'reopening' events, and showed that the mutants had differing effects on activity of some clients. Here we used those mutants to dissect Hsp90 and cochaperone interactions. Overexpression or deletion of HCH1 had dramatically opposing effects on the growth of cells expressing different mutants, with a phenotypic shift coinciding with formation of the closed conformation. Hch1 appears to destabilize Hsp90-nucleotide interaction, hindering formation of the closed conformation, whereas Cpr6 counters the effects of Hch1 by stabilizing the closed conformation. Hch1 and the homologous Aha1 share some functions, but the role of Hch1 in inhibiting progression through the early stages of the folding cycle is unique. Sensitivity to the Hsp90 inhibitor NVP-AUY922 also correlates with the conformational cycle, with mutants defective in the loading phase being most sensitive and those defective in the reopening phase being most resistant to the drug. Overall, our results indicate that the timing of transition into and out of the closed conformation is tightly regulated by cochaperones. Further analysis will help elucidate additional steps required for progression through the Hsp90 folding cycle and may lead to new strategies for modulating Hsp90 function.

It remains unknown to what extent gene-gene interactions contribute to complex traits. Here, we introduce a new approach using predicted gene expression to perform exhaustive transcriptome-wide interaction studies (TWISs) for multiple traits across all pairs of genes expressed in several tissue types. Using imputed transcriptomes, we simultaneously reduce the computational challenge and improve interpretability and statistical power. We discover (in the UK Biobank) and replicate (in independent cohorts) several interaction associations, and find several hub genes with numerous interactions. We also demonstrate that TWIS can identify novel associated genes because genes with many or strong interactions have smaller single-locus model effect sizes. Finally, we develop a method to test gene set enrichment of TWIS associations (E-TWIS), finding numerous pathways and networks enriched in interaction associations. Epistasis is may be widespread, and our procedure represents a tractable framework for beginning to explore gene interactions and identify novel genomic targets.

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features-referred to as canonical variables (CVs)-within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Pbp1 (poly(A)-binding protein-binding protein 1) is a cytoplasmic stress granule marker that is capable of forming condensates that function in the negative regulation of TORC1 signaling under respiratory conditions. Polyglutamine expansions in its mammalian ortholog ataxin-2 lead to spinocerebellar dysfunction due to toxic protein aggregation. Here, we show that loss of Pbp1 in S. cerevisiae leads to decreased amounts of mRNAs and mitochondrial proteins which are targets of Puf3, a member of the PUF (Pumilio and FBF) family of RNA-binding proteins. We found that Pbp1 supports the translation of Puf3-target mRNAs in respiratory conditions, such as those involved in the assembly of cytochrome c oxidase and subunits of mitochondrial ribosomes. We further show that Pbp1 and Puf3 interact through their respective low complexity domains, which is required for Puf3-target mRNA translation. Our findings reveal a key role for Pbp1-containing assemblies in enabling the translation of mRNAs critical for mitochondrial biogenesis and respiration. They may further explain prior associations of Pbp1/ataxin-2 with RNA, stress granule biology, mitochondrial function, and neuronal health.

Males have finite resources to spend on reproduction. Thus, males rely on a 'time investment strategy' to maximize their reproductive success. For example, male Drosophila melanogaster extends their mating duration when surrounded by conditions enriched with rivals. Here we report a different form of behavioral plasticity whereby male fruit flies exhibit a shortened duration of mating when they are sexually experienced; we refer to this plasticity as 'shorter-mating-duration (SMD)'. SMD is a plastic behavior and requires sexually dimorphic taste neurons. We identified several neurons in the male foreleg and midleg that express specific sugar and pheromone receptors. Using a cost-benefit model and behavioral experiments, we further show that SMD behavior exhibits adaptive behavioral plasticity in male flies. Thus, our study delineates the molecular and cellular basis of the sensory inputs required for SMD; this represents a plastic interval timing behavior that could serve as a model system to study how multisensory inputs converge to modify interval timing behavior for improved adaptation.

Integration of light and phytohormones is essential for plant growth and development. FAR-RED INSENSITIVE 219 (FIN219)/JASMONATE RESISTANT 1 (JAR1) participates in phytochrome A (phyA)-mediated far-red (FR) light signaling in Arabidopsis and is a jasmonate (JA)-conjugating enzyme for the generation of an active JA-isoleucine. Accumulating evidence indicates that FR and JA signaling integrate with each other. However, the molecular mechanisms underlying their interaction remain largely unknown. Here, the phyA mutant was hypersensitive to JA. The double mutant fin219-2phyA-211 showed a synergistic effect on seedling development under FR light. Further evidence revealed that FIN219 and phyA antagonized with each other in a mutually functional demand to modulate hypocotyl elongation and expression of light- and JA-responsive genes. Moreover, FIN219 interacted with phyA under prolonged FR light, and MeJA could enhance their interaction with CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1) in the dark and FR light. FIN219 and phyA interaction occurred mainly in the cytoplasm, and they regulated their mutual subcellular localization under FR light. Surprisingly, the fin219-2 mutant abolished the formation of phyA nuclear bodies under FR light. Overall, these data identified a vital mechanism of phyA-FIN219-COP1 association in response to FR light, and MeJA may allow the photoactivated phyA to trigger photomorphogenic responses.

We report three novel deletions involving the Multispecies Conserved Sequences (MCS) R2, also known as the Major Regulative Element (MRE), in patients showing the α-thalassemia phenotype. The three new rearrangements showed peculiar positions of the breakpoints. 1) The (αα)ES is a telomeric 110 kb deletion ending inside the MCS-R3 element. 2) The (αα)FG, 984 bp-long, ends 51 bp upstream to MCS-R2; both are associated with a severe α-thalassemia phenotype. 3) The (αα)CT, 5058 bp-long starts at position +93 of MCS-R2 and is the only one associated to a mild α-thalassemia phenotype. To understand the specific role of different segments of the MCS-R2 element and of its boundary regions we carried out transcriptional and expression analysis. Transcriptional analysis of patients' reticulocytes showed that (αα)ES was unable to produce α2-globin mRNA, while a high level of expression of the α2-globin genes (56%) was detected in (αα)CT deletion, characterized by the presence of the first 93 bp of MCS-R2. Expression analysis of constructs containing breakpoints and boundary regions of the deletions (αα)CT and (αα)FG, showed comparable activity both for MCS-R2 and the boundary region (-682/-8). Considering that the (αα)CT deletion, almost entirely removing MCS-R2, has a less severe phenotype than the (αα)FG α0thalassemia deletion, removing both MCS-R2 almost entirely and an upstream 679 bp, we infer for the first time that an enhancer element must exist in this region that helps to increase the expression of the α-globin genes. The genotype-phenotype relationship of other previously published MCS-R2 deletions strengthened our hypothesis.

[This corrects the article DOI: 10.1371/journal.pgen.1008873.].

Heterozygous variants in the glucocerebrosidase (GBA) gene are common and potent risk factors for Parkinson's disease (PD). GBA also causes the autosomal recessive lysosomal storage disorder (LSD), Gaucher disease, and emerging evidence from human genetics implicates many other LSD genes in PD susceptibility. We have systemically tested 86 conserved fly homologs of 37 human LSD genes for requirements in the aging adult Drosophila brain and for potential genetic interactions with neurodegeneration caused by α-synuclein (αSyn), which forms Lewy body pathology in PD. Our screen identifies 15 genetic enhancers of αSyn-induced progressive locomotor dysfunction, including knockdown of fly homologs of GBA and other LSD genes with independent support as PD susceptibility factors from human genetics (SCARB2, SMPD1, CTSD, GNPTAB, SLC17A5). For several genes, results from multiple alleles suggest dose-sensitivity and context-dependent pleiotropy in the presence or absence of αSyn. Homologs of two genes causing cholesterol storage disorders, Npc1a / NPC1 and Lip4 / LIPA, were independently confirmed as loss-of-function enhancers of αSyn-induced retinal degeneration. The enzymes encoded by several modifier genes are upregulated in αSyn transgenic flies, based on unbiased proteomics, revealing a possible, albeit ineffective, compensatory response. Overall, our results reinforce the important role of lysosomal genes in brain health and PD pathogenesis, and implicate several metabolic pathways, including cholesterol homeostasis, in αSyn-mediated neurotoxicity.

Mendelian randomization (MR) has been increasingly applied for causal inference with observational data by using genetic variants as instrumental variables (IVs). However, the current practice of MR has been largely restricted to investigating the total causal effect between two traits, while it would be useful to infer the direct causal effect between any two of many traits (by accounting for indirect or mediating effects through other traits). For this purpose we propose a two-step approach: we first apply an extended MR method to infer (i.e. both estimate and test) a causal network of total effects among multiple traits, then we modify a graph deconvolution algorithm to infer the corresponding network of direct effects. Simulation studies showed much better performance of our proposed method than existing ones. We applied the method to 17 large-scale GWAS summary datasets (with median N = 256879 and median #IVs = 48) to infer the causal networks of both total and direct effects among 11 common cardiometabolic risk factors, 4 cardiometabolic diseases (coronary artery disease, stroke, type 2 diabetes, atrial fibrillation), Alzheimer's disease and asthma, identifying some interesting causal pathways. We also provide an R Shiny app (https://zhaotongl.shinyapps.io/cMLgraph/) for users to explore any subset of the 17 traits of interest.

Glucose is the preferred carbon source for most eukaryotes, and the first step in its metabolism is phosphorylation to glucose-6-phosphate. This reaction is catalyzed by hexokinases or glucokinases. The yeast Saccharomyces cerevisiae encodes three such enzymes, Hxk1, Hxk2, and Glk1. In yeast and mammals, some isoforms of this enzyme are found in the nucleus, suggesting a possible moonlighting function beyond glucose phosphorylation. In contrast to mammalian hexokinases, yeast Hxk2 has been proposed to shuttle into the nucleus in glucose-replete conditions, where it reportedly moonlights as part of a glucose-repressive transcriptional complex. To achieve its role in glucose repression, Hxk2 reportedly binds the Mig1 transcriptional repressor, is dephosphorylated at serine 15 and requires an N-terminal nuclear localization sequence (NLS). We used high-resolution, quantitative, fluorescent microscopy of live cells to determine the conditions, residues, and regulatory proteins required for Hxk2 nuclear localization. Countering previous yeast studies, we find that Hxk2 is largely excluded from the nucleus under glucose-replete conditions but is retained in the nucleus under glucose-limiting conditions. We find that the Hxk2 N-terminus does not contain an NLS but instead is necessary for nuclear exclusion and regulating multimerization. Amino acid substitutions of the phosphorylated residue, serine 15, disrupt Hxk2 dimerization but have no effect on its glucose-regulated nuclear localization. Alanine substation at nearby lysine 13 affects dimerization and maintenance of nuclear exclusion in glucose-replete conditions. Modeling and simulation provide insight into the molecular mechanisms of this regulation. In contrast to earlier studies, we find that the transcriptional repressor Mig1 and the protein kinase Snf1 have little effect on Hxk2 localization. Instead, the protein kinase Tda1 regulates Hxk2 localization. RNAseq analyses of the yeast transcriptome dispels the idea that Hxk2 moonlights as a transcriptional regulator of glucose repression, demonstrating that Hxk2 has a negligible role in transcriptional regulation in both glucose-replete and limiting conditions. Our studies define a new model of cis- and trans-acting regulators of Hxk2 dimerization and nuclear localization. Based on our data, the nuclear translocation of Hxk2 in yeast occurs in glucose starvation conditions, which aligns well with the nuclear regulation of mammalian orthologs. Our results lay the foundation for future studies of Hxk2 nuclear activity.

There are many mechanisms that give rise to genomic change: while point mutations are often emphasized in genomic analyses, evolution acts upon many other types of genetic changes that can result in less subtle perturbations. Changes in chromosome structure, DNA copy number, and novel transposon insertions all create large genomic changes, which can have correspondingly large impacts on phenotypes and fitness. In this study we investigate the spectrum of adaptive mutations that arise in a population under consistently fluctuating nitrogen conditions. We specifically contrast these adaptive alleles and the mutational mechanisms that create them, with mechanisms of adaptation under batch glucose limitation and constant selection in low, non-fluctuating nitrogen conditions to address if and how selection dynamics influence the molecular mechanisms of evolutionary adaptation. We observe that retrotransposon activity accounts for a substantial number of adaptive events, along with microhomology-mediated mechanisms of insertion, deletion, and gene conversion. In addition to loss of function alleles, which are often exploited in genetic screens, we identify putative gain of function alleles and alleles acting through as-of-yet unclear mechanisms. Taken together, our findings emphasize that how selection (fluctuating vs. non-fluctuating) is applied also shapes adaptation, just as the selective pressure (nitrogen vs. glucose) does itself. Fluctuating environments can activate different mutational mechanisms, shaping adaptive events accordingly. Experimental evolution, which allows a wider array of adaptive events to be assessed, is thus a complementary approach to both classical genetic screens and natural variation studies to characterize the genotype-to-phenotype-to-fitness map.

The rice blast fungus Magnaporthe oryzae forms specialized infectious structures called appressoria that breach host cells to initiate infection. Previous studies demonstrated that the regulator of G-protein signaling (RGS)-like protein MoRgs7 undergoes endocytosis upon fungal sensing of hydrophobic environmental cues to activate cAMP signaling required for appressorium formation. However, the mechanism by which MoRgs7 internalizes and its fate remains undetermined. We here show that MoSep1, a conserved protein kinase of Mitotic Exit Network (MEN), phosphorylates MoRgs7 to regulate its function. MoRgs7 phosphorylation determines its interaction with MoCrn1, a coronin-like actin-binding protein homolog that also modulates the internalization of MoRgs7. Importantly, the endocytic transport of MoRgs7 is critical for its GTPase-activating protein (GAP) function important in cAMP signaling. Together, our findings revealed a novel mechanism by which M. oryzae activates MoRgs7-mediated hydrophobic cue-sensing signal transduction involving protein phosphorylation and endocytic transport to govern appressorium formation and fungal pathogenicity.

An essential process during Danio rerio's left-right organizer (Kupffer's Vesicle, KV) formation is the formation of a motile cilium by developing KV cells which extends into the KV lumen. Beating of motile cilia within the KV lumen directs fluid flow to establish the embryo's left-right axis. However, the timepoint at which KV cells start to form cilia and how cilia formation is coordinated with KV lumen formation have not been examined. We identified that nascent KV cells form cilia at their centrosomes at random intracellular positions that then move towards a forming apical membrane containing cystic fibrosis transmembrane conductance regulator (CFTR). Using optogenetic clustering approaches, we found that Rab35 positive membranes recruit Rab11 to modulate CFTR delivery to the apical membrane, which is required for lumen opening, and subsequent cilia extension into the lumen. Once the intracellular cilia reach the CFTR positive apical membrane, Arl13b-positive cilia extend and elongate in a Rab8 dependent manner into the forming lumen once the lumen reaches an area of 300 μm2. These studies demonstrate the need to acutely coordinate Rab8, Rab11, and Rab35-mediated membrane trafficking events to ensure appropriate timing in lumen and cilia formation during KV development.

Plasmids are diverse extrachromosomal elements significantly that contribute to interspecies dissemination of antimicrobial resistance (AMR) genes. However, within clinically important bacteria, plasmids can exhibit unexpected narrow host ranges, a phenomenon that has scarcely been examined. Here we show that pConj is largely restricted to the human-specific pathogen, Neisseria gonorrhoeae. pConj can confer tetracycline resistance and is central to the dissemination of other AMR plasmids. We tracked pConj evolution from the pre-antibiotic era 80 years ago to the modern day and demonstrate that, aside from limited gene acquisition and loss events, pConj is remarkably conserved. Notably, pConj has remained prevalent in gonococcal populations despite cessation of tetracycline use, thereby demonstrating pConj adaptation to its host. Equally, pConj imposes no measurable fitness costs and is stably inherited by the gonococcus. Its maintenance depends on the co-operative activity of plasmid-encoded Toxin:Antitoxin (TA) and partitioning systems rather than host factors. An orphan VapD toxin encoded on pConj forms a split TA with antitoxins expressed from an ancestral co-resident plasmid or a horizontally-acquired chromosomal island, potentially explaining pConj's limited distribution. Finally, ciprofloxacin can induce loss of this highly stable plasmid, reflecting epidemiological evidence of transient reduction in pConj prevalence when fluoroquinolones were introduced to treat gonorrhoea.

Curli amyloid fibers are a major constituent of the extracellular biofilm matrix formed by bacteria of the Enterobacteriaceae family. Within Escherichia coli biofilms, curli gene expression is limited to a subpopulation of bacteria, leading to heterogeneity of extracellular matrix synthesis. Here we show that bimodal activation of curli gene expression also occurs in well-mixed planktonic cultures of E. coli, resulting in all-or-none stochastic differentiation into distinct subpopulations of curli-positive and curli-negative cells at the entry into the stationary phase of growth. Stochastic curli activation in individual E. coli cells could further be observed during continuous growth in a conditioned medium in a microfluidic device, which further revealed that the curli-positive state is only metastable. In agreement with previous reports, regulation of curli gene expression by the second messenger c-di-GMP via two pairs of diguanylate cyclase and phosphodiesterase enzymes, DgcE/PdeH and DgcM/PdeR, modulates the fraction of curli-positive cells. Unexpectedly, removal of this regulatory network does not abolish the bimodality of curli gene expression, although it affects dynamics of activation and increases heterogeneity of expression levels among individual cells. Moreover, the fraction of curli-positive cells within an E. coli population shows stronger dependence on growth conditions in the absence of regulation by DgcE/PdeH and DgcM/PdeR pairs. We thus conclude that, while not required for the emergence of bimodal curli gene expression in E. coli, this c-di-GMP regulatory network attenuates the frequency and dynamics of gene activation and increases its robustness to cellular heterogeneity and environmental variation.

Vibrio cholerae is a facultative pathogen that primarily occupies marine environments. In this niche, V. cholerae commonly interacts with the chitinous shells of crustacean zooplankton. As a chitinolytic microbe, V. cholerae degrades insoluble chitin into soluble oligosaccharides. Chitin oligosaccharides serve as both a nutrient source and an environmental cue that induces a strong transcriptional response in V. cholerae. Namely, these oligosaccharides induce the chitin sensor, ChiS, to activate the genes required for chitin utilization and horizontal gene transfer by natural transformation. Thus, interactions with chitin impact the survival of V. cholerae in marine environments. Chitin is a complex carbon source for V. cholerae to degrade and consume, and the presence of more energetically favorable carbon sources can inhibit chitin utilization. This phenomenon, known as carbon catabolite repression (CCR), is mediated by the glucose-specific Enzyme IIA (EIIAGlc) of the phosphoenolpyruvate-dependent phosphotransferase system (PTS). In the presence of glucose, EIIAGlc becomes dephosphorylated, which inhibits ChiS transcriptional activity by an unknown mechanism. Here, we show that dephosphorylated EIIAGlc interacts with ChiS. We also isolate ChiS suppressor mutants that evade EIIAGlc-dependent repression and demonstrate that these alleles no longer interact with EIIAGlc. These findings suggest that EIIAGlc must interact with ChiS to exert its repressive effect. Importantly, the ChiS suppressor mutations we isolated also relieve repression of chitin utilization and natural transformation by EIIAGlc, suggesting that CCR of these behaviors is primarily regulated through ChiS. Together, our results reveal how nutrient conditions impact the fitness of an important human pathogen in its environmental reservoir.

Neurite pruning and regrowth are important mechanisms to adapt neural circuits to distinct developmental stages. Neurite regrowth after pruning often depends on differential regulation of growth signaling pathways, but their precise mechanisms of action during regrowth are unclear. Here, we show that the PI3K/TORC1 pathway is required for dendrite regrowth after pruning in Drosophila peripheral neurons during metamorphosis. TORC1 impinges on translation initiation, and our analysis of 5' untranslated regions (UTRs) of remodeling factor mRNAs linked to actin suggests that TOR selectively stimulates the translation of regrowth over pruning factors. Furthermore, we find that dendrite regrowth also requires the GTPase RalA and the exocyst complex as regulators of polarized secretion, and we provide evidence that this pathway is also regulated by TOR. We propose that TORC1 coordinates dendrite regrowth after pruning by coordinately stimulating the translation of regrowth factors involved in cytoskeletal regulation and secretion.

Stem cell differentiation is a highly dynamic process involving pervasive changes in gene expression. The large majority of existing studies has characterized differentiation at the level of individual molecular profiles, such as the transcriptome or the proteome. To obtain a more comprehensive view, we measured protein, mRNA and microRNA abundance during retinoic acid-driven differentiation of mouse embryonic stem cells. We found that mRNA and protein abundance are typically only weakly correlated across time. To understand this finding, we developed a hierarchical dynamical model that allowed us to integrate all data sets. This model was able to explain mRNA-protein discordance for most genes and identified instances of potential microRNA-mediated regulation. Overexpression or depletion of microRNAs identified by the model, followed by RNA sequencing and protein quantification, were used to follow up on the predictions of the model. Overall, our study shows how multi-omics integration by a dynamical model could be used to nominate candidate regulators.

Protein degradation is an essential biological process that regulates protein abundance and removes misfolded and damaged proteins from cells. In eukaryotes, most protein degradation occurs through the stepwise actions of two functionally distinct entities, the ubiquitin system and the proteasome. Ubiquitin system enzymes attach ubiquitin to cellular proteins, targeting them for degradation. The proteasome then selectively binds and degrades ubiquitinated substrate proteins. Genetic variation in ubiquitin system genes creates heritable differences in the degradation of their substrates. However, the challenges of measuring the degradative activity of the proteasome independently of the ubiquitin system in large samples have limited our understanding of genetic influences on the proteasome. Here, using the yeast Saccharomyces cerevisiae, we built and characterized reporters that provide high-throughput, ubiquitin system-independent measurements of proteasome activity. Using single-cell measurements of proteasome activity from millions of genetically diverse yeast cells, we mapped 15 loci across the genome that influence proteasomal protein degradation. Twelve of these 15 loci exerted specific effects on the degradation of two distinct proteasome substrates, revealing a high degree of substrate-specificity in the genetics of proteasome activity. Using CRISPR-Cas9-based allelic engineering, we resolved a locus to a causal variant in the promoter of RPT6, a gene that encodes a subunit of the proteasome's 19S regulatory particle. The variant increases RPT6 expression, which we show results in increased proteasome activity. Our results reveal the complex genetic architecture of proteasome activity and suggest that genetic influences on the proteasome may be an important source of variation in the many cellular and organismal traits shaped by protein degradation.

Transposable elements constitute nearly half of the mammalian genome and play important roles in genome evolution. While a multitude of both transcriptional and post-transcriptional mechanisms exist to silence transposable elements, control of transposition in vivo remains poorly understood. MOV10, an RNA helicase, is an inhibitor of mobilization of retrotransposons and retroviruses in cell culture assays. Here we report that MOV10 restricts LINE1 retrotransposition in mice. Although MOV10 is broadly expressed, its loss causes only incomplete penetrance of embryonic lethality, and the surviving MOV10-deficient mice are healthy and fertile. Biochemically, MOV10 forms a complex with UPF1, a key component of the nonsense-mediated mRNA decay pathway, and primarily binds to the 3' UTR of somatically expressed transcripts in testis. Consequently, loss of MOV10 results in an altered transcriptome in testis. Analyses using a LINE1 reporter transgene reveal that loss of MOV10 leads to increased LINE1 retrotransposition in somatic and reproductive tissues from both embryos and adult mice. Moreover, the degree of LINE1 retrotransposition inhibition is dependent on the Mov10 gene dosage. Furthermore, MOV10 deficiency reduces reproductive fitness over successive generations. Our findings demonstrate that MOV10 attenuates LINE1 retrotransposition in a dosage-dependent manner in mice.