Search

This article discusses the concept of pangenomics and its relevance to the study of human fungal pathogens. Pangenomes are the total collection of genes within a given phylogenetic group, consisting of a core genome and an accessory genome. The accessory genome is hypothesized to be involved in fungal adaptation to the host or environment, as well as in communication, pathogenicity, and antifungal resistance. Pangenomes can help analyze genomic differences between strains of the same species and provide a reference for experimental studies. They also reduce bias introduced by single reference genomes and elucidate phenotypic variation within a species. The article highlights the importance of functional annotation and experimental characterization of accessory genes, as well as the challenges of genomic data availability and quality, and the lack of a methodological "gold standard" and data curation. Overall, fungal pangenomes have the potential to reveal novel insights into the evolution, pathogenesis, and phenotypic heterogeneity of human fungal pathogens. [Extracted from the article]

Homologous recombination is a key process that governs the stability of eukaryotic genomes during DNA replication and repair. Multiple auxiliary factors regulate the choice of homologous recombination pathway in response to different types of replication stress. Using Schizosaccharomyces pombe we have previously suggested the role of DNA translocases Rrp1 and Rrp2, together with Srs2 helicase, in the common synthesis-dependent strand annealing sub-pathway of homologous recombination. Here we show that all three proteins are important for completion of replication after hydroxyurea exposure and provide data comparing the effect of overproduction of Srs2 with Rrp1 and Rrp2. We demonstrate that Srs2 localises to rDNA region and is required for proper replication of rDNA arrays. Upregulation of Srs2 protein levels leads to enhanced replication stress, chromosome instability and viability loss, as previously reported for Rrp1 and Rrp2. Interestingly, our data suggests that dysregulation of Srs2, Rrp1 and Rrp2 protein levels differentially affects checkpoint response: overproduction of Srs2 activates simultaneously DNA damage and replication stress response checkpoints, while cells overproducing Rrp1 mainly launch DNA damage checkpoint. On the other hand, upregulation of Rrp2 primarily leads to replication stress response checkpoint activation. Overall, we propose that Srs2, Rrp1 and Rrp2 have important and at least partially independent functions in the maintenance of distinct difficult to replicate regions of the genome. [ABSTRACT FROM AUTHOR]

Exceptions to Mendelian inheritance often highlight novel chromosomal behaviors. The maize Pl1-Rhoades allele conferring plant pigmentation can display inheritance patterns deviating from Mendelian expectations in a behavior known as paramutation. However, the chromosome features mediating such exceptions remain unknown. Here we show that small RNA production reflecting RNA polymerase IV function within a distal downstream set of five tandem repeats is coincident with meiotically-heritable repression of the Pl1-Rhoades transcription unit. A related pl1 haplotype with three, but not one with two, repeat units also displays the trans-homolog silencing typifying paramutations. 4C interactions, CHD3a-dependent small RNA profiles, nuclease sensitivity, and polyadenylated RNA levels highlight a repeat subregion having regulatory potential. Our comparative and mutant analyses show that transcriptional repression of Pl1-Rhoades correlates with 24-nucleotide RNA production and cytosine methylation at this subregion indicating the action of a specific DNA-dependent RNA polymerase complex. These findings support a working model in which pl1 paramutation depends on trans-chromosomal RNA-directed DNA methylation operating at a discrete cis-linked and copy-number-dependent transcriptional regulatory element. Author summary: Although switching genes "on" and "off" in a coordinated fashion is essential to organismal development and homeostasis, some genes ‐ or specific forms thereof ‐ display unusual switching behaviors where the "off" state is sexually transmitted to offspring. Further, only "off" states persist when both "on" and "off" states are combined in heterozygous condition. Examples of this unusual behavior, known as paramutation, exist for certain transgenes and/or endogenous genes in both plants and animals but the chromosome features triggering this switching and the responsible cellular machinery are poorly understood. Here, we used long-read DNA sequencing, comparative genomics, DNA -DNA interaction assays, mutant analyses, and cytosine methylation profiling to identify a switchable regulatory sequence associated with paramutation at the maize purple plant1 (pl1) locus. These key sequences are embedded within a 5-copy tandem repeat array of mostly transposable elements ~14 kilobases distal to the pl1 gene. The "off" state is associated with these sequences producing so-called "small interfering RNAs" (siRNAs) that may direct the recruitment of chromatin-modifying enzymes. Our findings support a molecular model for paramutation and highlight the regulatory importance of repeated sequences found in large eukaryotic genomes. [ABSTRACT FROM AUTHOR]

DNA methylation is a key regulator of eukaryote genomes, and is of particular relevance in the regulation of gene expression on the sex chromosomes, with a key role in dosage compensation in mammalian XY systems. In the case of birds, dosage compensation is largely absent, with it being restricted to two small Male Hyper-Methylated (MHM) regions on the Z chromosome. To investigate how variation in DNA methylation is regulated on the Z chromosome we utilised a wild x domestic advanced intercross in the chicken, with both hypothalamic methylomes and transcriptomes assayed in 124 individuals. The relatively large numbers of individuals allowed us to identify additional genomic MHM regions on the Z chromosome that were significantly differentially methylated between the sexes. These regions appear to down-regulate local gene expression in males, but not remove it entirely (unlike the lncRNAs identified in the initial MHM regions). These MHM regions were further tested and the most balanced genes appear to show decreased expression in males, whilst methylation appeared to be far more correlated with gene expression in the less balanced, as compared to the most balanced genes. In addition, trans effect hotspots were also identified that were based on the autosomes but affected the Z, and also that were based on the Z chromosome but that affected autosomal DNA methylation regulation. In addition, quantitative trait loci (QTL) that regulate variation in methylation on the Z chromosome, and those loci that regulate methylation on the autosomes that derive from the Z chromosome were mapped. Trans-effect hotspots were also identified that were based on the autosomes but affected the Z, and also one that was based on the Z chromosome but that affected both autosomal and sex chromosome DNA methylation regulation. We show that both cis and trans loci that originate from the Z chromosome never exhibit an interaction with sex, whereas trans loci originating from the autosomes but affecting the Z chromosome always display such an interaction. Our results highlight how additional MHM regions are actually present on the Z chromosome, and they appear to have smaller-scale effects on gene expression in males. Quantitative variation in methylation is also regulated both from the autosomes to the Z chromosome, and from the Z chromosome to the autosomes. Author summary: DNA methylation is a key regulator of eukaryote genomes, and is of particular relevance in the regulation of gene expression on the sex chromosomes, with a key role in dosage compensation in mammalian XY systems. In the case of birds, dosage compensation is largely absent, with it being restricted to two small Male Hyper-Methylated (MHM) regions on the Z chromosome. We utilised a wild x domestic advanced intercross in the chicken, with both hypothalamic methylomes and transcriptomes assayed in 124 individuals, to investigate the role that methylation plays in regulating gene expression on the Z chromosome. Our results highlight how additional MHM regions are actually present on the Z chromosome, and they appear to have smaller-scale effects on gene expression in males. Quantitative variation in methylation is also regulated both from the autosomes to the Z chromosome, and from the Z chromosome to the autosomes. In addition, these MHM regions were further tested and the most balanced genes appear to show decreased expression in males, whilst methylation appeared to be far more correlated with gene expression in the less balanced, as compared to the most balanced genes. [ABSTRACT FROM AUTHOR]

Recent advances in long-read sequencing technologies have enabled the complete assembly of eukaryotic genomes from telomere to telomere by allowing repeated regions to be fully sequenced and assembled, thus filling the gaps left by previous short-read sequencing methods. Furthermore, long-read sequencing can also help characterizing structural variants, with applications in the fields of genome evolution or cancer genomics. For many organisms, the main bottleneck to sequence long reads remains the lack of robust methods to obtain high-molecular-weight (HMW) DNA. For this purpose, we developed an optimized protocol to extract DNA suitable for long-read sequencing from the unicellular green alga Chlamydomonas reinhardtii, based on CTAB/phenol extraction followed by a size selection step for long DNA molecules. We provide validation results for the extraction protocol, as well as statistics obtained with Oxford Nanopore Technologies sequencing. [ABSTRACT FROM AUTHOR]

Analysis of eukaryotic genomes requires the detection and classification of transposable elements (TEs), a crucial but complex and time-consuming task. To improve the performance of tools that accomplish these tasks, Machine Learning approaches (ML) that leverage computer resources, such as GPUs (Graphical Processing Unit) and multiple CPU (Central Processing Unit) cores, have been adopted. However, until now, the use of ML techniques has mostly been limited to classification of TEs. Herein, a detection-classification strategy (named YORO) based on convolutional neural networks is adapted from computer vision (YOLO) to genomics. This approach enables the detection of genomic objects through the prediction of the position, length, and classification in large DNA sequences such as fully sequenced genomes. As a proof of concept, the internal protein-coding domains of LTR-retrotransposons are used to train the proposed neural network. Precision, recall, accuracy, F1-score, execution times and time ratios, as well as several graphical representations were used as metrics to measure performance. These promising results open the door for a new generation of Deep Learning tools for genomics. YORO architecture is available at https://github.com/simonorozcoarias/YORO. [ABSTRACT FROM AUTHOR]

Pervasive transcription of eukaryotic genomes generates non-coding transcripts with regulatory potential. We examined the effects of non-coding antisense transcription on the regulation of expression of the yeast PHO5 gene, a paradigmatic case for gene regulation through promoter chromatin remodeling. A negative role for antisense transcription at the PHO5 gene locus was demonstrated by leveraging the level of overlapping antisense transcription through specific mutant backgrounds, expression from a strong promoter in cis, and use of the CRISPRi system. Furthermore, we showed that enhanced elongation of PHO5 antisense leads to a more repressive chromatin conformation at the PHO5 gene promoter, which is more slowly remodeled upon gene induction. The negative effect of antisense transcription on PHO5 gene transcription is mitigated upon inactivation of the histone deacetylase Rpd3, showing that PHO5 antisense RNA acts via histone deacetylation. This regulatory pathway leads to Rpd3-dependent decreased recruitment of the RSC chromatin remodeling complex to the PHO5 gene promoter upon induction of antisense transcription. Overall, the data in this work reveal an additional level in the complex regulatory mechanism of PHO5 gene expression by showing antisense transcription-mediated repression at the level of promoter chromatin structure remodeling. Author summary: Non-coding transcripts synthesized by RNA Poll are short-lived because they are rapidly terminated by the NNS complex and degraded by the RNA exosome. However, some non-coding transcripts can regulate the expression of coding genes at whose loci they are initiated. In this work, we show that a non-coding antisense transcript regulates the expression of the yeast S. cerevisiae PHO5 gene that encodes a periplasmic acid phosphatase. Furthermore, we show that the repressive effect of antisense transcription at the PHO5 gene locus is mediated by chromatin-modifying complexes, such as the histone deacetylase Rrp3 and the chromatin-remodeling complex RSC. Consequently, elongation of the antisense transcript through the PHO5 promoter region results in a chromatin conformation that is more difficult to remodel upon gene induction and therefore in slower kinetics of gene expression. Overall, this work paints a more complete picture of the regulatory process involved in repression of the yeast PHO5 promoter, a well-studied model for gene regulation through chromatin structure remodeling. [ABSTRACT FROM AUTHOR]

