Variation within and among species in gene expression: raw material for evolution

Transcriptomic response to nitrogen availability reveals signatures of adaptive plasticity during tetraploid wheat domestication

The domestication of crops, coupled with agroecosystem development, is associated with major environmental changes and provides an ideal model of phenotypic plasticity. Here, we examined 32 genotypes of three tetraploid wheat (Triticum turgidum L.) subspecies, wild emmer, emmer, and durum wheat, which are representative of the key stages in the domestication of tetraploid wheat. We developed a pipeline that integrates RNA-Seq data and population genomics to assess gene expression plasticity and identify selection signatures under diverse nitrogen availability conditions. Our analysis revealed differing gene expression responses to nitrogen availability across primary (wild emmer to emmer) and secondary (emmer to durum wheat) domestication. Notably, nitrogen triggered the expression of twice as many genes in durum wheat compared to that in emmer and wild emmer. Unique selection signatures were identified at each stage: primary domestication mainly influenced genes related to biotic interactions, whereas secondary domestication affected genes related to amino acid metabolism, in particular lysine. Selection signatures were found in differentially expressed genes (DEGs), notably those associated with nitrogen metabolism, such as the gene encoding glutamate dehydrogenase (GDH). Overall, our study highlights the pivotal role of nitrogen availability in the domestication and adaptive responses of a major food crop, with varying effects across different traits and growth conditions.

Tyrosyl-DNA phosphodiesterase 2 (Tdp2) repairs DNA-protein crosslinks and protects against double strand breaks in vivo

DNA-protein crosslinks pose a significant challenge to genome stability and cell viability. Efficient repair of DPCs is crucial for preserving genomic integrity and preventing the accumulation of DNA damage. Despite recent advances in our understanding of DPC repair, many aspects of this process, especially at the organismal level, remain elusive. In this study, we used zebrafish as a model organism to investigate the role of TDP2 (Tyrosyl-DNA phosphodiesterase 2) in DPC repair. We characterized the two tdp2 orthologs in zebrafish using phylogenetic, syntenic and expression analysis and investigated the phenotypic consequences of tdp2 silencing in zebrafish embryos. We then quantified the effects of tdp2a and tdp2b silencing on cellular DPC levels and DSB accumulation in zebrafish embryos. Our findings revealed that tdp2b is the main ortholog during embryonic development, while both orthologs are ubiquitously present in adult tissues. Notably, the tdp2b ortholog is phylogenetically closer to human TDP2. Silencing of tdp2b, but not tdp2a, resulted in the loss of Tdp2 activity in zebrafish embryos, accompanied by the accumulation of DPCs and DSBs. Our findings contribute to a more comprehensive understanding of DPC repair at the organismal level and underscore the significance of TDP2 in maintaining genome stability.

Expression of Survivin, CK7, ASH1, HMGB3, L587S, and CLCA2 in Peripheral Blood of Lung Cancer Patients by Real-Time Polymerase Chain Reaction

Introduction The objective of the present study was to identify gene expression in peripheral blood by a real-time polymerase chain reaction (PCR) technique in patients who have lung carcinoma. Material and methods Peripheral blood samples of patients with non-small cell and small cell lung cancer were collected. Target genes included survivin, CK7, ASH1, HMGB3, L587S, and CLCA2. β-Actin was the reference gene. If the mean CT (threshold cycle) value for a target gene is ≥40, the gene expression is considered undetectable. Results Fifty patients with lung carcinoma were included and 30 healthy controls. Out of the six genes, survivin showed 26.8 times fold change as compared to controls; ASH1 and L587S were 0.54 and 0.06, respectively; and HMGB3, CLCA2, and CK7 had non-significant fold change in comparison to controls. The overall detection rate of the six target genes examined in lung cancer was 84%, with 42 out of 50 patients testing positive. Higher stages and ASH1 (p = 0.031), CK7 (p = <0.001), and HMGB3, p = 0.011 were associated significantly. CLCA2 had higher expression in patients without adrenal metastases (p = 0.044). Conclusions Lifestyle and geographical variation might be a probable cause of variable gene expression as compared to other studies. However, further research is needed to determine the clinical implication of these markers, especially in larger groups of early-stage patients.

Transcriptomics of differences in thermal plasticity associated with selection for an exaggerated male sexual trait

The information about the magnitude of differences in thermal plasticity both between and within populations, as well as identification of the underlying molecular mechanisms are key to understanding the evolution of thermal plasticity. In particular, genes underlying variation in the physiological response to temperature can provide raw material for selection acting on plastic traits. Using RNAseq, we investigate the transcriptional response to temperature in males and females from bulb mite populations selected for the increased frequency of one of two discrete male morphs (fighter- and scrambler-selected populations) that differ in relative fitness depending on temperature. We show that different mechanisms underlie the divergence in thermal response between fighter- and scrambler-selected populations at decreased vs. increased temperature. Temperature decrease to 18 °C was associated with higher transcriptomic plasticity of males with more elaborate armaments, as indicated by a significant selection-by-temperature interaction effect on the expression of 40 genes, 38 of which were upregulated in fighter-selected populations in response to temperature decrease. In response to 28 °C, no selection-by-temperature interaction in gene expression was detected. Hence, differences in phenotypic response to temperature increase likely depended on genes associated with their distinct morph-specific thermal tolerance. Selection of males also drove gene expression patterns in females. These patterns could be associated with temperature-dependent fitness differences between females from fighter- vs. scrambler-selected populations reported in previous studies. Our study shows that selection for divergent male sexually selected morphologies and behaviors has a potential to drive divergence in metabolic pathways underlying plastic response to temperature in both sexes.

Synergistic response to climate stressors in coral is associated with genotypic variation in baseline expression

As environments are rapidly reshaped due to climate change, phenotypic plasticity plays an important role in the ability of organisms to persist and is considered an especially important acclimatization mechanism for long-lived sessile organisms such as reef-building corals. Often, this ability of a single genotype to display multiple phenotypes depending on the environment is modulated by changes in gene expression, which can vary in response to environmental changes via two mechanisms: baseline expression and expression plasticity. We used transcriptome-wide expression profiling of eleven genotypes of common-gardened Acropora cervicornis to explore genotypic variation in the expression response to thermal and acidification stress, both individually and in combination. We show that the combination of these two stressors elicits a synergistic gene expression response, and that both baseline expression and expression plasticity in response to stress show genotypic variation. Additionally, we demonstrate that frontloading of a large module of coexpressed genes is associated with greater retention of algal symbionts under combined stress. These results illustrate that variation in the gene expression response of individuals to climate change stressors can persist even when individuals have shared environmental histories, affecting their performance under future climate change scenarios.

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Sexual conflict, heterochrony and tissue specificity as evolutionary problems of adaptive plasticity in development

