The solute carrier (SLC) complement of the human genome: phylogenetic classification reveals four major families

Solute carriers (SLCs) is the largest group of transporters, embracing transporters for inorganic ions, amino acids, neurotransmitters, sugars, purines and fatty acids among other substrates. We mined the finished assembly of the human genome using Hidden Markov Models (HMMs) obtaining a total of 384 unique SLC sequences. Detailed clustering and phylogenetic analysis of the entire SLC family showed that 15 of the families place into four large phylogenetic clusters with the largest containing eight SLC families, suggesting that many of the distinct families of SLCs have a common evolutionary origin. This study represents the first overall genomic roadmap of the SLCs providing large sequence sets and clarifies the phylogenetic relationships among the families of the second largest group of membrane proteins.     

Long evolutionary conservation and considerable tissue specificity of several atypical solute carrier transporters

The superfamily of Solute Carriers (SLCs) has around 384 members in the human genome grouped into at least 48 families. While many of these transporters have been well characterized with established important biological functions, there are few recently identified genes that are not studied regarding tissue distribution or evolutionary origin. Here we study 14 of these recently discovered SLC genes (HIAT1, HIATL1, MFSD1, MFSD5, MFSD6, MFSD9, MFSD10, SLC7A14, SLC7A15, SLC10A6, SLC15A5, SLC16A12, SLC30A10 and SLC21A21) with the purpose to give much better picture over the sequence relationship and tissue expression of the diverse SLC gene family. We used a range of bioinformatic methods to classify each of these genes into the different SLC gene families. We found that 9 of the 14 atypical SLCs are distant members of the Major Facilitator Superfamily (MFS) clan while the others belong to the APC clan, the DMT clan, the CPA_AT clan and the IT clan. We found most of the genes to be highly evolutionary conserved, likely to be present in most bilateral species, except for SLC21A21 that we found only present in mammals. Several of these transporter genes have highly specific tissue expression profile while it is notable that most are expressed in the CNS with the exception of SLC21A21 and SLC15A5. This work provides fundamental information on 14 transporters that previously have not received much attention enabling a more comprehensive view over the SLC superfamily.     

Evolutionary genomics of the Fox genes: Origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.     

Mito-nuclear genetic comparison in a Wolbachia infected weevil: insights on reproductive mode, infection age and evolutionary forces shaping genetic variation

Background

The composition and expression of vertebrate gene families is shaped by species specific gene loss in combination with a number of gene and genome duplication events (R1, R2 in all vertebrates, R3 in teleosts) and depends on the ecological and evolutionary context. In this study we analyzed the evolutionary history of the solute carrier 1 (SLC1) gene family. These genes are supposed to be under strong selective pressure (purifying selection) due to their important role in the timely removal of glutamate at the synapse.

Results

In a genomic survey where we manually annotated and analyzing sequences from more than 300 SLC1 genes (from more than 40 vertebrate species), we found evidence for an interesting evolutionary history of this gene family. While human and mouse genomes contain 7 SLC1 genes, in prototheria, sauropsida, and amphibia genomes up to 9 and in actinopterygii up to 13 SLC1 genes are present. While some of the additional slc1 genes in ray-finned fishes originated from R3, the increased number of SLC1 genes in prototheria, sauropsida, and amphibia genomes originates from specific genes retained in these lineages. Phylogenetic comparison and microsynteny analyses of the SLC1 genes indicate, that theria genomes evidently lost several SLC1 genes still present in the other lineage. The genes lost in theria group into two new subfamilies of the slc1 gene family which we named slc1a8/eaat6 and slc1a9/eaat7.

Conclusions

The phylogeny of the SLC1/EAAT gene family demonstrates how multiple genome reorganization and duplication events can influence the number of active genes. Inactivation and preservation of specific SLC1 genes led to the complete loss of two subfamilies in extant theria, while other vertebrates have retained at least one member of two newly identified SLC1 subfamilies.     

Fourteen novel human members of mitochondrial solute carrier family 25 (SLC25) widely expressed in the central nervous system

Members of the solute carrier family 25 (SLC25) are known to transport molecules over the mitochondrial membrane. In this paper we present 14 novel members of SLC25 family in human. These were provided with following gene symbols by the HGNC: SLC25A32, SLC25A33, SLC25A34, SLC25A35, SLC25A37, SLC25A38, SLC25A39, SLC25A40, SLC25A41, SLC25A42, SLC25A43, SLC25A44, SLC25A45, and SLC25A46. We also identified the orthologues for these genes in rat and mouse. Moreover, we found yeast orthologues for 9 of these genes and show that the predicted substrate binding residues are highly conserved in the human and yeast proteins. We performed a comprehensive tissue localization study for 9 of these genes on a panel of 30 rat tissues with quantitative real-time polymerse chain reaction. We detected their mRNA in a wide number of tissues, both in brain and in periphery. This study provides an overall roadmap of the repertoire of the SLC25 family in mammals, showing that there are at least 46 genes in the human genome coding for mitochondrial transporters.     

