Gene Family Size Research Articles

Amino acid transporter expansions associated with the evolution of obligate endosymbiosis in sap-feeding insects (Hemiptera: sternorrhyncha).

BackgroundMutualistic obligate endosymbioses shape the evolution of endosymbiont genomes, but their impact on host genomes remains unclear. Insects of the sub-order Sternorrhyncha (Hemiptera) depend on bacterial endosymbionts for essential amino acids present at low abundances in their phloem-based diet. This obligate dependency has been proposed to explain why multiple amino acid transporter genes are maintained in the genomes of the insect hosts. We implemented phylogenetic comparative methods to test whether amino acid transporters have proliferated in sternorrhynchan genomes at rates grater than expected by chance.ResultsBy applying a series of methods to reconcile gene and species trees, inferring the size of gene families in ancestral lineages, and simulating the null process of birth and death in multi-gene families, we uncovered a 10-fold increase in duplication rate in the AAAP family of amino acid transporters within Sternorrhyncha. This gene family expansion was unmatched in other closely related clades lacking endosymbionts that provide essential amino acids.ConclusionsOur findings support the influence of obligate endosymbioses on host genome evolution by both inferring significant expansions of gene families involved in symbiotic interactions, and discovering increases in the rate of duplication associated with multiple emergences of obligate symbiosis in Sternorrhyncha.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0315-3) contains supplementary material, which is available to authorized users.

Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families.

The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.

Evolutionary Perspectives on Human Fungal Pathogens.

Owing to their small size and paucity of phenotypic characters, progress in the evolutionary biology of microbes in general, and human pathogenic fungi in particular, has been linked to a series of advances in DNA sequencing over the past quarter century. Phylogenetics was the first area to benefit, with the achievement of a basic understanding of fungal phylogeny. Population genetics was the next advance, finding cryptic species everywhere, and recombination in species previously thought to be asexual. Comparative genomics saw the next advance, in which variation in gene content and changes in gene family size were found to be important sources of variation. Fungal population genomics is showing that gene flow among closely related populations and species provides yet another source of adaptive, genetic variation. Now, two means to associate genetic variation with phenotypic variation, "reverse ecology" for adaptive phenotypes, and genome-wide association of any phenotype, are letting evolutionary biology make a profound contribution to molecular developmental biology of pathogenic fungi.

Origin of a novel regulatory module by duplication and degeneration of an ancient plant transcription factor

It is commonly believed that gene duplications provide the raw material for morphological evolution. Both the number of genes and size of gene families have increased during the diversification of land plants. Several small proteins that regulate transcription factors have recently been identified in plants, including the LITTLE ZIPPER (ZPR) proteins. ZPRs are post-translational negative regulators, via heterodimerization, of class III Homeodomain Leucine Zipper (C3HDZ) proteins that play a key role in directing plant form and growth. We show that ZPR genes originated as a duplication of a C3HDZ transcription factor paralog in the common ancestor of euphyllophytes (ferns and seed plants). The ZPRs evolved by degenerative mutations resulting in loss all of the C3HDZ functional domains, except the leucine zipper that modulates dimerization. ZPRs represent a novel regulatory module of the C3HDZ network unique to the euphyllophyte lineage, and their origin correlates to a period of rapid morphological changes and increased complexity in land plants. The origin of the ZPRs illustrates the significance of gene duplications in creating developmental complexity during land plant evolution that likely led to morphological evolution.

Gene duplication followed by exon structure divergence substitutes for alternative splicing in zebrafish.

In this study we report novel findings regarding the evolutionary relationship between gene duplication and alternative splicing, two processes that increase proteomic diversity. By studying teleost fish, we find that gene duplication followed by exon structure divergence between paralogs, but not gene duplication alone, leads to a significant reduction in alternative splicing, as measured by both the proportion of genes that undergo alternative splicing as well as mean number of transcripts per gene. Additionally, we show that this effect is independent of gene family size and gene function. Furthermore, we provide evidence that the reduction in alternative splicing may be due to the partitioning of ancestral splice forms among the duplicate genes - a form of subfunctionalization. Taken together these results indicate that exon structure evolution subsequent to gene duplication may be a common substitute for alternative splicing.

Strong links between genomic and anatomical diversity in both mammalian olfactory chemosensory systems