Differential gene expression represents a fundamental cause and manifestation of phenotypic plasticity. Adaptive phenotypic plasticity in gene expression as a trait evolves when alleles that mediate gene regulation serve to increase organismal fitness by improving the alignment of variation in gene expression with variation in circumstances. Among the diverse circumstances that a gene encounters are distinct cell types, developmental stages and sexes, as well as an organism's extrinsic ecological environments. Consequently, adaptive phenotypic plasticity provides a common framework to consider diverse evolutionary problems by considering the shared implications of alleles that produce context-dependent gene expression. From this perspective, adaptive plasticity represents an evolutionary resolution to conflicts of interest that arise from any negatively pleiotropic effects of expression of a gene across ontogeny, among tissues, between the sexes, or across extrinsic environments. This view highlights shared properties within the general relation of fitness, trait expression and context that may nonetheless differ substantively in the grain of selection within and among generations to influence the likelihood of adaptive plasticity as an evolutionary response. Research programmes that historically have focused on these separate issues may use the insights from one another by recognizing their shared dependence on context-dependent gene regulatory evolution.

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii: Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST-QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis-environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

Population transcriptogenomics highlights impaired metabolism and small population sizes in tree frogs living in the Chernobyl Exclusion Zone

BACKGROUND

Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant.

RESULTS

We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes.

CONCLUSIONS

Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Stochastic Modeling of Gene Expression Evolution Uncovers Tissue- and Sex-Specific Properties of Expression Evolution in the Drosophila Genus

Gene expression evolution is typically modeled with the stochastic Ornstein-Uhlenbeck (OU) process. It has been suggested that the estimation of within-species variations using replicated data can increase the predictive power of such models, but this hypothesis has not been fully tested. We developed EvoGeneX, a computationally efficient implementation of the OU-based method that models within-species variation. Using extensive simulations, we show that modeling within-species variations and appropriate selection of species improve the performance of the model. Further, to facilitate a comparative analysis of expression evolution, we introduce a formal measure of evolutionary expression divergence for a group of genes using the rate and the asymptotic level of divergence. With these tools in hand, we performed the first-ever analysis of the evolution of gene expression across different body-parts, species, and sexes of the Drosophila genus. We observed that genes with adaptive expression evolution tend to be body-part specific, whereas the genes with constrained evolution tend to be shared across body-parts. Among the neutrally evolving gene expression patterns, gonads in both sexes have higher expression divergence relative to other tissues and the male gonads have even higher divergence than the female gonads. Among the evolutionarily constrained genes, the gonads show different divergence patterns, where the male gonads, and not the female gonads, show less constrained divergence than other body-parts. Finally, we show interesting examples of adaptive expression evolution, including adaptation of odor binding proteins.

Placing human gene families into their evolutionary context

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

Novel-miR-310 mediated response mechanism to Cry1Ac protoxin in Plutella xylostella (L.)

The diamondback moth (DBM), Plutella xylostella (L.), has evolved resistance to multiple insecticides including Bacillus thuringiensis (Bt). ATP-binding cassette (ABC) transporters are a class of transmembrane protein families, involved in multiple physiological processes and pesticide resistances in insects. However, the role and regulatory mechanism of ABC transporter in mediating the response to Bt Cry1Ac toxin remain unclear. Here, we characterized a MAPK signaling pathway-enriched ABCG subfamily gene PxABCG20 from DBM, and found it was differentially expressed in the Cry1Ac-resistant and Cry1Ac-susceptible strains. RNAi knockdown of PxABCG20 increased the tolerance of DBM to Cry1Ac protoxin. To explore the regulatory mechanism of PxABCG20 expression, we predicted the potential miRNAs targeting PxABCG20 using two target prediction algorithms. Luciferase reporter assay confirmed that novel-miR-310 was able to down-regulate PxABCG20 expression in HEK293T cells. Furthermore, injection of novel-miR-310 agomir markedly inhibited PxABCG20 expression, resulting in increased tolerance to Cry1Ac protoxin in susceptible strain, while injection of novel-miR-310 antagomir markedly induced the expression of PxABCG20, leading to decreased tolerance to Cry1Ac protoxin. Our work provides theoretical basis for exploring novel targets for the DBM response to Cry1Ac toxin and expands the understanding of miRNA role in mediating the susceptibility of insect pest to Cry1Ac toxin.

Environment-driven shifts in inter-individual variation and phenotypic integration within subnetworks of the mussel transcriptome and proteome

The environment can alter the magnitude of phenotypic variation among individuals, potentially influencing evolutionary trajectories. However, environmental influences on variation are complex and remain understudied. Populations in heterogeneous environments might exhibit more variation, the amount of variation could differ between benign and stressful conditions, and/or variation might manifest in different ways among stages of the gene‐to‐protein expression cascade or among physiological functions. Here, we explore these three issues by quantifying patterns of inter‐individual variation in both transcript and protein expression levels among California mussels, Mytilus californianus Conrad. Mussels were exposed to five ecologically relevant treatments that varied in the mean and interindividual heterogeneity of body temperature. To target a diverse set of physiological functions, we assessed variation within 19 expression subnetworks, including canonical stress‐response pathways and empirically derived coexpression clusters that represent a diffuse set of cellular processes. Variation in expression was particularly pronounced in the treatments with high mean and heterogeneous body temperatures. However, with few exceptions, environment‐dependent shifts of variation in the transcriptome were not reflected in the proteome. A metric of phenotypic integration provided evidence for a greater degree of constraint on relative expression levels (i.e., stronger correlation) within expression subnetworks in benign, homogeneous environments. Our results suggest that environments that are more stressful on average – and which also tend to be more heterogeneous – can relax these expression constraints and reduce phenotypic integration within biochemical subnetworks. Context‐dependent “unmasking” of functional variation may contribute to interindividual differences in physiological phenotype and performance in stressful environments.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the allele frequency spectrum of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a superabundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of coexpressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

Long-lasting effects of chronic exposure to chemical pollution on the hologenome of the Manila clam

Chronic exposure to pollutants affects natural populations, creating specific molecular and biochemical signatures. In the present study, we tested the hypothesis that chronic exposure to pollutants might have substantial effects on the Manila clam hologenome long after removal from contaminated sites. To reach this goal, a highly integrative approach was implemented, combining transcriptome, genetic and microbiota analyses with the evaluation of biochemical and histological profiles of the edible Manila clam Ruditapes philippinarum, as it was transplanted for 6 months from the polluted area of Porto Marghera (PM) to the clean area of Chioggia (Venice lagoon, Italy). One month post-transplantation, PM clams showed several modifications to its resident microbiota, including an overrepresentation of the opportunistic pathogen Arcobacter spp. This may be related to the upregulation of several immune genes in the PM clams, potentially representing a host response to the increased abundance of deleterious bacteria. Six months after transplantation, PM clams demonstrated a lower ability to respond to environmental/physiological stressors related to the summer season, and the hepatopancreas-associated microbiota still showed different compositions among PM and CH clams. This study confirms that different stressors have predictable effects in clams at different biological levels and demonstrates that chronic exposure to pollutants leads to long-lasting effects on the animal hologenome. In addition, no genetic differentiation between samples from the two areas was detected, confirming that PM and CH clams belong to a single population. Overall, the obtained responses were largely reversible and potentially related to phenotypic plasticity rather than genetic adaptation. The results here presented will be functional for the assessment of the environmental risk imposed by chemicals on an economically important bivalve species.

