The origin of introns and their role in eukaryogenesis: a compromise solution to the introns-early versus introns-late debate?

Study of the RNA splicing kinetics via in vivo 5-EU labeling

Splicing is an important step of gene expression in all eukaryotes. Splice sites might be used with different efficiency, giving rise to alternative splicing products. At the same time, splice sites might be used at a variable rate. We used 5-ethynyl uridine labeling to sequence a nascent transcriptome of HeLa cells and deduced the rate of splicing for each donor and acceptor splice site. The following correlation analysis showed a correspondence of primary transcript features with the rate of splicing. Some dependencies we revealed were anticipated, such as a splicing rate decrease with a decreased complementarity of the donor splice site to U1 and acceptor sites to U2 snRNAs. Other dependencies were more surprising, like a negative influence of a distance to the 5' end on the rate of the acceptor splicing site utilization, or the differences in splicing rate between long, short, and RBM17-dependent introns. We also observed a deceleration of last intron splicing with an increase of the distance to the poly(A) site, which might be explained by the cooperativity of the splicing and polyadenylation. Additional analysis of splicing kinetics of SF3B4 knockdown cells suggested the impairment of a U2 snRNA recognition step. As a result, we deconvoluted the effects of several examined features on the splicing rate into a single regression model. The data obtained here are useful for further studies in the field, as they provide general splicing rate dependencies as well as help to justify the existence of slowly removed splice sites.

Evolution and Comparative Genomics of the Transforming Growth Factor-β-Related Proteins in Nile Tilapia

The members of the transforming growth factor β (TGF-β) family of cell signaling polypeptides have garnered a great deal of interest due to its capacity from nematodes to mammals to regulate cell-based activities which control the growth of embryos and sustain tissue homeostasis. The current study designed a computational analysis of the TGF-β protein family for understanding these proteins at the molecular level. This study determined the genomic structure of TGF-β gene family in Nile tilapia for the first time. We chose 33 TGF-β genes for identification and divided them into two subgroups, TGF-like and BMP-like. Moreover, the subcellular localization of the Nile tilapia TGF-β proteins have showed that majority of the members of TGF-β proteins family are present into extracellular matrix and plasma except BMP6, BMP7, and INHAC. All TGF-β proteins were thermostable excluding BMP1. Each protein exhibited basic nature, excluding of BMP1, BMP2, BMP7, BMP10, GDF2, GDF8, GDF11, AMH, INHA, INHBB, and NODAL M. All proteins gave impression of being unstable depending on the instability index, having values exceeding 40 excluding BMP1 and BMP2. Each TGF-β protein was found to be hydrophobic with lowered values of GRAVY. Moreover, every single one of the discovered TGF-β genes had a consistent evolutionary pattern. The TGF-β gene family had eight segmental duplications, and the Ka/Ks ratio demonstrated that purifying selection had an impact on the duplicated gene pairs which have experienced selection pressure. This study highlights important functionality of TGF-β and depicts the demand for further investigation to better understand the role and mechanism of transforming growth factor β in fishes and other species.

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

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.

Copper-binding proteins and exonic splicing enhancers and silencers

Eukaryotic DNA codes not only for proteins but contains a wealth of information required for accurate splicing of messenger RNA precursors and inclusion of constitutively or alternatively spliced exons in mature transcripts. This "auxiliary" splicing code has been characterized as exonic splicing enhancers and silencers (ESE and ESS). The exact interplay between protein and splicing codes is, however, poorly understood. Here, we show that exons encoding copper-coordinating amino acids in human cuproproteins lack ESEs and/or have an excess of ESSs, yet RNA sequencing and expressed sequence tags data show that they are more efficiently included in mature transcripts by the splicing machinery than average exons. Their largely constitutive inclusion in messenger RNA is facilitated by stronger splice sites, including polypyrimidine tracts, consistent with an important role of the surrounding intron architecture in ensuring high expression of metal-binding residues during evolution. ESE/ESS profiles of codons and entire exons that code for copper-coordinating residues were very similar to those encoding residues that coordinate zinc but markedly different from those that coordinate calcium. Together, these results reveal how the traditional and auxiliary splicing motifs responded to constraints of metal coordination in proteins.

Comparative transcriptional analysis of Persea americana MYB, WRKY and AP2/ERF transcription factors following Phytophthora cinnamomi infection

Plant cells undergo extensive transcriptional reprogramming following pathogen infection, with these reprogramming patterns becoming more complex when pathogens, such as hemibiotrophs, exhibit different lifestyles. These transcriptional changes are often orchestrated by MYB, WRKY and AP2/ERF transcription factors (TFs), which modulate both growth and defence-related gene expression. Transcriptional analysis of defence-related genes in avocado (Persea americana) infected with Phytophthora cinnamomi indicated differential immune response activation when comparing a partially resistant and susceptible rootstock. This study identified 226 MYB, 82 WRKY, and 174 AP2/ERF TF-encoding genes in avocado, using a genome-wide approach. Phylogenetic analysis revealed substantial sequence conservation within TF groups underscoring their functional significance. RNA-sequencing analysis in a partially resistant and susceptible avocado rootstock infected with P. cinnamomi was indicative of an immune response switch occurring in either rootstock after 24 and 6 h post-inoculation, respectively. Different clusters of co-expressed TF genes were observed at these times, suggesting the activation of necrotroph-related immune responses at varying intervals between the two rootstocks. This study aids our understanding of avocado immune response activation following P. cinnamomi infection, and the role of the TFs therein, elucidating the transcriptional reprogramming disparities between partially resistant and susceptible rootstocks.

On the origin of the nucleus: a hypothesis

SUMMARYIn this hypothesis article, we explore the origin of the eukaryotic nucleus. In doing so, we first look afresh at the nature of this defining feature of the eukaryotic cell and its core functions-emphasizing the utility of seeing the eukaryotic nucleoplasm and cytoplasm as distinct regions of a common compartment. We then discuss recent progress in understanding the evolution of the eukaryotic cell from archaeal and bacterial ancestors, focusing on phylogenetic and experimental data which have revealed that many eukaryotic machines with nuclear activities have archaeal counterparts. In addition, we review the literature describing the cell biology of representatives of the TACK and Asgardarchaeaota - the closest known living archaeal relatives of eukaryotes. Finally, bringing these strands together, we propose a model for the archaeal origin of the nucleus that explains much of the current data, including predictions that can be used to put the model to the test.

Proteome-wide structural analysis quantifies structural conservation across distant species

