Further examples of evolution by gene duplication revealed through DNA sequence comparisons

A degenerative process underlying hierarchic transitions in evolution

This paper describes an evolutionary process likely involved in hierarchic transitions in biological evolution at many levels, from genetics to social organization. It is related to the evolutionary process described as contingent neutral evolution (CNE). It involves a sequence of stages initiated by the spontaneous appearance of functional redundancy. This redundancy can be the result of gene duplication, symbiosis, cell-cell interactions, environmental supports, etc. The availability of redundant sources of biological functionality relaxes purifying selection and allows degenerative changes to accumulate in one or more of the duplicates, potentially degrading or otherwise fractionating its function. This degeneration will be effectively neutral so long as another maintains functional integrity. Sexual recombination can potentially sample different combinations of these sub functional alternatives, with the result that favorable synergistic interactions between independently degenerate duplicates will have a non-negligible probability of being uncovered. The expression of such a synergistic combinatorial effect will result in the irreversible degradation of any remaining autonomous functionality, thereby initiating selection to prevent breakup of co-dependency. This becomes relevant to the evolution of hierarchic transitions when two or more organisms reciprocally duplicate functions that each other requires. If the resulting relaxation of selection reliably persists for an extended evolutionary period it will tend to produce complementary degenerative effects in each organism, leading to their irreversible codependency and purifying selection to avoid loss of integrity of their higher order functional unity. This provides a partial inversion of Darwinian logic that explains how the potential costs of the loss of organism autonomy can be mitigated, enabling the incremental transition to a synergistic higher order unit of evolution.

New Transcriptome-Based SNP Markers for Noug (Guizotia abyssinica) and Their Conversion to KASP Markers for Population Genetics Analyses

The development and use of genomic resources are essential for understanding the population genetics of crops for their efficient conservation and enhancement. Noug (Guizotia abyssinica) is an economically important oilseed crop in Ethiopia and India. The present study sought to develop new DNA markers for this crop. Transcriptome sequencing was conducted on two genotypes and 628 transcript sequences containing 959 single nucleotide polymorphisms (SNPs) were developed. A competitive allele-specific PCR (KASP) assay was developed for the SNPs and used for genotyping of 24 accessions. A total of 554 loci were successfully genotyped across the accessions, and 202 polymorphic loci were used for population genetics analyses. Polymorphism information content (PIC) of the loci varied from 0.01 to 0.37 with a mean of 0.24, and about 49% of the loci showed significant deviation from the Hardy-Weinberg equilibrium. The mean expected heterozygosity was 0.27 suggesting moderately high genetic variation within accessions. Low but significant differentiation existed among accessions (FST = 0.045, p < 0.0001). Landrace populations from isolated areas may have useful mutations and should be conserved and used in breeding this crop. The genomic resources developed in this study were shown to be useful for population genetics research and can also be used in, e.g., association genetics.

Contrasting Patterns of Nucleotide Substitution Rates Provide Insight into Dynamic Evolution of Plastid and Mitochondrial Genomes of Geranium

Geraniaceae have emerged as a model system for investigating the causes and consequences of variation in plastid and mitochondrial genomes. Incredible structural variation in plastid genomes (plastomes) and highly accelerated evolutionary rates have been reported in selected lineages and functional groups of genes in both plastomes and mitochondrial genomes (mitogenomes), and these phenomena have been implicated in cytonuclear incompatibility. Previous organelle genome studies have included limited sampling of Geranium, the largest genus in the family with over 400 species. This study reports on rates and patterns of nucleotide substitutions in plastomes and mitogenomes of 17 species of Geranium and representatives of other Geraniaceae. As detected across other angiosperms, substitution rates in the plastome are 3.5 times higher than the mitogenome in most Geranium. However, in the branch leading to Geranium brycei/Geranium incanum mitochondrial genes experienced significantly higher dN and dS than plastid genes, a pattern that has only been detected in one other angiosperm. Furthermore, rate accelerations differ in the two organelle genomes with plastomes having increased dN and mitogenomes with increased dS. In the Geranium phaeum/Geranium reflexum clade, duplicate copies of clpP and rpoA genes that experienced asymmetric rate divergence were detected in the single copy region of the plastome. In the case of rpoA, the branch leading to G. phaeum/G. reflexum experienced positive selection or relaxation of purifying selection. Finally, the evolution of acetyl-CoA carboxylase is unusual in Geraniaceae because it is only the second angiosperm family where both prokaryotic and eukaryotic ACCases functionally coexist in the plastid.

Origin and evolution of GATA2a and GATA2b in teleosts: insights from tongue sole, Cynoglossus semilaevis

Background. Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication occurred in the lineage of teleost fish and is considered to be responsible for much of the biological diversification within the lineage. GATA2, a member of GATA family of transcription factors, is an important regulator of gene expression in hematopoietic cell in mammals, yet the role of this gene or its putative paralogs in ray-finned fishes remains relatively unknown.

Methods. In this study, we attempted to identify GATA2 sequences from the transcriptomes and genomes of multiple teleosts using the bioinformatic tools MrBayes, MEME, and PAML. Following identification, comparative analysis of genome structure, molecular evolution rate, and expression by real-time qPCR were used to predict functional divergence of GATA2 paralogs and their relative transcription in organs of female and male tongue soles (Cynoglossus semilaevis).

Results. Two teleost GATA2 genes were identified in the transcriptomes of tongue sole and Japanese flounder (Paralichthysolivaceus). Synteny and phylogenetic analysis confirmed that the two genes likely originated from the teleost-specific genome duplication . Additionally, selection pressure analysis predicted these gene duplicates to have undergone purifying selection and possible divergent new functions. This was supported by differential expression pattern of GATA2a and GATA2b observed in organs of female and male tongue soles.

Discussion. Our results indicate that two GATA2 genes originating from the first teleost-specific genome duplication have remained transcriptionally active in some fish species and have likely undergone neofunctionalization. This knowledge provides novel insights into the evolution of the teleost GATA2 genes and constituted important groundwork for further research on the GATA gene family.

Evolutionary Analyses and Natural Selection of Betaine-Homocysteine S-Methyltransferase (BHMT) and BHMT2 Genes