Chromatin insulators are responsible for orchestrating long-range interactions between enhancers and promoters throughout the genome and align with the boundaries of Topologically Associating Domains (TADs). Here, we demonstrate an association between gypsy insulator proteins and the phosphorylated histone variant H2Av (γH2Av), normally a marker of DNA double strand breaks. Gypsy insulator components colocalize with γH2Av throughout the genome, in polytene chromosomes and in diploid cells in which Chromatin IP data shows it is enriched at TAD boundaries. Mutation of insulator components su(Hw) and Cp190 results in a significant reduction in γH2Av levels in chromatin and phosphatase inhibition strengthens the association between insulator components and γH2Av and rescues γH2Av localization in insulator mutants. We also show that γH2Av, but not H2Av, is a component of insulator bodies, which are protein condensates that form during osmotic stress. Phosphatase activity is required for insulator body dissolution after stress recovery. Together, our results implicate the H2A variant with a novel mechanism of insulator function and boundary formation. Author summary: The DNA in eukaryotic genomes is folded into domains called Topologically Associating Domains (TADS), which promote gene specific transcription regulation. Insulator proteins are DNA binding proteins that bind at the boundaries between adjacent TADs. Loop extrusion is a mechanism by which insulators promote TAD formation in vertebrate genomes, but the mechanism by which insulator proteins facilitate the formation of boundaries in Drosophila is not well understood. In this wok we show that there is an association between Drosophila Gypsy insulator proteins and the phosphorylated version of the histone variant H2Av (γH2Av). γH2Av has been traditionally linked to the mechanism of DNA repair, but our data shows that Gypsy insulator components colocalize with γH2Av throughout the genome, and that γH2Av is also enriched at TAD boundaries. Mutation of genes encoding insulator proteins su(Hw) and Cp190 results in a significant reduction in γH2Av levels, and inhibition of the phosphatase activity that removes phosphate from γH2Av strengthens the association between insulator proteins and γH2Av. We also show that γH2Av is a component of insulator bodies, which are protein condensates that form during osmotic stress. Together, our results implicate the H2A variant with a novel mechanism of insulator function and boundary formation. [ABSTRACT FROM AUTHOR]

Rapid advancements in automated genomic technologies have uncovered many unique findings about the turtle genome and its associated features including olfactory gene expansions and duplications of toll-like receptors. However, despite the advent of large-scale sequencing, assembly, and annotation, about 40–50% of genes in eukaryotic genomes are left without functional annotation, severely limiting our knowledge of the biological information of genes. Additionally, these automated processes are prone to errors since draft genomes consist of several disconnected scaffolds whose order is unknown; erroneous draft assemblies may also be contaminated with foreign sequences and propagate to cause errors in annotation. Many of these automated annotations are thus incomplete and inaccurate, highlighting the need for functional annotation to link gene sequences to biological identity. In this study, we have functionally annotated two genes of the red-bellied short-neck turtle (Emydura subglobosa), a member of the relatively understudied pleurodire lineage of turtles. We improved upon initial ab initio gene predictions through homology-based evidence and generated refined consensus gene models. Through functional, localization, and structural analyses of the predicted proteins, we discovered conserved putative genes encoding mitochondrial proteins that play a role in C21-steroid hormone biosynthetic processes and fatty acid catabolism—both of which are distantly related by the tricarboxylic acid (TCA) cycle and share similar metabolic pathways. Overall, these findings further our knowledge about the genetic features underlying turtle physiology, morphology, and longevity, which have important implications for the treatment of human diseases and evolutionary studies. [ABSTRACT FROM AUTHOR]

Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps. Author summary: Nucleosomes are the conservative building blocks of the chromatin, and their array regularity and positioning level correlate with transcription activity, but their underlying distributions in eukaryote genomes have not been comprehensively studied due to the poorly positioned ones. Recently available high-throughput enhanced chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO provide information of nucleosome-level chromatin proximity, including the positions of mononucleosomes and the spacing between nucleosome pairs, enabling identifying nucleosomes in poorly positioned regions. Here, we present NucleoMap, a nucleosome identification approach from ultra-high resolution chromatin contact maps. In this paper, we provide an overview of NucleoMap's workflow and capabilities. We benchmark NucleoMap with popular baseline methods in multiple datasets. Next, we rigorously characterize genome-wide association between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions in eukaryotes. Based on the identified nucleosomes, we also construct more precise and more interpretable nucleosome contact maps, which preserve the inter-nucleosome distances. [ABSTRACT FROM AUTHOR]

Eukaryotic genomes often encode multiple transporters for the same nutrient. For example, budding yeast has 17 hexose transporters (HXTs), all of which potentially transport glucose. Using mathematical modelling, we show that transporters that use either facilitated diffusion or symport can have a rate-affinity tradeoff, where an increase in the maximal rate of transport decreases the transporter's apparent affinity. These changes affect the import flux non-monotonically, and for a given concentration of extracellular nutrient there is one transporter, characterised by its affinity, that has a higher import flux than any other. Through encoding multiple transporters, cells can therefore mitigate the tradeoff by expressing those transporters with higher affinities in lower concentrations of nutrients. We verify our predictions using fluorescent tagging of seven HXT genes in budding yeast and follow their expression over time in batch culture. Using the known affinities of the corresponding transporters, we show that their regulation in glucose is broadly consistent with a rate-affinity tradeoff: as glucose falls, the levels of the different transporters peak in an order that mostly follows their affinity for glucose. More generally, evolution is constrained by tradeoffs. Our findings indicate that one such tradeoff often occurs in the cellular transport of nutrients. Author summary: From yeast to humans, cells often express multiple different types of transporters for the same nutrient, and it is puzzling why a single high-affinity transporter is not expressed instead. Here we initially use mathematical modelling to demonstrate that transporters facilitating diffusion and those powered by the proton motive force can both exhibit a rate-affinity tradeoff, for quite general conditions. A transporter with a higher affinity necessarily has a lower rate, and vice versa. The tradeoff implies that there is a range of nutrient concentrations for which a transporter, characterised by its affinity, has a higher import flux than any other transporter with a different affinity. To mitigate the tradeoff, genomes may therefore encode multiple different transporters, and cells that express each transporter in the concentrations where it imports best will uptake nutrients at higher rates. Consistently, we show that as cells of budding yeast consume glucose, they express five types of hexose transporters in an order that follows the transporters' affinities. [ABSTRACT FROM AUTHOR]

Phylogenetic profiling in eukaryotes is of continued interest to study and predict the functional relationships between proteins. This interest is likely driven by the increased number of available diverse genomes and computational methods to infer orthologies. The evaluation of phylogenetic profiles has mainly focussed on reference genome selection in prokaryotes. However, it has been proven to be challenging to obtain high prediction accuracies in eukaryotes. As part of our recent comparison of orthology inference methods for eukaryotic genomes, we observed a surprisingly high performance for predicting interacting orthologous groups. This high performance, in turn, prompted the question of what factors influence the success of phylogenetic profiling when applied to eukaryotic genomes. Here we analyse the effect of species, orthologous group and interactome selection on protein interaction prediction using phylogenetic profiles. We select species based on the diversity and quality of the genomes and compare this supervised selection with randomly generated genome subsets. We also analyse the effect on the performance of orthologous groups defined to be in the last eukaryotic common ancestor of eukaryotes to that of orthologous groups that are not. Finally, we consider the effects of reference interactome set filtering and reference interactome species. In agreement with other studies, we find an effect of genome selection based on quality, less of an effect based on genome diversity, but a more notable effect based on the amount of information contained within the genomes. Most importantly, we find it is not merely selecting the correct genomes that is important for high prediction performance. Other choices in meta parameters such as orthologous group selection, the reference species of the interaction set, and the quality of the interaction set have a much larger impact on the performance when predicting protein interactions using phylogenetic profiles. These findings shed light on the differences in reported performance amongst phylogenetic profiles approaches, and reveal on a more fundamental level for which types of protein interactions this method has most promise when applied to eukaryotes. [ABSTRACT FROM AUTHOR]

Pangenome analyses are increasingly being utilized to study the evolution of eukaryotic organisms. While pangenomes can provide insight into polymorphic gene content, inferences about the ecological and adaptive potential of such organisms also need to be accompanied by additional supportive genomic analyses. In this study we constructed a pangenome of Claviceps purpurea from 24 genomes and examined the positive selection and recombination landscape of an economically important fungal organism for pharmacology and agricultural research. Together, these analyses revealed that C. purpurea has a relatively large accessory genome (~ 38%), high recombination rates (ρ = 0.044), and transposon mediated gene duplication. However, due to observations of relatively low transposable element (TE) content (8.8%) and a lack of variability in genome sizes, prolific TE expansion may be controlled by frequent recombination. We additionally identified that within the ergoline biosynthetic cluster the lpsA1 and lpsA2 were the result of a recombination event. However, the high recombination rates observed in C. purpurea may be influencing an overall trend of purifying selection across the genome. These results showcase the use of selection and recombination landscapes to identify mechanisms contributing to pangenome structure and primary factors influencing the evolution of an organism. [ABSTRACT FROM AUTHOR]

The cellular energy metabolism is one of the most conserved processes, as it is present in all living organisms. Mitochondria are providing the eukaryotic cell with energy and thus their genome and gene expression has been of broad interest for a long time. Mitochondrial gene expression changes under different conditions and is regulated by genes encoded in the nucleus of the cell. In this context, little is known about non-model organisms and we provide the first large-scaled gene expression analysis of mitochondrial-linked genes in laying hens. We analysed 28 mitochondrial and nuclear genes in 100 individuals in the context of five life-stages and strain differences among five tissues. Our study showed that mitochondrial gene expression increases during the productive life span, and reacts tissue and strain specific. In addition, the strains react different to potential increased oxidative stress, resulting from the increase in mitochondrial gene expression. The results suggest that the cellular energy metabolism as part of a complex regulatory system is strongly affected by the productive life span in laying hens and thus partly comparable to model organisms. This study provides a starting point for further analyses in this field on non-model organisms, especially in laying-hens. [ABSTRACT FROM AUTHOR]

