More genes in vertebrates?

With the acquisition of complete genome sequences from several animals, there is renewed interest in the pattern of genome evolution on our own lineage. One key question is whether gene number increased during chordate or vertebrate evolution. It is argued here that comparing the total number of genes between a fly, a nematode and human is not appropriate to address this question. Extensive gene loss after duplication is one complication; another is the problem of comparing taxa that are phylogenetically very distant. Amphioxus and tunicates are more appropriate animals for comparison to vertebrates. Comparisons of clustered homeobox genes, where gene loss can be identified, reveals a one to four mode of evolution for Hox and ParaHox genes. Analyses of other gene families in amphioxus and vertebrates confirm that gene duplication was very widespread on the vertebrate lineage. These data confirm that vertebrates have more genes than their closest invertebrate relatives, acquired through gene duplication. abbreviations IHGSC, International Human Genome Sequencing Consortium; TCESC, The C. elegans Sequencing Consortium.     

Evolution of the F-Box Gene Family in Euarchontoglires: Gene Number Variation and Selection Patterns

F-box proteins are substrate adaptors used by the SKP1–CUL1–F-box protein (SCF) complex, a type of E3 ubiquitin ligase complex in the ubiquitin proteasome system (UPS). SCF-mediated ubiquitylation regulates proteolysis of hundreds of cellular proteins involved in key signaling and disease systems. However, our knowledge of the evolution of the F-box gene family in Euarchontoglires is limited. In the present study, 559 F-box genes and nine related pseudogenes were identified in eight genomes. Lineage-specific gene gain and loss events occurred during the evolution of Euarchontoglires, resulting in varying F-box gene numbers ranging from 66 to 81 among the eight species. Both tandem duplication and retrotransposition were found to have contributed to the increase of F-box gene number, whereas mutation in the F-box domain was the main mechanism responsible for reduction in the number of F-box genes, resulting in a balance of expansion and contraction in the F-box gene family. Thus, the Euarchontoglire F-box gene family evolved under a birth-and-death model. Signatures of positive selection were detected in substrate-recognizing domains of multiple F-box proteins, and adaptive changes played a role in evolution of the Euarchontoglire F-box gene family. In addition, single nucleotide polymorphism (SNP) distributions were found to be highly non-random among different regions of F-box genes in 1092 human individuals, with domain regions having a significantly lower number of non-synonymous SNPs.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Out of the testis,into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome‐scale sequence and multi‐tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well‐known “out of the testis” origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary‐expressed genes. Therefore, I propose a novel “into the ovary” hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined “out of the testis, into the ovary” evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex‐specific adaptation in Drosophila.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Exon duplication from a fork head to a homeodomain protein

The evolution of complex organisms such as animals requires a large expansion of the number of genes controlling developmental events. In addition, it is thought that domains are shuffled between genes to further increase the complexity and generate new types of genes and functions. Working with the Caenorhabditis elegans homeobox gene ceh-43, the orthologue of fly Distal-less (Dll), we observed sequence similarity to the C. elegans gene fkh-1. Now, with the complete genomic sequence available, we examined this similarity in detail. The region of similarity is confined essentially to one exon in the carboxy terminus of the two genes. Based on the gene structure, we think that an exon of fkh-1 was duplicated to the carboxy terminus of ceh-43, where it was incorporated as the last exon. This duplication event seems to have happened recently since the similarity on the nucleotide level is higher than the sequence similarity between fkh-1 of C. elegans and C. briggsae. Potentially the duplication event was mediated via a short region of sequence similarity between the two open reading frames of the genes. This duplication event clear shows that a part of a gene can successfully be juxtaposed to another gene. These events may perhaps not be rare. Received: 28 March 1999 / Accepted: 8 June 1999     

Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth

