Effects of cis and trans regulatory variations on the expression divergence of heat shock response genes between yeast strains

Phenotypic variation among individuals in a population can be due to DNA sequence variation in protein coding regions or in regulatory elements. Recently, many studies have indicated that mutations in regulatory elements may be the major cause of phenotypic evolution. However, the mechanisms for evolutionary changes in gene expression are still not well understood. Here, we studied the relative roles of cis and trans regulatory changes in Saccharomyces cerevisiae cells to cope with heat stress. It has been found that the expression level of ~ 300 genes was induced at least two fold and that of ~ 500 genes was repressed at least two fold in response to heat shock. From the former set of genes, we randomly selected 65 genes that showed polymorphism(s) between the BY and RM strains for pyrosequencing analysis to explore the relative contributions of cis and trans regulatory variations to the expression divergence between BY and RM. Our data indicated that the expression divergence between BY and RM was mainly due to trans regulatory variations under either the normal condition or the heat stress condition. However, the relative contribution of trans regulatory variation was decreased from 76.9% to 61.5% after the heat shock stress. These results indicated that the cis regulatory variation may play an important role in the adaption to heat stress. In our data, 43.1% (28 genes) of the 65 genes showed the same trend of cis or trans variation effect after the heat shock stress, 35.4% (23 genes) showed an increased cis variation effect and 21.5% (14 genes) showed an increased trans variation effect after the heat shock stress. Thus, our data give insights into the relative roles of cis and trans variations in response to heat shock in yeast.     

Comparative genomic analysis of the Sm gene family in rice and maize

Sm proteins are a group of ubiquitous ring-shaped oligomers that function in multiple aspects of RNA metabolism. However, until this study, no comprehensive study incorporating phylogeny, chromosomal location, gene organization, adaptive evolution, expression profiling and functional networks has been reported for rice and maize. In this study, twenty-five and thirty-three Sm genes have been identified in rice and maize, respectively. Phylogenetic analyses identified eighteen gene groups. Results by gene locations indicated that segmental duplication contributes to the expansion of this gene family in rice and maize. Gene organization and motif compositions of the Sm members are highly conserved in each group, indicative of their functional conservation. Expression profiles have provided insights into the possible functional divergence among members of the Sm gene family. Adaptive evolution analyses suggested that purifying selection was the main force driving Sm evolution, but some critical sites might be responsible for functional divergence. In addition, four hundred and seventy-nine interactions were identified by functional network analyses, and most of which were associated with binding, cellular macromolecule biosynthesis, pre-mRNA processing and transferase activity. Overall, the data contribute to a better understanding of the complexity of Sm gene family in rice and maize and will provide a solid foundation for future functional studies.     

Cloning and expression of boule and dazl in the Nile tilapia (Oreochromis niloticus)

The Deleted in Azoospermia (DAZ) family of RNA binding proteins consists of highly conserved genes boule, daz and daz-like (dazl) essential for germ cell development. boule is known for its unisexual meiotic expression in invertebrates and mammals, but meiotic-specific female expression plus meiosis-preferential male expression in trout, and meiosis-preferential bisexual expression in medaka. dazl shows highly conserved bisexual expression throughout gametogenesis in diverse species. Here we report the cloning and expression of boule and dazl in the Nile tilapia (Oreochromis niloticus), an important aquaculture fish. Molecular cloning and sequence analysis led to the identification of tilapia boule and dazl cDNAs. The predicted partial Boule contains a conserved RRM motif and Dazl has the C-terminal sequence. On a phylogenetic tree, tilapia Boule and Dazl are in separate clades of Boule and Dazl homologs from other species, indicating their divergence during early vertebrate evolution. By RT-PCR analysis, boule and dazl showed bisexual gonad-specific expression. By in situ hybridization analysis, both boule and dazl RNAs were restricted to female and male germ cells of adult gonads but absent in gonadal soma. In the ovary, boule and dazl RNAs were abundant in oocytes. In the testis, boule and dazl RNAs were prominent in meiotic spermatocytes but barely detectable in meiotic products. These data show that boule and dazl are expressed bisexually in germ cells and provide useful markers to study gametogenesis in the adult tilapia.     

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.     

