The strength and pattern of natural selection on gene expression in rice.

Levels of gene expression underpin organismal phenotypes1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown1,2. Here we assayed gene expression in rice (Oryza sativa)3, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic6 (0.04%), and transcripts that display lower levels of expression and stochastic noise7-9 and higher levels of plasticity9 are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity9. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes4,5, but that the efficacy of selection is genetically constrained under drought conditions10. Drought selected for earlier flowering11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering13-marking this gene as a drought-escape gene11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

Evolutionary systems biology reveals patterns of rice adaptation to drought-prone agro-ecosystems.

Rice (Oryza sativa) was domesticated around 10,000 years ago and has developed into a staple for half of humanity. The crop evolved and is currently grown in stably wet and intermittently dry agro-ecosystems, but patterns of adaptation to differences in water availability remain poorly understood. While previous field studies have evaluated plant developmental adaptations to water deficit, adaptive variation in functional and hydraulic components, particularly in relation to gene expression, has received less attention. Here, we take an evolutionary systems biology approach to characterize adaptive drought resistance traits across roots and shoots. We find that rice harbors heritable variation in molecular, physiological, and morphological traits that is linked to higher fitness under drought. We identify modules of co-expressed genes that are associated with adaptive drought avoidance and tolerance mechanisms. These expression modules showed evidence of polygenic adaptation in rice subgroups harboring accessions that evolved in drought-prone agro-ecosystems. Fitness-linked expression patterns allowed us to identify the drought-adaptive nature of optimizing photosynthesis and interactions with arbuscular mycorrhizal fungi. Taken together, our study provides an unprecedented, integrative view of rice adaptation to water-limited field conditions.

Rapid evolutionary changes in gene expression in response to climate fluctuations.

There is now abundant evidence of rapid evolution in natural populations, but the genetic mechanisms of these changes remain unclear. One possible route to rapid evolution is through changes in the expression of genes that influence traits under selection. We examined contemporary evolutionary gene expression changes in plant populations responding to environmental fluctuations. We compared genome-wide gene expression, using RNA-seq, in two populations of Brassica rapa collected over four time points between 1997 and 2014, during which precipitation in southern California fluctuated dramatically and phenotypic and genotypic changes occurred. By combining transcriptome profiling with the resurrection approach, we directly examined evolutionary changes in gene expression over time. For both populations, we found a substantial number of differentially expressed genes between generations, indicating rapid evolution in the expression of many genes. Using existing gene annotations, we found that many changes occurred in genes involved in regulating stress responses and flowering time. These appeared related to the fluctuations in precipitation and were potentially adaptive. However, the evolutionary changes in gene expression differed across generations within and between populations, indicating largely independent evolutionary trajectories across populations and over time. Our study provides strong evidence for rapid evolution in gene expression, and indicates that changes in gene expression can be one mechanism of rapid evolutionary responses to selection episodes. This study also illustrates that combining resurrection studies with transcriptomics is a powerful approach for investigating evolutionary changes at the gene regulatory level, and will provide new insights into the genetic basis of contemporary evolution.

MT1-MMP Cooperates with TGF-ß Receptor-Mediated Signaling to Trigger SNAIL and Induce Epithelial-to-Mesenchymal-like Transition in U87 Glioblastoma Cells.

Epithelial-to-mesenchymal transition (EMT) recapitulates metastasis and can be induced in vitro through transforming growth factor (TGF)-ß signaling. A role for MMP activity in glioblastoma multiforme has been ascribed to EMT, but the molecular crosstalk between TGF-ß signaling and membrane type 1 MMP (MT1-MMP) remains poorly understood. Here, the expression of common EMT biomarkers, induced through TGF-ß and the MT1-MMP inducer concanavalin A (ConA), was explored using RNA-seq analysis and differential gene arrays in human U87 glioblastoma cells. TGF-ß triggered SNAIL and fibronectin expressions in 2D-adherent and 3D-spheroid U87 glioblastoma cell models. Those inductions were antagonized by the TGF-ß receptor kinase inhibitor galunisertib, the JAK/STAT inhibitors AG490 and tofacitinib, and by the diet-derived epigallocatechin gallate (EGCG). Transient gene silencing of MT1-MMP prevented the induction of SNAIL by ConA and abrogated TGF-ß-induced cell chemotaxis. Moreover, ConA induced STAT3 and Src phosphorylation, suggesting these pathways to be involved in the MT1-MMP-mediated signaling axis that led to SNAIL induction. Our findings highlight a new signaling axis linking MT1-MMP to TGF-ß-mediated EMT-like induction in glioblastoma cells, the process of which can be prevented by the diet-derived EGCG.

