Decoupled evolution of coding region and mRNA expression patterns after gene duplication: implications for the neutralist-selectionist debate

The link between gene duplication and divergent patterns of gene expression across a complex life cycle

The diversification of many lineages throughout natural history has frequently been associated with evolutionary changes in life cycle complexity. However, our understanding of the processes that facilitate differentiation in the morphologies and functions expressed by organisms throughout their life cycles is limited. Theory suggests that the expression of traits is decoupled across life stages, thus allowing for their evolutionary independence. Although trait decoupling between stages is well established, explanations of how said decoupling evolves have seldom been considered. Because the different phenotypes expressed by organisms throughout their life cycles are coded for by the same genome, trait decoupling must be mediated through divergence in gene expression between stages. Gene duplication has been identified as an important mechanism that enables divergence in gene function and expression between cells and tissues. Because stage transitions across life cycles require changes in tissue types and functions, we investigated the potential link between gene duplication and expression divergence between life stages. To explore this idea, we examined the temporal changes in gene expression across the monarch butterfly (Danaus plexippus) metamorphosis. We found that within homologous groups, more phylogenetically diverged genes exhibited more distinct temporal expression patterns. This relationship scaled such that more phylogenetically diverse homologous groups showed more diverse patterns of gene expression. Furthermore, we found that duplicate genes showed increased stage-specificity relative to singleton genes. Overall, our findings suggest an important link between gene duplication and the evolution of complex life cycles.

Comparative genomics of sirenians reveals evolution of filaggrin and caspase-14 upon adaptation of the epidermis to aquatic life

The mammalian epidermis has evolved to protect the body in a dry environment. Genes of the epidermal differentiation complex (EDC), such as FLG (filaggrin), are implicated in the barrier function of the epidermis. Here, we investigated the molecular evolution of the EDC in sirenians (manatees and dugong), which have adapted to fully aquatic life, in comparison to the EDC of terrestrial mammals and aquatic mammals of the clade Cetacea (whales and dolphins). We show that the main subtypes of EDC genes are conserved or even duplicated, like late cornified envelope (LCE) genes of the dugong, whereas specific EDC genes have undergone inactivating mutations in sirenians. FLG contains premature stop codons in the dugong, and the ortholog of human CASP14 (caspase-14), which proteolytically processes filaggrin, is pseudogenized in the same species. As FLG and CASP14 have also been lost in whales, these mutations represent convergent evolution of skin barrier genes in different lineages of aquatic mammals. In contrast to the dugong, the manatee has retained functional FLG and CASP14 genes. FLG2 (filaggrin 2) is truncated in both species of sirenians investigated. We conclude that the land-to-water transition of sirenians was associated with modifications of the epidermal barrier at the molecular level.

Cryptic genetic variation shapes the fate of gene duplicates in a protein interaction network

Paralogous genes are often redundant for long periods of time before they diverge in function. While their functions are preserved, paralogous proteins can accumulate mutations that, through epistasis, could impact their fate in the future. By quantifying the impact of all single-amino acid substitutions on the binding of two myosin proteins to their interaction partners, we find that the future evolution of these proteins is highly contingent on their regulatory divergence and the mutations that have silently accumulated in their protein binding domains. Differences in the promoter strength of the two paralogs amplify the impact of mutations on binding in the lowly expressed one. While some mutations would be sufficient to non-functionalize one paralog, they would have minimal impact on the other. Our results reveal how functionally equivalent protein domains could be destined to specific fates by regulatory and cryptic coding sequence changes that currently have little to no functional impact.

Retention of duplicated genes in evolution

Gene duplication is a prevalent phenomenon across the tree of life. The processes that lead to the retention of duplicated genes are not well understood. Functional genomics approaches in model organisms, such as yeast, provide useful tools to test the mechanisms underlying retention with functional redundancy and divergence of duplicated genes, including fates associated with neofunctionalization, subfunctionalization, back-up compensation, and dosage amplification. Duplicated genes may also be retained as a consequence of structural and functional entanglement. Advances in human gene editing have enabled the interrogation of duplicated genes in the human genome, providing new tools to evaluate the relative contributions of each of these factors to duplicate gene retention and the evolution of genome structure.

Molecular evolution of bumble bee vitellogenin and vitellogenin-like genes

Vitellogenin (Vg), a storage protein, has been significantly studied for its egg yolk precursor role in oviparous animals. Recent studies found that vitellogenin and its Vg-like homologs were fundamentally involved in many other biological processes in social insects such as female caste differences and oxidative stress resilience. In this study, we conducted the first large-scale molecular evolutionary analyses of vitellogenin coding genes (Vg) and Vg-like genes of bumble bees, a primitively eusocial insect belonging to the genus Bombus. We obtained sequences for each of the four genes (Vg, Vg-like-A, Vg-like-B, and Vg-like-C) from 27 bumble bee genomes (nine were newly sequenced in this study), and sequences from the two closest clades of Bombus, including five Apis species and five Tetragonula species. Our molecular evolutionary analyses show that in bumble bee, the conventional Vg experienced strong positive selection, while the Vg-like genes showed overall relaxation of purifying selection. In Apis and Tetragonula; however, all four genes were found under purifying selection. Furthermore, the conventional Vg showed signs of strong positive selection in most subgenera in Bombus, apart from the obligate parasitic subgenus Psithyrus which has no caste differentiation. Together, these results indicate that the conventional Vg, a key pleiotropic gene in social insects, is the most rapidly evolving copy, potentially due to its multiple known social functions for both worker and queen castes. This study shows that concerted evolution and purifying selection shaped the evolution of the Vg gene family following their ancient gene duplication and may be the leading forces behind the evolution of new potential protein function enabling functional social pleiotropy.

Genome-Wide Identification of Flowering-Time Genes in Brassica Species and Reveals a Correlation between Selective Pressure and Expression Patterns of Vernalization-Pathway Genes in Brassica napus

Flowering time is a key agronomic trait, directly influencing crop yield and quality. Many flowering-time genes have been identified and characterized in the model plant Arabidopsis thaliana; however, these genes remain uncharacterized in many agronomically important Brassica crops. In this study, we identified 1064, 510, and 524 putative orthologs of A. thaliana flowering-time genes from Brassica napus, Brassica rapa, and Brassica oleracea, respectively, and found that genes involved in the aging and ambient temperature pathways were fewer than those in other flowering pathways. Flowering-time genes were distributed mostly on chromosome C03 in B. napus and B. oleracea, and on chromosome A09 in B. rapa. Calculation of non-synonymous (Ka)/synonymous substitution (Ks) ratios suggested that flowering-time genes in vernalization pathways experienced higher selection pressure than those in other pathways. Expression analysis showed that most vernalization-pathway genes were expressed in flowering organs. Approximately 40% of these genes were highly expressed in the anther, whereas flowering-time integrator genes were expressed in a highly organ-specific manner. Evolutionary selection pressures were negatively correlated with the breadth and expression levels of vernalization-pathway genes. These findings provide an integrated framework of flowering-time genes in these three Brassica crops and provide a foundation for deciphering the relationship between gene expression patterns and their evolutionary selection pressures in Brassica napus.

Population genomics of parallel evolution in gene expression and gene sequence during ecological adaptation