Gene expression responses to thermal shifts in the endangered lichen Lobaria pulmonaria

Anthropogenic climate change has led to unprecedented shifts in temperature across many ecosystems. In a context of rapid environmental changes, acclimation is an important process as it may influence the capacity of organisms to survive under novel thermal conditions. Mechanisms of acclimation could involve upregulation of stress response genes involved in protein folding, DNA damage repair and the regulation of signal transduction genes, along with a simultaneous downregulation of genes involved in growth or the cell cycle, in order to maintain cellular functions and equilibria. We transplanted Lobaria pulmonaria lichens originating from different forests to determine the relative effects of long-term acclimation and genetic factors on the variability in expression of mycobiont and photobiont genes. We found a strong response of the mycobiont and photobiont to high temperatures, regardless of sample origin. The green-algal photobiont had an overall lower response than the mycobiont. Gene expression of both symbionts was also influenced by acclimation to transplantation sites and by genetic factors. L. pulmonaria seems to have evolved powerful molecular pathways to deal with environmental fluctuations and stress and can acclimate to new habitats by transcriptomic convergence. Although L. pulmonaria has the molecular machinery to counteract short-term thermal stress, survival of lichens such as L. pulmonaria depends mostly on their long-term positive carbon balance, which can be compromised by higher temperatures and reduced precipitation, and both these outcomes have been predicted for Central Europe in connection with global climate change.

A cis-Acting Mutation in the PxABCG1 Promoter Is Associated with Cry1Ac Resistance in Plutella xylostella (L.)

The molecular mechanisms of insect resistance to Cry toxins generated from the bacterium Bacillus thuringiensis (Bt) urgently need to be elucidated to enable the improvement and sustainability of Bt-based products. Although downregulation of the expression of midgut receptor genes is a pivotal mechanism of insect resistance to Bt Cry toxins, the underlying transcriptional regulation of these genes remains elusive. Herein, we unraveled the regulatory mechanism of the downregulation of the ABC transporter gene PxABCG1 (also called Pxwhite), a functional midgut receptor of the Bt Cry1Ac toxin in Plutella xylostella. The PxABCG1 promoters of Cry1Ac-susceptible and Cry1Ac-resistant strains were cloned and analyzed, and they showed clear differences in activity. Subsequently, a dual-luciferase reporter assay, a yeast one-hybrid (Y1H) assay, and RNA interference (RNAi) experiments demonstrated that a cis-mutation in a binding site of the Hox transcription factor Antennapedia (Antp) decreased the promoter activity of the resistant strain and eliminated the binding and regulation of Antp, thereby enhancing the resistance of P. xylostella to the Cry1Ac toxin. These results advance our knowledge of the roles of cis- and trans-regulatory variations in the regulation of midgut Cry receptor genes and the evolution of Bt resistance, contributing to a more complete understanding of the Bt resistance mechanism.

Variation of microRNA expression in the human placenta driven by population identity and sex of the newborn

BACKGROUND

Analysis of lymphocyte cell lines revealed substantial differences in the expression of mRNA and microRNA (miRNA) among human populations. The extent of such population-associated differences in actual human tissues remains largely unexplored. The placenta is one of the few solid human tissues that can be collected in substantial numbers in a controlled manner, enabling quantitative analysis of transient biomolecules such as RNA transcripts. Here, we analyzed microRNA (miRNA) expression in human placental samples derived from 36 individuals representing four genetically distinct human populations: African Americans, European Americans, South Asians, and East Asians. All samples were collected at the same hospital following a unified protocol, thus minimizing potential biases that might influence the results.

RESULTS

Sequence analysis of the miRNA fraction yielded 938 annotated and 70 novel miRNA transcripts expressed in the placenta. Of them, 82 (9%) of annotated and 11 (16%) of novel miRNAs displayed quantitative expression differences among populations, generally reflecting reported genetic and mRNA-expression-based distances. Several co-expressed miRNA clusters stood out from the rest of the population-associated differences in terms of miRNA evolutionary age, tissue-specificity, and disease-association characteristics. Among three non-environmental influenced demographic parameters, the second largest contributor to miRNA expression variation after population was the sex of the newborn, with 32 miRNAs (3% of detected) exhibiting significant expression differences depending on whether the newborn was male or female. Male-associated miRNAs were evolutionarily younger and correlated inversely with the expression of target mRNA involved in neuron-related functions. In contrast, both male and female-associated miRNAs appeared to mediate different types of hormonal responses. Demographic factors further affected reported imprinted expression of 66 placental miRNAs: the imprinting strength correlated with the mother's weight, but not height.

CONCLUSIONS

Our results showed that among 12 assessed demographic variables, population affiliation and fetal sex had a substantial influence on miRNA expression variation among human placental samples. The effect of newborn-sex-associated miRNA differences further led to expression inhibition of the target genes clustering in specific functional pathways. By contrast, population-driven miRNA differences might mainly represent neutral changes with minimal functional impacts.

Magnitude and Predictability of pH Fluctuations Shape Plastic Responses to Ocean Acidification

AbstractPhenotypic plasticity is expected to facilitate the persistence of natural populations as global change progresses. The attributes of fluctuating environments that favor the evolution of plasticity have received extensive theoretical investigation, yet empirical validation of these findings is still in its infancy. Here, we combine high-resolution environmental data with a laboratory-based experiment to explore the influence of habitat pH fluctuation dynamics on the plasticity of gene expression in two populations of the Mediterranean mussel, Mytilus galloprovincialis. We linked differences in the magnitude and predictability of pH fluctuations in two habitats to population-specific gene expression profiles in ambient and stressful pH treatments. Our results demonstrate population-based differentiation in gene expression plasticity, whereby mussels native to a habitat exhibiting a large magnitude of pH fluctuations with low predictability display reduced phenotypic plasticity between experimentally imposed pH treatments. This work validates recent theoretical findings on evolution in fluctuating environments, suggesting that the predictability of fluctuating selection pressures may play a predominant role in shaping the phenotypic variation observed across natural populations.

A framework for understanding gene expression plasticity and its influence on stress tolerance

