Mosaic evolution of ruminant stomach lysozyme genes

The genomes of ruminant artiodactyls, such as cow and sheep, have approximately 10 lysozyme genes, 4 of which are expressed in the stomach. Most of the duplications of the lysozyme genes occurred 40-50 million years ago, before the divergence of cow and sheep. Despite this, the coding regions of stomach lysozyme genes within a species (e.g., cow, sheep, or deer) are more similar to each other than to lysozyme genes in other ruminants. This observation suggests that the coding regions of the stomach lysozyme genes have evolved in a concerted fashion. Our previous characterization of 3 cow stomach lysozyme genes suggested that it was only the coding exons that had participated in concerted evolution. To determine whether the introns and flanking regions of ruminant stomach lysozyme genes are evolving in a concerted or a divergent fashion, we have isolated and characterized 2 sheep stomach lysozyme genes. Comparison of the sequences of the sheep and cow stomach lysozyme genes clearly shows that the introns and flanking regions have evolved, like the 3' untranslated region of the mRNAs, in a divergent manner. Thus, if the four coding exons are evolving by concerted evolution, then a mosaic pattern of concerted and divergent evolution is occurring in these genes. The independent concerted evolution of coding exons of the ruminant stomach lysozyme gene may have assisted in the accelerated adaptive evolution of the lysozyme to new function in the early ruminant.     

Evolutionary conservation of the chromosomal configuration and regulation of amylase genes among eight species of the Drosophila melanogaster species subgroup

Nuclear DNA was extracted from each of the eight species comprising the Drosophila melanogaster species subgroup. Southern hybridization of this DNA by using a molecular probe specific for the alpha-amylase coding region showed that the duplicated structure of the amylase locus, first found in D. melanogaster, is conserved among all species of the melanogaster subgroup. Evidence is also presented for the concerted evolution of the duplicated genes within each species. In addition, it is shown that the glucose repression of amylase gene expression, which has been extensively studied in D. melanogaster, is not confined to this species but occurs in all eight members of the species subgroup. Thus, both the duplicated gene structure and the glucose repression of Drosophila amylase gene activity are stable over extended periods of evolutionary time.      

Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Molecular Characterization and Evolution of the Amylase Multigene Family of Drosophila ananassae

Drosophila ananassae is known to produce numerous alpha-amylase variants. We have cloned seven different Amy genes in an African strain homozygous for the AMY1,2,3,4 electrophoretic pattern. These genes are organized as two main clusters: the first one contains three intronless copies on the 2L chromosome arm, two of which are tandemly arranged. The other cluster, on the 3L arm, contains two intron-bearing copies. The amylase variants AMY1 and AMY2 have been assigned to the intronless cluster, and AMY3 and AMY4 to the second one. The divergence of coding sequences between clusters is moderate (6.1% in amino acids), but the flanking regions are very different, which could explain their differential regulation. Within each cluster, coding and noncoding regions are conserved. Two very divergent genes were also cloned, both on chromosome 3L, but very distant from each other and from the other genes. One is the Amyrel homologous (41% divergent), the second one, Amyc1 (21.6% divergent) is unknown outside the D. ananassae subgroup. These two genes have unknown functions. Received: 30 May 2000 / Accepted: 17 July 2000     

Molecular Evolution of the Duplicated Amy Locus in the Drosophila Melanogaster Species Subgroup: Concerted Evolution Only in the Coding Region and an Excess of Nonsynonymous Substitutions in Speciation

From the analysis of restriction maps of the Amy region in eight sibling species belonging to the Drosophila melanogaster species subgroup, we herein show that the patterns of duplication of the Amy gene are almost the same in all species. This indicates that duplication occurred before speciation within this species subgroup. From the nucleotide sequence data, we show a strong within-species similarity between the duplicated loci in the Amy coding region. This is in contrast to a strong similarity in the 5' and 3' flanking regions within each locus (proximal or distal) throughout the species subgroup. This means that concerted evolution occurred only in the Amy coding region and that differentiated evolution between the duplication occurred in the flanking regions. Moreover, when comparing the species, we also found a significant excess of nonsynonymous substitutions. In particular, all the fixed substitutions specific to D. erecta were found to be nonsynonymous. We thus conclude that adaptive protein evolution occurred in the lineage of D. erecta that is a ``specialist' species for host plants and probably also occurs in the process of speciation in general.     

