A role for gene duplication and natural variation of gene expression in the evolution of metabolism

Background

Most eukaryotic genomes have undergone whole genome duplications during their evolutionary history. Recent studies have shown that the function of these duplicated genes can diverge from the ancestral gene via neo- or sub-functionalization within single genotypes. An additional possibility is that gene duplicates may also undergo partitioning of function among different genotypes of a species leading to genetic differentiation. Finally, the ability of gene duplicates to diverge may be limited by their biological function.

Methodology/Principal Findings

To test these hypotheses, I estimated the impact of gene duplication and metabolic function upon intraspecific gene expression variation of segmental and tandem duplicated genes within Arabidopsis thaliana. In all instances, the younger tandem duplicated genes showed higher intraspecific gene expression variation than the average Arabidopsis gene. Surprisingly, the older segmental duplicates also showed evidence of elevated intraspecific gene expression variation albeit typically lower than for the tandem duplicates. The specific biological function of the gene as defined by metabolic pathway also modulated the level of intraspecific gene expression variation. The major energy metabolism and biosynthetic pathways showed decreased variation, suggesting that they are constrained in their ability to accumulate gene expression variation. In contrast, a major herbivory defense pathway showed significantly elevated intraspecific variation suggesting that it may be under pressure to maintain and/or generate diversity in response to fluctuating insect herbivory pressures.

Conclusion

These data show that intraspecific variation in gene expression is facilitated by an interaction of gene duplication and biological activity. Further, this plays a role in controlling diversity of plant metabolism.     

Early evolutionary history and genomic features of gene duplicates in the human genome

Background

Human gene duplicates have been the focus of intense research since the development of array-based and targeted next-generation sequencing approaches in the last decade. These studies have primarily concentrated on determining the extant copy-number variation from a population-genomic perspective but lack a robust evolutionary framework to elucidate the early structural and genomic characteristics of gene duplicates at emergence and their subsequent evolution with increasing age.

Results

We analyzed 184 gene duplicate pairs comprising small gene families in the draft human genome with 10 % or less synonymous sequence divergence. Human gene duplicates primarily originate from DNA-mediated events, taking up genomic residence as intrachromosomal copies in direct or inverse orientation. The distribution of paralogs on autosomes follows random expectations in contrast to their significant enrichment on the sex chromosomes. Furthermore, human gene duplicates exhibit a skewed gradient of distribution along the chromosomal length with significant clustering in pericentromeric regions. Surprisingly, despite the large average length of human genes, the majority of extant duplicates (83 %) are complete duplicates, wherein the entire ORF of the ancestral copy was duplicated. The preponderance of complete duplicates is in accord with an extremely large median duplication span of 36 kb, which enhances the probability of capturing ancestral ORFs in their entirety. With increasing evolutionary age, human paralogs exhibit declines in (i) the frequency of intrachromosomal paralogs, and (ii) the proportion of complete duplicates. These changes may reflect lower survival rates of certain classes of duplicates and/or the role of purifying selection. Duplications arising from RNA-mediated events comprise a small fraction (11.4 %) of all human paralogs and are more numerous in older evolutionary cohorts of duplicates.

Conclusions

The degree of structural resemblance, genomic location and duplication span appear to influence the long-term maintenance of paralogs in the human genome. The median duplication span in the human genome far exceeds that in C. elegans and yeast and likely contributes to the high prevalence of complete duplicates relative to structurally heterogeneous duplicates (partial and chimeric). The relative roles of regulatory sequence versus exon-intron structure changes in the acquisition of novel function by human paralogs remains to be determined.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1827-3) contains supplementary material, which is available to authorized users.     

Characterization of the caleosin gene family in the Triticeae

Background

The caleosin genes encode proteins with a single conserved EF hand calcium-binding domain and comprise small gene families found in a wide range of plant species. Some members of the gene family have been shown to be upregulated by environmental stresses including low water availability and high salinity. Caleosin 3 from wheat has been shown to interact with the α-subunit of the heterotrimeric G proteins, and to act as a GTPase activating protein (GAP). This study characterizes the size and diversity of the gene family in wheat and related species and characterizes the differential tissue-specific expression of members of the gene family.