Phenotypic plasticity can serve as a stepping stone towards adaptation. Recently, studies have shown that gene expression contributes to emergent stress responses such as thermal tolerance, with tolerant and susceptible populations showing distinct transcriptional profiles. However, given the dynamic nature of gene expression, interpreting transcriptomic results in a way that elucidates the functional connection between gene expression and the observed stress response is challenging. Here, we present a conceptual framework to guide interpretation of gene expression reaction norms in the context of stress tolerance. We consider the evolutionary and adaptive potential of gene expression reaction norms and discuss the influence of sampling timing, transcriptomic resilience, as well as complexities related to life history when interpreting gene expression dynamics and how these patterns relate to host tolerance. We highlight corals as a case study to demonstrate the value of this framework for non-model systems. As species face rapidly changing environmental conditions, modulating gene expression can serve as a mechanistic link from genetic and cellular processes to the physiological responses that allow organisms to thrive under novel conditions. Interpreting how or whether a species can employ gene expression plasticity to ensure short-term survival will be critical for understanding the global impacts of climate change across diverse taxa.

Signals of selection beyond bottlenecks between exotic populations of the bull-headed dung beetle, Onthophagus taurus

Colonization of new environments can lead to population bottlenecks and rapid phenotypic evolution that could be due to neutral and selective processes. Exotic populations of the bull-headed dung beetle (Onthophagus taurus) have differentiated in opposite directions from native beetles in male horn-to-body size allometry and female fecundity. Here we test for genetic and transcriptional differences among two exotic and one native O. taurus populations after three generations in common garden conditions. We sequenced RNA from 24 individuals for each of the three populations including both sexes, and spanning four developmental stages for the two exotic, differentiated populations. Identifying 270,400 high-quality single nucleotide polymorphisms, we revealed a strong signal of genetic differentiation between the three populations, and evidence of recent bottlenecks within and an excess of outlier loci between exotic populations. Differences in gene expression between populations were greatest in prepupae and early adult life stages, stages during which differences in male horn development and female fecundity manifest. Finally, genes differentially expressed between exotic populations also had greater genetic differentiation and performed functions related to chitin biosynthesis and nutrient sensing, possibly underlying allometry and fecundity trait divergences. Our results suggest that beyond bottlenecks, recent introductions have led to genetic and transcriptional differences in genes correlated with observed phenotypic differences.

Nonadditive Transcriptomic Signatures of Genotype-by-Genotype Interactions during the Initiation of Plant-Rhizobium Symbiosis

Rhizobia are ecologically important, facultative plant-symbiotic microbes. In nature, there is a large variability in the association of rhizobial strains and host plants of the same species. Here, we evaluated whether plant and rhizobial genotypes influence the initial transcriptional response of rhizobium following perception of a host plant. RNA sequencing of the model rhizobium Sinorhizobium meliloti exposed to root exudates or luteolin (an inducer of nod genes, involved in the early steps of symbiotic interaction) was performed on a combination of three S. meliloti strains and three alfalfa varieties as host plants. The response to root exudates involved hundreds of changes in the rhizobium transcriptome. Of the differentially expressed genes, 35% were influenced by the strain genotype, 16% were influenced by the plant genotype, and 29% were influenced by strain-by-host plant genotype interactions. We also examined the response of a hybrid S. meliloti strain in which the symbiotic megaplasmid (∼20% of the genome) was mobilized between two of the above-mentioned strains. Dozens of genes were upregulated in the hybrid strain, indicative of nonadditive variation in the transcriptome. In conclusion, this study demonstrated that transcriptional responses of rhizobia upon perception of legumes are influenced by the genotypes of both symbiotic partners and their interaction, suggesting a wide spectrum of genetic determinants involved in the phenotypic variation of plant-rhizobium symbiosis.IMPORTANCE A sustainable way for meeting the need of an increased global food demand should be based on a holobiont perspective, viewing crop plants as intimately associated with their microbiome, which helps improve plant nutrition, tolerance to pests, and adverse climate conditions. However, the genetic repertoire needed for efficient association with plants by the microbial symbionts is still poorly understood. The rhizobia are an exemplary model of facultative plant symbiotic microbes. Here, we evaluated whether genotype-by-genotype interactions could be identified in the initial transcriptional response of rhizobium perception of a host plant. We performed an RNA sequencing study to analyze the transcriptomes of different rhizobial strains elicited by root exudates of three alfalfa varieties as a proxy of an early step of the symbiotic interaction. The results indicated strain- and plant variety-dependent variability in the observed transcriptional changes, providing fundamentally novel insights into the genetic basis of rhizobium-plant interactions. Our results provide genetic insights and perspective to aid in the exploitation of natural rhizobium variation for improvement of legume growth in agricultural ecosystems.

Nonadditive Transcriptomic Signatures of Genotype-by-Genotype Interactions during the Initiation of Plant-Rhizobium Symbiosis

Rhizobia are ecologically important, facultative plant-symbiotic microbes. In nature, there is a large variability in the association of rhizobial strains and host plants of the same species. Here, we evaluated whether plant and rhizobial genotypes influence the initial transcriptional response of rhizobium following perception of a host plant. RNA sequencing of the model rhizobium Sinorhizobium meliloti exposed to root exudates or luteolin (an inducer of nod genes, involved in the early steps of symbiotic interaction) was performed on a combination of three S. meliloti strains and three alfalfa varieties as host plants. The response to root exudates involved hundreds of changes in the rhizobium transcriptome. Of the differentially expressed genes, 35% were influenced by the strain genotype, 16% were influenced by the plant genotype, and 29% were influenced by strain-by-host plant genotype interactions. We also examined the response of a hybrid S. meliloti strain in which the symbiotic megaplasmid (∼20% of the genome) was mobilized between two of the above-mentioned strains. Dozens of genes were upregulated in the hybrid strain, indicative of nonadditive variation in the transcriptome. In conclusion, this study demonstrated that transcriptional responses of rhizobia upon perception of legumes are influenced by the genotypes of both symbiotic partners and their interaction, suggesting a wide spectrum of genetic determinants involved in the phenotypic variation of plant-rhizobium symbiosis.IMPORTANCE A sustainable way for meeting the need of an increased global food demand should be based on a holobiont perspective, viewing crop plants as intimately associated with their microbiome, which helps improve plant nutrition, tolerance to pests, and adverse climate conditions. However, the genetic repertoire needed for efficient association with plants by the microbial symbionts is still poorly understood. The rhizobia are an exemplary model of facultative plant symbiotic microbes. Here, we evaluated whether genotype-by-genotype interactions could be identified in the initial transcriptional response of rhizobium perception of a host plant. We performed an RNA sequencing study to analyze the transcriptomes of different rhizobial strains elicited by root exudates of three alfalfa varieties as a proxy of an early step of the symbiotic interaction. The results indicated strain- and plant variety-dependent variability in the observed transcriptional changes, providing fundamentally novel insights into the genetic basis of rhizobium-plant interactions. Our results provide genetic insights and perspective to aid in the exploitation of natural rhizobium variation for improvement of legume growth in agricultural ecosystems.

