Identifying the mechanisms of intron gain: progress and trends

Transposon-derived introns as an element shaping the structure of eukaryotic genomes

The widely accepted hypothesis postulates that the first spliceosomal introns originated from group II self-splicing introns. However, it is evident that not all spliceosomal introns in the nuclear genes of modern eukaryotes are inherited through vertical transfer of intronic sequences. Several phenomena contribute to the formation of new introns but their most common origin seems to be the insertion of transposable elements. Recent analyses have highlighted instances of mass gains of new introns from transposable elements. These events often coincide with an increase or change in the spliceosome's tolerance to splicing signals, including the acceptance of noncanonical borders. Widespread acquisitions of transposon-derived introns occur across diverse evolutionary lineages, indicating convergent processes. These events, though independent, likely require a similar set of conditions. These conditions include the presence of transposon elements with features enabling their removal at the RNA level as introns and/or the existence of a splicing mechanism capable of excising unusual sequences that would otherwise not be recognized as introns by standard splicing machinery. Herein we summarize those mechanisms across different eukaryotic lineages.

Discovering Intron Gain Events in Humans through Large-Scale Evolutionary Comparisons

The rapid growth in the number of sequenced genomes makes it possible to search for the appearance of entirely new introns in the human lineage. In this study, we compared the genomic sequences for 19,120 human protein-coding genes to a collection of 3493 vertebrate genomes, mapping the patterns of intron alignments onto a phylogenetic tree. This mapping allowed us to trace many intron gain events to precise locations in the tree, corresponding to distinct points in evolutionary history. We discovered 584 intron gain events, all of them relatively recent, in 514 distinct human genes. Among these events, we explored the hypothesis that intronization was the mechanism responsible for intron gain. Intronization events were identified by locating instances where human introns correspond to exonic sequences in homologous vertebrate genes. Although apparently rare, we found three compelling cases of intronization, and for each of those we compared the human protein sequence and structure to homologous genes that lack the introns.

Genome-wide identification of WRKY transcription factor family members in Miscanthus sinensis (Miscanthus sinensis Anderss)

Miscanthus is an emerging sustainable bioenergy crop whose growing environment is subject to many abiotic and biological stresses. WRKY transcription factors play an important role in stress response and growth of biotic and abiotic. To clarify the distribution and expression of the WRKY genes in Miscanthus, it is necessary to classify and phylogenetically analyze the WRKY genes in Miscanthus. The v7.1 genome assembly of Miscanthus was analyzed by constructing an evolutionary tree. In Miscanthus, there are 179 WRKY genes were identified. The 179 MsWRKYs were classified into three groups with conserved gene structure and motif composition. The tissue expression profile of the WRKY genes showed that MsWRKY genes played an essential role in all growth stages of plants. At the early stage of plant development, the MsWRKY gene is mainly expressed in the rhizome of plants. In the middle stage, it is mainly expressed in the leaf. At the end stage, mainly in the stem. According to the results, it showed significant differences in the expression of the MsWRKY in different stages of Miscanthus sinensis. The results of the study contribute to a better understanding of the role of the MsWRKY gene in the growth and development of Miscanthus.

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

The spread of the first introns in proto-eukaryotic paralogs

Spliceosomal introns are a unique feature of eukaryotic genes. Previous studies have established that many introns were present in the protein-coding genes of the last eukaryotic common ancestor (LECA). Intron positions shared between genes that duplicated before LECA could in principle provide insight into the emergence of the first introns. In this study we use ancestral intron position reconstructions in two large sets of duplicated families to systematically identify these ancient paralogous intron positions. We found that 20-35% of introns inferred to have been present in LECA were shared between paralogs. These shared introns, which likely preceded ancient duplications, were wide spread across different functions, with the notable exception of nuclear transport. Since we observed a clear signal of pervasive intron loss prior to LECA, it is likely that substantially more introns were shared at the time of duplication than we can detect in LECA. The large extent of shared introns indicates an early origin of introns during eukaryogenesis and suggests an early origin of a nuclear structure, before most of the other complex eukaryotic features were established.

New Genomic Signals Underlying the Emergence of Human Proto-Genes

De novo genes are novel genes which emerge from non-coding DNA. Until now, little is known about de novo genes’ properties, correlated to their age and mechanisms of emergence. In this study, we investigate four related properties: introns, upstream regulatory motifs, 5′ Untranslated regions (UTRs) and protein domains, in 23,135 human proto-genes. We found that proto-genes contain introns, whose number and position correlates with the genomic position of proto-gene emergence. The origin of these introns is debated, as our results suggest that 41% of proto-genes might have captured existing introns, and 13.7% of them do not splice the ORF. We show that proto-genes which emerged via overprinting tend to be more enriched in core promotor motifs, while intergenic and intronic genes are more enriched in enhancers, even if the TATA motif is most commonly found upstream in these genes. Intergenic and intronic 5′ UTRs of proto-genes have a lower potential to stabilise mRNA structures than exonic proto-genes and established human genes. Finally, we confirm that proteins expressed by proto-genes gain new putative domains with age. Overall, we find that regulatory motifs inducing transcription and translation of previously non-coding sequences may facilitate proto-gene emergence. Our study demonstrates that introns, 5′ UTRs, and domains have specific properties in proto-genes. We also emphasize that the genomic positions of de novo genes strongly impacts these properties.

A new comprehensive annotation of leucine-rich repeat-containing receptors in rice

Oryza sativa (rice) plays an essential food security role for more than half of the world's population. Obtaining crops with high levels of disease resistance is a major challenge for breeders, especially today, given the urgent need for agriculture to be more sustainable. Plant resistance genes are mainly encoded by three large leucine-rich repeat (LRR)-containing receptor (LRR-CR) families: the LRR-receptor-like kinase (LRR-RLK), LRR-receptor-like protein (LRR-RLP) and nucleotide-binding LRR receptor (NLR). Using lrrprofiler, a pipeline that we developed to annotate and classify these proteins, we compared three publicly available annotations of the rice Nipponbare reference genome. The extended discrepancies that we observed for LRR-CR gene models led us to perform an in-depth manual curation of their annotations while paying special attention to nonsense mutations. We then transferred this manually curated annotation to Kitaake, a cultivar that is closely related to Nipponbare, using an optimized strategy. Here, we discuss the breakthrough achieved by manual curation when comparing genomes and, in addition to 'functional' and 'structural' annotations, we propose that the community adopts this approach, which we call 'comprehensive' annotation. The resulting data are crucial for further studies on the natural variability and evolution of LRR-CR genes in order to promote their use in breeding future resilient varieties.

