The utility of shallow RNA-Seq for documenting differential gene expression in genes with high and low levels of expression

Supplementation of the Plant Conditioner ELICE Vakcina® Product with β-Aminobutyric Acid and Salicylic Acid May Lead to Trans-Priming Signaling in Barley (Hordeum vulgare)

Plant immunological memory, priming, is a defense mechanism that can be triggered by external stimuli, leading to the activation of biochemical pathways and preparing plants for disease resistance. Plant conditioners improve yield and crop quality through nutrient efficiency and abiotic stress tolerance, which is enhanced by the addition of resistance- and priming-induced compounds. Based on this hypothesis, this study aimed to investigate plant responses to priming actives of different natures, including salicylic acid and beta-aminobutyric acid, in combination with the plant conditioning agent ELICE Vakcina^®. Phytotron experiments and RNA-Seq analyses of differentially expressed genes using the combinations of these three investigated compounds were performed in a barley culture to investigate possible synergistic relationships in the genetic regulatory network. The results indicated a strong regulation of defense responses, which was enhanced by supplemental treatments; however, both synergistic and antagonistic effects were enhanced with one or two components, depending on the supplementation. The overexpressed transcripts were functionally annotated to assess their involvement in jasmonic acid and salicylic acid signaling; however, their determinant genes were highly dependent on the supplemental treatments. Although the effects overlapped, the potential effects of trans-priming the two supplements tested could be largely separated.

Locomotion and brain gene expression exhibit sex-specific non-monotonic dose-response to HFPO-DA during Drosophila melanogaster lifespan

BACKGROUND

Legacy per- and polyfluoroalkyl substances (PFAS), known for their environmental persistence and bio-accumulative properties, have been phased out in the U.S. due to public health concerns. A newer polymerization aid used in the manufacture of some fluoropolymers, hexafluoropropylene oxide-dimer acid (HFPO-DA), has lower reported bioaccumulation and toxicity, but is a potential neurotoxicant implicated in dopaminergic neurodegeneration.

OBJECTIVE

We investigated HFPO-DA's bio-accumulative potential and sex-specific effects on lifespan, locomotion, and brain gene expression in fruit flies.

METHODS

We quantified bioaccumulation of HFPO-DA in fruit flies exposed to 8.7×10⁴µg/L of HFPO-DA in the fly media for 14 days via UHPLC-MS. Long-term effect on lifespan was determined by exposing both sexes to 8.7×10² - 8.7×10⁵µg/L of HFPO-DA in media. Locomotion was measured following 3, 7, and 14 days of exposures at 8.7×10¹ - 8.7×10⁵µg/L of HFPO-DA in media, and high-throughput 3'-end RNA-sequencing was used to quantify gene expression in fly brains across the same time points.

RESULTS

Bioaccumulation of HFPO-DA in fruit flies was not detected. HFPO-DA-induced effects on lifespan, locomotion, and brain gene expression, and lowest adverse effect level (LOAEL) showed sexually dimorphic patterns. Locomotion scores significantly decreased in at least one dose at all time points for females and only at 3-day exposure for males, while brain gene expression exhibited non-monotonic dose-response. Differentially expressed genes correlated to locomotion scores revealed sex-specific numbers of positively and negatively correlated genes per functional category.

CONCLUSION

Although HFPO-DA effects on locomotion and survival were significant at doses higher than the US EPA reference dose, the brain transcriptomic profiling reveals sex-specific changes and neurological molecular targets; gene enrichments highlight disproportionately affected categories, including immune response: female-specific co-upregulation suggests potential neuroinflammation. Consistent sex-specific exposure effects necessitate blocking for sex in experimental design during HFPO-DA risk assessment.

Toxicity assessment of hexafluoropropylene oxide-dimer acid on morphology, heart physiology, and gene expression during zebrafish (Danio rerio) development

