Origin of new genes: evidence from experimental and computational analyses

Exon shuffling is an essential molecular mechanism for the formation of new genes. Many cases of exon shuffling have been reported in vertebrate genes. These discoveries revealed the importance of exon shuffling in the origin of new genes. However, only a few cases of exon shuffling were reported from plants and invertebrates, which gave rise to the assertion that the intron-mediated recombination mechanism originated very recently. We focused on the origin of new genes by exon shuffling and retroposition. We will first summarize our experimental work, which revealed four new genes in Drosophila, plants, and humans. These genes are 10⁶ to 10⁸ million years old. The recency of these genes allows us to directly examine the origin and evolution of genes in detail. These observations show firstly the importance of exon shuffling and retroposition in the rapid creation of new gene structures. They also show that the resultant chimerical structures appearing as mosaic proteins or as retroposed coding structures with novel regulatory systems, often confer novel functions. Furthermore, these newly created genes appear to have been governed by positive Darwinian selection throughout their history, with rapid changes of amino acid sequence and gene structure in very short periods of evolution. We further analyzed the distribution of intron phases in three non-vertebrate species, Drosophila melanogaster, Caenorhabditis elegans, and Arabidosis thaliana, as inferred from their genome sequences. As in the case of vertebrate genes, we found that intron phases in these species are unevenly distributed with an excess of phase zero introns and a significant excess of symmetric exons. Both findings are consistent with the requirements for the molecular process of exon shuffling. Thus, these non-vertebrate genomes may have also been strongly impacted by exon shuffling in general.     

Evolution of Prokaryotic Genes by Shift of Stop Codons

De novo origin of coding sequence remains an obscure issue in molecular evolution. One of the possible paths for addition (subtraction) of DNA segments to (from) a gene is stop codon shift. Single nucleotide substitutions can destroy the existing stop codon, leading to uninterrupted translation up to the next stop codon in the gene’s reading frame, or create a premature stop codon via a nonsense mutation. Furthermore, short indels-caused frameshifts near gene’s end may lead to premature stop codons or to translation past the existing stop codon. Here, we describe the evolution of the length of coding sequence of prokaryotic genes by change of positions of stop codons. We observed cases of addition of regions of 3′UTR to genes due to mutations at the existing stop codon, and cases of subtraction of C-terminal coding segments due to nonsense mutations upstream of the stop codon. Many of the observed stop codon shifts cannot be attributed to sequencing errors or rare deleterious variants segregating within bacterial populations. The additions of regions of 3′UTR tend to occur in those genes in which they are facilitated by nearby downstream in-frame triplets which may serve as new stop codons. Conversely, subtractions of coding sequence often give rise to in-frame stop codons located nearby. The amino acid composition of the added region is significantly biased, compared to the overall amino acid composition of the genes. Our results show that in prokaryotes, shift of stop codon is an underappreciated contributor to functional evolution of gene length.     

Repetitive element-mediated recombination as a mechanism for new gene origination in Drosophila

Previous studies of repetitive elements (REs) have implicated a mechanistic role in generating new chimerical genes. Such examples are consistent with the classic model for exon shuffling, which relies on non-homologous recombination. However, recent data for chromosomal aberrations in model organisms suggest that ectopic homology-dependent recombination may also be important. Lack of a dataset comprising experimentally verified young duplicates has hampered an effective examination of these models as well as an investigation of sequence features that mediate the rearrangements. Here we use approximately 7,000 cDNA probes (approximately 112,000 primary images) to screen eight species within the Drosophila melanogaster subgroup and identify 17 duplicates that were generated through ectopic recombination within the last 12 mys. Most of these are functional and have evolved divergent expression patterns and novel chimeric structures. Examination of their flanking sequences revealed an excess of repetitive sequences, with the majority belonging to the transposable element DNAREP1 family, associated with the new genes. Our dataset strongly suggests an important role for REs in the generation of chimeric genes within these species.     

Evolution of the intron-exon structure of eukaryotic genes