Results

A total of 34 gene family members that belong to eleven paralogous groups of caleosins were identified in the hexaploid bread wheat, T. aestivum. Each group was represented by three homeologous copies of the gene located on corresponding homeologous chromosomes, except the caleosin 10, which has four gene copies. Ten gene family members were identified in diploid barley, Hordeum vulgare, and in rye, Secale cereale, seven in Brachypodium distachyon, and six in rice, Oryza sativa. The analysis of gene expression was assayed in triticale and rye by RNA-Seq analysis of 454 sequence sets and members of the gene family were found to have diverse patterns of gene expression in the different tissues that were sampled in rye and in triticale, the hybrid hexaploid species derived from wheat and rye. Expression of the gene family in wheat and barley was also previously determined by microarray analysis, and changes in expression during development and in response to environmental stresses are presented.

Conclusions

The caleosin gene family had a greater degree of expansion in the Triticeae than in the other monocot species, Brachypodium and rice. The prior implication of one member of the gene family in the stress response and heterotrimeric G protein signaling, points to the potential importance of the caleosin gene family. The complexity of the family and differential expression in various tissues and under conditions of abiotic stress suggests the possibility that caleosin family members may play diverse roles in signaling and development that warrants further investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-239) contains supplementary material, which is available to authorized users.     

Relationship between gene duplicability and diversifiability in the topology of biochemical networks

Background

Selective gene duplicability, the extensive expansion of a small number of gene families, is universal. Quantitatively, the number of genes (P_(K)) with K duplicates in a genome decreases precipitously as K increases, and often follows a power law (P_(k)∝k^-α). Functional diversification, either neo- or sub-functionalization, is a major evolution route for duplicate genes.

Results

Using three lines of genomic datasets, we studied the relationship between gene duplicability and diversifiability in the topology of biochemical networks. First, we explored scenario where two pathways in the biochemical networks antagonize each other. Synthetic knockout of respective genes for the two pathways rescues the phenotypic defects of each individual knockout. We identified duplicate gene pairs with sufficient divergences that represent this antagonism relationship in the yeast S. cerevisiae. Such pairs overwhelmingly belong to large gene families, thus tend to have high duplicability. Second, we used distances between proteins of duplicate genes in the protein interaction network as a metric of their diversification. The higher a gene’s duplicate count, the further the proteins of this gene and its duplicates drift away from one another in the networks, which is especially true for genetically antagonizing duplicate genes. Third, we computed a sequence-homology-based clustering coefficient to quantify sequence diversifiability among duplicate genes – the lower the coefficient, the more the sequences have diverged. Duplicate count (K) of a gene is negatively correlated to the clustering coefficient of its duplicates, suggesting that gene duplicability is related to the extent of sequence divergence within the duplicate gene family.

Conclusion

Thus, a positive correlation exists between gene diversifiability and duplicability in the context of biochemical networks – an improvement of our understanding of gene duplicability.     

Genome-Wide Analysis of the Expansin Gene Superfamily Reveals Grapevine-Specific Structural and Functional Characteristics

Background

Expansins are proteins that loosen plant cell walls in a pH-dependent manner, probably by increasing the relative movement among polymers thus causing irreversible expansion. The expansin superfamily (EXP) comprises four distinct families: expansin A (EXPA), expansin B (EXPB), expansin-like A (EXLA) and expansin-like B (EXLB). There is experimental evidence that EXPA and EXPB proteins are required for cell expansion and developmental processes involving cell wall modification, whereas the exact functions of EXLA and EXLB remain unclear. The complete grapevine (Vitis vinifera) genome sequence has allowed the characterization of many gene families, but an exhaustive genome-wide analysis of expansin gene expression has not been attempted thus far.