Traditional evolutionary biology research mainly relies on sequence information to infer evolutionary relationships between genes or proteins. In contrast, protein structural information has long been overlooked, although structures are more conserved and closely linked to the functions than the sequences. To address this gap, we conducted a proteome-wide structural analysis using experimental and computed protein structures for organisms from the three distinct domains, including Homo sapiens (eukarya), Escherichia coli (bacteria), and Methanocaldococcus jannaschii (archaea). We reveal the distribution of structural similarity and sequence identity at the genomic level and characterize the twilight zone, where signals obtained from sequence alignment are blurred and evolutionary relationships cannot be inferred unambiguously. We find that structurally similar homologous protein pairs in the twilight zone account for ∼0.004%-0.021% of all possible protein pair combinations, which translates to ∼8%-32% of the protein-coding genes, depending on the species under comparison. In addition, by comparing the structural homologs, we show that human proteins involved in the energy supply are more similar to their E. coli homologs, whereas proteins relating to the central dogma are more similar to their M. jannaschii homologs. We also identify a bacterial GPCR homolog in the E. coli proteome that displays distinctive domain architecture. Our results shed light on the characteristics of the twilight zone and the origin of different pathways from a protein structure perspective, highlighting an exciting new frontier in evolutionary biology.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

Spliceosome assembly and regulation: insights from analysis of highly reduced spliceosomes

Premessenger RNA splicing is catalyzed by the spliceosome, a multimegadalton RNA-protein complex that assembles in a highly regulated process on each intronic substrate. Most studies of splicing and spliceosomes have been carried out in human or S. cerevisiae model systems. There exists, however, a large diversity of spliceosomes, particularly in organisms with reduced genomes, that suggests a means of analyzing the essential elements of spliceosome assembly and regulation. In this review, we characterize changes in spliceosome composition across phyla, describing those that are most frequently observed and highlighting an analysis of the reduced spliceosome of the red alga Cyanidioschyzon merolae We used homology modeling to predict what effect splicing protein loss would have on the spliceosome, based on currently available cryo-EM structures. We observe strongly correlated loss of proteins that function in the same process, for example, in interacting with the U1 snRNP (which is absent in C. merolae), regulation of Brr2, or coupling transcription and splicing. Based on our observations, we predict splicing in C. merolae to be inefficient, inaccurate, and post-transcriptional, consistent with the apparent trend toward its elimination in this lineage. This work highlights the striking flexibility of the splicing pathway and the spliceosome when viewed in the context of eukaryotic diversity.

Transposon-encoded nucleases use guide RNAs to selfishly bias their inheritance

Insertion sequences (IS) are compact and pervasive transposable elements found in bacteria, which encode only the genes necessary for their mobilization and maintenance. IS 200 /IS 605 elements undergo 'peel-and-paste' transposition catalyzed by a TnpA transposase, but intriguingly, they also encode diverse, TnpB- and IscB-family proteins that are evolutionarily related to the CRISPR-associated effectors Cas12 and Cas9, respectively. Recent studies demonstrated that TnpB-family enzymes function as RNA-guided DNA endonucleases, but the broader biological role of this activity has remained enigmatic. Here we show that TnpB/IscB are essential to prevent permanent transposon loss as a consequence of the TnpA transposition mechanism. We selected a family of related IS elements from Geobacillus stearothermophilus that encode diverse TnpB/IscB orthologs, and showed that a single TnpA transposase was active for transposon excision. The donor joints formed upon religation of IS-flanking sequences were efficiently targeted for cleavage by RNA-guided TnpB/IscB nucleases, and co-expression of TnpB together with TnpA led to significantly greater transposon retention, relative to conditions in which TnpA was expressed alone. Remarkably, TnpA and TnpB/IscB recognize the same AT-rich transposon-adjacent motif (TAM) during transposon excision and RNA-guided DNA cleavage, respectively, revealing a striking convergence in the evolution of DNA sequence specificity between collaborating transposase and nuclease proteins. Collectively, our study reveals that RNA-guided DNA cleavage is a primal biochemical activity that arose to bias the selfish inheritance and spread of transposable elements, which was later co-opted during the evolution of CRISPR-Cas adaptive immunity for antiviral defense.

Recognition of transcription terminators during retrotransposition: How to keep a group II intron quiet

Chung et al. recently presented the structure of a primitive group IIC intron with its DNA target, which reveals the structural requirements that this class of intron uses to recognize a transcription terminator stem loop at the DNA level for insertion during retrotransposition.

Casposons - silent heroes of the CRISPR-Cas systems evolutionary history

Many archaeal and bacterial organisms possess an adaptive immunity system known as CRISPR-Cas. Its role is to recognize and degrade foreign DNA showing high similarity to repeats within the CRISPR array. In recent years computational techniques have been used to identify cas1 genes that are not associated with CRISPR systems, named cas1-solo. Often, cas1-solo genes are present in a conserved neighborhood of PolB-like polymerase genes, which is a characteristic feature of self-synthesizing, eukaryotic transposons of the Polinton class. Nearly all cas1-polB genomic islands are flanked by terminal inverted repeats and direct repeats which correspond to target site duplications. Considering the patchy taxonomic distribution of the identified islands in archaeal and bacterial genomes, they were characterized as a new superfamily of mobile genetic elements and called casposons. Here, we review recent experiments on casposons' mobility and discuss their discovery, classification, and evolutionary relationship with the CRISPR-Cas systems.

Mitochondrial Dynamics during Development

Mitochondria are dynamic membrane-bound organelles in eukaryotic cells. These are important for the generation of chemical energy needed to power various cellular functions and also support metabolic, energetic, and epigenetic regulation in various cells. These organelles are also important for communication with the nucleus and other cellular structures, to maintain developmental sequences and somatic homeostasis, and for cellular adaptation to stress. Increasing information shows mitochondrial defects as an important cause of inherited disorders in different organ systems. In this article, we provide an extensive review of ontogeny, ultrastructural morphology, biogenesis, functional dynamics, important clinical manifestations of mitochondrial dysfunction, and possibilities for clinical intervention. We present information from our own clinical and laboratory research in conjunction with information collected from an extensive search in the databases PubMed, EMBASE, and Scopus.

The evolution of pre-mRNA splicing and its machinery revealed by reduced extremophilic red algae