Phylogeny of the macaques (Cercopithecidae: Macaca) based on Alu elements

Genus Macaca (Cercopithecidae: Papionini) is one of the most successful primate radiations. Despite previous studies on morphology and mitochondrial DNA analysis, a number of issues regarding the details of macaque evolution remain unsolved. Alu elements are a class of non-autonomous retroposons belonging to short interspersed elements that are specific to the primate lineage. Because retroposon insertions show very little homoplasy, and because the ancestral state (absence of the SINE) is known, Alu elements are useful genetic markers and have been utilized for analyzing primate phylogenentic relationships and human population genetic relationships. Using PCR display methodology, 298 new Alu insertions have been identified from ten species of macaques. Together with 60 loci reported previously, a total of 358 loci are used to infer the phylogenetic relationships of genus Macaca. With regard to earlier unresolved issues on the macaque evolution, the topology of our tree suggests that: 1) genus Macaca contains four monophyletic species groups; 2) within the Asian macaques, the silenus group diverged first, and members of the sinica and fascicularis groups share a common ancestor; 3) Macaca arctoides are classified in the sinica group. Our results provide a robust molecular phylogeny for genus Macaca with stronger statistical support than previous studies. The present study also illustrates that SINE-based approaches are a powerful tool in primate phylogenetic studies and can be used to successfully resolve evolutionary relationships between taxa at scales from the ordinal level to closely related species within one genus.     

Evolved Lactococcus lactis Strains for Enhanced Expression of Recombinant Membrane Proteins

The production of complex multidomain (membrane) proteins is a major hurdle in structural genomics and a generic approach for optimizing membrane protein expression is still lacking. We have devised a selection method to isolate mutant strains with improved functional expression of recombinant membrane proteins. By fusing green fluorescent protein and an erythromycin resistance marker (ErmC) to the C-terminus of a target protein, one simultaneously selects for variants with enhanced expression (increased erythromycin resistance) and correct folding (green fluorescent protein fluorescence). Three evolved hosts, displaying 2- to 8-fold increased expression of a plethora of proteins, were fully sequenced and shown to carry single-site mutations in the nisK gene. NisK is the sensor protein of a two-component regulatory system that directs nisin-A-mediated expression. The levels of recombinant membrane proteins were increased in the evolved strains, and in some cases their folding states were improved. The generality and simplicity of our approach allow rapid improvements of protein production yields by directed evolution in a high-throughput way.     

Divergence of the mitochondrial electron transport chains from the green alga Chlamydomonas reinhardtii and its colorless close relative Polytomella sp.

Compelling evidence exists that the colorless algae of the genus Polytomella arose from a green Chlamydomonas-like ancestor by losing its functional photosynthetic apparatus. Due to the close relationship between the colorless and the green chlorophyte, Polytomella sp. appeared as a useful indicative framework for structural studies of Chlamydomonas reinhardtii mitochondria. However, comparative studies reported here unexpectedly revealed significant differences between the mitochondrial respiratory systems of the two algae. Two-dimensional blue native/SDS-PAGE of isolated mitochondria indicated that cytochrome-containing respiratory complexes III and IV in the two chlorophytes contrast in size, subunit composition and relative abundance. Complex IV in Polytomella is smaller than its counterpart in C. reinhardtii and occurs in two forms that differ presumably in the presence of subunit COXIII. The cytochrome c and the iron-sulfur Rieske protein of both chlorophytes revealed structural differences on the amino acid sequence level. Under comparable culture conditions, the colorless alga exhibits lower levels of cytochrome c and complex IV but a higher respiratory activity than the green alga. Cytochrome c levels were also found to be differently regulated by the growth conditions in both algae. The divergence between the respiratory systems in the two related chlorophytes can be viewed as a consequence of the loss of photosynthetic activity and/or of the adaptation to the environment via the acquisition of a more flexible, heterotrophic metabolism. Our understanding of mitochondrial function and evolution is expected to be greatly enhanced via further parallel studies of photosynthetic/non-photosynthetic algae, for which this study forms an incentive.     