Gene duplication plays important roles in organismal evolution, because duplicate genes provide raw materials for the evolution of mechanisms controlling physiological and/or morphological novelties. Gene duplication can occur via several mechanisms, including segmental duplication, tandem duplication and retroposition. Although segmental and tandem duplications have been found to be important for the expansion of a number of multigene families, the contribution of retroposition is not clear. Here we show that plant SKP1 genes have evolved by multiple duplication events from a single ancestral copy in the most recent common ancestor (MRCA) of eudicots and monocots, resulting in 19 ASK (Arabidopsis SKP1-like) and 28 OSK (Oryza SKP1-like) genes. The estimated birth rates are more than ten times the average rate of gene duplication, and are even higher than that of other rapidly duplicating plant genes, such as type I MADS box genes, R genes, and genes encoding receptor-like kinases. Further analyses suggest that a relatively large proportion of the duplication events may be explained by tandem duplication, but few, if any, are likely to be due to segmental duplication. In addition, by mapping the gain/loss of a specific intron on gene phylogenies, and by searching for the features that characterize retrogenes/retrosequences, we show that retroposition is an important mechanism for expansion of the plant SKP1 gene family. Specifically, we propose that two and three ancient retroposition events occurred in lineages leading to Arabidopsis and rice, respectively, followed by repeated tandem duplications and chromosome rearrangements. Our study represents a thorough investigation showing that retroposition can play an important role in the evolution of a plant gene family whose members do not encode mobile elements.     

Duplication and diversification of trehalase confers evolutionary advantages on lepidopteran insects

Gene duplication provides a major source of new genes for evolutionary novelty and ecological adaptation. However, the maintenance of duplicated genes and their relevance to adaptive evolution has long been debated. Insect trehalase (Treh) plays key roles in energy metabolism, growth, and stress recovery. Here, we show that the duplication of Treh in Lepidoptera (butterflies and moths) is linked with their adaptation to various environmental stresses. Generally, two Treh genes are present in insects: Treh1 and Treh2. We report three distinct forms of Treh in lepidopteran insects, where Treh1 was duplicated into two gene clusters (Treh1a and Treh1b). These gene clusters differ in gene expression patterns, enzymatic properties, and subcellular localizations, suggesting that the enzymes probably underwent sub‐ and/or neofunctionalization in the lepidopteran insects. Interestingly, selective pressure analysis provided significant evidence of positive selection on duplicate Treh1b gene in lepidopteran insect lineages. Most positively selected sites were located in the alpha‐helical region, and several sites were close to the trehalose binding and catalytic sites. Subcellular adaptation of duplicate Treh1b driven by positive selection appears to have occurred as a result of selected changes in specific sequences, allowing for rapid reprogramming of duplicated Treh during evolution. Our results suggest that gene duplication of Treh and subsequent functional diversification could increase the survival rate of lepidopteran insects through various regulations of intracellular trehalose levels, facilitating their adaptation to diverse habitats. This study provides evidence regarding the mechanism by which gene family expansion can contribute to species adaptation through gene duplication and subsequent functional diversification.     

Evolutionary relationship between human and mouse immunoglobulin kappa light chain variable region genes

Similar to the Igh-V multigene family, the human or mouse Igk-V repertoirer is a distorted continuum of homologous genes that may be grouped into families displaying >80% nucleic acid sequence similarity among their members. systematic interspecies sequence comparisons reveal that most human Igk-V gene families exhibit clear homology to mouse Ogk-V families (sequence similarity >74%). A hypothetical phylogenetic tree of Igk-V genes predicts that a minimum of seven Igk-V genes/families predate mammalian radiation. In two cases, several interrelated mouse Igk-V families exhibit phylogenetic equidistance with just one human Igk-V family, implying a more pronounced divergence for the elevated number of Igk-V gene families in the mouse. Mouse-human Igk-V comaprisons, moreover, illustrate how expansion, contraction, and perhaps deletion of Igk-V gene families shape the Igk-V repertoire during mammalian evolution.     

Massive gene rearrangements of mitochondrial genomes and implications for the phylogeny of Trichoptera (Insecta)

Mitochondrial genomes have been widely used for phylogenetic reconstruction and evolutionary analysis in various groups of Insecta. Gene rearrangements in the mitogenome can be informative characters for phylogenetic reconstruction and adaptive evolution. Trichoptera is one of the most important groups of aquatic insects. Prior to this study, complete mitogenomes from Trichoptera were restricted to eight families, resulting in a biased view of their mitogenome structure and evolution. Here, we assemble new mitogenomes for 66 species by high-throughput sequencing. The mitogenomes of 19 families and 47 genera are documented for the first time. Combined with 16 previously published mitogenomes of Trichoptera, we find 14 kinds of gene rearrangement patterns novel for Trichoptera, including rearrangement of protein-coding genes, tRNAs and control regions. Simultaneously, we provide evidence for the occurrence of tandem duplication and non-random loss events in the mitogenomes of three families. Phylogenetic analyses show that Hydroptilidae was recovered as a sister group to Annulipalpia. The increased nucleotide substitution rate and adaptive evolution may have affected the mitochondrial gene rearrangements in Trichoptera. Our study offers new insights into the mechanisms and patterns of mitogenome rearrangements in Insecta at large and into the usefulness of mitogenomic gene order as a phylogenetic marker within Trichoptera.     