Genomic Analyses of Phenotypic Differences Between Native and Invasive Populations of Diffuse Knapweed (Centaurea diffusa)

Invasive species represent excellent opportunities to study the evolutionary potential of traits important to success in novel environments. Although some ecologically important traits have been identified in invasive species, little is typically known about the genetic mechanisms that underlie invasion success in non-model species. Here, we use a genome-wide association (GWAS) approach to identify the genetic basis of trait variation in the non-model, invasive, diffuse knapweed [Centaurea diffusa Lam. (Asteraceae)]. To assist with this analysis, we have assembled the first draft genome reference and fully annotated plastome assembly for this species, and one of the first from this large, weedy, genus, which is of major ecological and economic importance. We collected phenotype data from 372 individuals from four native and four invasive populations of C. diffusa grown in a common environment. Using these individuals, we produced reduced-representation genotype-by-sequencing (GBS) libraries and identified 7,058 SNPs. We identify two SNPs associated with leaf width in these populations, a trait which significantly varies between native and invasive populations. In this rosette forming species, increased leaf width is a major component of increased biomass, a common trait in invasive plants correlated with increased fitness. Finally, we use annotations from Arabidopsis thaliana to identify 98 candidate genes that are near the associated SNPs and highlight several good candidates for leaf width variation.

Phylostratic Shift of Whole-Genome Duplications in Normal Mammalian Tissues towards Unicellularity Is Driven by Developmental Bivalent Genes and Reveals a Link to Cancer

Tumours were recently revealed to undergo a phylostratic and phenotypic shift to unicellularity. As well, aggressive tumours are characterized by an increased proportion of polyploid cells. In order to investigate a possible shared causation of these two features, we performed a comparative phylostratigraphic analysis of ploidy-related genes, obtained from transcriptomic data for polyploid and diploid human and mouse tissues using pairwise cross-species transcriptome comparison and principal component analysis. Our results indicate that polyploidy shifts the evolutionary age balance of the expressed genes from the late metazoan phylostrata towards the upregulation of unicellular and early metazoan phylostrata. The up-regulation of unicellular metabolic and drug-resistance pathways and the downregulation of pathways related to circadian clock were identified. This evolutionary shift was associated with the enrichment of ploidy with bivalent genes (p < 10⁻¹⁶). The protein interactome of activated bivalent genes revealed the increase of the connectivity of unicellulars and (early) multicellulars, while circadian regulators were depressed. The mutual polyploidy-c-MYC-bivalent genes-associated protein network was organized by gene-hubs engaged in both embryonic development and metastatic cancer including driver (proto)-oncogenes of viral origin. Our data suggest that, in cancer, the atavistic shift goes hand-in-hand with polyploidy and is driven by epigenetic mechanisms impinging on development-related bivalent genes.

Differential immune gene expression associated with contemporary range expansion in two invasive rodents in Senegal

Biological invasions are major anthropogenic changes associated with threats to biodiversity and health. However, what determines the successful establishment and spread of introduced populations remains unclear. Here, we explore several hypotheses linking invasion success and immune phenotype traits, including those based on the evolution of increased competitive ability concept. We compared gene expression profiles between anciently and recently established populations of two major invading species, the house mouse Mus musculus domesticus and the black rat Rattus rattus, in Senegal (West Africa). Transcriptome analyses identified differential expression between anciently and recently established populations for 364 mouse genes and 83 rat genes. All immune-related genes displaying differential expression along the mouse invasion route were overexpressed at three of the four recently invaded sites studied. Complement activation pathway genes were overrepresented among these genes. By contrast, no particular immunological process was found to be overrepresented among the differentially expressed genes of black rat. Changes in transcriptome profiles were thus observed along invasion routes, but with different specific patterns between the two invasive species. These changes may be driven by increases in infection risks at sites recently invaded by the house mouse, and by stochastic events associated with colonization history for the black rat. These results constitute a first step toward the identification of immune eco-evolutionary processes potentially involved in the invasion success of these two rodent species.

Abiotic and past climatic conditions drive protein abundance variation among natural populations of the caddisfly Crunoecia irrorata

Deducing impacts of environmental change on species and the populations they form in nature is an important goal in contemporary ecology. Achieving this goal is hampered by our limited understanding of the influence of naturally occurring environmental variation on the molecular systems of ecologically relevant species, as the pathways underlying fitness-affecting plastic responses have primarily been studied in model organisms and under controlled laboratory conditions. Here, to test the hypothesis that proteome variation systematically relates to variation in abiotic conditions, we establish such relationships by profiling the proteomes of 24 natural populations of the spring-dwelling caddisfly Crunoecia irrorata. We identified protein networks whose abundances correlated with environmental (abiotic) gradients such as in situ pH, oxygen- and nitrate concentrations but also climatic data such as past thermal minima and temperature seasonality. Our analyses suggest that variations in abiotic conditions induce discrete proteome responses such as the differential abundance of proteins associated with cytoskeletal function, heat-shock proteins and proteins related to post-translational modification. Identifying these drivers of proteome divergence characterizes molecular "noise", and positions it as a background against which molecular signatures of species' adaptive responses to stressful conditions can be identified.

Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions

Phenotypic plasticity has been increasingly recognized for its importance in adaptation to novel environments, and initial rapid plastic response to acute stresses usually serves as the stepping stone for future adaptation. Differential gene expression and alternative splicing have been proposed as two underlying mechanisms for rapid plastic response to environmental stresses. Here, we used an invasive model species, Ciona savignyi, to investigate the temporary plastic changes under temperature stresses on gene expression and alternative splicing. Our results revealed rapid and highly dynamic gene expression reprogramming and alternative splicing switch under acute stresses. Distinct transcriptional response profiles were triggered by two types of temperature stresses, showing resilience recovery and increasing divergence under heat and cold challenges, respectively. Interestingly, alternative exons were more inclined to be skipped under both heat and cold stresses, leading to shorter isoforms but with maintained Open Reading Frames (ORFs). Although similar response patterns were observed between differential gene expression and alternative splicing, low overlap between Differentially Expressed Genes (DEGs) and Differentially Alternative Spliced Genes (DASGs) suggests that distinct gene sets and associated functions should be involved in temperature challenges. Thus, alternative splicing should offer an additional layer of plastic response to environmental challenges. Finally, we identified key plastic genes involved in both gene expression regulation and alternative splicing. The results obtained here shed light on adaptation and accommodation mechanisms during biological invasions, particularly for acute environmental changes at early stages of biological invasions such as transport and introduction.

Transcriptome profiling revealed potentially important roles of defensive gene expression in the divergence of insect biotypes: a case study with the cereal aphid Sitobion avenae

Background

