Divergence of Gene Body DNA Methylation and Evolution of Plant Duplicate Genes

It has been shown that gene body DNA methylation is associated with gene expression. However, whether and how deviation of gene body DNA methylation between duplicate genes can influence their divergence remains largely unexplored. Here, we aim to elucidate the potential role of gene body DNA methylation in the fate of duplicate genes. We identified paralogous gene pairs from Arabidopsis and rice (Oryza sativa ssp. japonica) genomes and reprocessed their single-base resolution methylome data. We show that methylation in paralogous genes nonlinearly correlates with several gene properties including exon number/gene length, expression level and mutation rate. Further, we demonstrated that divergence of methylation level and pattern in paralogs indeed positively correlate with their sequence and expression divergences. This result held even after controlling for other confounding factors known to influence the divergence of paralogs. We observed that methylation level divergence might be more relevant to the expression divergence of paralogs than methylation pattern divergence. Finally, we explored the mechanisms that might give rise to the divergence of gene body methylation in paralogs. We found that exonic methylation divergence more closely correlates with expression divergence than intronic methylation divergence. We show that genomic environments (e.g., flanked by transposable elements and repetitive sequences) of paralogs generated by various duplication mechanisms are associated with the methylation divergence of paralogs. Overall, our results suggest that the changes in gene body DNA methylation could provide another avenue for duplicate genes to develop differential expression patterns and undergo different evolutionary fates in plant genomes.     

Low rates of silent substitution in nuclear genes of two distantly related Scrophulariaceae (Antirrhinum and Verbascum)

Low levels of genetic diversity and divergence at nuclear loci have previously been observed for cycloidea and fil1-like genes within and between several Antirrhinum species, and divergence at these loci is also low between species in genera at different levels of relatedness in the former family Scrophulariaceae (Digitalis and Verbascum). The low divergence values are surprising, because (based on the sequences of chloroplast loci) the Scrophulariaceae are thought to be polyphyletic, with two anciently diverged clades, and the species we compared belonged to the two different clades. Here, we extend our studies of sequence divergence to more nuclear genes: fil2, far, globosa, and ADH: Detailed studies revealed that in Antirrhinum these genes belong to gene families. Low levels of divergence between Antirrhinum and Verbascum were observed for four of the loci studied, fil2-1, fil2-2, far-L, and globosa, similar to our previous observations. We discuss hypotheses to explain these low synonymous divergence values. For Adh, no cases of very similar sequences were found, but, rather, our sequences from the three different genera (Antirrhinum, Digitalis, and Verbascum) were all very diverged. Repeated gene duplication and loss of elements in the Adh gene family is likely in these lineages, making it impossible to determine orthology of the Adh genes.     

Likelihood Analysis of the Chalcone Synthase Genes Suggests the Role of Positive Selection in Morning Glories (<Emphasis Type="Italic">Ipomoea</Emphasis>)

Chalcone synthase (CHS) is a key enzyme in the biosynthesis of flavonoides, which are important for the pigmentation of flowers and act as attractants to pollinators. Genes encoding CHS constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. In morning glories (Ipomoea), five functional CHS genes (A–E) have been described. Phylogenetic analysis of the Ipomoea CHS gene family revealed that CHS A, B, and C experienced accelerated rates of amino acid substitution relative to CHS D and E. To examine whether the CHS genes of the morning glories underwent adaptive evolution, maximum-likelihood models of codon substitution were used to analyze the functional sequences in the Ipomoea CHS gene family. These models used the nonsynonymous/synonymous rate ratio ( = d_N/d_S) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio test suggested significant variation in selection pressure among amino acid sites, with a small proportion of them detected to be under positive selection along the branches ancestral to CHS A, B, and C. Positive Darwinian selection appears to have promoted the divergence of subfamily ABC and subfamily DE and is at least partially responsible for a rate increase following gene duplication.     

Heterogeneity of the glutathione transferase genes encoding enzymes responsible for insecticide degradation in the housefly

