Metabolic pathways are known to sense the environmental stimuli and result in physiological adjustments. The responding processes need to be tightly controlled. Here, we show that upon encountering P. aeruginosa, C. elegans upregulate the transcription factor ets-4, but this upregulation is attenuated by the ribonuclease, rege-1. As such, mutants with defective REGE-1 ribonuclease activity undergo ets-4-dependent early death upon challenge with P. aeruginosa. Furthermore, mRNA-seq analysis revealed associated global changes in two key metabolic pathways, the IIS (insulin/IGF signaling) and TOR (target of rapamycin) kinase signaling pathways. In particular, failure to degrade ets-4 mRNA in activity-defective rege-1 mutants resulted in upregulation of class II longevity genes, which are suppressed during longevity, and activation of TORC1 kinase signaling pathway. Genetic inhibition of either pathway way was sufficient to abolish the poor survival phenotype in rege-1 worms. Further analysis of ETS-4 ChIP data from ENCODE and characterization of one upregulated class II gene, ins-7, support that the Class II genes are activated by ETS-4. Interestingly, deleting an upregulated Class II gene, acox-1.5, a peroxisome β-oxidation enzyme, largely rescues the fat lost phenotype and survival difference between rege-1 mutants and wild-types. Thus, rege-1 appears to be crucial for animal survival due to its tight regulation of physiological responses to environmental stimuli. This function is reminiscent of its mammalian ortholog, Regnase-1, which modulates the intestinal mTORC1 signaling pathway.

The double stranded RNA binding protein Adad1 (adenosine deaminase domain containing 1) is a member of the adenosine deaminase acting on RNAs (Adar) protein family with germ cell-specific expression. In mice, Adad1 is necessary for sperm differentiation, however its function outside of mammals has not been investigated. Here, through an N-ethyl-N-nitrosourea (ENU) based forward genetic screen, we identified an adad1 mutant zebrafish line that develops as sterile males. Further histological examination revealed complete lack of germ cells in adult mutant fish, however germ cells populated the gonad, proliferated, and entered meiosis in larval and juvenile fish. Although meiosis was initiated in adad1 mutant testes, the spermatocytes failed to progress beyond the zygotene stage. Thus, Adad1 is essential for meiosis and germline maintenance in zebrafish. We tested if spermatogonial stem cells were affected using nanos2 RNA FISH and a label retaining cell (LRC) assay, and found that the mutant testes had fewer LRCs and nanos2-expressing cells compared to wild-type siblings, suggesting that failure to maintain the spermatogonial stem cells resulted in germ cell loss by adulthood. To identify potential molecular processes regulated by Adad1, we sequenced bulk mRNA from mutants and wild-type testes and found mis-regulation of genes involved in RNA stability and modification, pointing to a potential broader role in post-transcriptional regulation. Our findings suggest that the RNA regulatory protein Adad1 is required for fertility through regulation of spermatogonial stem cell maintenance in zebrafish.

Reshuffling of genetic variation occurs both by independent assortment of chromosomes and by homologous recombination. Such reshuffling can generate novel allele combinations and break linkage between advantageous and deleterious variants which increases both the potential and the efficacy of natural selection. Here we used high-density linkage maps to characterize global and regional recombination rate variation in two populations of the wood white butterfly (Leptidea sinapis) that differ considerably in their karyotype as a consequence of at least 27 chromosome fissions and fusions. The recombination data were compared to estimates of genetic diversity and measures of selection to assess the relationship between chromosomal rearrangements, crossing over, maintenance of genetic diversity and adaptation. Our data show that the recombination rate is influenced by both chromosome size and number, but that the difference in the number of crossovers between karyotypes is reduced as a consequence of a higher frequency of double crossovers in larger chromosomes. As expected from effects of selection on linked sites, we observed an overall positive association between recombination rate and genetic diversity in both populations. Our results also revealed a significant effect of chromosomal rearrangements on the rate of intergenic diversity change between populations, but limited effects on polymorphisms in coding sequence. We conclude that chromosomal rearrangements can have considerable effects on the recombination landscape and consequently influence both maintenance of genetic diversity and efficiency of selection in natural populations.

The integration of individually replicating genes into a primitive chromosome is a key evolutionary transition in the development of life, allowing the simultaneous inheritance of genes. However, how this transition occurred is unclear because the extended size of primitive chromosomes replicate slower than unlinked genes. Theoretical studies have suggested that a primitive chromosome can evolve in the presence of cell-like compartments, as the physical linkage prevents the stochastic loss of essential genes upon division, but experimental support for this is lacking. Here, we demonstrate the evolution of a chromosome-like RNA from two cooperative RNA replicators encoding replication and metabolic enzymes. Through their long-term replication in cell-like compartments, linked RNAs emerged with the two cooperative RNAs connected end-to-end. The linked RNAs had different mutation patterns than the two unlinked RNAs, suggesting that they were maintained as partially distinct lineages in the population. Our results provide experimental evidence supporting the plausibility of the evolution of a primitive chromosome from unlinked gene fragments, an important step in the emergence of complex biological systems.