Many insects can develop differential biotypes on variable host plants, but the underlying molecular factors and mechanisms are not well understood. To address this issue, transcriptome profiling analyses were conducted for two biotypes of the cereal aphid, Sitobion avenae (Fabricius), on both original and alternative plants.

Results

Comparisons between both biotypes generated 4174 differentially expressed unigenes (DEGs). In their response to host plant shift, 39 DEGs were shared by both biotypes, whereas 126 and 861 DEGs occurred only in biotypes 1 and 3, respectively. MMC (modulated modularity clustering) analyses showed that specific DEGs of biotypes 1 and 3 clustered into five and nine transcriptional modules, respectively. Among these DEGs, defense-related genes underwent intensive expression restructuring in both biotypes. However, biotype 3 was found to have relatively lower gene transcriptional plasticity than biotype 1. Gene enrichment analyses of the abovementioned modules showed functional divergence in defensive DEGs for the two biotypes in response to host transfer. The expression plasticity for some defense related genes was showed to be directly related to fecundity of S. avenae biotypes on both original and alternative plants, suggesting that expression plasticity of key defensive genes could have significant impacts on the adaptive potential and differentiation of S. avenae biotypes on different plants.

Conclusions

The divergence patterns of transcriptional plasticity in defense related genes may play important roles in the phenotypic evolution and differentiation of S. avenae biotypes. Our results can provide insights into the role of gene expression plasticity in the divergence of insect biotypes and adaptive evolution of insect populations.

Ecological divergence in sympatry causes gene misexpression in hybrids

Ecological speciation occurs when reproductive isolation evolves as a byproduct of adaptive divergence between populations. Selection favouring gene regulatory divergence between species could result in transgressive levels of gene expression in F1 hybrids that may lower hybrid fitness. We combined 58 resequenced genomes with 124 transcriptomes to identify patterns of hybrid gene misexpression that may be driven by adaptive regulatory divergence within a young radiation of Cyprinodon pupfishes, which consists of a dietary generalist and two trophic specialists-a molluscivore and a scale-eater. We found more differential gene expression between closely related sympatric specialists than between allopatric generalist populations separated by 1,000 km. Intriguingly, 9.6% of genes that were differentially expressed between sympatric species were also misexpressed in F1 hybrids. A subset of these genes were in highly differentiated genomic regions and enriched for functions important for trophic specialization, including head, muscle and brain development. These regions also included genes that showed evidence of hard selective sweeps and were significantly associated with oral jaw length-the most rapidly diversifying skeletal trait in this radiation. Our results indicate that divergent ecological selection in sympatry can contribute to hybrid gene misexpression which may act as a reproductive barrier between nascent species.

Evolutionary Patterns of Sex-Biased Genes in Three Species of Haplodiploid Insects

Females and males often differ obviously in morphology and behavior, and the differences between sexes are the result of natural selection and/or sexual selection. To a great extent, the differences between the two sexes are the result of differential gene expression. In haplodiploid insects, this phenomenon is obvious, since males develop from unfertilized zygotes and females develop from fertilized zygotes. Whiteflies of the Bemisia tabaci species complex are typical haplodiploid insects, and some species of this complex are important pests of many crops worldwide. Here, we report the transcriptome profiles of males and females in three species of this whitefly complex. Between-species comparisons revealed that non-sex-biased genes display higher variation than male-biased or female-biased genes. Sex-biased genes evolve at a slow rate in protein coding sequences and gene expression and have a pattern of evolution that differs from those of social haplodiploid insects and diploid animals. Genes with high evolutionary rates are more related to non-sex-biased traits-such as nutrition, immune system, and detoxification-than to sex-biased traits, indicating that the evolution of protein coding sequences and gene expression has been mainly driven by non-sex-biased traits.

Genome-Wide Genotype-Expression Relationships Reveal Both Copy Number and Single Nucleotide Differentiation Contribute to Differential Gene Expression between Stickleback Ecotypes

Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.

Linking stress coping styles with brain mRNA abundance of selected transcripts for Senegalese sole (Solea senegalensis) juveniles

In fish, proactive and reactive individual stress copying styles (SCS) have been used to resolve variation in molecular expression data. Stress coping styles have been previously described in several stages of Solea senegalensis by validating for the species the use of standard behavioural screening tests. The present study aimed to link behavioural SCS tests with brain transcript abundance in early Senegalese sole juveniles in order to observe the natural variation in a molecular pathway in this species. A total of 50 juveniles were subjected to three individual behavioural (Restraining, New environment and Confinement) and one group (Risk-taking) screening tests. The fish were classified in SCS categories by applying a hierarchical cluster to the variable "Total activity" (the total activity time that the fish was moving in each individual test). Three categories were defined, proactive, intermediate and reactive sole. Six transcripts were chosen and tested, one related to basic metabolism (gapdh-2), three to feeding behaviour (per1, igf-Ia, pparß) and two to the stress response (crh-BP and hsp90aa) in 30 juveniles (10 individuals per SCS category) using rt-qPCR to observe differences in the abundance of those transcripts among SCS. Four transcripts were differentially expressed (DETs) among them. The transcript gapdh-2 showed up-regulation for proactive and intermediate SCS sole while reactive individuals showed down-regulation. Target mRNAs per1, igf-Ia and pparß, showed different levels of up-regulation for proactive and reactive fish while intermediates were highly down-regulated. Surprisingly no differences in stress related transcripts were observed. Correlations were found between variation in coping styles and variation in the abundance of mRNAs involved in important biological functions in Senegalese sole. These results are the first evidence of the relationship between the behavioural individual variation and the fluctuation in brain transcripts abundance in Senegalese sole.

Drift and Directional Selection Are the Evolutionary Forces Driving Gene Expression Divergence in Eye and Brain Tissue of Heliconius Butterflies

Investigating gene expression evolution over micro- and macroevolutionary timescales will expand our understanding of the role of gene expression in adaptation and speciation. In this study, we characterized the evolutionary forces acting on gene expression levels in eye and brain tissue of five Heliconius butterflies with divergence times of ∼5-12 MYA. We developed and applied Brownian motion (BM) and Ornstein-Uhlenbeck (OU) models to identify genes whose expression levels are evolving through drift, stabilizing selection, or a lineage-specific shift. We found that 81% of the genes evolve under genetic drift. When testing for branch-specific shifts in gene expression, we detected 368 (16%) shift events. Genes showing a shift toward upregulation have significantly lower gene expression variance than those genes showing a shift leading toward downregulation. We hypothesize that directional selection is acting in shifts causing upregulation, since transcription is costly. We further uncovered through simulations that parameter estimation of OU models is biased when using small phylogenies and only becomes reliable with phylogenies having ≥ 50 taxa. Therefore, we developed a new statistical test based on BM to identify highly conserved genes (i.e., evolving under strong stabilizing selection), which comprised 3% of the orthoclusters. In conclusion, we found that drift is the dominant evolutionary force driving gene expression evolution in eye and brain tissue in Heliconius Nevertheless, the higher proportion of genes evolving under directional than under stabilizing selection might reflect species-specific selective pressures on vision and the brain that are necessary to fulfill species-specific requirements.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