Hexafluoropropylene oxide-dimer acid (HFPO-DA) is one of the emerging replacements for the "forever" carcinogenic and toxic long-chain PFAS. HFPO-DA is a polymerization aid used for manufacturing fluoropolymers, whose global distribution and undetermined toxic properties are a concern regarding human and ecological health. To assess embryotoxic potential, zebrafish embryos were exposed to HFPO-DA at concentrations of 0.5-20,000 mg/L at 24-, 48-, and 72-h post-fertilization (hpf). Heart rate increased significantly in embryos exposed to 2 mg/L and 10 mg/L HFPO-DA across all time points. Spinal deformities and edema phenotypes were evident among embryos exposed to 1000-16,000 mg/L HFPO-DA at 72 hpf. A median lethal concentration (LC₅₀) was derived as 7651 mg/L at 72 hpf. Shallow RNA sequencing analysis of 9465 transcripts identified 38 consistently differentially expressed genes at 0.5 mg/L, 1 mg/L, 2 mg/L, and 10 mg/L HFPO-DA exposures. Notably, seven downregulated genes were associated with visual response, and seven upregulated genes were expressed in or regulated the cardiovascular system. This study identifies biological targets and molecular pathways affected during animal development by an emerging, potentially problematic, and ubiquitous industrial chemical.

Natural immunity stimulation using ELICE16INDURES® plant conditioner in field culture of soybean

Recently, climate change has had an increasing impact on the world. Innate defense mechanisms operating in plants - such as PAMP-triggered Immunity (PTI) - help to reduce the adverse effects caused by various abiotic and biotic stressors. In this study, the effects of ELICE16INDURES® plant conditioner for organic farming, developed by the Research Institute for Medicinal Plants and Herbs Ltd. Budakalász Hungary, were studied in a soybean population in Northern Hungary. The active compounds and ingredients of this product were selected in such a way as to facilitate the triggering of general plant immunity without the presence and harmful effects of pathogens, thereby strengthening the healthy plant population and preparing it for possible stress effects. In practice, treatments of this agent were applied at two different time points and two concentrations. The conditioning effect was well demonstrated by using agro-drone and ENDVI determination in the soybean field. The genetic background of healthier plants was investigated by NGS sequencing, and by the expression levels of genes encoding enzymes involved in the catalysis of metabolic pathways regulating PTI. The genome-wide transcriptional profiling resulted in 13 contigs related to PAMP-triggered immunity and activated as a result of the treatments. Further analyses showed 16 additional PTI-related contigs whose gene expression changed positively as a result of the treatments. The gene expression values of genes encoded in these contigs were determined by in silico mRNA quantification and validated by RT-qPCR. Both - relatively low and high treatments - showed an increase in gene expression of key genes involving AOC, IFS, MAPK4, MEKK, and GST. Transcriptomic results indicated that the biosyntheses of jasmonic acid (JA), salicylic acid (SA), phenylpropanoid, flavonoid, phytoalexin, and cellular detoxification processes were triggered in the appropriate molecular steps and suggested that plant immune reactions may be activated also artificially, and innate immunity can be enhanced with proper plant biostimulants.

Physiological and transcriptomic effects of hexafluoropropylene oxide dimer acid in Caenorhabditis elegans during development

Per- and polyfluoroalkyl substances (PFAS) are chemicals resistant to degradation. While such a feature is desirable in consumer and industrial products, some PFAS, including perfluorooctanoic acid (PFOA), are toxic and bioaccumulate. Hexafluoropropylene oxide dimer acid (HFPO-DA), an emerging PFAS developed to replace PFOA, has not been extensively studied. To evaluate the potential toxicity of HFPO-DA with a cost- and time-efficient approach, we exposed C. elegans larvae for 48 h to 4 × 10^-9-4 g/L HFPO-DA in liquid media and measured developmental, behavioral, locomotor, and transcriptional effects at various exposure levels. Worms exposed to 1.5-4 g/L HFPO-DA were developmentally delayed, and progeny production was significantly delayed (p < 0.05) in worms exposed to 2-4 g/L HFPO-DA. Statistically significant differential gene expression was identified in all fourteen HFPO-DA exposure groups ranging from 1.25 × 10^-5 to 4 g/L, except for 6.25 × 10^-5 g/L. Among 10298 analyzed genes, 2624 differentially expressed genes (DEGs) were identified in the developmentally delayed 4 g/L group only, and 78 genes were differentially expressed in at least one of the thirteen groups testing 1.25 × 10^-5-2 g/L HFPO-DA exposures. Genes encoding for detoxification enzymes including cytochrome P450 and UDP glucuronosyltransferases were upregulated in 0.25-4 g/L acute exposure groups. DEGs were also identified in lower exposure level groups, though they did not share biological functions except for six ribosomal protein-coding genes. While our transcriptional data is inconclusive to infer mechanisms of toxicity, the significant gene expression differences at 1.25 × 10^-5 g/L, the lowest concentration tested for transcriptional changes, calls for further targeted analyses of low-dose HFPO-DA exposure effects.