Methodology/Principal Findings

We identified 29 EXP superfamily genes in the grapevine genome, representing all four EXP families. Members of the same EXP family shared the same exon–intron structure, and phylogenetic analysis confirmed a closer relationship between EXP genes from woody species, i.e. grapevine and poplar (Populus trichocarpa), compared to those from Arabidopsis thaliana and rice (Oryza sativa). We also identified grapevine-specific duplication events involving the EXLB family. Global gene expression analysis confirmed a strong correlation among EXP genes expressed in mature and green/vegetative samples, respectively, as reported for other gene families in the recently-published grapevine gene expression atlas. We also observed the specific co-expression of EXLB genes in woody organs, and the involvement of certain grapevine EXP genes in berry development and post-harvest withering.

Conclusion

Our comprehensive analysis of the grapevine EXP superfamily confirmed and extended current knowledge about the structural and functional characteristics of this gene family, and also identified properties that are currently unique to grapevine expansin genes. Our data provide a model for the functional characterization of grapevine gene families by combining phylogenetic analysis with global gene expression profiling.     

Genome-wide and molecular evolution analysis of the subtilase gene family in Vitis vinifera

Background

Vitis vinifera (grape) is one of the most economically significant fruit crops in the world. The availability of the recently released grape genome sequence offers an opportunity to identify and analyze some important gene families in this species. Subtilases are a group of subtilisin-like serine proteases that are involved in many biological processes in plants. However, no comprehensive study incorporating phylogeny, chromosomal location and gene duplication, gene organization, functional divergence, selective pressure and expression profiling has been reported so far for the grape.

Results

In the present study, a comprehensive analysis of the subtilase gene family in V. vinifera was performed. Eighty subtilase genes were identified. Phylogenetic analyses indicated that these subtilase genes comprised eight groups. The gene organization is considerably conserved among the groups. Distribution of the subtilase genes is non-random across the chromosomes. A high proportion of these genes are preferentially clustered, indicating that tandem duplications may have contributed significantly to the expansion of the subtilase gene family. Analyses of divergence and adaptive evolution show that while purifying selection may have been the main force driving the evolution of grape subtilases, some of the critical sites responsible for the divergence may have been under positive selection. Further analyses of real-time PCR data suggested that many subtilase genes might be important in the stress response and functional development of plants.

Conclusions

Tandem duplications as well as purifying and positive selections have contributed to the functional divergence of subtilase genes in V. vinifera. The data may contribute to a better understanding of the grape subtilase gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1116) contains supplementary material, which is available to authorized users.     

Phylogeny of Oedogoniales, Chaetophorales and Chaetopeltidales (Chlorophyceae): inferences from sequence-structure analysis of ITS2

Background and Aims

The green algal class Chlorophyceae comprises five orders (Chlamydomonadales, Sphaeropleales, Chaetophorales, Chaetopeltidales and Oedogoniales). Attempts to resolve the relationships among these groups have met with limited success. Studies of single genes (18S rRNA, 26S rRNA, rbcL or atpB) have largely failed to unambiguously resolve the relative positions of Oedogoniales, Chaetophorales and Chaetopeltidales (the OCC taxa). In contrast, recent genomics analyses of plastid data from OCC exemplars provided a robust phylogenetic analysis that supports a monophyletic OCC alliance.

Methods

An ITS2 data set was assembled to independently test the OCC hypothesis and to evaluate the performance of these data in assessing green algal phylogeny at the ordinal or class level. Sequence-structure analysis designed for use with ITS2 data was employed for phylogenetic reconstruction.

Key Results

Results of this study yielded trees that were, in general, topologically congruent with the results from the genomic analyses, including support for the monophyly of the OCC alliance.

Conclusions

Not all nodes from the ITS2 analyses exhibited robust support, but our investigation demonstrates that sequence-structure analyses of ITS2 provide a taxon-rich means of testing phylogenetic hypotheses at high taxonomic levels. Thus, the ITS2 data, in the context of sequence-structure analysis, provide an economical supplement or alternative to the single-marker approaches used in green algal phylogeny.     