Complex phylogeographic history of central African forest elephants and its implications for taxonomy

Background

The mitochondrial genomes of snakes are characterized by an overall evolutionary rate that appears to be one of the most accelerated among vertebrates. They also possess other unusual features, including short tRNAs and other genes, and a duplicated control region that has been stably maintained since it originated more than 70 million years ago. Here, we provide a detailed analysis of evolutionary dynamics in snake mitochondrial genomes to better understand the basis of these extreme characteristics, and to explore the relationship between mitochondrial genome molecular evolution, genome architecture, and molecular function. We sequenced complete mitochondrial genomes from Slowinski's corn snake (Pantherophis slowinskii) and two cottonmouths (Agkistrodon piscivorus) to complement previously existing mitochondrial genomes, and to provide an improved comparative view of how genome architecture affects molecular evolution at contrasting levels of divergence.

Results

We present a Bayesian genetic approach that suggests that the duplicated control region can function as an additional origin of heavy strand replication. The two control regions also appear to have different intra-specific versus inter-specific evolutionary dynamics that may be associated with complex modes of concerted evolution. We find that different genomic regions have experienced substantial accelerated evolution along early branches in snakes, with different genes having experienced dramatic accelerations along specific branches. Some of these accelerations appear to coincide with, or subsequent to, the shortening of various mitochondrial genes and the duplication of the control region and flanking tRNAs.

Conclusion

Fluctuations in the strength and pattern of selection during snake evolution have had widely varying gene-specific effects on substitution rates, and these rate accelerations may have been functionally related to unusual changes in genomic architecture. The among-lineage and among-gene variation in rate dynamics observed in snakes is the most extreme thus far observed in animal genomes, and provides an important study system for further evaluating the biochemical and physiological basis of evolutionary pressures in vertebrate mitochondria.     

Concerted evolution of the primate immunoglobulin alpha-gene through gene conversion.

We determined four nucleotide sequences of the hominoid immunoglobulin alpha (C alpha) genes (chimpanzee C alpha 2, gorilla C alpha 2, and gibbon C alpha 1 and C alpha 2 genes), which made possible the examination of gene conversions in all hominoid C alpha genes. The following three methods were used to detect gene conversions: 1) phenetic tree construction; 2) detection of a DNA segment with extremely low variability between duplicated C alpha genes; and 3) a site by site search of shared nucleotide changes between duplicated C alpha genes. Results obtained from method 1 indicated a concerted evolution of the duplicated C alpha genes in the human, chimpanzee, gorilla, and gibbon lineages, while results obtained from method 2 suggested gene conversions in the human, gorilla, and gibbon C alpha genes. With method 3 we identified clusters of shared nucleotide changes between duplicated C alpha genes in human, chimpanzee, gorilla, and gibbon lineages, and in their hypothetical ancestors. In the present study converted regions were identified over the entire C alpha gene region excluding a few sites in the coding region which have escaped from gene conversion. This indicates that gene conversion is a general phenomenon in evolution, that can be clearly observed in non-functional regions.     

Evolution of genes for allelic and isotypic forms of immunoglobulin kappa chains and of the genes for T-cell receptor beta chains in rabbits