The Cyanidiales are a group of mostly thermophilic and acidophilic red algae that thrive near volcanic vents. Despite their phylogenetic relationship, the reduced genomes of Cyanidioschyzon merolae and Galdieria sulphuraria are strikingly different with respect to pre-mRNA splicing, a ubiquitous eukaryotic feature. Introns are rare and spliceosomal machinery is extremely reduced in C. merolae, in contrast to G. sulphuraria. Previous studies also revealed divergent spliceosomes in the mesophilic red alga Porphyridium purpureum and the red algal derived plastid of Guillardia theta (Cryptophyta), along with unusually high levels of unspliced transcripts. To further examine the evolution of splicing in red algae, we compared C. merolae and G. sulphuraria, investigating splicing levels, intron position, intron sequence features, and the composition of the spliceosome. In addition to identifying 11 additional introns in C. merolae, our transcriptomic analysis also revealed typical eukaryotic splicing in G. sulphuraria, whereas most transcripts in C. merolae remain unspliced. The distribution of intron positions within their host genes was examined to provide insight into patterns of intron loss in red algae. We observed increasing variability of 5' splice sites and branch donor regions with increasing intron richness. We also found these relationships to be connected to reductions in and losses of corresponding parts of the spliceosome. Our findings highlight patterns of intron and spliceosome evolution in related red algae under the pressures of genome reduction.

Exon Elongation Added Intrinsically Disordered Regions to the Encoded Proteins and Facilitated the Emergence of the Last Eukaryotic Common Ancestor

Most prokaryotic proteins consist of a single structural domain (SD) with little intrinsically disordered regions (IDRs) that by themselves do not adopt stable structures, whereas the typical eukaryotic protein comprises multiple SDs and IDRs. How eukaryotic proteins evolved to differ from prokaryotic proteins has not been fully elucidated. Here, we found that the longer the internal exons are, the more frequently they encode IDRs in eight eukaryotes including vertebrates, invertebrates, a fungus, and plants. Based on this observation, we propose the "small bang" model from the proteomic viewpoint: the protoeukaryotic genes had no introns and mostly encoded one SD each, but a majority of them were subsequently divided into multiple exons (step 1). Many exons unconstrained by SDs elongated to encode IDRs (step 2). The elongated exons encoding IDRs frequently facilitated the acquisition of multiple SDs to make the last common ancestor of eukaryotes (step 3). One prediction of the model is that long internal exons are mostly unconstrained exons. Analytical results of the eight eukaryotes are consistent with this prediction. In support of the model, we identified cases of internal exons that elongated after the rat-mouse divergence and discovered that the expanded sections are mostly in unconstrained exons and preferentially encode IDRs. The model also predicts that SDs followed by long internal exons tend to have other SDs downstream. This prediction was also verified in all the eukaryotic species analyzed. Our model accounts for the dichotomy between prokaryotic and eukaryotic proteins and proposes a selective advantage conferred by IDRs.

Repeat sequences limit the effectiveness of lateral gene transfer and favored the evolution of meiotic sex in early eukaryotes

The transition from prokaryotic lateral gene transfer to eukaryotic meiotic sex is poorly understood. Phylogenetic evidence suggests that it was tightly linked to eukaryogenesis, which involved an unprecedented rise in both genome size and the density of genetic repeats. Expansion of genome size raised the severity of Muller's ratchet, while limiting the effectiveness of lateral gene transfer (LGT) at purging deleterious mutations. In principle, an increase in recombination length combined with higher rates of LGT could solve this problem. Here, we show using a computational model that this solution fails in the presence of genetic repeats prevalent in early eukaryotes. The model demonstrates that dispersed repeat sequences allow ectopic recombination, which leads to the loss of genetic information and curtails the capacity of LGT to prevent mutation accumulation. Increasing recombination length in the presence of repeat sequences exacerbates the problem. Mutational decay can only be resisted with homology along extended sequences of DNA. We conclude that the transition to homologous pairing along linear chromosomes was a key innovation in meiotic sex, which was instrumental in the expansion of eukaryotic genomes and morphological complexity.

Imaging Intron Evolution

Intron evolution may be readily imaged through the combined use of the “dot plot” function of the NCBI BLAST, aligning two sequences at a time, and the Vertebrate “Multiz” alignment and conservation tool of the UCSC Genome Browser. With the NCBI BLAST, an ideal alignment of two highly conserved sequences generates a diagonal straight line in the plot from the lower left corner to the upper right corner. Gaps in this line correspond to non-conserved sections. In addition, the dot plot of the alignment of a sequence with the same sequence after the removal of the Transposable Elements (TEs) can be observed along the diagonal gaps that correspond to the sites of TE insertion. The UCSC Genome Browser can graph, along the entire sequence of a single gene, the level of overall conservation in vertebrates. This level can be compared with the conservation level of the gene in one or more selected vertebrate species. As an example, we show the graphic analysis of the intron conservation in two genes: the mitochondrial solute carrier 21 (SLC25A21) and the growth hormone receptor (GHR), whose coding sequences are conserved through vertebrates, while their introns show dramatic changes in nucleotide composition and even length. In the SLC25A21, a few short but significant nucleotide sequences are conserved in zebrafish, Xenopus and humans, and the rate of conservation steadily increases from chicken/human to mouse/human alignments. In the GHR, a less conserved gene, the earlier indication of intron conservation is a small signal in chicken/human alignment. The UCSC tool may simultaneously display the conservation level of a gene in different vertebrates, with reference to the level of overall conservation in Vertebrates. It is shown that, at least in SLC25A21, the sites of higher conservation are not always coincident in chicken and zebrafish nor are the sites of higher vertebrate conservation.

Can Introns Stabilize Gene Duplication?

Gene duplication is considered one of the most important events that determine the evolution of genomes. However, the neo-duplication condition of a given gene is particularly unstable due to recombination events. Several mechanisms have been proposed to justify this step. In this "opinion article" we propose a role for intron sequences in stabilizing gene duplication by limiting and reducing the identity of the gene sequence between the two duplicated copies. A review of the topic and a detailed hypothesis are presented.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

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.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Independent Evolution of the MYB Family in Brown Algae

Myeloblastosis (MYB) proteins represent one of the largest families of eukaryotic transcription factors and regulate important processes in growth and development. Studies on MYBs have mainly focused on animals and plants; however, comprehensive analysis across other supergroups such as SAR (stramenopiles, alveolates, and rhizarians) is lacking. This study characterized the structure, evolution, and expression of MYBs in four brown algae, which comprise the biggest multicellular lineage of SAR. Subfamily 1R-MYB comprised heterogeneous proteins, with fewer conserved motifs found outside the MYB domain. Unlike the SHAQKY subgroup of plant 1R-MYB, THAQKY comprised the largest subgroup of brown algal 1R-MYBs. Unlike the expansion of 2R-MYBs in plants, brown algae harbored more 3R-MYBs than 2R-MYBs. At least ten 2R-MYBs, fifteen 3R-MYBs, and one 6R-MYB orthologs existed in the common ancestor of brown algae. Phylogenetic analysis showed that brown algal MYBs had ancient origins and a diverged evolution. They showed strong affinity with stramenopile species, while not with red algae, green algae, or animals, suggesting that brown algal MYBs did not come from the secondary endosymbiosis of red and green plastids. Sequence comparison among all repeats of the three types of MYB subfamilies revealed that the repeat of 1R-MYBs showed higher sequence identity with the R3 of 2R-MYBs and 3R-MYBs, which supports the idea that 1R-MYB was derived from loss of the first and second repeats of the ancestor MYB. Compared with other species of SAR, brown algal MYB proteins exhibited a higher proportion of intrinsic disordered regions, which might contribute to multicellular evolution. Expression analysis showed that many MYB genes are responsive to different stress conditions and developmental stages. The evolution and expression analyses provided a comprehensive analysis of the phylogeny and functions of MYBs in brown algae.