Intergenic transcription is a common feature of eukaryotic genomes and performs important and diverse cellular functions. Here, we investigate the iab-8 ncRNA from the Drosophila Bithorax Complex and show that this RNA is able to repress the transcription of genes located at its 3' end by a sequence-independent, transcriptional interference mechanism. Although this RNA is expressed in the early epidermis and CNS, we find that its repressive activity is limited to the CNS, where, in wild-type embryos, it acts on the Hox gene, abd-A, located immediately downstream of it. The CNS specificity is achieved through a 3' extension of the transcript, mediated by the neuronal-specific, RNA-binding protein, ELAV. Loss of ELAV activity eliminates the 3' extension and results in the ectopic activation of abd-A. Thus, a tissue-specific change in the length of a ncRNA is used to generate a precise pattern of gene expression in a higher eukaryote. Author summary: Although all of the cells making up complex organisms contain the same genetic material, they are nevertheless able to create the diverse tissues of the body. They do this by changing the genes they express. Thus, understanding how genes are controlled in a tissue-specific fashion is one of the primary interests of molecular genetics. Within the bithorax homeotic complex of the fruit fly Drosophila melanogaster, we, and others, previously showed that a >92 kb-long non-coding RNA, called the iab-8 ncRNA, downregulates many important developmental genes, including its genomic downstream neighbor, the homeotic gene abd-A. This downregulation is important as its loss is linked to female sterility. Interestingly, we find that the iab-8 ncRNA regulates abd-A through a mechanism called transcriptional interference, where one gene downregulates a target gene by transcribing over it. In the case of iab-8, this process is limited to the posterior central nervous system, where the iab-8 ncRNA is specifically extended into the abd-A gene by the action of the neuronal-specific RNA binding protein, ELAV. Overall, our work highlights a largely unexplored mechanism by which tissue-specific gene regulation is achieved. [ABSTRACT FROM AUTHOR]

The cost of sequencing the genome is dropping at a much faster rate compared to assembling and finishing the genome. The use of lightly sampled genomes (genome-skims) could be transformative for genomic ecology, and results using k-mers have shown the advantage of this approach in identification and phylogenetic placement of eukaryotic species. Here, we revisit the basic question of estimating genomic parameters such as genome length, coverage, and repeat structure, focusing specifically on estimating the k-mer repeat spectrum. We show using a mix of theoretical and empirical analysis that there are fundamental limitations to estimating the k-mer spectra due to ill-conditioned systems, and that has implications for other genomic parameters. We get around this problem using a novel constrained optimization approach (Spline Linear Programming), where the constraints are learned empirically. On reads simulated at 1X coverage from 66 genomes, our method, REPeat SPECTra Estimation (RESPECT), had < 1.5% error in length estimation compared to 34% error previously achieved. In shotgun sequenced read samples with contaminants, RESPECT length estimates had median error 4%, in contrast to other methods that had median error 80%. Together, the results suggest that low-pass genomic sequencing can yield reliable estimates of the length and repeat content of the genome. The RESPECT software will be publicly available at https://urldefense.proofpoint.com/v2/url?u=https-3A__github.com_shahab-2Dsarmashghi_RESPECT.git&d=DwIGAw&c=-35OiAkTchMrZOngvJPOeA&r=ZozViWvD1E8PorCkfwYKYQMVKFoEcqLFm4Tg49XnPcA&m=f-xS8GMHKckknkc7Xpp8FJYw_ltUwz5frOw1a5pJ81EpdTOK8xhbYmrN4ZxniM96&s=717o8hLR1JmHFpRPSWG6xdUQTikyUjicjkipjFsKG4w&e=. Author summary: The cost of sequencing the genome is dropping at a much faster rate compared to assembling and finishing the genome. The use of lightly sampled genomes (genome skims) could be transformative for genomic ecology. Analyzing genome skims, mostly based on statistics of small oligomers, remains challenging, but recent results have shown the advantage of this approach for the identification and phylogenetic placement of eukaryotic species. In this paper, we present a method, RESPECT, to estimate genomic properties such as genome length and repetitiveness from low-coverage genome skims. We trained RESPECT using assembled genomes and tested it on low-coverage simulated and real reads. Benchmarking results reveal that RESPECT has excellent accuracy in estimating the genome length compared to other methods, and can provide critical information regarding the repeat structure of the genome. [ABSTRACT FROM AUTHOR]

Beyond their canonical function in nucleocytoplasmic exchanges, nuclear pore complexes (NPCs) regulate the expression of protein-coding genes. Here, we have implemented transcriptomic and molecular methods to specifically address the impact of the NPC on retroelements, which are present in multiple copies in genomes. We report a novel function for the Nup84 complex, a core NPC building block, in specifically restricting the transcription of LTR-retrotransposons in yeast. Nup84 complex-dependent repression impacts both Copia and Gypsy Ty LTR-retrotransposons, all over the S. cerevisiae genome. Mechanistically, the Nup84 complex restricts the transcription of Ty1, the most active yeast retrotransposon, through the tethering of the SUMO-deconjugating enzyme Ulp1 to NPCs. Strikingly, the modest accumulation of Ty1 RNAs caused by Nup84 complex loss-of-function is sufficient to trigger an important increase of Ty1 cDNA levels, resulting in massive Ty1 retrotransposition. Altogether, our study expands our understanding of the complex interactions between retrotransposons and the NPC, and highlights the importance for the cells to keep retrotransposon under tight transcriptional control. Author summary: Retroelements, which replicate by reverse transcription of their RNA into a cDNA that is integrated into the host genetic material, play an important role in the plasticity of eukaryotic genomes. The study of yeast retrotransposons has led to the identification of host factors that limit or facilitate retroelement mobility, including components of the nuclear pore complex (NPC), most of them still awaiting mechanistic characterization. Here, we investigated the contribution of the Nup84 complex, a core NPC scaffold, to retrotransposon biology in budding yeast. Our findings uncover that the Nup84 complex restricts the transcription of phylogenetically-distinct Ty retroelements. By focusing on Ty1 retrotransposons, we provide evidence that repression by the Nup84 complex depends on the maintenance at the NPC of the SUMO-protease Ulp1, an essential enzyme of the SUMO pathway with multiple targets in the transcription machinery. We finally show that this transcriptional control is critical for genome dynamics, since a small increase in the accumulation of Ty1 RNAs leads to massive retrotransposition. Our data reveal that although relatively abundant, Ty transcripts are limiting for retrotransposition, underscoring the importance of having their expression tightly controlled. We also characterize a new non-canonical function of NPCs, confirming their connection with genome expression and stability. [ABSTRACT FROM AUTHOR]

Transposable elements (TEs) represent a major portion of most eukaryotic genomes, yet little is known about their mutation rates or how their activity is shaped by other evolutionary forces. Here, we compare short- and long-term patterns of genome-wide mutation accumulation (MA) of TEs among 9 genotypes from three populations of Daphnia magna from across a latitudinal gradient. While the overall proportion of the genome comprised of TEs is highly similar among genotypes from Finland, Germany, and Israel, populations are distinguishable based on patterns of insertion site polymorphism. Our direct rate estimates indicate TE movement is highly variable (net rates ranging from -11.98 to 12.79 x 10−5 per copy per generation among genotypes), differing both among populations and TE families. Although gains outnumber losses when selection is minimized, both types of events appear to be highly deleterious based on their low frequency in control lines where propagation is not limited to random, single-progeny descent. With rate estimates 4 orders of magnitude higher than base substitutions, TEs clearly represent a highly mutagenic force in the genome. Quantifying patterns of intra- and interspecific variation in TE mobility with and without selection provides insight into a powerful mechanism generating genetic variation in the genome. Author summary: Transposable elements (TEs) are a significant portion of most eukaryotic genomes, yet our understanding of their rates of mobility and their patterns of accumulation remain very limited. Here, we estimate genome-wide rates of gain and loss of TEs in Daphnia magna, a well-studied model organism in ecology, and compare these rates of mutation to the long-term accumulation of TEs in the genome. Rates vary remarkably among genotypes and populations, and between different types of TEs within the same lineage. Despite this variation, over long time periods, TE content in the genome is extremely similar across genotypes within the species. We compare our results to the few estimates available from other taxa, and argue that TEs are an important source of mutagenesis in the genome worthy of further investigation. [ABSTRACT FROM AUTHOR]

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity. Author summary: Transposable elements (TEs) are pieces of DNA that can jump to new positions in the genome. When they land at a new location, they generate a mutation. Such mutations in genes affecting kernel and plant pigmentation allowed the discovery of TEs in maize in the 1940's. Since, we have learned that TEs are a ubiquitous feature of eukaryotic genomes, and that TEs make up over 85% of all the DNA in a maize genome. Here, we investigate the roles of individual TE copies and TE species interacting within the maize genome, and how these relationships are analogous to ecological communities. The community of transposable elements within a single genome represents a rich diversity of ecological strategies for survival in the complex, hostile environment of a host genome. [ABSTRACT FROM AUTHOR]

Regulatory elements control gene expression through transcription initiation (promoters) and by enhancing transcription at distant regions (enhancers). Accurate identification of regulatory elements is fundamental for annotating genomes and understanding gene expression patterns. While there are many attempts to develop computational promoter and enhancer identification methods, reliable tools to analyze long genomic sequences are still lacking. Prediction methods often perform poorly on the genome-wide scale because the number of negatives is much higher than that in the training sets. To address this issue, we propose a dynamic negative set updating scheme with a two-model approach, using one model for scanning the genome and the other one for testing candidate positions. The developed method achieves good genome-level performance and maintains robust performance when applied to other vertebrate species, without re-training. Moreover, the unannotated predicted regulatory regions made on the human genome are enriched for disease-associated variants, suggesting them to be potentially true regulatory elements rather than false positives. We validated high scoring "false positive" predictions using reporter assay and all tested candidates were successfully validated, demonstrating the ability of our method to discover novel human regulatory regions. Author summary: Identification of regulatory elements (promoters and enhancers) is important for understanding gene expression patterns. The set of promoters and enhancers is not complete for non-model organisms and even for the human genome there are still unannotated regions, such as alternative promoters for the known genes or promoters that are only expressed in a small fraction of cells or under specific conditions. Despite the development of experimental techniques, the regulatory regions annotation remains expensive and laborious and computational methods can speed up this process by providing candidates for the validation. We developed an easy-to-use tool capable of regulatory regions annotation in eukaryotic genomes. The developed method reduces the number of false positives made by including difficult samples in the training set. The method consists of two deep learning models, where one model scans the genome and identifies putative regulatory regions while the other model pinpoints the Transcription Start Site (TSS) location within the identified region. The predicted regions were validated using reporter assay, finding previously unknown regulatory regions in the human genome. The trained model achieved good genome-wide performance and was supported by meaningful extracted biological features. [ABSTRACT FROM AUTHOR]