Cyanolichens can have both cyanobacteria and green algae in a common layer as major contributors to photosynthesis

Background and Aims

Cyanolichens are usually stated to be bipartite (mycobiont plus cyanobacterial photobiont). Analyses revealed green algal carbohydrates in supposedly cyanobacterial lichens (in the genera Pseudocyphellaria, Sticta and Peltigera). Investigations were carried out to determine if both cyanobacteria and green algae were present in these lichens and, if so, what were their roles.

Methods

The types of photobiont present were determined by light and fluorescence microscopy. Small carbohydrates were analysed to detect the presence of green algal metabolites. Thalli were treated with selected strengths of Zn²⁺ solutions that stop cyanobacterial but not green algal photosynthesis. CO₂ exchange was measured before and after treatment to determine the contribution of each photobiont to total thallus photosynthesis. Heterocyst frequencies were determined to clarify whether the cyanobacteria were modified for increased nitrogen fixation (high heterocyst frequencies) or were normal, vegetative cells.

Key Results

Several cyanobacterial lichens had green algae present in the photosynthetic layer of the thallus. The presence of the green algal transfer carbohydrate (ribitol) and the incomplete inhibition of thallus photosynthesis upon treatment with Zn²⁺ solutions showed that both photobionts contributed to the photosynthesis of the lichen thallus. Low heterocyst frequencies showed that, despite the presence of adjacent green algae, the cyanobacteria were not altered to increase nitrogen fixation.

Conclusions

These cyanobacterial lichens are a tripartite lichen symbiont combination in which the mycobiont has two primarily photosynthetic photobionts, ‘co-primary photobionts’, a cyanobacterium (dominant) and a green alga. This demonstrates high flexibility in photobiont choice by the mycobiont in the Peltigerales. Overall thallus appearance does not change whether one or two photobionts are present in the cyanobacterial thallus. This suggests that, if there is a photobiont effect on thallus structure, it is not specific to one or the other photobiont.     

Endogenous cytokinin in developing kiwifruit is implicated in maintaining fruit flesh chlorophyll levels

Background and Aims

Green kiwifruit (Actinidia deliciosa) retain high concentrations of chlorophyll in the fruit flesh, whereas in gold-fleshed kiwifruit (A. chinensis) chlorophyll is degraded to colourless catabolites during fruit development, leaving yellow carotenoids visible. The plant hormone group the cytokinins has been implicated in the delay of senescence, and so the aim of this work was to investigate the link between cytokinin levels in ripening fruit and chlorophyll de-greening.

Methods

The expression of genes related to cytokinin metabolism and signal transduction and the concentration of cytokinin metabolites were measured. The regulation of gene expression was assayed using transient activation of the promoter of STAY-GREEN2 (SGR2) by cytokinin response regulators.

Key Results

While the total amount of cytokinin increased in fruit of both species during maturation and ripening, a high level of expression of two cytokinin biosynthetic gene family members, adenylate isopentenyltransferases, was only detected in green kiwifruit fruit during ripening. Additionally, high levels of O-glucosylated cytokinins were detected only in green kiwifruit, as was the expression of the gene for zeatin O-glucosyltransferase, the enzyme responsible for glucosylating cytokinin into a storage form. Season to season variation in gene expression was seen, and some de-greening of the green kiwifruit fruit occurred in the second season, suggesting environmental effects on the chlorophyll degradation pathway. Two cytokinin-related response regulators, RRA17 and RRB120, showed activity against the promoter of kiwifruit SGR2.

Conclusions

The results show that in kiwifruit, levels of cytokinin increase markedly during fruit ripening, and that cytokinin metabolism is differentially regulated in the fruit of the green and gold species. However, the causal factor(s) associated with the maintenance or loss of chlorophyll in kiwifruit during ripening remains obscure.     