Insights into the seasonal adaptive mechanisms of Chinese alligators (Alligator sinensis) from transcriptomic analyses

The Chinese alligator (Alligator sinensis) is an endemic and rare species in China, and is considered to be one of the most endangered vertebrates in the world. It is known to hibernate, an energy-saving strategy against cold temperatures and food deprivation. Changes in gene expression during hibernation remain largely unknown. To understand these complex seasonal adaptive mechanisms, we performed a comprehensive survey of differential gene expression in heart, skeletal muscle, and kidney of hibernating and active Chinese alligators using RNA-Sequencing. In total, we identified 4780 genes differentially expressed between the active and hibernating periods. GO and KEGG pathway analysis indicated the likely role of these differentially expressed genes (DEGs). The upregulated DEGs in the active Chinese alligator, CSRP3, MYG and PCKGC, may maintain heart and skeletal muscle contraction, transport and storage of oxygen, and enhance the body’s metabolism, respectively. The upregulated DEGs in the dormant Chinese alligator, ADIPO, CIRBP and TMM27, may improve insulin sensitivity and glucose/lipid metabolism, protect cells against harmful effects of cold temperature and hypoxia, regulate amino acid transport and uptake, and stimulate the proliferation of islet cells and the secretion of insulin. These results provide a foundation for understanding the molecular mechanisms of the seasonal adaptation required for hibernation in Chinese alligators, as well as effective information for other non-model organisms research.

Gene duplication and the evolution of phenotypic diversity in insect societies

Gene duplication is an important evolutionary process thought to facilitate the evolution of phenotypic diversity. We investigated if gene duplication was associated with the evolution of phenotypic differences in a highly social insect, the honeybee Apis mellifera. We hypothesized that the genetic redundancy provided by gene duplication could promote the evolution of social and sexual phenotypes associated with advanced societies. We found a positive correlation between sociality and rate of gene duplications across the Apoidea, indicating that gene duplication may be associated with sociality. We also discovered that genes showing biased expression between A. mellifera alternative phenotypes tended to be found more frequently than expected among duplicated genes than singletons. Moreover, duplicated genes had higher levels of caste-, sex-, behavior-, and tissue-biased expression compared to singletons, as expected if gene duplication facilitated phenotypic differentiation. We also found that duplicated genes were maintained in the A. mellifera genome through the processes of conservation, neofunctionalization, and specialization, but not subfunctionalization. Overall, we conclude that gene duplication may have facilitated the evolution of social and sexual phenotypes, as well as tissue differentiation. Thus this study further supports the idea that gene duplication allows species to evolve an increased range of phenotypic diversity.

Virus and dsRNA-triggered transcriptional responses reveal key components of honey bee antiviral defense

Recent high annual losses of honey bee colonies are associated with many factors, including RNA virus infections. Honey bee antiviral responses include RNA interference and immune pathway activation, but their relative roles in antiviral defense are not well understood. To better characterize the mechanism(s) of honey bee antiviral defense, bees were infected with a model virus in the presence or absence of dsRNA, a virus associated molecular pattern. Regardless of sequence specificity, dsRNA reduced virus abundance. We utilized next generation sequencing to examine transcriptional responses triggered by virus and dsRNA at three time-points post-infection. Hundreds of genes exhibited differential expression in response to co-treatment of dsRNA and virus. Virus-infected bees had greater expression of genes involved in RNAi, Toll, Imd, and JAK-STAT pathways, but the majority of differentially expressed genes are not well characterized. To confirm the virus limiting role of two genes, including the well-characterized gene, dicer, and a probable uncharacterized cyclin dependent kinase in honey bees, we utilized RNAi to reduce their expression in vivo and determined that virus abundance increased, supporting their involvement in antiviral defense. Together, these results further our understanding of honey bee antiviral defense, particularly the role of a non-sequence specific dsRNA-mediated antiviral pathway.