piRNAs are small non-coding RNAs that guide the silencing of transposons and other targets in animal gonads. In Drosophila female germline, many piRNA source loci dubbed "piRNA clusters" lack hallmarks of active genes and exploit an alternative path for transcription, which relies on the Rhino-Deadlock-Cutoff (RDC) complex. RDC was thought to be absent in testis, so it remains to date unknown how piRNA cluster transcription is regulated in the male germline. We found that components of RDC complex are expressed in male germ cells during early spermatogenesis, from germline stem cells (GSCs) to early spermatocytes. RDC is essential for expression of dual-strand piRNA clusters and transposon silencing in testis; however, it is dispensable for expression of Y-linked Suppressor of Stellate piRNAs and therefore Stellate silencing. Despite intact Stellate repression, males lacking RDC exhibited compromised fertility accompanied by germline DNA damage and GSC loss. Thus, piRNA-guided repression is essential for normal spermatogenesis beyond Stellate silencing. While RDC associates with multiple piRNA clusters in GSCs and early spermatogonia, its localization changes in later stages as RDC concentrates on a single X-linked locus, AT-chX. Dynamic RDC localization is paralleled by changes in piRNA cluster expression, indicating that RDC executes a fluid piRNA program during different stages of spermatogenesis. These results disprove the common belief that RDC is dispensable for piRNA biogenesis in testis and uncover the unexpected, sexually dimorphic and dynamic behavior of a core piRNA pathway machinery. Author summary: Large fractions of eukaryotic genomes are occupied by mobile genetic elements called transposons. Active transposons can move in the genome causing DNA damage and mutations, while inactive copies can contribute to chromosome organization and regulation of gene expression. Host cells employ several mechanisms to discriminate transposons from other genes and repress transposon activities. In germ cells, a conserved class of short RNAs called Piwi-interacting (pi)RNAs recognize target RNAs in both the nucleus and cytoplasm and then guide transposon repression by preventing their transcription and destroying their RNAs. piRNAs are encoded in extended genomic regions dubbed piRNA clusters. Previously, composition and regulation of piRNA clusters were studied in the female germline of fruit flies, where a nuclear protein complex, the RDC complex, was shown to promote non-canonical transcription of these regions. However, RDC was believed to be dispensable in males. Here, we showed that RDC is essential for transposon repression in males, and males lacking RDC exhibit compromised fertility and loss of germ cells. We found that RDC binds multiple piRNA clusters in early germ cells but concentrates on a single locus at later stages. Our results indicate dynamic regulation of loci that produce piRNAs and, therefore, piRNA targets throughout spermatogenesis. [ABSTRACT FROM AUTHOR]

In every kingdom of life, GC->AT transitions occur more frequently than any other type of mutation due to the spontaneous deamination of cytidine. In eukaryotic genomes, this slow loss of GC base pairs is counteracted by biased gene conversion which increases genomic GC content as part of the recombination process. However, this type of biased gene conversion has not been observed in bacterial genomes, so we hypothesized that GC->AT transitions cause a reduction of genomic GC content in prokaryotic genomes on an evolutionary time scale. To test this hypothesis, we used a phylogenetic approach to analyze triplets of closely related genomes representing a wide range of the bacterial kingdom. The resulting data indicate that genomic GC content is drifting downward in bacterial genomes where GC base pairs comprise 40% or more of the total genome. In contrast, genomes containing less than 40% GC base pairs have fewer opportunities for GC->AT transitions to occur so genomic GC content is relatively stable or actually increasing. It should be noted that this observed change in genomic GC content is the net change in shared parts of the genome and does not apply to parts of the genome that have been lost or acquired since the genomes being compared shared common ancestor. However, a more detailed analysis of two Caulobacter genomes revealed that the acquisition of mobile elements by the two genomes actually reduced the total genomic GC content as well. [ABSTRACT FROM AUTHOR]

We present a model of the evolution of control systems in a genome under environmental constraints. The model conceptually follows the Jacob and Monod model of gene control. Genes contain control elements which respond to the internal state of the cell as well as the environment to control expression of a coding region. Control and coding regions evolve to maximize a fitness function between expressed coding sequences and the environment. The model was run 118 times to an average of 1.4∙106 'generations' each with a range of starting parameters probed the conditions under which genomes evolved a 'default style' of control. Unexpectedly, the control logic that evolved was not significantly correlated to the complexity of the environment. Genetic logic was strongly correlated with genome complexity and with the fraction of genes active in the cell at any one time. More complex genomes correlated with the evolution of genetic controls in which genes were active ('default on'), and a low fraction of genes being expressed correlated with a genetic logic in which genes were biased to being inactive unless positively activated ('default off' logic). We discuss how this might relate to the evolution of the complex eukaryotic genome, which operates in a 'default off' mode. [ABSTRACT FROM AUTHOR]

The human zinc finger antiviral protein (ZAP) recognizes RNA by binding to CpG dinucleotides. Mammalian transcriptomes are CpG-poor, and ZAP may have evolved to exploit this feature to specifically target non-self viral RNA. Phylogenetic analyses reveal that ZAP and its paralogue PARP12 share an ancestral gene that arose prior to extensive eukaryote divergence, and the ZAP lineage diverged from the PARP12 lineage in tetrapods. Notably, the CpG content of modern eukaryote genomes varies widely, and ZAP-like genes arose subsequent to the emergence of CpG-suppression in vertebrates. Human PARP12 exhibited no antiviral activity against wild type and CpG-enriched HIV-1, but ZAP proteins from several tetrapods had antiviral activity when expressed in human cells. In some cases, ZAP antiviral activity required a TRIM25 protein from the same or related species, suggesting functional co-evolution of these genes. Indeed, a hypervariable sequence in the N-terminal domain of ZAP contributed to species-specific TRIM25 dependence in antiviral activity assays. Crosslinking immunoprecipitation coupled with RNA sequencing revealed that ZAP proteins from human, mouse, bat and alligator exhibit a high degree of CpG-specificity, while some avian ZAP proteins appear more promiscuous. Together, these data suggest that the CpG- rich RNA directed antiviral activity of ZAP-related proteins arose in tetrapods, subsequent to the onset of CpG suppression in certain eukaryote lineages, with subsequent species-specific adaptation of cofactor requirements and RNA target specificity. Author summary: To control viral infections, cells have evolved a variety of mechanisms that detect, modify and sometimes eliminate viral components. One of such mechanism is the Zinc Finger Antiviral Protein (ZAP) which binds RNA sequences that are rich in elements composed of a cytosine followed by a guanine. Selection of viral RNA can only be achieved because such elements are sparse in RNAs encoded by human genes. Here, we traced the molecular evolution of ZAP. We found that ZAP and a closely related gene, PARP12, originated from the same ancestral gene that existed in a predecessor of vertebrates and invertebrates. We found that ZAP proteins from mammals, birds and reptiles have antiviral activity but only in the presence of a co-factor, TRIM25, from the same species. ZAP proteins from birds were particularly interesting since they demonstrated a broader antiviral activity, primarily driven by relaxed requirement for cytosine-guanine. Our findings suggest that viruses that infect birds–which are important vectors for human diseases–are under differential selective pressures and this property may influence the outcome of interspecies transmission. [ABSTRACT FROM AUTHOR]

The overarching trend in mitochondrial genome evolution is functional streamlining coupled with gene loss; therefore, gene acquisition by mitochondria is considered to be exceedingly rare. Selfish elements in the form of self-splicing introns occur in many organellar genomes, but the wider diversity of selfish elements, and how they persist in the DNA of organelles, has not been explored. In the mitochondrial genome of a marine heterotrophic katablepharid protist, we identify a functional type II restriction modification (RM) system originating from a horizontal gene transfer (HGT) event involving bacteria related to flavobacteria. This RM system consists of an HpaII-like endonuclease and a cognate cytosine methyltransferase (CM). We demonstrate that these proteins are functional by heterologous expression in both bacterial and eukaryotic cells. These results suggest that a mitochondrial-encoded RM system can function as a toxin–antitoxin selfish element and that such elements could be co-opted by eukaryotic genomes to drive biased organellar inheritance. This study reveals that a functional type II restriction modification system of flavobacterial ancestry has been horizontally transferred into the mitochondrion of a marine protist and is capable of encoding potent function, perhaps allowing it to play a role in inter-organellar warfare or protection against further integration of foreign DNA. [ABSTRACT FROM AUTHOR]

Long regulatory elements (LREs), such as CpG islands, polydA:dT tracts or AU-rich elements, are thought to play key roles in gene regulation but, as opposed to conventional binding sites of transcription factors, few methods have been proposed to formally and automatically characterize them. We present here a computational approach named DExTER (Domain Exploration To Explain gene Regulation) dedicated to the identification of candidate LREs (cLREs) and apply it to the analysis of the genomes of P. falciparum and other eukaryotes. Our analyses show that all tested genomes contain several cLREs that are somewhat conserved along evolution, and that gene expression can be predicted with surprising accuracy on the basis of these long regions only. Regulation by cLREs exhibits very different behaviours depending on species and conditions. In P. falciparum and other Apicomplexan organisms as well as in Dictyostelium discoideum, the process appears highly dynamic, with different cLREs involved at different phases of the life cycle. For multicellular organisms, the same cLREs are involved in all tissues, but a dynamic behavior is observed along embryonic development stages. In P. falciparum, whose genome is known to be strongly depleted of transcription factors, cLREs are predictive of expression with an accuracy above 70%, and our analyses show that they are associated with both transcriptional and post-transcriptional regulation signals. Moreover, we assessed the biological relevance of one LRE discovered by DExTER in P. falciparum using an in vivo reporter assay. The source code (python) of DExTER is available at https://gite.lirmm.fr/menichelli/DExTER. Author summary: Gene expression is regulated at different levels and by different mechanisms in Eukaryotes. At the DNA level, transcription factors (TFs) are supposed to play a key role by binding short motifs in promoters or enhancers. In Plasmodium falciparum, the causative agent of severe malaria in humans, different levels of gene regulation are also present, but very few TFs have been identified and validated so far. We propose here a computational method for the identification of a new type of regulatory elements called long regulatory elements (LRE). Contrary to TF motifs, that are usually 6-12bp long, LREs may span dozens or hundreds of base pairs. Moreover, no computational method have been specifically dedicated to their identification until now. We show with our method that, depending on species and conditions, LREs may play important role in gene regulation. For P. falciparum, these elements appear to determine a very large part of gene expression variation in all stages of the parasite life cycle. [ABSTRACT FROM AUTHOR]