Phylogenetic and Genomic Analyses Resolve the Origin of Important Plant Genes Derived from Transposable Elements.

Once perceived as merely selfish, transposable elements (TEs) are now recognized as potent agents of adaptation. One way TEs contribute to evolution is through TE exaptation, a process whereby TEs, which persist by replicating in the genome, transform into novel host genes, which persist by conferring phenotypic benefits. Known exapted TEs (ETEs) contribute diverse and vital functions, and may facilitate punctuated equilibrium, yet little is known about this process. To better understand TE exaptation, we designed an approach to resolve the phylogenetic context and timing of exaptation events and subsequent patterns of ETE diversification. Starting with known ETEs, we search in diverse genomes for basal ETEs and closely related TEs, carefully curate the numerous candidate sequences, and infer detailed phylogenies. To distinguish TEs from ETEs, we also weigh several key genomic characteristics including repetitiveness, terminal repeats, pseudogenic features, and conserved domains. Applying this approach to the well-characterized plant ETEs MUG and FHY3, we show that each group is paraphyletic and we argue that this pattern demonstrates that each originated in not one but multiple exaptation events. These exaptations and subsequent ETE diversification occurred throughout angiosperm evolution including the crown group expansion, the angiosperm radiation, and the primitive evolution of angiosperms. In addition, we detect evidence of several putative novel ETE families. Our findings support the hypothesis that TE exaptation generates novel genes more frequently than is currently thought, often coinciding with key periods of evolution.

Diversity and evolution of transposable elements in Arabidopsis.

Transposable elements are mobile genetic elements that have successfully populated eukaryotic genomes and show diversity in their structure and transposition mechanisms. Although first viewed solely as selfish, transposable elements are now known as important vectors to drive the adaptation and evolution of their host genome. Transposable elements can affect host gene structures, gene copy number, gene expression, and even as a source for novel genes. For example, a number of transposable element sequences have been co-opted to contribute to evolutionary innovation, such as the mammalian placenta and the vertebrate immune system. In plants, the need to adapt rapidly to changing environmental conditions is essential and is reflected, as will be discussed, by genome plasticity and an abundance of diverse, active transposon families. This review focuses on transposable elements in plants, particularly those that have beneficial effects on the host. We also emphasize the importance of having proper tools to annotate and classify transposons to better understand their biology.

A gene family derived from transposable elements during early angiosperm evolution has reproductive fitness benefits in Arabidopsis thaliana.

The benefits of ever-growing numbers of sequenced eukaryotic genomes will not be fully realized until we learn to decipher vast stretches of noncoding DNA, largely composed of transposable elements. Transposable elements persist through self-replication, but some genes once encoded by transposable elements have, through a process called molecular domestication, evolved new functions that increase fitness. Although they have conferred numerous adaptations, the number of such domesticated transposable element genes remains unknown, so their evolutionary and functional impact cannot be fully assessed. Systematic searches that exploit genomic signatures of natural selection have been employed to identify potential domesticated genes, but their predictions have yet to be experimentally verified. To this end, we investigated a family of domesticated genes called MUSTANG (MUG), identified in a previous bioinformatic search of plant genomes. We show that MUG genes are functional. Mutants of Arabidopsis thaliana MUG genes yield phenotypes with severely reduced plant fitness through decreased plant size, delayed flowering, abnormal development of floral organs, and markedly reduced fertility. MUG genes are present in all flowering plants, but not in any non-flowering plant lineages, such as gymnosperms, suggesting that the molecular domestication of MUG may have been an integral part of early angiosperm evolution. This study shows that systematic searches can be successful at identifying functional genetic elements in noncoding regions and demonstrates how to combine systematic searches with reverse genetics in a fruitful way to decipher eukaryotic genomes.

Systems genomics of salinity stress response in rice.

Populations can adapt to stressful environments through changes in gene expression. However, the role of gene regulation in mediating stress response and adaptation remains largely unexplored. Here, we use an integrative field dataset obtained from 780 plants of Oryza sativa ssp. indica (rice) grown in a field experiment under normal or moderate salt stress conditions to examine selection and evolution of gene expression variation under salinity stress conditions. We find that salinity stress induces increased selective pressure on gene expression. Further, we show that trans-eQTLs rather than cis-eQTLs are primarily associated with rice's gene expression under salinity stress, potentially via a few master-regulators. Importantly, and contrary to the expectations, we find that cis-trans reinforcement is more common than cis-trans compensation which may be reflective of rice diversification subsequent to domestication. We further identify genetic fixation as the likely mechanism underlying this compensation/reinforcement. Additionally, we show that cis- and trans-eQTLs are under different selection regimes, giving us insights into the evolutionary dynamics of gene expression variation. By examining genomic, transcriptomic, and phenotypic variation across a rice population, we gain insights into the molecular and genetic landscape underlying adaptive salinity stress responses, which is relevant for other crops and other stresses.