In silico identification and Bayesian phylogenetic analysis of multiple new mammalian kallikrein gene families

Kallikrein gene families have been identified previously in genomes of the human, the mouse, and the rat, and individual kallikrein-like genes have been found in many more species. This study presents the in silico identification of kallikrein gene families in the recently sequenced genomes of four additional mammalian species, the chimpanzee, the dog, the pig, and the opossum. Phylogenies were constructed with gene sequences from all seven mammalian families, using Bayesian analysis, which clarified the evolutionary relationships between these genes. Individual gene sequences, as well as concatenated constructs of multiple sequences, were used. Fifteen kallikrein genes were located in the chimpanzee (Pan troglodytes) genome, while only 14 were identified in the canine (Canis familiaris) genome as no orthologue to human KLK3 was found. Thirteen genes were identified from the pig (Sus scrofa) genome, which lacked homologues to KLK2 and KLK3, and 11 genes, orthologous to human KLK5 through KLK15, were found in the opossum (Monodelphis domestica) genome. No kallikrein genes were identified from the available genome sequences of the chicken (Gallus gallus) or African clawed frog (Xenopus tropicalis). Within the family of kallikreins several subfamilies were suggested by phylogenetic analysis. One consisted of KLK4, KLK5, and KLK14; another of KLK9, KLK11, and KLK15; a third of KLK10 and KLK12; a fourth of KLK6 and KLK13; and finally one of KLK8 and the classical kallikreins (KLK1, KLK2, and KLK3).     

The repertoire of solute carriers of family 6: identification of new human and rodent genes

Tremendous amount of primary sequence information has been made available from the genome sequencing projects, although a complete annotation and identification of all genes is still far from being complete. Here, we present the identification of two new human genes from the pharmacologically important family of transporter proteins, solute carriers family 6 (SLC6). These were named SLC6A17 and SLC6A18 by HUGO. The human repertoire of SLC6 proteins now consists of 19 functional members and four pseudogenes. We also identified the corresponding orthologues and additional genes from mouse and rat genomes. Detailed phylogenetic analysis of the entire family of SLC6 proteins in mammals shows that this family can be divided into four subgroups. We used Hidden Markov Models for these subgroups and identified in total 430 unique SLC6 proteins from 10 animal, one plant, two fungi, and 196 bacterial genomes. It is evident that SLC6 proteins are present in both animals and bacteria, and that three of the four subfamilies of mammalian SLC6 proteins are present in Caenorhabditis elegans, showing that these subfamilies are evolutionary very ancient. Moreover, we performed tissue localization studies on the entire family of SLC6 proteins on a panel of 15 rat tissues and further, the expression of three of the new genes was studied using quantitative real-time PCR showing expression in multiple central and peripheral tissues. This paper presents an overall overview of the gene repertoire of the SLC6 gene family and its expression profile in rats.     

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.     

Poxvirus orthologous clusters: toward defining the minimum essential poxvirus genome

Increasingly complex bioinformatic analysis is necessitated by the plethora of sequence information currently available. A total of 21 poxvirus genomes have now been completely sequenced and annotated, and many more genomes will be available in the next few years. First, we describe the creation of a database of continuously corrected and updated genome sequences and an easy-to-use and extremely powerful suite of software tools for the analysis of genomes, genes, and proteins. These tools are available free to all researchers and, in most cases, alleviate the need for using multiple Internet sites for analysis. Further, we describe the use of these programs to identify conserved families of genes (poxvirus orthologous clusters) and have named the software suite POCs, which is available at www.poxvirus.org. Using POCs, we have identified a set of 49 absolutely conserved gene families-those which are conserved between the highly diverged families of insect-infecting entomopoxviruses and vertebrate-infecting chordopoxviruses. An additional set of 41 gene families conserved in chordopoxviruses was also identified. Thus, 90 genes are completely conserved in chordopoxviruses and comprise the minimum essential genome, and these will make excellent drug, antibody, vaccine, and detection targets. Finally, we describe the use of these tools to identify necessary annotation and sequencing updates in poxvirus genomes. For example, using POCs, we identified 19 genes that were widely conserved in poxviruses but missing from the vaccinia virus strain Tian Tan 1998 GenBank file. We have reannotated and resequenced fragments of this genome and verified that these genes are conserved in Tian Tan. The results for poxvirus genes and genomes are discussed in light of evolutionary processes.     