Loading of the MCM replicative helicase at origins of replication is a highly regulated process that precedes DNA replication in all eukaryotes. The stoichiometry of MCM loaded at origins has been proposed to be a key determinant of when those origins initiate replication during S phase. Nevertheless, the genome-wide regulation of MCM loading stoichiometry and its direct effect on replication timing remain unclear. In order to investigate why some origins load more MCM than others, we perturbed MCM levels in budding yeast cells and, for the first time, directly measured MCM levels and replication timing in the same experiment. Reduction of MCM levels through degradation of Mcm4, one of the six obligate components of the MCM complex, slowed progression through S phase and increased sensitivity to replication stress. Reduction of MCM levels also led to differential loading at origins during G1, revealing origins that are sensitive to reductions in MCM and others that are not. Sensitive origins loaded less MCM under normal conditions and correlated with a weak ability to recruit the origin recognition complex (ORC). Moreover, reduction of MCM loading at specific origins of replication led to a delay in their replication during S phase. In contrast, overexpression of MCM had no effects on cell cycle progression, relative MCM levels at origins, or replication timing, suggesting that, under optimal growth conditions, cellular MCM levels are not limiting for MCM loading. Our results support a model in which the loading capacity of origins is the primary determinant of MCM stoichiometry in wild-type cells, but that stoichiometry is controlled by origins' ability to recruit ORC and compete for MCM when MCM becomes limiting. Author summary: The coordinated replication of complex eukaryotic genomes is primarily regulated at the level of replication initiation. Therefore, how initiation is regulated in time and space is a active field of research. A leading hypothesis is that the stoichiometry of MCM, the replicative helicase, at replication origins is a significant determinant of initiation timing. However, how MCM stoichiometry at origins is regulated is unknown. We tested two hypothetical mechanisms for the control of MCM stoichiometry: Origin capacity, in which MCM stoichiometry is constrained by the varying capacity of individual origins to load MCM, and the ORC activity, in which MCM stoichiometry is regulated by the varying activity of ORC, the MCM loader, at individual origins. We find that neither mechanism is consistent with our model, but that a hybrid of the two mechanisms is. [ABSTRACT FROM AUTHOR]

Low-cost, high-throughput sequencing has led to an enormous increase in the number of sequenced microbial genomes, with well over 100,000 genomes in public archives today. Automatic genome annotation tools are integral to understanding these organisms, yet older gene finding methods must be retrained on each new genome. We have developed a universal model of prokaryotic genes by fitting a temporal convolutional network to amino-acid sequences from a large, diverse set of microbial genomes. We incorporated the new model into a gene finding system, Balrog (Bacterial Annotation by Learned Representation Of Genes), which does not require genome-specific training and which matches or outperforms other state-of-the-art gene finding tools. Balrog is freely available under the MIT license at https://github.com/salzberg-lab/Balrog. Author summary: Annotating the protein-coding genes in a newly sequenced prokaryotic genome is a critical part of describing their biological function. Relative to eukaryotic genomes, prokaryotic genomes are small and structurally simple, with 90% of their DNA typically devoted to protein-coding genes. Current computational gene finding tools are therefore able to achieve close to 99% sensitivity to known genes using species-specific gene models. Though highly sensitive at finding known genes, all current prokaryotic gene finders also predict large numbers of additional genes, which are labelled as "hypothetical protein" in GenBank and other annotation databases. Many hypothetical gene predictions likely represent true protein-coding sequence, but it is not known how many of them represent false positives. Additionally, all current gene finding tools must be trained specifically for each genome as a preliminary step in order to achieve high sensitivity. This requirement limits their ability to detect genes in fragmented sequences commonly seen in metagenomic samples. We took a data-driven approach to prokaryotic gene finding, relying on the large and diverse collection of already-sequenced genomes. By training a single, universal model of bacterial genes on protein sequences from many different species, we were able to match the sensitivity of current gene finders while reducing the overall number of gene predictions. Our model does not need to be refit on any new genome. Balrog (Bacterial Annotation by Learned Representation of Genes) represents a fundamentally different yet effective method for prokaryotic gene finding. [ABSTRACT FROM AUTHOR]

Pantoea agglomerans is a Gram-negative facultative anaerobic bacillus causing a wide range of opportunistic infections in humans including septicemia, pneumonia, septic arthritis, wound infections and meningitis. To date, the determinants of virulence, antibiotic resistance, metabolic features conferring survival and host-associated pathogenic potential of this bacterium remain largely underexplored. In this study, we sequenced and assembled the whole-genome of P. agglomerans KM1 isolated from kimchi in South Korea. The genome contained one circular chromosome of 4,039,945 bp, 3 mega plasmids, and 2 prophages. The phage-derived genes encoded integrase, lysozyme and terminase. Six CRISPR loci were identified within the bacterial chromosome. Further in-depth analysis showed that the genome contained 13 antibiotic resistance genes conferring resistance to clinically important antibiotics such as penicillin G, bacitracin, rifampicin, vancomycin, and fosfomycin. Genes involved in adaptations to environmental stress were also identified which included factors providing resistance to osmotic lysis, oxidative stress, as well as heat and cold shock. The genomic analysis of virulence factors led to identification of a type VI secretion system, hemolysin, filamentous hemagglutinin, and genes involved in iron uptake and sequestration. Finally, the data provided here show that, the KM1 isolate exerted strong immunostimulatory properties on RAW 264.7 macrophages in vitro. Stimulated cells produced Nitric Oxide (NO) and pro-inflammatory cytokines TNF-α, IL-6 and the anti-inflammatory cytokine IL-10. The upstream signaling for production of TNF-α, IL-6, IL-10, and NO depended on TLR4 and TLR1/2. While production of TNF-α, IL-6 and NO involved solely activation of the NF-κB, IL-10 secretion was largely dependent on NF-κB and to a lesser extent on MAPK Kinases. Taken together, the analysis of the whole-genome and immunostimulatory properties provided in-depth characterization of the P. agglomerans KM1 isolate shedding a new light on determinants of virulence that drive its interactions with the environment, other microorganisms and eukaryotic hosts [ABSTRACT FROM AUTHOR]

Histone variants expand chromatin functions in eukaryote genomes. H2A.B genes are testis-expressed short histone H2A variants that arose in placental mammals. Their biological functions remain largely unknown. To investigate their function, we generated a knockout (KO) model that disrupts all 3 H2A.B genes in mice. We show that H2A.B KO males have globally altered chromatin structure in postmeiotic germ cells. Yet, they do not show impaired spermatogenesis or testis function. Instead, we find that H2A.B plays a crucial role postfertilization. Crosses between H2A.B KO males and females yield embryos with lower viability and reduced size. Using a series of genetic crosses that separate parental and zygotic contributions, we show that the H2A.B status of both the father and mother, but not of the zygote, affects embryonic viability and growth during gestation. We conclude that H2A.B is a novel parental-effect gene, establishing a role for short H2A histone variants in mammalian development. We posit that parental antagonism over embryonic growth drove the origin and ongoing diversification of short histone H2A variants in placental mammals. The unusual short histone variant H2A.B is a novel parental-effect gene that plays an important role in early mammalian development. Parental antagonism over embryonic growth resource allocation may have driven the origin and ongoing diversification of short histone H2A variants in placental mammals. [ABSTRACT FROM AUTHOR]

Transposable elements (TEs) are mobile genetic elements that parasitize basically all eukaryotic species genomes. Due to their complexity, an in-depth TE characterization is only available for a handful of model organisms. In the present study, we performed a de novo and homology-based characterization of TEs in the genomes of 24 mosquito species and investigated their mode of inheritance. More than 40% of the genome of Aedes aegypti, Aedes albopictus, and Culex quinquefasciatus is composed of TEs, varying substantially among Anopheles species (0.13%–19.55%). Class I TEs are the most abundant among mosquitoes and at least 24 TE superfamilies were found. Interestingly, TEs have been extensively exchanged by horizontal transfer (172 TE families of 16 different superfamilies) among mosquitoes in the last 30 million years. Horizontally transferred TEs represents around 7% of the genome in Aedes species and a small fraction in Anopheles genomes. Most of these horizontally transferred TEs are from the three ubiquitous LTR superfamilies: Gypsy, Bel-Pao and Copia. Searching more than 32,000 genomes, we also uncovered transfers between mosquitoes and two different Phyla—Cnidaria and Nematoda—and two subphyla—Chelicerata and Crustacea, identifying a vector, the worm Wuchereria bancrofti, that enabled the horizontal spread of a Tc1-mariner element among various Anopheles species. These data also allowed us to reconstruct the horizontal transfer network of this TE involving more than 40 species. In summary, our results suggest that TEs are frequently exchanged by horizontal transfers among mosquitoes, influencing mosquito's genome size and variability. Author summary: A large fraction of most eukaryotic genomes is derived from DNA fragments that can multiply by inserting themselves in other regions of the genome, generating genomic instability and variability. These fragments are called Transposable Elements (TEs). Since they are a constituent part of the eukaryote genomes, these DNA pieces are normally inherited vertically along with other genes to the host offspring. Host genomes have an arsenal of molecular mechanisms that were naturally selected to control TEs replication avoiding most deleterious damage to the genome and hence leading to long-term TEs inactivation and degeneration. However, TEs sometimes get out of control and invade other species through horizontal transfer, the transfer of genetic information through other means than sexual intercourse. This phenomenon is poorly studied in mosquitoes, a highly important Insect taxa that transmit many human pathogens. We characterized the whole TE content of 24 mosquito genomes and investigated the TE inheritance mode uncovering hundreds of horizontal transfer events among these species and distantly related ones. We also identify a mosquito parasitic filarial worm, that is involved as a transfer vector between many species of mosquitoes and that horizontally transferred TEs contributed significantly to mosquito's genome expansion and variation. [ABSTRACT FROM AUTHOR]