Internally Symmetrical Stwintrons and Related Canonical Introns in Hypoxylaceae Species

Spliceosomal introns are pervasive in eukaryotes. Intron gains and losses have occurred throughout evolution, but the origin of new introns is unclear. Stwintrons are complex intervening sequences where one of the sequence elements (5′-donor, lariat branch point element or 3′-acceptor) necessary for excision of a U2 intron (external intron) is itself interrupted by a second (internal) U2 intron. In Hypoxylaceae, a family of endophytic fungi, we uncovered scores of donor-disrupted stwintrons with striking sequence similarity among themselves and also with canonical introns. Intron–exon structure comparisons suggest that these stwintrons have proliferated within diverging taxa but also give rise to proliferating canonical introns in some genomes. The proliferated (stw)introns have integrated seamlessly at novel gene positions. The recently proliferated (stw)introns appear to originate from a conserved ancestral stwintron characterised by terminal inverted repeats (45–55 nucleotides), a highly symmetrical structure that may allow the formation of a double-stranded intron RNA molecule. No short tandem duplications flank the putatively inserted intervening sequences, which excludes a DNA transposition-based mechanism of proliferation. It is tempting to suggest that this highly symmetrical structure may have a role in intron proliferation by (an)other mechanism(s).

Exon Shuffling Played a Decisive Role in the Evolution of the Genetic Toolkit for the Multicellular Body Plan of Metazoa

Division of labor and establishment of the spatial pattern of different cell types of multicellular organisms require cell type-specific transcription factor modules that control cellular phenotypes and proteins that mediate the interactions of cells with other cells. Recent studies indicate that, although constituent protein domains of numerous components of the genetic toolkit of the multicellular body plan of Metazoa were present in the unicellular ancestor of animals, the repertoire of multidomain proteins that are indispensable for the arrangement of distinct body parts in a reproducible manner evolved only in Metazoa. We have shown that the majority of the multidomain proteins involved in cell-cell and cell-matrix interactions of Metazoa have been assembled by exon shuffling, but there is no evidence for a similar role of exon shuffling in the evolution of proteins of metazoan transcription factor modules. A possible explanation for this difference in the intracellular and intercellular toolkits is that evolution of the transcription factor modules preceded the burst of exon shuffling that led to the creation of the proteins controlling spatial patterning in Metazoa. This explanation is in harmony with the temporal-to-spatial transition hypothesis of multicellularity that proposes that cell differentiation may have predated spatial segregation of cell types in animal ancestors.

Dbr1 functions in mRNA processing, intron turnover and human diseases

Pre-mRNA processing and mRNA stability play direct roles in controlling protein abundance in a cell. Before the mRNA can be translated into a protein, the introns in the pre-mRNA transcripts need to be removed by splicing, such that exons can be ligated together and can code for a protein. In this process, the function of the RNA lariat debranching enzyme or Dbr1 provides a rate-limiting step in the intron turnover process and possibly regulating the production of translation competent mRNAs. Surprising new roles of Dbr1 are emerging in cellular metabolism which extends beyond intron turnover processes, ranging from splicing regulation to translational control. In this review, we highlight the importance of the Dbr1 enzyme, its structure and how anomalies in its function could relate to various human diseases.

The Globin Gene Family in Arthropods: Evolution and Functional Diversity

Globins are small heme-proteins that reversibly bind oxygen. Their most prominent roles in vertebrates are the transport and storage of O₂ for oxidative energy metabolism, but recent research has suggested alternative, non-respiratory globin functions. In the species-rich and ecologically highly diverse taxon of arthropods, the copper-containing hemocyanin is considered the main respiratory protein. However, recent studies have suggested the presence of globin genes and their proteins in arthropod taxa, including model species like Drosophila. To systematically assess the taxonomic distribution, evolution and diversity of globins in arthropods, we systematically searched transcriptome and genome sequence data and found a conserved, widespread occurrence of three globin classes in arthropods: hemoglobin-like (HbL), globin X (GbX), and globin X-like (GbXL) protein lineages. These globin types were previously identified in protostome and deuterostome animals including vertebrates, suggesting their early ancestry in Metazoa. The HbL genes show multiple, lineage-specific gene duplications in all major arthropod clades. Some HbL genes (e.g., Glob2 and 3 of Drosophila) display particularly fast substitution rates, possibly indicating the evolution of novel functions, e.g., in spermatogenesis. In contrast, arthropod GbX and GbXL globin genes show high evolutionary stability: GbXL is represented by a single-copy gene in all arthropod groups except Brachycera, and representatives of the GbX clade are present in all examined taxa except holometabolan insects. GbX and GbXL both show a brain-specific expression. Most arthropod GbX and GbXL proteins, but also some HbL variants, include sequence motifs indicative of potential N-terminal acylation (i.e., N-myristoylation, 3C-palmitoylation). All arthropods except for the brachyceran Diptera harbor at least one such potentially acylated globin copy, confirming the hypothesis of an essential, conserved globin function associated with the cell membrane. In contrast to other animals, the fourth ancient globin lineage, represented by neuroglobin, appears to be absent in arthropods, and the putative arthropod orthologs of the fifth metazoan globin lineage, androglobin, lack a recognizable globin domain. Thus, the remarkable evolutionary stability of some globin variants is contrasted by occasional dynamic gene multiplication or even loss of otherwise strongly conserved globin lineages in arthropod phylogeny.

Alternative Splicing by NOVA Factors: From Gene Expression to Cell Physiology and Pathology