Trichoderma spp. are ubiquitous rhizosphere fungi capable of producing several classes of secondary metabolites that can modify the dynamics of the plant-associated microbiome. However, the bacterial-fungal mechanisms that mediate these interactions have not been fully characterized. Here, a random barcode transposon-site sequencing (RB-TnSeq) approach was employed to identify bacterial genes important for fitness in the presence of Trichoderma atroviride exudates. We selected three rhizosphere bacteria with RB-TnSeq mutant libraries that can promote plant growth: the nitrogen fixers Klebsiella michiganensis M5aI and Herbaspirillum seropedicae SmR1, and Pseudomonas simiae WCS417. As a non-rhizosphere species, Pseudomonas putida KT2440 was also included. From the RB-TnSeq data, nitrogen-fixing bacteria competed mainly for iron and required the siderophore transport system TonB/Exb for optimal fitness in the presence of T. atroviride exudates. In contrast, P. simiae and P. putida were highly dependent on mechanisms associated with membrane lipid modification that are required for resistance to cationic antimicrobial peptides (CAMPs). A mutant in the Hog1-MAP kinase (Δtmk3) gene of T. atroviride showed altered expression patterns of many nonribosomal peptide synthetase (NRPS) biosynthetic gene clusters with potential antibiotic activity. In contrast with exudates from wild-type T. atroviride, bacterial mutants containing lesions in genes associated with resistance to antibiotics did not show fitness defects when RB-TnSeq libraries were exposed to exudates from the Δtmk3 mutant. Unexpectedly, exudates from wild-type T. atroviride and the Δtmk3 mutant rescued purine auxotrophic mutants of H. seropedicae, K. michiganensis and P. simiae. Metabolomic analysis on exudates from wild-type T. atroviride and the Δtmk3 mutant showed that both strains excrete purines and complex metabolites; functional Tmk3 is required to produce some of these metabolites. This study highlights the complex interplay between Trichoderma-metabolites and soil bacteria, revealing both beneficial and antagonistic effects, and underscoring the intricate and multifaceted nature of this relationship.

The emergence of antibiotic resistance under treatment depends on the availability of resistance alleles and their establishment in the population. Novel resistance alleles are encoded either in chromosomal or extrachromosomal genetic elements; both types may be present in multiple copies within the cell. However, the effect of polyploidy on the emergence of antibiotic resistance remains understudied. Here we show that the establishment of resistance alleles in microbial populations depends on the ploidy level. Evolving bacterial populations under selection for antibiotic resistance, we demonstrate that resistance alleles in polyploid elements are lost frequently in comparison to alleles in monoploid elements due to segregational drift. Integrating the experiments with a mathematical model, we find a remarkable agreement between the theoretical and empirical results, confirming our understanding of the allele segregation process. Using the mathematical model, we further show that the effect of polyploidy on the establishment probability of beneficial alleles is strongest for low replicon copy numbers and plateaus for high replicon copy numbers. Our results suggest that the distribution of fitness effects for mutations that are eventually fixed in a population depends on the replicon ploidy level. Our study indicates that the emergence of antibiotic resistance in bacterial pathogens depends on the pathogen ploidy level.

Rapid depletion of cellular ATP can occur by oxidative stress induced by reactive oxygen species (ROS). Maintaining energy homeostasis requires the key molecular components AMP-activated protein kinase (AMPK) and arginine kinase (AK), an invertebrate orthologue of the mammalian creatine kinase (CK). Here, we deciphered two independent and synergistic pathways of AMPK acting on AK by using the beetle Tribolium castaneum as a model system. First, AMPK acts on transcriptional factor forkhead box O (FOXO) leading to phosphorylation and nuclear translocation of the FOXO. The phospho-FOXO directly promotes the expression of AK upon oxidative stress. Concomitantly, AMPK directly phosphorylates the AK to switch the direction of enzymatic catalysis for rapid production of ATP from the phosphoarginine-arginine pool. Further in vitro assays revealed that Sf9 cells expressing phospho-deficient AK mutants displayed the lower ATP/ADP ratio and cell viability under paraquat-induced oxidative stress conditions when compared with Sf9 cells expressing wild-type AKs. Additionally, the AMPK-FOXO-CK pathway is also involved in the restoration of ATP homeostasis under oxidative stress in mammalian HEK293 cells. Overall, we provide evidence that two distinct AMPK-AK pathways, transcriptional and post-translational regulations, are coherent responders to acute oxidative stresses and distinguished from classical AMPK-mediated long-term metabolic adaptations to energy challenge.