Natural selection often produces parallel phenotypic changes in response to a similar adaptive challenge. However, the extent to which parallel gene expression differences and genomic divergence underlie parallel phenotypic traits and whether they are decoupled or not remains largely unexplored. We performed a population genomic study of parallel ecological adaptation among replicate ecotype pairs of the rough periwinkle (Littorina saxatilis) at a regional geographical scale (NW Spain). We show that genomic changes underlying parallel phenotypic divergence followed a complex pattern of both repeatable differences and of differences unique to specific ecotype pairs, in which parallel changes in expression or sequence are restricted to a limited set of genes. Yet, the majority of divergent genes were divergent either for gene expression or coding sequence, but not for both simultaneously. Overall, our findings suggest that divergent selection significantly contributed to the process of parallel molecular differentiation among ecotype pairs, and that changes in expression and gene sequence underlying phenotypic divergence could, at least to a certain extent, be considered decoupled processes.

Identification and analysis of glutathione S-transferase gene family in sweet potato reveal divergent GST-mediated networks in aboveground and underground tissues in response to abiotic stresses

BACKGROUND

Sweet potato, a hexaploid species lacking a reference genome, is one of the most important crops in many developing countries, where abiotic stresses are a primary cause of reduction of crop yield. Glutathione S-transferases (GSTs) are multifunctional enzymes that play important roles in oxidative stress tolerance and cellular detoxification.

RESULTS

A total of 42 putative full-length GST genes were identified from two local transcriptome databases and validated by molecular cloning and Sanger sequencing. Sequence and intraspecific phylogenetic analyses revealed extensive differentiation in their coding sequences and divided them into eight subfamilies. Interspecific phylogenetic and comparative analyses indicated that most examined GST paralogs might originate and diverge before the speciation of sweet potato. Results from large-scale RNA-seq and quantitative real-time PCR experiments exhibited extensive variation in gene-expression profiles across different tissues and varieties, which implied strong evolutionary divergence in their gene-expression regulation. Moreover, we performed five manipulated stress experiments and uncovered highly divergent stress-response patterns of sweet potato GST genes in aboveground and underground tissues.

CONCLUSIONS

Our study identified a large number of sweet potato GST genes, systematically investigated their evolutionary diversification, and provides new insights into the GST-mediated stress-response mechanisms in this worldwide crop.

Underlying Principles of Natural Selection in Network Evolution: Systems Biology Approach

Systems biology is a rapidly expanding field that integrates diverse areas of science such as physics, engineering, computer science, mathematics, and biology toward the goal of elucidating the underlying principles of hierarchical metabolic and regulatory systems in the cell, and ultimately leading to predictive understanding of cellular response to perturbations. Because post-genomics research is taking place throughout the tree of life, comparative approaches offer a way for combining data from many organisms to shed light on the evolution and function of biological networks from the gene to the organismal level. Therefore, systems biology can build on decades of theoretical work in evolutionary biology, and at the same time evolutionary biology can use the systems biology approach to go in new uncharted directions. In this study, we present a review of how the post-genomics era is adopting comparative approaches and dynamic system methods to understand the underlying design principles of network evolution and to shape the nascent field of evolutionary systems biology. Finally, the application of evolutionary systems biology to robust biological network designs is also discussed from the synthetic biology perspective.

Genome-Wide Characterization and Expression Profiling of Sugar Transporter Family in the Whitefly, Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae)

Sugar transporters (STs) play pivotal roles in the growth, development, and stress responses of phloem-sucking insects, such as the whitefly, Bemisia tabaci. In this study, 137 sugar transporters (STs) were identified based on analysis of the genome and transcriptome of B. tabaci MEAM1. B. tabaci MEAM1 encodes a larger number of STs than other selected insects. Phylogenetic and molecular evolution analysis showed that the 137 STs formed three expanded clades and that the genes in Sternorrhyncha expanded clades had accelerated rates of evolution. B. tabaci sugar transporters (BTSTs) were divided into three groups based on their expression profiles across developmental stages; however, no host-specific BTST was found in B. tabaci fed on different host plants. Feeding of B. tabaci adults with feeding diet containing dsRNA significantly reduced the transcript level of the target genes in B. tabaci and mortality was significantly improved in B. tabaci fed on dsRNA compared to the control, which indicates the sugar transporters may be used as potential RNAi targets for B. tabaci bio-control. These results provide a foundation for further studies of STs in B. tabaci.

Bacilli glutamate dehydrogenases diverged via coevolution of transcription and enzyme regulation

The linkage between regulatory elements of transcription, such as promoters, and their protein products is central to gene function. Promoter-protein coevolution is therefore expected, but rarely observed, and the manner by which these two regulatory levels are linked remains largely unknown. We study glutamate dehydrogenase-a hub of carbon and nitrogen metabolism. In Bacillus subtilis, two paralogues exist: GudB is constitutively transcribed whereas RocG is tightly regulated. In their active, oligomeric states, both enzymes show similar enzymatic rates. However, swaps of enzymes and promoters cause severe fitness losses, thus indicating promoter-enzyme coevolution. Characterization of the proteins shows that, compared to RocG, GudB's enzymatic activity is highly dependent on glutamate and pH Promoter-enzyme swaps therefore result in excessive glutamate degradation when expressing a constitutive enzyme under a constitutive promoter, or insufficient activity when both the enzyme and its promoter are tightly regulated. Coevolution of transcriptional and enzymatic regulation therefore underlies paralogue-specific spatio-temporal control, especially under diverse growth conditions.

The Phenotypic Plasticity of Duplicated Genes in Saccharomyces cerevisiae and the Origin of Adaptations

Gene and genome duplication are the major sources of biological innovations in plants and animals. Functional and transcriptional divergence between the copies after gene duplication has been considered the main driver of innovations . However, here we show that increased phenotypic plasticity after duplication plays a more major role than thought before in the origin of adaptations. We perform an exhaustive analysis of the transcriptional alterations of duplicated genes in the unicellular eukaryote Saccharomyces cerevisiae when challenged with five different environmental stresses. Analysis of the transcriptomes of yeast shows that gene duplication increases the transcriptional response to environmental changes, with duplicated genes exhibiting signatures of adaptive transcriptional patterns in response to stress. The mechanism of duplication matters, with whole-genome duplicates being more transcriptionally altered than small-scale duplicates. The predominant transcriptional pattern follows the classic theory of evolution by gene duplication; with one gene copy remaining unaltered under stress, while its sister copy presents large transcriptional plasticity and a prominent role in adaptation. Moreover, we find additional transcriptional profiles that are suggestive of neo- and subfunctionalization of duplicate gene copies. These patterns are strongly correlated with the functional dependencies and sequence divergence profiles of gene copies. We show that, unlike singletons, duplicates respond more specifically to stress, supporting the role of natural selection in the transcriptional plasticity of duplicates. Our results reveal the underlying transcriptional complexity of duplicated genes and its role in the origin of adaptations.

Widespread and Adaptive Alterations in Genome-Wide Gene Expression Associated with Ecological Divergence of Two Oryza Species