NOVA1 and NOVA2, the two members of the NOVA family of alternative splicing factors, bind YCAY clusters of pre-mRNAs and assemble spliceosomes to induce the maintenance/removal of introns and exons, thus governing the development of mRNAs. Members of other splicing families operate analogously. Activity of NOVAs accounts for up to 700 alternative splicing events per cell, taking place both in the nucleus (co-transcription of mRNAs) and in the cytoplasm. Brain neurons express high levels of NOVAs, with NOVA1 predominant in cerebellum and spinal cord, NOVA2 in the cortex. Among brain physiological processes NOVAs play critical roles in axon pathfinding and spreading, structure and function of synapses, as well as the regulation of surface receptors and voltage-gated channels. In pathology, NOVAs contribute to neurodegenerative diseases and epilepsy. In vessel endothelial cells, NOVA2 is essential for angiogenesis, while in adipocytes, NOVA1 contributes to regulation of thermogenesis and obesity. In many cancers NOVA1 and also NOVA2, by interacting with specific miRNAs and by additional mechanisms, activate oncogenic roles promoting cell proliferation, colony formation, migration, and invasion. In conclusion, NOVAs regulate cell functions of physiological and pathological nature. Single cell identification and distinction, and new therapies addressed to NOVA targets might be developed in the near future.

Complex intron generation in the yeast genus Lipomyces

In primary transcripts of eukaryotic nuclear genes, coding sequences are often interrupted by U2-type introns. Such intervening sequences can constitute complex introns excised by consecutive splicing reactions. The origin of spliceosomal introns is a vexing problem. Sequence variation existent across fungal taxa provides means to study their structure and evolution. In one class of complex introns called [D] stwintrons, an (internal) U2 intron is nested within the 5'-donor element of another (external) U2 intron. In the gene for a reticulon-like protein in species of the ascomycete yeast genus Lipomyces, the most 5' terminal intron position is occupied by one of three complex intervening sequences consistent of differently nested U2 intron units, as demonstrated in L. lipofer, L. suomiensis, and L. starkeyi. In L. starkeyi, the donor elements of the constituent introns are abutting and the complex intervening sequence can be excised alternatively either with one standard splicing reaction or, as a [D] stwintron, by two consecutive reactions. Our work suggests how [D] stwintrons could emerge by the appearance of new functional splice sites within an extant intron. The stepwise stwintronisation mechanism may involve duplication of the functional intron donor element of the ancestor intron.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

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.

The protein coding information in eukaryotic genes is broken by intervening sequences called introns that are removed from RNA during transcription by a large protein-RNA complex called the spliceosome. Where introns come from and how the spliceosome contributes to genome evolution are open questions. In this study, we find more than 150 new places in the yeast genome that are recognized by the spliceosome and spliced out as introns. Since they appear to have arisen very recently in evolution by sequence drift and do not appear to contribute to gene expression or its regulation, we call these protointrons. Protointrons are found in both protein-coding and non-coding RNAs and are not efficiently removed by the splicing machinery. Although most protointrons are not conserved and will likely disappear as evolution proceeds, a few are spliced more efficiently, and are located where they might begin to play functional roles in gene expression, as predicted by the proposed process of intronization. The challenge now is to understand how spontaneously appearing splicing events like protointrons might contribute to the creation of new genes, new genetic controls, and new protein isoforms as genomes evolve.

Rapidly evolving protointrons in Saccharomyces genomes revealed by a hungry spliceosome

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.

The Divergent Genomes of Teleosts

Boasting nearly 30,000 species, teleosts account for half of all extant vertebrates and approximately 98% of all ray-finned fish species (Actinopterygii). Teleosts are also the largest and most diverse group of vertebrates, exhibiting an astonishing level of morphological, physiological, and behavioral diversity. Previous studies had indicated that the teleost lineage has experienced an additional whole-genome duplication event. Recent comparative genomic analyses of teleosts and other bony vertebrates using spotted gar (a nonteleost ray-finned fish) and elephant shark (a cartilaginous fish) as outgroups have revealed several divergent features of teleost genomes. These include an accelerated evolutionary rate of protein-coding and nucleotide sequences, a higher rate of intron turnover, loss of many potential cis-regulatory elements and shorter conserved syntenic blocks. A combination of these divergent genomic features might have contributed to the evolution of the amazing phenotypic diversity and morphological innovations of teleosts.

Microsporidian Introns Retained against a Background of Genome Reduction: Characterization of an Unusual Set of Introns

Spliceosomal introns are ubiquitous features of eukaryotic genomes, but the mechanisms responsible for their loss and gain are difficult to identify. Microsporidia are obligate intracellular parasites that have significantly reduced genomes and, as a result, have lost many if not all of their introns. In the microsporidian Encephalitozoon cuniculi, a relatively long intron was identified and was spliced at higher levels than the remaining introns. This long intron is part of a set of unique introns in two unrelated genes that show high levels of sequence conservation across diverse microsporidia. The introns possess a unique internal conserved region, which overlaps with a shared, predicted stem–loop structure. The unusual similarity and retention of these long introns in reduced microsporidian genomes could indicate that these introns function similarly, are homologous, or both. Regardless, the significant genome reduction in microsporidia provides a rare opportunity to understand intron evolution.

Origin and Spread of Spliceosomal Introns: Insights from the Fungal Clade Zymoseptoria

Spliceosomal introns are a key feature of eukaryote genome architecture and have been proposed to originate from selfish group II introns from an endosymbiotic bacterium, that is, the ancestor of mitochondria. However, the mechanisms underlying the wide spread of spliceosomal introns across eukaryotic genomes have been obscure. In this study, we characterize the dynamic evolution of spliceosomal introns in the fungal genus Zymoseptoria at different evolutionary scales, that is, within a genome, among conspecific strains within species, and between different species. Within the genome, spliceosomal introns can proliferate in unrelated genes and intergenic regions. Among conspecific strains, spliceosomal introns undergo rapid turnover (gains and losses) and frequent sequence exchange between geographically distinct strains. Furthermore, spliceosomal introns could undergo introgression between distinct species, which can further promote intron invasion and proliferation. The dynamic invasion and proliferation processes of spliceosomal introns resemble the life cycles of mobile selfish (group I/II) introns, and these intron movements, at least in part, account for the dramatic processes of intron gain and intron loss during eukaryotic evolution.

Tracing the evolution of the heterotrimeric G protein α subunit in Metazoa

BACKGROUND

Heterotrimeric G proteins are fundamental signaling proteins composed of three subunits, Gα and a Gβγ dimer. The role of Gα as a molecular switch is critical for transmitting and amplifying intracellular signaling cascades initiated by an activated G protein Coupled Receptor (GPCR). Despite their biochemical and therapeutic importance, the study of G protein evolution has been limited to the scope of a few model organisms. Furthermore, of the five primary Gα subfamilies, the underlying gene structure of only two families has been thoroughly investigated outside of Mammalia evolution. Therefore our understanding of Gα emergence and evolution across phylogeny remains incomplete.