Suppression of transposable elements (TEs) is paramount to maintain genomic integrity and organismal fitness. In D. melanogaster, the flamenco locus is a master suppressor of TEs, preventing the mobilization of certain endogenous retrovirus-like TEs from somatic ovarian support cells to the germline. It is transcribed by Pol II as a long (100s of kb), single-stranded, primary transcript, and metabolized into ~24-32 nt Piwi-interacting RNAs (piRNAs) that target active TEs via antisense complementarity. flamenco is thought to operate as a trap, owing to its high content of recent horizontally transferred TEs that are enriched in antisense orientation. Using newly-generated long read genome data, which is critical for accurate assembly of repetitive sequences, we find that flamenco has undergone radical transformations in sequence content and even copy number across simulans clade Drosophilid species. Drosophila simulans flamenco has duplicated and diverged, and neither copy exhibits synteny with D. melanogaster beyond the core promoter. Moreover, flamenco organization is highly variable across D. simulans individuals. Next, we find that D. simulans and D. mauritiana flamenco display signatures of a dual-stranded cluster, with ping-pong signals in the testis and/or embryo. This is accompanied by increased copy numbers of germline TEs, consistent with these regions operating as functional dual-stranded clusters. Overall, the physical and functional diversity of flamenco orthologs is testament to the extremely dynamic consequences of TE arms races on genome organization, not only amongst highly related species, but even amongst individuals.

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

Polo-like kinase 1 (PLK1) is a serine/threonine kinase required for mitosis and cytokinesis. As cancer cells are often hypersensitive to partial PLK1 inactivation, chemical inhibitors of PLK1 have been developed and tested in clinical trials. However, these small molecule inhibitors alone are not completely effective. PLK1 promotes numerous molecular and cellular events in the cell division cycle and it is unclear which of these events most crucially depend on PLK1 activity. We used a CRISPR-based genome-wide screening strategy to identify genes whose inactivation enhances cell proliferation defects upon partial chemical inhibition of PLK1. Genes identified encode proteins that are functionally linked to PLK1 in multiple ways, most notably factors that promote centromere and kinetochore function. Loss of the kinesin KIF18A or the outer kinetochore protein SKA1 in PLK1-compromised cells resulted in mitotic defects, activation of the spindle assembly checkpoint and nuclear reassembly defects. We also show that PLK1-dependent CENP-A loading at centromeres is extremely sensitive to partial PLK1 inhibition. Our results suggest that partial inhibition of PLK1 compromises the integrity and function of the centromere/kinetochore complex, rendering cells hypersensitive to different kinetochore perturbations. We propose that KIF18A is a promising target for combinatorial therapies with PLK1 inhibitors.

Transcription of ribosomal RNA (rRNA) by RNA Polymerase (Pol) I in the nucleolus is necessary for ribosome biogenesis, which is intimately tied to cell growth and proliferation. Perturbation of ribosome biogenesis results in tissue specific disorders termed ribosomopathies in association with alterations in nucleolar structure. However, how rRNA transcription and ribosome biogenesis regulate nucleolar structure during normal development and in the pathogenesis of disease remains poorly understood. Here we show that homozygous null mutations in Pol I subunits required for rRNA transcription and ribosome biogenesis lead to preimplantation lethality. Moreover, we discovered that Polr1a-/-, Polr1b-/-, Polr1c-/- and Polr1d-/- mutants exhibit defects in the structure of their nucleoli, as evidenced by a decrease in number of nucleolar precursor bodies and a concomitant increase in nucleolar volume, which results in a single condensed nucleolus. Pharmacological inhibition of Pol I in preimplantation and midgestation embryos, as well as in hiPSCs, similarly results in a single condensed nucleolus or fragmented nucleoli. We find that when Pol I function and rRNA transcription is inhibited, the viscosity of the granular compartment of the nucleolus increases, which disrupts its phase separation properties, leading to a single condensed nucleolus. However, if a cell progresses through mitosis, the absence of rRNA transcription prevents reassembly of the nucleolus and manifests as fragmented nucleoli. Taken together, our data suggests that Pol I function and rRNA transcription are required for maintaining nucleolar structure and integrity during development and in the pathogenesis of disease.

The conserved transcription factor Myc regulates cell growth, proliferation and apoptosis, and its deregulation has been associated with human pathologies. Although specific miRNAs have been identified as fundamental components of the Myc tumorigenic program, how Myc regulates miRNA biogenesis remains controversial. Here we showed that Myc functions as an important regulator of miRNA biogenesis in Drosophila by influencing both miRNA gene expression and processing. Through the analysis of ChIP-Seq datasets, we discovered that nearly 56% of Drosophila miRNA genes show dMyc binding, exhibiting either the canonical or non-canonical E-box sequences within the peak region. Consistently, reduction of dMyc levels resulted in widespread downregulation of miRNAs gene expression. dMyc also modulates miRNA processing and activity by controlling Drosha and AGO1 levels through direct transcriptional regulation. By using in vivo miRNA activity sensors we demonstrated that dMyc promotes miRNA-mediated silencing in different tissues, including the wing primordium and the fat body. We also showed that dMyc-dependent expression of miR-305 in the fat body modulates Dmp53 levels depending on nutrient availability, having a profound impact on the ability of the organism to respond to nutrient stress. Indeed, dMyc depletion in the fat body resulted in extended survival to nutrient deprivation which was reverted by expression of either miR-305 or a dominant negative version of Dmp53. Our study reveals a previously unrecognized function of dMyc as an important regulator of miRNA biogenesis and suggests that Myc-dependent expression of specific miRNAs may have important tissue-specific functions.