Ecological speciation is a common mechanism by which new species arise. Despite great efforts, the role of gene expression in ecological divergence and speciation is poorly understood. Here, we conducted a genome-wide gene expression investigation of two Oryza species that are evolutionarily young and distinct in ecology and morphology. Using digital gene expression technology and the paired-end RNA sequencing method, we obtained 21,415 expressed genes across three reproduction-related tissues. Of them, approximately 8% (1,717) differed significantly in expression levels between the two species and these differentially expressed genes are randomly distributed across the genome. Moreover, 62% (1,064) of the differentially expressed genes exhibited a signature of directional selection in at least one species. Importantly, the genes with differential expression between species evolved more rapidly at the 5' flanking sequences than the genes without differential expression relative to coding sequences, suggesting that cis-regulatory changes are likely adaptive and play an important role in the ecological divergence of the two species. Finally, we showed evidence of significant differentiation between species in phenotype traits and observed that genes with differential expression were overrepresented with functional terms involving phenotypic and ecological differentiation between the two species, including reproduction- and stress-related characteristics. Our findings demonstrate that ecological speciation is associated with widespread and adaptive alterations in genome-wide gene expression and provide new insights into the importance of regulatory evolution in ecological speciation in plants.

Key Role of Amino Acid Repeat Expansions in the Functional Diversification of Duplicated Transcription Factors

The high regulatory complexity of vertebrates has been related to two rounds of whole genome duplication (2R-WGD) that occurred before the divergence of the major vertebrate groups. Following these events, many developmental transcription factors (TFs) were retained in multiple copies and subsequently specialized in diverse functions, whereas others reverted to their singleton state. TFs are known to be generally rich in amino acid repeats or low-complexity regions (LCRs), such as polyalanine or polyglutamine runs, which can evolve rapidly and potentially influence the transcriptional activity of the protein. Here we test the hypothesis that LCRs have played a major role in the diversification of TF gene duplicates. We find that nearly half of the TF gene families originated during the 2R-WGD contains LCRs. The number of gene duplicates with LCRs is 155 out of 550 analyzed (28%), about twice as many as the number of single copy genes with LCRs (15 out of 115, 13%). In addition, duplicated TFs preferentially accumulate certain LCR types, the most prominent of which are alanine repeats. We experimentally test the role of alanine-rich LCRs in two different TF gene families, PHOX2A/PHOX2B and LHX2/LHX9. In both cases, the presence of the alanine-rich LCR in one of the copies (PHOX2B and LHX2) significantly increases the capacity of the TF to activate transcription. Taken together, the results provide strong evidence that LCRs are important driving forces of evolutionary change in duplicated genes.

Maintenance and Loss of Duplicated Genes by Dosage Subfunctionalization

Whole-genome duplications (WGDs) have contributed to gene-repertoire enrichment in many eukaryotic lineages. However, most duplicated genes are eventually lost and it is still unclear why some duplicated genes are evolutionary successful whereas others quickly turn to pseudogenes. Here, we show that dosage constraints are major factors opposing post-WGD gene loss in several Paramecium species that share a common ancestral WGD. We propose a model where a majority of WGD-derived duplicates preserve their ancestral function and are retained to produce enough of the proteins performing this same ancestral function. Under this model, the expression level of individual duplicated genes can evolve neutrally as long as they maintain a roughly constant summed expression, and this allows random genetic drift toward uneven contributions of the two copies to total expression. Our analysis suggests that once a high level of imbalance is reached, which can require substantial lengths of time, the copy with the lowest expression level contributes a small enough fraction of the total expression that selection no longer opposes its loss. Extension of our analysis to yeast species sharing a common ancestral WGD yields similar results, suggesting that duplicated-gene retention for dosage constraints followed by divergence in expression level and eventual deterministic gene loss might be a universal feature of post-WGD evolution.

New insights into the evolutionary history of plant sorbitol dehydrogenase

BACKGROUND

Sorbitol dehydrogenase (SDH, EC 1.1.1.14) is the key enzyme involved in sorbitol metabolism in higher plants. SDH genes in some Rosaceae species could be divided into two groups. L-idonate-5-dehydrogenase (LIDH, EC 1.1.1.264) is involved in tartaric acid (TA) synthesis in Vitis vinifera and is highly homologous to plant SDHs. Despite efforts to understand the biological functions of plant SDH, the evolutionary history of plant SDH genes and their phylogenetic relationship with the V. vinifera LIDH gene have not been characterized.

RESULTS

A total of 92 SDH genes were identified from 42 angiosperm species. SDH genes have been highly duplicated within the Rosaceae family while monocot, Brassicaceae and most Asterid species exhibit singleton SDH genes. Core Eudicot SDHs have diverged into two phylogenetic lineages, now classified as SDH Class I and SDH Class II. V. vinifera LIDH was identified as a Class II SDH. Tandem duplication played a dominant role in the expansion of plant SDH family and Class II SDH genes were positioned in tandem with Class I SDH genes in several plant genomes. Protein modelling analyses of V. vinifera SDHs revealed 19 putative active site residues, three of which exhibited amino acid substitutions between Class I and Class II SDHs and were influenced by positive natural selection in the SDH Class II lineage. Gene expression analyses also demonstrated a clear transcriptional divergence between Class I and Class II SDH genes in V. vinifera and Citrus sinensis (orange).

CONCLUSIONS

Phylogenetic, natural selection and synteny analyses provided strong support for the emergence of SDH Class II by positive natural selection after tandem duplication in the common ancestor of core Eudicot plants. The substitutions of three putative active site residues might be responsible for the unique enzyme activity of V. vinifera LIDH, which belongs to SDH Class II and represents a novel function of SDH in V. vinifera that may be true also of other Class II SDHs. Gene expression analyses also supported the divergence of SDH Class II at the expression level. This study will facilitate future research into understanding the biological functions of plant SDHs.

Spatial features for Escherichia coli genome organization

Background

In bacterial genomes, the compactly encoded genes and operons are well organized, with genes in the same biological pathway or operons in the same regulon close to each other on the genome sequence. In addition, the linearly close genes have a higher probability of co-expression and their protein products tend to form protein–protein interactions. However, the organization features of bacterial genomes in a three-dimensional space remain elusive. The DNA interaction data of Escherichia coli, measured by the genome conformation capture (GCC) technique, have recently become available, which allowed us to investigate the spatial features of bacterial genome organization.

Results

By renormalizing the GCC data, we compared the interaction frequency of operon pairs in the same regulon with that of random operon pairs. The results showed that arrangements of operons in the E. coli genome tend to minimize the spatial distance between operons in the same regulon. A similar global organization feature exists for genes in biological pathways of E. coli. In addition, the genes close to each other spatially (even if they are far from each other on the genome sequence) tend to be co-expressed and form protein–protein interactions. These results provided new insights into the organization principles of bacterial genomes and support the notion of transcription factory.

Conclusions

This study revealed the organization features of Escherichia coli genomic functional units in the 3D space and furthered our understanding of the link between the three-dimensional structure of chromosomes and biological function.

Electronic supplementary material

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

Not only for egg yolk--functional and evolutionary insights from expression, selection, and structural analyses of Formica ant vitellogenins