Evolutionary stability of the immunoglobulin heavy chain variable region gene families in teleost

 A comparison between related species would allow us to study the evolutionary changes in complex gene families. To investigate the evolution of immunoglobulin VH gene families in lower vertebrates, we compared cDNA VH clones from two related teleost fish species, Arctic charr (Salvelinus alpinus) and rainbow trout (Oncorhynchus mykiss), which are separated from their common ancestor by 12–20 million years (MY). The results showed that randomly isolated charr VH genes could be closely grouped to known VH genes of rainbow trout, suggesting that the VH family structure is stable during 12–20 MY and that the total number of VH families changes only gradually over a longer period. This finding also led us to define eight VH gene families of Arctic charr, designated Salalp VH I, VH II, and so on. The presence of species-specific amino acids suggests that non-reciprocal genetic exchanges (e.g., gene duplication) play an important role in shaping the evolution of the V gene family. Received: 23 July 1997 / Revised: 6 October 1997     

Genome-wide identification and comparative analysis of lipocalin families in Lepidoptera with an emphasis on Bombyx mori

Lipocalins exhibit functional diversity, including roles in retinol transport, invertebrate cryptic coloration, and stress response. However, genome-wide identification and characterization of lipocalin in the insect lineage have not been thoroughly explored. Here, we found that a lineage-specific expansion of the lipocalin genes in Lepidoptera occurred in large part due to tandem duplication events and several lipocalin genes involving insect coloration were expanded more via tandem duplication in butterflies. A comparative analysis of conserved motifs showed both conservation and divergence of lepidopteran lipocalin family protein structures during evolution. We observe dynamic changes in tissue expression preference of paralogs in Bombyx mori, suggesting differential contribution of paralogs to specific organ functions during evolution. Subcellular localization experiments revealed that lipocalins localize to the cytoplasm, nuclear membrane, or nucleus in BmN cells. Moreover, several lipocalin genes exhibited divergent responses to abiotic and biotic stresses, and 1 lipocalin gene was upregulated by 300 fold in B. mori. These results suggest that lipocalins act as signaling components in defense responses by mediating crosstalk between abiotic and biotic stress responses. This study deepens our understanding of the comprehensive characteristics of lipocalins in insects.     

Evolutionary genomics analysis reveals gene expansion and functional diversity of arylalkylamine N-acetyltransferases in the Culicinae subfamily of mosquitoes

Arylalkylamine N-acetyltransferase (aaNAT), considered a potential new insecticide target, catalyzes the acetylation of arylalkylamine substrates such as serotonin and dopamine and, hence, mediates diverse functions in insects. However, the origin of insect aaNATs (iaaNATs) and the evolutionary process that generates multiple aaNATs in mosquitoes remain largely unknown. Here, we have analyzed the genomes of 33 species to explore and expand our understanding of the molecular evolution of this gene family in detail. We show that aaNAT orthologs are present in Bacteria, Cephalochordata, Chondrichthyes, Cnidaria, Crustacea, Mammalia, Placozoa, and Teleoste, as well as those from a number of insects, but are absent in some species of Annelida, Echinozoa, and Mollusca as well as Arachnida. Particularly, more than 10 aaNATs were detected in the Culicinae subfamily of mosquitoes. Molecular evolutionary analysis of aaNAT/aaNAT-like genes in mosquitoes reveals that tandem duplication events led to gene expansion in the Culicinae subfamily of mosquitoes more than 190 million years ago. Further selection analysis demonstrates that mosquito aaNATs evolved under strongly positive pressures that generated functional diversity following gene duplication events. Overall, this study may provide novel insights into the molecular evolution of the aaNAT family in mosquitoes.     

Early vertebrate chromosome duplications and the evolution of the neuropeptide Y receptor gene regions

Background  

One of the many gene families that expanded in early vertebrate evolution is the neuropeptide (NPY) receptor family of G-protein coupled receptors. Earlier work by our lab suggested that several of the NPY receptor genes found in extant vertebrates resulted from two genome duplications before the origin of jawed vertebrates (gnathostomes) and one additional genome duplication in the actinopterygian lineage, based on their location on chromosomes sharing several gene families. In this study we have investigated, in five vertebrate genomes, 45 gene families with members close to the NPY receptor genes in the compact genomes of the teleost fishes Tetraodon nigroviridis and Takifugu rubripes. These correspond to Homo sapiens chromosomes 4, 5, 8 and 10.     