The mammalian cochlea is composed of sensory hair cells as well as multiple different types of non-sensory supporting cells. Pillar cells are one type of supporting cell that form the tunnel of Corti and include two morphologically and functionally distinct subtypes: inner pillar cells (IPCs) and outer pillar cells (OPCs). The processes of specification and differentiation of inner versus outer pillar cells are still unclear. Here, we show that β-Catenin is required for establishing IPC identity in the mammalian cochlea. To differentiate the transcriptional and adhesion roles of β-Catenin in establishing IPC identity, we examined two different models of β-Catenin deletion; one that deletes both transcriptional and structural functions and one which retains cell adhesion function but lacks transcriptional function. Here, we show that cochleae lacking β-Catenin transcriptional function lost IPCs and displayed extranumerary OPCs, indicating its requirement for establishing IPC identity. Overexpression of β-Catenin induced proliferation within IPCs but not ectopic IPCs. Single-cell transcriptomes of supporting cells lacking β-Catenin transcriptional function show a loss of the IPC and gain of OPC signatures. Finally, targeted deletion of β-Catenin in IPCs also led to the loss of IPC identity, indicating a cell autonomous role of β-Catenin in establishing IPC identity. As IPCs have the capacity to regenerate sensory hair cells in the postnatal cochlea, our results will aid in future IPC-based hair cell regeneration strategies.

Localization of oskar mRNA to the posterior of the Drosophila oocyte is essential for abdominal patterning and germline development. oskar localization is a multi-step process involving temporally and mechanistically distinct transport modes. Numerous cis-acting elements and trans-acting factors have been identified that mediate earlier motor-dependent transport steps leading to an initial accumulation of oskar at the posterior. Little is known, however, about the requirements for the later localization phase, which depends on cytoplasmic flows and results in the accumulation of large oskar ribonucleoprotein granules, called founder granules, by the end of oogenesis. Using super-resolution microscopy, we show that founder granules are agglomerates of smaller oskar transport particles. In contrast to the earlier kinesin-dependent oskar transport, late-phase localization depends on the sequence as well as on the structure of the spliced oskar localization element (SOLE), but not on the adjacent exon junction complex deposition. Late-phase localization also requires the oskar 3' untranslated region (3' UTR), which targets oskar to founder granules. Together, our results show that 3' UTR-mediated targeting together with SOLE-dependent agglomeration leads to accumulation of oskar in large founder granules at the posterior of the oocyte during late stages of oogenesis. In light of previous work showing that oskar transport particles are solid-like condensates, our findings indicate that founder granules form by a process distinct from that of well-characterized ribonucleoprotein granules like germ granules, P bodies, and stress granules. Additionally, they illustrate how an individual mRNA can be adapted to exploit different localization mechanisms depending on the cellular context.

Striated muscle laminopathies caused by missense mutations in the nuclear lamin gene LMNA are characterized by cardiac dysfunction and often skeletal muscle defects. Attempts to predict which LMNA variants are pathogenic and to understand their physiological effects lag behind variant discovery. We created Caenorhabditis elegans models for striated muscle laminopathies by introducing pathogenic human LMNA variants and variants of unknown significance at conserved residues within the lmn-1 gene. Severe missense variants reduced fertility and/or motility in C. elegans. Nuclear morphology defects were evident in the hypodermal nuclei of many lamin variant strains, indicating a loss of nuclear envelope integrity. Phenotypic severity varied within the two classes of missense mutations involved in striated muscle disease, but overall, variants associated with both skeletal and cardiac muscle defects in humans lead to more severe phenotypes in our model than variants predicted to disrupt cardiac function alone. We also identified a separation of function allele, lmn-1(R204W), that exhibited normal viability and swimming behavior but had a severe nuclear migration defect. Thus, we established C. elegans avatars for striated muscle laminopathies and identified LMNA variants that offer insight into lamin mechanisms during normal development.

Quantitative traits may be controlled by many loci, many alleles at each locus, and subject to genotype-by-environment interactions, making them difficult to map. One example of such a complex trait is shoot branching in the model plant Arabidopsis, and its plasticity in response to nitrate. Here, we use artificial selection under contrasting nitrate supplies to dissect the genetic architecture of this complex trait, where loci identified by association mapping failed to explain heritability estimates. We found a consistent response to selection for high branching, with correlated responses in other traits such as plasticity and flowering time. Genome-wide scans for selection and simulations suggest that at least tens of loci control this trait, with a distinct genetic architecture between low and high nitrate treatments. While signals of selection could be detected in the populations selected for high branching on low nitrate, there was very little overlap in the regions selected in three independent populations. Thus the regulatory network controlling shoot branching can be tuned in different ways to give similar phenotypes.