RESULTS

We have computationally identified the presence and absence of every Gα gene (GNA-) across all major branches of Deuterostomia and evaluated the conservation of the underlying exon-intron structures across these phylogenetic groups. We provide evidence of mutually exclusive exon inclusion through alternative splicing in specific lineages. Variations of splice site conservation and isoforms were found for several paralogs which coincide with conserved, putative motifs of DNA-/RNA-binding proteins. In addition to our curated gene annotations, within Primates, we identified 15 retrotranspositions, many of which have undergone pseudogenization. Most importantly, we find numerous deviations from previous findings regarding the presence and absence of individual GNA- genes, nuanced differences in phyla-specific gene copy numbers, novel paralog duplications and subsequent intron gain and loss events.

CONCLUSIONS

Our curated annotations allow us to draw more accurate inferences regarding the emergence of all Gα family members across Metazoa and to present a new, updated theory of Gα evolution. Leveraging this, our results are critical for gaining new insights into the co-evolution of the Gα subunit and its many protein binding partners, especially therapeutically relevant G protein - GPCR signaling pathways which radiated in Vertebrata evolution.

The Role of Transposable Elements in Emergence of Metazoa

Systems initially emerged for protecting genomes against insertions of transposable elements and represented by mechanisms of splicing regulation, RNA-interference, and epigenetic factors have played a key role in the evolution of animals. Many studies have shown inherited transpositions of mobile elements in embryogenesis and preservation of their activities in certain tissues of adult organisms. It was supposed that on the emergence of Metazoa the self-regulation mechanisms of transposons related with the gene networks controlling their activity could be involved in intercellular cell coordination in the cascade of successive divisions with differentiated gene expression for generation of tissues and organs. It was supposed that during evolution species-specific features of transposons in the genomes of eukaryotes could form the basis for creation of dynamically related complexes of systems for epigenetic regulation of gene expression. These complexes could be produced due to the influence of noncoding transposon-derived RNAs on DNA methylation, histone modifications, and processing of alternative splicing variants, whereas the mobile elements themselves could be directly involved in the regulation of gene expression in cis and in trans. Transposons are widely distributed in the genomes of eukaryotes; therefore, their activation can change the expression of specific genes. In turn, this can play an important role in cell differentiation during ontogenesis. It is supposed that transposons can form a species-specific pattern for control of gene expression, and that some variants of this pattern can be favorable for adaptation. The presented data indicate the possible influence of transposons in karyotype formation. It is supposed that transposon localization relative to one another and to protein-coding genes can influence the species-specific epigenetic regulation of ontogenesis.

Emergence and loss of spliceosomal twin introns

BACKGROUND

In the primary transcript of nuclear genes, coding sequences-exons-usually alternate with non-coding sequences-introns. In the evolution of spliceosomal intron-exon structure, extant intron positions can be abandoned and new intron positions can be occupied. Spliceosomal twin introns ("stwintrons") are unconventional intervening sequences where a standard "internal" intron interrupts a canonical splicing motif of a second, "external" intron. The availability of genome sequences of more than a thousand species of fungi provides a unique opportunity to study spliceosomal intron evolution throughout a whole kingdom by means of molecular phylogenetics.

RESULTS

A new stwintron was encountered in Aspergillus nidulans and Aspergillus niger. It is present across three classes of Leotiomyceta in the transcript of a well-conserved gene encoding a putative lipase (lipS). It occupies the same position as a standard intron in the orthologue gene in species of the early divergent classes of the Pezizomycetes and the Orbiliomycetes, suggesting that an internal intron has appeared within a pre-extant intron. On the other hand, the stwintron has been lost from certain taxa in Leotiomycetes and Eurotiomycetes at several occasions, most likely by a mechanism involving reverse transcription and homologous recombination. Another ancient stwintron present across whole Pezizomycotina orders-in the transcript of the bifunctional biotin biosynthesis gene bioDA-occurs at the same position as a standard intron in many species of non-Dikarya. Nevertheless, also the bioDA stwintron has disappeared from certain lineages within the taxa where it occurs, i.e., Sordariomycetes and Botryosphaeriales. Intriguingly, only the internal intron was lost from the Sordariomycetes bioDA stwintron at all but one occasion, leaving a standard intron in the same position, while where the putative lipase stwintron was lost, no intronic sequences remain.

CONCLUSIONS

Molecular phylogeny of the peptide product was used to monitor the existence and fate of a stwintron in the transcripts of two neatly defined fungal genes, encoding well conserved proteins. Both defining events-stwintron emergence and loss-can be explained with extant models for intron insertion and loss. We thus demonstrate that stwintrons can serve as model systems to study spliceosomal intron evolution.

Ancestry & molecular evolutionary analyses of heat shock protein 47 kDa (HSP47/SERPINH1)

HSP47/SERPINH1 is key-regulator for collagen biosynthesis and its structural assembly. To date, there is no comprehensive study on the phylogenetic history of HSP47. Herein we illustrate the evolutionary history of HSP47/SERPINH1 along with sequence, structural and syntenic traits for HSP47/SERPINH1. We have identified ancestral HSP47/SERPINH1 locus in Japanese lamprey (Lethenteron japonicum). This gene remains on the same or similar locus for ~500 million years (MY), but chromosomal duplication was observed in ray-finned fishes, leading into three sets of three sets (I-III) of HSP47/SERPINH1. Two novel introns were inserted at the positions 36b and 102b in the first exon of only HSP47_1 gene from the selected ray-finned fishes. On the evolutionary time scale, the events of HSP47 duplications took placed between 416–360 MY ago (MYA) while intron insertion dates back to 231–190 MYA after early divergence of ray-finned fishes.

Uncovering missing pieces: duplication and deletion history of arrestins in deuterostomes

BACKGROUND

The cytosolic arrestin proteins mediate desensitization of activated G protein-coupled receptors (GPCRs) via competition with G proteins for the active phosphorylated receptors. Arrestins in active, including receptor-bound, conformation are also transducers of signaling. Therefore, this protein family is an attractive therapeutic target. The signaling outcome is believed to be a result of structural and sequence-dependent interactions of arrestins with GPCRs and other protein partners. Here we elucidated the detailed evolution of arrestins in deuterostomes.

