The Tic20 gene family: Phylogenetic analysis and evolutionary considerations

Tic20 is a polytopic protein of the inner envelope membrane of chloroplasts, and it is proposed to act as a translocation channel during chloroplast protein import. By analyzing 29 sequences from diverse organisms, it was evident that Tic20-related proteins form two distinct clades, termed Group 1 and Group 2. The former group includes canonical Tic20 proteins that are essential for chloroplast development, while members of the latter are of unknown function. An increased evolutionary rate, in connection with adaptation to terrestrial life, was detected in Group 1. Interestingly, the sub-cellular (genomic) localization of genes coding for Group 1 proteins differs between evolutionary lineages.Key words: Arabidopsis, chloroplast protein import, plastids, protein targeting, Tic20     

Phylogenetic and evolutionary analysis of the PLUNC gene family

The PLUNC family of human proteins are candidate host defense proteins expressed in the upper airways. The family subdivides into short (SPLUNC) and long (LPLUNC) proteins, which contain domains predicted to be structurally similar to one or both of the domains of bactericidal/permeability-increasing protein (BPI), respectively. In this article we use analysis of the human, mouse, and rat genomes and other sequence data to examine the relationships between the PLUNC family proteins from humans and other species, and between these proteins and members of the BPI family. We show that PLUNC family clusters exist in the mouse and rat, with the most significant diversification in the locus occurring for the short PLUNC family proteins. Clear orthologous relationships are established for the majority of the proteins, and ambiguities are identified. Completion of the prediction of the LPLUNC4 proteins reveals that these proteins contain approximately a 150-residue insertion encoded by an additional exon. This insertion, which is predicted to be largely unstructured, replaces the structure homologous to the 40s hairpin of BPI. We show that the exon encoding this region is anomalously variable in size across the LPLUNC proteins, suggesting that this region is key to functional specificity. We further show that the mouse and human PLUNC family orthologs are evolving rapidly, which supports the hypothesis that these proteins are involved in host defense. Intriguingly, this rapid evolution between the human and mouse sequences is replaced by intense purifying selection in a large portion of the N-terminal domain of LPLUNC4. Our data provide a basis for future functional studies of this novel protein family.     

Cloning and evolutionary analysis of tobacco MAPK gene family

The mitogen-activated protein (MAP) kinase cascade is an important signaling module which is involved in biotic and abiotic stress responses as well as plant growth and development. In this study, we identified 17 tobacco MAPKs including 11 novel tobacco MAPK genes that have not been identified before. Comparative analysis with MAPK gene families from other plants, such as Athaliana thaliana, rice and poplar, suggested that tobacco MAPKs (such as NtMPK1, NtMPK3 and NtMPK8) might play similar functions in response to abiotic and biotic stresses. QRT-PCR analysis revealed that a total of 14 NtMPKs were regulated by SA and/or MeJA, suggesting their potential roles involved in plant defense response. In addition, 6 NtMPKs were induced by drought treatment, implying their roles in response to drought stress. Our results indicated that most of tobacco MAPK might be involved in plant defense response, which provides the basis for further analysis on physiological functions of tobacco MAPKs.     

Characterization and comparative analysis of the EGLN gene family

Rat Sm-20 is a homologue of the Caenorhabditis elegans gene egl-9 and has been implicated in the regulation of growth, differentiation and apoptosis in muscle and nerve cells. Null mutants in egl-9 result in a complete tolerance to an otherwise lethal toxin produced by Pseudomonas aeruginosa. This study describes the conserved Egl-Nine (EGLN) gene family of which rat SM-20 and C. elegans Egl-9 are members and characterizes the mouse and human homologues. Each of the human genes (EGLN1, EGLN2 and EGLN3) are of a conserved genomic structure consisting of five coding exons. Phylogenetic analysis and domain organization show that EGLN1 represents the ancestral form of the gene family and that EGLN3 is the human orthologue of rat Sm-20. The previously observed mitochondrial targeting of rat SM-20 is unlikely to be a general feature of the protein family and may be a feature specific to rats. An EGLN gene is unexpectedly found in the genome of P. aeruginosa, a bacterium known to produce a toxin that acts through the Egl-9 protein. The pathogenic bacterium Vibrio cholerae is also shown to have an EGLN gene suggesting that it is an important pathogenicity factor. These results provide new insights into host-pathogen interactions and a basis for further functional characterization of the gene family and resolve discrepancies in annotation between gene family members.     