Genome-Wide Identification, Characterization and Function Analysis of Lineage-Specific Genes in the Tea Plant Camellia sinensis

Genes that have no homologous sequences with other species are called lineage-specific genes (LSGs), are common in living organisms, and have an important role in the generation of new functions, adaptive evolution and phenotypic alteration of species. Camellia sinensis var. sinensis (CSS) is one of the most widely distributed cultivars for quality green tea production. The rich catechins in tea have antioxidant, free radical elimination, fat loss and cancer prevention potential. To further understand the evolution and utilize the function of LSGs in tea, we performed a comparative genomics approach to identify Camellia-specific genes (CSGs). Our result reveals that 1701 CSGs were identified specific to CSS, accounting for 3.37% of all protein-coding genes. The majority of CSGs (57.08%) were generated by gene duplication, and the time of duplication occurrence coincide with the time of two genome-wide replication (WGD) events that happened in CSS genome. Gene structure analysis revealed that CSGs have shorter gene lengths, fewer exons, higher GC content and higher isoelectric point. Gene expression analysis showed that CSG had more tissue-specific expression compared to evolutionary conserved genes (ECs). Weighted gene co-expression network analysis (WGCNA) showed that 18 CSGs are mainly associated with catechin synthesis-related pathways, including phenylalanine biosynthesis, biosynthesis of amino acids, pentose phosphate pathway, photosynthesis and carbon metabolism. Besides, we found that the expression of three CSGs (CSS0030246, CSS0002298, and CSS0030939) was significantly down-regulated in response to both types of stresses (salt and drought). Our study first systematically identified LSGs in CSS, and comprehensively analyzed the features and potential functions of CSGs. We also identified key candidate genes, which will provide valuable assistance for further studies on catechin synthesis and provide a molecular basis for the excavation of excellent germplasm resources.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

Evolution of plant telomerase RNAs: farther to the past, deeper to the roots

The enormous sequence heterogeneity of telomerase RNA (TR) subunits has thus far complicated their characterization in a wider phylogenetic range. Our recent finding that land plant TRs are, similarly to known ciliate TRs, transcribed by RNA polymerase III and under the control of the type-3 promoter, allowed us to design a novel strategy to characterize TRs in early diverging Viridiplantae taxa, as well as in ciliates and other Diaphoretickes lineages. Starting with the characterization of the upstream sequence element of the type 3 promoter that is conserved in a number of small nuclear RNAs, and the expected minimum TR template region as search features, we identified candidate TRs in selected Diaphoretickes genomes. Homologous TRs were then used to build covariance models to identify TRs in more distant species. Transcripts of the identified TRs were confirmed by transcriptomic data, RT-PCR and Northern hybridization. A templating role for one of our candidates was validated in Physcomitrium patens. Analysis of secondary structure demonstrated a deep conservation of motifs (pseudoknot and template boundary element) observed in all published TRs. These results elucidate the evolution of the earliest eukaryotic TRs, linking the common origin of TRs across Diaphoretickes, and underlying evolutionary transitions in telomere repeats.

An Update on Trichoderma Mitogenomes: Complete De Novo Mitochondrial Genome of the Fungal Biocontrol Agent Trichoderma harzianum (Hypocreales, Sordariomycetes), an Ex-Neotype Strain CBS 226.95, and Tracing the Evolutionary Divergences of Mitogenomes in Trichoderma

Members of the genus Trichoderma (Hypocreales), widely used as biofungicides, biofertilizers, and as model fungi for the industrial production of CAZymes, have actively been studied for the applications of their biological functions. Recently, the study of the nuclear genomes of Trichoderma has expanded in the directions of adaptation and evolution to gain a better understanding of their ecological traits. However, Trichoderma's mitochondria have received much less attention despite mitochondria being the most necessary element for sustaining cell life. In this study, a mitogenome of the fungus Trichoderma harzianum CBS 226.95 was assembled de novo. A 27,632 bp circular DNA molecule was revealed with specific features, such as the intronless of all core PCGs, one homing endonuclease, and a putative overlapping tRNA, on a closer phylogenetic relationship with T. reesei among hypocrealean fungi. Interestingly, the mitogenome of T. harzianum CBS 226.95 was predicted to have evolved earlier than those of other Trichoderma species and also assumed with a selection pressure in the cox3. Considering the bioavailability, both for the ex-neotype strain of the T. harzianum species complex and the most globally representative commercial fungal biocontrol agent, our results on the T. harzianum CBS 226.95 mitogenome provide crucial information which will be helpful criteria in future studies on Trichoderma.

The Coevolution of Fungal Mitochondrial Introns and Their Homing Endonucleases (GIY-YIG and LAGLIDADG)

Fungal mitochondrial (mt) genomes exhibit great diversity in size which is partially attributed to their variable intergenic regions and most importantly to the inclusion of introns within their genes. These introns belong to group I or II, and both of them are self-splicing. The majority of them carry genes encoding homing endonucleases, either LAGLIDADG or GIY-YIG. In this study, it was found that these intronic homing endonucleases genes (HEGs) may originate from mt free-standing open reading frames which can be found nowadays in species belonging to Early Diverging Fungi as “living fossils.” A total of 487 introns carrying HEGs which were located in the publicly available mt genomes of representative species belonging to orders from all fungal phyla was analyzed. Their distribution in the mt genes, their insertion target sequence, and the phylogenetic analyses of the HEGs showed that these introns along with their HEGs form a composite structure in which both selfish elements coevolved. The invasion of the ancestral free-standing HEGs in the introns occurred through a perpetual mechanism, called in this study as “aenaon” hypothesis. It is based on recombination, transpositions, and horizontal gene transfer events throughout evolution. HEGs phylogenetically clustered primarily according to their intron hosts and secondarily to the mt genes carrying the introns and their HEGs. The evolutionary models created revealed an “intron-early” evolution which was enriched by “intron-late” events through many different independent recombinational events which resulted from both vertical and horizontal gene transfers.

U5 snRNA Interactions With Exons Ensure Splicing Precision