Lateral gene transfer (LGT) has impacted prokaryotic genome evolution, yet the extent to which LGT compromises vertical evolution across individual genes and individual phyla is unknown, as are the factors that govern LGT frequency across genes. Estimating LGT frequency from tree comparisons is problematic when thousands of genomes are compared, because LGT becomes difficult to distinguish from phylogenetic artefacts. Here we report quantitative estimates for verticality across all genes and genomes, leveraging a well-known property of phylogenetic inference: phylogeny works best at the tips of trees. From terminal (tip) phylum level relationships, we calculate the verticality for 19,050,992 genes from 101,422 clusters in 5,655 prokaryotic genomes and rank them by their verticality. Among functional classes, translation, followed by nucleotide and cofactor biosynthesis, and DNA replication and repair are the most vertical. The most vertically evolving lineages are those rich in ecological specialists such as Acidithiobacilli, Chlamydiae, Chlorobi and Methanococcales. Lineages most affected by LGT are the α-, β-, γ-, and δ- classes of Proteobacteria and the Firmicutes. The 2,587 eukaryotic clusters in our sample having prokaryotic homologues fail to reject eukaryotic monophyly using the likelihood ratio test. The low verticality of α-proteobacterial and cyanobacterial genomes requires only three partners—an archaeal host, a mitochondrial symbiont, and a plastid ancestor—each with mosaic chromosomes, to directly account for the prokaryotic origin of eukaryotic genes. In terms of phylogeny, the 100 most vertically evolving prokaryotic genes are neither representative nor predictive for the remaining 97% of an average genome. In search of factors that govern LGT frequency, we find a simple but natural principle: Verticality correlates strongly with gene distribution density, LGT being least likely for intruding genes that must replace a preexisting homologue in recipient chromosomes. LGT is most likely for novel genetic material, intruding genes that encounter no competing copy. Author summary: Because multicellular life is a latecomer in Earth history, most of evolutionary history is microbial evolution. Scientists investigate microbial evolution by studying the evolution of genes. One of the main surprises of the genomic era is the amount of lateral gene transfer that has gone on in prokaryote genome evolution. Gene transfer clouds evolutionary history, but by how much: How lateral and how vertical is the microbial evolutionary process across genes, genomes and lineages? We introduce measures of verticality in genome evolution that permit a ranking of genes and lineages according to their degree of verticality. We show that genes already present in genomes are less likely to be replaced by a newly introduced copy than genes that offer new evolutionary opportunities for the recipient, providing a simple and natural mechanism that limits and promotes lateral gene transfer frequency. Only a very small minority of prokaryotic genes evolve vertically. While the 100 genes that are most widely used to describe the phylogenetic relationships of microbes are indeed the most vertical, they are not at all representative for the evolution of other genes. These findings have broad implications for how we understand the evolutionary process as inferred from gene trees. [ABSTRACT FROM AUTHOR]

Functional analyses of genes are crucial for unveiling biological responses, genetic engineering, and developing new medicines. However, functional analyses have largely been restricted to model organisms, representing a major hurdle for functional studies and industrial applications. To resolve this, comparative genome analyses can be used to provide clues to gene functions as well as their evolutionary history. To this end, we present Prometheus, a web-based omics portal that contains more than 17,215 sequences from prokaryotic and eukaryotic genomes. This portal supports interkingdom comparative analyses via a domain architecture-based gene identification system and Gene Search, and users can easily and rapidly identify single or entire gene sets in specific pathways. Bioinformatics tools for further analyses are provided in Prometheus or through Bio-Express, a cloud-based bioinformatics analysis platform. Prometheus is a new paradigm for comparative analyses of large amounts of genomic information. [ABSTRACT FROM AUTHOR]

Transposable elements (TEs) are mobile genetic elements in eukaryotic genomes. Recent research highlights the important role of TEs in the embryogenesis, neurodevelopment, and immune functions. However, there is a lack of a one-stop and easy to use computational pipeline for expression analysis of both genes and locus-specific TEs from RNA-Seq data. Here, we present GeneTEFlow, a fully automated, reproducible and platform-independent workflow, for the comprehensive analysis of gene and locus-specific TEs expression from RNA-Seq data employing Nextflow and Docker technologies. This application will help researchers more easily perform integrated analysis of both gene and TEs expression, leading to a better understanding of roles of gene and TEs regulation in human diseases. GeneTEFlow is freely available at https://github.com/zhongw2/GeneTEFlow. [ABSTRACT FROM AUTHOR]

Achieving complete and precise genome duplication requires that each genomic segment be replicated only once per cell division cycle. Protecting large eukaryotic genomes from re-replication requires an overlapping set of molecular mechanisms that prevent the first DNA replication step, the DNA loading of MCM helicase complexes to license replication origins, after S phase begins. Previous reports have defined many such origin licensing inhibition mechanisms, but the temporal relationships among them are not clear, particularly with respect to preventing re-replication in G2 and M phases. Using a combination of mutagenesis, biochemistry, and single cell analyses in human cells, we define a new mechanism that prevents re-replication through hyperphosphorylation of the essential MCM loading protein, Cdt1. We demonstrate that Cyclin A/CDK1 can hyperphosphorylate Cdt1 to inhibit MCM re-loading in G2 phase. The mechanism of inhibition is to block Cdt1 binding to MCM independently of other known Cdt1 inactivation mechanisms such as Cdt1 degradation during S phase or Geminin binding. Moreover, our findings suggest that Cdt1 dephosphorylation at the mitosis-to-G1 phase transition re-activates Cdt1. We propose that multiple distinct, non-redundant licensing inhibition mechanisms act in a series of sequential relays through each cell cycle phase to ensure precise genome duplication. Author summary: The initial step of DNA replication is loading the DNA helicase, MCM, onto DNA during the first phase of the cell division cycle. If MCM loading occurs inappropriately onto DNA that has already been replicated, then cells risk DNA re-replication, a source of endogenous DNA damage and genome instability. How mammalian cells prevent any sections of their very large genomes from re-replicating is still not fully understood. We found that the Cdt1 protein, one of the critical MCM loading factors, is inhibited specifically in late cell cycle stages through a mechanism involving protein phosphorylation. This phosphorylation prevents Cdt1 from binding MCM; when Cdt1 cannot be phosphorylated MCM is inappropriately re-loaded onto DNA and cells are prone to re-replication. When cells divide and transition into G1 phase, Cdt1 is then dephosphorylated to re-activate it for MCM loading. Based on these findings we assert that the different mechanisms that cooperate to avoid re-replication are not redundant. Instead, different cell cycle phases are dominated by different re-replication control mechanisms. These findings have implications for understanding how genomes are duplicated precisely once per cell cycle and shed light on how that process is perturbed by changes in Cdt1 levels or phosphorylation activity. [ABSTRACT FROM AUTHOR]

Chromatin regulation of eukaryotic genomes depends on the formation of nucleosome complexes between histone proteins and DNA. Histone variants, which are diversified by sequence or expression pattern, can profoundly alter chromatin properties. While variants in histone H2A and H3 families are well characterized, the extent of diversification of histone H2B proteins is less understood. Here, we report a systematic analysis of the histone H2B family in plants, which have undergone substantial divergence during the evolution of each major group in the plant kingdom. By characterising Arabidopsis H2Bs, we substantiate this diversification and reveal potential functional specialization that parallels the phylogenetic structure of emergent clades in eudicots. In addition, we identify a new class of highly divergent H2B variants, H2B.S, that specifically accumulate during chromatin compaction of dry seed embryos in multiple species of flowering plants. Our findings thus identify unsuspected diverse properties among histone H2B proteins in plants that has manifested into potentially novel groups of histone variants. Author summary: In addition to well-studied variants from core histones families H2A and H3, we report that land plants diversified their H2B family, leading to specialized H2B variants with specific patterns of expression, genomic distributions and properties. [ABSTRACT FROM AUTHOR]

Phylogenetic profiling is a computational method to predict genes involved in the same biological process by identifying protein families which tend to be jointly lost or retained across the tree of life. Phylogenetic profiling has customarily been more widely used with prokaryotes than eukaryotes, because the method is thought to require many diverse genomes. There are now many eukaryotic genomes available, but these are considerably larger, and typical phylogenetic profiling methods require at least quadratic time as a function of the number of genes. We introduce a fast, scalable phylogenetic profiling approach entitled HogProf, which leverages hierarchical orthologous groups for the construction of large profiles and locality-sensitive hashing for efficient retrieval of similar profiles. We show that the approach outperforms Enhanced Phylogenetic Tree, a phylogeny-based method, and use the tool to reconstruct networks and query for interactors of the kinetochore complex as well as conserved proteins involved in sexual reproduction: Hap2, Spo11 and Gex1. HogProf enables large-scale phylogenetic profiling across the three domains of life, and will be useful to predict biological pathways among the hundreds of thousands of eukaryotic species that will become available in the coming few years. HogProf is available at https://github.com/DessimozLab/HogProf. Author summary: Genes that are involved in the same biological process tend to co-evolve. This property is exploited by the technique of phylogenetic profiling, which identifies co-evolving (and therefore likely functionally related) genes through patterns of correlated gene retention and loss in evolution and across species. However, conventional methods to computing and clustering these correlated genes do not scale with increasing numbers of genomes. HogProf is a novel phylogenetic profiling tool built on probabilistic data structures. It allows the user to construct searchable databases containing the evolutionary history of hundreds of thousands of protein families. Such fast detection of coevolution takes advantage of the rapidly increasing amount of genomic data publicly available, and can uncover unknown biological networks and guide in-vivo research and experimentation. We have applied our tool to describe the biological networks underpinning sexual reproduction in eukaryotes. [ABSTRACT FROM AUTHOR]