The eta-globin gene. Its long evolutionary history in the beta-globin gene family of mammals

In phylogenetic reconstructions by the parsimony method, utilizing 62 sequenced globin genes and pseudogenes (including 34 of the beta-globin gene family from eutherian orders Primates, Lagomorpha, Artiodactyla and Rodentia), the branch of primate psi beta pseudogenes and the goat embryonically expressed epsilon II gene group monophyletically together as orthologues of a common ancestral gene (labelled eta) distinct from orthologues of epsilon, gamma, delta and beta. This primate psi eta-goat eta branch is cladistically closer to epsilon and gamma than to delta and beta branches. In each eutherian order gene conversions replaced portions of delta by beta sequences, whereas in descent of Primates epsilon, gamma and eta mostly retained their separate ancient identities predating the radiation of Eutheria in all their exons and non-coding regions. The loci of the ancestral beta-globin gene cluster in basal eutherians and proto-primates, as deduced from beta-clusters representing the four eutherian orders, were linked 5'-epsilon-gamma-eta-delta-beta-3' with epsilon, gamma and eta being embryonically expressed genes, and delta and beta ontogenetically later expressed genes. Through deletions gamma was lost in artiodactyl evolution, eta in lagomorph and rodent evolution, and all DNA between exon 2 3' boundaries of eta and delta in prosimian lemuriform evolution (lemur having the hybrid pseudogene psi eta delta). Simian primates retained intact the five loci of the ancestral cluster. Not only did eta, after it became a pseudogene in the basal primates, persist intact in descent to present-day simians but in the line to hominoids it evolved during the last 40 million years at the decelerated rate of 1 X 10(-9) substitutions/site per year which is one-fifth the expected neutral rate. The possibility is suggested that the psi eta locus situated between fetal and adult chromosomal domains of the simian beta-globin gene cluster might play some role in a mechanism for ontogenetic switches of globin gene expression. However, not enough sequence data on genes and intergenic regions in DNA of species of primates and other mammals as yet exist to know if the slow rate of 1 X 10(-9) reflects the rate of a conserved functional gene or primarily reflects a decelerated neutral rate of hominoid DNA evolution, conceivably from enhanced DNA repair and longer generation times in hominoids. The further possibility is raised that gene correction (repair of damaged DNA that prevents emergence of new alleles) and gene conversion both more often involve strand copying of conserved than of rapidly evolving DNA.     

Immunoglobulin CH gene family in hominoids and its evolutionary history.

The organization of the human immunoglobulin CH gene suggests that a gene duplication involving the C gamma-C gamma-C epsilon-C alpha region has occurred during evolution. We previously showed that both chimpanzee and gorilla have two 5'-C epsilon-C alpha-3', as in human, and that orangutan, gibbon, and Old World monkeys have one C epsilon gene and one, two, and one C alpha gene(s), respectively. In addition to these clustered CH genes, there is one processed C epsilon pseudogene in each species. The present study revealed that orangutan and crab-eating macaque (an Old World monkey) both have one 5'-C epsilon-C alpha-3' and that gibbon has two 5'-C epsilon-C alpha-3', one C epsilon gene of which is completely deleted. By Southern analysis, the number of C gamma genes in all the nonhuman hominoids was estimated to be four to five, as in human, in comparison with two for crab-eating macaque. The C mu and C delta genes were estimated to be present as single copies in both hominoids and crab-eating macaque. Furthermore, it was proved that there are two copies of the C epsilon 5'-flanking region in both the orangutan and the gibbon genomes. These results show that gene duplication including the C gamma-C gamma-C epsilon-C alpha genes occurred in the common ancestor of hominoids and that subsequent deletion of the C epsilon gene (in orangutan, including one of the C alpha genes) occurred independently in each hominoid species.     

