Evolution of the vertebrate ABC gene family: analysis of gene birth and death

Vertebrate evolution has been largely driven by the duplication of genes that allow for the acquisition of new functions. The ATP-binding cassette (ABC) proteins constitute a large and functionally diverse family of membrane transporters. The members of this multigene family are found in all cellular organisms, most often engaged in the translocation of a wide variety of substrates across lipid membranes. Because of the diverse function of these genes, their large size, and the large number of orthologs, ABC genes represent an excellent tool to study gene family evolution. We have identified ABC proteins from the sea squirt (Ciona intestinalis), zebrafish (Danio rerio), and chicken (Gallus gallus) and, using phylogenetic analysis, identified those genes with a one-to-one orthologous relationship to human ABC proteins. All ABC protein subfamilies found in Ciona and zebrafish correspond to the human subfamilies, with the exception of a single ABCH subfamily gene found only in zebrafish. Multiple gene duplication and deletion events were identified in different lineages, indicating an ongoing process of gene evolution. As many ABC genes are involved in human genetic diseases, and important drug transport phenotypes, the understanding of ABC gene evolution is important to the development of animal models and functional studies.     

Evolutionary History of Eukaryotic α-Glucosidases from the α-Amylase Family

Although some α-glucosidases from the α-amylase family (glycoside hydrolase family GH13) have been studied extensively, their exact number, organization on the chromosome, and orthology/paralogy relationship were unknown. This was true even for important disease vectors where gut α-glucosidase is known to be receptor for the Bin toxin used to control the population of some mosquito species. In some cases orthologs from related species were studied intensively, while potentially important paralogs were omitted. We have, therefore, used a bioinformatics approach to identify all family GH13 α-glucosidases from the selected species from Metazoa (including three mosquito species: Aedes aegypti, Anopheles gambiae, and Culex quinquefasciatus) as well as from Fungi in an effort to characterize their arrangement on the chromosome and evolutionary relationships among orthologs and among paralogs. We also searched for pseudogenes and genes coding for enzymatically inactive proteins with a possible new function. We have found GH13 α-glucosidases mostly in Arthropoda and Fungi where they form gene families, as a result of multiple lineage-specific gene duplications. In mosquito species we have identified 14 α-glucosidase (Aglu) genes of which only five have been biochemically characterized so far, two are putative pseudogenes and the rest remains uncharacterized. We also revealed quite a complex evolutionary history of the eukaryotic α-glucosidases probably involving multiple losses of genes or horizontal gene transfer from bacteria.     

Origin of the murine implantation serine proteinase subfamily

The S1 serine protease family is one of the largest gene families known. Within this family there are several subfamilies that have been grouped together as a result of sequence comparisons and substrate identification. The grouping of related genes allows for the speculation of function for newly found members by comparison and for novel subfamilies by contrast. Analysis of the evolutionary patterns of genes indicates whether or not orthologs are likely to be identified in other species as well as potentially indicating that hypothesized orthologs are in fact not. Looking at subtle differences between subfamily members can reveal intricacies about function and expression. Previously, we have described genes encoding two novel serine proteinases, ISP1 and ISP2, which are most closely related to tryptases. The ISP1 gene encodes the embryo-derived enzyme strypsin, which is necessary for blastocyst hatching and invasion in vitro. Additionally both ISP1 and ISP2 are co-expressed in the endometrial gland during the time of hatching, suggesting that they may also both participate in zona lysis from within the uterine lumen. Here, we demonstrate that the ISPs are tandemly linked within the tryptase cluster on 17A3.3. We suggest that remarkable similarities within the 5'-untranslated and first intron regions of ISP1 and ISP2 may explain their intimate co-regulation in uterus. We also suggest that ISP genes have evolved through gene duplication and that the ISP1 gene has also begun to adopt an additional new function in the murine preimplantation embryo.     

Atg8: an autophagy-related ubiquitin-like protein family

