The power and promise of RNA-seq in ecology and evolution

Accelerating crop improvement via integration of transcriptome-based network biology and genome editing

MAIN CONCLUSION

Big data and network biology infer functional coupling between genes. In combination with machine learning, network biology can dramatically accelerate the pace of gene discovery using modern transcriptomics approaches and be validated via genome editing technology for improving crops to stresses. Unlike other living things, plants are sessile and frequently face various environmental challenges due to climate change. The cumulative effects of combined stresses can significantly influence both plant growth and yields. In navigating the complexities of climate change, ensuring the nourishment of our growing population hinges on implementing precise agricultural systems. Conventional breeding methods have been commonly employed; however, their efficacy has been impeded by limitations in terms of time, cost, and infrastructure. Cutting-edge tools focussing on big data are being championed to usher in a new era in stress biology, aiming to cultivate crops that exhibit enhanced resilience to multifactorial stresses. Transcriptomics, combined with network biology and machine learning, is proving to be a powerful approach for identifying potential genes to target for gene editing, specifically to enhance stress tolerance. The integration of transcriptomic data with genome editing can yield significant benefits, such as gaining insights into gene function by modifying or manipulating of specific genes in the target plant. This review provides valuable insights into the use of transcriptomics platforms and the application of biological network analysis and machine learning in the discovery of novel genes, thereby enhancing the understanding of plant responses to combined or sequential stress. The transcriptomics as a forefront omics platform and how it is employed through biological networks and machine learning that lead to novel gene discoveries for producing multi-stress-tolerant crops, limitations, and future directions have also been discussed.

Comparative transcriptomics reveals potential regulators of climate adaptation in Santalum album L. (Indian Sandalwood)

Santalum album L. (Indian Sandalwood), a valued tree species known for its fragrant heartwood and essential oil is facing increasing threat due to severe anthropogenic pressures compounded by climate change which has resulted in depletion of its adaptive gene pool. The present study investigates the transcriptome-level responses of nine sandalwood genotypes sourced from diverse climatic zones to identify adaptive genes in the species. Comparative transcriptomics predicted 727, 1141 and 479 differentially expressed transcripts (DETs) across wet vs. dry; monsoon vs. dry and wet vs. monsoon conditions, respectively, and majority of DETs were up-regulated in samples sourced from high rainfall areas. Transcripts including heat shock proteins, Zinc finger binding protein, ribosomal proteins, transcription factors and protein kinase were identified as probable regulators of climate adaptation in S. album. The expression changes of eight selected transcripts were further validated by real-time quantitative PCR. Protein-protein interaction analysis revealed key hub transcripts involved in climate response, while alternative splicing events in transcripts such as SURP and G-patch domain-containing protein 1-like protein, G-type lectin S-receptor-like serine/threonine protein kinase B120, Tetraspanin-3 and ARM repeat superfamily protein indicated the probable role of alternate splicing in increasing the transcript diversity during adaptation. This study presents the first insight into the molecular mechanisms of climate adaptation in the species and can form the basis for specific interventions such as selective breeding, genetic manipulation, and habitat management for conservation and long-term survival of sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-025-04218-4.

Metabolic differentiation of brushtail possum populations resistant and susceptible to plant toxins revealed via differential gene expression

The Australian brushtail possum (Trichosurus vulpecula) is adapted to a wide range of food plants across its range and is exposed to numerous physiological challenges. Populations that are resistant to the plant toxin sodium fluoroacetate are of particular interest as this compound has been used since the 1940s for vertebrate pest management around the world. Candidate gene identification is an important first step in understanding how spatial populations have responded to local selection resulting in local physiological divergence. We employ differential gene expression of liver samples from wild-caught brushtail possums from toxin-resistant and toxin-susceptible populations to identify candidate genes that might be involved in metabolic pathways associated with toxin-resistance. This allowed us to identify genetic pathways involved in resistance to the plant toxin sodium fluoroacetate in Western Australian possums but not those originally from south eastern Australia. We identified differentially expressed genes in the liver that are associated with cell signalling, encapsulating structure, cell mobility, and tricarboxylic acid cycle. The gene expression differences detected indicate which metabolic pathways are most likely to be associated with sodium fluoroacetate resistance in these marsupials and we provide a comprehensive list of candidate genes and pathways to focus on for future studies.

Thermal variation influences the transcriptome of the major malaria vector Anopheles stephensi

The distribution and abundance of ectothermic mosquitoes are strongly affected by temperature, but mechanisms remain unexplored. We describe the effect of temperature on the transcriptome of Anopheles stephensi, an invasive vector of human malaria. Adult females were maintained across a range of mean temperatures (20 °C, 24 °C and 28 °C), with daily fluctuations of +5 °C and -4 °C at each mean temperature. Transcriptomes were described up to 19 days post-blood meal. Of the >3100 differentially expressed genes, we observed shared temporal expression profiles across all temperatures, suggesting their indispensability to mosquito life history. Tolerance to 20 and 28 ( + 5°C/-4°C) was associated with larger and more diverse transcriptomes compared to 24 ( + 5 °C/-4 °C). Finally, we identified two distinct trends in gene expression in response to blood meal ingestion, oxidative stress, and reproduction. Our work has implications for mosquitoes' response to thermal variation, mosquito immune-physiology, mosquito-malaria interactions and the development of vector control tools.

Experimental reporting of fish transcriptomic responses in environmental toxicology and ecotoxicology

Due to its increasing affordability and efforts to understand transcriptional responses of organisms to biotic and abiotic stimuli, transcriptomics has become an important tool with significant impact on toxicological investigations and hazard and risk assessments, especially during development and application of new approach methodologies (NAMs). Data generated using transcriptomic methodologies have directly informed adverse outcome pathway frameworks, chemical and biological read across, and aided in the identification of points of departure. Using data reporting frameworks for transcriptomics data offers improved transparency and reproducibility of research and an opportunity to identify barriers to adoption of these NAMs, especially in environmental toxicology and ecotoxicology with aquatic models. Improved reporting also allows for reexamination of existing data, limiting needs for experiment replication and further reducing animal experimentation. Here, we use a standardized form of data reporting for omics-based studies, the Organisation for Economic Co-operation and Development omics reporting framework, which specifically reports on a list of parameters that should be included in transcriptomics studies used in a regulatory context. We focused specifically on fish studies using RNA- Sequencing (Seq)/microarray technologies within a toxicology context. Inconsistencies in reporting and methodologies among the experimental designs (toxicology vs. molecular characterization) were observed in addition to foundational differences in reporting of sample concentration or preparation or quality assessments, which can affect reproducibility and read across, confidence in results, and contribute substantially to understanding molecular mechanisms of toxicants and toxins. Our findings present an opportunity for improved research reporting. We also provide several recommendations as logical steps to reduce barriers to adoption of transcriptomics within environmental toxicology and ecotoxicology.

Gene expression comparisons between captive and wild shrew brains reveal captivity effects