Divergence of the vertebrate sp1A/ryanodine receptor domain and SOCS box-containing (Spsb) gene family and its expression and regulation within the mouse brain

The Spsb family of genes encode well-conserved proteins of unknown function. Mammalian Spsb genes are likely the result of three separate duplication and divergence events during vertebrate evolution. The phylogenetic relationship along with expression and regulation of Spsb genes may offer insight into the evolution and function of this gene family in vertebrates. We have established that Spsb genes are expressed in numerous tissues, however their pattern and level of expression is tissue-dependent. Further, only Spsb1 is responsive to stress caused by ethanol exposure in the mouse brain, which suggests that Spsb genes have acquired different regulatory mechanisms. Analysis of cis-regulatory elements supports this, but also reveals some common regulatory modules involved in cell proliferation and stress response. Our results contribute to the growing body of data on the expression and function of Spsb genes, which serve as a model for studies on the origin, divergence and specialization of eukaryotic gene families.     

The Determinants of Stability and Folding in Evolutionarily Diverged Cytochromes c

Cytochrome c has served as a paradigm for the study of protein stability, folding, and molecular evolution, but it remains unclear how these aspects of the protein are related. For example, while the bovine and equine cytochromes c are known to have different stabilities, and possibly different folding mechanisms, it is not known how these differences arise from just three amino acid substitutions introduced during divergence. Using site-selectively incorporated carbon-deuterium bonds, we show that like the equine protein, bovine cytochrome c is induced to unfold by guanidine hydrochloride via a stepwise mechanism, but it does not populate an intermediate as is observed with the equine protein. The increased stability also results in more similar free energies of unfolding observed at different sites within the protein, giving the appearance of a more concerted mechanism. Furthermore, we show that the differences in stability and folding appear to result from a single amino acid substitution that stabilizes a helix by allowing for increased solvation of its N-terminus.     

Transfection of Eimeria and Toxoplasma using heterologous regulatory sequences

Eimeriatenella and Toxoplasmagondii are Apicomplexan protozoa and share many similarities in biology and genomics. While the latter parasites are easily cultured in vitro and genetically manipulated, many Eimeria species are difficult to grow in vitro. We hypothesised that molecular tools for the genetic manipulation of T. gondii could be applied to the study of Eimeria parasites. Here we show that three different promoter sequences originating from E. tenella could function effectively not only in other species of the Eimeria genus (histone H4) but also in T. gondii (histone H4, actin and tubulin). Similarly, promoters of the “housekeeping” gene (tubulin) and differentially regulated gene (surface antigen gene, sag1) of T. gondii were effective in driving the expression of the yellow fluorescent protein (YFP) maker gene in E. tenella. The transfection efficiency with heterologous regulatory sequences was similar to that with homologous promoters; while the promoter strength of heterologous vectors is slightly weaker than the homologous vectors in both E. tenella and T. gondii. The results suggest that 5′ regulatory sequences are functionally conserved not only among the Eimeria species, but also between T. gondii and E. tenella, and that T. gondii could be used as a novel transfection check system for Eimeria-rooted vectors, accelerating the development of reverse genetics in Eimeria spp.     

Complete chloroplast genome sequences of Praxelis (Eupatorium catarium Veldkamp), an important invasive species

Praxelis (Eupatorium catarium Veldkamp) is a new hazardous invasive plant species that has caused serious economic losses and environmental damage in the Northern hemisphere tropical and subtropical regions. Although previous studies focused on detecting the biological characteristics of this plant to prevent its expansion, little effort has been made to understand the impact of Praxelis on the ecosystem in an evolutionary process. The genetic information of Praxelis is required for further phylogenetic identification and evolutionary studies. Here, we report the complete Praxelis chloroplast (cp) genome sequence. The Praxelis chloroplast genome is 151,410 bp in length including a small single-copy region (18,547 bp) and a large single-copy region (85,311 bp) separated by a pair of inverted repeats (IRs; 23,776 bp). The genome contains 85 unique and 18 duplicated genes in the IR region. The gene content and organization are similar to other Asteraceae tribe cp genomes. We also analyzed the whole cp genome sequence, repeat structure, codon usage, contraction of the IR and gene structure/organization features between native and invasive Asteraceae plants, in order to understand the evolution of organelle genomes between native and invasive Asteraceae. Comparative analysis identified the 14 markers containing greater than 2% parsimony-informative characters, indicating that they are potential informative markers for barcoding and phylogenetic analysis. Moreover, a sister relationship between Praxelis and seven other species in Asteraceae was found based on phylogenetic analysis of 28 protein-coding sequences. Complete cp genome information is useful for plant phylogenetic and evolutionary studies within this invasive species and also within the Asteraceae family.     