An integrated RNAseq-1H NMR metabolomics approach to understand soybean primary metabolism regulation in response to Rhizoctonia foliar blight disease

BACKGROUND

Rhizoctonia solani AG1-IA is a devastating phytopathogen causing Rhizoctonia foliar blight (RFB) of soybean worldwide with yield losses reaching 60%. Plant defense mechanisms are complex and information from different metabolic pathways is required to thoroughly understand plant defense regulation and function. Combining information from different "omics" levels such as transcriptomics, metabolomics, and proteomics is required to gain insights into plant metabolism and its regulation. As such, we studied fluctuations in soybean metabolism in response to R. solani infection at early and late disease stages using an integrated transcriptomics-metabolomics approach, focusing on the regulation of soybean primary metabolism and oxidative stress tolerance.

RESULTS

Transcriptomics (RNAseq) and metabolomics (1H NMR) data were analyzed individually and by integration using bidirectional orthogonal projections to latent structures (O2PLS) to reveal possible links between the metabolome and transcriptome during early and late infection stages. O2PLS analysis detected 516 significant transcripts, double that reported in the univariate analysis, and more significant metabolites than detected in partial least squares discriminant analysis. Strong separation of treatments based on integration of the metabolomes and transcriptomes of the analyzed soybean leaves was revealed, similar trends as those seen in analyses done on individual datasets, validating the integration method being applied. Strong fluctuations of soybean primary metabolism occurred in glycolysis, the TCA cycle, photosynthesis and photosynthates in response to R. solani infection. Data were validated using quantitative real-time PCR on a set of specific markers as well as randomly selected genes. Significant increases in transcript and metabolite levels involved in redox reactions and ROS signaling, such as peroxidases, thiamine, tocopherol, proline, L-alanine and GABA were also recorded. Levels of ethanol increased 24 h post-infection in soybean leaves, and alcohol dehydrogenase (ADH) loss-of-function mutants of Arabidopsis thaliana had higher necrosis than wild type plants.

CONCLUSIONS

As a proof-of-concept, this study offers novel insights into the biological correlations and identification of candidate genes and metabolites that can be used in soybean breeding for resistance to R. solani AG1-IA infection. Additionally, these findings imply that alcohol and its associated gene product ADH may have important roles in plant resistance to R. solani AG1-IA causing foliar blight.

The reference transcriptome of the adult female biting midge (Culicoides sonorensis) and differential gene expression profiling during teneral, blood, and sucrose feeding conditions

Unlike other important vectors such as mosquitoes and sandflies, genetic and genomic tools for Culicoides biting midges are lacking, despite the fact that they vector a large number of arboviruses and other pathogens impacting humans and domestic animals world-wide. In North America, female Culicoides sonorensis midges are important vectors of bluetongue virus (BTV) and epizootic hemorrhagic disease virus (EHDV), orbiviruses that cause significant disease in livestock and wildlife. Libraries of tissue-specific transcripts expressed in response to feeding and oral orbivirus challenge in C. sonorensis have previously been reported, but extensive genome-wide expression profiling in the midge has not. Here, we successfully used deep sequencing technologies to construct the first adult female C. sonorensis reference transcriptome, and utilized genome-wide expression profiling to elucidate the genetic response to blood and sucrose feeding over time. The adult female midge unigene consists of 19,041 genes, of which less than 7% are differentially expressed during the course of a sucrose meal, while up to 52% of the genes respond significantly in blood-fed midges, indicating hematophagy induces complex physiological processes. Many genes that were differentially expressed during blood feeding were associated with digestion (e.g. proteases, lipases), hematophagy (e.g., salivary proteins), and vitellogenesis, revealing many major metabolic and biological factors underlying these critical processes. Additionally, key genes in the vitellogenesis pathway were identified, which provides the first glimpse into the molecular basis of anautogeny for C. sonorensis. This is the first extensive transcriptome for this genus, which will serve as a framework for future expression studies, RNAi, and provide a rich dataset contributing to the ultimate goal of informing a reference genome assembly and annotation. Moreover, this study will serve as a foundation for subsequent studies of genome-wide expression analyses during early orbivirus infection and dissecting the molecular mechanisms behind vector competence in midges.