Evolutionary Patterns and Selective Pressures of Odorant/Pheromone Receptor Gene Families in Teleost Fishes

Background

Teleost fishes do not have a vomeronasal organ (VNO), and their vomeronasal receptors (V1Rs, V2Rs) are expressed in the main olfactory epithelium (MOE), as are odorant receptors (ORs) and trace amine-associated receptors (TAARs). In this study, to obtain insights into the functional distinction among the four chemosensory receptor families in teleost fishes, their evolutionary patterns were examined in zebrafish, medaka, stickleback, fugu, and spotted green pufferfish.

Methodology/Principal Findings

Phylogenetic analysis revealed that many lineage-specific gene gains and losses occurred in the teleost fish TAARs, whereas only a few gene gains and losses have taken place in the teleost fish vomeronasal receptors. In addition, synonymous and nonsynonymous nucleotide substitution rate ratios (K_A/K_S) in TAARs tended to be higher than those in ORs and V2Rs.

Conclusions/Significance

Frequent gene gains/losses and high K_A/K_S in teleost TAARs suggest that receptors in this family are used for detecting some species-specific chemicals such as pheromones. Conversely, conserved repertoires of V1R and V2R families in teleost fishes may imply that receptors in these families perceive common odorants for teleosts, such as amino acids. Teleost ORs showed intermediate evolutionary pattern between TAARs and vomeronasal receptors. Many teleost ORs seem to be used for common odorants, but some ORs may have evolved to recognize lineage-specific odors.     

Oxidative stress induces distinct physiological responses in the two Trebouxia phycobionts of the lichen Ramalina farinacea

Background and Aims

Most lichens form associations with Trebouxia phycobionts and some of them simultaneously include genetically different algal lineages. In other symbiotic systems involving algae (e.g. reef corals), the relative abundances of different endosymbiotic algal clades may change over time. This process seems to provide a mechanism allowing the organism to respond to environmental stress. A similar mechanism may operate in lichens with more than one algal lineage, likewise protecting them against environmental stresses. Here, the physiological responses to oxidative stress of two distinct Trebouxia phycobionts (provisionally named TR1 and TR9) that coexist within the lichen Ramalina farinacea were analysed.

Methods

Isolated phycobionts were exposed to oxidative stress through the reactive oxygen species propagator cumene hydroperoxide (CuHP). Photosynthetic pigments and proteins, photosynthesis (through modulated chlorophyll fluorescence), the antioxidant enzymes superoxide dismutase (SOD) and glutathione reductase (GR), and the stress-related protein HSP70 were analysed.

Key Results

Photosynthetic performance was severely impaired by CuHP in phycobionts, as indicated by decreases in the maximal PSII photochemical efficiency (F_v/F_m), the quantum efficiency of PSII (Φ_PSII) and the non-photochemical dissipation of energy (NPQ). However, the CuHP-dependent decay in photosynthesis was significantly more severe in TR1, which also showed a lower NPQ and a reduced ability to preserve chlorophyll a, carotenoids and D1 protein. Additionally, differences were observed in the capacities of the two phycobionts to modulate antioxidant activities and HPS70 levels when exposed to oxidative stress. In TR1, CuHP significantly diminished HSP70 and GR but did not change SOD activities. In contrast, in TR9 the levels of both antioxidant enzymes and those of HSP70 increased in response to CuHP.

Conclusions

The better physiological performance of TR9 under oxidative conditions may reflect its greater capacity to undertake key metabolic adjustments, including increased non-photochemical quenching, higher antioxidant protection and the induction of repair mechanisms.     

The Arabidopsis Kinome: phylogeny and evolutionary insights into functional diversification

Background

Protein kinases constitute a particularly large protein family in Arabidopsis with important functions in cellular signal transduction networks. At the same time Arabidopsis is a model plant with high frequencies of gene duplications. Here, we have conducted a systematic analysis of the Arabidopsis kinase complement, the kinome, with particular focus on gene duplication events. We matched Arabidopsis proteins to a Hidden-Markov Model of eukaryotic kinases and computed a phylogeny of 942 Arabidopsis protein kinase domains and mapped their origin by gene duplication.