Vitellogenin (Vg), a storage protein, has been extensively studied for its egg-yolk precursor role, and it has been suggested to be fundamentally involved in caste differences in social insects. More than one Vg copy has been reported in several oviparous species, including ants. However, the number and function of different Vgs, their phylogenetic relatedness, and their role in reproductive queens and nonreproductive workers have been studied in few species only. We studied caste-biased expression of Vgs in seven Formica ant species. Only one copy of conventional Vg was identified in Formica species, and three Vg homologs, derived from ancient duplications, which represent yet undiscovered Vg-like genes. We show that each of these Vg-like genes is present in all studied Hymenoptera and some of them in other insects as well. We show that after each major duplication event, at least one of the Vg-like genes has experienced a period of positive selection. This, combined with the observation that the Vg-like genes have acquired or lost specific protein domains suggests sub- or neofunctionalization between Vg and the duplicated genes. In contrast to earlier studies, Vg was not consistently queen biased in its expression, and the caste bias of the three Vg-like genes was highly variable among species. Furthermore, a truncated and Hymenoptera-specific Vg-like gene, Vg-like-C, was consistently worker biased. Multispecies comparisons are essential for Vg expression studies, and for gene expression studies in general, as we show that expression and also, putative functions cannot be generalized even among closely related species.

The evolutionary pattern and the regulation of stearoyl-CoA desaturase genes

Stearoyl-CoA desaturase (SCD) is a key enzyme that converts saturated fatty acids (SFAs) to monounsaturated fatty acids (MUFAs) in the biosynthesis of fat. To date, two isoforms of scd gene (scd1 and scd5) have been found widely existent in most of the vertebrate animals. However, the evolutionary patterns of both isofoms and the function of scd5 are poorly understandable. Herein, we aim to characterize the evolutionary pattern of scd genes and further predict the function differentiation of scd genes. The sequences of scd genes were highly conserved among eukaryote. Phylogenetic analysis identified two duplications of scd gene early in vertebrate evolution. The relative rate ratio test, branch-specific dN/dS ratio tests, and branch-site dN/dS ratio tests all suggested that the scd genes were evolved at a similar rate. The evolution of scd genes among eukaryote was under strictly purifying selection though several sites in scd1 and scd5 were undergone a relaxed selection pressure. The variable binding sites by transcriptional factors at the 5′-UTR and by miRNAs at 3′-UTR of scd genes suggested that the regulators of scd5 may be different from that of scd1. This study promotes our understanding of the evolutionary patterns and function of SCD genes in eukaryote.

Genome-wide analysis and evolutionary study of sucrose non-fermenting 1-related protein kinase 2 (SnRK2) gene family members in Arabidopsis and Oryza

The over-expression of plant specific SnRK2 gene family members by hyperosmotic stress and some by abscisic acid is well established. In this report, we have analyzed the evolution of SnRK2 gene family in different plant lineages including green algae, moss, lycophyte, dicot and monocot. Our results provide some evidences to indicate that the natural selection pressure had considerable influence on cis-regulatory promoter region and coding region of SnRK2 members in Arabidopsis and Oryza independently through time. Observed degree of sequence/motif conservation amongst SnRK2 homolog in all the analyzed plant lineages strongly supported their inclusion as members of this family. The chromosomal distributions of duplicated SnRK2 members have also been analyzed in Arabidopsis and Oryza. Massively Parallel Signature Sequencing (MPSS) database derived expression data and the presence of abiotic stress related promoter elements within the 1 kb upstream promoter region of these SnRK2 family members further strengthen the observations of previous workers. Additionally, the phylogenetic relationships of SnRK2 have been studied in all plant lineages along with their respective exon-intron structural patterns. Our results indicate that the ancestral SnRK2 gene of land plants gradually evolved by duplication and diversification and modified itself through exon-intron loss events to survive under environmental stress conditions.

In silico characterization and evolutionary analyses of CCAAT binding proteins in the lycophyte plant Selaginella moellendorffii genome: a growing comparative genomics resource

NF-Y transcription factors encoded by HAP gene family, composed of three subunits (HAP2/NF-YA, HAP3/NF-YB and HAP5/NF-YC), are capable of transcriptional regulation of target genes with high specificity by binding to the CCAAT-containing promoter sequences. Here, we have characterized duplicated HAP genes in Selaginella moellendorffii and explored some features that might be involved in the regulation of gene expression and their function. Subsequently, the evolutionary relationships of LEC1-type of HAP3 genes have been studied starting from lycophytes to angiosperm to reveal the details of conservation and diversification of these genes during plant evolution. Computational analyses demonstrated the variation in length of cis-regulatory region of HAP3 duplicates in S. moellendorffii containing three thermodynamically stable and evolutionarily conserved RNA secondary structures. The homology modeling of NF-Y proteins, secondary structural details, DNA binding large positive patches, binding affinity of H2A-H2B interactive residues of NF-YC subunits on the duplicated NF-YB subunits, conserved domain analyses and protein structural alignments indicated that gene duplication process of HAP genes in S. moellendorffii, followed by structural diversification, provide specific hints about their functional specificity under various circumstances for the survival of this lycophytic plant. We have identified several conserved motifs in LEC1 proteins among all plant lineages during evolution.

Accelerated evolution after gene duplication: a time-dependent process affecting just one copy

Gene duplication is widely regarded as a major mechanism modeling genome evolution and function. However, the mechanisms that drive the evolution of the two, initially redundant, gene copies are still ill defined. Many gene duplicates experience evolutionary rate acceleration, but the relative contribution of positive selection and random drift to the retention and subsequent evolution of gene duplicates, and for how long the molecular clock may be distorted by these processes, remains unclear. Focusing on rodent genes that duplicated before and after the mouse and rat split, we find significantly increased sequence divergence after duplication in only one of the copies, which in nearly all cases corresponds to the novel daughter copy, independent of the mechanism of duplication. We observe that the evolutionary rate of the accelerated copy, measured as the ratio of nonsynonymous to synonymous substitutions, is on average 5-fold higher in the period spanning 4-12 My after the duplication than it was before the duplication. This increase can be explained, at least in part, by the action of positive selection according to the results of the maximum likelihood-based branch-site test. Subsequently, the rate decelerates until purifying selection completely returns to preduplication levels. Reversion to the original rates has already been accomplished 40.5 My after the duplication event, corresponding to a genetic distance of about 0.28 synonymous substitutions per site. Differences in tissue gene expression patterns parallel those of substitution rates, reinforcing the role of neofunctionalization in explaining the evolution of young gene duplicates.

Expression divergence between Escherichia coli and Salmonella enterica serovar Typhimurium reflects their lifestyles

Escherichia coli K12 is a commensal bacteria and one of the best-studied model organisms. Salmonella enterica serovar Typhimurium, on the other hand, is a facultative intracellular pathogen. These two prokaryotic species can be considered related phylogenetically, and they share a large amount of their genetic material, which is commonly termed the "core genome." Despite their shared core genome, both species display very different lifestyles, and it is unclear to what extent the core genome, apart from the species-specific genes, plays a role in this lifestyle divergence. In this study, we focus on the differences in expression domains for the orthologous genes in E. coli and S. Typhimurium. The iterative comparison of coexpression methodology was used on large expression compendia of both species to uncover the conservation and divergence of gene expression. We found that gene expression conservation occurs mostly independently from amino acid similarity. According to our estimates, at least more than one quarter of the orthologous genes has a different expression domain in E. coli than in S. Typhimurium. Genes involved with key cellular processes are most likely to have conserved their expression domains, whereas genes showing diverged expression are associated with metabolic processes that, although present in both species, are regulated differently. The expression domains of the shared "core" genome of E. coli and S. Typhimurium, consisting of highly conserved orthologs, have been tuned to help accommodate the differences in lifestyle and the pathogenic potential of Salmonella.