A chromosome-level Amaranthus cruentus genome assembly highlights gene family evolution and biosynthetic gene clusters that may underpin the nutritional value of this traditional crop

Traditional crops have historically provided accessible and affordable nutrition to millions of rural dwellers but have been neglected, with most modern agricultural systems over-reliant on a small number of internationally traded crops. Traditional crops are typically well-adapted to local agro-ecological conditions and many are nutrient-dense. They can play a vital role in local food systems through enhanced nutrition (particularly where diets are dominated by starch crops), food security and livelihoods for smallholder farmers, and a climate-resilient and biodiverse agriculture. Using short-read, long-read and phased sequencing technologies, we generated a high-quality chromosome-level genome assembly for Amaranthus cruentus, an under-researched crop with micronutrient- and protein-rich leaves and gluten-free seed, but lacking improved varieties, with respect to productivity and quality traits. The 370.9 Mb genome demonstrates a shared whole genome duplication with a related species, Amaranthus hypochondriacus. Comparative genome analysis indicates chromosomal loss and fusion events following genome duplication that are common to both species, as well as fission of chromosome 2 in A. cruentus alone, giving rise to a haploid chromosome number of 17 (versus 16 in A. hypochondriacus). Genomic features potentially underlying the nutritional value of this crop include two A. cruentus-specific genes with a likely role in phytic acid synthesis (an anti-nutrient), expansion of ion transporter gene families, and identification of biosynthetic gene clusters conserved within the amaranth lineage. The A. cruentus genome assembly will underpin much-needed research and global breeding efforts to develop improved varieties for economically viable cultivation and realization of the benefits to global nutrition security and agrobiodiversity.     

Selection and Methionine Accumulation in the Fat Body Protein 2 Gene (FBP2), a Duplicate of the Drosophila Alcohol Dehydrogenase (ADH) Gene

The Drosophila fat body protein 2 gene (Fbp2) is an ancient duplication of the alcohol dehydrogenase gene (Adh) which encodes a protein that differs substantially from ADH in its methionine content. In D. melanogaster, there is one methionine in ADH, while there are 51 (20% of all amino acids) in FBP2. Methionine is involved in 46% of amino acid replacements when Fbp2 DNA sequences are compared between D. melanogaster and D. pseudoobscura. Methionine accumulation does not affect conserved residues of the ADH-ADH^r-FBP2 multigene family. The multigene family has evolved by replacement of mildly hydrophobic amino acids by methionine with no apparent reversion. Its short-term evolution was compared between two Drosophila species, while its long-term evolution was compared between two genera belonging respectively to acalyptrate and calyptrate Diptera, Drosophila and Sarcophaga. The pattern of nucleotide substitution was consistent with an independent accumulation of methionines at the Fbp2 locus in each lineage. Under a steady-state model, the rate of methionine accumulation was constant in the lineage leading to Drosophila, and was twice as fast as that in the calyptrate lineage. Substitution rates were consistent with a slight positive selective advantage for each methionine change in about one-half of amino acid sites in Drosophila. This shows that selection can potentially account for a large proportion of amino acid replacements in the molecular evolution of proteins. Received: 12 December 1994 / Accepted: 15 April 1996     

Duplications and deletions of Vh genes in inbred strains of mice

evolution of variable region (Vh) gene family copy number and polymorphism was investigated by the analysis of the immunoglobulin heavy chain variable region (Igh-V) locus in 74 inbred strains and substrains of mice. Several strains were found to have slight differences from Igh-V haplotypes previously identified, usually of a single Vh gene family. These results indicate that the evolution of copy number in the mouse Igh-V locus proceeds largely by the accumulation of incremental changes, reflecting the clustered organization of the mouse Igh-V locus. We have found no evidence of very large or frequent duplication or deletion events indicative of rapid expansion or contraction processes. The existence of one or more particularly large Vh gene families most likely reflects random copy number variation, rather than selection for the amplification of their members. The identification of strains with recombinant Vh gene arrays demonstrates that recombination, both within and between haplotypes, appears to be the predominant mechanism generating the high restriction fragment length polymorphism in the Igh-V locus.Abbreviations used in this paper Igh-V immunoglobulin heavy chain variable region locus - Vh heavy chain variable region gene - Dh heavy chain diversity region gene - V immunoglobulin kappa light chain variable region gene - V T-cell receptor beta chain variable region gene