Molecular evolution of the ankyrin gene family

Ankyrins are membrane adaptor molecules that play important roles in coupling integral membrane proteins to the spectrin-based cytoskeleton network. Human mutations of ankyrin genes lead to severe genetic diseases such as fatal cardiac arrhythmias and hereditary spherocytosis. To elucidate the evolutionary history of ankyrins, we have identified novel ankyrin sequences in insect, fish, frog, chicken, dog, and chimpanzee genomes and explored the phylogenetic relationships of the ankyrin gene family. Our data demonstrate that duplication of ankyrin genes occurred at two different stages. The first duplication resulted from an independent evolution event specific in Arthropoda after its divergence from Chordata. Following the separation from Urochordata, expansion of ankyrins in vertebrates involved ancestral genome duplications. We did not find evidence of coordinated arrangements of gene families of ankyrin-associated membrane proteins on paralogous chromosomes. In addition, evolution of the 24 ANK-repeats strikingly correlated with the exon boundary sites of ankyrin genes, which might have occurred before its duplication in vertebrates. Such correlation is speculated to bring functional diversity and complexity. Moreover, based on the phylogenetic analysis of the ANK-repeat domain, we put forward a novel model for the putative primordial ankyrin that contains the fourth six-ANK-repeat subdomain and the spectrin-binding domain. These findings will provide guides for future studies concerning structure, function, evolutionary origins of ankyrins, and possibly other cytoskeletal proteins.     

Genome Wide Identification,Phylogeny and Expression of Zinc Transporter Genes in Common Carp

Background

Zinc is an essential trace element in organisms, which serves as a cofactor for hundreds of enzymes that are involved in many pivotal biological processes including growth, development, reproduction and immunity. Therefore, the homeostasis of zinc in the cell is fundamental. The zinc transporter gene family is a large gene family that encodes proteins which regulate the movement of zinc across cellular and intracellular membranes. However, studies on teleost zinc transporters are mainly limited to model species.

Methodology/Principal Findings

We identified a set of 37 zinc transporters in common carp genome, including 17 from SLC30 family (ZnT), and 20 from SLC39 family (ZIP). Phylogenetic and syntenic analysis revealed that most of the zinc transporters are highly conserved, though recent gene duplication and gene losses do exist. Through examining the copy number of zinc transporter genes across several vertebrate genomes, thirteen zinc transporters in common carp are found to have undergone the gene duplications, including SLC30A1, SLC30A2, SLC30A5, SLC30A7, SLC30A9, SLC30A10, SLC39A1, SLC39A3, SLC39A4, SLC39A5, SLC39A6, SLC39A7 and SLC39A9. The expression patterns of all zinc transporters were established in various tissues, including blood, brain, gill, heart, intestine, liver, muscle, skin, spleen and kidney, and showed that most of the zinc transporters were ubiquitously expressed, indicating the critical role of zinc transporters in common carp.

Conclusions

To some extent, examination of gene families with detailed phylogenetic or orthology analysis could verify the authenticity and accuracy of assembly and annotation of the recently published common carp whole genome sequences. The gene families are also considered as a unique source for evolutionary studies. Moreover, the whole set of common carp zinc transporters provides an important genomic resource for future biochemical, toxicological and physiological studies of zinc in teleost.     

Gene recruitment--a common mechanism in the evolution of transfer RNA gene families

The evolution of alloacceptor transfer RNAs (tRNAs) has been traditionally thought to occur vertically and reflect the evolution of the genetic code. Yet there have been several indications that a tRNA gene could evolve horizontally, from a copy of an alloacceptor tRNA gene in the same genome. Earlier, we provided the first unambiguous evidence for the occurrence of such "tRNA gene recruitment" in nature--in the mitochondrial (mt) genome of the demosponge Axinella corrugata. Yet the extent and the pattern of this process in the evolution of tRNA gene families remained unclear. Here we analyzed tRNA genes from 21 mt genomes of demosponges as well as nuclear genomes of rhesus macaque, chimpanzee and human. We found four new cases of alloacceptor tRNA gene recruitment in mt genomes and eleven cases in the nuclear genomes. In most of these cases we observed a single nucleotide substitution at the middle position of the anticodon, which resulted in the change of not only the tRNA's amino-acid identity but also the class of the amino-acyl tRNA synthetases (aaRSs) involved in amino-acylation. We hypothesize that the switch to a different class of aaRSs may have prevented the conflict between anticodon and amino-acid identities of recruited tRNAs. Overall our results suggest that gene recruitment is a common phenomenon in tRNA multigene family evolution and should be taken into consideration when tRNA evolutionary history is reconstructed.     