Persons living with HIV (PLWH) have an increased risk for tuberculosis (TB). After prolonged and repeated exposure, some PLWH never develop TB and show no evidence of immune sensitization to Mycobacterium tuberculosis (Mtb) as defined by persistently negative tuberculin skin tests (TST) and interferon gamma release assays (IGRA). This group has been identified and defined as HIV+ persistently TB, tuberculin and IGRA negative (HITTIN). To investigate potential innate mechanisms unique to individuals with the HITTIN phenotype we compared their neutrophil Mtb infection response to that of PLWH, with no TB history, but who test persistently IGRA positive, and tuberculin positive (HIT). Neutrophil samples from 17 HITTIN (PMNHITTIN) and 11 HIT (PMNHIT) were isolated and infected with Mtb H37Rv for 1h and 6h. RNA was extracted and used for RNAseq analysis. Since there was no significant differential transcriptional response at 1h between infected PMNHITTIN and PMNHIT, we focused on the 6h timepoint. When compared to uninfected PMN, PMNHITTIN displayed 3106 significantly upregulated and 3548 significantly downregulated differentially expressed genes (DEGs) (absolute cutoff of a log2FC of 0.2, FDR < 0.05) whereas PMNHIT demonstrated 3816 significantly upregulated and 3794 significantly downregulated DEGs following 6h Mtb infection. Contrasting the log2FC 6h infection response to Mtb from PMNHITTIN against PMNHIT, 2285 genes showed significant differential response between the two groups. Overall PMNHITTIN had a lower fold change response to Mtb infection compared to PMNHIT. According to pathway enrichment, Apoptosis and NETosis were differentially regulated between HITTIN and HIT PMN responses after 6h Mtb infection. To corroborate the blunted NETosis transcriptional response measured among HITTIN, fluorescence microscopy revealed relatively lower neutrophil extracellular trap formation and cell loss in PMNHITTIN compared to PMNHIT, showing that PMNHITTIN have a distinct response to Mtb.

Regulation of the microtubule cytoskeleton is crucial for the development and maintenance of neuronal architecture, and recent studies have highlighted the significance of regulated RNA processing in the establishment and maintenance of neural circuits. In a genetic screen conducted using mechanosensory neurons of C. elegans, we identified a mutation in muscleblind-1/mbl-1 as a suppressor of loss of kinesin-13 family microtubule destabilizing factor klp-7. Muscleblind-1(MBL-1) is an RNA-binding protein that regulates the splicing, localization, and stability of RNA. Our findings demonstrate that mbl-1 is required cell-autonomously for axon growth and proper synapse positioning in the posterior lateral microtubule (PLM) neuron. Loss of mbl-1 leads to increased microtubule dynamics and mixed orientation of microtubules in the anterior neurite of PLM. These defects are also accompanied by abnormal axonal transport of the synaptic protein RAB-3 and reduction of gentle touch sensation in mbl-1 mutant. Our data also revealed that mbl-1 is genetically epistatic to mec-7 (β tubulin) and mec-12 (α tubulin) in regulating axon growth. Furthermore, mbl-1 is epistatic to sad-1, an ortholog of BRSK/Brain specific-serine/threonine kinase and a known regulator of synaptic machinery, for synapse formation at the correct location of the PLM neurite. Notably, the immunoprecipitation of MBL-1 resulted in the co-purification of mec-7, mec-12, and sad-1 mRNAs, suggesting a direct interaction between MBL-1 and these transcripts. Additionally, mbl-1 mutants exhibited reduced levels and stability of mec-7 and mec-12 transcripts. Our study establishes a previously unknown link between RNA-binding proteins and cytoskeletal machinery, highlighting their crucial roles in the development and maintenance of the nervous system.

The mitochondrial C-to-U RNA editing factor PPR56 of the moss Physcomitrium patens is an RNA-binding pentatricopeptide repeat protein equipped with a terminal DYW-type cytidine deaminase domain. Transferred into Escherichia coli, PPR56 works faithfully on its two native RNA editing targets, nad3eU230SL and nad4eU272SL, and also converts cytidines into uridines at over 100 off-targets in the bacterial transcriptome. Accordingly, PPR56 is attractive for detailed mechanistic studies in the heterologous bacterial setup, allowing for scoring differential RNA editing activities of many target and protein variants in reasonable time. Here, we report (i) on the effects of numerous individual and combined PPR56 protein and target modifications, (ii) on the spectrum of off-target C-to-U editing in the bacterial background transcriptome for PPR56 and two variants engineered for target re-direction and (iii) on combinations of targets in tandem or separately at the 5'- and 3'-ends of large mRNAs. The latter experimentation finds enhancement of RNA editing at weak targets in many cases, including cox3eU290SF as a new candidate mitogenome target. We conclude that C-to-U RNA editing can be much enhanced by transcript features also outside the region ultimately targeted by PPRs of a plant editing factor, possibly facilitated by its enrichment or scanning along transcripts.