The Evolution of Gene Expression Underlying Vision Loss in Cave Animals

Dissecting the evolutionary genetic processes underlying eye reduction and vision loss in obligate cave-dwelling organisms has been a long-standing challenge in evolutionary biology. Independent vision loss events in related subterranean organisms can provide critical insight into these processes as well as into the nature of convergent loss of complex traits. Advances in evolutionary developmental biology have illuminated the significant role of heritable gene expression variation in the evolution of new forms. Here, we analyze gene expression variation in adult eye tissue across the freshwater crayfish, representing four independent vision-loss events in caves. Species and individual expression patterns cluster by eye function rather than phylogeny, suggesting convergence in transcriptome evolution in independently blind animals. However, this clustering is not greater than what is observed in surface species with conserved eye function after accounting for phylogenetic expectations. Modeling expression evolution suggests that there is a common increase in evolutionary rates in the blind lineages, consistent with a relaxation of selective constraint maintaining optimal expression levels. This is evidence for a repeated loss of expression constraint in the transcriptomes of blind animals and that convergence occurs via a similar trajectory through genetic drift.

Genome-Wide Microarray Analysis Suggests Transcriptomic Response May Not Play a Major Role in High- to Low-Altitude Acclimation in Harvest Mouse (Micromys minutus)

Simple Summary

Micromys minutus is a small rodent species that has a wide range of vertical distribution in Taiwan. By comparing the gene expression profile of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site, the Tnfrsf12a gene was demonstrated to have a differential expression pattern. Although this finding may be correlated with the altitude acclimation, the observation of only one gene transcript with significant alteration leads us to suggest that genetic response may not play a major role in altitude acclimation in M. minutus. Future comparative functional genomics studies involving other organ systems (in addition to skeletal muscles) and alarger sample size are warranted for better insight into the altitude acclimation of this small rodent species.

Abstract

The harvest mouse (Micromys minutus) is a small rodent species with a wide range of vertical distribution in Taiwan, extending from the sea level to 3100 m altitude. This species has recently suffered from habitat loss in high-altitude areas due to orchard cultivation, which may have resulted in mouse migration from high to low altitude. To investigate whether there is any physiological mechanism involved in altitude acclimation, rat cDNA microarray was used to compare transcriptomic patterns of the skeletal muscle tissues taken from individuals native to the high-altitude environment and those transferred to the low-altitude captive site. Of the 23,188 genes being analyzed, 47 (33 up-regulated and 14 down-regulated) were found to have differential expression (fold change > 4 or < −4, ANOVA p < 0.05). However, after multiple testing correction with a false discovery rate (FDR), only the result for Tnfrsf12a was found to be statistically significant (fold change = 13, FDR p < 0.05). The result was confirmed by quantitative polymerase chain reaction (q-PCR). The expression of Tnfrsf12a possibly relates to the skeletal muscle biology and thus can be correlated with altitude acclimation. However, finding only one gene transcript with significant alteration suggests that transcriptomic response may not play a major role in high- to low-altitude acclimation in harvest mouse.

Phototransduction Gene Expression and Evolution in Cave and Surface Crayfishes

In the absence of light in caves, animals have repeatedly evolved reduced eyes and visual systems. Whether the underlying genetic components remain intact in blind species remains unanswered across taxa. The freshwater crayfish have evolved to live in caves multiple times throughout their history; therefore, this system provides an opportunity to probe the genetic patterns and processes underlying repeated vision loss. Using transcriptomic data from the eyes of 14 species of cave and surface crayfishes, we identify the expression of 17 genes putatively related to visual phototransduction. We find a similarly complete repertoire of phototransduction gene families expressed in cave and surface species, but that the expression levels of those transcripts are consistently lower in cave species. We find statistical support for episodic positive selection, increased and decreased selection strength in caves, depending on the gene family. Analyses of gene expression evolution suggest convergent and possibly adaptive downregulation of these genes across eye-reduction events. Our results reveal a combination of evolutionary processes acting on the sequences and gene expression levels of vision-related genes underlying the loss of vision in caves.

Evolution of DNA Methylation in Papio Baboons

Changes in gene regulation have long been thought to play an important role in primate evolution. However, although a number of studies have compared genome-wide gene expression patterns across primate species, fewer have investigated the gene regulatory mechanisms that underlie such patterns, or the relative contribution of drift versus selection. Here, we profiled genome-scale DNA methylation levels in blood samples from five of the six extant species of the baboon genus Papio (4-14 individuals per species). This radiation presents the opportunity to investigate DNA methylation divergence at both shallow and deeper timescales (0.380-1.4 My). In contrast to studies in human populations, but similar to studies in great apes, DNA methylation profiles clearly mirror genetic and geographic structure. Divergence in DNA methylation proceeds fastest in unannotated regions of the genome and slowest in regions of the genome that are likely more constrained at the sequence level (e.g., gene exons). Both heuristic approaches and Ornstein-Uhlenbeck models suggest that DNA methylation levels at a small set of sites have been affected by positive selection, and that this class is enriched in functionally relevant contexts, including promoters, enhancers, and CpG islands. Our results thus indicate that the rate and distribution of DNA methylation changes across the genome largely mirror genetic structure. However, at some CpG sites, DNA methylation levels themselves may have been a target of positive selection, pointing to loci that could be important in connecting sequence variation to fitness-related traits.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.

Variable gene transcription underlies phenotypic convergence of hypoxia tolerance in sculpins

Background

The degree by which mechanisms underlying phenotypic convergence are similar among taxa depends on the number of evolutionary paths available for selection to act upon. Likelihood of convergence will be influenced by an interplay of factors such as genetic architecture, phylogenetic history and population demography. To determine if there is convergence or divergence in mechanisms underlying phenotypic similarity, we assessed whether gene transcription patterns differed among species with similar levels of hypoxia tolerance.

Results

Three species of marine fish from the superfamily Cottoidea (smoothhead sculpin [Artedius lateralis], sailfin sculpin [Nautichthys oculofasciatus] and Pacific staghorn sculpin [Leptocottus armatus]), all of which have previously been shown to share the same level of hypoxia tolerance, were exposed to short-(8 h) and longer-term (72 h) hypoxia and mRNA transcripts were assessed using a custom microarray. We examined hypoxia-induced transcription patterns in metabolic and protein production pathways and found that a high proportion of genes associated with these biological processes showed significant differences among the species. Specifically, the data suggest that the smoothhead sculpin, unlike the sailfin sculpin and the Pacific staghorn sculpin, relied on amino acid degradation rather than glycolysis or fatty acid oxidation to generate ATP during hypoxia exposure. There was also variation across the species in the transcription of genes involved in protein production (e.g. mRNA processing and protein translation), such that it increased in the smoothhead sculpin, decreased in the sailfin sculpin and was variable in the Pacific staghorn sculpin.