Embryonic expression and evolutionary analysis of the amphioxus Dickkopf and Kremen family genes

The secreted Wnt signaling inhibitor Dickkopf1(Dkk1)plays key role in vertebrate head induction.Its receptor Kremen synergizes with Dkkl in Wnt inhibition.Here we have carried out expression and functional studies of the Dkk and Kremen genes in amphioxus(Branchiostoma belcheri).During embryonic and larval development,BbDkk1/2/4 is expressed in the posterior mesoendoderm,anterior somatic mesoderm and the pharyngeal regions.Its expression becomes restricted to the pharyngeal region on the left side at larval stages.In 45 h larvae,BbDkk1/2/4 is expressed specifically in the cerebral vesicle.BbDkk3 was only detected at larval stages in the mid-intestine region.Seven Kremen related genes were identified in the genome of the Florida amphioxus(Branchiostoma floridae),clustered in 4scaffolds,and are designated Kremen1-4 and Kremen-like 1-3,respectively.In B.belcheri,Kremenl is strongly expressed in the mesoendoderm during early development and Kremen3 is expressed asymmetrically in spots in the larval pharyngeal region.In luciferase reporter assays,BbDkk1/2/4 can strongly inhibit Writ signaling,while BbDkk3,BbKremen1 and BbKremen3 can not.No co-operative effect was observed between amphioxus Dkk1/2/4 and Kremens,suggesting that the interaction between Dkk and Kremen likely originated later during evolution.     

The extracellular Leucine-Rich Repeat superfamily; a comparative survey and analysis of evolutionary relationships and expression patterns

Correction to Dolan J, Walshe K, Alsbury S, Hokamp K, O'Keeffe S, Okafuji T, Miller SF, Tear G, Mitchell KJ: The extracellular leucine-rich repeat superfamily; a comparative survey and analysis of evolutionary relationships and expression patterns. BMC Genomics 2007, 8:320.     

Genome-wide comparative gene family classification

Correct classification of genes into gene families is important for understanding gene function and evolution. Although gene families of many species have been resolved both computationally and experimentally with high accuracy, gene family classification in most newly sequenced genomes has not been done with the same high standard. This project has been designed to develop a strategy to effectively and accurately classify gene families across genomes. We first examine and compare the performance of computer programs developed for automated gene family classification. We demonstrate that some programs, including the hierarchical average-linkage clustering algorithm MC-UPGMA and the popular Markov clustering algorithm TRIBE-MCL, can reconstruct manual curation of gene families accurately. However, their performance is highly sensitive to parameter setting, i.e. different gene families require different program parameters for correct resolution. To circumvent the problem of parameterization, we have developed a comparative strategy for gene family classification. This strategy takes advantage of existing curated gene families of reference species to find suitable parameters for classifying genes in related genomes. To demonstrate the effectiveness of this novel strategy, we use TRIBE-MCL to classify chemosensory and ABC transporter gene families in C. elegans and its four sister species. We conclude that fully automated programs can establish biologically accurate gene families if parameterized accordingly. Comparative gene family classification finds optimal parameters automatically, thus allowing rapid insights into gene families of newly sequenced species.     

Combining gene expression QTL mapping and phenotypic spectrum analysis to uncover gene regulatory relationships