Compared with their free-ranging counterparts, wild animals in captivity experience different conditions with lasting physiological and behavioural effects. Although shifts in gene expression are expected to occur upstream of these phenotypes, we found no previous gene expression comparisons of captive versus free-ranging mammals. We assessed gene expression profiles of three brain regions (cortex, olfactory bulb and hippocampus) of wild shrews (Sorex araneus) compared with shrews kept in captivity for two months and undertook sample dropout to examine robustness given limited sample sizes. Consistent with captivity effects, we found hundreds of differentially expressed genes in all three brain regions, 104 overlapping across all three, that enriched pathways associated with neurodegenerative disease, oxidative phosphorylation and genes encoding ribosomal proteins. In the shrew, transcriptomic changes detected under captivity resemble responses in several human pathologies, including major depressive disorder and neurodegeneration. While interpretations of individual genes are tempered by small sample sizes, we propose captivity influences brain gene expression and function and can confound analyses of natural processes in wild individuals under captive conditions.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

The roles of different gene expression regulators in acoustic variation in the intermediate horseshoe bat revealed by long-read and short-read RNA sequencing data

Gene expression changes contribute greatly to phenotypic variations in nature. Studying patterns of regulators of gene expression is important to fully understand the molecular mechanism underlying phenotypic variations. In horseshoe bats, the cochleae are finely tuned to echoes of call frequency. Here, using 2 recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis hainanus and R. a. himalayanus) with great acoustic variations as the system, we aim to explore relative roles of different regulators of gene expression (differential gene expression, alternative splicing (AS) and long non-coding RNAs (lncRNAs)) in phenotypic variation with a combination of Illumina short-read and Nanopore long-read RNA-seq data from the cochlea. Compared to R. a. hainanus, R. a. himalayanus exhibited much more upregulated differentially expressed genes (DEGs) and multiple of them may play important roles in the maintenance and damage repair of auditory hair cells. We identified 411 differentially expressed lncRNAs and their target DEGs upregulated in R. a. himalayanus were also mainly involved in a protective mechanism for auditory hair cells. Using 3 different methods of AS analysis, we identified several candidate alternatively spliced genes (ASGs) that expressed different isoforms which may be associated with acoustic divergence of the 2 subspecies. We observed significantly less overlap than expected between DEGs and ASGs, supporting complementary roles of differential gene expression and AS in generating phenotypic variations. Overall, our study highlights the importance of a combination of short-read and long-read RNA-seq data in examining the regulation of gene expression changes responsible for phenotypic variations.

Repeated marine heatwaves aggravate the adverse effects of nano-TiO2 on physiological metabolism of the thick-shelled mussel Mytilus coruscus

Global climate change is a major trigger of unexpected temperature fluctuations. The impacts of marine heatwaves (MHWs) and nano-titanium dioxide (nano-TiO2) on marine organisms have been extensively investigated. However, the potential mechanisms underlying their interactive effects on physiological processes and metabolism remain poorly understood, especially regarding periodic MHWs in real-world conditions. In this study, the effects of nano-TiO2 (at concentrations of 0, 25, and 250 μg/L) and periodic MHWs on the condition index (CI) and underlying metabolic mechanisms were investigated in mussels (Mytilus coruscus). The results showed that mussels try to upregulate their respiration rate (RR) to enhance aerobic metabolism (indicated by elevated succinate dehydrogenase) under short-term nano-TiO2 exposure. However, even at ambient concentration (25 μg/L), prolonged nano-TiO2 exposure inhibited ingestion ability (decreased clearance rate) and glycolysis (inhibited pyruvate kinase, hexokinase, and phosphofructokinase activities), which led to an insufficient energy supply (decreased triglyceride, albumin, and ATP contents). Repeated thermal scenarios caused more severe physiological damage, demonstrating that mussels are fragile to periodic MHWs. MHWs decreased the zeta potential of the nano-TiO2 particles but increased the hydrodynamic diameter. Additionally, exposure to nano-TiO2 and periodic MHWs further affected aerobic respiration (inhibited lactate dehydrogenase and succinate dehydrogenase activities), metabolism (decreased RR, activities of respiratory metabolism-related enzymes, and expressions of PEPCK, PPARγ, and ACO), and overall health condition (decreased ATP and CI). These findings indicate that the combined stress of these two stressors exerts more detrimental impact on the physiological performance and energy metabolism of mussels, and periodic MHWs exacerbate the toxicological effects of ambient concentration nano-TiO2. Given the potential worsening of nanoparticle pollution and the increase in extreme heat events in the future, the well-being of mussels in the marine environment may face further threats.

Expression-environment associations in transcriptomic heat stress responses for a global plant lineage

The increasing frequency and severity of heatwaves will intensify stress on plants. Given regional variation in heatwave exposure and expected differences in thermal tolerance between species it is unlikely that all plant species will be affected equally by climate change. However, little is currently known about variation in the responses of plants to heat stress, or how those responses differ among closely related species adapted to different environments. Here we quantify the response of 17 Acacia species (175 RNA-seq libraries), from across Australia's diverse biomes, to a multi-day experimental heatwave treatment to identify variation in transcriptomic and physiological responses to heat stress. Genes with known heat response functions showed consistent responses across Acacia species. Up to 10% of all genes and over 100 gene families showed significant clinal variation in the magnitude of their expression plasticity across species. Specifically, gene families linked to the temperature stress response were overrepresented among significant relationships with home range temperature conditions. Gene expression responses seen on the first day of the heatwave were more frequently associated with home range climates, while expression responses by day four were more commonly related to photosystem II acclimation. Comparative transcriptomics on non-model species has the potential to provide key information on stress response plasticity, especially when linked with our understanding of model species. Our study indicates that the pressing challenge to identifying potentially vulnerable species to climate change could be benefited by the further exploration of clinal variation in transcriptome plasticity.

Transcriptomic changes in the posterior pallium of male zebra finches associated with social niche conformance

Animals plastically adjust their physiological and behavioural phenotypes to conform to their social environment-social niche conformance. The degree of sexual competition is a critical part of the social environment to which animals adjust their phenotypes, but the underlying genetic mechanisms are poorly understood. We conducted a study to investigate how differences in sperm competition risk affect the gene expression profiles of the testes and two brain areas (posterior pallium and optic tectum) in breeding male zebra finches (Taeniopygia castanotis). In this pre-registered study, we investigated a large sample of 59 individual transcriptomes. We compared two experimental groups: males held in single breeding pairs (low sexual competition) versus those held in two pairs (elevated sexual competition) per breeding cage. Using weighted gene co-expression network analysis (WGCNA), we observed significant effects of the social treatment in all three tissues. However, only the treatment effects found in the pallium were confirmed by an additional randomisation test for statistical robustness. Likewise, the differential gene expression analysis revealed treatment effects only in the posterior pallium (ten genes) and optic tectum (six genes). No treatment effects were found in the testis at the single gene level. Thus, our experiments do not provide strong evidence for transcriptomic adjustment specific to manipulated sperm competition risk. However, we did observe transcriptomic adjustments to the manipulated social environment in the posterior pallium. These effects were polygenic rather than based on few individual genes with strong effects. Our findings are discussed in relation to an accompanying paper using the same animals, which reports behavioural results consistent with the results presented here.

Metabolomic and transcriptomic analyses highlight metabolic regulatory networks of Salvia miltiorrhiza in response to replant disease

BACKGROUND

Salvia miltiorrhiza, a well-known traditional Chinese medicine, frequently suffers from replant diseases that adversely affect its quality and yield. To elucidate S. miltiorrhiza's metabolic adaptations to replant disease, we analyzed its metabolome and transcriptome, comparing normal and replant diseased plants for the first time.