Autophagy-related (Atg) proteins are eukaryotic factors participating in various stages of the autophagic process. Thus far 34 Atgs have been identified in yeast, including the key autophagic protein Atg8. The Atg8 gene family encodes ubiquitin-like proteins that share a similar structure consisting of two amino-terminal α helices and a ubiquitin-like core. Atg8 family members are expressed in various tissues, where they participate in multiple cellular processes, such as intracellular membrane trafficking and autophagy. Their role in autophagy has been intensively studied. Atg8 proteins undergo a unique ubiquitin-like conjugation to phosphatidylethanolamine on the autophagic membrane, a process essential for autophagosome formation. Whereas yeast has a single Atg8 gene, many other eukaryotes contain multiple Atg8 orthologs. Atg8 genes of multicellular animals can be divided, by sequence similarities, into three subfamilies: microtubule-associated protein 1 light chain 3 (MAP1LC3 or LC3), γ-aminobutyric acid receptor-associated protein (GABARAP) and Golgi-associated ATPase enhancer of 16 kDa (GATE-16), which are present in sponges, cnidarians (such as sea anemones, corals and hydras) and bilateral animals. Although genes from all three subfamilies are found in vertebrates, some invertebrate lineages have lost the genes from one or two subfamilies. The amino terminus of Atg8 proteins varies between the subfamilies and has a regulatory role in their various functions. Here we discuss the evolution of Atg8 proteins and summarize the current view of their function in intracellular trafficking and autophagy from a structural perspective.     

Elucidation of subfamily segregation and intramolecular coevolution of the olfactomedin-like proteins by comprehensive phylogenetic analysis and gene expression pattern assessment

The categorization of genes by structural distinctions relevant to biological characteristics is very important for understanding of gene functions and predicting functional implications of uncharacterized genes. It was absolutely necessary to deploy an effective and efficient strategy to deal with the complexity of the large olfactomedin-like (OLF) gene family sharing sequence similarity but playing diversified roles in many important biological processes, as the simple highest-hit homology analysis gave incomprehensive results and led to inappropriate annotation for some uncharacterized OLF members. In light of evolutionary information that may facilitate the classification of the OLF family and proper association of novel OLF genes with characterized homologs, we performed phylogenetic analysis on all 116 OLF proteins currently available, including two novel members cloned by our group. The OLF family segregated into seven subfamilies and members with similar domain compositions or functional properties all fell into relevant subfamilies. Furthermore, our Northern blot analysis and previous studies revealed that the typical human OLF members in each subfamily exhibited tissue-specific expression patterns, which in turn supported the segregation of the OLF subfamilies with functional divergence. Interestingly, the phylogenetic tree topology for the OLF domains alone was almost identical with that of the full-length tree representing the unique phylogenetic feature of full-length OLF proteins and their particular domain compositions. Moreover, each of the major functional domains of OLF proteins kept the same phylogenetic feature in defining similar topology of the tree. It indicates that the OLF domain and the various domains in flanking non-OLF regions have coevolved and are likely to be functionally interdependent. Expanded by a plausible gene duplication and domain couplings scenario, the OLF family comprises seven evolutionarily and functionally distinct subfamilies, in which each member shares similar structural and functional characteristics including the composition of coevolved and interdependent domains. The phylogenetically classified and preliminarily assessed subfamily framework may greatly facilitate the studying on the OLF proteins. Furthermore, it also demonstrated a feasible and reliable strategy to categorize novel genes and predict the functional implications of uncharacterized proteins based on the comprehensive phylogenetic classification of the subfamilies and their relevance to preliminary functional characteristics.     

Phylogenetic analysis of the eukaryotic RNA (cytosine-5)-methyltransferases

RNA (cytosine-5)-methyltransferases (RCMTs) have been characterized both in prokaryotic and eukaryotic organisms. The RCMT family, however, remains largely uncharacterized, as opposed to the family of DNA (cytosine-5)-methyltransferases which has been studied in depth. In the present study, an in silico identification of the putative 5-methylcytosine RNA-generating enzymes in the eukaryotic genomes was performed. A comprehensive phylogenetic analysis of the putative eukaryotic RCMT-related proteins has been performed in order to redefine subfamilies within the RCMT family. Five distinct eukaryotic subfamilies were identified, including the three already known (NOP2, NCL1 and YNL022c), one novel subfamily (RCMT9) and a fifth one which hitherto was considered to exist exclusively in prokaryotes (Fmu). The potential evolutionary relationships among the different eukaryotic RCMT subfamilies were also investigated. Furthermore, the results of this study add further support to a previous hypothesis that RCMTs represent evolutionary intermediates of RNA (uridine-5)-methyltransferases and DNA (cytosine-5)-methyltransferases.     