Gene expression QTL (eQTL) mapping can suggest candidate regulatory relationships between genes. Recent advances in mammalian phenotype annotation such as mammalian phenotype ontology (MPO) enable systematic analysis of the phenotypic spectrum subserved by many genes. In this study we combined eQTL mapping and phenotypic spectrum analysis to predict gene regulatory relationships. Five pairs of genes with similar phenotypic effects and potential regulatory relationships suggested by eQTL mapping were identified. Lines of evidence supporting some of the predicted regulatory relationships were obtained from biological literature. A particularly notable example is that promoter sequence analysis and real-time PCR assays support the predicted regulation of protein kinase C epsilon (Prkce) by cAMP responsive element binding protein 1 (Creb1). Our results show that the combination of gene eQTL mapping and phenotypic spectrum analysis may provide a valuable approach to uncovering gene regulatory relations underlying mammalian phenotypes.     

Sequence features and evolutionary mechanisms in the chicken avidin gene family.

The chicken avidin gene family comprises the avidin gene (avd) and several homologous avidin-related genes (avrs). The sequences of the avr genes are nearly identical to each other but exhibit nonrandomly distributed, frequently nonsynonymous nucleotide substitutions compared to avd. In this study, we determined the genetic distances and the phylogeny of the avd and avr genes and found differences between different exons and introns. Our results suggest the involvement of biased gene conversion in the evolution of the genes. Furthermore, one of the genes was identified as a putative fusion gene. The occurrence of both gene conversion and recombination supports the models suggesting a common initiation mechanism for conversion and crossing-over. The existence of avidin-related proteins (AVRs) is currently unknown, but the putative AVRs are expected to bind biotin similarly to avidin. However, the observed sequence differences may affect the stability and glycosylation patterns of the putative AVR proteins.     

The rabbitCD1 and the evolutionary conservation of theCD1 gene family

A comparison of the genes encoding the CD1 leucocyte differentiation antigens in man and mouse shows important differences which prompted us to analyze theCD1 genes of the rabbit. We have found that the rabbit genome contains multipleCD1 loci. Upon cloning and sequencing, one of these loci was found to encode the known rabbit CD1-like antigen (R-Ta) and to be closely related to the humanCD1b gene, which is absent in the mouse, while a second rabbit gene is closely related to both the humanR3 and the mouseCD1 genes. The data reinforce the notion of the existence of two classes ofCD1 genes, one of which is conserved in all species, while the other, albeit also evolutionarily old, has been deleted in mice as well as in other rodents. The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M26248 and M26249.     

The SAUR gene family in coffee: genome-wide identification and gene expression analysis during somatic embryogenesis

Molecular Biology Reports - Small auxin-up RNA (SAUR) genes form a wide family supposedly involved in different physiological and developmental processes in plants such as leaf senescence, auxin...     

葡萄TST基因家族鉴定及表达分析

【目的】液泡膜糖转运蛋白（tonoplast sugar transporter,TST或者TMT）是植物发育过程中发挥重要作用的一种糖转运蛋白。为探究该基因家族在葡萄生长发育中的作用,并进一步为阐明TST基因功能提供坚实的基础。【方法】通过同源分析法从葡萄基因组中鉴定出13个TST基因,对基因的结构和编码蛋白质进行生物信息学分析;利用qRT-PCR技术分析‘鄞红’葡萄发育过程中不同组织的TST表达水平,并与不同时期葡萄果肉中可溶性糖含量进行相关性分析。【结果】结果表明：该基因家族分布在6条染色体上,存在3对片段重复和3对串联重复基因;根据系统发育将其分为3个亚家族,各亚族家族成员结构相似;顺式作用元件表明TST基因含有大量与激素、光和胁迫响应相关的顺式作用元件;蛋白质结构均显示该家族由α-螺旋和无规则卷曲组成,各亚族蛋白模型相似;qRT-PCR结果显示,VvTST在不同组织中均有表达,并存在时空表达特异性;对葡萄果肉中基因表达量变化与可溶性糖含量变化进行相关性分析发现,4个VvTST（VIT_18s0001g12560、VIT_18s0122g00850、VIT_04s0023g01860和VIT_03s0038g03940）的表达水平与葡萄果肉中可溶性糖的积累呈现相似趋势。【结论】上述研究结果表明,VvTST可能对葡萄果肉中可溶性糖积累起到至关重要的作用。