Conclusions

Changes in metabolic and protein production pathways are part of the key responses of fishes to exposures to environmental hypoxia. Yet, species with similar overall hypoxia tolerance exhibited different transcriptional responses in these pathways, indicating flexibility and complexity of interactions in the evolution of the mechanisms underlying the hypoxia tolerance phenotype. The variation in the hypoxia-induced transcription of genes across species with similar hypoxia tolerance suggests that similar whole-animal phenotypes can emerge from divergent evolutionary paths that may affect metabolically important functions.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1275-1) contains supplementary material, which is available to authorized users.

Comparative cochlear transcriptomics of echolocating bats provides new insights into different nervous activities of CF bat species

The molecular mechanisms used by echolocating bats to deal with different ultrasonic signals remain to be revealed. Here, we utilised RNA-Seq data to conduct comparative cochlear transcriptomics to assess the variation of gene expression among bats with three types of echolocation: constant-frequency (CF) bats, frequency-modulated (FM) bats and click bats. Our results suggest larger differences in gene expression between CF and click bats than between CF and FM bats and small differences between FM and click bats. We identified 426 and 1,504 differentially expressed genes (DEGs) by the different methods as functionally important for CF bats, in that they showed consistent upregulation in the cochlea of two CF bats, relative to the levels in click and FM bats. Subsequently, downstream GO and KEGG enrichment analyses indicated that both the 426 and 1,504 gene sets were associated with changes in nervous activities in the cochleae of CF bats. In addition, another set of 1,764 DEGs were identified to have crucial hearing related physiological functions for laryngeally echolocating bats. Our study provides a comprehensive overview of the genetic basis of differences among echolocating bats, revealing different nervous system activities during auditory perception in the cochlea particularly in CF bats.

Endocrine-immune signaling as a predictor of survival: A prospective study in developing songbird chicks

Immune function varies with an animal's endocrine physiology and energy reserves, as well as its abiotic and biotic environment. This context-dependency is thought to relate to adaptive trade-off resolution that varies from one context to the next; however, it is less clear how state- and environmentally-dependent differences in endocrine-immune signaling relate to survival in natural populations. We begin to address this question in a prospective study on a free-living passerine bird, the tree swallow (Tachycineta bicolor), by capitalizing upon naturally-occurring variation in ectoparasitism in 12-day old chicks. We measured body mass, hematological gene expression of the pro-inflammatory cytokine interleukin-6 (IL-6) as well as corticosterone (CORT) secretion at baseline and in response to 30 min of handling. We found that chicks with ectoparasites had smaller body mass and higher levels of IL-6 gene expression at this critical stage of post-natal growth and development. Mass and IL-6 were positively correlated, but only among parasitized chicks, suggesting that larger chicks mount stronger immune responses when necessary, i.e. in the presence of ectoparasites that are known to induce inflammation. IL-6 mRNA expression was negatively correlated with stress-induced CORT levels, suggesting that this proxy of inflammation may be co-regulated with or coordinated by glucocorticoids. More importantly, these endocrine-immune parameters predicted survival to fledging, which was positively associated with IL-6 mRNA abundance and, to a lesser degree, CORT reactivity. These results suggest a link between endocrine-immune interactions and performance in nature, and as a consequence, they shed light on the potentially adaptive, context-dependent interplay between body mass, immunity, and endocrine physiology during development.

Somatic polyploidy is associated with the upregulation of c-MYC interacting genes and EMT-like signature

The dependence of cancer on overexpressed c-MYC and its predisposition for polyploidy represents a double puzzle. We address this conundrum by cross-species transcription analysis of c-MYC interacting genes in polyploid vs. diploid tissues and cells, including human vs. mouse heart, mouse vs. human liver and purified 4n vs. 2n mouse decidua cells. Gene-by-gene transcriptome comparison and principal component analysis indicated that c-MYC interactants are significantly overrepresented among ploidy-associated genes. Protein interaction networks and gene module analysis revealed that the most upregulated genes relate to growth, stress response, proliferation, stemness and unicellularity, as well as to the pathways of cancer supported by MAPK and RAS coordinated pathways. A surprising feature was the up-regulation of epithelial-mesenchymal transition (EMT) modules embodied by the N-cadherin pathway and EMT regulators from SNAIL and TWIST families. Metabolic pathway analysis also revealed the EMT-linked features, such as global proteome remodeling, oxidative stress, DNA repair and Warburg-like energy metabolism. Genes associated with apoptosis, immunity, energy demand and tumour suppression were mostly down-regulated. Noteworthy, despite the association between polyploidy and ample features of cancer, polyploidy does not trigger it. Possibly it occurs because normal polyploidy does not go that far in embryonalisation and linked genome destabilisation. In general, the analysis of polyploid transcriptome explained the evolutionary relation of c-MYC and polyploidy to cancer.

Differential gene expression during early development in recently evolved and sympatric Arctic charr morphs

Phenotypic differences between closely related taxa or populations can arise through genetic variation or be environmentally induced, leading to altered transcription of genes during development. Comparative developmental studies of closely related species or variable populations within species can help to elucidate the molecular mechanisms related to evolutionary divergence and speciation. Studies of Arctic charr (Salvelinus alpinus) and related salmonids have revealed considerable phenotypic variation among populations and in Arctic charr many cases of extensive variation within lakes (resource polymorphism) have been recorded. One example is the four Arctic charr morphs in the ∼10,000 year old Lake Thingvallavatn, which differ in numerous morphological and life history traits. We set out to investigate the molecular and developmental roots of this polymorphism by studying gene expression in embryos of three of the morphs reared in a common garden set-up. We performed RNA-sequencing, de-novo transcriptome assembly and compared gene expression among morphs during an important timeframe in early development, i.e., preceding the formation of key trophic structures. Expectedly, developmental time was the predominant explanatory variable. As the data were affected by some form of RNA-degradation even though all samples passed quality control testing, an estimate of 3'-bias was the second most common explanatory variable. Importantly, morph, both as an independent variable and as interaction with developmental time, affected the expression of numerous transcripts. Transcripts with morph effect, separated the three morphs at the expression level, with the two benthic morphs being more similar. However, Gene Ontology analyses did not reveal clear functional enrichment of transcripts between groups. Verification via qPCR confirmed differential expression of several genes between the morphs, including regulatory genes such as AT-Rich Interaction Domain 4A (arid4a) and translin (tsn). The data are consistent with a scenario where genetic divergence has contributed to differential expression of multiple genes and systems during early development of these sympatric Arctic charr morphs.