Phylogenomic analysis reveals dynamic evolutionary history of the Drosophila heterochromatin protein 1 (HP1) gene family

Heterochromatin is the gene-poor, satellite-rich eukaryotic genome compartment that supports many essential cellular processes. The functional diversity of proteins that bind and often epigenetically define heterochromatic DNA sequence reflects the diverse functions supported by this enigmatic genome compartment. Moreover, heterogeneous signatures of selection at chromosomal proteins often mirror the heterogeneity of evolutionary forces that act on heterochromatic DNA. To identify new such surrogates for dissecting heterochromatin function and evolution, we conducted a comprehensive phylogenomic analysis of the Heterochromatin Protein 1 gene family across 40 million years of Drosophila evolution. Our study expands this gene family from 5 genes to at least 26 genes, including several uncharacterized genes in Drosophila melanogaster. The 21 newly defined HP1s introduce unprecedented structural diversity, lineage-restriction, and germline-biased expression patterns into the HP1 family. We find little evidence of positive selection at these HP1 genes in both population genetic and molecular evolution analyses. Instead, we find that dynamic evolution occurs via prolific gene gains and losses. Despite this dynamic gene turnover, the number of HP1 genes is relatively constant across species. We propose that karyotype evolution drives at least some HP1 gene turnover. For example, the loss of the male germline-restricted HP1E in the obscura group coincides with one episode of dramatic karyotypic evolution, including the gain of a neo-Y in this lineage. This expanded compendium of ovary- and testis-restricted HP1 genes revealed by our study, together with correlated gain/loss dynamics and chromosome fission/fusion events, will guide functional analyses of novel roles supported by germline chromatin.     

Primate evolution of an olfactory receptor cluster: diversification by gene conversion and recent emergence of pseudogenes.

The olfactory receptor (OR) subgenome harbors the largest known gene family in mammals, disposed in clusters on numerous chromosomes. We have carried out a comparative evolutionary analysis of the best characterized genomic OR gene cluster, on human chromosome 17p13. Fifteen orthologs from chimpanzee (localized to chromosome 19p15), as well as key OR counterparts from other primates, have been identified and sequenced. Comparison among orthologs and paralogs revealed a multiplicity of gene conversion events, which occurred exclusively within OR subfamilies. These appear to lead to segment shuffling in the odorant binding site, an evolutionary process reminiscent of somatic combinatorial diversification in the immune system. We also demonstrate that the functional mammalian OR repertoire has undergone a rapid decline in the past 10 million years: while for the common ancestor of all great apes an intact OR cluster is inferred, in present-day humans and great apes the cluster includes nearly 40% pseudogenes.     

Coalescent methods for estimating species trees from phylogenomic data

Genome-scale sequence data have become increasingly available in the phylogenetic studies for understanding the evolutionary histories of species. However, it is challenging to develop probabilistic models to account for heterogeneity of phylogenomic data. The multispecies coalescent model describes gene trees as independent random variables generated from a coalescence process occurring along the lineages of the species tree. Since the multispecies coalescent model allows gene trees to vary across genes, coalescent-based methods have been popularly used to account for heterogeneous gene trees in phylogenomic data analysis. In this paper, we summarize and evaluate the performance of coalescent-based methods for estimating species trees from genome-scale sequence data. We investigate the effects of deep coalescence and mutation on the performance of species tree estimation methods. We found that the coalescent-based methods perform well in estimating species trees for a large number of genes, regardless of the degree of deep coalescence and mutation. The performance of the coalescent methods is negatively correlated with the lengths of internal branches of the species tree.     