The origin and evolution of intron-exon structures continue to be controversial topics. Two alternative theories, the ‘exon theory of genes’ and the ‘insertional theory of introns’, debate the presence or absence of introns in primordial genes. Both sides of the argument have focused on the positions of introns with respect to protein and gene structures. A new approach has emerged in the study of the evolution of intron-exon structures: a population analysis of genes. One example is the statistical analysis of intron phases — the position of introns within or between codons. This analysis detected a significant signal of exon shuffling in the DNA sequence database containing both ancient and modern exon sequences: intron phase correlations, that is, the association together within genes of introns of the same phase. The results of this analysis suggest that exon shuffling played an important role in the origin of both ancient and modern genes.     

Lineage-specific diversification of killer cell Ig-like receptors in the owl monkey,a New World primate

Killer cell Ig-like receptors (KIRs) modulate the cytotoxic effects of natural killer cells. In primates, the KIRs are highly diverse as a consequence of variation in gene content, alternative domain composition, and loci polymorphism. We analyzed a bacterial artificial chromosome (BAC) clone draft sequence spanning the owl monkey KIR cluster. The draft sequence had seven ordered yet unconnected contigs containing six full-length and two partial gene models, flanked by the LILRB and FcAR framework genes. Gene models were predicted to encode KIRs with inhibitory, activating, or dual functionality. Four gene models encoded three Ig domain receptors, while three others encoded molecules with four Ig domains. The additional domain resulted from an insertion in tandem of a 2,101 bp fragment containing the last 289 bp of intron 2, exon 3, and intron 3, resulting in molecules with two D0 domains. Re-screening of the owl monkey BAC library and sequencing of partial cDNAs from an owl monkey yielded five additional KIRs, four of which encoded receptors with short cytoplasmic domains with premature stop codons due to either a single nucleotide substitution or deletion or the absence of exon 8. Phylogenetic analysis by domains showed that owl monkey KIRs were monophyletic, clustering independently from other primate KIR lineages. Retroelements found in introns, however, were shared by KIRs from different primate lineages. This suggests that the owl monkey inherited a KIR cluster with a rich history of exon shuffling upon which positive selection for ligand binding operated to diversify the receptors in a lineage-specific fashion. Nucleotide sequence data reported are available in the GenBank database under accession numbers FJ154791–5.     

Nonsense-mediated mRNA decay in Drosophila: at the intersection of the yeast and mammalian pathways

The nonsense-mediated mRNA decay (NMD) pathway promotes the rapid degradation of mRNAs containing premature stop codons (PTCs). In Caenorhabditis elegans, seven genes (smg1-7) playing an essential role in NMD have been identified. Only SMG2-4 (known as UPF1-3) have orthologs in Saccharomyces cerevisiae. Here we show that the Drosophila orthologs of UPF1-3, SMG1, SMG5 and SMG6 are required for the degradation of PTC-containing mRNAs, but that there is no SMG7 ortholog in this organism. In contrast, orthologs of SMG5-7 are encoded by the human genome and all three are required for NMD. In human cells, exon boundaries have been shown to play a critical role in defining PTCs. This role is mediated by components of the exon junction complex (EJC). Contrary to expectation, however, we show that the components of the EJC are dispensable for NMD in Drosophila cells. Consistently, PTC definition occurs independently of exon boundaries in DROSOPHILA: Our findings reveal that despite conservation of the NMD machinery, different mechanisms have evolved to discriminate premature from natural stop codons in metazoa.     

Formation and Longevity of Chimeric and Duplicate Genes in Drosophila melanogaster

Historically, duplicate genes have been regarded as a major source of novel genetic material. However, recent work suggests that chimeric genes formed through the fusion of pieces of different genes may also contribute to the evolution of novel functions. To compare the contribution of chimeric and duplicate genes to genome evolution, we measured their prevalence and persistence within Drosophila melanogaster. We find that ~80.4 duplicates form per million years, but most are rapidly eliminated from the genome, leaving only 4.1% to be preserved by natural selection. Chimeras form at a comparatively modest rate of ~11.4 per million years but follow a similar pattern of decay, with ultimately only 1.4% of chimeras preserved. We propose two mechanisms of chimeric gene formation, which rely entirely on local, DNA-based mutations to explain the structure and placement of the youngest chimeric genes observed. One involves imprecise excision of an unpaired duplication during large-loop mismatch repair, while the other invokes a process akin to replication slippage to form a chimeric gene in a single event. Our results paint a dynamic picture of both chimeras and duplicate genes within the genome and suggest that chimeric genes contribute substantially to genomic novelty.     