RESULTS

We identified 1,269 metabolites, 257 of which were differentially accumulated metabolites, and identified 217 differentially expressed genes. Integrated transcriptomic and metabolomic analyses revealed a significant up-regulation and co-expression of metabolites and genes associated with plant hormone signal transduction and flavonoid biosynthesis pathways in replant diseases. Within plant hormone signal transduction pathway, plants afflicted with replant disease markedly accumulated indole-3-acetic acid and abscisic acid, correlating with high expression of their biosynthesis-related genes (SmAmidase, SmALDH, SmNCED, and SmAAOX3). Simultaneously, changes in hormone concentrations activated plant hormone signal transduction pathways. Moreover, under replant disease, metabolites in the local flavonoid metabolite biosynthetic pathway were significantly accumulated, consistent with the up-regulated gene (SmHTC1 and SmHTC2). The qRT-PCR analysis largely aligned with the transcriptomic results, confirming the trends in gene expression. Moreover, we identified 10 transcription factors co-expressed with differentially accumulated metabolites.

CONCLUSIONS

Overall, we revealed the key genes and metabolites of S. miltiorrhiza under replant disease, establishing a robust foundation for future inquiries into the molecular responses to combat replant stress.

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

Effect of RNA-Seq data normalization on protein interactome mapping for Alzheimer's disease

High throughput RNA sequencing brings new perspective to the elucidation of molecular mechanisms of diseases. Normalization is the first and most important step for RNA-Seq data, and it can differ based on the purpose of the analysis. Within-sample normalization methods (eg. TPM) are preferred when genes in a sample are compared with each other, and between-sample normalization methods (eg. deseq2, TMM, Voom) are used when the samples in a dataset are compared. Normalization approaches rescale the data, and, therefore, they affect the results of the analysis. Here, we selected two most commonly used Alzheimer's disease RNA-Seq datasets from ROSMAP and Mayo Clinic cohorts and mapped the differentially expressed genes on human protein interactome to discover disease-specific subnetworks. To this end, the raw count data were first processed with four different, commonly used RNA-Seq normalization methods (deseq2, TMM, Voom and TPM). Then, covariate adjustment was applied to the normalized data for gender, age of death and post-mortem interval. Each normalized dataset was separately mapped on the human protein-protein interaction network either in covariate-adjusted or non-adjusted form. Capturing known Alzheimer's disease genes and genes associated with the disease-related functional terms in the discovered subnetworks were the criteria to compare different normalization methods. Based on our results, applying covariate adjustment has a positive effect on normalization by removing the confounder effects. Covariate-adjusted TMM and covariate-adjusted deseq2 methods performed better in both transcriptome datasets.

Molecular adaptations underlying high-frequency hearing in the brain of CF bats species

BACKGROUND

The majority of bat species have developed remarkable echolocation ability, especially for the laryngeally echolocating bats along with high-frequency hearing. Adaptive evolution has been widely detected for the cochleae in the laryngeally echolocating bats, however, limited understanding for the brain which is the central to echolocation signal processing in the auditory perception system, the laryngeally echolocating bats brain may also undergo adaptive changes.

RESULT

In order to uncover the molecular adaptations related with high-frequency hearing in the brain of laryngeally echolocating bats, the genes expressed in the brain of Rhinolophus ferrumequinum (CF bat) and Myotis pilosus (FM bat) were both detected and also compared. A total of 346,891 genes were detected and the signal transduction mechanisms were annotated by the most abundant genes, followed by the transcription. In hence, there were 3,088 DEGs were found between the two bat brains, with 1,426 highly expressed in the brain of R. ferrumequinum, which were significantly enriched in the neuron and neurodevelopmental processes. Moreover, we found a key candidate hearing gene, ADCY1, playing an important role in the R. ferrumequinum brain and undergoing adaptive evolution in CF bats.

CONCLUSIONS

Our study provides a new insight to the molecular bases of high-frequency hearing in two laryngeally echolocating bats brain and revealed different nervous system activities during auditory perception in the brain of CF bats.

Hybridization and gene expression: Beyond differentially expressed genes

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

Integration of Physiological, Transcriptomic, and Metabolomic Analyses Reveal Molecular Mechanisms of Salt Stress in Maclura tricuspidata

Salt stress is a universal abiotic stress that severely affects plant growth and development. Understanding the mechanisms of Maclura tricuspidate's adaptation to salt stress is crucial for developing salt-tolerant plant varieties. This article discusses the integration of physiology, transcriptome, and metabolome to investigate the mechanism of salt adaptation in M. tricuspidata under salt stress conditions. Overall, the antioxidant enzyme system (SOD and POD) of M. tricuspidata exhibited higher activities compared with the control, while the content of soluble sugar and concentrations of chlorophyll a and b were maintained during salt stress. KEGG analysis revealed that deferentially expressed genes were primarily involved in plant hormone signal transduction, phenylpropanoid and flavonoid biosynthesis, alkaloids, and MAPK signaling pathways. Differential metabolites were enriched in amino acid metabolism, the biosynthesis of plant hormones, butanoate, and 2-oxocarboxylic acid metabolism. Interestingly, glycine, serine, and threonine metabolism were found to be important both in the metabolome and transcriptome-metabolome correlation analyses, suggesting their essential role in enhancing the salt tolerance of M. tricuspidata. Collectively, our study not only revealed the molecular mechanism of salt tolerance in M. tricuspidata, but also provided a new perspective for future salt-tolerant breeding and improvement in salt land for this species.

Two avian Plasmodium species trigger different transcriptional responses on their vector Culex pipiens