Phylogenetic relationships within the Opisthokonta based on phylogenomic analyses of conserved single-copy protein domains

Many of the eukaryotic phylogenomic analyses published to date were based on alignments of hundreds to thousands of genes. Frequently, in such analyses, the most realistic evolutionary models currently available are often used to minimize the impact of systematic error. However, controversy remains over whether or not idiosyncratic gene family dynamics (i.e., gene duplications and losses) and incorrect orthology assignments are always appropriately taken into account. In this paper, we present an innovative strategy for overcoming orthology assignment problems. Rather than identifying and eliminating genes with paralogy problems, we have constructed a data set comprised exclusively of conserved single-copy protein domains that, unlike most of the commonly used phylogenomic data sets, should be less confounded by orthology miss-assignments. To evaluate the power of this approach, we performed maximum likelihood and Bayesian analyses to infer the evolutionary relationships within the opisthokonts (which includes Metazoa, Fungi, and related unicellular lineages). We used this approach to test 1) whether Filasterea and Ichthyosporea form a clade, 2) the interrelationships of early-branching metazoans, and 3) the relationships among early-branching fungi. We also assessed the impact of some methods that are known to minimize systematic error, including reducing the distance between the outgroup and ingroup taxa or using the CAT evolutionary model. Overall, our analyses support the Filozoa hypothesis in which Ichthyosporea are the first holozoan lineage to emerge followed by Filasterea, Choanoflagellata, and Metazoa. Blastocladiomycota appears as a lineage separate from Chytridiomycota, although this result is not strongly supported. These results represent independent tests of previous phylogenetic hypotheses, highlighting the importance of sophisticated approaches for orthology assignment in phylogenomic analyses.     

Evolutionary Genomics and Adaptive Evolution of the Hedgehog Gene Family (Shh,Ihh and Dhh) in Vertebrates

The Hedgehog (Hh) gene family codes for a class of secreted proteins composed of two active domains that act as signalling molecules during embryo development, namely for the development of the nervous and skeletal systems and the formation of the testis cord. While only one Hh gene is found typically in invertebrate genomes, most vertebrates species have three (Sonic hedgehog – Shh; Indian hedgehog – Ihh; and Desert hedgehog – Dhh), each with different expression patterns and functions, which likely helped promote the increasing complexity of vertebrates and their successful diversification. In this study, we used comparative genomic and adaptive evolutionary analyses to characterize the evolution of the Hh genes in vertebrates following the two major whole genome duplication (WGD) events. To overcome the lack of Hh-coding sequences on avian publicly available databases, we used an extensive dataset of 45 avian and three non-avian reptilian genomes to show that birds have all three Hh paralogs. We find suggestions that following the WGD events, vertebrate Hh paralogous genes evolved independently within similar linkage groups and under different evolutionary rates, especially within the catalytic domain. The structural regions around the ion-binding site were identified to be under positive selection in the signaling domain. These findings contrast with those observed in invertebrates, where different lineages that experienced gene duplication retained similar selective constraints in the Hh orthologs. Our results provide new insights on the evolutionary history of the Hh gene family, the functional roles of these paralogs in vertebrate species, and on the location of mutational hotspots.     

New insight into the evolution of aquaporins from flowering plants and vertebrates: orthologous identification and functional transfer is possible

Aquaporins (AQPs) represent a family of channel proteins that transport water and/or small solutes across cell membranes in the three domains of life. In all previous phylogenetic analysis of aquaporin, trees constructed using proteins with very low amino acid identity (<15%) were incongruent with rRNA data. In this work, restricting the evolutionary study of aquaporins to proteins with high amino acid identity (>25%), we showed congruence between AQPs and organismal trees. On the basis of this analysis, we defined 19 orthologous gene clusters in flowering plant species (3 PIP-like, 7 TIP-like, 6 NIP-like and 3 SIP-like). We described specific conserved motifs for each subfamily and each cluster, which were used to develop a method for automatic classification. Analysis of amino acid identity between orthologous monocotyledon and dicotyledon AQPs from each cluster, suggested that PIPs are under high evolutionary constraint. The phylogenetic analysis allowed us the assignment of orthologous aquaporins for very distant animal lineages (tetrapods-fishes). We also demonstrated that the location of all vertebrate AQPs in the ortholog clusters could be predicted by comparing their amino acid identity with human AQPs. We defined four AQP subfamilies in animals: AQP1-like, AQP8-like, AQP3-like and AQP11-like. Phylogenetic analysis showed that the four animal AQPs subfamilies are related with PIP-like, TIP-like, NIP-like and SIP-like subfamilies, respectively. Thus, this analysis would allow the prediction of individual AQPs function on the basis of orthologous genes from Arabidopsis thaliana and Homo sapiens.     

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

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