RESULTS

Identity and number of arrestin paralogs were determined searching deuterostome genomes and gene expression data. In contrast to standard gene prediction methods, our strategy first detects exons situated on different scaffolds and then solves the problem of assigning them to the correct gene. This increases both the completeness and the accuracy of the annotation in comparison to conventional database search strategies applied by the community. The employed strategy enabled us to map in detail the duplication- and deletion history of arrestin paralogs including tandem duplications, pseudogenizations and the formation of retrogenes. The two rounds of whole genome duplications in the vertebrate stem lineage gave rise to four arrestin paralogs. Surprisingly, visual arrestin ARR3 was lost in the mammalian clades Afrotheria and Xenarthra. Duplications in specific clades, on the other hand, must have given rise to new paralogs that show signatures of diversification in functional elements important for receptor binding and phosphate sensing.

CONCLUSION

The current study traces the functional evolution of deuterostome arrestins in unprecedented detail. Based on a precise re-annotation of the exon-intron structure at nucleotide resolution, we infer the gain and loss of paralogs and patterns of conservation, co-variation and selection.

Exchange of Genetic Sequences Between Viruses and Hosts

Although genetic transfer between viruses and vertebrate hosts occurs less frequently than gene flow between bacteriophages and prokaryotes, it is extensive and has affected the evolution of both parties. With retroviruses, the integration of proviral DNA into chromosomal DNA can result in the activation of adjacent host gene expression and in the transduction of host transcripts into retroviral genomes as oncogenes. Yet in contrast to lysogenic phage, there is little evidence that viral oncogenes persist in a chain of natural transmission or that retroviral transduction is a significant driver of the horizontal spread of host genes. Conversely, integration of proviruses into the host germ line has generated endogenous retroviral genomes (ERV) in all vertebrate genomes sequenced to date. Some of these genomes retain potential infectivity and upon reactivation may transmit to other host species. During mammalian evolution, sequences of retroviral origin have been repurposed to serve host functions, such as the viral envelope glycoproteins crucial to the development of the placenta. Beyond retroviruses, DNA viruses with complex genomes have acquired numerous genes of host origin which influence replication, pathogenesis and immune evasion, while host species have accumulated germline sequences of both DNA and RNA viruses. A codicil is added on lateral transmission of cancer cells between hosts and on migration of host mitochondria into cancer cells.

Mechanism for DNA transposons to generate introns on genomic scales

The discovery of introns four decades ago was one of the most unexpected findings in molecular biology. Introns are sequences interrupting genes that must be removed as part of messenger RNA production. Genome sequencing projects have shown that most eukaryotic genes contain at least one intron, and frequently many. Comparison of these genomes reveals a history of long evolutionary periods during which few introns were gained, punctuated by episodes of rapid, extensive gain. However, although several detailed mechanisms for such episodic intron generation have been proposed, none has been empirically supported on a genomic scale. Here we show how short, non-autonomous DNA transposons independently generated hundreds to thousands of introns in the prasinophyte Micromonas pusilla and the pelagophyte Aureococcus anophagefferens. Each transposon carries one splice site. The other splice site is co-opted from the gene sequence that is duplicated upon transposon insertion, allowing perfect splicing out of the RNA. The distributions of sequences that can be co-opted are biased with respect to codons, and phasing of transposon-generated introns is similarly biased. These transposons insert between pre-existing nucleosomes, so that multiple nearby insertions generate nucleosome-sized intervening segments. Thus, transposon insertion and sequence co-option may explain the intron phase biases and prevalence of nucleosome-sized exons observed in eukaryotes. Overall, the two independent examples of proliferating elements illustrate a general DNA transposon mechanism that can plausibly account for episodes of rapid, extensive intron gain during eukaryotic evolution.

Intron gain by tandem genomic duplication: a novel case in a potato gene encoding RNA-dependent RNA polymerase

The origin and subsequent accumulation of spliceosomal introns are prominent events in the evolution of eukaryotic gene structure. However, the mechanisms underlying intron gain remain unclear because there are few proven cases of recently gained introns. In an RNA-dependent RNA polymerase (RdRp) gene, we found that a tandem duplication occurred after the divergence of potato and its wild relatives among other Solanum plants. The duplicated sequence crosses the intron-exon boundary of the first intron and the second exon. A new intron was detected at this duplicated region, and it includes a small previously exonic segment of the upstream copy of the duplicated sequence and the intronic segment of the downstream copy of the duplicated sequence. The donor site of this new intron was directly obtained from the small previously exonic segment. Most of the splicing signals were inherited directly from the parental intron/exon structure, including a putative branch site, the polypyrimidine tract, the 3' splicing site, two putative exonic splicing enhancers, and the GC contents differed between the intron and exon. In the widely cited model of intron gain by tandem genomic duplication, the duplication of an AGGT-containing exonic segment provides the GT and AG splicing sites for the new intron. Our results illustrate that the tandem duplication model of intron gain should be diverse in terms of obtaining the proper splicing signals.

Genome Sequencing of the Phytoseiid Predatory Mite Metaseiulus occidentalis Reveals Completely Atomized Hox Genes and Superdynamic Intron Evolution

Metaseiulus occidentalis is an eyeless phytoseiid predatory mite employed for the biological control of agricultural pests including spider mites. Despite appearances, these predator and prey mites are separated by some 400 Myr of evolution and radically different lifestyles. We present a 152-Mb draft assembly of the M. occidentalis genome: Larger than that of its favored prey, Tetranychus urticae, but considerably smaller than those of many other chelicerates, enabling an extremely contiguous and complete assembly to be built-the best arachnid to date. Aided by transcriptome data, genome annotation cataloged 18,338 protein-coding genes and identified large numbers of Helitron transposable elements. Comparisons with other arthropods revealed a particularly dynamic and turbulent genomic evolutionary history. Its genes exhibit elevated molecular evolution, with strikingly high numbers of intron gains and losses, in stark contrast to the deer tick Ixodes scapularis Uniquely among examined arthropods, this predatory mite's Hox genes are completely atomized, dispersed across the genome, and it encodes five copies of the normally single-copy RNA processing Dicer-2 gene. Examining gene families linked to characteristic biological traits of this tiny predator provides initial insights into processes of sex determination, development, immune defense, and how it detects, disables, and digests its prey. As the first reference genome for the Phytoseiidae, and for any species with the rare sex determination system of parahaploidy, the genome of the western orchard predatory mite improves genomic sampling of chelicerates and provides invaluable new resources for functional genomic analyses of this family of agriculturally important mites.