Imperfect conservation of human pre-mRNA splice sites is necessary to produce alternative isoforms. This flexibility is combined with the precision of the message reading frame. Apart from intron-termini GU_AG and the branchpoint A, the most conserved are the exon-end guanine and +5G of the intron start. Association between these guanines cannot be explained solely by base-pairing with U1 snRNA in the early spliceosome complex. U6 succeeds U1 and pairs +5G in the pre-catalytic spliceosome, while U5 binds the exon end. Current U5 snRNA reconstructions by CryoEM cannot explain the conservation of the exon-end G. Conversely, human mutation analyses show that guanines of both exon termini can suppress splicing mutations. Our U5 hypothesis explains the mechanism of splicing precision and the role of these conserved guanines in the pre-catalytic spliceosome. We propose: (1) optimal binding register for human exons and U5-the exon junction positioned at U5Loop1 C₃₉|C₃₈; (2) common mechanism for base-pairing of human U5 snRNA with diverse exons and bacterial Ll.LtrB intron with new loci in retrotransposition-guided by base pair geometry; and (3) U5 plays a significant role in specific exon recognition in the pre-catalytic spliceosome. Statistical analyses showed increased U5 Watson-Crick pairs with the 5'exon in the absence of +5G at the intron start. In 5'exon positions -3 and -5, this effect is specific to U5 snRNA rather than U1 snRNA of the early spliceosome. Increased U5 Watson-Crick pairs with 3'exon position +1 coincide with substitutions of the conserved -3C at the intron 3'end. Based on mutation and X-ray evidence, we propose that -3C pairs with U2 G₃₁ juxtaposing the branchpoint and the 3'intron end. The intron-termini pair, formed in the pre-catalytic spliceosome to be ready for transition after branching, and the early involvement of the 3'intron end ensure that the 3'exon contacts U5 in the pre-catalytic complex. We suggest that splicing precision is safeguarded cooperatively by U5, U6, and U2 snRNAs that stabilize the pre-catalytic complex by Watson-Crick base pairing. In addition, our new U5 model explains the splicing effect of exon-start +1G mutations: U5 Watson-Crick pairs with exon +2C/+3G strongly promote exon inclusion. We discuss potential applications for snRNA therapeutics and gene repair by reverse splicing.

Extraordinary diversity of telomeres, telomerase RNAs and their template regions in Saccharomycetaceae

Telomerase RNA (TR) carries the template for synthesis of telomere DNA and provides a scaffold for telomerase assembly. Fungal TRs are long and have been compared to higher eukaryotes, where they show considerable diversity within phylogenetically close groups. TRs of several Saccharomycetaceae were recently identified, however, many of these remained uncharacterised in the template region. Here we show that this is mainly due to high variability in telomere sequence. We predicted the telomere sequences using Tandem Repeats Finder and then we identified corresponding putative template regions in TR candidates. Remarkably long telomere units and the corresponding putative TRs were found in Tetrapisispora species. Notably, variable lengths of the annealing sequence of the template region (1–10 nt) were found. Consequently, species with the same telomere sequence may not harbour identical TR templates. Thus, TR sequence alone can be used to predict a template region and telomere sequence, but not to determine these exactly. A conserved feature of telomere sequences, tracts of adjacent Gs, led us to test the propensity of individual telomere sequences to form G4. The results show highly diverse values of G4-propensity, indicating the lack of ubiquitous conservation of this feature across Saccharomycetaceae.

Analysis of Different Size Fractions Provides a More Complete Perspective of Viral Diversity in a Freshwater Embayment

Inspired by recent discoveries of the prevalence of large viruses in the environment, we reassessed the longstanding approach of filtering water through small-pore-size filters to separate viruses from cells before metagenomic analysis. We collected samples from three sites in Hamilton Harbour, an embayment of Lake Ontario, and studied 6 data sets derived from <0.45-μm- and >0.45-μm-size fractions to compare the diversity of viruses in these fractions. At the level of virus order/family, we observed highly diverse and distinct virus communities in the >0.45-μm-size fractions, whereas the <0.45-μm-size fractions were composed primarily of Caudovirales The relative abundances of Caudovirales for which hosts could be inferred varied widely between size fractions, with higher relative abundances of cyanophages in the >0.45-μm-size fractions, potentially indicating replication within cells during ongoing infections. Many viruses of eukaryotes, such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, were detected exclusively in the often-disregarded >0.45-μm-size fractions. In addition to observing unique virus communities associated with each size fraction from every site we examined, we detected viruses common to both fractions, suggesting that these are candidates for further exploration because they could be the product of ongoing or recent lytic events. Most importantly, our observations indicate that analysis of either fraction alone provides only a partial perspective of double-stranded DNA (dsDNA) viruses in the environment, highlighting the need for more comprehensive approaches for analyzing virus communities inferred from metagenomic sequencing.IMPORTANCE Most studies of aquatic virus communities analyze DNA sequences derived from the smaller-size "free-virus" fraction. Our study demonstrates that analysis of virus communities using only the smaller-size fraction can lead to erroneously low diversity estimates for many of the larger viruses such as Mimiviridae, Phycodnaviridae, Iridoviridae, and Poxviridae, whereas analyzing only the larger->0.45-μm-size fraction can lead to underestimates of Caudovirales diversity and relative abundance. Similarly, our data show that examining only the smaller-size fraction can lead to underestimations of virophage and cyanophage relative abundances that could, in turn, cause researchers to assume their limited ecological importance. Given the considerable differences we observed in this study, we recommend cautious interpretations of environmental virus community assemblages and dynamics when based on metagenomic data derived from different size fractions.

Mutations in cis that affect mRNA synthesis, processing and translation

Genetic mutations that cause hereditary diseases usually affect the composition of the transcribed mRNA and its encoded protein, leading to instability of the mRNA and/or the protein. Sometimes, however, such mutations affect the synthesis, the processing or the translation of the mRNA, with similar disastrous effects. We here present an overview of mRNA synthesis, its posttranscriptional modification and its translation into protein. We then indicate which elements in these processes are known to be affected by pathogenic mutations, but we restrict our review to mutations in cis, in the DNA of the gene that encodes the affected protein. These mutations can be in enhancer or promoter regions of the gene, which act as binding sites for transcription factors involved in pre-mRNA synthesis. We also describe mutations in polyadenylation sequences and in splice site regions, exonic and intronic, involved in intron removal. Finally, we include mutations in the Kozak sequence in mRNA, which is involved in protein synthesis. We provide examples of genetic diseases caused by mutations in these DNA regions and refer to databases to help identify these regions. The over-all knowledge of mRNA synthesis, processing and translation is essential for improvement of the diagnosis of patients with genetic diseases.