Structural biology of solute carrier (SLC) membrane transport proteins

The human solute carriers (SLCs) comprise over 400 different transporters, organized into 65 families (http://slc.bioparadigms.org/) based on their sequence homology and transport function. SLCs are responsible for transporting extraordinarily diverse solutes across biological membranes, including inorganic ions, amino acids, lipids, sugars, neurotransmitters and drugs. Most of these membrane proteins function as coupled symporters (co-transporters) utilizing downhill ion (H⁺ or Na⁺) gradients as the driving force for the transport of substrate against its concentration gradient into cells. Other members work as antiporters (exchangers) that typically contain a single substrate-binding site with an alternating access mode of transport, while a few members exhibit channel-like properties. Dysfunction of SLCs is correlated with numerous human diseases and therefore they are potential therapeutic drug targets. In this review, we identified all of the SLC crystal structures that have been determined, most of which are from prokaryotic species. We further sorted all the SLC structures into four main groups with different protein folds and further discuss the well-characterized MFS (major facilitator superfamily) and LeuT (leucine transporter) folds. This review provides a systematic analysis of the structure, molecular basis of substrate recognition and mechanism of action in different SLC family members.     

A comprehensive expression analysis of all members of a gene family encoding cell-wall enzymes allowed us to predict cis-regulatory regions involved in cell-wall construction in specific organs of Arabidopsis.

The Arabidopsis thaliana genome sequencing project has revealed that multigene families, such as those generated by genome duplications, are more abundant among plant genomes than among animal genomes. To gain insight into the evolutionary implications of the multigene families in higher plants, we examined the XTH gene family, a group of genes encoding xyloglucan endotransglucosylase/hydrolase, which are responsible for cell-wall construction in plants. Expression analysis of all members (33 genes) of this family, using quantitative real-time RT-PCR, revealed that most members exhibit distinct expression profiles in terms of tissue specificity and responses to hormonal signals, with some members exhibiting similar expression patterns. By comparing the flanking sequences of individual genes, we identified four sets of large-segment duplications and two sets of solitary gene duplications. In each set of gene duplicates, long nucleotide sequences, ranging from one to two hundred base pairs, are conserved. Furthermore, gene duplicates exhibit similar organ-specific expression profiles. These facts allowed us to predict putative cis-regulatory regions, particularly those responsible for cell-wall construction, and hence for morphogenesis, that are specific for certain organs or tissues in plants.     

Using paleogenomics to study the evolution of gene families: origin and duplication history of the relaxin family hormones and their receptors

Recent progress in the analysis of whole genome sequencing data has resulted in the emergence of paleogenomics, a field devoted to the reconstruction of ancestral genomes. Ancestral karyotype reconstructions have been used primarily to illustrate the dynamic nature of genome evolution. In this paper, we demonstrate how they can also be used to study individual gene families by examining the evolutionary history of relaxin hormones (RLN/INSL) and relaxin family peptide receptors (RXFP). Relaxin family hormones are members of the insulin superfamily, and are implicated in the regulation of a variety of primarily reproductive and neuroendocrine processes. Their receptors are G-protein coupled receptors (GPCR's) and include members of two distinct evolutionary groups, an unusual characteristic. Although several studies have tried to elucidate the origins of the relaxin peptide family, the evolutionary origin of their receptors and the mechanisms driving the diversification of the RLN/INSL-RXFP signaling systems in non-placental vertebrates has remained elusive. Here we show that the numerous vertebrate RLN/INSL and RXFP genes are products of an ancestral receptor-ligand system that originally consisted of three genes, two of which apparently trace their origins to invertebrates. Subsequently, diversification of the system was driven primarily by whole genome duplications (WGD, 2R and 3R) followed by almost complete retention of the ligand duplicates in most vertebrates but massive loss of receptor genes in tetrapods. Interestingly, the majority of 3R duplicates retained in teleosts are potentially involved in neuroendocrine regulation. Furthermore, we infer that the ancestral AncRxfp3/4 receptor may have been syntenically linked to the AncRln-like ligand in the pre-2R genome, and show that syntenic linkages among ligands and receptors have changed dynamically in different lineages. This study ultimately shows the broad utility, with some caveats, of incorporating paleogenomics data into understanding the evolution of gene families.     

Reconstructing ancestral genome content based on symmetrical best alignments and Dollo parsimony

MOTIVATION: Gene duplications and losses (GDLs) are important events in genome evolution. They result in expansion or contraction of gene families, with a likely role in phenotypic evolution. As more genomes become available and their annotations are improved, software programs capable of rapidly and accurately identifying the content of ancestral genomes and the timings of GDLs become necessary to understand the unique evolution of each lineage. RESULTS: We report EvolMAP, a new algorithm and software that utilizes a species tree-based gene clustering method to join all-to-all symmetrical similarity comparisons of multiple gene sets in order to infer the gene composition of multiple ancestral genomes. The algorithm further uses Dollo parsimony-based comparison of the inferred ancestral genes to pinpoint the timings of GDLs onto evolutionary intervals marked by speciation events. Using EvolMAP, first we analyzed the expansion of four families of G-protein coupled receptors (GPCRs) within animal lineages. Additional to demonstrating the unique expansion tree for each family, results also show that the ancestral eumetazoan genome contained many fewer GPCRs than modern animals, and these families expanded through concurrent lineage-specific duplications. Second, we analyzed the history of GDLs in mammalian genomes by comparing seven proteomes. In agreement with previous studies, we report that the mammalian gene family sizes have changed drastically through their evolution. Interestingly, although we identified a potential source of duplication for 75% of the gained genes, remaining 25% did not have clear-cut sources, revealing thousands of genes that have likely gained their distinct sequence identities within the descent of mammals. AVAILABILITY: Query server, source code and executable are available at http://kosik-web.mcdb.ucsb.edu/evolmap/index.htm .     

The Secretin GPCRs descended from the family of Adhesion GPCRs

The Adhesion G-protein-coupled receptors (GPCRs) are the most complex gene family among GPCRs with large genomic size, multiple introns, and a fascinating flora of functional domains, though the evolutionary origin of this family has been obscure. Here we studied the evolution of all class B (7tm2)-related genes, including the Adhesion, Secretin, and Methuselah families of GPCRs with a focus on nine genomes. We found that the cnidarian genome of Nematostella vectensis has a remarkably rich set of Adhesion GPCRs with a broad repertoire of N-terminal domains although this genome did not have any Secretin GPCRs. Moreover, the single-celled and colony-forming eukaryotes Monosiga brevicollis and Dictyostelium discoideum contain Adhesion-like GPCRs although these genomes do not have any Secretin GPCRs suggesting that the Adhesion types of GPCRs are the most ancient among class B GPCRs. Phylogenetic analysis found Adhesion group V (that contains GPR133 and GPR144) to be the closest relative to the Secretin family in the Adhesion family. Moreover, Adhesion group V sequences in N. vectensis share the same splice site setup as the Secretin GPCRs. Additionally, one of the most conserved motifs in the entire Secretin family is only found in group V of the Adhesion family. We suggest therefore that the Secretin family of GPCRs could have descended from group V Adhesion GPCRs. We found a set of unique Adhesion-like GPCRs in N. vectensis that have long N-termini containing one Somatomedin B domain each, which is a domain configuration similar to that of a set of Adhesion-like GPCRs found in Branchiostoma floridae. These sequences show slight similarities to Methuselah sequences found in insects. The extended class B GPCRs have a very complex evolutionary history with several species-specific expansions, and we identified at least 31 unique N-terminal domains originating from other protein classes. The overall N-terminal domain structure, however, concurs with the phylogenetic analysis of the transmembrane domains, thus enabling us to track the origin of most of the subgroups.     

Poxvirus genome evolution by gene gain and loss

The poxviruses (Poxviridae) are a family of viruses with double-stranded DNA genomes and substantial numbers (often >200) of genes per genome. We studied the patterns of gene gain and loss over the evolutionary history of 17 poxvirus complete genomes. A phylogeny based on gene family presence/absence showed good agreement with families based on concatenated amino acid sequences of conserved single-copy genes. Gene duplications in poxviruses were often lineage specific, and the most extensively duplicated viral gene families were found in only a few of the genomes analyzed. A total of 34 gene families were found to include a member in at least one of the poxvirus genomes analyzed and at least one animal genome; in 16 (47%) of these families, there was evidence of recent horizontal gene transfer (HGT) from host to virus. Gene families with evidence of HGT included several involved in host immune defense mechanisms (the MHC class I, interleukin-10, interleukin-24, interleukin-18, the interferon gamma receptor, and tumor necrosis factor receptor II) and others (glutaredoxin and glutathione peroxidase) involved in resistance of cells to oxidative stress. Thus "capture" of host genes by HGT has been a recurrent feature of poxvirus evolution and has played an important role in adapting the virus to survive host antiviral defense mechanisms.