New insights into the evolution of the families of genes encoding immunoglobulins and T-cell receptors of rabbits (Oryctolagus cuniculus) have come from molecular genetic studies. In contrast to human and mouse, rabbits were shown to have two genes for the constant region of immunoglobulin light chains (C kappa 1 and C kappa 2 isotypes) and complex allelic variants of K1 (allotypes). Although K1 allotype protein sequences differed at up to 41% of the amino acid positions, 3' untranslated, 5', and 3' flanking regions were conserved, and in the coding regions 78-80% of the codons with differences had replacement changes. Proportions of silent changes and changes in noncoding regions were comparable. Thus, in spite of their markedly different protein sequences, the K1b4, b5, and b9 allotypes appeared to be products of allelic genes. Molecular genetic analyses suggested that they may have undergone rapid divergence after an ancestral K2-like gene duplicated. Some rabbits were found to have two similar T-cell receptor C beta genes as do humans and many strains of mice, but others appeared to have three different C beta. In addition, we found allotypic forms of C beta. Some of the C beta allotypic differences occurred at positions where analogous C kappa allotypic differences were found. We also found V beta in mouse and human that were more similar to rabbit V beta than closely linked rabbit genes were to each other. This contrasts with rabbit immunoglobulin VH gene sequences that reflect concerted evolution. The data suggested that T-cell receptor V beta genes duplicated prior to mammalian radiation.     

Concerted evolution of duplicated mitochondrial control regions in three related seabird species

Background  

Many population genetic and phylogenetic analyses of mitochondrial DNA (mtDNA) assume that mitochondrial genomes do not undergo recombination. Recently, concerted evolution of duplicated mitochondrial control regions has been documented in a range of taxa. Although the molecular mechanism that facilitates concerted evolution is unknown, all proposed mechanisms involve mtDNA recombination.     

The effect of gene conversion on the divergence between duplicated genes

Nonindependent evolution of duplicated genes is called concerted evolution. In this article, we study the evolutionary process of duplicated regions that involves concerted evolution. The model incorporates mutation and gene conversion: the former increases d, the divergence between two duplicated regions, while the latter decreases d. It is demonstrated that the process consists of three phases. Phase I is the time until d reaches its equilibrium value, d(0). In phase II d fluctuates around d(0), and d increases again in phase III. Our simulation results demonstrate that the length of concerted evolution (i.e., phase II) is highly variable, while the lengths of the other two phases are relatively constant. It is also demonstrated that the length of phase II approximately follows an exponential distribution with mean tau, which is a function of many parameters including gene conversion rate and the length of gene conversion tract. On the basis of these findings, we obtain the probability distribution of the level of divergence between a pair of duplicated regions as a function of time, mutation rate, and tau. Finally, we discuss potential problems in genomic data analysis of duplicated genes when it is based on the molecular clock but concerted evolution is common.     

Nucleotide variation of the duplicated amylase genes in Drosophila kikkawai

We examined levels and patterns of the nucleotide polymorphism of the Amylase genes with a head-to-head duplication in Drosophila kikkawai. The levels of variation in D. kikkawai were comparable to those in Drosophila melanogaster. Tajima's test, Fu and Li's test, HKA test, and MK test did not show significant departure from neutrality. We found an excess of replacement changes in the within-locus class, representing polymorphism in one of the duplicated genes, compared with the between-locus class, representing polymorphism shared between the duplicated genes. Most replacement changes in the within-locus class were singletons. These results suggest that most replacement changes are deleterious. A contrasting evolutionary pattern, involving concerted evolution in the coding regions but differential evolution in the 5'-flanking regions, was observed. However, unlike the duplicated Amy genes of D. melanogaster, the coding regions of the duplicated genes in D. kikkawai tended to diverge. Using Ohta's model of the small multigene family, we found that recombination (interchromosomal equal crossing-over) rate was one order higher than gene conversion (unequal crossing-over) rate, resulting in a considerable but incomplete homogenization of the duplicated coding regions. Linkage disequilibria were found in the intron as well as within and around the regulatory cis-element sequences of one of the duplicated genes (Amy1). The possible causes of these linkage disequilibria were discussed.     

Identification of a two-domain antifreeze protein gene in Antarctic eelpout <Emphasis Type="Italic">Lycodichthys dearborni</Emphasis>