Mammalian olfaction comprises two chemosensory systems: the odorant-detecting main olfactory system (MOS) and the pheromone-detecting vomeronasal system (VNS). Mammals are diverse in their anatomical and genomic emphases on olfactory chemosensation, including the loss or reduction of these systems in some orders. Despite qualitative evidence linking the genomic evolution of the olfactory systems to specific functions and phenotypes, little work has quantitatively tested whether the genomic aspects of the mammalian olfactory chemosensory systems are correlated to anatomical diversity. We show that the genomic and anatomical variation in these systems is tightly linked in both the VNS and the MOS, though the signature of selection is different in each system. Specifically, the MOS appears to vary based on absolute organ and gene family size while the VNS appears to vary according to the relative proportion of functional genes and relative anatomical size and complexity. Furthermore, there is little evidence that these two systems are evolving in a linked fashion. The relationships between genomic and anatomical diversity strongly support a role for natural selection in shaping both the anatomical and genomic evolution of the olfactory chemosensory systems in mammals.

Messages from Powdery Mildew DNA: How the Interplay with a Host Moulds Pathogen Genomes

The genomes of the barley, Arabidopsis and pea powdery mildew are significantly larger than those of related fungi. This is due to an extraordinary expansion of retro-trasposons that are evident as repetitive elements in the sequence. The protein coding genes are fewer than expected due to an overall reduction in the size of gene families, a reduction in the number of paralogs and because of the loss of certain metabolic pathways. Many of these changes have also been observed in the genomes of other taxonomically unrelated obligate biotrophic pathogens. The only group of genes that bucks the trend of gene loss, are those encoding small secreted proteins that bear the hall marks of effectors.

Increased brain size in mammals is associated with size variations in gene families with cell signalling, chemotaxis and immune-related functions

Genomic determinants underlying increased encephalization across mammalian lineages are unknown. Whole genome comparisons have revealed large and frequent changes in the size of gene families, and it has been proposed that these variations could play a major role in shaping morphological and physiological differences among species. Using a genome-wide comparative approach, we examined changes in gene family size (GFS) and degree of encephalization in 39 fully sequenced mammalian species and found a significant over-representation of GFS variations in line with increased encephalization in mammals. We found that this relationship is not accounted for by known correlates of brain size such as maximum lifespan or body size and is not explained by phylogenetic relatedness. Genes involved in chemotaxis, immune regulation and cell signalling-related functions are significantly over-represented among those gene families most highly correlated with encephalization. Genes within these families are prominently expressed in the human brain, particularly the cortex, and organized in co-expression modules that display distinct temporal patterns of expression in the developing cortex. Our results suggest that changes in GFS associated with encephalization represent an evolutionary response to the specific functional requirements underlying increased brain size in mammals.

Abundance and Distribution of Transposable Elements in Two Drosophila QTL Mapping Resources

Here we present computational machinery to efficiently and accurately identify transposable element (TE) insertions in 146 next-generation sequenced inbred strains of Drosophila melanogaster. The panel of lines we use in our study is composed of strains from a pair of genetic mapping resources: the Drosophila Genetic Reference Panel (DGRP) and the Drosophila Synthetic Population Resource (DSPR). We identified 23,087 TE insertions in these lines, of which 83.3% are found in only one line. There are marked differences in the distribution of elements over the genome, with TEs found at higher densities on the X chromosome, and in regions of low recombination. We also identified many more TEs per base pair of intronic sequence and fewer TEs per base pair of exonic sequence than expected if TEs are located at random locations in the euchromatic genome. There was substantial variation in TE load across genes. For example, the paralogs derailed and derailed-2 show a significant difference in the number of TE insertions, potentially reflecting differences in the selection acting on these loci. When considering TE families, we find a very weak effect of gene family size on TE insertions per gene, indicating that as gene family size increases the number of TE insertions in a given gene within that family also increases. TEs are known to be associated with certain phenotypes, and our data will allow investigators using the DGRP and DSPR to assess the functional role of TE insertions in complex trait variation more generally. Notably, because most TEs are very rare and often private to a single line, causative TEs resulting in phenotypic differences among individuals may typically fail to replicate across mapping panels since individual elements are unlikely to segregate in both panels. Our data suggest that “burden tests” that test for the effect of TEs as a class may be more fruitful.