Organ evolution in angiosperms driven by correlated divergences of gene sequences and expression patterns

The evolution of a species involves changes in its genome and its transcriptome. Divergence in expression patterns may be more important than divergence in sequences for determining phenotypic changes, particularly among closely related species. We examined the relationships between organ evolution, sequence evolution, and expression evolution in Arabidopsis thaliana, rice (Oryza sativa), and maize (Zea mays). We found correlated divergence of gene sequences and expression patterns, with distinct divergence rates that depend on the organ types in which a gene is expressed. For instance, genes specifically expressed in reproductive organs (i.e., stamen) evolve more quickly than those specifically expressed in vegetative organs (e.g., root). The different rates in organ evolution may be due to different degrees of functional constraint associated with the different physiological functions of plant organs. Additionally, the evolutionary rate of a gene sequence is correlated with the breadth of its expression in terms of the number of tissues, the number of coregulation modules, and the number of species in which the gene is expressed, as well as the number of genes with which it may interact. This linkage supports the hypothesis that constitutively expressed genes may experience higher levels of functional constraint accumulated from multiple tissues than do tissue-specific genes.

Social plasticity in fish: integrating mechanisms and function

Social plasticity is a ubiquitous feature of animal behaviour. Animals must adjust the expression of their social behaviour to the nuances of daily social life and to the transitions between life-history stages, and the ability to do so affects their Darwinian fitness. Here, an integrative framework is proposed for understanding the proximate mechanisms and ultimate consequences of social plasticity. According to this framework, social plasticity is achieved by rewiring or by biochemically switching nodes of the neural network underlying social behaviour in response to perceived social information. Therefore, at the molecular level, it depends on the social regulation of gene expression, so that different brain genomic and epigenetic states correspond to different behavioural responses and the switches between states are orchestrated by signalling pathways that interface the social environment and the genotype. At the evolutionary scale, social plasticity can be seen as an adaptive trait that can be under positive selection when changes in the environment outpace the rate of genetic evolutionary change. In cases when social plasticity is too costly or incomplete, behavioural consistency can emerge by directional selection that recruits gene modules corresponding to favoured behavioural states in that environment. As a result of this integrative approach, how knowledge of the proximate mechanisms underlying social plasticity is crucial to understanding its costs, limits and evolutionary consequences is shown, thereby highlighting the fact that proximate mechanisms contribute to the dynamics of selection. The role of teleosts as a premier model to study social plasticity is also highlighted, given the diversity and plasticity that this group exhibits in terms of social behaviour. Finally, the proposed integrative framework to social plasticity also illustrates how reciprocal causation analysis of biological phenomena (i.e. considering the interaction between proximate factors and evolutionary explanations) can be a more useful approach than the traditional proximate-ultimate dichotomy, according to which evolutionary processes can be understood without knowledge on proximate causes, thereby black-boxing developmental and physiological mechanisms.

Gene duplication in the evolution of sexual dimorphism

Males and females share most of the same genes, so selection in one sex will typically produce a correlated response in the other sex. Yet, the sexes have evolved to differ in a multitude of behavioral, morphological, and physiological traits. How did this sexual dimorphism evolve despite the presence of a common underlying genome? We investigated the potential role of gene duplication in the evolution of sexual dimorphism. Because duplication events provide extra genetic material, the sexes each might use this redundancy to facilitate sex-specific gene expression, permitting the evolution of dimorphism. We investigated this hypothesis at the genome-wide level in Drosophila melanogaster, using the presence of sex-biased expression as a proxy for the sex-specific specialization of gene function. We expected that if sexually antagonistic selection is a potent force acting upon individual genes, duplication will result in paralog families whose members differ in sex-biased expression. Gene members of the same duplicate family can have different expression patterns in males versus females. In particular, duplicate pairs containing a male-biased gene are found more frequently than expected, in agreement with previous studies. Furthermore, when the singleton ortholog is unbiased, duplication appears to allow one of the paralog copies to acquire male-biased expression. Conversely, female-biased expression is not common among duplicates; fewer duplicate genes are expressed in the female-soma and ovaries than in the male-soma and testes. Expression divergence exists more in older than in younger duplicates pairs, but expression divergence does not correlate with protein sequence divergence. Finally, genomic proximity may have an effect on whether paralogs differ in sex-biased expression. We conclude that the data are consistent with a role of gene duplication in fostering male-biased, but not female-biased, gene expression, thereby aiding the evolution of sexual dimorphism.

Genome-wide approaches to the study of adaptive gene expression evolution: systematic studies of evolutionary adaptations involving gene expression will allow many fundamental questions in evolutionary biology to be addressed

The role of gene expression in evolutionary adaptation has been a subject of debate for over 40 years. cis-regulation of transcription has been proposed to be the primary source of morphological novelty in evolution, though this is based on only a handful of examples. Recently the first genome-wide studies of gene expression adaptation have been published, giving us an initial global view of this process. Systematic studies such as these will allow a number of key questions currently facing the field of gene expression evolution to be addressed.

Evolutionary constraint in flanking regions of avian genes

An important comprehension from comparative genomic analysis is that sequence conservation beyond neutral expectations is frequently found outside protein-coding regions, indicating important functional roles of noncoding DNA. Understanding the causes of constraint on noncoding sequence evolution forms an important area of research, not least in light of the importance for understanding the evolution of gene expression. We aligned all orthologous genes of chicken and zebra finch together with 5 kb of their upstream and downstream noncoding sequences, to study the evolution of gene flanking sequences in the avian genome. Using ancestral repeats as a neutral reference, we detected significant evolutionary constraint in the 3' flanking region, highest directly after termination (60%) and then gradually decreasing to about 20% 5 kb downstream. Constraint was higher in annotated 3' untranslated regions (UTRs) than in non-UTRs at the same distance from the stop codon and higher in sequences annotated as microRNA (miRNA)-binding sites than in non-miRNA-binding sites within 3' UTRs. Constraint was also higher when estimated for a smaller data set of genes from more closely related songbird species, indicating turnover of functional elements during avian evolution. On the 5' flanking side constraint was readily seen within the first 125 bp immediately upstream of the start codon (34%) and was about 10% for remaining sequence within 5 kb upstream. Analysis of chicken polymorphism data gave further support for the highest constraint directly before and after the translated region. Finally, we found that genes evolving under the highest constraint measured by d(N)/d(S) also had the highest level of constraint in the 3' flanking region. This study broadens the insights into gene flanking sequence evolution by adding new findings from a vertebrate lineage other than mammals.

A single enhancer regulating the differential expression of duplicated red-sensitive opsin genes in zebrafish