One of the four glutathione-S-transferases (GST) that is overproduced in the insecticide-resistant Cornell-R strain of the housefly (Musca domestica) produces an activity that degrades the insecticide dimethyl parathion and conjugates glutathione to lindane. In earlier work, it was shown that the resistant Cornell-R carries an amplification, probably a duplication, of one or more of its GST loci and that this amplification is directly related to resistance. Using polymerase chain reaction (PCR) amplification with genomic DNA, multiple copies of the gene encoding the parathion-degrading activity (called MdGst-3) were subcloned from both the ancestral, insecticide-susceptible strain BPM and from the insecticide-resistant Cornell-R. In BPM, three different MdGst-3 genes were identified while in Cornell-R, 12 different MdGst-3 sequences were found that, though closely related to ancestral genes, had diverged by a few nucleotides. This diversity in MdGst-3 genomic sequences in Cornell-R is reflected in the expressed sequences, as sampled through a cDNA bank. Population heterozygosity cannot account for these multiple GST genes. We suggest that selection for resistance to insecticides has resulted in not only amplification of the MdGst-3 genes but also in the divergence of sequence between the amplified copies. Received: 22 November 1995 / Accepted: 23 February 1996     

Not born equal: increased rate asymmetry in relocated and retrotransposed rodent gene duplicates

Duplicated genes frequently evolve at different rates. This asymmetry is evidence of natural selection's ability to discriminate between the 2 copies, subjecting them to different levels of purifying selection or even permitting adaptive evolution of one or both copies. However, if gene duplication creates pairs of protein-coding sequences that are initially identical, this raises the question of how selection tells the 2 copies apart. Here, we investigated asymmetric sequence divergence of recently duplicated genes in rodents and related this to 2 possible sources of such asymmetry: gene relocation as a consequence of duplication and retrotransposition as a mechanism of gene duplication. We found that most young rodent duplicates that have been relocated were created by retrotransposition. The degree of rate asymmetry in gene pairs where one copy has been relocated (either by retrotransposition or DNA-based duplication) is greater than in pairs formed by local DNA-based duplication events. Furthermore, by considering the direction of transposition for distant duplicates, we found a consistent tendency for retrogenes to undergo accelerated protein evolution relative to their static paralogs, whereas DNA-based transpositions showed no such tendency. Finally, we demonstrate that the faster sequence evolution of retrogenes correlates with the profound alteration of their expression pattern that is precipitated by retrotransposition.     

Characterization of duplicated Dunaliella viridis SPT1 genes provides insights into early gene divergence after duplication

The sodium-dependent phosphate transporter gene from unicellular green algae Dunaliella viridis, DvSPT1, shares similarity with members of Pi transporter family. Sequencing analysis of D. viridis BAC clone containing the DvSPT1 gene revealed two inverted duplicated copies of this gene (DvSPT1 and DvSPT1-2 respectively). The duplication covered most of both genes except for their 3' downstream region. The duplicated genomic sequences exhibited 97.9% identity with a synonymous divergence of Ks=0.0126 in the coding region. This data indicated very recent gene duplication in D. viridis genome, providing an excellent opportunity to investigate sequence and expression divergence of duplicated genes at an early stage. Scatted point mutations and length polymorphism of simple sequence repeats (SSRs) were predominant among the sequence divergence soon after gene duplication. Due to sequence divergence in the 5' regulatory regions and a swap of the entire 3' downstream regions (3'-UTR), DvSPT1 and DvSPT1-2 showed expression divergence in response to extra-cellular NaCl concentration changes. According to their expression patterns, the two diverged gene copies would provide better adaptation to a broader range of extra-cellular NaCl concentration. Furthermore, Southern blot analysis indicated that there might be a large phosphate transporter gene family in D. viridis.     

Genome structure analysis of molluscs revealed whole genome duplication and lineage specific repeat variation