Spliceosomal intronogenesis

The presence of intervening sequences, termed introns, is a defining characteristic of eukaryotic nuclear genomes. Once transcribed into pre-mRNA, these introns must be removed within the spliceosome before export of the processed mRNA to the cytoplasm, where it is translated into protein. Although intron loss has been demonstrated experimentally, several mysteries remain regarding the origin and propagation of introns. Indeed, documented evidence of gain of an intron has only been suggested by phylogenetic analyses. We report the use of a strategy that detects selected intron gain and loss events. We have experimentally verified, to our knowledge, the first demonstrations of intron transposition in any organism. From our screen, we detected two separate intron gain events characterized by the perfect transposition of a reporter intron into the yeast genes RPL8B and ADH2, respectively. We show that the newly acquired introns are able to be removed from their respective pre-mRNAs by the spliceosome. Additionally, the novel allele, RPL8Bint, is functional when overexpressed within the genome in a strain lacking the Rpl8 paralogue RPL8A, demonstrating that the gene targeted for intronogenesis is functional.

Deltex1 Polymorphisms Are Associated with Hepatitis B Vaccination Non-Response in Southwest China

Background

Hepatitis B vaccination is the most important tool available for preventing hepatitis B virus (HBV) infection and reducing the prevalence of infection. However, epidemiological studies have demonstrated that morethan 5% of patients exhibit a non- or hypo-response to the HBV vaccine. Genetic variations associated with T cell immunity contribute to the immune response to HBV vaccination. The deltex 1 (DTX1) gene is involved in T cell anergy, which may also be associated with the immune response to the HBV vaccination.

Methods

We detected 10 single nucleotide polymorphisms (SNPs) in or around the DTX1 gene in 601 infants out of a population from Southwest China, including 299 high responders(HRs; HBsAb > 100 mIU/mL) and 302 non-responders (NRs; HBsAb < 10 mIU/mL). An additional validation study was performed, comprising 230 adult patients(135 HRs and 95 NRs) from Southwest China.

Results

This study found that the minor allele ‘G’ of rs2384077 (adjusted p = 2.63E^-04,) and the minor allele ‘C’ of rs10744794 (adjusted p = 3.69E^-04) in the first intron of the DTX1 gene were remarkably associated with the immune response to HBV vaccination in both infant and adult populations. Moreover, a subsequent analysis indicated that haplotypes (A-T, G-C) of the two SNPs were significantly associated with the immune response to HBV vaccination.

Conclusions

Two SNPs (rs2384077 and rs10744794) in an intron of DTX1 and the linkage disequilibrium (LD) block are significantly associated with the immune response to HBV vaccination. The functional element annotation of the LD block between the two SNPs contains four transcriptional regulatory elements. The results suggest that these two SNPs may be involved in the immune response to HBV vaccination.

Evaluation of the mechanisms of intron loss and gain in the social amoebae Dictyostelium

Background

Spliceosomal introns are a common feature of eukaryotic genomes. To approach a comprehensive understanding of intron evolution on Earth, studies should look beyond repeatedly studied groups such as animals, plants, and fungi. The slime mold Dictyostelium belongs to a supergroup of eukaryotes not covered in previous studies.

Results

We found 441 precise intron losses in Dictyostelium discoideum and 202 precise intron losses in Dictyostelium purpureum. Consistent with these observations, Dictyostelium discoideum was found to have significantly more copies of reverse transcriptase genes than Dictyostelium purpureum. We also found that the lost introns are significantly further from the 5′ end of genes than the conserved introns. Adjacent introns were prone to be lost simultaneously in Dictyostelium discoideum. In both Dictyostelium species, the exonic sequences flanking lost introns were found to have a significantly higher GC content than those flanking conserved introns. Together, these observations support a reverse-transcription model of intron loss in which intron losses were caused by gene conversion between genomic DNA and cDNA reverse transcribed from mature mRNA. We also identified two imprecise intron losses in Dictyostelium discoideum that may have resulted from genomic deletions. Ninety-eight putative intron gains were also observed. Consistent with previous studies of other lineages, the source sequences were found in only a small number of cases, with only two instances of intron gain identified in Dictyostelium discoideum.

Conclusions

Although they diverged very early from animals and fungi, Dictyostelium species have similar mechanisms of intron loss.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0567-y) contains supplementary material, which is available to authorized users.

A Recently Evolved Alternative Splice Site in the BRANCHED1a Gene Controls Potato Plant Architecture

Amplification and diversification of transcriptional regulators that control development is a driving force of morphological evolution. A major source of protein diversity is alternative splicing, which leads to the generation of different isoforms from a single gene. The mechanisms and timing of intron evolution nonetheless remain unclear, and the functions of alternative splicing-generated protein isoforms are rarely studied. In Solanum tuberosum, the BRANCHED1a (BRC1a) gene encodes a TCP transcription factor that controls lateral shoot outgrowth. Here, we report the recent evolution in Solanum of an alternative splice site in BRC1a that leads to the generation of two BRC1a protein isoforms with distinct C-terminal regions, BRC1a(Long) and BRC1a(Short), encoded by unspliced and spliced mRNA, respectively. The BRC1a(Long) C-terminal region has a strong activation domain, whereas that of BRC1a(S) lacks an activation domain and is predicted to form an amphipathic helix, the H domain, which prevents protein nuclear targeting. BRC1a(Short) is thus mainly cytoplasmic, while BRC1a(Long) is mainly nuclear. BRC1a(Long) functions as a transcriptional activator, whereas BRC1a(Short) appears to have no transcriptional activity. Moreover, BRC1a(Short) can heterodimerize with BRC1a(Long) and act as a dominant-negative factor; it increases BRC1a(Long) concentration in cytoplasm and reduces its transcriptional activity. This alternative splicing mechanism is regulated by hormones and external stimuli that control branching. The evolution of a new alternative splicing site and a novel protein domain in Solanum BRC1a led to a multi-level mechanism of post-transcriptional and post-translational BRC1a regulation that effectively modulates its branch suppressing activity in response to environmental and endogenous cues.