The family of Dof transcription factors in pine

The Dof proteins are members of a major family of plant transcription factors that contain a single C2/C2 class of zinc-finger DNA-binding domain. The size of gene families and the role of transcription factors in trees have only been partially addressed. We studied the size and structure of the Dof gene family in loblolly and maritime pines. The presence of the conserved Dof domain was a typical feature in the ten Dof genes identified in both pine species. Three of the genes lack intron structures and the other seven have a single intron. An analysis of putative cis elements in the promoter region of the genes suggests that the members of the family play different roles throughout the life cycle of the trees in a way that is unrelated to the subfamily to which they belong. Using the information available in the maritime pine database, we have cloned the corresponding open reading frames and initiated structural and functional studies. The gene expression of maritime pine Dof factors was examined at different developmental stages: from embryo to adult tree. This study has revealed a differential pattern of gene expression based on developmental stage and tissue specificity. The interaction both in vitro and in vivo between Dof factors has been addressed and our data showed a negative effect of Dof 2 in the Dof 5-regulated expression of a target gene in maritime pine. These findings suggest that both genes interact in the regulation of gene expression.

Smallest bitter taste receptor (T2Rs) gene repertoire in carnivores.

Bitter taste reception is presumably associated with dietary selection, preventing animals from ingesting potentially harmful compounds. Accordingly, carnivores, who encounter these toxic substances less often, should have fewer genes associated with bitter taste reception compared with herbivores and omnivores. To investigate the genetic basis of bitter taste reception, we confirmed bitter taste receptor (T2R) genes previously found in the genome sequences of two herbivores (cow and horse), two omnivores (mouse and rat) and one carnivore (dog). We also identified, for the first time, the T2R repertoire from the genome of other four carnivore species (ferret, giant panda, polar bear and cat) and detected 17-20 bitter receptor genes from the five carnivore genomes, including 12-16 intact genes, 0-1 partial but putatively functional genes, and 3-8 pseudogenes. Both the intact T2R genes and the total T2R gene number among carnivores were the smallest among the tested species, supporting earlier speculations that carnivores have fewer T2R genes, herbivores an intermediate number, and omnivores the largest T2R gene repertoire. To further explain the genetic basis for this disparity, we constructed a phylogenetic tree, which showed most of the T2R genes from the five carnivores were one-to-one orthologs across the tree, suggesting that carnivore T2Rs were conserved among mammals. Similarly, the small carnivore T2R family size was likely due to rare duplication events. Collectively, these results strengthen arguments for the connection between T2R gene family size, diet and habit.

The pattern of Phosphate transporter 1 genes evolutionary divergence in Glycine maxL.

BackgroundThe Phosphate transporter 1 (PHT1) gene family has crucial roles in phosphate uptake, translocation, remobilization, and optimization of metabolic processes using of Pi. Gene duplications expand the size of gene families, and subfunctionalization of paralog gene pairs is a predominant tendency after gene duplications. To date, experimental evidence for the evolutionary relationships among different paralog gene pairs of a given gene family in soybean is limited.ResultsAll potential Phosphate transporter 1 genes in Glycine max L. (GmPHT1) were systematically analyzed using both bioinformatics and experimentation. The soybean PHT1 genes originated from four distinct ancestors prior to the Gamma WGT and formed 7 paralog gene pairs and a singleton gene. Six of the paralog gene pairs underwent subfunctionalization, and while GmPHT1;4 paralog gene experienced pseudogenization. Examination of long-term evolutionary changes, six GmPHT1 paralog gene pairs diverged at multiple levels, in aspects of spatio-temporal expression patterns and/or quanta, phosphates affinity properties, subcellular localization, and responses to phosphorus stress.ConclusionsThese characterized divergences occurred in tissue- and/or development-specific modes, or conditional modes. Moreover, they have synergistically shaped the evolutionary rate of GmPHT1 family, as well as maintained phosphorus homeostasis at cells and in the whole plant.

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A.thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the 'core' proteome of green plants; (iii) 70% of A.thaliana genes could be traced back to 450million years ago; and (iv) intriguingly, A.thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome-wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A.thaliana.

Analysing complex Triticeae genomes – concepts and strategies