Type III antifreeze proteins (AFP III) in the Antarctic eelpout Lycodichthys dearborni contain at least two size variants—a 7-kDa protein family and a specific 14-kDa isoform composed of two 7-kDa domains linked in tandem. We report the characterization of a two-domain AFP III gene from L. dearborni, and propose that the two-domain AFP III gene arose from a single-domain AFP III gene through duplication and degeneration. AT-rich regions played an important role in the degeneration of the duplicated AFP III gene that resulted in the concatenation of two originally separated 7-kDa AFP-coding exons into a single gene. We also identified a pseudo-AFP III gene interrupted at an AT-rich coding region, supporting AT-rich regions as hotspots for DNA recombination in AFP III gene evolution. Interestingly, study of AFP III genes in the related Antarctic eelpout Pachycara brachycephalum showed absence of two- and multi-domain AFP III genes, indicating that modes of AFP III gene family evolution are specific within species. Nucleotide sequences have been deposited into NCBI Genbank under Accession Numbers: EU627165, EU627166.     

Complete sequence of the mitochondrial genome of Tetrahymena thermophila and comparative methods for identifying highly divergent genes

The complete sequence of the mitochondrial genome of Tetrahymena thermophila has been determined and compared with the mitochondrial genome of Tetrahymena pyriformis. The sequence similarity clearly indicates homology of the entire T.thermophila and T.pyriformis mitochondrial genomes. The T.thermophila genome is very compact, most of the intergenic regions are short (only three are longer than 63 bp) and comprise only 3.8% of the genome. The nad9 gene is tandemly duplicated in T.thermophila. Long terminal inverted repeats and the nad9 genes are undergoing concerted evolution. There are 55 putative genes: three ribosomal RNA genes, eight transfer RNA genes, 22 proteins with putatively assigned functions and 22 additional open reading frames of unknown function. In order to extend indications of homology beyond amino acid sequence similarity we have examined a number of physico-chemical properties of the mitochondrial proteins, including theoretical pI, molecular weight and particularly the predicted transmembrane spanning regions. This approach has allowed us to identify homologs to ymf58 (nad4L), ymf62 (nad6) and ymf60 (rpl6).     

Mapping and recombination analysis of two moth colour mutations,Black moth and Wild wing spot,in the silkworm Bombyx mori

Many lepidopteran insects exhibit body colour variations, where the high phenotypic diversity observed in the wings and bodies of adults provides opportunities for studying adaptive morphological evolution. In the silkworm Bombyx mori, two genes responsible for moth colour mutation, Bm and Ws, have been mapped to 0.0 and 14.7 cM of the B. mori genetic linkage group 17; however, these genes have not been identified at the molecular level. We performed positional cloning of both genes to elucidate the molecular mechanisms that underlie the moth wing- and body-colour patterns in B. mori. We successfully narrowed down Bm and Ws to ~2-Mb-long and 100-kb-long regions on the same scaffold Bm_scaf33. Gene prediction analysis of this region identified 77 candidate genes in the Bm region, whereas there were no candidate genes in the Ws region. Fluorescence in-situ hybridisation analysis in Bm mutant detected chromosome inversion, which explains why there are no recombination in the corresponding region. The comparative genomic analysis demonstrated that the candidate regions of both genes shared synteny with a region associated with wing- and body-colour variations in other lepidopteran species including Biston betularia and Heliconius butterflies. These results suggest that the genes responsible for wing and body colour in B. mori may be associated with similar genes in other Lepidoptera.     

Both Positive and Negative Selection Pressures Contribute to the Polymorphism Pattern of the Duplicated Human CYP21A2 Gene

The human steroid 21-hydroxylase gene (CYP21A2) participates in cortisol and aldosterone biosynthesis, and resides together with its paralogous (duplicated) pseudogene in a multiallelic copy number variation (CNV), called RCCX CNV. Concerted evolution caused by non-allelic gene conversion has been described in great ape CYP21 genes, and the same conversion activity is responsible for a serious genetic disorder of CYP21A2, congenital adrenal hyperplasia (CAH). In the current study, 33 CYP21A2 haplotype variants encoding 6 protein variants were determined from a European population. CYP21A2 was shown to be one of the most diverse human genes (HHe=0.949), but the diversity of intron 2 was greater still. Contrary to previous findings, the evolution of intron 2 did not follow concerted evolution, although the remaining part of the gene did. Fixed sites (different fixed alleles of sites in human CYP21 paralogues) significantly accumulated in intron 2, indicating that the excess of fixed sites was connected to the lack of effective non-allelic conversion and concerted evolution. Furthermore, positive selection was presumably focused on intron 2, and possibly associated with the previous genetic features. However, the positive selection detected by several neutrality tests was discerned along the whole gene. In addition, the clear signature of negative selection was observed in the coding sequence. The maintenance of the CYP21 enzyme function is critical, and could lead to negative selection, whereas the presumed gene regulation altering steroid hormone levels via intron 2 might help fast adaptation, which broadly characterizes the genes of human CNVs responding to the environment.     