Escherichia coli is both a highly prevalent commensal and a major opportunistic pathogen causing bloodstream infections (BSI). A systematic analysis characterizing the genomic determinants of extra-intestinal pathogenic vs. commensal isolates in human populations, which could inform mechanisms of pathogenesis, diagnostic, prevention and treatment is still lacking. We used a collection of 912 BSI and 370 commensal E. coli isolates collected in France over a 17-year period (2000-2017). We compared their pangenomes, genetic backgrounds (phylogroups, STs, O groups), presence of virulence-associated genes (VAGs) and antimicrobial resistance genes, finding significant differences in all comparisons between commensal and BSI isolates. A machine learning linear model trained on all the genetic variants derived from the pangenome and controlling for population structure reveals similar differences in VAGs, discovers new variants associated with pathogenicity (capacity to cause BSI), and accurately classifies BSI vs. commensal strains. Pathogenicity is a highly heritable trait, with up to 69% of the variance explained by bacterial genetic variants. Lastly, complementing our commensal collection with an older collection from 1980, we predict that pathogenicity continuously increased through 1980, 2000, to 2010. Together our findings imply that E. coli exhibit substantial genetic variation contributing to the transition between commensalism and pathogenicity and that this species evolved towards higher pathogenicity.

In the mammalian cerebral cortex, the hippocampal primordium (Hcp) occupies a discrete position in the dorsal telencephalic neuroepithelium adjacent to the neocortical primordium (Ncp). We examined transcriptomic and chromatin-level features that distinguish the Hcp from the Ncp in the mouse during the early neurogenic period, embryonic day (E)12.5. ATAC-seq revealed that the Hcp was more accessible than the Ncp at this stage. Motif analysis of the differentially accessible loci in these tissues revealed LHX2 as a candidate transcription factor for modulating gene regulatory networks (GRNs). We analyzed LHX2 occupancy profiles and compared these with transcriptomic data from control and Lhx2 mutant Hcp and Ncp at E12.5. Our results revealed that LHX2 directly regulates distinct genes in the Hcp and Ncp within a set of common pathways that control fundamental aspects of development namely pluripotency, axon pathfinding, Wnt, and Hippo signaling. Loss of Lhx2 caused a decrease in accessibility, specifically in hippocampal chromatin, suggesting that this factor may play a unique role in hippocampal development. We identified 14 genes that were preferentially enriched in the Hcp, for which LHX2 regulates both chromatin accessibility and mRNA expression, which have not thus far been examined in hippocampal development. Together, these results provide mechanistic insight into how LHX2 function in the Hcp may contribute to the process by which the hippocampus acquires features distinct from the neocortex.

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Assessment of the genetic similarity between two phenotypes can provide insight into a common genetic aetiology and inform the use of pleiotropy-informed, cross-phenotype analytical methods to identify novel genetic associations. The genetic correlation is a well-known means of quantifying and testing for genetic similarity between traits, but its estimates are subject to comparatively large sampling error. This makes it unsuitable for use in a small-sample context. We discuss the use of a previously published nonparametric test of genetic similarity for application to GWAS summary statistics. We establish that the null distribution of the test statistic is modelled better by an extreme value distribution than a transformation of the standard exponential distribution. We show with simulation studies and real data from GWAS of 18 phenotypes from the UK Biobank that the test is to be preferred for use with small sample sizes, particularly when genetic effects are few and large, outperforming the genetic correlation and another nonparametric statistical test of independence. We find the test suitable for the detection of genetic similarity in the rare disease context.

N-terminal ends of polypeptides are critical for the selective co-translational recruitment of N-terminal modification enzymes. However, it is unknown whether specific N-terminal signatures differentially regulate protein fate according to their cellular functions. In this work, we developed an in-silico approach to detect functional preferences in cellular N-terminomes, and identified in S. cerevisiae more than 200 Gene Ontology terms with specific N-terminal signatures. In particular, we discovered that Mitochondrial Targeting Sequences (MTS) show a strong and specific over-representation at position 2 of hydrophobic residues known to define potential substrates of the N-terminal acetyltransferase NatC. We validated mitochondrial precursors as co-translational targets of NatC by selective purification of translating ribosomes, and found that their N-terminal signature is conserved in Saccharomycotina yeasts. Finally, systematic mutagenesis of the position 2 in a prototypal yeast mitochondrial protein confirmed its critical role in mitochondrial protein import. Our work highlights the hydrophobicity of MTS N-terminal residues and their targeting by NatC as important features for the definition of the mitochondrial proteome, providing a molecular explanation for mitochondrial defects observed in yeast or human NatC-depleted cells. Functional mapping of N-terminal residues thus has the potential to support the discovery of novel mechanisms of protein regulation or targeting.