The genomic sequences of many important Triticeae crop species are hard to assemble and analyse due to their large genome sizes, (in part) polyploid genomes and high repeat content. Recently, the draft genomes of barley and bread wheat were reported thanks to cost-efficient and fast NGS technologies. The genome of barley is estimated to be 5 Gb in size whereas the genome of bread wheat accounts for 17 Gb and harbours an allo-hexaploid genome. Direct assembly of the sequence reads and access to the gene content is hampered by the repeat content. As a consequence, novel strategies and data analysis concepts had to be developed to provide much-needed whole genome sequence surveys and access to the gene repertoires. Here we describe some analytical strategies that now enable structuring of massive NGS data generated and pave the way towards structured and ordered sequence data and gene order. Specifically we report on the GenomeZipper, a synteny driven approach to order and structure NGS survey sequences of grass genomes that lack a physical map. In addition, to access and analyse the gene repertoire of allo-hexaploid bread wheat from the raw sequence reads, a reference-guided approach was developed utilizing representative genes from rice, Brachypodium distachyon, sorghum and barley. Stringent sub-assembly on the reference genes prevented collapsing of homeologous wheat genes and allowed to estimate gene retention rate and determine gene family sizes. Genomic sequences from the wheat sub-genome progenitors enabled to discriminate a large number of sub-assemblies between the wheat A, B or D sub-genome using machine learning algorithms. Many of the concepts outlined here can readily be applied to other complex plant and non-plant genomes.

The Chordate Proteome History Database