Scientists are sequencing new genomes at an increasing rate with the goal of associating genome contents with phenotypic traits. After a new genome is sequenced and assembled, structural gene annotation is often the first step in analysis. Despite advances in computational gene prediction algorithms, most eukaryotic genomes still benefit from manual gene annotation. This requires access to good genome browsers to enable annotators to visualize and evaluate multiple lines of evidence (e.g., sequence similarity, RNA sequencing [RNA-Seq] results, gene predictions, repeats) and necessitates many volunteers to participate in the work. To address the technical barriers to creating genome browsers, the Genomics Education Partnership (GEP; https://gep.wustl.edu/) has partnered with the Galaxy Project (https://galaxyproject.org) to develop G-OnRamp (http://g-onramp.org), a web-based platform for creating UCSC Genome Browser Assembly Hubs and JBrowse genome browsers. G-OnRamp also converts a JBrowse instance into an Apollo instance for collaborative genome annotations in research and educational settings. The genome browsers produced can be transferred to the CyVerse Data Store for long-term access. G-OnRamp enables researchers to easily visualize their experimental results, educators to create Course-based Undergraduate Research Experiences (CUREs) centered on genome annotation, and students to participate in genomics research. In the process, students learn about genes/genomes and about how to utilize large datasets. Development of G-OnRamp was guided by extensive user feedback. Sixty-five researchers/educators from >40 institutions participated through in-person workshops, which produced >20 genome browsers now available for research and education. Genome browsers generated for four parasitoid wasp species have been used in a CURE engaging students at 15 colleges and universities. Our assessment results in the classroom demonstrate that the genome browsers produced by G-OnRamp are effective tools for engaging undergraduates in research and in enabling their contributions to the scientific literature in genomics. Expansion of such genomics research/education partnerships will be beneficial to researchers, faculty, and students alike. Author summary: Major projects now underway aim to sequence most of the multicellular organisms on earth (e.g., the Earth Biogenome Project). But obtaining this data is only the beginning. To understand these organisms and how they relate to each other, we need to annotate their genomes (i.e., identify the genes and other features). While computers are essential for this process, most annotation tasks still require or benefit from human analyses. Genome browsers allow annotators to quickly visualize and evaluate multiple lines of evidence to create the best gene models. Hence, annotation of large number of eukaryotic species requires efficient generation of genome browsers and recruitment of many volunteers to participate. We have previously developed a web-based platform (G-OnRamp) to reduce the technical barriers for creating genome browsers. Using the G-OnRamp browsers, we engaged 15 faculty and their students in a Course-based Undergraduate Research Experience (CURE) focused on genome annotation of parasitoid wasp species. We find that G-OnRamp browsers work well in the classroom, and these efforts are beneficial for students and researchers. Students gain research experience, learn about genes and genomes, and learn how to work with large datasets. Researchers obtain high-quality datasets that could not be generated in any other way. [ABSTRACT FROM AUTHOR]

Background: Translationally controlled tumor protein (TCTP) is a conserved, multifunctional protein involved in numerous cellular processes in eukaryotes. Although the functions of TCTP have been investigated sporadically in animals, invertebrates, and plants, few lineage-specific activities of this molecule, have been reported. An exception is in Arabidopsis thaliana, in which TCTP (AtTCTP1) functions in stomatal closuer by regulating microtubule stability. Further, although the development of next-generation sequencing technologies has facilitated the analysis of many eukaryotic genomes in public databases, inter-kingdom comparative analyses using available genome information are comparatively scarce. Methodology: To carry out inter-kingdom comparative analysis of TCTP, TCTP genes were identified from 377 species. Then phylogenetic analysis, prediction of protein structure, molecular docking simulation and molecular dynamics analysis were performed to investigate the evolution of TCTP genes and their binding proteins. Results: A total of 533 TCTP genes were identified from 377 eukaryotic species, including protozoa, fungi, invertebrates, vertebrates, and plants. Phylogenetic and secondary structure analyses reveal lineage-specific evolution of TCTP, and inter-kingdom comparisons highlight the lineage-specific emergence of, or changes in, secondary structure elements in TCTP proteins from different kingdoms. Furthermore, secondary structure comparisons between TCTP proteins within each kingdom, combined with measurements of the degree of sequence conservation, suggest that TCTP genes have evolved to conserve protein secondary structures in a lineage-specific manner. Additional tertiary structure analysis of TCTP-binding proteins and their interacting partners and docking simulations between these proteins further imply that TCTP gene variation may influence the tertiary structures of TCTP-binding proteins in a lineage-specific manner. Conclusions: Our analysis suggests that TCTP has undergone lineage-specific evolution and that structural changes in TCTP proteins may correlate with the tertiary structure of TCTP-binding proteins and their binding partners in a lineage-specific manner. [ABSTRACT FROM AUTHOR]

Membraneless pericentromeric heterochromatin (PCH) domains play vital roles in chromosome dynamics and genome stability. However, our current understanding of 3D genome organization does not include PCH domains because of technical challenges associated with repetitive sequences enriched in PCH genomic regions. We investigated the 3D architecture of Drosophila melanogaster PCH domains and their spatial associations with the euchromatic genome by developing a novel analysis method that incorporates genome-wide Hi-C reads originating from PCH DNA. Combined with cytogenetic analysis, we reveal a hierarchical organization of the PCH domains into distinct "territories." Strikingly, H3K9me2-enriched regions embedded in the euchromatic genome show prevalent 3D interactions with the PCH domain. These spatial contacts require H3K9me2 enrichment, are likely mediated by liquid-liquid phase separation, and may influence organismal fitness. Our findings have important implications for how PCH architecture influences the function and evolution of both repetitive heterochromatin and the gene-rich euchromatin. Author summary: The three dimensional (3D) organization of genomes in cell nuclei can influence a wide variety of genome functions. However, most of our understanding of this critical architecture has been limited to the gene-rich euchromatin, and largely ignores the gene-poor and repeat-rich pericentromeric heterochromatin, or PCH. PCH comprises a large part of most eukaryotic genomes, forms 3D membraneless PCH domains in nuclei, and plays a vital role in chromosome dynamics and genome stability. In this study, we developed a new method that overcomes the technical challenges imposed by the highly repetitive PCH DNA, and generated a comprehensive picture of its 3D organization. Combined with image analyses, we reveal a hierarchical organization of the PCH domains. Surprisingly, we showed that distant euchromatic regions enriched for repressive epigenetic marks also dynamically interact with the main PCH domains. These 3D interactions are likely mediated by liquid-liquid phase separation (similar to how oil and vinegar separate in salad dressing) and the resulting liquid-like fusion events, and can influence the fitness of individuals. Our discoveries have strong implications for how seemingly "junk" DNA could impact functions in the gene-rich euchromatin. [ABSTRACT FROM AUTHOR]

Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs). Intronic TEs often trigger the formation of transcriptionally repressive heterochromatin, even within transcription-permissive chromatin environments. Although TE-bearing introns are widely observed in eukaryotic genomes, their epigenetic states, impacts on gene regulation and function, and their contributions to genetic diversity and evolution, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain intronic heterochromatin, most of which are associated with TEs and repetitive sequences. These heterochromatic introns are longer and highly enriched in promoter-proximal positions. On the other hand, introns also accumulate hypomethylated short TEs. Genes with heterochromatic introns are implicated in various biological functions. Transcription of genes bearing intronic heterochromatin is regulated by an epigenetic mechanism involving the conserved factor OsIBM2, mutation of which results in severe developmental and reproductive defects. Furthermore, we found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study demonstrates that heterochromatin is a common epigenetic feature associated with actively transcribed genes in the rice genome. Author summary: Intronic regions of eukaryotic genomes accumulate many Transposable Elements (TEs) and repeats. These intronic repeats are often targeted by epigenetic silencing mechanisms and form a repressive heterochromatin structure, even within transcriptionally active genes. However, the distribution of TEs in the intragenic regions, and their contributions to genetic diversity and evolution in plant genomes, remain poorly understood. In this study, we investigated the genome-wide distribution of intronic TEs and their epigenetic states in the Oryza sativa genome, where TEs comprise 35% of the genome. We found that over 10% of rice genes contain introns associated with repressive heterochromatin. Genes with heterochromatic introns implicated in various biological functions. The conserved protein OsIBM2 is required for proper transcription of a group of heterochromatin-containing genes. We also found that heterochromatic introns evolve rapidly compared to non-heterochromatic introns. Our study indicates that heterochromatin is a common feature in transcribed genes in the rice genome. [ABSTRACT FROM AUTHOR]

Transposable elements constitute a large fraction of most eukaryotic genomes. Insertion of mobile DNA sequences typically has deleterious effects on host fitness, and thus diverse mechanisms have evolved to control mobile element proliferation. Mobility of the Ty1 retrotransposon in Saccharomyces yeasts is regulated by copy number control (CNC) mediated by a self-encoded restriction factor derived from the Ty1 gag capsid gene that inhibits virus-like particle function. Here, we survey a panel of wild and human-associated strains of S. cerevisiae and S. paradoxus to investigate how genomic Ty1 content influences variation in Ty1 mobility. We observe high levels of mobility for a tester element with a gag sequence from the canonical Ty1 subfamily in permissive strains that either lack full-length Ty1 elements or only contain full-length copies of the Ty1' subfamily that have a divergent gag sequence. In contrast, low levels of canonical Ty1 mobility are observed in restrictive strains carrying full-length Ty1 elements containing a canonical gag sequence. Phylogenomic analysis of full-length Ty1 elements revealed that Ty1' is the ancestral subfamily present in wild strains of S. cerevisiae, and that canonical Ty1 in S. cerevisiae is a derived subfamily that acquired gag from S. paradoxus by horizontal transfer and recombination. Our results provide evidence that variation in the ability of S. cerevisiae and S. paradoxus strains to repress canonical Ty1 transposition via CNC is regulated by the genomic content of different Ty1 subfamilies, and that self-encoded forms of transposon control can spread across species boundaries by horizontal transfer. Author summary: Mobile DNA elements dominate the genomes of many eukaryotes, although their evolutionary dynamics are often not well understood. Here, we describe how the retrovirus-like Ty1 retrotransposon diversified in Saccharomyces yeasts in the face of a potent restriction factor encoded by the Ty1 gag gene. We show that mobility of a Ty1 tester element from S. cerevisiae with a gag sequence from the "canonical" Ty1 subfamily varies substantially across diverse yeast strains. Strains with full-length Ty1 elements carrying a canonical gag sequence in their genomes restrict tester Ty1 mobility while strains containing only elements from the Ty1' subfamily do not, suggesting functional divergence among Ty1 subfamilies. Permissive strains lacking Ty1 elements with canonical gag can be converted to restrictive strains by introducing Ty1 elements with canonical gag into permissive genomes. Phylogenetic analyses show that Ty1' is the ancestral subfamily in wild S. cerevisiae lineages, whereas canonical Ty1 is a derived subfamily found only in human-associated strains. Furthermore, the canonical Ty1 gag arose by recombination between an ancestral Ty1' element and a Ty1 element horizontally transferred from S. paradoxus. Our results highlight how exchange of retrotransposon sequences across species can impact the repertoire and mobility of mobile DNA elements in eukaryotic genomes. [ABSTRACT FROM AUTHOR]