Refolding of the integral membrane protein light-harvesting complex II monitored by pulse EPR.

The major light-harvesting chlorophyll a/b complex (LHCII) of the photosynthetic apparatus in plants self-organizes in vitro. The recombinant apoprotein, denatured in dodecyl sulfate, spontaneously folds when it is mixed with its pigments, chlorophylls, and carotenoids in detergent solution, and assembles into structurally authentic LHCII in the course of several minutes. Pulse EPR techniques, specifically double-electron-electron resonance (DEER), have been used to analyze protein folding during this process. Pairs of nitroxide labels were introduced site-specifically into recombinant LHCII and shown not to affect the stability and function of the pigment-protein complex. Interspin distance distributions between two spin pairs were measured at various time points, one pair located on either end of the second transmembrane helix (helix 3), the other one located near the luminal ends of the intertwined transmembrane helices 1 and 4. In the dodecyl sulfate-solubilized apoprotein, both distance distributions were consistent with a random-coil protein structure. A rapid freeze-quench experiment on the latter spin pair indicated that 1 s after initiating reconstitution the protein structure is virtually unchanged. Subsequently, both distance distributions monitored protein folding in the same time range in which the assembly of chlorophylls into the complex had been observed. The positioning of the spin pair spanning the hydrophobic core of LHCII clearly preceded the juxtaposition of the spin pair on the luminal side of the complex. This indicates that superhelix formation of helices 1 and 4 is a late step in LHCII assembly.

The influence of genetic architecture on responses to selection under drought in rice.

Accurately predicting responses to selection is a major goal in biology and important for successful crop breeding in changing environments. However, evolutionary responses to selection can be constrained by such factors as genetic and cross-environment correlations, linkage, and pleiotropy, and our understanding of the extent and impact of such constraints is still developing. Here, we conducted a field experiment to investigate potential constraints to selection for drought resistance in rice (Oryza sativa) using phenotypic selection analysis and quantitative genetics. We found that traits related to drought response were heritable, and some were under selection, including selection for earlier flowering, which could allow drought escape. However, patterns of selection generally were not opposite under wet and dry conditions, and we did not find individual or closely linked genes that influenced multiple traits, indicating a lack of evidence that antagonistic pleiotropy, linkage, or cross-environment correlations would constrain selection for drought resistance. In most cases, genetic correlations had little influence on responses to selection, with direct and indirect selection largely congruent. The exception to this was seed mass under drought, which was predicted to evolve in the opposite direction of direct selection due to correlations. Because of this indirect effect on selection on seed mass, selection for drought resistance was not accompanied by a decrease in seed mass, and yield increased with fecundity. Furthermore, breeding lines with high fitness and yield under drought also had high fitness and yield under wet conditions, indicating that there was no evidence for a yield penalty on drought resistance. We found multiple genes in which expression influenced both water use efficiency (WUE) and days to first flowering, supporting a genetic basis for the trade-off between drought escape and avoidance strategies. Together, these results can provide helpful guidance for understanding and managing evolutionary constraints and breeding stress-resistant crops.

An inferred fitness consequence map of the rice genome.

The extent to which sequence variation impacts plant fitness is poorly understood. High-resolution maps detailing the constraint acting on the genome, especially in regulatory sites, would be beneficial as functional annotation of noncoding sequences remains sparse. Here, we present a fitness consequence (fitCons) map for rice (Oryza sativa). We inferred fitCons scores (ρ) for 246 inferred genome classes derived from nine functional genomic and epigenomic datasets, including chromatin accessibility, messenger RNA/small RNA transcription, DNA methylation, histone modifications and engaged RNA polymerase activity. These were integrated with genome-wide polymorphism and divergence data from 1,477 rice accessions and 11 reference genome sequences in the Oryzeae. We found ρ to be multimodal, with ~9% of the rice genome falling into classes where more than half of the bases would probably have a fitness consequence if mutated. Around 2% of the rice genome showed evidence of weak negative selection, frequently at candidate regulatory sites, including a novel set of 1,000 potentially active enhancer elements. This fitCons map provides perspective on the evolutionary forces associated with genome diversity, aids in genome annotation and can guide crop breeding programs.

Exaptation of transposable element coding sequences.