Results

The phylogeny showed two major clades of receptor kinases and soluble kinases, each of which was divided into functional subclades. Based on this phylogeny, association of yet uncharacterized kinases to families was possible which extended functional annotation of unknowns. Classification of gene duplications within these protein kinases revealed that representatives of cytosolic subfamilies showed a tendency to maintain segmentally duplicated genes, while some subfamilies of the receptor kinases were enriched for tandem duplicates. Although functional diversification is observed throughout most subfamilies, some instances of functional conservation among genes transposed from the same ancestor were observed. In general, a significant enrichment of essential genes was found among genes encoding for protein kinases.

Conclusions

The inferred phylogeny allowed classification and annotation of yet uncharacterized kinases. The prediction and analysis of syntenic blocks and duplication events within gene families of interest can be used to link functional biology to insights from an evolutionary viewpoint. The approach undertaken here can be applied to any gene family in any organism with an annotated genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-548) contains supplementary material, which is available to authorized users.     

Expulsion of Symbiotic Algae during Feeding by the Green Hydra – a Mechanism for Regulating Symbiont Density?

Background

Algal-cnidarian symbiosis is one of the main factors contributing to the success of cnidarians, and is crucial for the maintenance of coral reefs. While loss of the symbionts (such as in coral bleaching) may cause the death of the cnidarian host, over-proliferation of the algae may also harm the host. Thus, there is a need for the host to regulate the population density of its symbionts. In the green hydra, Chlorohydra viridissima, the density of symbiotic algae may be controlled through host modulation of the algal cell cycle. Alternatively, Chlorohydra may actively expel their endosymbionts, although this phenomenon has only been observed under experimentally contrived stress conditions.

Principal Findings

We show, using light and electron microscopy, that Chlorohydra actively expel endosymbiotic algal cells during predatory feeding on Artemia. This expulsion occurs as part of the apocrine mode of secretion from the endodermal digestive cells, but may also occur via an independent exocytotic mechanism.

Significance

Our results demonstrate, for the first time, active expulsion of endosymbiotic algae from cnidarians under natural conditions. We suggest this phenomenon may represent a mechanism whereby cnidarians can expel excess symbiotic algae when an alternative form of nutrition is available in the form of prey.     

Nitrogen stress response of a hybrid species: a gene expression study

Background and Aims

Low soil fertility limits growth and productivity in many natural and agricultural systems, where the ability to sense and respond to nutrient limitation is important for success. Helianthus anomalus is an annual sunflower of hybrid origin that is adapted to desert sand-dune substrates with lower fertility than its parental species, H. annuus and H. petiolaris. Previous studies have shown that H. anomalus has traits generally associated with adaptation to low-fertility habitats, including a lower inherent relative growth rate and longer leaf lifetime.

Methods

Here, a cDNA microarray is used to identify gene expression differences that potentially contribute to increased tolerance of low fertility of the hybrid species by comparing the nitrogen stress response of all three species with high- and low-nutrient treatments.

Key Results

Relative to the set of genes on the microarray, the genes showing differential expression in the hybrid species compared with its parents are enriched in stress-response genes, developmental genes, and genes involved in responses to biotic or abiotic stimuli. After a correction for multiple comparisons, five unique genes show a significantly different response to nitrogen limitation in H. anomalus compared with H. petiolaris and H. annuus. The Arabidopsis thaliana homologue of one of the five genes, catalase 1, has been shown to affect the timing of leaf senescence, and thus leaf lifespan.

Conclusions

The five genes identified in this analysis will be examined further as candidate genes for the adaptive stress response in H. anomalus. Genes that improve growth and productivity under nutrient stress could be used to improve crops for lower soil fertility which is common in marginal agricultural settings.