Comparative genome structure analysis allows us to identify novel genes, repetitive sequences and gene duplications. To explore lineage-specific genomic changes of the molluscs that is good model for development of nervous system in invertebrate, we conducted comparative genome structure analyses of three molluscs, pygmy squid, nautilus and scallops using partial genome shotgun sequencing. Most effective elements on the genome structural changes are repetitive elements (REs) causing expansion of genome size and whole genome duplication producing large amount of novel functional genes. Therefore, we investigated variation and proportion of REs and whole genome duplication. We, first, identified variations of REs in the three molluscan genomes by homology-based and de novo RE detection. Proportion of REs were 9.2%, 4.0%, and 3.8% in the pygmy squid, nautilus and scallop, respectively. We, then, estimated genome size of the species as 2.1, 4.2 and 1.8 Gb, respectively, with 2× coverage frequency and DNA sequencing theory. We also performed a gene duplication assay based on coding genes, and found that large-scale duplication events occurred after divergence from the limpet Lottia, an out-group of the three molluscan species. Comparison of all the results suggested that RE expansion did not relate to the increase in genome size of nautilus. Despite close relationships to nautilus, the squid has the largest portion of REs and smaller genome size than nautilus. We also identified lineage-specific RE and gene-family expansions, possibly relate to acquisition of the most complicated eye and brain systems in the three species.     

Duplication and adaptive evolution of the COR15 genes within the highly cold-tolerant Draba lineage (Brassicaceae)

Plants have evolved diverse adaptive mechanisms that enable them to tolerate abiotic stresses, to varying degrees, and such stresses may have strongly influenced evolutionary changes at levels ranging from molecular to morphological. Previous studies on these phenomena have focused on the adaptive evolution of stress-related orthologous genes in specific lineages. However, heterogenetic evolution of the paralogous genes following duplication has only been examined in a very limited number of stress-response gene families. The COR15 gene encodes a low molecular weight protein that plays an important role in protecting plants from cold stresses. Although two different copies of this gene have been found in the model species, Arabidopsis thaliana, evolutionary patterns of this small gene family in plants have not been previously explored. In this study, we cloned COR15-like sequences and performed evolutionary analyses of these sequences (including those previously reported) in the highly cold-tolerant Draba lineage and related lineages of Brassicaceae. Our phylogenetic analyses indicate that all COR15-like sequences clustered into four clades that corresponded well to the morphological lineages. Gene conversions were found to have probably occurred before/during the divergence of Brassica and Draba lineage. However, repeated, independent duplications of this gene have occurred in different lineages of Brassicaceae. Further comparisons of all sequences suggest that there have been significant inter-lineage differences in evolutionary rates between the duplicated and original genes. We assessed the likelihood that the differences between two well-supported gene subfamilies that appear to have originated from a single duplication, COR15a and COR15b, within the Draba lineage have been driven by adaptive evolution. Comparisons of their non-synonymous/synonymous substitution ratios and rates of predicted amino acid changes indicate that these two gene groups are evolving under different selective pressures and may be functionally divergent. This functional divergence was confirmed by comparing site-specific shifts in evolution indexes of the two groups of predicted proteins. The evidence of differential selection and possible functional divergence suggests that the duplication may be of adaptive significance, with possible implications for the explosive diversification of the Draba lineage during the cooling Quaternary stages and the following worldwide colonization of arid alpine and artic regions.     

Comparative genomic analysis of dinucleotide repeats in Tritryps