Peony is a famous ornamental and medicinal plant in China, and peony hybrid breeding is an important means of germplasm innovation. However, research on the genome of this species is limited, thereby hindering the genetic and breeding research on peony. In the present study, simple sequence repeat (SSR) locus analysis was performed on expressed sequence tags obtained by the transcriptome sequencing of Paeonia using Microsatellite software. Primers with polymorphism were obtained via polymerase chain reaction amplification and electrophoresis. As a result, a total of 86,195 unigenes were obtained by assembling the transcriptome data of Paeonia. Functional annotations were obtained in seven functional databases including 49,172 (Non-Redundant Protein Sequence Database: 57.05%), 38,352 (Nucleotide Sequence Database: 44.49%), 36,477 (Swiss Prot: 42.32%), 38,905 (Clusters of Orthologous Groups for Eukaryotic Complete Genomes: 45.14%), 37,993 (Kyoto Encyclopedia of Genes and Genomes: 44.08%), 26,832 (Gene Ontology: 31.13%) and 37,758 (Pfam: 43.81%) unigenes. Meanwhile, 21,998 SSR loci were distributed in 17,567 unigenes containing SSR sequences, and the SSR distribution frequency was 25.52%, with an average of one SSR sequence per 4.66 kb. Mononucleotide, dinucleotide, and trinucleotide were the main repeat types, accounting for 55.74%, 25.58%, and 13.21% of the total repeat times, respectively. Forty-five pairs of the 100 pairs of primers selected randomly could amplify clear polymorphic bands. The polymorphic primers of these 45 pairs were used to cluster and analyze 16 species of peony. The new SSR molecular markers can be useful for the study of genetic diversity and marker-assisted breeding of peony. [ABSTRACT FROM AUTHOR]

Interspecific hybridization is the process where closely related species mate and produce offspring with admixed genomes. The genomic revolution has shown that hybridization is common, and that it may represent an important source of novel variation. Although most interspecific hybrids are sterile or less fit than their parents, some may survive and reproduce, enabling the transfer of adaptive variants across the species boundary, and even result in the formation of novel evolutionary lineages. There are two main variants of hybrid species genomes: allopolyploid, which have one full chromosome set from each parent species, and homoploid, which are a mosaic of the parent species genomes with no increase in chromosome number. The establishment of hybrid species requires the development of reproductive isolation against parental species. Allopolyploid species often have strong intrinsic reproductive barriers due to differences in chromosome number, and homoploid hybrids can become reproductively isolated from the parent species through assortment of genetic incompatibilities. However, both types of hybrids can become further reproductively isolated, gaining extrinsic isolation barriers, by exploiting novel ecological niches, relative to their parents. Hybrids represent the merging of divergent genomes and thus face problems arising from incompatible combinations of genes. Thus hybrid genomes are highly dynamic and undergo rapid evolutionary change, including genome stabilization in which selection against incompatible combinations results in fixation of compatible ancestry block combinations within the hybrid species. The potential for rapid adaptation or speciation makes hybrid genomes a particularly exciting subject of in evolutionary biology. Here we summarize how introgressed alleles or hybrid species can establish and how the resulting hybrid genomes evolve. [ABSTRACT FROM AUTHOR]

The sophistication of gene prediction algorithms and the abundance of RNA-based evidence for the maize genome may suggest that manual curation of gene models is no longer necessary. However, quality metrics generated by the MAKER-P gene annotation pipeline identified 17,225 of 130,330 (13%) protein-coding transcripts in the B73 Reference Genome V4 gene set with models of low concordance to available biological evidence. Working with eight graduate students, we used the Apollo annotation editor to curate 86 transcript models flagged by quality metrics and a complimentary method using the Gramene gene tree visualizer. All of the triaged models had significant errors–including missing or extra exons, non-canonical splice sites, and incorrect UTRs. A correct transcript model existed for about 60% of genes (or transcripts) flagged by quality metrics; we attribute this to the convention of elevating the transcript with the longest coding sequence (CDS) to the canonical, or first, position. The remaining 40% of flagged genes resulted in novel annotations and represent a manual curation space of about 10% of the maize genome (~4,000 protein-coding genes). MAKER-P metrics have a specificity of 100%, and a sensitivity of 85%; the gene tree visualizer has a specificity of 100%. Together with the Apollo graphical editor, our double triage provides an infrastructure to support the community curation of eukaryotic genomes by scientists, students, and potentially even citizen scientists. [ABSTRACT FROM AUTHOR]

DNA methylation and epigenetic silencing play important roles in the regulation of transposable elements (TEs) in many eukaryotic genomes. A majority of the maize genome is derived from TEs that can be classified into different orders and families based on their mechanism of transposition and sequence similarity, respectively. TEs themselves are highly methylated and it can be tempting to view them as a single uniform group. However, the analysis of DNA methylation profiles in flanking regions provides evidence for distinct groups of chromatin properties at different TE families. These differences among TE families are reproducible in different tissues and different inbred lines. TE families with varying levels of DNA methylation in flanking regions also show distinct patterns of chromatin accessibility and modifications within the TEs. The differences in the patterns of DNA methylation flanking TE families arise from a combination of non-random insertion preferences of TE families, changes in DNA methylation triggered by the insertion of the TE and subsequent selection pressure. A set of nearly 70,000 TE polymorphisms among four assembled maize genomes were used to monitor the level of DNA methylation at haplotypes with and without the TE insertions. In many cases, TE families with high levels of DNA methylation in flanking sequence are enriched for insertions into highly methylated regions. The majority of the >2,500 TE insertions into unmethylated regions result in changes in DNA methylation in haplotypes with the TE, suggesting the widespread potential for TE insertions to condition altered methylation in conserved regions of the genome. This study highlights the interplay between TEs and the methylome of a major crop species. [ABSTRACT FROM AUTHOR]

Introns are a prevalent feature of eukaryotic genomes, yet their origins and contributions to genome function and evolution remain mysterious. In budding yeast, repression of the highly transcribed intron-containing ribosomal protein genes (RPGs) globally increases splicing of non-RPG transcripts through reduced competition for the spliceosome. We show that under these “hungry spliceosome” conditions, splicing occurs at more than 150 previously unannotated locations we call protointrons that do not overlap known introns. Protointrons use a less constrained set of splice sites and branchpoints than standard introns, including in one case AT-AC in place of GT-AG. Protointrons are not conserved in all closely related species, suggesting that most are not under positive selection and are fated to disappear. Some are found in non-coding RNAs (e. g. CUTs and SUTs), where they may contribute to the creation of new genes. Others are found across boundaries between noncoding and coding sequences, or within coding sequences, where they offer pathways to the creation of new protein variants, or new regulatory controls for existing genes. We define protointrons as (1) nonconserved intron-like sequences that are (2) infrequently spliced, and importantly (3) are not currently understood to contribute to gene expression or regulation in the way that standard introns function. A very few protointrons in S. cerevisiae challenge this classification by their increased splicing frequency and potential function, consistent with the proposed evolutionary process of “intronization”, whereby new standard introns are created. This snapshot of intron evolution highlights the important role of the spliceosome in the expansion of transcribed genomic sequence space, providing a pathway for the rare events that may lead to the birth of new eukaryotic genes and the refinement of existing gene function. [ABSTRACT FROM AUTHOR]

Group II introns are ancient retroelements that significantly shaped the origin and evolution of contemporary eukaryotic genomes. These self-splicing ribozymes share a common ancestor with the telomerase enzyme, the spliceosome machinery as well as the highly abundant spliceosomal introns and non-LTR retroelements. More than half of the human genome thus consists of various elements that evolved from ancient group II introns, which altogether significantly contribute to key functions and genetic diversity in eukaryotes. Similarly, group II intron-related elements in bacteria such as abortive phage infection (Abi) retroelements, diversity generating retroelements (DGRs) and some CRISPR-Cas systems have evolved to confer important functions to their hosts. In sharp contrast, since bacterial group II introns are scarce, irregularly distributed and frequently spread by lateral transfer, they have mainly been considered as selfish retromobile elements with no beneficial function to their host. Here we unveil a new group II intron function that generates genetic diversity at the RNA level in bacterial cells. We demonstrate that Ll.LtrB, the model group II intron from Lactococcus lactis, recognizes specific sequence motifs within cellular mRNAs by base pairing, and invades them by reverse splicing. Subsequent splicing of ectopically inserted Ll.LtrB, through circularization, induces a novel trans-splicing pathway that generates exon 1-mRNA and mRNA-mRNA intergenic chimeras. Our data also show that recognition of upstream alternative circularization sites on intron-interrupted mRNAs release Ll.LtrB circles harboring mRNA fragments of various lengths at their splice junction. Intergenic trans-splicing and alternative circularization both produce novel group II intron splicing products with potential new functions. Overall, this work describes new splicing pathways in bacteria that generate, similarly to the spliceosome in eukaryotes, genetic diversity at the RNA level while providing additional functional and evolutionary links between group II introns, spliceosomal introns and the spliceosome. [ABSTRACT FROM AUTHOR]

The article offers information on the genome size evolution in maize. Topics discussed include component of eukaryotic genomes such as transposable elements; effects of DNA content on organismal phenotype; and the test for natural selection on genome size along with usage of a quantitative genetics framework to conduct a test of local adaptation, recognizing that genome size.

Eukaryotic translation initiation factor 2-alpha kinase (EIF2AK) proteins inhibit protein synthesis at translation initiation level, in response to various stress conditions, including oxidative stress, heme deficiency, osmotic shock, and heat shock. Origin and functional diversification of EIF2AK sequences remain ambiguous. Here we determine the origin and molecular evolution of EIF2AK proteins in lower eukaryotes and studied the molecular basis of divergence among sub-family sequences. Present work emphasized primitive origin of EIF2AK4 sub-family gene in lower eukaryotes of protozoan lineage. Phylogenetic analysis supported common origin and sub-family based classification of EIF2AKs. Functional divergence studies across sub-families revealed several putative amino acid sites, which assist in altered protein interactions of kinase domains. The data can facilitate designing site-directed experimental studies aiming at elucidating diverse functional aspects of kinase domains regarding down-regulation of protein synthesis. [ABSTRACT FROM AUTHOR]