Transposable elements (TEs) are mobile genetic elements that were once perceived as merely selfish, but are now recognized as potent agents of adaptation. One way TEs contribute to genome evolution is through TE exaptation, a process whereby TEs, which usually persist by replicating in the genome, transform into novel host genes, which thereafter persist by conferring phenotypic benefits. Exapted TEs are known to contribute diverse and vital functions, and may facilitate punctuated equilibrium, yet we have little understanding about the process of TE exaptation. In order to facilitate our understanding of how TE coding sequences may become exapted, here we incorporate the findings of recent publications into a framework and six-step model.

Abiotic Stress Phenotypes Are Associated with Conserved Genes Derived from Transposable Elements.

Plant phenomics offers unique opportunities to accelerate our understanding of gene function and plant response to different environments, and may be particularly useful for studying previously uncharacterized genes. One important type of poorly characterized genes is those derived from transposable elements (TEs), which have departed from a mobility-driven lifestyle to attain new adaptive roles for the host (exapted TEs). We used phenomics approaches, coupled with reverse genetics, to analyze T-DNA insertion mutants of both previously reported and novel protein-coding exapted TEs in the model plant Arabidopsis thaliana. We show that mutations in most of these exapted TEs result in phenotypes, particularly when challenged by abiotic stress. We built statistical multi-dimensional phenotypic profiles and compared them to wild-type and known stress responsive mutant lines for each particular stress condition. We found that these exapted TEs may play roles in responses to phosphate limitation, tolerance to high salt concentration, freezing temperatures, and arsenic toxicity. These results not only experimentally validate a large set of putative functional exapted TEs recently discovered through computational analysis, but also uncover additional novel phenotypes for previously well-characterized exapted TEs in A. thaliana.

Developing maps of fitness consequences for plant genomes.

Predicting the fitness consequences of mutations, and their concomitant impacts on molecular and cellular function as well as organismal phenotypes, is an important challenge in biology that has new relevance in an era when genomic data is readily available. The ability to construct genomewide maps of fitness consequences in plant genomes is a recent development that has profound implications for our ability to predict the fitness effects of mutations and discover functional elements. Here we highlight approaches to building fitness consequence maps to infer regions under selection. We emphasize computational methods applied primarily to the study of human disease that translate physical maps of within-species genome variation into maps of fitness effects of individual natural mutations. Maps of fitness consequences in plants, combined with traditional genetic approaches, could accelerate discovery of functional elements such as regulatory sequences in non-coding DNA and genetic polymorphisms associated with key traits, including agronomically-important traits such as yield and environmental stress responses.

Reproducing on time when temperature varies: shifts in the timing of courtship by fiddler crabs.

Many species reproduce when conditions are most favorable for the survival of young. Numerous intertidal fish and invertebrates release eggs or larvae during semilunar, large amplitude, nocturnal tides when these early life stages are best able to escape predation by fish that feed near the shore during the day. Remarkably, some species, including the fiddler crabs Uca terpsichores and Uca deichmanni, maintain this timing throughout the year as temperature, and thus the rate of embryonic development, vary. The mechanisms that allow such precision in the timing of the production of young are poorly known. A preliminary study suggested that when temperature decreases, U. terpsichores mate earlier in the tidal amplitude cycle such that larvae are released at the appropriate time. We tested this idea by studying the timing of courtship in U. terpsichores and U. deichmanni as temperature varied annually during two years, at 5 locations that differed in the temperature of the sediment where females incubate their eggs. Uca terpsichores courted earlier at locations where sediment temperature declined seasonally but not where sediment temperature remained elevated throughout the year. In contrast, clear shifts in courtship timing were not observed for U. deichmanni despite variation in sediment temperature. We discuss other mechanisms by which this species may maintain reproductive timing. These two species are likely to be affected differently by changes in the frequency and intensity of cold periods that are expected to accompany climate change.

An atlas of over 90,000 conserved noncoding sequences provides insight into crucifer regulatory regions.

Despite the central importance of noncoding DNA to gene regulation and evolution, understanding of the extent of selection on plant noncoding DNA remains limited compared to that of other organisms. Here we report sequencing of genomes from three Brassicaceae species (Leavenworthia alabamica, Sisymbrium irio and Aethionema arabicum) and their joint analysis with six previously sequenced crucifer genomes. Conservation across orthologous bases suggests that at least 17% of the Arabidopsis thaliana genome is under selection, with nearly one-quarter of the sequence under selection lying outside of coding regions. Much of this sequence can be localized to approximately 90,000 conserved noncoding sequences (CNSs) that show evidence of transcriptional and post-transcriptional regulation. Population genomics analyses of two crucifer species, A. thaliana and Capsella grandiflora, confirm that most of the identified CNSs are evolving under medium to strong purifying selection. Overall, these CNSs highlight both similarities and several key differences between the regulatory DNA of plants and other species.