A fundamental step in the evolution of the visual system is the gene duplication of visual opsins and differentiation between the duplicates in absorption spectra and expression pattern in the retina. However, our understanding of the mechanism of expression differentiation is far behind that of spectral tuning of opsins. Zebrafish (Danio rerio) have two red-sensitive cone opsin genes, LWS-1 and LWS-2. These genes are arrayed in a tail-to-head manner, in this order, and are both expressed in the long member of double cones (LDCs) in the retina. Expression of the longer-wave sensitive LWS-1 occurs later in development and is thus confined to the peripheral, especially ventral-nasal region of the adult retina, whereas expression of LWS-2 occurs earlier and is confined to the central region of the adult retina, shifted slightly to the dorsal-temporal region. In this study, we employed a transgenic reporter assay using fluorescent proteins and P1-artificial chromosome (PAC) clones encompassing the two genes and identified a 0.6-kb “LWS-activating region” (LAR) upstream of LWS-1, which regulates expression of both genes. Under the 2.6-kb flanking upstream region containing the LAR, the expression pattern of LWS-1 was recapitulated by the fluorescent reporter. On the other hand, when LAR was directly conjugated to the LWS-2 upstream region, the reporter was expressed in the LDCs but also across the entire outer nuclear layer. Deletion of LAR from the PAC clones drastically lowered the reporter expression of the two genes. These results suggest that LAR regulates both LWS-1 and LWS-2 by enhancing their expression and that interaction of LAR with the promoters is competitive between the two genes in a developmentally restricted manner. Sharing a regulatory region between duplicated genes could be a general way to facilitate the expression differentiation in duplicated visual opsins.

Among vertebrates, fish may have the most advanced color vision. They have greatly varied repertoires of color sensors called visual opsins, possibly reflecting evolutionary adaptation to their diverse photic environments in water, and are an excellent model to study the evolution of vertebrate color vision. This is achieved by multiplying opsin genes and differentiating their absorption light spectra and expression patterns. However, little is understood regarding how the opsin genes are regulated to achieve the differential expression pattern. In this study, we focused on the duplicated red-sensitive opsin genes of zebrafish to tackle this problem. We discovered an “enhancer” region near the two red opsin genes that plays a crucial role in their differential expression pattern. Our results suggest that the two red opsin genes interact with the enhancer competitively in a developmentally restricted manner. Sharing a regulatory region could be a general way to facilitate the expression differentiation in duplicated visual opsin genes.

Gene expression in a paleopolyploid: a transcriptome resource for the ciliate Paramecium tetraurelia

Background

The genome of Paramecium tetraurelia, a unicellular model that belongs to the ciliate phylum, has been shaped by at least 3 successive whole genome duplications (WGD). These dramatic events, which have also been documented in plants, animals and fungi, are resolved over evolutionary time by the loss of one duplicate for the majority of genes. Thanks to a low rate of large scale genome rearrangement in Paramecium, an unprecedented large number of gene duplicates of different ages have been identified, making this organism an outstanding model to investigate the evolutionary consequences of polyploidization. The most recent WGD, with 51% of pre-duplication genes still in 2 copies, provides a snapshot of a phase of rapid gene loss that is not accessible in more ancient polyploids such as yeast.

Results

We designed a custom oligonucleotide microarray platform for P. tetraurelia genome-wide expression profiling and used the platform to measure gene expression during 1) the sexual cycle of autogamy, 2) growth of new cilia in response to deciliation and 3) biogenesis of secretory granules after massive exocytosis. Genes that are differentially expressed during these time course experiments have expression patterns consistent with a very low rate of subfunctionalization (partition of ancestral functions between duplicated genes) in particular since the most recent polyploidization event.

Conclusions

A public transcriptome resource is now available for Paramecium tetraurelia. The resource has been integrated into the ParameciumDB model organism database, providing searchable access to the data. The microarray platform, freely available through NimbleGen Systems, provides a robust, cost-effective approach for genome-wide expression profiling in P. tetraurelia. The expression data support previous studies showing that at short evolutionary times after a whole genome duplication, gene dosage balance constraints and not functional change are the major determinants of gene retention.

Comparison of computational models for assessing conservation of gene expression across species

Assessing conservation/divergence of gene expression across species is important for the understanding of gene regulation evolution. Although advances in microarray technology have provided massive high-dimensional gene expression data, the analysis of such data is still challenging. To date, assessing cross-species conservation of gene expression using microarray data has been mainly based on comparison of expression patterns across corresponding tissues, or comparison of co-expression of a gene with a reference set of genes. Because direct and reliable high-throughput experimental data on conservation of gene expression are often unavailable, the assessment of these two computational models is very challenging and has not been reported yet. In this study, we compared one corresponding tissue based method and three co-expression based methods for assessing conservation of gene expression, in terms of their pair-wise agreements, using a frequently used human-mouse tissue expression dataset. We find that 1) the co-expression based methods are only moderately correlated with the corresponding tissue based methods, 2) the reliability of co-expression based methods is affected by the size of the reference ortholog set, and 3) the corresponding tissue based methods may lose some information for assessing conservation of gene expression. We suggest that the use of either of these two computational models to study the evolution of a gene's expression may be subject to great uncertainty, and the investigation of changes in both gene expression patterns over corresponding tissues and co-expression of the gene with other genes is necessary.

Expression divergence of tandemly arrayed genes in human and mouse

Tandemly arrayed genes (TAGs) account for about one third of the duplicated genes in eukaryotic genomes, yet there has not been any systematic study of their gene expression patterns. Taking advantage of recently published large-scale microarray data sets, we studied the expression divergence of 361 two-member TAGs in human and 212 two-member TAGs in mouse and examined the effect of sequence divergence, gene orientation, and chromosomal proximity on the divergence of TAG expression patterns. Our results show that there is a weak negative correlation between sequence divergence of TAG members and their expression similarity. There is also a weak negative correlation between chromosomal proximity of TAG members and their expression similarity. We did not detect any significant relationship between gene orientation and expression similarity. We also found that downstream TAG members do not show significantly narrower expression breadth than upstream members, contrary to what we predict based on TAG expression divergence hypothesis that we propose. Finally, we show that both chromosomal proximity and expression correlation in TAGs do not differ significantly from their neighboring non-TAG gene pairs, suggesting that tandem duplication is unlikely to be the cause for the higher-than-random expression association between neighboring genes on a chromosome in human and mouse.

The developmental expression dynamics of Drosophila melanogaster transcription factors

BACKGROUND

Site-specific transcription factors (TFs) are coordinators of developmental and physiological gene expression programs. Their binding to cis-regulatory modules of target genes mediates the precise cell- and context-specific activation and repression of genes. The expression of TFs should therefore reflect the core expression program of each cell.

RESULTS

We studied the expression dynamics of about 750 TFs using the available genomics resources in Drosophila melanogaster. We find that 95% of these TFs are expressed at some point during embryonic development, with a peak roughly between 10 and 12 hours after egg laying, the core stages of organogenesis. We address the differential utilization of DNA-binding domains in different developmental programs systematically in a spatio-temporal context, and show that the zinc finger class of TFs is predominantly early expressed, while Homeobox TFs exhibit later expression in embryogenesis.

CONCLUSIONS

Previous work, dissecting cis-regulatory modules during Drosophila development, suggests that TFs are deployed in groups acting in a cooperative manner. In contrast, we find that there is rapid exchange of co-expressed partners amongst the fly TFs, at rates similar to the genome-wide dynamics of co-expression clusters. This suggests there may also be a high level of combinatorial complexity of TFs at cis-regulatory modules.