Diabetic retinopathy (DR) is a common complication of diabetes. Approximately 20% of DR patients have diabetic macular edema (DME) characterized by fluid leakage into the retina. There is a genetic component to DR and DME risk, but few replicable loci. Because not all DR cases have DME, we focused on DME to increase power, and conducted a multi-ancestry GWAS to assess DME risk in a total of 1,502 DME patients and 5,603 non-DME controls in discovery and replication datasets. Two loci reached GWAS significance (p<5x10-8). The strongest association was rs2239785, (K150E) in APOL1. The second finding was rs10402468, which co-localized to PLVAP and ANKLE1 in vascular / endothelium tissues. We conducted multiple sensitivity analyses to establish that the associations were specific to DME status and did not reflect diabetes status or other diabetic complications. Here we report two novel loci for risk of DME which replicated in multiple clinical trial and biobank derived datasets. One of these loci, containing the gene APOL1, is a risk factor in African American DME and DKD patients, indicating that this locus plays a broader role in diabetic complications for multiple ancestries. Trial Registration: NCT00473330, NCT00473382, NCT03622580, NCT03622593, NCT04108156.

Evidence of interbreeding between archaic hominins and humans comes from methods that infer the locations of segments of archaic haplotypes, or 'archaic coverage' using the genomes of people living today. As more estimates of archaic coverage have emerged, it has become clear that most of this coverage is found on the autosomes- very little is retained on chromosome X. Here, we summarize published estimates of archaic coverage on autosomes and chromosome X from extant human samples. We find on average 7 times more archaic coverage on autosomes than chromosome X, and identify broad continental patterns in this ratio: greatest in European samples, and least in South Asian samples. We also perform extensive simulation studies to investigate how the amount of archaic coverage, lengths of coverage, and rates of purging of archaic coverage are affected by sex-bias caused by an unequal sex ratio within the archaic introgressors. Our results generally confirm that, with increasing male sex-bias, less archaic coverage is retained on chromosome X. Ours is the first study to explicitly model such sex-bias and its potential role in creating the dearth of archaic coverage on chromosome X.

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Lipid droplets (LDs) are ubiquitous organelles that facilitate neutral lipid storage in cells, including energy-dense triglycerides. They are found in all investigated metazoan embryos where they are thought to provide energy for development. Intriguingly, early embryos of diverse metazoan species asymmetrically allocate LDs amongst cellular lineages, a process which can involve massive intracellular redistribution of LDs. However, the biological reason for asymmetric lineage allocation is unknown. To address this issue, we utilize the Drosophila embryo where the cytoskeletal mechanisms that drive allocation are well characterized. We disrupt allocation by two different means: Loss of the LD protein Jabba results in LDs adhering inappropriately to glycogen granules; loss of Klar alters the activities of the microtubule motors that move LDs. Both mutants cause the same dramatic change in LD tissue inheritance, shifting allocation of the majority of LDs to the yolk cell instead of the incipient epithelium. Embryos with such mislocalized LDs do not fully consume their LDs and are delayed in hatching. Through use of a dPLIN2 mutant, which appropriately localizes a smaller pool of LDs, we find that failed LD transport and a smaller LD pool affect embryogenesis in a similar manner. Embryos of all three mutants display overlapping changes in their transcriptome and proteome, suggesting that lipid deprivation results in a shared embryonic response and a widespread change in metabolism. Excitingly, we find abundant changes related to redox homeostasis, with many proteins related to glutathione metabolism upregulated. LD deprived embryos have an increase in peroxidized lipids and rely on increased utilization of glutathione-related proteins for survival. Thus, embryos are apparently able to mount a beneficial response upon lipid stress, rewiring their metabolism to survive. In summary, we demonstrate that early embryos allocate LDs into specific lineages for subsequent optimal utilization, thus protecting against oxidative stress and ensuring punctual development.

The Tup1-Cyc8 complex in Saccharomyces cerevisiae was one of the first global co-repressors of gene transcription discovered. However, despite years of study, a full understanding of the contribution of Tup1p and Cyc8p to complex function is lacking. We examined TUP1 and CYC8 single and double deletion mutants and show that CYC8 represses more genes than TUP1, and that there are genes subject to (i) unique repression by TUP1 or CYC8, (ii) redundant repression by TUP1 and CYC8, and (iii) there are genes at which de-repression in a cyc8 mutant is dependent upon TUP1, and vice-versa. We also reveal that Tup1p and Cyc8p can make distinct contributions to commonly repressed genes most likely via specific interactions with different histone deacetylases. Furthermore, we show that Tup1p and Cyc8p can be found independently of each other to negatively regulate gene transcription and can persist at active genes to negatively regulate on-going transcription. Together, these data suggest that Tup1p and Cyc8p can associate with active and inactive genes to mediate distinct negative and positive regulatory roles when functioning within, and possibly out with the complex.