The evolution and functional divergence of the beta-carotene oxygenase gene family in teleost fish—Exemplified by Atlantic salmon

In mammals, two carotenoid cleaving oxygenases are known; beta-carotene 15,15′-monooxygenase (BCMO1) and beta-carotene 9′,10′-oxygenase (BCO2). BCMO1 is a key enzyme in vitamin A synthesis by symmetrically cleaving beta-carotene into 2 molecules of all-trans-retinal, while BCO2 is responsible for asymmetric cleavage of a broader range of carotenoids. Here, we show that the Atlantic salmon beta-carotene oxygenase (bco) gene family contains 5 members, three bco2 and two bcmo1 paralogs. Using public sequence databases, multiple bco genes were also found in several additional teleost species. Phylogenetic analysis indicates that bco2a and bco2b originate from the teleost fish specific genome duplication (FSGD or 3R), while the third and more distant paralog, bco2 like, might stem from a prior duplication event in the teleost lineage. The two bcmo1 paralogs (bcmo1 and bcmo1 like) appear to be the result of an ancient duplication event that took place before the divergence of ray-finned (Actinopterygii) and lobe-finned fish (Sarcopterygii), with subsequent nonfunctionalization and loss of one Sarcopterygii paralog. Gene expression analysis of the bcmo1 and bco2 paralogs in Atlantic salmon reveals regulatory divergence with tissue specific expression profiles, suggesting that the beta-carotene oxygenase subtypes have evolved functional divergences. We suggest that teleost fish have evolved and maintained an extended repertoire of beta-carotene oxygenases compared to the investigated Sarcopterygii species, and hypothesize that the main driver behind this functional divergence is the exposure to a diverse set of carotenoids in the aquatic environment.     

Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants

Background and Aims

The OVATE gene encodes a nuclear-localized regulatory protein belonging to a distinct family of plant-specific proteins known as the OVATE family proteins (OFPs). OVATE was first identified as a key regulator of fruit shape in tomato, with nonsense mutants displaying pear-shaped fruits. However, the role of OFPs in plant development has been poorly characterized.

Methods

Public databases were searched and a total of 265 putative OVATE protein sequences were identified from 13 sequenced plant genomes that represent the major evolutionary lineages of land plants. A phylogenetic analysis was conducted based on the alignment of the conserved OVATE domain from these 13 selected plant genomes. The expression patterns of tomato SlOFP genes were analysed via quantitative real-time PCR. The pattern of OVATE gene duplication resulting in the expansion of the gene family was determined in arabidopsis, rice and tomato.

Key Results

Genes for OFPs were found to be present in all the sampled land plant genomes, including the early-diverged lineages, mosses and lycophytes. Phylogenetic analysis based on the amino acid sequences of the conserved OVATE domain defined 11 sub-groups of OFPs in angiosperms. Different evolutionary mechanisms are proposed for OVATE family evolution, namely conserved evolution and divergent expansion. Characterization of the AtOFP family in arabidopsis, the OsOFP family in rice and the SlOFP family in tomato provided further details regarding the evolutionary framework and revealed a major contribution of tandem and segmental duplications towards expansion of the OVATE gene family.

Conclusions

This first genome-wide survey on OFPs provides new insights into the evolution of the OVATE protein family and establishes a solid base for future functional genomics studies on this important but poorly characterized regulatory protein family in plants.