Duplication mechanism and disruptions in flanking regions determine the fate of Mammalian gene duplicates

Here we identify duplicated genes in five mammalian genomes and classify these duplicates based on the mechanisms by which they were generated. Retrotransposition accounts for at least half of all predicted duplicate genes in these genomes, with tandem and interspersed DNA-mediated duplicates comprising the other half. Estimation of the evolutionary rates in each class revealed greater rate asymmetry between retrotransposed and interspersed DNA duplicate pairs than between tandem duplicates, suggesting that retrotransposed and interspersed DNA duplicates are diverging more quickly. In an attempt to understand the basis of this asymmetry, we identified disruption of flanking DNA as an indicator of new duplicate fate-loss of local synteny accelerates the asymmetry of divergence of interspersed DNA duplicates. We also show that intact retrogenes are enriched in intergenic regions and indel purified regions of the human genome. Moreover, intact retrogenes closest to annotated genes show the greatest levels of purifying selective pressure. Together, these findings suggest that the differential evolution of duplicate genes may be significantly influenced by changes in local genome architecture.

From evidence to inference: probing the evolution of protein interaction networks

The evolutionary mechanisms by which protein interaction networks grow and change are beginning to be appreciated as a major factor shaping their present-day structures and properties. Starting with a consideration of the biases and errors inherent in our current views of these networks, we discuss the dangers of constructing evolutionary arguments from naïve analyses of network topology. We argue that progress in understanding the processes of network evolution is only possible when hypotheses are formulated as plausible evolutionary models and compared against the observed data within the framework of probabilistic modeling. The value of such models is expected to be greatly enhanced as they incorporate more of the details of the biophysical properties of interacting proteins, gene phylogeny, and measurement error and as more advanced methodologies emerge for model comparison and the inference of ancestral network states.

Variation in gene duplicates with low synonymous divergence in Saccharomyces cerevisiae relative to Caenorhabditis elegans

Differences between yeast and worm duplicates result from differences in mechanisms of duplication and effective population size.

Background

The direct examination of large, unbiased samples of young gene duplicates in their early stages of evolution is crucial to understanding the origin, divergence and preservation of new genes. Furthermore, comparative analysis of multiple genomes is necessary to determine whether patterns of gene duplication can be generalized across diverse lineages or are species-specific. Here we present results from an analysis comprising 68 duplication events in the Saccharomyces cerevisiae genome. We partition the yeast duplicates into ohnologs (generated by a whole-genome duplication) and non-ohnologs (from small-scale duplication events) to determine whether their disparate origins commit them to divergent evolutionary trajectories and genomic attributes.

Results

We conclude that, for the most part, ohnologs tend to appear remarkably similar to non-ohnologs in their structural attributes (specifically the relative composition frequencies of complete, partial and chimeric duplicates), the discernible length of the duplicated region (duplication span) as well as genomic location. Furthermore, we find notable differences in the features of S. cerevisiae gene duplicates relative to those of another eukaryote, Caenorhabditis elegans, with respect to chromosomal location, extent of duplication and the relative frequencies of complete, partial and chimeric duplications.

Conclusions

We conclude that the variation between yeast and worm duplicates can be attributed to differing mechanisms of duplication in conjunction with the varying efficacy of natural selection in these two genomes as dictated by their disparate effective population sizes.

A comprehensive analysis of gene expression evolution between humans and mice

Evolutionary changes in gene expression account for most phenotypic differences between species. Advances in microarray technology have made the systematic study of gene expression evolution possible. In this study, gene expression patterns were compared between human and mouse genomes using two published methods. Specifically, we studied how gene expression evolution was related to GO terms and tried to decode the relationship between promoter evolution and gene expression evolution. The results showed that (1) the significant enrichment of biological processes in orthologs of expression conservation reveals functional significance of gene expression conservation. The more conserved gene expression in some biological processes than is expected in a purely neutral model reveals negative selection on gene expression. However, fast evolving genes mainly support the neutrality of gene expression evolution, and (2) gene expression conservation is positively but only slightly correlated with promoter conservation based on a motif-count score of the promoter alignment. Our results suggest a neutral model with negative selection for gene expression evolution between humans and mice, and promoter evolution could have some effects on gene expression evolution.

Comparative transcriptomics in the Triticeae

Background

Barley and particularly wheat are two grass species of immense agricultural importance. In spite of polyploidization events within the latter, studies have shown that genotypically and phenotypically these species are very closely related and, indeed, fertile hybrids can be created by interbreeding. The advent of two genome-scale Affymetrix GeneChips now allows studies of the comparison of their transcriptomes.

Results

We have used the Wheat GeneChip to create a "gene expression atlas" for the wheat transcriptome (cv. Chinese Spring). For this, we chose mRNA from a range of tissues and developmental stages closely mirroring a comparable study carried out for barley (cv. Morex) using the Barley1 GeneChip. This, together with large-scale clustering of the probesets from the two GeneChips into "homologous groups", has allowed us to perform a genomic-scale comparative study of expression patterns in these two species. We explore the influence of the polyploidy of wheat on the results obtained with the Wheat GeneChip and quantify the correlation between conservation in gene sequence and gene expression in wheat and barley. In addition, we show how the conservation of expression patterns can be used to elucidate, probeset by probeset, the reliability of the Wheat GeneChip.

Conclusion

While there are many differences in expression on the level of individual genes and tissues, we demonstrate that the wheat and barley transcriptomes appear highly correlated. This finding is significant not only because given small evolutionary distance between the two species it is widely expected, but also because it demonstrates that it is possible to use the two GeneChips for comparative studies. This is the case even though their probeset composition reflects rather different design principles as well as, of course, the present incomplete knowledge of the gene content of the two species. We also show that, in general, the Wheat GeneChip is not able to distinguish contributions from individual homoeologs. Furthermore, the comparison between the two species leads us to conclude that the conservation of both gene sequence as well as gene expression is positively correlated with absolute expression levels, presumably reflecting increased selection pressure on genes coding for proteins present at high levels. In addition, the results indicate the presence of a correlation between sequence and expression conservation within the Triticeae.

Coding region structural heterogeneity and turnover of transcription start sites contribute to divergence in expression between duplicate genes

BACKGROUND

Gene expression divergence is one manifestation of functional differences between duplicate genes. Although rapid accumulation of expression divergence between duplicate gene copies has been observed, the driving mechanisms behind this phenomenon have not been explored in detail.

RESULTS

We examine which factors influence expression divergence between human duplicate genes, utilizing the latest genome-wide data sets. We conclude that the turnover of transcription start sites between duplicate genes occurs rapidly after gene duplication and that gene pairs with shared transcription start sites have significantly higher expression similarity than those without shared transcription start sites. Moreover, we find that most (55%) duplicate gene pairs do not retain the same coding sequence structure between the two duplicate copies and this also contributes to divergence in their expression. Furthermore, the proportion of aligned sequences in cis-regulatory regions between the two copies is positively correlated with expression similarity. Surprisingly, we find no effect of copy-specific transposable element insertions on the divergence of duplicate gene expression.

CONCLUSIONS

Our results suggest that turnover of transcription start sites, structural heterogeneity of coding sequences, and divergence of cis-regulatory regions between copies play a pivotal role in determining the expression divergence of duplicate genes.