Malaria is a mosquito-borne disease caused by protozoans of the genus Plasmodium that affects both humans and wildlife. The fitness consequences of infections by avian malaria are well known in birds, however, little information exists on its impact on mosquitoes. Here we study how Culex pipiens mosquitoes transcriptionally respond to infection by two different Plasmodium species, P. relictum and P. cathemerium, differing in their virulence (mortality rate) and transmissibility (parasite presence in exposed mosquitoes' saliva). We studied the mosquito response to the infection at three critical stages of parasite development: the formation of ookinetes at 24 h post-infection (hpi), the release of sporozoites into the hemocoel at 10 days post-infection (dpi), and the storage of sporozoites in the salivary glands at 21 dpi. For each time point, we characterized the gene expression of mosquitoes infected with each P. relictum and P. cathemerium and mosquitoes fed on an uninfected bird and, subsequently, compared their transcriptomic responses. Differential gene expression analysis showed that most transcriptomic changes occurred during the early infection stage (24 hpi), especially when comparing P. relictum and P. cathemerium-infected mosquitoes. Differentially expressed genes in mosquitoes infected with each species were related mainly to the metabolism of the immune response, trypsin, and other serine-proteases. We conclude that these differences in response may partly play a role in the differential virulence and transmissibility previously observed between P. relictum and P. cathemerium in Cx. pipiens.

Antagonistic Activity of Streptomyces alfalfae 11F against Fusarium Wilt of Watermelon and Transcriptome Analysis Provides Insights into the Synthesis of Phenazine-1-Carboxamide

Streptomyces alfalfa strain 11F has inhibitory effects on many phytopathogenic fungi and improves the establishment and biomass yield of switchgrass. However, the antagonistic effects of strain 11F on Fusarium wilt of watermelon and its secondary metabolites that contribute to its biocontrol activity are poorly understood. We evaluated the antagonistic and growth-promoting effects of strain 11F and conducted a transcriptome analysis to identify the metabolites contributing to antifungal activity. Strain 11F had marked inhibitory effects on six fungal pathogens. The incidence of Fusarium wilt of watermelon seedlings was decreased by 46.02%, while watermelon seedling growth was promoted, as indicated by plant height (8.7%), fresh weight (23.1%), and dry weight (60.0%). Clean RNA-sequencing data were annotated with 7553 functional genes. The 2582 differentially expressed genes (DEGs) detected in the Control vs. Case 2 comparison were divided into 42 subcategories of the biological process, cellular component, and molecular function Gene Ontology categories. Seven hundred and forty functional genes (55.47% of the DEGs) were assigned to Kyoto Encyclopedia of Genes and Genomes metabolic pathways, reflecting the complexity of the strain 11F metabolic regulatory system. The expression level of the gene phzF, which encodes an enzyme essential for phenazine-1-carboxylic acid (PCA) synthesis, was downregulated 3.7-fold between the 24 h and 48 h fermentation time points, suggesting that strain 11F can produce phenazine compounds. A phenazine compound from 11F was isolated and identified as phenazine-1-carboxamide (PCN), which contributed to the antagonistic activity against Fusarium oxysporum f. sp. niveum. PCA was speculated to be the synthetic precursor of PCN. The downregulation in phzF expression might be associated with the decrease in PCA accumulation and the increase in PCN synthesis in strain 11F from 24 to 48 h. Streptomyces alfalfae 11F protects watermelon seedlings from Fusarium wilt of watermelon and promotes seedling growth. The transcriptome analysis of strain 11F provides insights into the synthesis of PCN, which has antifungal activity against F. oxysporum f. sp. niveum of watermelon.

Potential Effects of High Temperature and Heat Wave on Nanorana pleskei Based on Transcriptomic Analysis

In the context of climate change, understanding how indigenous amphibians of the Qinghai-Tibet plateau react to stresses and their coping mechanisms could be crucial for predicting their fate and successful conservation. A liver transcriptome for Nanorana pleskei was constructed using high-throughput RNA sequencing, and its gene expression was compared with frogs acclimated under either room temperature or high temperature and also heat wave exposed ones. A total of 126,465 unigenes were produced, with 66,924 (52.92%) of them being annotated. Up to 694 genes were found to be differently regulated as a result of abnormal temperature acclimatization. Notably, genes belonging to the heat shock protein (HSP) family were down-regulated in both treated groups. Long-term exposure to high-temperature stress may impair the metabolic rate of the frog and trigger the body to maintain a hypometabolic state in an effort to survive challenging times. During heat waves, unlike the high-temperature group, mitochondrial function was not impaired, and the energy supply was largely normal to support the highly energy-consuming metabolic processes. Genes were more transcriptionally suppressed when treated with high temperatures than heat waves, and the body stayed in low-energy states for combating these long-term adverse environments to survive. It might be strategic to preserve initiation to executive protein activity under heat wave stress. Under both stress conditions, compromising the protection of HSP and sluggish steroid activity occurred in frogs. Frogs were more affected by high temperatures than by heat waves.

Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3' RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

Optimization and application of non-native Phragmites australis transcriptome assemblies

Phragmites australis (common reed) has a cosmopolitan distribution and has been suggested as a model organism for the study of invasive plant species. In North America, the non-native subspecies (ssp. australis) is widely distributed across the contiguous 48 states in the United States and large parts of Canada. Even though millions of dollars are spent annually on Phragmites management, insufficient knowledge of P. australis impeded the efficiency of management. To solve this problem, transcriptomic information generated from multiple types of tissue could be a valuable resource for future studies. Here, we constructed forty-nine P. australis transcriptomes assemblies via different assembly tools and multiple parameter settings. The optimal transcriptome assembly for functional annotation and downstream analyses was selected among these transcriptome assemblies by comprehensive assessments. For a total of 422,589 transcripts assembled in this transcriptome assembly, 319,046 transcripts (75.5%) have at least one functional annotation. Within the transcriptome assembly, we further identified 1,495 transcripts showing tissue-specific expression pattern, 10,828 putative transcription factors, and 72,165 candidates for simple sequence repeats markers. The identification and analyses of predicted transcripts related to herbicide- and salinity-resistant genes were shown as two applications of the transcriptomic information to facilitate further research on P. australis. Transcriptome assembly and selection would be important for the transcriptome annotation. With this optimal transcriptome assembly and all relative information from downstream analyses, we have helped to establish foundations for future studies on the mechanisms underlying the invasiveness of non-native P. australis subspecies.

Co-expression Gene Networks and Machine-learning Algorithms Unveil a Core Genetic Toolkit for Reproductive Division of Labour in Rudimentary Insect Societies

The evolution of eusociality requires that individuals forgo some or all their own reproduction to assist the reproduction of others in their group, such as a primary egg-laying queen. A major open question is how genes and genetic pathways sculpt the evolution of eusociality, especially in rudimentary forms of sociality-those with smaller cooperative nests when compared with species such as honeybees that possess large societies. We lack comprehensive comparative studies examining shared patterns and processes across multiple social lineages. Here we examine the mechanisms of molecular convergence across two lineages of bees and wasps exhibiting such rudimentary societies. These societies consist of few individuals and their life histories range from facultative to obligately social. Using six species across four independent origins of sociality, we conduct a comparative meta-analysis of publicly available transcriptomes. Standard methods detected little similarity in patterns of differential gene expression in brain transcriptomes among reproductive and non-reproductive individuals across species. By contrast, both supervised machine learning and consensus co-expression network approaches uncovered sets of genes with conserved expression patterns among reproductive and non-reproductive phenotypes across species. These sets overlap substantially, and may comprise a shared genetic "toolkit" for sociality across the distantly related taxa of bees and wasps and independently evolved lineages of sociality. We also found many lineage-specific genes and co-expression modules associated with social phenotypes and possible signatures of shared life-history traits. These results reveal how taxon-specific molecular mechanisms complement a core toolkit of molecular processes in sculpting traits related to the evolution of eusociality.

Histological and MRI brain atlas of the common shrew, Sorex araneus, with brain region-specific gene expression profiles

The common shrew, Sorex araneus, is a small mammal of growing interest in neuroscience research, as it exhibits dramatic and reversible seasonal changes in individual brain size and organization (a process known as Dehnel's phenomenon). Despite decades of studies on this system, the mechanisms behind the structural changes during Dehnel's phenomenon are not yet understood. To resolve these questions and foster research on this unique species, we present the first combined histological, magnetic resonance imaging (MRI), and transcriptomic atlas of the common shrew brain. Our integrated morphometric brain atlas provides easily obtainable and comparable anatomic structures, while transcriptomic mapping identified distinct expression profiles across most brain regions. These results suggest that high-resolution morphological and genetic research is pivotal for elucidating the mechanisms underlying Dehnel's phenomenon while providing a communal resource for continued research on a model of natural mammalian regeneration. Morphometric and NCBI Sequencing Read Archive are available at https://doi.org/10.17617/3.HVW8ZN.

Metformin modulates the gut microbiome in a mice model of high-fat diet-induced glycolipid metabolism disorder

INTRODUCTION

Metformin (MET) can regulate glucose and lipid levels, and the gut microbiota may be involved in the control of metabolism. We hypothesized that MET alleviates glucolipid metabolism disorder by modulating gut microbiota and microbial metabolites.

RESEARCH DESIGN AND METHODS

A total of 24 male C57BL/6 J mice were equally divided into three groups (normal control, model control (MC), and MET-treated groups). Model mice were established by feeding a high-fat diet for 6 weeks. The MET-treated group was administered MET solution (2.5 g/100 mL, 250 mg/kg). Fecal samples were collected to characterize the microbiota system using metagenomic shotgun sequencing and gas chromatography-time of flight-mass spectrometry analysis. Phenotypic and biochemical indices were obtained for further correlation analysis.

RESULTS

Compared with the MC group, MET reduced the levels of weight, glucose, areas under the glucose curve in the glucose tolerance test, triglyceride (TG), and total cholesterol (TC). A decreasing abundance of bacteria, including Parabacteroides distasonis, and an increasing abundance of bacteria, including Bacteroides vulgatus, were observed in the MET-treated group. The 2-deoxytetronic acid declined after MET intervention and was positively correlated with species over-represented in the MC group and negatively correlated with species enriched in the MET-treated group. Additionally, species enriched in the MET-treated group negatively correlated with glucose, areas under the glucose curve in the glucose tolerance test, and TGs. Further, the correlation between the differential metabolites, which decreased after MET intervention, and the phenotypic indices was positive.

CONCLUSIONS

MET-induced restoration of intestinal homeostasis correlates with the amelioration of host glucolipid metabolism.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

Transcriptomic analysis of pre-ovipositional embryonic arrest in a non-squamate reptile (Chelonia mydas)

After gastrulation, oviductal hypoxia maintains turtle embryos in an arrested state prior to oviposition. Subsequent exposure to atmospheric oxygen upon oviposition initiates recommencement of embryonic development. Arrest can be artificially extended for several days after oviposition by incubation of the egg under hypoxic conditions, with development recommencing in an apparently normal fashion after subsequent exposure to normoxia. To examine the transcriptomic events associated with embryonic arrest in green sea turtles (Chelonia mydas), RNA‐sequencing analysis was performed on embryos from freshly laid eggs and eggs incubated in either normoxia (oxygen tension ~159 mmHg) or hypoxia (<8 mmHg) for 36 h after oviposition (n = 5 per group). The patterns of gene expression differed markedly among the three experimental groups. Normal embryonic development in normoxia was associated with upregulation of genes involved in DNA replication, the cell cycle, and mitosis, but these genes were commonly downregulated after incubation in hypoxia. Many target genes of hypoxia inducible factors, including the gene encoding insulin‐like growth factor binding protein 1 (igfbp1), were downregulated by normoxic incubation but upregulated by incubation in hypoxia. Notably, some of the transcriptomic effects of hypoxia in green turtle embryos resembled those reported to be associated with hypoxia‐induced embryonic arrest in diverse taxa, including the nematode Caenorhabditis elegans and zebrafish (Danio rerio). Hypoxia‐induced preovipositional embryonic arrest appears to be a unique adaptation of turtles. However, our findings accord with the proposition that the mechanisms underlying hypoxia‐induced embryonic arrest per se are highly conserved across diverse taxa.

Genomics and Informatics, Conjoined Tools Vital for Understanding and Protecting Plant Health

Genomics' impact on crop production continuously expands. The number of sequenced plant and microbial species and strains representing diverse populations of individual species rapidly increases thanks to the advent of next-generation sequencing technologies. Their genomic blueprints revealed candidate genes involved in various functions and processes crucial for crop health and helped in understanding how the sequenced organisms have evolved at the genome level. Functional genomics quickly translates these blueprints into a detailed mechanistic understanding of how such functions and processes work and are regulated; this understanding guides and empowers efforts to protect crops from diverse biotic and abiotic threats. Metagenome analyses help identify candidate microbes crucial for crop health and uncover how microbial communities associated with crop production respond to environmental conditions and cultural practices, presenting opportunities to enhance crop health by judiciously configuring microbial communities. Efficient conversion of disparate types of massive genomics data into actionable knowledge requires a robust informatics infrastructure supporting data preservation, analysis, and sharing. This review starts with an overview of how genomics came about and has quickly transformed life science. We illuminate how genomics and informatics can be applied to investigate various crop health-related problems using selected studies. We end the review by noting why community empowerment via crowdsourcing is crucial to harnessing genomics to protect global food and nutrition security without continuously expanding the environmental footprint of crop production.

Genomic reaction norms inform predictions of plastic and adaptive responses to climate change

Genomic reaction norms represent the range of gene expression phenotypes (usually mRNA transcript levels) expressed by a genotype along an environmental gradient. Reaction norms derived from common‐garden experiments are powerful approaches for disentangling plastic and adaptive responses to environmental change in natural populations. By treating gene expression as a phenotype in itself, genomic reaction norms represent invaluable tools for exploring causal mechanisms underlying organismal responses to climate change across multiple levels of biodiversity.

Our goal is to provide the context, framework and motivation for applying genomic reaction norms to study the responses of natural populations to climate change.

Here, we describe the utility of integrating genomics with common‐garden‐gradient experiments under a reaction norm analytical framework to answer fundamental questions about phenotypic plasticity, local adaptation, their interaction (i.e. genetic variation in plasticity) and future adaptive potential.

An experimental and analytical framework for constructing and analysing genomic reaction norms is presented within the context of polygenic climate change responses of structured populations with gene flow. Intended for a broad eco‐evo readership, we first briefly review adaptation with gene flow and the importance of understanding the genomic basis and spatial scale of adaptation for conservation and management of structured populations under anthropogenic change. Then, within a high‐dimensional reaction norm framework, we illustrate how to distinguish plastic, differentially expressed (difference in reaction norm intercepts) and differentially plastic (difference in reaction norm slopes) genes, highlighting the areas of opportunity for applying these concepts.

We conclude by discussing how genomic reaction norms can be incorporated into a holistic framework to understand the eco‐evolutionary dynamics of climate change responses from molecules to ecosystems. We aim to inspire researchers to integrate gene expression measurements into common‐garden experimental designs to investigate the genomics of climate change responses as sequencing costs become increasingly accessible.

Unraveling patterns of disrupted gene expression across a complex tissue

Whole tissue RNASeq is the standard approach for studying gene expression divergence in evolutionary biology and provides a snapshot of the comprehensive transcriptome for a given tissue. However, whole tissues consist of diverse cell types differing in expression profiles, and the cellular composition of these tissues can evolve across species. Here, we investigate the effects of different cellular composition on whole tissue expression profiles. We compared gene expression from whole testes and enriched spermatogenesis populations in two species of house mice, Mus musculus musculus and M. m. domesticus, and their sterile and fertile F1 hybrids, which differ in both cellular composition and regulatory dynamics. We found that cellular composition differences skewed expression profiles and differential gene expression in whole testes samples. Importantly, both approaches were able to detect large-scale patterns such as disrupted X chromosome expression, although whole testes sampling resulted in decreased power to detect differentially expressed genes. We encourage researchers to account for histology in RNASeq and consider methods that reduce sample complexity whenever feasible. Ultimately, we show that differences in cellular composition between tissues can modify expression profiles, potentially altering inferred gene ontological processes, insights into gene network evolution, and processes governing gene expression evolution.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.

Athermal, Chemically Triggered Release of RNA from Thioester Nucleic Acids

An athermal approach to mRNA enrichment from total RNA using a self-immolative thioester linked nucleic acids (TENA) is described. Oligo(thymine) (oT) TENA has a six-atom spacing between bases which allowed TENA to selectively base-pair with polyadenine RNA. As a result of the neutral backbone of TENA and the hydrophobicity of the octanethiol end group, oT TENA is water insoluble and efficiently pulled down 93±2 % of EGFP mRNA at a concentration of 10 ng μL^-1 . Self-immolative degradation of TENA upon ambient temperature exposure to nucleophilic buffer components (Tris, DTT) allowed recovery of 55±27 ng of mRNA from 3.1 μg of total RNA, which was not statistically different from the amount recovered using Dynabeads® mRNA DIRECT Kit (89±24 ng). Gene expression as measured by RT-qPCR was comparable for both enrichment methods, suggesting that the mild conditions required for enrichment of mRNA using oT TENA are compatible with RT-qPCR and other downstream molecular biology applications.

Prediction of adverse effects of effluents containing phenolic compounds in the Ba River on the ovary of fish (Hemiculter leucisculus) using transcriptomic and metabolomic analyses

The aim of this work was to evaluate the endocrine disrupting effects on the ovarian development of sharpbelly (Hemiculter leucisculus) caused by effluents containing phenolic compounds. This was achieved using integrated transcriptomic and metabolomic analyses, along with histopathological examinations. Sharpbelly, an indigenous freshwater fish widely distributed in East Asia, were collected by pole fishing from three sampling sites in the Ba River. These sampling sites include a mid-stream site near a wastewater outfall and a reference site located upstream and a far field comparison site located downstream. In sharpbelly collected near the wastewater discharge, the oocyte development was activated, compared to the other two sites. Histopathological alterations in the fish ovaries were likely due to the upregulated steroid hormone biosynthesis process, as suggested by the differentially expressed genes (e.g., hsd3b, hsd17b1) and differentially accumulated metabolites (e.g., pregnenolone). Additionally, under the stress of effluents containing phenolic compounds, genes related to the signaling pathways for oxidative phosphorylation and leukocyte transendothelial migration were dysregulated, suggesting the potential induction of inflammation and several ovarian diseases. Overall, these findings suggest that effluents containing phenolic compounds influence ovary development and reproductive function of female sharpbelly. Whether there is any resulting dysfunction of folliculogenesis, abnormality of ovulation, production of premature eggs and/or potential induction of ovarian cancers remains to be determined by further studies, for a better evaluation on effluents containing phenolic compounds to the fish fertility and the health of their offspring, and even the stability of the wild fish population. Notably, the integration of transcriptomics and metabolomics can complement the routine chemical analysis to comprehensively monitor the effects of wastewater treatment plant effluents on the health of wild fish.

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 novel TFL1 gene induces flowering in the mast seeding alpine snow tussock, Chionochloa pallens (Poaceae)

Masting, the synchronous, highly variable flowering across years by a population of perennial plants, has been reported to be precipitated by various factors including nitrogen levels, drought conditions, and spring and summer temperatures. However, the molecular mechanism leading to the initiation of flowering in masting plants in particular years remains largely unknown, despite the potential impact of climate change on masting phenology. We studied genes controlling flowering in the alpine snow tussock Chionochloa pallens (Poaceae), a strongly masting perennial grass. We used a range of in situ and manipulated plants to obtain leaf samples from tillers (shoots) which subsequently remained vegetative or flowered. Here, we show that a novel orthologue of TERMINAL FLOWER 1 (TFL1; normally a repressor of flowering in other species) promotes the induction of flowering in C. pallens (hence Anti-TFL1), a conclusion supported by structural, functional and expression analyses. Global transcriptomic analysis indicated differential expression of CpTPS1, CpGA20ox1, CpREF6 and CpHDA6, emphasizing the role of endogenous cues and epigenetic regulation in terms of responsiveness of plants to initiate flowering. Our molecular-based study provides insights into the cellular mechanism of flowering in masting plants and will supplement ecological and statistical models to predict how masting will respond to global climate change.

Evolutionary history of two rare endemic conifer species from the eastern Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Understanding the population genetics and evolutionary history of endangered species is urgently needed in an era of accelerated biodiversity loss. This knowledge is most important for regions with high endemism that are ecologically vulnerable, such as the Qinghai-Tibet Plateau (QTP).

METHODS

The genetic variation of 84 juniper trees from six populations of Juniperus microsperma and one population of Juniperus erectopatens, two narrow-endemic junipers from the QTP that are sister to each other, was surveyed using RNA-sequencing data. Coalescent-based analyses were used to test speciation, migration and demographic scenarios. Furthermore, positively selected and climate-associated genes were identified, and the genetic load was assessed for both species.

KEY RESULTS

Analyses of 149 052 single nucleotide polymorphisms showed that the two species are well differentiated and monophyletic. They diverged around the late Pliocene, but interspecific gene flow continued until the Last Glacial Maximum. Demographic reconstruction by Stairway Plot detected two severe bottlenecks for J. microsperma but only one for J. erectopatens. The identified positively selected genes and climate-associated genes revealed habitat adaptation of the two species. Furthermore, although J. microsperma had a much wider geographical distribution than J. erectopatens, the former possesses lower genetic diversity and a higher genetic load than the latter.

CONCLUSIONS

This study sheds light on the evolution of two endemic juniper species from the QTP and their responses to Quaternary climate fluctuations. Our findings emphasize the importance of speciation and demographic history reconstructions in understanding the current distribution pattern and genetic diversity of threatened species in mountainous regions.

Admixture mapping reveals loci for carcass mass in red deer x sika hybrids in Kintyre, Scotland

We deployed admixture mapping on a sample of 386 deer from a hybrid swarm between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon) sampled in Kintyre, Scotland to search for quantitative trait loci (QTLs) underpinning phenotypic differences between the species. These two species are highly diverged genetically [Fst between pure species, based on 50K single nucleotide polymorphism (SNPs) = 0.532] and phenotypically: pure red have on average twice the carcass mass of pure sika in our sample (38.7 kg vs 19.1 kg). After controlling for sex, age, and population genetic structure, we found 10 autosomal genomic locations with QTL for carcass mass. Effect sizes ranged from 0.191 to 1.839 kg and as expected, in all cases the allele derived from sika conferred lower carcass mass. The sika population was fixed for all small carcass mass alleles, whereas the red deer population was typically polymorphic. GO term analysis of genes lying in the QTL regions are associated with oxygen transport. Although body mass is a likely target of selection, none of the SNPs marking QTL are introgressing faster or slower than expected in either direction.

Can functional genomic diversity provide novel insights into mechanisms of community assembly? A pilot study from an invaded alpine streambed

An important focus of community ecology, including invasion biology, is to investigate functional trait diversity patterns to disentangle the effects of environmental and biotic interactions. However, a notable limitation is that studies usually rely on a small and easy-to-measure set of functional traits, which might not immediately reflect ongoing ecological responses to changing abiotic or biotic conditions, including those that occur at a molecular or physiological level. We explored the potential of using the diversity of expressed genes-functional genomic diversity (FGD)-to understand ecological dynamics of a recent and ongoing alpine invasion. We quantified FGD based on transcriptomic data measured for 26 plant species occurring along adjacent invaded and pristine streambeds. We used an RNA-seq approach to summarize the overall number of expressed transcripts and their annotations to functional categories, and contrasted this with functional trait diversity (FTD) measured from a suite of characters that have been traditionally considered in plant ecology. We found greater FGD and FTD in the invaded community, independent of differences in species richness. However, the magnitude of functional dispersion was greater from the perspective of FGD than from FTD. Comparing FGD between congeneric alien-native species pairs, we did not find many significant differences in the proportion of genes whose annotations matched functional categories. Still, native species with a greater relative abundance in the invaded community compared with the pristine tended to express a greater fraction of genes at significant levels in the invaded community, suggesting that changes in FGD may relate to shifts in community composition. Comparisons of diversity patterns from the community to the species level offer complementary insights into processes and mechanisms driving invasion dynamics. FGD has the potential to illuminate cryptic changes in ecological diversity, and we foresee promising avenues for future extensions across taxonomic levels and macro-ecosystems.

Dietary Diversification and Specialization in Neotropical Bats Facilitated by Early Molecular Evolution

Dietary adaptation is a major feature of phenotypic and ecological diversification, yet the genetic basis of dietary shifts is poorly understood. Among mammals, Neotropical leaf-nosed bats (family Phyllostomidae) show unmatched diversity in diet; from a putative insectivorous ancestor, phyllostomids have radiated to specialize on diverse food sources including blood, nectar, and fruit. To assess whether dietary diversification in this group was accompanied by molecular adaptations for changing metabolic demands, we sequenced 89 transcriptomes across 58 species and combined these with published data to compare ∼13,000 protein coding genes across 66 species. We tested for positive selection on focal lineages, including those inferred to have undergone dietary shifts. Unexpectedly, we found a broad signature of positive selection in the ancestral phyllostomid branch, spanning genes implicated in the metabolism of all major macronutrients, yet few positively selected genes at the inferred switch to plantivory. Branches corresponding to blood- and nectar-based diets showed selection in loci underpinning nitrogenous waste excretion and glycolysis, respectively. Intriguingly, patterns of selection in metabolism genes were mirrored by those in loci implicated in craniofacial remodeling, a trait previously linked to phyllostomid dietary specialization. Finally, we show that the null model of the widely-used branch-site test is likely to be misspecified, with the implication that the test is too conservative and probably under-reports true cases of positive selection. Our findings point to a complex picture of adaptive radiation, in which the evolution of new dietary specializations has been facilitated by early adaptations combined with the generation of new genetic variation.

Transcriptomes of Saussurea (Asteraceae) Provide Insights into High-Altitude Adaptation

Understanding how species adapt to extreme environments is an extension of the main goals of evolutionary biology. While alpine plants are an ideal system for investigating the genetic basis of high-altitude adaptation, genomic resources in these species are still limited. In the present study, we generated reference-level transcriptomic data of five Saussurea species through high-throughput sequencing and de novo assembly. Three of them are located in the highland of the Qinghai-Tibet Plateau (QTP), and the other two are close relatives distributed in the lowland. A series of comparative and evolutionary genomics analyses were conducted to explore the genetic signatures of adaptive evolution to high-altitude environments. Estimation of divergence time using single-copy orthologs revealed that Saussurea species diversified during the Miocene, a period with extensive tectonic movement and climatic fluctuation on the QTP. We characterized gene families specific to the alpine species, including genes involved in oxidoreductase activity, pectin catabolic process, lipid transport, and polysaccharide metabolic process, which may play important roles in defense of hypoxia and freezing temperatures of the QTP. Furthermore, in a phylogenetic context with the branch model, we identified hundreds of genes with signatures of positive selection. These genes are involved in DNA repair, membrane transport, response to UV-B and hypoxia, and reproductive processes, as well as some metabolic processes associated with nutrient intake, potentially responsible for Saussurea adaptation to the harsh environments of high altitude. Overall, our study provides valuable genomic resources for alpine species and gained helpful insights into the genomic basis of plants adapting to extreme environments.

Comparative Analyses of Gene Co-expression Networks: Implementations and Applications in the Study of Evolution

Similarities and differences in the associations of biological entities among species can provide us with a better understanding of evolutionary relationships. Often the evolution of new phenotypes results from changes to interactions in pre-existing biological networks and comparing networks across species can identify evidence of conservation or adaptation. Gene co-expression networks (GCNs), constructed from high-throughput gene expression data, can be used to understand evolution and the rise of new phenotypes. The increasing abundance of gene expression data makes GCNs a valuable tool for the study of evolution in non-model organisms. In this paper, we cover motivations for why comparing these networks across species can be valuable for the study of evolution. We also review techniques for comparing GCNs in the context of evolution, including local and global methods of graph alignment. While some protein-protein interaction (PPI) bioinformatic methods can be used to compare co-expression networks, they often disregard highly relevant properties, including the existence of continuous and negative values for edge weights. Also, the lack of comparative datasets in non-model organisms has hindered the study of evolution using PPI networks. We also discuss limitations and challenges associated with cross-species comparison using GCNs, and provide suggestions for utilizing co-expression network alignments as an indispensable tool for evolutionary studies going forward.

Contributions of cis- and trans-Regulatory Evolution to Transcriptomic Divergence across Populations in the Drosophila mojavensis Larval Brain

Natural selection on gene expression was originally predicted to result primarily in cis- rather than trans-regulatory evolution, due to the expectation of reduced pleiotropy. Despite this, numerous studies have ascribed recent evolutionary divergence in gene expression predominantly to trans-regulation. Performing RNA-seq on single isofemale lines from genetically distinct populations of the cactophilic fly Drosophila mojavensis and their F₁ hybrids, we recapitulated this pattern in both larval brains and whole bodies. However, we demonstrate that improving the measurement of brain expression divergence between populations by using seven additional genotypes considerably reduces the estimate of trans-regulatory contributions to expression evolution. We argue that the finding of trans-regulatory predominance can result from biases due to environmental variation in expression or other sources of noise, and that cis-regulation is likely a greater contributor to transcriptional evolution across D. mojavensis populations. Lastly, we merge these lines of data to identify several previously hypothesized and intriguing novel candidate genes, and suggest that the integration of regulatory and population-level transcriptomic data can provide useful filters for the identification of potentially adaptive genes.

Multiomic analysis of the Arabian camel (Camelus dromedarius) kidney reveals a role for cholesterol in water conservation

The Arabian camel (Camelus dromedarius) is the most important livestock animal in arid and semi-arid regions and provides basic necessities to millions of people. In the current context of climate change, there is renewed interest in the mechanisms that enable camelids to survive in arid conditions. Recent investigations described genomic signatures revealing evolutionary adaptations to desert environments. We now present a comprehensive catalogue of the transcriptomes and proteomes of the dromedary kidney and describe how gene expression is modulated as a consequence of chronic dehydration and acute rehydration. Our analyses suggested an enrichment of the cholesterol biosynthetic process and an overrepresentation of categories related to ion transport. Thus, we further validated differentially expressed genes with known roles in water conservation which are affected by changes in cholesterol levels. Our datasets suggest that suppression of cholesterol biosynthesis may facilitate water retention in the kidney by indirectly facilitating the AQP2-mediated water reabsorption.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Big Data in Conservation Genomics: Boosting Skills, Hedging Bets, and Staying Current in the Field

A current challenge in the fields of evolutionary, ecological, and conservation genomics is balancing production of large-scale datasets with additional training often required to handle such datasets. Thus, there is an increasing need for conservation geneticists to continually learn and train to stay up-to-date through avenues such as symposia, meetings, and workshops. The ConGen meeting is a near-annual workshop that strives to guide participants in understanding population genetics principles, study design, data processing, analysis, interpretation, and applications to real-world conservation issues. Each year of ConGen gathers a diverse set of instructors, students, and resulting lectures, hands-on sessions, and discussions. Here, we summarize key lessons learned from the 2019 meeting and more recent updates to the field with a focus on big data in conservation genomics. First, we highlight classical and contemporary issues in study design that are especially relevant to working with big datasets, including the intricacies of data filtering. We next emphasize the importance of building analytical skills and simulating data, and how these skills have applications within and outside of conservation genetics careers. We also highlight recent technological advances and novel applications to conservation of wild populations. Finally, we provide data and recommendations to support ongoing efforts by ConGen organizers and instructors-and beyond-to increase participation of underrepresented minorities in conservation and eco-evolutionary sciences. The future success of conservation genetics requires both continual training in handling big data and a diverse group of people and approaches to tackle key issues, including the global biodiversity-loss crisis.

Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Resilience to climate change depends on a species' adaptive potential and phenotypic plasticity. The latter can enhance survival of individual organisms during short periods of extreme environmental perturbations, allowing genetic adaptation to take place over generations. Along the U.S. East Coast, estuarine-dependent spotted seatrout (Cynoscion nebulosus) populations span a steep temperature gradient that provides an ideal opportunity to explore the molecular basis of phenotypic plasticity. Genetically distinct spotted seatrout sampled from a northern and a southern population were exposed to acute cold and heat stress (5 biological replicates in each treatment and control group), and their transcriptomic responses were compared using RNA-sequencing (RNA-seq). The southern population showed a larger transcriptomic response to acute cold stress, whereas the northern population showed a larger transcriptomic response to acute heat stress compared with their respective population controls. Shared transcripts showing significant differences in expression levels were predominantly enriched in pathways that included metabolism, transcriptional regulation, and immune response. In response to heat stress, only the northern population significantly upregulated genes in the apoptosis pathway, which could suggest greater vulnerability to future heat waves in this population as compared to the southern population. Genes showing population-specific patterns of expression, including hpt, acot, hspa5, and hsc71, are candidates for future studies aiming to monitor intraspecific differences in temperature stress responses in spotted seatrout. Our findings contribute to the current understanding of phenotypic plasticity and provide a basis for predicting the response of a eurythermal fish species to future extreme temperatures.

Enhanced Biosynthesis of Chlorogenic Acid and Its Derivatives in Methyl-Jasmonate-Treated Gardenia jasminoides Cells: A Study on Metabolic and Transcriptional Responses of Cells

Chlorogenic acid and its derivatives (CQAs) are considered as important bioactive secondary metabolites in Gardenia jasminoides Ellis (G. jasminoides). However, few studies have investigated the biosynthesis and regulation of CQAs in G. jasminoides. In this study, methyl jasmonate (MeJA) was used to enhance CQAs accumulation in cultured G. jasminoides cells. Moreover, the possible molecular mechanism of MeJA-mediated accumulation of CQAs is also explored. To this end, time-course transcriptional profiles of G. jasminoides cells responding to MeJA were used to investigate the mechanism from different aspects, including jasmonate (JAs) biosynthesis, signal transduction, biosynthesis of precursor, CQAs biosynthesis, transporters, and transcription factors (TFs). A total of 57,069 unigenes were assembled from the clean reads, in which 80.7% unigenes were successfully annotated. Furthermore, comparative transcriptomic results indicated that differentially expressed genes (DEGs) were mainly involved in JAs biosynthesis and signal transduction (25 DEGs), biosynthesis of precursor for CQAs (18 DEGs), CQAs biosynthesis (19 DEGs), and transporters (29 DEGs). Most of these DEGs showed continuously upregulated expressions over time, which might activate the jasmonic acid (JA) signal transduction network, boost precursor supply, and ultimately stimulate CQAs biosynthesis. Additionally, various TFs from different TF families also responded to MeJA elicitation. Interestingly, 38 DEGs from different subgroups of the MYB family might display positive or negative regulations on phenylpropanoids, especially on CQAs biosynthesis. Conclusively, our results provide insight into the possible molecular mechanism of regulation on CQAs biosynthesis, which led to a high CQAs yield in the G. jasminoides cells under MeJA treatment.

Lipid metabolism in Calanus finmarchicus is sensitive to variations in predation risk and food availability

Late developmental stages of the marine copepods in the genus Calanus can spend extended periods in a dormant stage (diapause) that is preceded by the accumulation of large lipid stores. We assessed how lipid metabolism during development from the C4 stage to adult is altered in response to predation risk and varying food availability, to ultimately understand more of the metabolic processes during development in Calanus copepods. We used RNA sequencing to assess if perceived predation risk in combination with varied food availability affects expression of genes associated with lipid metabolism and diapause preparation in C. finmarchicus. The lipid metabolism response to predation risk differed depending on food availability, time and life stage. Predation risk caused upregulation of lipid catabolism with high food, and downregulation with low food. Under low food conditions, predation risk disrupted lipid accumulation. The copepods showed no clear signs of diapause preparation, supporting earlier observations of the importance of multiple environmental cues in inducing diapause in C. finmarchicus. This study demonstrates that lipid metabolism is a sensitive endpoint for the interacting environmental effects of predation pressure and food availability. As diapause may be controlled by lipid accumulation, our findings may contribute towards understanding processes that can ultimately influence diapause timing.

Messenger RNA enrichment using synthetic oligo(T) click nucleic acids

Enrichment of mRNA is a key step in a number of molecular biology techniques, particularly in the rapidly growing field of transcriptomics. Currently, mRNA is isolated using oligo(thymine) DNA (oligo(dT)) immobilized on solid supports, which binds to the poly(A) tail of mRNA to pull the mRNA out of solution through the use of magnets or centrifugal filters. Here, a simple method to isolate mRNA by complexing it with synthetic click nucleic acids (CNAs) is described. Oligo(T) CNA bound efficiently to mRNA, and because of the insolubility of CNA in water, >90% of mRNA was readily removed from solution using this method. Simple washing, buffer exchange, and heating steps enabled mRNA's enrichment from total RNA, with a yield of 3.1 ± 1.5% of the input total RNA by mass, comparable to the yield from commercially available mRNA enrichment beads. Further, the integrity and activity of mRNA after CNA-facilitated pulldown and release was evaluated through two assays. In vitro translation of EGFP mRNA confirmed the translatability of mRNA into functional protein and RT-qPCR was used to amplify enriched mRNA from total RNA extracts and compare gene expression to results obtained using commercially available products.