An Intron of Invertebrate Microphthalmia Transcription Factor Gene Is Evolved from a Longer Ancestral Sequence

Introns are highly variable in number and size. Sequence simulation is an effective method to elucidate intron evolution patterns. Previously, we have reported that introns are more likely to evolve through mutation-and-deletion (MD) rather than through mutation-and-insertion (MI). In the present study, we further studied evolution models by allowing insertion in the MD model and by allowing deletion in the MI model at various frequencies. It was found that all deletion-biased models with proper parameter settings could generate sequences with attributes matchable to 16 invertebrate introns from the microphthalmia transcription factor gene, whereas all insertion-biased models with any parameter settings failed to generate such sequences. We conclude that the examined invertebrate introns may have evolved from a longer ancestral sequence in a deletion-biased pattern. The constructed models are useful for studying the evolution of introns from other genes and/or from other taxonomic groups. (C++ scripts of all deletion- and insertion-biased models are available upon request.)

Composition and Function of Telomerase-A Polymerase Associated with the Origin of Eukaryotes

The canonical DNA polymerases involved in the replication of the genome are unable to fully replicate the physical ends of linear chromosomes, called telomeres. Chromosomal termini thus become shortened in each cell cycle. The maintenance of telomeres requires telomerase—a specific RNA-dependent DNA polymerase enzyme complex that carries its own RNA template and adds telomeric repeats to the ends of chromosomes using a reverse transcription mechanism. Both core subunits of telomerase—its catalytic telomerase reverse transcriptase (TERT) subunit and telomerase RNA (TR) component—were identified in quick succession in Tetrahymena more than 30 years ago. Since then, both telomerase subunits have been described in various organisms including yeasts, mammals, birds, reptiles and fish. Despite the fact that telomerase activity in plants was described 25 years ago and the TERT subunit four years later, a genuine plant TR has only recently been identified by our group. In this review, we focus on the structure, composition and function of telomerases. In addition, we discuss the origin and phylogenetic divergence of this unique RNA-dependent DNA polymerase as a witness of early eukaryotic evolution. Specifically, we discuss the latest information regarding the recently discovered TR component in plants, its conservation and its structural features.

Thaumatin-Like Protein (TLP) Gene Family in Barley: Genome-Wide Exploration and Expression Analysis during Germination

Thaumatin-like Proteins (TLPs) are known to play a vital role in plant defense, developmental processes and seed germination. We identified 19 TLP genes from the reference genome of barley and 37, 28 and 35 TLP genes from rice, Brachypodium and sorghum genomes, respectively. Comparative phylogenetic analysis classified the TLP family into nine groups. Localized gene duplications with diverse exon/intron structures contributed to the expansion of the TLP gene family in cereals. Most of the barley TLPs were localized on chromosome 5H. The spatiotemporal expression pattern of HvTLP genes indicated their predominant expression in the embryo, developing grains, root and shoot tissues. Differential expression of HvTLP14, HvTLP17 and HvTLP18 in the malting variety (Morex) over 16-96 h of grain germination revealed their possible role in malting. This study provides a description of the TLP gene family in barley and their possible involvement in seed germination and the malting process.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.

Division of labour in a matrix, rather than phagocytosis or endosymbiosis, as a route for the origin of eukaryotic cells

Abstract

Two apparently irreconcilable models dominate research into the origin of eukaryotes. In one model, amitochondrial proto-eukaryotes emerged autogenously from the last universal common ancestor of all cells. Proto-eukaryotes subsequently acquired mitochondrial progenitors by the phagocytic capture of bacteria. In the second model, two prokaryotes, probably an archaeon and a bacterial cell, engaged in prokaryotic endosymbiosis, with the species resident within the host becoming the mitochondrial progenitor. Both models have limitations. A search was therefore undertaken for alternative routes towards the origin of eukaryotic cells. The question was addressed by considering classes of potential pathways from prokaryotic to eukaryotic cells based on considerations of cellular topology. Among the solutions identified, one, called here the “third-space model”, has not been widely explored. A version is presented in which an extracellular space (the third-space), serves as a proxy cytoplasm for mixed populations of archaea and bacteria to “merge” as a transitionary complex without obligatory endosymbiosis or phagocytosis and to form a precursor cell. Incipient nuclei and mitochondria diverge by division of labour. The third-space model can accommodate the reorganization of prokaryote-like genomes to a more eukaryote-like genome structure. Nuclei with multiple chromosomes and mitosis emerge as a natural feature of the model. The model is compatible with the loss of archaeal lipid biochemistry while retaining archaeal genes and provides a route for the development of membranous organelles such as the Golgi apparatus and endoplasmic reticulum. Advantages, limitations and variations of the “third-space” models are discussed.

Reviewers

This article was reviewed by Damien Devos, Buzz Baum and Michael Gray.

Comprehensive genomic analyses with 115 plastomes from algae to seed plants: structure, gene contents, GC contents, and introns

BACKGROUND

Chloroplasts are a common character in plants. The chloroplasts in each plant lineage have shaped their own genomes, plastomes, by structural changes and transferring many genes to nuclear genomes during plant evolution. Some plastid genes have introns that are mostly group II introns.

OBJECTIVE

This study aimed to get genomic and evolutionary insights on the plastomes from green algae to flowering plants.

METHODS

Plastomes of 115 species from green algae, bryophytes, pteridophytes (spore bearing vascular plants), gymnosperms, and angiosperms were mined from NCBI organelle genome database. Plastome structure, gene contents and GC contents were analyzed by the in-house developed Phyton code. Intronic features including presence/absence, length, intron phases were analyzed by manually in the annotated information in NCBI.

RESULTS

The canonical quadripartite structures were retained in most plastomes except of a few plastomes that had lost an invert repeat (IR). Expansion or reduction or deletion of IRs resulted in the length variation of the plastomes. The number of protein coding genes ranged from 40 to 92 with an average 79.43 ± 5.84 per plastome and gene losses were apparent in specific lineages. The number of trn genes ranged from 13 to 33 with an average 21.19 ± 2.42 per plastome. Ribosomal RNA genes, rrn, were located in the IRs so that they were present in a duplicate except of the species that had lost one of the IR. GC contents were variable from 24.9 to 51.0% with an average 38.21 ± 3.27%, indicating bias to high AT contents. Plastid introns were present in 18 protein coding genes, six trn genes, and one rrn gene. Intron losses occurred among the orthologous genes in different plant lineages. The plastid introns were long compared with the nuclear introns, which might be related with the spliceosome nuclear introns and self-splicing group II plastid introns. The trnK-UUU intron contained the maturase encoding matK gene except in the chlorophyte algae and monilophyte ferns in which the trnK-UUU was lost, but matK retained. There were many annotation artefacts in the intron positions in the NCBI database. In the analysis of intron phases, phase 0 introns were more frequent than those of phase 2 and 3 introns. Phase polymorphism was observed in the introns of clpP which was derived from nucleotide insertion. Plastid trn introns were long compared to the archaeal or eukaryotic nuclear tRNA introns. Of the six plastid trn introns, one was at the D loop and other five were at the anticodon loop. The insertion sites were conserved among the trn genes in archaea, eukaryotic nuclear and plastid tRNA genes.