The chordate proteome history database (http://ioda.univ-provence.fr) comprises some 20,000 evolutionary analyses of proteins from chordate species. Our main objective was to characterize and study the evolutionary histories of the chordate proteome, and in particular to detect genomic events and automatic functional searches. Firstly, phylogenetic analyses based on high quality multiple sequence alignments and a robust phylogenetic pipeline were performed for the whole protein and for each individual domain. Novel approaches were developed to identify orthologs/paralogs, and predict gene duplication/gain/loss events and the occurrence of new protein architectures (domain gains, losses and shuffling). These important genetic events were localized on the phylogenetic trees and on the genomic sequence. Secondly, the phylogenetic trees were enhanced by the creation of phylogroups, whereby groups of orthologous sequences created using OrthoMCL were corrected based on the phylogenetic trees; gene family size and gene gain/loss in a given lineage could be deduced from the phylogroups. For each ortholog group obtained from the phylogenetic or the phylogroup analysis, functional information and expression data can be retrieved. Database searches can be performed easily using biological objects: protein identifier, keyword or domain, but can also be based on events, eg, domain exchange events can be retrieved. To our knowledge, this is the first database that links group clustering, phylogeny and automatic functional searches along with the detection of important events occurring during genome evolution, such as the appearance of a new domain architecture.

Family Size and Turnover Rates among Several Classes of Small Non–Protein-Coding RNA Genes in Caenorhabditis Nematodes

It is important to understand the forces that shape the size and evolutionary histories of gene families. Here, we investigated the evolution of non–protein-coding RNA genes in the genomes of Caenorhabditis nematodes. We specifically focused on nested arrangements, that is, cases in which an RNA gene is entirely contained in an intron of another gene. Comparing these arrangements between species simplifies the inference of orthology and, therefore, of evolutionary fates of nested genes. Two distinct patterns are evident in the data. Genes encoding small nuclear RNAs (snRNAs) and transfer RNAs form large families, which have persisted since before the common ancestor of Metazoa. Yet, individual genes die relatively rapidly, with few orthologs having survived since the divergence of Caenorhabditis elegans and Caenorhabditis briggsae. In contrast, genes encoding small nucleolar RNAs (snoRNAs) are either single-copy or form small families. Individual snoRNAs turn over at a relatively slow rate—most C. elegans genes have clearly identifiable orthologs in C. briggsae. We also found that in Drosophila, genes from larger snRNA families die at a faster rate than their counterparts from single-gene families. These results suggest that a relationship between family size and the rate of gene turnover may be a general feature of genome evolution.

BadiRate: estimating family turnover rates by likelihood-based methods

The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset. BadiRate software and documentation is available from http://www.ub.edu/softevol/badirate.

A Bayesian model for gene family evolution

BackgroundA birth and death process is frequently used for modeling the size of a gene family that may vary along the branches of a phylogenetic tree. Under the birth and death model, maximum likelihood methods have been developed to estimate the birth and death rate and the sizes of ancient gene families (numbers of gene copies at the internodes of the phylogenetic tree). This paper aims to provide a Bayesian approach for estimating parameters in the birth and death model.ResultsWe develop a Bayesian approach for estimating the birth and death rate and other parameters in the birth and death model. In addition, a Bayesian hypothesis test is developed to identify the gene families that are unlikely under the birth and death process. Simulation results suggest that the Bayesian estimate is more accurate than the maximum likelihood estimate of the birth and death rate. The Bayesian approach was applied to a real dataset of 3517 gene families across genomes of five yeast species. The results indicate that the Bayesian model assuming a constant birth and death rate among branches of the phylogenetic tree cannot adequately explain the observed pattern of the sizes of gene families across species. The yeast dataset was thus analyzed with a Bayesian heterogeneous rate model that allows the birth and death rate to vary among the branches of the tree. The unlikely gene families identified by the Bayesian heterogeneous rate model are different from those given by the maximum likelihood method.ConclusionsCompared to the maximum likelihood method, the Bayesian approach can produce more accurate estimates of the parameters in the birth and death model. In addition, the Bayesian hypothesis test is able to identify unlikely gene families based on Bayesian posterior p-values. As a powerful statistical technique, the Bayesian approach can effectively extract information from gene family data and thereby provide useful information regarding the evolutionary process of gene families across genomes.

Origin and evolution of gene families in Bacteria and Archaea

BackgroundComparison of complete genomes of Bacteria and Archaea shows that gene content varies considerably and that genomes evolve quite rapidly via gene duplication and deletion and horizontal gene transfer. We analyze a diverse set of 92 Bacteria and 79 Archaea in order to investigate the processes governing the origin and evolution of families of related genes within genomes.ResultsGenes were clustered into related groups using similarity criteria derived from BLAST. Most clusters contained genes from only one or a small number of genomes, and relatively few core clusters were found that spanned all genomes. Gene clusters found in larger numbers of genomes tended to have larger numbers of genes per genome; however, clusters with unusually large numbers of genes per genome were found among both narrowly and widely distributed clusters. Larger genomes were found to have larger mean gene family sizes and a greater proportion of families of very large size. We used a model of birth, death, and innovation to predict the distribution of gene family sizes. The key parameter is r, the ratio of duplications to deletions. It was found that the model can give a good fit to the observed distribution only if there are several classes of genes with different values of r. The preferred model in most cases had three classes of genes.ConclusionsThere appears to be a rapid rate of origination of new gene families within individual genomes. Most of these gene families are deleted before they spread to large numbers of genomes, which suggests that they may not be generally beneficial to the organisms. The family size distribution is best described by a large fraction of families that tend to have only one or two genes and a small fraction of families of multi-copy genes that are highly prone to duplication. Larger families occur more frequently in larger genomes, indicating higher r in these genomes, possibly due to a greater tolerance for non-beneficial gene duplicates. The smallest genomes contain very few multi-copy families, suggesting a high rate of deletion of all but the most beneficial genes in these genomes.

A New Genomic Evolutionary Model for Rearrangements, Duplications, and Losses that Applies across Eukaryotes and Prokaryotes

Genomic rearrangements have been studied since the beginnings of modern genetics and models for such rearrangements have been the subject of many papers over the last 10 years. However, none of the extant models can predict the evolution of genomic organization into circular unichromosomal genomes (as in most prokaryotes) and linear multichromosomal genomes (as in most eukaryotes). Very few of these models support gene duplications and losses--yet these events may be more common in evolutionary history than rearrangements and themselves cause apparent rearrangements. We propose a new evolutionary model that integrates gene duplications and losses with genome rearrangements and that leads to genomes with either one (or a very few) circular chromosome or a collection of linear chromosomes. Our model is based on existing rearrangement models and inherits their linear-time algorithms for pairwise distance computation (for rearrangement only). Moreover, our model predictions fit observations about the evolution of gene family sizes and agree with the existing predictions about the growth in the number of chromosomes in eukaryotic genomes.