Aberrantly up-regulated activity of the type II transmembrane protease Matriptase-1 has been associated with the development and progression of a range of epithelial-derived carcinomas, and a variety of signaling pathways can mediate Matriptase-dependent tumorigenic events. During mammalian carcinogenesis, gain of Matriptase activity often results from imbalanced ratios between Matriptase and its cognate transmembrane inhibitor Hai1. Similarly, in zebrafish, unrestrained Matriptase activity due to loss of hai1a results in epidermal pre-neoplasms already during embryogenesis. Here, based on our former findings of a similar tumor-suppressive role for the Na+/K+-pump beta subunit ATP1b1a, we identify epithelial polarity defects and systemic hypotonic stress as another mode of aberrant Matriptase activation in the embryonic zebrafish epidermis in vivo. In this case, however, a different oncogenic pathway is activated which contains PI3K, AKT and NFkB, rather than EGFR and PLD (as in hai1a mutants). Strikingly, epidermal pre-neoplasm is only induced when epithelial polarity defects in keratinocytes (leading to disturbed Matriptase subcellular localization) occur in combination with systemic hypotonic stress (leading to increased proteolytic activity of Matriptase). A similar combinatorial effect of hypotonicity and loss of epithelial polarity was also obtained for the activity levels of Matriptase-1 in human MCF-10A epithelial breast cells. Together, this is in line with the multi-factor concept of carcinogenesis, with the notion that such factors can even branch off from one and the same initiator (here ATP1a1b) and can converge again at the level of one and the same mediator (here Matriptase). In sum, our data point to tonicity and epithelial cell polarity as evolutionarily conserved regulators of Matriptase activity that upon de-regulation can constitute an alternative mode of Matriptase-dependent carcinogenesis in vivo.

Protein kinases are central components of almost all signaling pathways that control cellular activities. In the model organism Saccharomyces cerevisiae, the paralogous protein kinases Ypk1 and Ypk2, which control membrane lipid homeostasis, are essential for viability, and previous studies strongly indicated that this is also the case for their single ortholog Ypk1 in the pathogenic yeast Candida albicans. Here, using FLP-mediated inducible gene deletion, we reveal that C. albicans ypk1Δ mutants are viable but slow-growing, explaining prior failures to obtain null mutants. Phenotypic analyses of the mutants showed that the functions of Ypk1 in regulating sphingolipid biosynthesis and cell membrane lipid asymmetry are conserved, but the consequences of YPK1 deletion are milder than in S. cerevisiae. Mutational studies demonstrated that the highly conserved PDK1 phosphorylation site T548 in its activation loop is essential for Ypk1 function, whereas the TORC2 phosphorylation sites S687 and T705 at the C-terminus are important for Ypk1-dependent resistance to membrane stress. Unexpectedly, Pkh1, the single C. albicans orthologue of Pkh1/Pkh2, which mediate Ypk1 phosphorylation at the PDK1 site in S. cerevisiae, was not required for normal growth of C. albicans under nonstressed conditions, and Ypk1 phosphorylation at T548 was only slightly reduced in pkh1Δ mutants. We found that another protein kinase, Pkh3, whose ortholog in S. cerevisiae cannot substitute Pkh1/2, acts redundantly with Pkh1 to activate Ypk1 in C. albicans. No phenotypic effects were observed in cells lacking Pkh3 alone, but pkh1Δ pkh3Δ double mutants had a severe growth defect and Ypk1 phosphorylation at T548 was completely abolished. These results establish that Ypk1 is not essential for viability in C. albicans and that, despite its generally conserved function, the Ypk1 signaling pathway is rewired in this pathogenic yeast and includes a novel upstream kinase to activate Ypk1 by phosphorylation at the PDK1 site.

Establishment of a proper DNA methylation landscape in mammalian oocytes is important for maternal imprinting and embryonic development. De novo DNA methylation in oocytes is mediated by the DNA methyltransferase DNMT3A, which has an ATRX-DNMT3-DNMT3L (ADD) domain that interacts with histone H3 tail unmethylated at lysine-4 (H3K4me0). The domain normally blocks the methyltransferase domain via intramolecular interaction and binding to histone H3K4me0 releases the autoinhibition. However, H3K4me0 is widespread in chromatin and the role of the ADD-histone interaction has not been studied in vivo. We herein show that amino-acid substitutions in the ADD domain of mouse DNMT3A cause dwarfism. Oocytes derived from homozygous females show mosaic loss of CG methylation and almost complete loss of non-CG methylation. Embryos derived from such oocytes die in mid-to-late gestation, with stochastic and often all-or-none-type CG-methylation loss at imprinting control regions and misexpression of the linked genes. The stochastic loss is a two-step process, with loss occurring in cleavage-stage embryos and regaining occurring after implantation. These results highlight an important role for the ADD domain in efficient, and likely processive, de novo CG methylation and pose a model for stochastic inheritance of epigenetic perturbations in germ cells to the next generation.