Molecular phylogenetics of the family Cyprinidae (Actinopterygii: Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichowensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages.     

The ABC gene family in arthropods: Comparative genomics and role in insecticide transport and resistance

About a 100 years ago, the Drosophila white mutant marked the birth of Drosophila genetics. The white gene turned out to encode the first well studied ABC transporter in arthropods. The ABC gene family is now recognized as one of the largest transporter families in all kingdoms of life. The majority of ABC proteins function as primary-active transporters that bind and hydrolyze ATP while transporting a large diversity of substrates across lipid membranes. Although extremely well studied in vertebrates for their role in drug resistance, less is known about the role of this family in the transport of endogenous and exogenous substances in arthropods. The ABC families of five insect species, a crustacean and a chelicerate have been annotated in some detail. We conducted a thorough phylogenetic analysis of the seven arthropod and human ABC protein subfamilies, to infer orthologous relationships that might suggest conserved function. Most orthologous relationships were found in the ABCB half transporter, ABCD, ABCE and ABCF subfamilies, but specific expansions within species and lineages are frequently observed and discussed. We next surveyed the role of ABC transporters in the transport of xenobiotics/plant allelochemicals and their involvement in insecticide resistance. The involvement of ABC transporters in xenobiotic resistance in arthropods is historically not well documented, but an increasing number of studies using unbiased differential gene expression analysis now points to their importance. We give an overview of methods that can be used to link ABC transporters to resistance. ABC proteins have also recently been implicated in the mode of action and resistance to Bt toxins in Lepidoptera. Given the enormous interest in Bt toxicology in transgenic crops, such findings will provide an impetus to further reveal the role of ABC transporters in arthropods.     

Evolution of Chemical Diversity in Echinocandin Lipopeptide Antifungal Metabolites

The echinocandins are a class of antifungal drugs that includes caspofungin, micafungin, and anidulafungin. Gene clusters encoding most of the structural complexity of the echinocandins provided a framework for hypotheses about the evolutionary history and chemical logic of echinocandin biosynthesis. Gene orthologs among echinocandin-producing fungi were identified. Pathway genes, including the nonribosomal peptide synthetases (NRPSs), were analyzed phylogenetically to address the hypothesis that these pathways represent descent from a common ancestor. The clusters share cooperative gene contents and linkages among the different strains. Individual pathway genes analyzed in the context of similar genes formed unique echinocandin-exclusive phylogenetic lineages. The echinocandin NRPSs, along with the NRPS from the inp gene cluster in Aspergillus nidulans and its orthologs, comprise a novel lineage among fungal NRPSs. NRPS adenylation domains from different species exhibited a one-to-one correspondence between modules and amino acid specificity that is consistent with models of tandem duplication and subfunctionalization. Pathway gene trees and Ascomycota phylogenies are congruent and consistent with the hypothesis that the echinocandin gene clusters have a common origin. The disjunct Eurotiomycete-Leotiomycete distribution appears to be consistent with a scenario of vertical descent accompanied by incomplete lineage sorting and loss of the clusters from most lineages of the Ascomycota. We present evidence for a single evolutionary origin of the echinocandin family of gene clusters and a progression of structural diversification in two fungal classes that diverged approximately 290 to 390 million years ago. Lineage-specific gene cluster evolution driven by selection of new chemotypes contributed to diversification of the molecular functionalities.