CONCLUSIONS

Current study refurbrished the previous findings of structural variations, gene contents, and GC contents of the chloroplast genomes from green algae to flowering plants. The study also included some noble findings and discussions on the plastome introns including their length variations and phase variation. We also presented and corrected some false annotations on the introns in protein coding and tRNA genes in the genome database, which might be confirmed by the chloroplast transcriptome analysis in the future.

Origin, Diversity, and Evolution of Telomere Sequences in Plants

Telomeres are basic structures of eukaryote genomes. They distinguish natural chromosome ends from double-stranded breaks in DNA and protect chromosome ends from degradation or end-to-end fusion with other chromosomes. Telomere sequences are usually tandemly arranged minisatellites, typically following the formula (TxAyGz)n. Although they are well conserved across large groups of organisms, recent findings in plants imply that their diversity has been underestimated. Changes in telomeres are of enormous evolutionary importance as they can affect whole-genome stability. Even a small change in the telomere motif of each repeat unit represents an important interference in the system of sequence-specific telomere binding proteins. Here, we provide an overview of telomere sequences, considering the latest phylogenomic evolutionary framework of plants in the broad sense (Archaeplastida), in which new telomeric sequences have recently been found in diverse and economically important families such as Solanaceae and Amaryllidaceae. In the family Lentibulariaceae and in many groups of green algae, deviations from the typical plant telomeric sequence have also been detected recently. Ancestry and possible homoplasy in telomeric motifs, as well as extant gaps in knowledge are discussed. With the increasing availability of genomic approaches, it is likely that more telomeric diversity will be uncovered in the future. We also discuss basic methods used for telomere identification and we explain the implications of the recent discovery of plant telomerase RNA on further research about the role of telomerase in eukaryogenesis or on the molecular causes and consequences of telomere variability.

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

Telomeres ensure chromosome length homeostasis and protection from catastrophic end-to-end chromosome fusions. All eukaryotes require this essential, strictly conserved telomere-dependent genome preservation. However, recent evolutionary analyses of mammals, plants, and flies report pervasive rapid evolution of telomere proteins. The causes of this paradoxical observation - that unconserved machinery underlies an essential, conserved function - remain enigmatic. Indeed, these fast-evolving telomere proteins bind, extend, and protect telomeric DNA, which itself evolves slowly in most systems. We hypothesize that the universally fast-evolving subtelomere - the telomere-adjacent, repetitive sequence - is a primary driver of the 'telomere paradox'. Under this model, radical sequence changes in the subtelomere perturb subtelomere-dependent, telomere functions. Compromised telomere function then spurs adaptation of telomere proteins to maintain telomere length homeostasis and protection. We propose an experimental framework that leverages both protein divergence and subtelomeric sequence divergence to test the hypothesis that subtelomere sequence evolution shapes recurrent innovation of telomere machinery.

Retroelement origins of pre-mRNA splicing

Recent cryo-EM structures of a group II intron caught in the process of invading DNA have given new insight into the mechanisms of both splicing and retrotransposition. Conformational dynamics involving the branch-site helix domain VI are responsible for substrate exchange between the two steps of splicing. These structural rearrangements have strong parallels with the movement of the branch-site helix in the spliceosome during catalysis. This is strong evidence for the spliceosome evolving from a group II intron ancestor. We observe other topological changes in the overall structure of the catalytic domain V that may occur in the spliceosome as well. Therefore, studying group II introns not only provides us with insight into the evolutionary origins of the spliceosome, but also may inform the design of experiments to further probe structure-function relationships in this eukaryotic splicing apparatus. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Telomerase RNAs in land plants

To elucidate the molecular nature of evolutionary changes of telomeres in the plant order Asparagales, we aimed to characterize telomerase RNA subunits (TRs) in these plants. The unusually long telomere repeat unit in Allium plants (12 nt) allowed us to identify TRs in transcriptomic data of representative species of the Allium genus. Orthologous TRs were then identified in Asparagales plants harbouring telomere DNA composed of TTAGGG (human type) or TTTAGGG (Arabidopsis-type) repeats. Further, we identified TRs across the land plant phylogeny, including common model plants, crop plants, and plants with unusual telomeres. Several lines of functional testing demonstrate the templating telomerase function of the identified TRs and disprove a functionality of the only previously reported plant telomerase RNA in Arabidopsis thaliana. Importantly, our results change the existing paradigm in plant telomere biology which has been based on the existence of a relatively conserved telomerase reverse transcriptase subunit (TERT) associating with highly divergent TRs even between closely related plant taxa. The finding of a monophyletic origin of genuine TRs across land plants opens the possibility to identify TRs directly in transcriptomic or genomic data and/or predict telomere sequences synthesized according to the respective TR template region.

An Integrated Model of Minor Intron Emergence and Conservation

Minor introns constitute <0.5% of the introns in the human genome and have remained an enigma since their discovery. These introns are removed by a distinct splicing complex, the minor spliceosome. Both are ancient, tracing back to the last eukaryotic common ancestor (LECA), which is reflected by minor intron enrichment in specific gene families, such as the mitogen activated-protein kinase kinases, voltage-gated sodium and calcium ion channels, and E2F transcription factors. Most minor introns occur as single introns in genes with predominantly major introns. Due to this organization, minor intron-containing gene (MIG) expression requires the coordinated action of two spliceosomes, which increases the probability of missplicing. Thus, one would expect loss of minor introns via purifying selection. This has resulted in complete minor intron loss in at least nine eukaryotic lineages. However, minor introns are highly conserved in land plants and metazoans, where their importance is underscored by embryonic lethality when the minor spliceosome is inactivated. Conditional inactivation of the minor spliceosome has shown that rapidly dividing progenitor cells are highly sensitive to minor spliceosome loss. Indeed, we found that MIGs were significantly enriched in a screen for genes essential for survival in 341 cycling cell lines. Here, we propose that minor introns inserted randomly into genes in LECA or earlier and were subsequently conserved in genes crucial for cycling cell survival. We hypothesize that the essentiality of MIGs allowed minor introns to endure through the unicellularity of early eukaryotic evolution. Moreover, we identified 59 MIGs that emerged after LECA, and that many of these are essential for cycling cell survival, reinforcing our essentiality model for MIG conservation. This suggests that minor intron emergence is dynamic across eukaryotic evolution, and that minor introns should not be viewed as molecular fossils. We also posit that minor intron splicing was co-opted in multicellular evolution as a regulatory switch for en masse control of MIG expression and the biological processes they regulate. Specifically, this mode of regulation could control cell proliferation and thus body size, an idea supported by domestication syndrome, wherein MIGs are enriched in common candidate animal domestication genes.