The protozoans Trypanosoma cruzi, Trypanosoma brucei and Leishmania major (Tritryps), are evolutionarily ancient eukaryotes which cause worldwide human parasitosis. They present unique biological features. Indeed, canonical DNA/RNA cis-acting elements remain mostly elusive. Repetitive sequences, originally considered as selfish DNA, have been lately recognized as potentially important functional sequence elements in cell biology. In particular, the dinucleotide patterns have been related to genome compartmentalization, gene evolution and gene expression regulation. Thus, we perform a comparative analysis of the occurrence, length and location of dinucleotide repeats (DRs) in the Tritryp genomes and their putative associations with known biological processes. We observe that most types of DRs are more abundant than would be expected by chance. Complementary DRs usually display asymmetrical strand distribution, favoring TT and GT repeats in the coding strands. In addition, we find that GT repeats are among the longest DRs in the three genomes. We also show that specific DRs are non-uniformly distributed along the polycistronic unit, decreasing toward its boundaries. Distinctive non-uniform density patterns were also found in the intergenic regions, with predominance at the vicinity of the ORFs. These findings further support that DRs may control genome structure and gene expression.     

On the nature of gene innovation: duplication patterns in microbial genomes

Gene duplication is considered a major force in gene family expansion and gene innovation. As gene copies assume novel functions, they must avoid periods of neutrality or be deleted from the genome. Current opinions state that copies avoid neutrality through gene dosage effects. These copies are therefore selected from an early stage. This study concentrates on the flow of copies from recent duplication to gene innovation. We have studied 21 microbial genomes using amino acid divergence to describe paralog evolution in the long-term perspective. Five of these were studied in closer detail using nucleotide divergence for a shorter perspective. It was found that rates of duplication and deletion are high, with only a small fraction of duplications retained and apparently selected. This leads to a steady accumulation of paralogs, which seems to be of a similar magnitude in most of the genomes. Furthermore, it is found that genes of high expression level, as measured by their codon bias, are strongly underrepresented among the most recent duplications. Based on these and other observations, it is suggested that gene innovation is driven by amplification of weak, ancillary functions rather than strong, established functions.     

Specificity of Repeat-Induced Point Mutation (Rip) in Neurospora: Sensitivity of Non-Neurospora Sequences, a Natural Diverged Tandem Duplication, and Unique DNA Adjacent to a Duplicated Region

The process designated RIP (repeat-induced point mutation) alters duplicated DNA sequences in the sexual cycle of Neurospora crassa. We tested whether non-Neurospora sequences are susceptible to RIP, explored the basis for the observed immunity to this process of a diverged tandem duplication that probably arose by a natural duplication followed by RIP (the Neurospora zeta-eta region), and investigated whether RIP extends at all into unique sequences bordering a duplicated region. Bacterial sequences of the plasmid pUC8 and of a gene conferring resistance to hygromycin B were sensitive to RIP in N. crassa when repeated in the genome. When the entire 1.6-kb zeta-eta region was duplicated, it was susceptible to RIP, but was affected by it to a lesser extent than other duplications. Only three of 62 progeny from crosses harboring unlinked duplications of the region showed evidence of changes. We attribute the low level of alterations to depletion of mutable sites. The stability of the zeta-eta region in strains having single copies of the region suggests that the 14% divergence of the tandem elements is sufficient to prevent RIP. DNA sequence analysis of unduplicated pUC8 sequences adjacent to a duplication revealed that RIP continued at least 180 bp beyond the boundary of the duplication. Three mutations occurred in the 200-bp segment of bordering sequences examined.     

Molecular trajectories leading to the alternative fates of duplicate genes

Gene duplication generates extra gene copies in which mutations can accumulate without risking the function of pre-existing genes. Such mutations modify duplicates and contribute to evolutionary novelties. However, the vast majority of duplicates appear to be short-lived and experience duplicate silencing within a few million years. Little is known about the molecular mechanisms leading to these alternative fates. Here we delineate differing molecular trajectories of a relatively recent duplication event between humans and chimpanzees by investigating molecular properties of a single duplicate: DNA sequences, gene expression and promoter activities. The inverted duplication of the Glutathione S-transferase Theta 2 (GSTT2) gene had occurred at least 7 million years ago in the common ancestor of African great apes and is preserved in chimpanzees (Pan troglodytes), whereas a deletion polymorphism is prevalent in humans. The alternative fates are associated with expression divergence between these species, and reduced expression in humans is regulated by silencing mutations that have been propagated between duplicates by gene conversion. In contrast, selective constraint preserved duplicate divergence in chimpanzees. The difference in evolutionary processes left a unique DNA footprint in which dying duplicates are significantly more similar to each other (99.4%) than preserved ones. Such molecular trajectories could provide insights for the mechanisms underlying duplicate life and death in extant genomes.     