Diverse transcriptional programs associated with environmental stress and hormones in the Arabidopsis receptor-like kinase gene family

The genome of Arabidopsis thaliana encodes more than 600 receptor-like kinase (RLK) genes, by far the dominant class of receptors found in land plants. Although similar to the mammalian receptor tyrosine kinases, plant RLKs are serine/threonine kinases that represent a novel signaling innovation unique to plants and, consequently, an excellent opportunity to understand how extracellular signaling evolved and functions in plants as opposed to animals. RLKs are predicted to be major components of the signaling pathways that allow plants to respond to environmental and developmental conditions. However, breakthroughs in identifying these processes have been limited to only a handful of individual RLKs. Here, we used a Syngenta custom Arabidopsis GeneChip array to compile a detailed profile of the transcriptional activity of 604 receptor-like kinase genes after exposure to a cross-section of known signaling factors in plants, including abiotic stresses, biotic stresses, and hormones. In the 68 experiments comprising the study, we found that 582 of the 604 RLK genes displayed a two-fold or greater change in expression to at least one of 12 types of treatments, thereby providing a large body of experimental evidence for targeted functional screens of individual RLK genes. We investigated whether particular subfamilies of RLK genes are responsive to specific types of signals and found that each subfamily displayed broad ranges of expression, as opposed to being targeted towards particular signal classes. Finally, by analyzing the divergence of sequence and gene expression among the RLK subfamilies, we present evidence as to the functional basis for the expansion of the RLKs and how this expansion may have affected conservation and divergences in their function. Taken as a whole, our study represents a preliminary, working model of processes and interactions in which the members of the RLK gene family may be involved, where such information has remained elusive for so many of its members.

Assessing the evolution of gene expression using microarray data

Classical studies of the evolution of gene function have predominantly focused on mutations within protein coding regions. With the advent of microarrays, however, it has become possible to evaluate the transcriptional activity of a gene as an additional characteristic of function. Recent studies have revealed an equally important role for gene regulation in the retention and evolution of duplicate genes. Here we review approaches to assessing the evolution of gene expression using microarray data, and discuss potential influences on expression divergence. Currently, there are no established standards on how best to identify and quantify instances of expression divergence. There have also been few efforts to date that incorporate suspected influences into mathematical models of expression divergence. Such developments will be crucial to a comprehensive understanding of the role gene duplications and expression evolution play in the emergence of complex traits and functional diversity. An integrative approach to gene family evolution, including both orthologous and paralogous genes, has the potential to bring strong predictive power both to the functional annotation of extant proteins and to the inference of functional characteristics of ancestral gene family members.

Weak correlation between sequence conservation in promoter regions and in protein-coding regions of human-mouse orthologous gene pairs

Background

Interspecies sequence comparison is a powerful tool to extract functional or evolutionary information from the genomes of organisms. A number of studies have compared protein sequences or promoter sequences between mammals, which provided many insights into genomics. However, the correlation between protein conservation and promoter conservation remains controversial.

Results

We examined promoter conservation as well as protein conservation for 6,901 human and mouse orthologous genes, and observed a very weak correlation between them. We further investigated their relationship by decomposing it based on functional categories, and identified categories with significant tendencies. Remarkably, the 'ribosome' category showed significantly low promoter conservation, despite its high protein conservation, and the 'extracellular matrix' category showed significantly high promoter conservation, in spite of its low protein conservation.

Conclusion

Our results show the relation of gene function to protein conservation and promoter conservation, and revealed that there seem to be nonparallel components between protein and promoter sequence evolution.

Evolution of gene sequence and gene expression are not correlated in yeast

We show that, in yeast, the divergence rate of gene expression is not correlated with that of its associated coding sequence. Gene essentiality influences both modes of evolution, but other properties related to protein structure or promoter composition are only correlated with coding-sequence divergence or gene expression divergence, respectively. Based on these findings, we discuss the possibilities of neutral evolution of gene expression and of different modes of evolution in unicellular versus multicellular organisms.

EPGD: a comprehensive web resource for integrating and displaying eukaryotic paralog/paralogon information

Gene duplication is common in all three domains of life, especially in eukaryotic genomes. The duplicates provide new material for the action of evolutionary forces such as selection or genetic drift. Here we describe a sophisticated procedure to extract duplicated genes (paralogs) from 26 available eukaryotic genomes, to pre-calculate several evolutionary indexes (evolutionary rate, synonymous distance/clock, transition redundant exchange clock, etc.) based on the paralog family, and to identify block or segmental duplications (paralogons). We also constructed an internet-accessible Eukaryotic Paralog Group Database (EPGD; http://epgd.biosino.org/EPGD/). The database is gene-centered and organized by paralog family. It focuses on paralogs and evolutionary duplication events. The paralog families and paralogons can be searched by text or sequence, and are downloadable from the website as plain text files. The database will be very useful for both experimentalists and bioinformaticians interested in the study of duplication events or paralog families.

Comparative analysis indicates regulatory neofunctionalization of yeast duplicates

Comparison of the expression profiles of S. cerevisiae duplicate pairs with that of their pre-duplication orthologs in C. albicans identified a class of genes that may present cases of regulatory neofunctionalization.

Background

Gene duplication provides raw material for the generation of new functions, but most duplicates are rapidly lost due to the initial redundancy in gene function. How gene function diversifies following duplication is largely unclear. Previous studies analyzed the diversification of duplicates by characterizing their coding sequence divergence. However, functional divergence can also be attributed to changes in regulatory properties, such as protein localization or expression, which require only minor changes in gene sequence.

Results

We developed a novel method to compare expression profiles from different organisms and applied it to analyze the expression divergence of yeast duplicated genes. The expression profiles of Saccharomyces cerevisiae duplicate pairs were compared with those of their pre-duplication orthologs in Candida albicans. Duplicate pairs were classified into two classes, corresponding to symmetric versus asymmetric rates of expression divergence. The latter class includes 43 duplicate pairs in which only one copy has a significant expression similarity to the C. albicans ortholog. These may present cases of regulatory neofunctionalization, as supported also by their dispensability and variability.

Conclusion

Duplicated genes may diversify through regulatory neofunctionalization. Notably, the asymmetry of gene sequence evolution and the asymmetry of gene expression evolution are only weakly correlated, underscoring the importance of expression analysis to elucidate the evolution of novel functions.

Underlying principles of natural selection in network evolution: systems biology approach

Systems biology is a rapidly expanding field that integrates diverse areas of science such as physics, engineering, computer science, mathematics, and biology toward the goal of elucidating the underlying principles of hierarchical metabolic and regulatory systems in the cell, and ultimately leading to predictive understanding of cellular response to perturbations. Because post-genomics research is taking place throughout the tree of life, comparative approaches offer a way for combining data from many organisms to shed light on the evolution and function of biological networks from the gene to the organismal level. Therefore, systems biology can build on decades of theoretical work in evolutionary biology, and at the same time evolutionary biology can use the systems biology approach to go in new uncharted directions. In this study, we present a review of how the post-genomics era is adopting comparative approaches and dynamic system methods to understand the underlying design principles of network evolution and to shape the nascent field of evolutionary systems biology. Finally, the application of evolutionary systems biology to robust biological network designs is also discussed from the synthetic biology perspective.