Betaine-homocysteine S-methyltransferase (BHMT) and BHMT2 convert homocysteine to methionine using betaine and S-methylmethionine, respectively, as methyl donor substrates. Increased levels of homocysteine in blood are associated with cardiovascular disease. Given their role in human health and nutrition, we identified BHMT and BHMT2 genes and proteins from 38 species of deuterostomes including human and non-human primates. We aligned the genes to look for signatures of selection, to infer evolutionary rates and events across lineages, and to identify the evolutionary timing of a gene duplication event that gave rise to two genes, BHMT and BHMT2. We found that BHMT was present in the genomes of the sea urchin, amphibians, reptiles, birds and mammals; BHMT2 was present only across mammals. BHMT and BHMT2 were present in tandem in the genomes of all monotreme, marsupial and placental species examined. Evolutionary rates were accelerated for BHMT2 relative to BHMT. Selective pressure varied across lineages, with the highest dN/dS ratios for BHMT and BHMT2 occurring immediately following the gene duplication event, as determined using GA Branch analysis. Nine codons were found to display signatures suggestive of positive selection; these contribute to the enzymatic or oligomerization domains, suggesting involvement in enzyme function. Gene duplication likely occurred after the divergence of mammals from other vertebrates but prior to the divergence of extant mammalian subclasses, followed by two deletions in BHMT2 that affect oligomerization and methyl donor specificity. The faster evolutionary rate of BHMT2 overall suggests that selective constraints were reduced relative to BHMT. The dN/dS ratios in both BHMT and BHMT2 was highest following the gene duplication, suggesting that purifying selection played a lesser role as the two paralogs diverged in function.

Effects of atrazine and chlorpyrifos on the mRNA levels of HSP70 and HSC70 in the liver, brain, kidney and gill of common carp (Cyprinus carpio L.)

Although the chaperone role of heat shock proteins (HSPs) has been demonstrated in invertebrates, the function of HSPs in vertebrates, especially in fish, remains unclear. In this study, relative changes in the mRNA abundance of the HSP70 gene were examined by real-time PCR in the muscle, spleen, head kidney, heart, liver, brain, kidney and gill of common carp. Results indicated that the highest and lowest levels of HSP70 expression were found in the heart and muscle, respectively, and the highest and lowest levels of HSC70 expression were found in the spleen and muscle, respectively. In addition, we investigated differential HSP70 gene expression in common carp after a 40-d exposure to chlorpyrifos (CPF) and atrazine (ATR), alone or in combination, and after a 20-d recovery. Results indicated that the expression of carp HSP70 and 70-kDa heat shock cognate protein (HSC70) with ATR and CPF treatment alone or in combination was significantly upregulated. The present results provide new insights into the mechanisms used by fish to adapt to stressful environments.

Differential evolution of MAGE genes based on expression pattern and selection pressure

Starting from publicly-accessible datasets, we have utilized comparative and phylogenetic genome analyses to characterize the evolution of the human MAGE gene family. Our characterization of genomic structures in representative genomes of primates, rodents, carnivora, and macroscelidea indicates that both Type I and Type II MAGE genes have undergone lineage-specific evolution. The restricted expression pattern in germ cells of Type I MAGE orthologs is observed throughout evolutionary history. Unlike Type II MAGEs that have conserved promoter sequences, Type I MAGEs lack promoter conservation, suggesting that epigenetic regulation is a central mechanism for controlling their expression. Codon analysis shows that Type I but not Type II MAGE genes have been under positive selection. The combination of genomic and expression analysis suggests that Type 1 MAGE promoters and genes continue to evolve in the hominin lineage, perhaps towards functional diversification or acquiring additional specific functions, and that selection pressure at codon level is associated with expression spectrum.

Novel polymorphism of internal transcribed spacers (ITS) and their utilization in phylogenetic analysis of Neanthes glandicincta (Annelida: Polychaeta: Nereididae)

Sequences of internal transcribed spacers (ITS1 and ITS2) are increasingly being used to infer phylogenetic relationships at or below species levels. Here we report a novel case of ITS polymorphism within Neanthes glandicincta (Annelida: Polychaeta: Nereididae). Two types of ITS sequence (Type I and Type II) were cloned and sequenced, which showed significant differences both in nucleotide composition and length. Variations of these two types sequences also differed from each other with Type I was highly divergent while Type II was highly conserved. Phylogenetic trees inferred from ITS1 and ITS2 sequences showed striking discrepancy in N. glandicincta. Non-concerted evolution of multi-gene is suggested to be responsible for the high degree of polymorphism in ITS regions. Due to the two divergent types of ITS presented within a single N. glandicincta individual, the utilization of ITS regions for delineation of population or closely related species cannot be substantiated. The finding of different types of ITS in a single individual also stresses the need for analyzing a large number of clones whenever ITS sequences obtained by PCR amplification and cloning are being used in phylogenetic reconstruction.

Colloquium paper: a role for relaxed selection in the evolution of the language capacity

Explaining the extravagant complexity of the human language and our competence to acquire it has long posed challenges for natural selection theory. To answer his critics, Darwin turned to sexual selection to account for the extreme development of language. Many contemporary evolutionary theorists have invoked incredibly lucky mutation or some variant of the assimilation of acquired behaviors to innate predispositions in an effort to explain it. Recent evodevo approaches have identified developmental processes that help to explain how complex functional synergies can evolve by Darwinian means. Interestingly, many of these developmental mechanisms bear a resemblance to aspects of Darwin's mechanism of natural selection, often differing only in one respect (e.g., form of duplication, kind of variation, competition/cooperation). A common feature is an interplay between processes of stabilizing selection and processes of relaxed selection at different levels of organism function. These may play important roles in the many levels of evolutionary process contributing to language. Surprisingly, the relaxation of selection at the organism level may have been a source of many complex synergistic features of the human language capacity, and may help explain why so much language information is "inherited" socially.

Molecular evolution of the duplicated TFIIAgamma genes in Oryzeae and its relatives

Background

Gene duplication provides raw genetic materials for evolutionary novelty and adaptation. The evolutionary fate of duplicated transcription factor genes is less studied although transcription factor gene plays important roles in many biological processes. TFIIAγ is a small subunit of TFIIA that is one of general transcription factors required by RNA polymerase II. Previous studies identified two TFIIAγ-like genes in rice genome and found that these genes either conferred resistance to rice bacterial blight or could be induced by pathogen invasion, raising the question as to their functional divergence and evolutionary fates after gene duplication.

Results

We reconstructed the evolutionary history of the TFIIAγ genes from main lineages of angiosperms and demonstrated that two TFIIAγ genes (TFIIAγ1 and TFIIAγ5) arose from a whole genome duplication that happened in the common ancestor of grasses. Likelihood-based analyses with branch, codon, and branch-site models showed no evidence of positive selection but a signature of relaxed selective constraint after the TFIIAγ duplication. In particular, we found that the nonsynonymous/synonymous rate ratio (ω = d_N/d_S) of the TFIIAγ1 sequences was two times higher than that of TFIIAγ5 sequences, indicating highly asymmetric rates of protein evolution in rice tribe and its relatives, with an accelerated rate of TFIIAγ1 gene. Our expression data and EST database search further indicated that after whole genome duplication, the expression of TFIIAγ1 gene was significantly reduced while TFIIAγ5 remained constitutively expressed and maintained the ancestral role as a subunit of the TFIIA complex.

Conclusion

The evolutionary fate of TFIIAγ duplicates is not consistent with the neofunctionalization model that predicts that one of the duplicated genes acquires a new function because of positive Darwinian selection. Instead, we suggest that subfunctionalization might be involved in TFIIAγ evolution in grasses. The fact that both TFIIAγ1 and TFIIAγ5 genes were effectively involved in response to biotic or abiotic factors might be explained by either Dykhuizen-Hartl effect or buffering hypothesis.