Bayesian phylogeny analysis of vertebrate serpins illustrates evolutionary conservation of the intron and indels based six groups classification system from lampreys for ∼500 MY

The serpin superfamily is characterized by proteins that fold into a conserved tertiary structure and exploits a sophisticated and irreversible suicide-mechanism of inhibition. Vertebrate serpins are classified into six groups (V1-V6), based on three independent biological features-genomic organization, diagnostic amino acid sites and rare indels. However, this classification system was based on the limited number of mammalian genomes available. In this study, several non-mammalian genomes are used to validate this classification system using the powerful Bayesian phylogenetic method. This method supports the intron and indel based vertebrate classification and proves that serpins have been maintained from lampreys to humans for about 500 MY. Lampreys have fewer than 10 serpins, which expand into 36 serpins in humans. The two expanding groups V1 and V2 have SERPINB1/SERPINB6 and SERPINA8/SERPIND1 as the ancestral serpins, respectively. Large clusters of serpins are formed by local duplications of these serpins in tetrapod genomes. Interestingly, the ancestral HCII/SERPIND1 locus (nested within PIK4CA) possesses group V4 serpin (A2APL1, homolog of α 2-AP/SERPINF2) of lampreys; hence, pointing to the fact that group V4 might have originated from group V2. Additionally in this study, details of the phylogenetic history and genomic characteristics of vertebrate serpins are revisited.