Introns: Good Day Junk Is Bad Day Treasure

Introns are ubiquitous in eukaryotic transcripts. They are often viewed as junk RNA but the huge energetic burden of transcribing, removing, and degrading them suggests a significant evolutionary advantage. Ostensibly, an intron functions within the host pre-mRNA to regulate its splicing, transport, and degradation. However, recent studies have revealed an entirely new class of trans-acting functions where the presence of intronic RNA in the cell impacts the expression of other genes in trans. Here, we review possible new mechanisms of intron functions, with a focus on the role of yeast introns in regulating the cell growth response to starvation.

Functional characterization of Gh_A08G1120 (GH3.5) gene reveal their significant role in enhancing drought and salt stress tolerance in cotton

BACKGROUND

Auxins play an important role in plant growth and development; the auxins responsive gene; auxin/indole-3-acetic acid (Aux/IAA), small auxin-up RNAs (SAUR) and Gretchen Hagen3 (GH3) control their mechanisms. The GH3 genes function in homeostasis by the catalytic activities in auxin conjugation and bounding free indole-3-acetic acid (IAA) to amino acids.

RESULTS

In our study, we identified the GH3 genes in three cotton species; Gossypium hirsutum, Gossypium arboreum and Gossypium raimondii, analyzed their chromosomal distribution, phylogenetic relationships, cis-regulatory element function and performed virus induced gene silencing of the novel Gh_A08G1120 (GH3.5) gene. The phylogenetic tree showed four clusters of genes with clade 1, 3 and 4 having mainly members of the GH3 of the cotton species while clade 2 was mainly members belonging to Arabidopsis. There were no paralogous genes, and few orthologous genes were observed between Gossypium and other species. All the GO terms were detected, but only 14 genes were found to have described GO terms in upland cotton, more biological functions were detected, as compared to the other functions. The GH3.17 subfamily harbored the highest number of the cis-regulatory elements, most having promoters towards dehydration-responsiveness. The RNA expression analysis revealed that 10 and 8 genes in drought and salinity stress conditions respectively were upregulated in G. hirsutum. All the genes that were upregulated in plants under salt stress conditions were also upregulated in drought stress; moreover, Gh_A08G1120 (GH3.5) exhibited a significant upregulation across the two stress factors. Functional characterization of Gh_A08G1120 (GH3.5) through virus-induced gene silencing (VIGS) revealed that the VIGS plants ability to tolerate drought and salt stresses was significantly reduced compared to the wild types. The chlorophyll content, relative leaf water content (RLWC), and superoxide dismutase (SOD) concentration level were reduced significantly while malondialdehyde concentration and ion leakage as a measure of cell membrane stability (CMS) increased in VIGS plants under drought and salt stress conditions.

CONCLUSION

This study revealed the significance of the GH3 genes in enabling the plant's adaptation to drought and salt stress conditions as evidenced by the VIGS results and RT-qPCR analysis.

Analysis of the Functional Relevance of Epigenetic Chromatin Marks in the First Intron Associated with Specific Gene Expression Patterns

We previously showed that the first intron of genes exhibits several interesting characteristics not seen in other introns: 1) it is the longest intron on average in almost all eukaryotes, 2) it presents the highest number of conserved sites, and 3) it exhibits the highest density of regulatory chromatin marks. Here, we expand on our previous study by integrating various multiomics data, leading to further evidence supporting the functionality of sites in the first intron. We first show that trait-associated single-nucleotide polymorphisms (TASs) are significantly enriched in the first intron. We also show that within the first intron, the density of epigenetic chromatin signals is higher near TASs than in distant regions. Furthermore, the distribution of several chromatin regulatory marks is investigated in relation to gene expression specificity (i.e., housekeeping vs. tissue-specific expression), essentiality (essential genes vs. nonessential genes), and levels of gene expression; housekeeping genes or essential genes contain greater proportions of active chromatin marks than tissue-specific genes or nonessential genes, and highly expressed genes exhibit a greater density of chromatin regulatory marks than genes with low expression. Moreover, we observe that genes carrying multiple first-intron TASs interact with each other within a large protein-protein interaction network, ultimately connecting to the UBC protein, a well-established protein involved in ubiquitination. We believe that our results shed light on the functionality of first introns as a genomic entity involved in gene expression regulation.

Enrichment of rare alleles within epigenetic chromatin marks in the first intron

In previous studies, we demonstrated that some sites in the first intron likely regulate gene expression. In the present work, we sought to further confirm the functional relevance of first intron sites by estimating the quantity of rare alleles in the first intron. A basic hypothesis posited herein is that genomic regions carrying more functionally important sites will have a higher proportion of rare alleles. We estimated the proportions of rare single nucleotide polymorphisms with a minor allele frequency < 0.01 located in several histone marks in the first introns of various genes, and compared them with those in other introns and those in 2-kb upstream regions. As expected, rare alleles were found to be significantly enriched in most of the regulatory sites located in the first introns. Meanwhile, transcription factor binding sites were significantly more enriched in the 2-kb upstream regions (i.e., the regions of putative promoters of genes) than in the first introns. These results strongly support our proposal that the first intron sites of genes may have important regulatory functions in gene expression independent of promoters.

Excised linear introns regulate growth in yeast

Spliceosomal introns are ubiquitous non-coding RNAs that are typically destined for rapid debranching and degradation. Here we describe 34 excised introns in Saccharomyces cerevisiae that-despite being rapidly degraded in log-phase growth-accumulate as linear RNAs under either saturated-growth conditions or other stresses that cause prolonged inhibition of TORC1, which is a key integrator of growth signalling. Introns that become stabilized remain associated with components of the spliceosome and differ from other spliceosomal introns in having a short distance between their lariat branch point and 3' splice site, which is necessary and sufficient for their stabilization. Deletion of these unusual introns is disadvantageous in saturated conditions and causes aberrantly high growth rates in yeast that are chronically challenged with the TORC1 inhibitor rapamycin. The reintroduction of native or engineered stable introns suppresses this aberrant rapamycin response. Thus, excised introns function within the TOR growth-signalling network of S. cerevisiae and, more generally, excised spliceosomal introns can have biological functions.