Adaptive evolution of recently duplicated accessory gland protein genes in desert Drosophila

The relationship between animal mating system variation and patterns of protein polymorphism and divergence is poorly understood. Drosophila provides an excellent system for addressing this issue, as there is abundant interspecific mating system variation. For example, compared to D. melanogaster subgroup species, repleta group species have higher remating rates, delayed sexual maturity, and several other interesting differences. We previously showed that accessory gland protein genes (Acp's) of Drosophila mojavensis and D. arizonae evolve more rapidly than Acp's in the D. melanogaster subgroup and that adaptive Acp protein evolution is likely more common in D. mojavensis/D. arizonae than in D. melanogaster/D. simulans. These findings are consistent with the idea that greater postcopulatory selection results in more adaptive evolution of seminal fluid proteins in the repleta group flies. Here we report another interesting evolutionary difference between the repleta group and the D. melanogaster subgroup Acp's. Acp gene duplications are present in D. melanogaster, but their high sequence divergence indicates that the fixation rate of duplicated Acp's has been low in this lineage. Here we report that D. mojavensis and D. arizonae genomes contain several very young duplicated Acp's and that these Acp's have experienced very rapid, adaptive protein divergence. We propose that rapid remating of female desert Drosophila generates selection for continuous diversification of the male Acp complement to improve male fertilization potential. Thus, mating system variation may be associated with adaptive protein divergence as well as with duplication of Acp's in Drosophila.

The alphaD-globin gene originated via duplication of an embryonic alpha-like globin gene in the ancestor of tetrapod vertebrates

Gene duplication is thought to play an important role in the co-option of existing protein functions to new physiological pathways. The globin superfamily of genes provides an excellent example of the kind of physiological versatility that can be attained through the functional and regulatory divergence of duplicated genes that encode different subunit polypeptides of the tetrameric hemoglobin protein. In contrast to prevailing views about the evolutionary history of the alpha-globin gene family, here we present phylogenetic evidence that the alpha(A)- and alpha(D)-globin genes are not the product of a single, tandem duplication of an ancestral globin gene with adult function in the common ancestor of extant birds, reptiles, and mammals. Instead, our analysis reveals that the alpha(D)-globin gene of amniote vertebrates arose via duplication of an embryonic alpha-like globin gene that predated the radiation of tetrapods. The important evolutionary implication is that the distinct biochemical properties of alpha(D)-hemoglobin (HbD) are not exclusively derived characters that can be attributed to a post-duplication process of neofunctionalization. Rather, many of the distinct biochemical properties of HbD are retained ancestral characters that reflect the fact that the alpha(D)-globin gene arose via duplication of a gene that had a larval/embryonic function. These insights into the evolutionary origin of HbD illustrate how adaptive modifications of physiological pathways may result from the retention and opportunistic co-option of ancestral protein functions.

The evolutionary fate of recently duplicated retrogenes in mice

Inferences about the evolutionary impact of gene duplications often rely on the analysis of their long-term outcome. The fate of the majority of them must, however, be decided shortly after duplication. Here we analysed the evolutionary pattern of 10 mouse genes very recently duplicated by retrotransposition, by sequencing the retroposed copy in five to 10 closely related mouse species. In all cases the retroposed copy experienced accelerated nonsynonymous evolution whereas the divergence pattern of the source copy appeared unaffected by the duplication, consistent with the neofunctionalization model. The analysis further revealed that most retrogenes, including pseudogenes, did not experience a period of relaxed neutral evolution, but have been submitted to purifying selection ever since their retroposition. We propose that these duplicates play a biochemical role but are not indispensable. Purifying selection prevents them from acquiring a negative role until they are lost or silenced. This period of unnecessary redundancy could in rare cases give the time for new functions to evolve.

Multiple copies of coding as well as pseudogene c-mos sequence exist in three lacertid species

The analysis of a 581 bp section of the nuclear gene c-mos revealed multiple copies of putative functional sequences as well as pseudogenes in three closely related lacertid species Lacerta laevis, L. kulzeri and L. cyanisparsa. A phylogenetic analysis of c-mos in comparison with a molecular phylogeny based on the mitochondrial cytochrome b gene supports our findings. The study also provides new insights into the phylogenetic relationships of L. cyanisparsa and L. laevis. Pseudogenes of the three species share 11 single-nucleotide substitutions, a 1 bp deletion and a premature stop codon but differ by group-specific mutations. This result suggests that the c-mos gene has become duplicated and subsequently silenced already in the common ancestor of the three species. Sequence divergence suggests that the duplication and the loss of function occurred in the late Miocene/early Pliocene, i.e., about 5 million years ago. Indications of gene conversion are discussed. We suggest that future studies using c-mos for phylogenetic studies should provide evidence for the orthology of the sequences compared.

Independent Duplications of the Acetylcholinesterase Gene Conferring Insecticide Resistance in the Mosquito Culex pipiens

Gene duplication is thought to be the main potential source of material for the evolution of new gene functions. Several models have been proposed for the evolution of new functions through duplication, most based on ancient events (Myr). We provide molecular evidence for the occurrence of several (at least 3) independent duplications of the ace-1 locus in the mosquito Culex pipiens, selected in response to insecticide pressure that probably occurred very recently (<40 years ago). This locus encodes the main target of several insecticides, the acetylcholinesterase. The duplications described consist of 2 alleles of ace-1, 1 susceptible and 1 resistant to insecticide, located on the same chromosome. These events were detected in different parts of the world and probably resulted from distinct mechanisms. We propose that duplications were selected because they reduce the fitness cost associated with the resistant ace-1 allele through the generation of persistent, advantageous heterozygosis. The rate of duplication of ace-1 in C. pipiens is probably underestimated, but seems to be rather high.

Sequencing and expression pattern of inducible heat shock gene products in the European flat oyster, Ostrea edulis

Heat shock proteins are a multigene family of polypeptides composed of the constitutively-expressed heat shock cognate (HSC) members and of the stress-inducible (HSP) proteins, whose expression is specifically induced by stress factors. Constitutive and inducible 70 kDa isoforms are reported in vertebrates and invertebrates. HSCs are expressed in all bivalve molluscs studied to date, while the occurrence of strictly heat-inducible HSPs seems a distinctive feature of oysters. To gain more insight into the molecular features of the Ostrea edulis HSP70, we have cloned and sequenced the gene product putatively encoding for the heat-inducible isoform HSP69 and examined the pattern of expression after heat exposure. Four different clones of approximately 1794 bp were obtained that share a high degree of homology with heat-inducible HSP70 from other bivalves. Amino acid sequence comparisons indicated that the main structural features of the heat-inducible HSP70 are highly conserved within the O. edulis HSP70 clones, while lower sequence homology occurred with respect to HSC70 transcripts. Northern blot analysis indicated that HSP69 mRNA was absent in control animals but induced after heat shock (1 h at 32 degrees C or higher). Induction was detectable immediately after heat shock, reaching a maximum after 2 to 3 h of post-stress recovery at 18 degrees C, and decreasing thereafter. A phylogenetic analysis of the HSP70 family members from oyster and other bivalves revealed a substantial conservation in the evolutionary pattern among constitutive and inducible gene products, from invertebrates to higher vertebrates.

Differential HSP70 gene expression in the Mediterranean mussel exposed to various stressors

HSP70 gene expression was studied by quantitative RT-PCR after cloning and sequencing of two different HSP70 gene fragments from the digestive gland of Mytilus galloprovincialis, called MgHSP70 and MgHSC70. Heat shock (1h at 35 degrees C) caused rapid induction of MgHSP70, while no change was observed for MgHSC70. Hg(2+) (150 microg/L for different time periods) significantly induced MgHSP70 expression that reached maximum levels after 24h, decreasing thereafter. MgHSC70 expression was inhibited after 1 day and induced after a 6-day exposure to Hg(2+). A 1-week exposure to Cr(6+) (1, 10, and 50 ng/L) induced and inhibited MgHSC70 and MgHSP70 transcript levels, respectively. MgHSC70 and MgHSP70 appear to play different roles in cell protection; the former is induced after acute stress and/or during the earlier phase of the response while the latter is induced by chronic stress. The present results provide new insights into mechanisms used by mussels to adapt to stressful environments.

Actin phylogeny and intron distribution in bangiophyte red algae(rhodoplantae)

The molecular phylogeny of red algal actin genes, with emphasis on the paraphyletic "Bangiophyceae," was examined and compared to the rhodophyte SSU rDNA phylogeny. Nineteen new genomic actin sequences and seven SSU rDNA sequences were obtained and subjected to diverse phylogenetic analyses (maximum likelihood, distance/neighbor-joining, maximum parsimony, Bayesian analyses, and, with respect to protein sequences, also quartet puzzling). The actin trees confirmed most of the major clades found in the SSU rDNA phylogenies, although with a lower resolution. An actin gene duplication in the florideophycean lineage is reported, presumably related to an increased complexity of sexual reproduction. In addition, the distribution and characteristics of spliceosomal introns found in some of the actin sequences were examined. Introns were found in almost all florideophycean actin genes, whereas only two bangiophyte sequences contained introns. One intron in the florideophycean actin genes was also found in metazoan, and, shifted by one or two nucleotides, in a glaucocystophyte, a cryptophyte, and two fungal actin genes, and thus may be an ancient intron.

Rapid evolution through gene duplication and subfunctionalization of the testes-specific alpha4 proteasome subunits in Drosophila

Gene duplication is an important mechanism for acquiring new genes and creating genetic novelty in organisms. Evidence suggests that duplicated genes are retained at a much higher rate than originally thought and that functional divergence of gene copies is a major factor promoting their retention in the genome. We find that two Drosophila testes-specific alpha4 proteasome subunit genes (alpha4-t1 and alpha4-t2) have a higher polymorphism within species and are significantly more diverged between species than the somatic alpha4 gene. Our data suggest that following gene duplication, the alpha4-t1 gene experienced relaxed selective constraints, whereas the alpha4-t2 gene experienced positive selection acting on several codons. We report significant heterogeneity in evolutionary rates among all three paralogs at homologous codons, indicating that functional divergence has coincided with genic divergence. Reproductive subfunctionalization may allow for a more rapid evolution of reproductive traits and a greater specialization of testes function. Our data add to the increasing evidence that duplicated genes experience lower selective constraints and in some cases positive selection following duplication. Newly duplicated genes that are freer from selective constraints may provide a mechanism for developing new interactions and a pathway for the evolution of new genes.

Birth-and-death evolution of the Cecropin multigene family in Drosophila

Cecropins are insect antibacterial peptides that are part of the insect humoral immune response and could, therefore, be potential targets of natural selection. In Drosophila, the Cec genes constitute a multigene family whose members are arranged in tandem. The complete Cec family was isolated in two obscura group species: D. subobscura and D. pseudoobscura. The chromosomal regions encompassing the Cec genes were subsequently sequenced and mapped by in situ hybridization. In D. pseudoobscura, as in species of the D. melanogaster complex and in D. virilis, the Cec genes constitute a single cluster with five genes. In D. subobscura, unlike in the rest of the species, the eight members of the family are split into two clusters located in different parts of the same chromosome. Remarkable differences in levels of divergence were observed between copies in both species. The genomic organization and the phylogenetic relationships among members of this family in the genus Drosophila indicate (i) that the presence of two clusters is the derived state of the family, (ii) repeated gene duplication within species, (iii) nonfunctionalization and loss of some Cec copies, and (iv) the presence of both highly divergent and highly similar copies within species. These features are better explained by the birth-and-death model of molecular evolution, which posits that the high instability of the Cec multigene family in Drosophila is simply determined by the duplication rate and the subsequent loss of duplicated loci.

Asymmetric evolution of duplicate genes encoding the CCAAT-binding factor NF-Y in plant genomes

NF-Y is a ubiquitous CCAAT-binding factor composed of NF-YA, NF-YB and NF-YC. Multiple genes encoding NF-Y subunits have been identified in plant genomes. It remains unclear whether the duplicate genes underwent different evolutionary patterns. Likelihood-ratio tests were used to examine whether the amino acid substitution rates are the same between duplicate genes. The influences of selection on evolution were evaluated by comparing the conservative and radical amino acid substitution rates, as well as maximum-likelihood analysis. Some NF-YB and NF-YC duplicates showed significant evidence of asymmetric evolution but not the NF-YA duplicates. Most amino acid replacements in the NF-YB and NF-YC duplicates result in changes in hydropathy, polar requirement and polarity. The physicochemical changes in the sequences of NF-YB seem to be coupled to asymmetric divergence in gene function. Plant NF-Y genes have evolved in different patterns. Relaxed selective constraints following gene duplication are most likely responsible for the unequal evolutionary rates and distinct divergence patterns of duplicate NF-Y genes. Positive selection may have promoted amino acid hydropathy changes in the NF-YC duplicates.

Rapid subfunctionalization accompanied by prolonged and substantial neofunctionalization in duplicate gene evolution

Gene duplication is the primary source of new genes. Duplicate genes that are stably preserved in genomes usually have divergent functions. The general rules governing the functional divergence, however, are not well understood and are controversial. The neofunctionalization (NF) hypothesis asserts that after duplication one daughter gene retains the ancestral function while the other acquires new functions. In contrast, the subfunctionalization (SF) hypothesis argues that duplicate genes experience degenerate mutations that reduce their joint levels and patterns of activity to that of the single ancestral gene. We here show that neither NF nor SF alone adequately explains the genome-wide patterns of yeast protein interaction and human gene expression for duplicate genes. Instead, our analysis reveals rapid SF, accompanied by prolonged and substantial NF in a large proportion of duplicate genes, suggesting a new model termed subneofunctionalization (SNF). Our results demonstrate that enormous numbers of new functions have originated via gene duplication.

Patterns of diversifying selection in the phytotoxin-like scr74 gene family of Phytophthora infestans

Phytophthora infestans, the organism responsible for the Irish famine, causes late blight, a re-emerging disease of potato and tomato. Little is known about the molecular evolution of P. infestans genes. To identify candidate effector genes (virulence or avirulence genes) that may have co-evolved with the host, we mined expressed sequence tag (EST) data from infection stages of P. infestans for secreted and potentially polymorphic genes. This led to the identification of scr74, a gene that encodes a predicted 74-amino acid secreted cysteine-rich protein with similarity to the Phytophthora cactorum phytotoxin PcF. The expression of scr74 was upregulated approximately 60-fold 2 to 4 days after inoculation of tomato and was also significantly induced during early stages of colonization of potato. The scr74 gene was found to belong to a highly polymorphic gene family within P. infestans with 21 different sequences identified. Using the approximate and maximum likelihood (ML) methods, we found that diversifying selection likely caused the extensive polymorphism observed within the scr74 gene family. Pairwise comparisons of 17 scr74 sequences revealed elevated ratios of nonsynonymous to synonymous nucleotide-substitution rates, particularly in the mature region of the proteins. Using ML, all 21 polymorphic amino acid sites were identified to be under diversifying selection. Of these 21 amino acids, 19 are located in the mature protein region, suggesting that selection may have acted on the functional portions of the proteins. Further investigation of gene copy number and organization revealed that the scr74 gene family comprises at least three copies located in a region of no more than 300 kb of the P. infestans genome. We found evidence that recombination contributed to sequence divergence within at least one gene locus. These results led us to propose an evolutionary model that involves gene duplication and recombination, followed by functional divergence of scr74 genes. This study provides support for using diversifying selection as a criterion for identifying candidate effector genes from sequence databases.

Mining EST databases to resolve evolutionary events in major crop species

Using plant EST collections, we obtained 1392 potential gene duplicates across 8 plant species: Zea mays, Oryza sativa, Sorghum bicolor, Hordeum vulgare, Solanum tuberosum, Lycopersicon esculentum, Medicago truncatula, and Glycine max. We estimated the synonymous and nonsynonymous distances between each gene pair and identified two to three mixtures of normal distributions corresponding to one to three rounds of genome duplication in each species. Within the Poaceae, we found a conserved duplication event among all four species that occurred approximately 50-60 million years ago (Mya); an event that probably occurred before the major radiation of the grasses. In the Solanaceae, we found evidence for a conserved duplication event approximately 50-52 Mya. A duplication in soybean occurred approximately 44 Mya and a duplication in Medicago about 58 Mya. Comparing synonymous and nonsynonymous distances allowed us to determine that most duplicate gene pairs are under purifying, negative selection. We calculated Pearson's correlation coefficients to provide us with a measure of how gene expression patterns have changed between duplicate pairs, and compared this across evolutionary distances. This analysis showed that some duplicates seemed to retain expression patterns between pairs, whereas others showed uncorrelated expression.

Gene expression and feeding ecology: evolution of piscivory in the venomous gastropod genus Conus

Differential expression of gene-family members is typically associated with the specific development of certain tissues and organs, but its importance in the ecological adaptation of organisms has rarely been investigated. Several specialized feeding modes have evolved within the predatory marine gastropod genus Conus, including molluscivory and piscivory. Based on phylogenetic investigations of Conus species, it has been concluded that piscivory arose at least twice in this genus. Moreover, molecular analyses of conotoxin mRNA transcripts reveal that piscivores from independent evolutionary lineages express the same subset of four-loop conotoxins, contrary to phylogenetic expectations. These results demonstrate that differential expression of gene-family members can play a key role in adaptive evolution, particularly during shifts to new ecological niches.

Evolutionary dynamics of the DNA-binding domains in putative R2R3-MYB genes identified from rice subspecies indica and japonica genomes

The molecular evolution of the R2R3-MYB gene family is of great interest because it is one of the most important transcription factor gene families in the plant kingdom. Comparative analyses of a gene family may reveal important adaptive changes at the protein level and thereby provide insights that relate structure to function. We have performed a range of comparative and bioinformatics analyses on R2R3-MYB genes identified from the rice (Oryza sativa subsp. japonica and indica) and Arabidopsis genome sequences. The study provides an initial framework to investigate how different evolutionary lineages in a gene family evolve new functions. Our results reveal a remarkable excess of non-synonymous substitutions, an indication of adaptive selection on protein structure that occurred during the evolution of both helix1 and helix2 of rice R2R3-MYB DNA-binding domains. These flexible alpha-helix regions associated with high frequencies of excess non-synonymous substitutions may play critical roles in the characteristic packing of R2R3-MYB DNA-binding domains and thereby modify the protein-DNA interaction process resulting in the recognition of novel DNA-binding sites. Furthermore, a co-evolutionary pattern is found between the second alpha-helix of the R2 domain and the second alpha-helix of the R3 domain by examining all the possible alpha-helix pairings in both the R2 and R3 domains. This points to the functional importance of pairing interactions between related secondary structures.

Rapid evolution of a pollen-specific oleosin-like gene family from Arabidopsis thaliana and closely related species

It has been shown in a variety of species that genes expressed in reproductive tissues evolve rapidly, which often appears to be the result of positive Darwinian selection. We investigated the evolution of a family of seven pollen-specific oleosin-like proteins (or oleopollenins) in Arabidopsis thaliana and two closely related species. More than 30 kb of a genomic region that harbors the complete, tandemly repeated oleopollenin cluster were sequenced from Arabidopsis lyrata ssp. lyrata, and Boechera drummondii. A phylogenetic analysis of the complete gene cluster from these three species and from Brassica oleracea confirmed its rapid evolution resulting from gene duplication and gene loss events, numerous amino acid substitutions, and insertions/deletions in the coding sequence. Independent duplications were inferred in the lineages leading to Arabidopsis and to Brassica, and gene loss was inferred in the lineage leading to B. drummondii. Comparisons of the ratio of nonsynonymous (d(N)) and synonymous (d(S)) divergence revealed that the oleopollenins are among the most rapidly evolving proteins currently known from Arabidopsis and that they may evolve under positive Darwinian selection. Reverse transcriptase polymerase chain reaction analysis demonstrated the expression of oleopollenins in flowers of the outcrossing A. lyrata, the selfing B. drummondii, and the apomictic Boechera holboellii, suggesting that oleopollenins play an important role in species with different breeding systems. These results are consistent with a putative function in species recognition, but further analyses of protein function and sequence variation in species with different breeding systems are necessary to reveal the underlying causes for the rapid evolution of oleopollenins.

The origin of new genes: glimpses from the young and old

Genome data have revealed great variation in the numbers of genes in different organisms, which indicates that there is a fundamental process of genome evolution: the origin of new genes. However, there has been little opportunity to explore how genes with new functions originate and evolve. The study of ancient genes has highlighted the antiquity and general importance of some mechanisms of gene origination, and recent observations of young genes at early stages in their evolution have unveiled unexpected molecular and evolutionary processes.

Excess non-synonymous substitutions suggest that positive selection episodes occurred during the evolution of DNA-binding domains in the Arabidopsis R2R3-MYB gene family

It has been suggested that evolutionary changes in regulatory genes may be the predominant molecular mechanism governing both physiological and morphological evolution. R2R3-AtMYB is one of the largest transcription factor gene families in Arabidopsis. Using inferred ancestral sequences we show that several lineages in the R2R3-AtMYB phylogeny experienced excess non-synonymous nucleotide substitution upon gene duplication, indicating episodes of positive selection driving adaptive shifts early in the evolution of this gene family. A noise reduction technique was then used to determine individual sites in DNA-binding domains (R2 domain and R3 domain) of R2R3-AtMYB protein sequence that were favored by frequent non-synonymous substitutions. The analyses reveal that the first helix (helix1) and the second helix (helix2) in both R2 and R3 domains are characterized by more frequent non-synonymous substitutions, and thus experienced significantly higher positive selection pressure than the third helix (helix3) in both domains. Previous MYB protein structure studies have suggested that helix1 and helix2 in both R2 and R3 domains are involved in the characteristic packing of R2R3-AtMYB DNA-binding domains. This suggests that excess non-synonymous substitutions in these helices could have resulted in MYB recognition of novel gene target sites.

Maximum likelihood methods for detecting adaptive evolution after gene duplication

The rapid accumulation of genomic sequences in public databases will finally allow large scale studies of gene family evolution, including evaluation of the role of positive Darwinian selection following a duplication event. This will be possible because recent statistical methods of comparing synonymous and nonsynonymous substitution rates permit reliable detection of positive selection at individual amino acid sites and along evolutionary lineages. Here, we summarize maximum-likelihood based methods, and present a framework for their application to analysis of gene families. Using these methods, we investigated the role of positive Darwinian selection in the ECP-EDN gene family of primates and the Troponin C gene family of vertebrates. We also comment on the limitations of these methods and discuss directions for further improvements.

Variant subunit specificity in the quaternary structure of Artemia hemoglobin

The brine shrimp Artemia has three extracellular hemoglobins (Hbs) that are developmentally expressed and exhibit distinct oxygen-binding characteristics (Heip, Moens, and Kondo 1978; Heip et al. 1978 ). These Hbs are composed of two polymers, each of which comprises nine covalently linked globin domains. Although the cDNA sequences of two nine-domain globins from Artemia have been published, there is evidence for the existence of further expressed globin genes (Manning, Trotman, and Tate 1990 ). In the present study extensive analysis at the cDNA and genomic levels was performed in order to determine the globin gene copy number in Artemia. Sequence and Southern analysis suggest that four Hb polymers (T1, T2, C1, and C2) are expressed in Artemia. In addition, there is also at least one globin pseudogene. Protein sequencing of the native Hbs revealed that there are limitations on which two polymers can associate. The composition of the Hbs has been determined to be: Hb I, C1C2; Hb II, C1T2; and Hb III, T1T2. These pairings allow the levels of the three Artemia Hbs to be regulated independently by polymer expression alone, therefore explaining the previously inconsistent developmental and hypoxia-induced expression patterns.

Evolution of novel genes

Much progress in understanding the evolution of new genes has been accomplished in the past few years. Molecular mechanisms such as illegitimate recombination and LINE element mediated 3' transduction underlying exon shuffling, a major process for generating new genes, are better understood. The identification of young genes in invertebrates and vertebrates has revealed a significant role of adaptive evolution acting on initially rudimentary gene structures created as if by evolutionary tinkers. New genes in humans and our primate relatives add a new component to the understanding of genetic divergence between humans and non-humans.

Accelerated regulatory gene evolution in an adaptive radiation

The disparity between rates of morphological and molecular evolution remains a key paradox in evolutionary genetics. A proposed resolution to this paradox has been the conjecture that morphological evolution proceeds via diversification in regulatory loci, and that phenotypic evolution may correlate better with regulatory gene divergence. This conjecture can be tested by examining rates of regulatory gene evolution in species that display rapid morphological diversification within adaptive radiations. We have isolated homologues to the Arabidopsis APETALA3 (ASAP3/TM6) and APETALA1 (ASAP1) floral regulatory genes and the CHLOROPHYLL A/B BINDING PROTEIN9 (ASCAB9) photosynthetic structural gene from species in the Hawaiian silversword alliance, a premier example of plant adaptive radiation. We have compared rates of regulatory and structural gene evolution in the Hawaiian species to those in related species of North American tarweeds. Molecular evolutionary analyses indicate significant increases in nonsynonymous relative to synonymous nucleotide substitution rates in the ASAP3/TM6 and ASAP1 regulatory genes in the rapidly evolving Hawaiian species. By contrast, no general increase is evident in neutral mutation rates for these loci in the Hawaiian species. An increase in nonsynonymous relative to synonymous nucleotide substitution rate is also evident in the ASCAB9 structural gene in the Hawaiian species, but not to the extent displayed in the regulatory loci. The significantly accelerated rates of regulatory gene evolution in the Hawaiian species may reflect the influence of allopolyploidy or of selection and adaptive divergence. The analyses suggest that accelerated rates of regulatory gene evolution may accompany rapid morphological diversification in adaptive radiations.

Positive and negative selection in the DAZ gene family

Because a microdeletion containing the DAZ gene is the most frequently observed deletion in infertile men, the DAZ gene was considered a strong candidate for the azoospermia factor. A recent evolutionary analysis, however, suggested that DAZ was free from functional constraints and consequently played little or no role in human spermatogenesis. The major evidence for this surprising conclusion is that the nonsynonymous substitution rate is similar to the synonymous rate and to the rate in introns. In this study, we reexamined the evolution of the DAZ gene family by using maximum-likelihood methods, which accommodate variable selective pressures among sites or among branches. The results suggest that DAZ is not free from functional constraints. Most amino acids in DAZ are under strong selective constraint, while a few sites are under diversifying selection with nonsynonymous/ synonymous rate ratios (d(N)/d(S)) well above 1. As a result, the average d(N)/d(S) ratio over sites is not a sensible measure of selective pressure on the protein. Lineage-specific analysis indicated that human members of this gene family were evolving by positive Darwinian selection, although the evidence was not strong.

Evolution of gene families

In eukaryote genomes, there are many kinds of gene families. Gene duplication and conversion are sources of the evolution of gene families, including those with uniform members and those with diverse functions. Population genetics theory on identity coefficients among gene members of a gene family shows that the balance between diversification by mutation, and homogenization by unequal crossing over and gene conversion, is important. Also, evolution of new functions is due to gene duplication followed by differentiation. Positive selection is necessary for the evolution of novel functions. However, many examples of current gene families suggest that both drift and selection are at work on their evolution.

Evolutionary diversification of multigene families: allelic selection of toxins in predatory cone snails

In order to investigate the evolution of conotoxin multigene families among two closely related vermivorous CONUS: species, we sequenced 104 four-loop conotoxin mRNAs from two individuals of CONUS: ebraeus and compared these with sequences already obtained from CONUS: abbreviatus. In contrast to the diversity of conotoxin sequences obtained from C. abbreviatus, only two common sequence variants were recovered from C. ebraeus. Segregation patterns of the variants in these two individuals and restriction digests of four-loop conotoxin amplification products from nine additional individuals suggest that the common variants are alleles from a single locus. These two putative alleles differ at nine positions that occur nonrandomly in the toxin-coding region of the sequences. Moreover, all substitutions are at nonsynonymous sites and are responsible for seven amino acid differences among the predicted amino acid sequences of the alleles. These results imply that conotoxin diversity is driven by strong diversifying selection and some form of frequency-dependent or overdominant selection at conotoxin loci, and they suggest that diverse conotoxin multigene families can originate from duplications at polymorphic loci. Furthermore, none of the sequences recovered from C. ebraeus appeared to be orthologs of loci from C. abbreviatus, and attempts to amplify orthologous sequences with locus-specific primers were unsuccessful among these species. These patterns suggest that venoms of closely related CONUS: species may differ due to the differential expression of conotoxin loci.

The origin of the Jingwei gene and the complex modular structure of its parental gene, yellow emperor, in Drosophila melanogaster

Jingwei (jgw) is the first gene found to be of sufficiently recent origin in Drosophila to offer insights into the origin of a gene. While its chimerical gene structure was partially resolved as including a retrosequence of alcohol dehydrogenase (ADH:), the structure of its non-ADH: parental gene, the donor of the N-terminal domain of jgw, is unclear. We characterized this non-ADH: parental locus, yellow emperor (ymp), by cloning it, mapping it onto the polytene chromosomes, sequencing the entire locus, and examining its expression patterns in Drosophila melanogaster. We show that ymp is located in the 96-E region; the N-terminal domain of ymp has donated the non-ADH: portion of jgw via a duplication. The similar 5' portions of the gene and its regulatory sequences give rise to similar testis-specific expression patterns in ymp and jgw in Drosophila teissieri. Furthermore, between-species comparison of ymp revealed purifying selection in the protein sequence, suggesting a functional constraint in ymp. While the structure of ymp provides clear information for the molecular origin of the new gene jgw, it unexpectedly casts a new light on the concept of genes. We found, for the first time, that the single locus of the ymp gene encompasses three major molecular mechanisms determining structure of eukaryotic genes: (1) the 5' exons of ymp are involved in an exon-shuffling event that has created the portion recruited by jgw; (2) using alternative cleavage sites and alternative splicing sites, the 3' exon groups of ymp produce two proteins with nonhomologous C-terminal domains, both exclusively in the testis; and (3) in the opposite strand of the third intron of ymp is an essential gene, musashi (msi), which encodes an RNA-binding protein. The composite gene structure of ymp manifests the complexity of the gene concept, which should be considered in genomic research, e.g., gene finding.

Reliable amplification of actin genes facilitates deep-level phylogeny

The gene for actin as a highly conserved and functionally essential genetic element is developing into a major tool for phylogenetic analysis within a broad organismic range. We therefore propose a set of universally applicable primers that allow reliable amplification of actin genes. For primer construction the amino acid sequences of 57 actin genes comprising fungi, animals, plants and protists were analysed, aligned and used for the definition of six well-conserved regions which are suitable as priming sites in PCR amplification experiments. Ten primers were designed for specific in vitro amplification of actin gene fragments from a wide range of microorganisms. The corresponding gene fragments provide a strong basis to isolate nearly complete actin genes for further molecular characterization and for establishing phylogenies based on actin gene trees.

Phylogeny, duplication, and intraspecific variation of Adh sequences in New World diploid cottons (Gossypium l., malvaceae)

The 13 "D-genome"cotton species are a monophyletic assemblage of morphologically diverse diploids that inhabit arid to semiarid regions in Mexico, with 1 disjunct species each in Peru and the Galapagos Islands and 1 species whose range extends northward into Arizona. While these species lack commercially significant fiber (i. e., cotton), they are important in that they represent one of the parental genomes of the cultivated tetraploid cottons. To assess phylogenetic relationships among these species, we sequenced and analyzed a region of a nuclear-encoded alcohol dehydrogenase gene (AdhA). Phylogenetic analysis resulted in a topology that is generally consistent with current taxonomic alignment of the species, although the phylogeny based on AdhA sequences conflicts with those inferred from cpDNA and ITS data sets, most notably in the position of the anomalous species Gossypium gossypioides. In one lineage, we detected both gene duplication and sequence polymorphisms that transcend species boundaries; sequences in this lineage formed a monophyletic clade, yet no taxon within the clade contained a monophyletic collection of sequences. Potential explanations for this latter phenomenon, including gene duplication, gene flow, and lineage sorting, are discussed.

Decoupled evolution of coding region and mRNA expression patterns after gene duplication: implications for the neutralist-selectionist debate

The neutralist perspective on molecular evolution maintains that the vast majority of mutations affecting gene function are neutral or deleterious. After a gene duplication where both genes are retained, it predicts that original and duplicate genes diverge at clock-like rates. This prediction is usually tested for coding sequences, but can also be applied to another important aspect of gene function, the genes' expression pattern. Moreover, if both sequence and expression pattern diverge at clock-like rates, a correlation between divergence in sequence and divergence in expression patterns is expected. Duplicate gene pairs with more highly diverged sequences should also show more highly diverged expression patterns. This prediction is tested for a large sample of duplicated genes in the yeast Saccharomyces cerevisiae, using both genome sequence and microarray expression data. Only a weak correlation is observed, suggesting that coding sequence and mRNA expression patterns of duplicate gene pairs evolve independently and at vastly different rates. Implications of this finding for the neutralist-selectionist debate are discussed.

The alpha-mannosidases: phylogeny and adaptive diversification

alpha-Mannosidase enzymes comprise a class of gylcoside hydrolases involved in the maturation and degradaton of glycoprotein-linked oligosaccharides. Various alpha-mannosidase enzymatic activities are encoded by an ancient and ubiquitous gene superfamily. A comparative sequence analysis was employed here to characterize the evolutionary relationships and dynamics of the alpha-mannosidase superfamily. A series of lineage-specific BLAST searches recovered the first ever recognized archaean and eubacterial alpha-mannosidase sequences, in addition to numerous eukaryotic sequences. Motif-based alignment and subsequent phylogenetic analysis of the entire superfamily revealed the presence of three well-supported monophyletic clades that represent discrete alpha-mannosidase families. The comparative method was used to evaluate the phylogenetic distribution of alpha-mannosidase functional variants within families. Results of this analysis demonstrate a pattern of functional diversification of alpha-mannosidase paralogs followed by conservation of function among orthologs. Nucleotide polymorphism among the most closely related pair of duplicated genes was analyzed to evaluate the role of natural selection in the functional diversification of alpha-mannosidase paralogs. Ratios of nonsynonymous and synonymous variation show an increase in the rate of nonsynonymous change after duplication and a relative excess of fixed nonsynonymous changes between the two groups of paralogs. These data point to a possible role for positive Darwinian selection in the evolution of alpha-mannosidase functional diversification following gene duplication.

Genome evolution in polyploids

Polyploidy is a prominent process in plants and has been significant in the evolutionary history of vertebrates and other eukaryotes. In plants, interdisciplinary approaches combining phylogenetic and molecular genetic perspectives have enhanced our awareness of the myriad genetic interactions made possible by polyploidy. Here, processes and mechanisms of gene and genome evolution in polyploids are reviewed. Genes duplicated by polyploidy may retain their original or similar function, undergo diversification in protein function or regulation, or one copy may become silenced through mutational or epigenetic means. Duplicated genes also may interact through inter-locus recombination, gene conversion, or concerted evolution. Recent experiments have illuminated important processes in polyploids that operate above the organizational level of duplicated genes. These include inter-genomic chromosomal exchanges, saltational, non-Mendelian genomic evolution in nascent polyploids, inter-genomic invasion, and cytonuclear stabilization. Notwithstanding many recent insights, much remains to be learned about many aspects of polyploid evolution, including: the role of transposable elements in structural and regulatory gene evolution; processes and significance of epigenetic silencing; underlying controls of chromosome pairing; mechanisms and functional significance of rapid genome changes; cytonuclear accommodation; and coordination of regulatory factors contributed by two, sometimes divergent progenitor genomes. Continued application of molecular genetic approaches to questions of polyploid genome evolution holds promise for producing lasting insight into processes by which novel genotypes are generated and ultimately into how polyploidy facilitates evolution and adaptation.

Genome evolution in polyploids

Polyploidy is a prominent process in plants and has been significant in the evolutionary history of vertebrates and other eukaryotes. In plants, interdisciplinary approaches combining phylogenetic and molecular genetic perspectives have enhanced our awareness of the myriad genetic interactions made possible by polyploidy. Here, processes and mechanisms of gene and genome evolution in polyploids are reviewed. Genes duplicated by polyploidy may retain their original or similar function, undergo diversification in protein function or regulation, or one copy may become silenced through mutational or epigenetic means. Duplicated genes also may interact through inter-locus recombination, gene conversion, or concerted evolution. Recent experiments have illuminated important processes in polyploids that operate above the organizational level of duplicated genes. These include inter-genomic chromosomal exchanges, saltational, non-Mendelian genomic evolution in nascent polyploids, inter-genomic invasion, and cytonuclear stabilization. Notwithstanding many recent insights, much remains to be learned about many aspects of polyploid evolution, including: the role of transposable elements in structural and regulatory gene evolution; processes and significance of epigenetic silencing; underlying controls of chromosome pairing; mechanisms and functional significance of rapid genome changes; cytonuclear accommodation; and coordination of regulatory factors contributed by two, sometimes divergent progenitor genomes. Continued application of molecular genetic approaches to questions of polyploid genome evolution holds promise for producing lasting insight into processes by which novel genotypes are generated and ultimately into how polyploidy facilitates evolution and adaptation.

Duplicated genes evolve independently after polyploid formation in cotton

Of the many processes that generate gene duplications, polyploidy is unique in that entire genomes are duplicated. This process has been important in the evolution of many eukaryotic groups, and it occurs with high frequency in plants. Recent evidence suggests that polyploidization may be accompanied by rapid genomic changes, but the evolutionary fate of discrete loci recently doubled by polyploidy (homoeologues) has not been studied. Here we use locus-specific isolation techniques with comparative mapping to characterize the evolution of homoeologous loci in allopolyploid cotton (Gossypium hirsutum) and in species representing its diploid progenitors. We isolated and sequenced 16 loci from both genomes of the allopolyploid, from both progenitor diploid genomes and appropriate outgroups. Phylogenetic analysis of the resulting 73.5 kb of sequence data demonstrated that for all 16 loci (14.7 kb/genome), the topology expected from organismal history was recovered. In contrast to observations involving repetitive DNAs in cotton, there was no evidence of interaction among duplicated genes in the allopolyploid. Polyploidy was not accompanied by an obvious increase in mutations indicative of pseudogene formation. Additionally, differences in rates of divergence among homoeologues in polyploids and orthologues in diploids were indistinguishable across loci, with significant rate deviation restricted to two putative pseudogenes. Our results indicate that most duplicated genes in allopolyploid cotton evolve independently of each other and at the same rate as those of their diploid progenitors. These indications of genic stasis accompanying polyploidization provide a sharp contrast to recent examples of rapid genomic evolution in allopolyploids.

Coding sequence evolution

Dramatic progress has been made in the past ten years in the development of statistical and experimental techniques for investigating features of molecular evolution. Applied to coding regions, these techniques have produced remarkable advances in our understanding of selection for codon usage but, ironically, have had little impact on our understanding of protein evolution. That may be about to change.

Origin of new genes and source for N-terminal domain of the chimerical gene, jingwei, in Drosophila

This paper deals with a general question posed by the origin of new processed chimerical genes: when a new retrosequence inserts into a new genome position, how does it become activated and acquire novel protein function by recruiting new functional domains and regulatory elements? Jingwei (jgw), a newly evolved functional gene with a chimerical structure in Drosophila, provides an opportunity to examine such questions. The source of its exon encoding C-terminal peptide has been identified as an Adh retrosequence, which extends the concept of exon shuffling from recombination to retroposition as a general molecular mechanism for the origin of a new gene. However, the origin of 5' exons remains unclear. We examined two hypotheses concerning the origin of these non-Adh-derived jgw exons: (i) these exons might originate from a unique genomic sequence that fortuitously evolved a standard intron-exon structure and regulatory sequence for jgw; (ii) these exons might be a duplicate of an unrelated previously existing gene. Genomic Southern analysis, in conjunction with construction and screening of a genomic bookshelf (sub-library), was conducted in a group of Drosophila species. The results demonstrated that there are duplicate genes containing the same structure as the recruited portion of jgw. We name this duplicate gene in Drosophila teissieri and Drosophila yakuba and its orthologous gene in Drosophila melanogaster as yellow-emperor (ymp). Thus, the 5' exons/introns originated from a previously existing gene that provided new modules with specific sub-function to create jgw.