四膜虫eRF1识别终止密码子的功能结构域

原生动物的一些纤毛虫中终止密码子发生重分配现象,将1个或2个终止密码子翻译为氨基酸.目前对这一现象的发生机制仍无合理解释．近年来,对蛋白质合成终止过程中肽链释放因子(eukaryotic polypeptide release factor, eRF)结构和功能的深入研究,为揭示终止密码子的重分配机制提供了重要的线索．本实验以具有终止密码子识别特异性的四膜虫Tt-eRF1为研究对象,将其中与密码子识别有关的GTx、NIKS和Y-C-F关键模体(motif) 引入识别通用终止密码子的酵母Sc=eRF1中,构建成各种嵌合体eRF1．利用双荧光素酶报告系统和细胞活性实验,分析关键模体及其周边的氨基酸对eRF1识别终止密码子性质的影响．结果表明,GTx和NIKS模体一定程度上决定eRF1识别终止密码子第1位碱基U和第2位碱基A;Y-C-F模体决定eRF1识别终止密码子UGA的第2位碱基G．模体内及其相邻氨基酸定点突变分析进一步支持以上结果．本研究推测,eRF1在进化过程中一些关键模体结构的改变决定其识别终止密码子的特异性,只能识别3个终止密码子中的1个或2个．随后,由于tRNA基因的突变产生阻抑性tRNA,促成终止密码子在原生动物纤毛虫中的重新分配．     

Positive selection in phytotoxic protein-encoding genes of Botrytis species

Evolutionary patterns of sequence divergence were analyzed in genes from the fungal genus Botrytis (Ascomycota), encoding phytotoxic proteins homologous to a necrosis and ethylene-inducing protein from Fusarium oxysporum. Fragments of two paralogous genes (designated NEP1 and NEP2) were amplified from all known Botrytis species and sequenced. NEP1 sequences of two Botrytis species contain premature stop codons, indicating that they may be non-functional. Both paralogs of all species encode proteins with a remarkably similar predicted secondary structure, however, they contain different types of post-translational modification motifs, which are conserved across the genus. While both NEP genes are, overall, under purifying selection, we identified a number of amino acids under positive selection based on inference using maximum likelihood models. Positively selected amino acids in NEP1 were not under selection in corresponding positions in NEP2. The biological significance of positively selected residues and the role of NEP proteins in pathogenesis remain to be resolved.     

Is Esterase-P Encoded by a Cryptic Pseudogene in Drosophila Melanogaster?

We have amplified and sequenced the gene encoding Esterase-P (Est-P) in 10 strains of Drosophila melanogaster. Three premature termination codons occur in the coding region of the gene in two strains. This observation, together with other indirect evidence, leads us to propose that Est-P may be a pseudogene in D. melanogaster. Est-P would be a ``cryptic' pseudogene, in the sense that it retains intact the coding sequence (without stop codons and other alterations usually observed in pseudogenes) in most D. melanogaster strains. We conjecture that the β-esterase cluster may consist in other Drosophila species of functional and nonfunctional genes. We also conjecture that the rarity of detected pseudogenes in Drosophila may be due to the difficulty of discovering them, because most of them are cryptic.     

Molecular Characterization of Two Silenced y-type Genes for Glu-B1 in Triticum aestivum ssp. yunnanese and ssp. tibetanum

The high molecular weight glutenln subunits (HMW-GSs) are a major class of common wheat storage proteins. The bread-making quality of common wheat flour is influenced by the composition of HMW-GSs. In the present study, two unexpressed 1By genes from Triticum aesitvum L.ssp.yunnanese AS332 and T. aesitvum sep.tibetanum AS908 were respectively cloned and characterized. The results indicated that both of the silenced 1By genes in AS332 and AS908 were 1By9. in contrast to previously reported mechanisms for silenced genes 1Ax and 1Ay, which was due to the insertion of transposon elements or the presence of premature stop codon via base substitution of C→T transition in tdnucleotides CAA or CAG, the silence of 1By9 genes was caused by premature stop codons via the deletion of base A in tdnucleotide CAA, which lead to frameshift mutation and indirectly produced several premature stop codons (TAG) downstream of the coding sequence.     