Molecular evolution of PISTILLATA-like genes in the dogwood genus Cornus (Cornaceae)

The MADS-box gene family encodes critical regulators determining floral organ development. Understanding evolutionary patterns and processes of MADS-box genes is an important step toward unraveling the molecular basis of floral morphological evolution. In this study, we investigated the evolution of PI-like genes of the MADS-box family in the dogwood genus Cornus (Cornaceae). Cornus is a eudicot lineage in the asterids clade, and is intriguing in evolving petaloid bract morphology in two major lineages within the genus. The gene genealogy reconstructed using genomic DNA and cDNA sequences suggests multiple PI-like gene duplication events in Cornus. An ancient duplication event resulted in two ancient paralogs, CorPI-A and CorPI-B, which have highly diverged intron regions. Duplication of CorPI-A further resulted in two paralogs in one subgroup of Cornus, the BW group that does not produce modified bracts. Most species analyzed were found to contain more than one copy of the PI-like gene with most copies derived recently within species. Estimation and comparison of dN/dS ratios revealed relaxed selection in the PI-like gene in Cornus in comparison with the gene in the closely related outgroups Alangium and Davidia, and in other flowering plants. Selection also differed among major gene copies, CorPI-A and CorPI-B, and among different morphological subgroups of Cornus. Variation in selection pressures may indicate functional changes in PI-like genes after gene duplication and among different lineages. Strong positive selection at three amino acid sites of CorPI was also detected from a region critical for dimerization activity. Total substitution rates of the CorPI gene also differ among lineages of Cornus, showing a trend similar to that found in dN/dS ratios. We also found that the CorPI-A copy contains informative phylogenetic information when compared across species of Cornus.     

Sporadic amplification of ID elements in rodents

ID sequences are members of a short interspersed element (SINE) repetitive DNA family within the rodent genome. The copy number of individual ID elements varies by up to three orders of magnitude between species. This amplification has been highly sporadic in the order Rodentia and does not follow any phylogenetic trend. Using library screening and dot-blot analysis, we estimate there are 25,000 copies of ID elements in the deer mouse, 1,500 copies in the gerbil (both cricetid rodents), and 60,000 copies of either ID or ID-like elements in a sciurid rodent (squirrel). By dot-blot analysis, we estimate there are 150,000, 4,000, 1,000, and 200 copies of ID elements in the rat, mouse, hamster, and guinea pig, respectively (which is consistent with previous reports) and 200 copies in the hystricognath rodent, nutria. Therefore, a rapid amplification took place not only after the divergence of rat and mouse but also following the deer mouse (Peromyscus) and hamster split, with no evidence of increased amplifications in hystricognath rodents. No notable variations of sequences from the BC1 genes of several myomorphic rodents were observed that would possibly explain the varied levels of ID amplification. We did observe subgenera and species-group-specific variation in the ID core sequence of the BC1 gene within the genus Peromyscus. Sequence analysis of cloned ID elements in Peromyscus show most ID elements in this genus arose prior to Peromyscus subgenus divergence. Correspondence of the consensus sequence of individual ID elements in gerbil and deer mouse further confirms BC1 as a master gene in ID amplification. Several possible mechanisms responsible for the quantitative variations are explored.The nucleotide sequences reported in this paper have been submitted to the GenBank/EMBL Data Bank with accession numbers: U33850, U33851, U33852 (BC1 sequences); and U33853, U33854, U33855, U33856, U33857, U33858, U33859, U33860, U33861, U33862, U33863, U33864, U33865, U33866, U33867 (ID sequences) Correspondence to: D.H. Kass