Intron Invasions Trace Algal Speciation and Reveal Nearly Identical Arctic and Antarctic Micromonas Populations

Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence–absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture (Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica—a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga.

A survey of transposable element classification systems--a call for a fundamental update to meet the challenge of their diversity and complexity

The increase of publicly available sequencing data has allowed for rapid progress in our understanding of genome composition. As new information becomes available we should constantly be updating and reanalyzing existing and newly acquired data. In this report we focus on transposable elements (TEs) which make up a significant portion of nearly all sequenced genomes. Our ability to accurately identify and classify these sequences is critical to understanding their impact on host genomes. At the same time, as we demonstrate in this report, problems with existing classification schemes have led to significant misunderstandings of the evolution of both TE sequences and their host genomes. In a pioneering publication Finnegan (1989) proposed classifying all TE sequences into two classes based on transposition mechanisms and structural features: the retrotransposons (class I) and the DNA transposons (class II). We have retraced how ideas regarding TE classification and annotation in both prokaryotic and eukaryotic scientific communities have changed over time. This has led us to observe that: (1) a number of TEs have convergent structural features and/or transposition mechanisms that have led to misleading conclusions regarding their classification, (2) the evolution of TEs is similar to that of viruses by having several unrelated origins, (3) there might be at least 8 classes and 12 orders of TEs including 10 novel orders. In an effort to address these classification issues we propose: (1) the outline of a universal TE classification, (2) a set of methods and classification rules that could be used by all scientific communities involved in the study of TEs, and (3) a 5-year schedule for the establishment of an International Committee for Taxonomy of Transposable Elements (ICTTE).

Origin and evolution of a new retained intron on the vulcan gene in Drosophila melanogaster subgroup species

Although numerous intron gains have been discovered, the mechanisms of intron creation have proven to be elusive. Previous study revealed that the vulcan gene of Drosophila melanogaster contained four exons in its coding region. In the current study, a newly created intron (Intron L) was identified on exon 2 of vulcan in D. melanogaster by comparing expression sequence tags. The RT-PCR experiment revealed that Intron L was associated with intron retention, in which two alternative transcripts of the gene differ by the inclusion or removal of an intron. It was found that Intron L was created by intronization of exonic sequence, and its donor and acceptor splice sites were created by synonymous mutation, leading to the origin of a new vulcan protein that is 22 amino acids shorter than the previously reported vulcan protein. Moreover, to track the origin of Intron L, 36 orthologous genes of species of Drosophila were cloned or annotated, and phylogenetic analysis was carried out. It indicated that the common ancestor of D. melangaster subgroup species created Intron L about 15 million years ago.

Intron evolution in Neurospora: the role of mutational bias and selection

We used comparative and population genomics to study intron evolutionary dynamics in the fungal model genus Neurospora. For our investigation, we used well-annotated genomes of N. crassa, N. discreta, and N. tetrasperma, and 92 resequenced genomes of N. tetrasperma from natural populations. By analyzing the four well-annotated genomes, we identified 9495 intron sites in 7619 orthologous genes. Our data supports nonhomologous end joining (NHEJ) and tandem duplication as mechanisms for intron gains in the genus and the RT-mRNA process as a mechanism for intron loss. We found a moderate intron gain rate (5.78–6.89 × 10⁻¹³ intron gains per nucleotide site per year) and a high intron loss rate (7.53–13.76 × 10⁻¹⁰ intron losses per intron sites per year) as compared to other eukaryotes. The derived intron gains and losses are skewed to high frequencies, relative to neutral SNPs, in natural populations of N. tetrasperma, suggesting that selection is involved in maintaining a high intron turnover. Furthermore, our analyses of the association between intron population-level frequency and genomic features suggest that selection is involved in shaping a 5′ intron position bias and a low intron GC content. However, intron sequence analyses suggest that the gained introns were not exposed to recent selective sweeps. Taken together, this work contributes to our understanding of the importance of mutational bias and selection in shaping the intron distribution in eukaryotic genomes.

Molecular evolution and functional divergence of X-intrinsic protein genes in plants

X-intrinsic proteins (XIPs) are a novel class of major intrinsic proteins found in diverse organisms. Recently, XIP genes have been reported to be involved in the transport of a wide range of hydrophobic solutes; however, the evolutionary forces driving their structural and functional divergence in plants are poorly understood. In the present study, comprehensive bioinformatics analyses were performed to gain insight into the molecular and evolutionary mechanisms driving this structural and functional diversification. Phylogenetic analyses have revealed the major lineage-specific expansions of XIP genes in plants. Within the eudicots, XIP genes have diverged into Asterid and Rosid-specific phylogenetic lineages and have also undergone several independent duplications during the course of evolution. Investigation of functional divergence at the protein level showed evidence for shifting evolutionary rate and/or altered constraints on the physiochemical properties of specific amino acid sites following gene duplication. Selection pressure analyses suggest that purifying selection is the predominant evolutionary force acting on the XIP gene subfamily, along with episodic positive selection. However, only a few amino acid sites were found to be subjected to such episodic positive selection. Furthermore, protein functional divergence analysis has identified critical amino acid residues, which must be validated by future experimental studies, that could provide new insights into the role of XIPs in transport of a wide range solutes of physiological importance.

The complex intron landscape and massive intron invasion in a picoeukaryote provides insights into intron evolution

Genes in pieces and spliceosomal introns are a landmark of eukaryotes, with intron invasion usually assumed to have happened early on in evolution. Here, we analyze the intron landscape of Micromonas, a unicellular green alga in the Mamiellophyceae lineage, demonstrating the coexistence of several classes of introns and the occurrence of recent massive intron invasion. This study focuses on two strains, CCMP1545 and RCC299, and their related individuals from ocean samplings, showing that they not only harbor different classes of introns depending on their location in the genome, as for other Mamiellophyceae, but also uniquely carry several classes of repeat introns. These introns, dubbed introner elements (IEs), are found at novel positions in genes and have conserved sequences, contrary to canonical introns. This IE invasion has a huge impact on the genome, doubling the number of introns in the CCMP1545 strain. We hypothesize that each IE class originated from a single ancestral IE that has been colonizing the genome after strain divergence by inserting copies of itself into genes by intron transposition, likely involving reverse splicing. Along with similar cases recently observed in other organisms, our observations in Micromonas strains shed a new light on the evolution of introns, suggesting that intron gain is more widespread than previously thought.

ESR1 rs9340799 is associated with endometriosis-related infertility and in vitro fertilization failure

Estrogen receptor alpha has a central role in human fertility by regulating estrogen action in all human reproductive tissues. Leukemia inhibitory factor (LIF) expression, a cytokine critical for blastocyst implantation, is mediated by estrogen signaling, so we hypothesized that ESR1 gene polymorphisms might be candidate risk markers for endometriosis-related infertility and in vitro fertilization (IVF) failure. We included 98 infertile women with endometriosis, 115 infertile women with at least one IVF failure and also 134 fertile women as controls. TaqMan SNP assays were used for genotyping LIF (rs929271), MDM2 (rs2279744), MDM4 (rs1563828), USP7 (rs1529916), and ESR1 (rs9340799 and rs2234693) polymorphisms. The SNP ESR1 rs9340799 was associated with endometriosis-related infertility (P < 0.001) and also with IVF failure (P = 0.018). After controlling for age, infertile women with ESR1 rs9340799 GG genotype presented 4-fold increased risk of endometriosis (OR 4.67, 95% CI 1.84–11.83, P = 0.001) and 3-fold increased risk of IVF failure (OR 3.33, 95% CI 1.38–8.03, P = 0.007). Our results demonstrate an association between ESR1 rs9340799 polymorphism and infertile women with endometriosis and also with women who were submitted to IVF procedures and had no blastocyst implantation.

Genetic variants and evolutionary analyses of heparin cofactor II

Heparin cofactor II (HCII) belongs to serpin superfamily and it acts as a thrombin inhibitor in the coagulation cascade, in a glycosaminoglycan-dependent pathway using the release of a sequestered hirudin-like N-terminal tail for interaction with thrombin. This serpin belongs to multiple member group V2 of vertebrate serpin classification. However, there is no comprehensive study illustrating the exact phylogenetic history of HCII, to date. Herein, we explored phylogenetic traits of HCII genes. Structures of HCII gene from selected ray-finned fishes and lamprey varied in exon I and II with insertions of novel introns of which one in core domain for ray-finned fishes in exon II at the position 241c. We found HCII remain nested in the largest intron of phosphatidylinositol (PI) 4-kinase (PIK4CA) gene (genetic variants of this gene cause schizophrenia) at the origin of vertebrates, dated about 500MY old. We found that sequence features such as two acidic repeats (AR1-II), GAG-binding helix-D, three serpin motifs and inhibitory reactive center loop (RCL) of HCII protein are highly conserved in 55 vertebrates analyzed. We identified 985 HCII variants by analysis of 1092 human genomes with top three variation classes belongs to SNPs (84.3%), insertion (7.1%) and deletion (5.0%). We identified 37 deleterious mutations in the human HCII protein and we have described these mutations in relation to HCII sequence-structure-function relationships. These understandings may have clinical and medical importance as well.

Distribution of conventional and nonconventional introns in tubulin (α and β) genes of euglenids

The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB, from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns--positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism--intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns--because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

Impact of an intragenic retrotransposon on the structural integrity and evolution of a major isoprenoid biosynthesis pathway gene in Hevea brasiliensis

Isoprenoids belong to a large family of structurally and functionally different natural compounds found universally from prokaryotes to higher animals and plants. In Hevea brasiliensis, the commercially important cis-polyisoprene (rubber) is synthesised as part of its defence mechanism in addition to other common isoprenoids like phytosterols, growth hormones etc. Farnesyl diphosphate synthase (FDPS) is a key enzyme in this process which catalyses the conversion of isoprene units into polyisoprene. Although prior sequence information is available, the structural variants of the FDPS gene presently existing in Hevea population are largely unknown. Since gene structure has a major role in gene regulation, extensive sequence analysis of this gene from different genotypes was carried out to identify the prevailing structural variants. We identified several SNPs and large indels which were associated with a partial transposable element (TE). Modification of key regulatory motifs and splice sites induced by the retroelement was also identified in the first intron. Screening of popular rubber clones, wild germplasm accessions and Hevea species revealed that the retroelement is responsible for the generation of new alleles with varying degrees of sequence homology. Segregation analysis of a progeny population confirmed that the alleles are not paralogs and are inherited in a Mendelian mode. Our findings suggest that the first intron of the FDPS gene has been subjected to various chromosomal rearrangements due to the interaction of a retrotransposon, resulting in novel alleles which may substantially contribute towards the evolution of this major gene in rubber. Moreover, the results indicate the possible existence of a retrotransposon-mediated epigenetic gene regulatory mechanism in Hevea.