Give-and-take: interactions between DNA transposons and their host plant genomes

Recent genome sequencing efforts have revealed how extensively transposable elements (TEs) have contributed to the shaping of present day plant genomes. DNA transposons associate preferentially with the euchromatic or genic component of plant genomes and have had the opportunity to interact intimately with the genes of the plant host. These interactions have resulted in TEs acquiring host sequences, forming chimeric genes through exon shuffling, replacing regulatory sequences, mobilizing genes around the genome, and contributing genes to the host. The close interaction of transposons with genes has also led to the evolution of intricate cellular mechanisms for silencing transposon activity. Transposons have thus become important subjects of study in understanding epigenetic regulation and, in cases where transposons have amplified to high numbers, how to escape that regulation.     

Frequent appearance of novel protein-coding sequences by frameshift translation

Genomic duplication, followed by divergence, contributes to organismal evolution. Several mechanisms, such as exon shuffling and alternative splicing, are responsible for novel gene functions, but they generate homologous domains and do not usually lead to drastic innovation. Major novelties can potentially be introduced by frameshift mutations and this idea can explain the creation of novel proteins. Here, we employ a strategy using simulated protein sequences and identify 470 human and 108 mouse frameshift events that originate new gene segments. No obvious interspecies overlap was observed, suggesting high rates of acquisition of evolutionary events. This inference is supported by a deficiency of TpA dinucleotides in the protein-coding sequences, which decreases the occurrence of translational termination, even on the complementary strand. Increased usage of the TGA codon as the termination signal in newer genes also supports our inference. This suggests that tolerated frameshift changes are a prevalent mechanism for the rapid emergence of new genes and that protein-coding sequences can be derived from existing or ancestral exons rather than from events that result in noncoding sequences becoming exons.     

A 17 bp deletion in the Bone Morphogenetic Protein 15 (BMP15) gene is associated to increased prolificacy in the Rasa Aragonesa sheep breed

Different mutations in the Bone Morphogenetic Protein 15 (BMP15) and the Growth Differentiation Factor 9 (GDF9) genes cause increased ovulation rate and infertility in a dosage-sensitive manner in sheep. They cause increased ovulation rate and twin and triplet births in heterozygotes, and complete primary ovarian failure in homozygotes resulting in total infertility. We are here presenting a novel mutation in the second exon of the ovine BMP15 gene, found in the Spanish breed Rasa Aragonesa. It consists of a 17 bp deletion resulting in displacement of the open reading frame and premature stop codons. As a consequence, nearly 85% of the sequence of the wild type aminoacidic chain in the second exon of the BMP15 pro-protein is modified or suppressed as only the first 45 amino acids are conserved of the 245 original. The mature peptide is lost. The ewes heterozygous for this deletion present very high prolificacy (2.66 lambs/birth) when compared to a mean flock prolificacy of 1.36 lambs. The deletion causes a complete lack of functionality of the second exon of BMP15, comparable to the effect of premature stop codons in other mutations. Therefore, homozygous females for the deletion are expected to present primary ovarian failure. DNA sequence analysis of the GDF9 coding regions detected only a synonymous Single Nucleotide Polymorphism (SNP), apparently not linked to changes in prolificacy.     

Exon size distribution and the origin of introns

Since it was first recognised that eukaryotic genes are fragmented into coding segments (exons) separated by non-coding segments (introns), the reason for this phenomenon has been debated. There are two dominant theories: that the piecewise arrangement of genes allows functional protein domains, represented by exons, to recombine by shuffling to form novel proteins with combinations of functions; or that introns represent parasitic DNA that can infest the eukaryotic genome because it does not interfere grossly with the fitness of its host. Differing distributions of exon lengths are predicted by these two theories. In this paper we examine distributions of exon lengths for six different organisms and find that they offer empirical evidence that both theories may in part be correct.