Mosaic Evolution of Silk Genes in Aliatypus Trapdoor Spiders (Mygalomorphae, Antrodiaetidae)

Spider silk genes are composed mostly of repetitive sequence that is flanked by non-repetitive terminal regions. Inferences about the evolutionary processes that influenced silk genes have largely been made from analyses using distantly related taxa and ancient silk gene duplicates. These studies have relied on comparisons across the conserved non-repetitive terminal regions to determine orthologous and paralogous relationships, as well as the influence of selection on silk genes. While the repetitive region heavily influences silk fiber mechanical properties, few molecular evolutionary analyses have been conducted on this region due to difficulty in determining homology. Here, we sample internal repetitive and carboxy terminal regions from all extant species of the trapdoor spider genus, Aliatypus. Aliatypus spiders are highly dispersal limited and rely on their silk lined burrow for protection. We determine positional homology across species for the carboxy terminal regions and relative positional homology for the internal repetitive regions. Gene trees based on each of these regions are in good agreement with the Aliatypus species tree, which indicates we sampled single spidroin orthologs in each species. In addition, we find that purifying selection and concerted evolution have acted to conserve Aliatypus spidroin internal repetitive regions. In contrast, selection testing identifies evidence of sites that evolved under positive selection and amino acid replacements that result in radical physicochemical changes in the carboxy terminal region. These findings indicate that comparison of spidroin orthologs across a comprehensive sample of congenerics reveal molecular evolutionary patterns obscured from studies using higher-level sampling of silk encoding genes.     

Phylogeny and the Evolution of the <Emphasis Type="Italic">Amylase</Emphasis> Multigenes in the <Emphasis Type="Italic">Drosophila montium</Emphasis> Species Subgroup

Abstract To investigate the phylogenetic relationships and molecular evolution of α-amylase (Amy) genes in the Drosophila montium species subgroup, we constructed the phylogenetic tree of the Amy genes from 40 species from the montium subgroup. On our tree the sequences of the auraria, kikkawai, and jambulina complexes formed distinct tight clusters. However, there were a few inconsistencies between the clustering pattern of the sequences and taxonomic classification in the kikkawai and jambulina complexes. Sequences of species from other complexes (bocqueti, bakoue, nikananu, and serrata) often did not cluster with their respective taxonomic groups. This suggests that relationships among the Amy genes may be different from those among species due to their particular evolution. Alternatively, the current taxonomy of the investigated species is unreliable. Two types of divergent paralogous Amy genes, the so-called Amy1- and Amy3-type genes, previously identified in the D. kikkawai complex, were common in the montium subgroup, suggesting that the duplication event from which these genes originate is as ancient as the subgroup or it could even predate its differentiation. Thc Amy1-type genes were closer to the Amy genes of D. melanogaster and D. pseudoobscura than to the Amy3-type genes. In the Amy1-type genes, the loss of the ancestral intron occurred independently in the auraria complex and in several Afrotropical species. The GC content at synonymous third codon positions (GC3s) of the Amy1-type genes was higher than that of the Amy3-type genes. Furthermore, the Amy1-type genes had more biased codon usage than the Amy3-type genes. The correlations between GC3s and GC content in the introns (GCi) differed between these two Amy-type genes. These findings suggest that the evolutionary forces that have affected silent sites of the two Amy-type genes in the montium species subgroup may differ.