Predicting yield of individual field-grown rapeseed plants from rosette-stage leaf gene expression.

In the plant sciences, results of laboratory studies often do not translate well to the field. To help close this lab-field gap, we developed a strategy for studying the wiring of plant traits directly in the field, based on molecular profiling and phenotyping of individual plants. Here, we use this single-plant omics strategy on winter-type Brassica napus (rapeseed). We investigate to what extent early and late phenotypes of field-grown rapeseed plants can be predicted from their autumnal leaf gene expression, and find that autumnal leaf gene expression not only has substantial predictive power for autumnal leaf phenotypes but also for final yield phenotypes in spring. Many of the top predictor genes are linked to developmental processes known to occur in autumn in winter-type B. napus accessions, such as the juvenile-to-adult and vegetative-to-reproductive phase transitions, indicating that the yield potential of winter-type B. napus is influenced by autumnal development. Our results show that single-plant omics can be used to identify genes and processes influencing crop yield in the field.

Impact of perfluoroalkyl substances (PFAS) and PFAS mixtures on lipid metabolism in differentiated HepaRG cells as a model for human hepatocytes.

Per- and polyfluoroalkyl substances (PFAS) are environmental contaminants with various adverse health effects in humans including disruption of lipid metabolism. Aim of the present study was to elucidate the molecular mechanisms of PFAS-mediated effects on lipid metabolism in human cells. Here, we examined the impact of a number of PFAS (PFOS, PFOA, PFNA, PFDA, PFHxA, PFBA, PFHxS, PFBS, HFPO-DA, and PMPP) and of some exposure-relevant PFAS mixtures being composed of PFOS, PFOA, PFNA and PFHxS on lipid metabolism in human HepaRG cells, an in vitro model for human hepatocytes. At near cytotoxic concentrations, the selected PFAS and PFAS mixtures induced triglyceride accumulation in HepaRG cells and consistently affected the expression of marker genes for steatosis, as well as PPARα target genes and genes related to lipid and cholesterol metabolism, pointing to common molecular mechanisms of PFAS in disrupting cellular lipid and cholesterol homeostasis. PPARα activation was examined by a transactivation assay in HEK293T cells, and synergistic effects were observed for the selected PFAS mixtures at sum concentrations higher than 25 µM, whereas additivity was observed at sum concentrations lower than 25 µM. Of note, any effect observed in the in vitro assays occurred at PFAS concentrations that were at least four to five magnitudes above real-life internal exposure levels of the general population.

A multi-omics approach to elucidate okadaic acid-induced changes in human HepaRG hepatocarcinoma cells.

Okadaic acid (OA), a prevalent marine biotoxin found in shellfish, is known for causing acute gastrointestinal symptoms. Despite its potential to reach the bloodstream and the liver, the hepatic effects of OA are not well understood, highlighting a significant research gap. This study aims to comprehensively elucidate the impact of OA on the liver by examining the transcriptome, proteome, and phosphoproteome alterations in human HepaRG liver cells exposed to non-cytotoxic OA concentrations. We employed an integrative multi-omics approach, encompassing RNA sequencing, shotgun proteomics, phosphoproteomics, and targeted DigiWest analysis. This enabled a detailed exploration of gene and protein expression changes, alongside phosphorylation patterns under OA treatment. The study reveals concentration- and time-dependent deregulation in gene and protein expression, with a significant down-regulation of xenobiotic and lipid metabolism pathways. Up-regulated pathways include actin crosslink formation and a deregulation of apoptotic pathways. Notably, our results revealed that OA, as a potent phosphatase inhibitor, induces alterations in actin filament organization. Phosphoproteomics data highlighted the importance of phosphorylation in enzyme activity regulation, particularly affecting proteins involved in the regulation of the cytoskeleton. OA's inhibition of PP2A further leads to various downstream effects, including alterations in protein translation and energy metabolism. This research expands the understanding of OA's systemic impact, emphasizing its role in modulating the phosphorylation landscape, which influences crucial cellular processes. The results underscore OA's multifaceted effects on the liver, particularly through PP2A inhibition, impacting xenobiotic metabolism, cytoskeletal dynamics, and energy homeostasis. These insights enhance our comprehension of OA's biological significance and potential health risks.

Characterization of hepatic zonation in mice by mass-spectrometric and antibody-based proteomics approaches.

Periportal and perivenous hepatocytes show zonal heterogeneity in metabolism and signaling. Here, hepatic zonation in mouse liver was analyzed by non-targeted mass spectrometry (MS) and by the antibody-based DigiWest technique, yielding a comprehensive overview of protein expression in periportal and perivenous hepatocytes. Targeted immunoaffinity-based proteomics were used to substantiate findings related to drug metabolism. 165 (MS) and 82 (DigiWest) zonated proteins were identified based on the selected criteria for statistical significance, including 7 (MS) and 43 (DigiWest) proteins not identified as zonated before. New zonated proteins especially comprised kinases and phosphatases related to growth factor-dependent signaling, with mainly periportal localization. Moreover, the mainly perivenous zonation of a large panel of cytochrome P450 enzymes was characterized. DigiWest data were shown to complement the MS results, substantially improving possibilities to bioinformatically identify zonated biological processes. Data mining revealed key regulators and pathways preferentially active in either periportal or perivenous hepatocytes, with ß-catenin signaling and nuclear xeno-sensing receptors as the most prominent perivenous regulators, and several kinase- and G-protein-dependent signaling cascades active mainly in periportal hepatocytes. In summary, the present data substantially broaden our knowledge of hepatic zonation in mouse liver at the protein level.

Proteomic analysis of hepatic effects of phenobarbital in mice with humanized liver.

Activation of the constitutive androstane receptor (CAR) may induce adaptive but also adverse effects in rodent liver, including the induction of drug-metabolizing enzymes, transient hepatocellular proliferation, and promotion of liver tumor growth. Human relevance of CAR-related adverse hepatic effects is controversially debated. Here, we used the chimeric FRG-KO mouse model with livers largely repopulated by human hepatocytes, in order to study human hepatocytes and their response to treatment with the model CAR activator phenobarbital (PB) in vivo. Mice received an intraperitoneal injection with 50 mg/kg body weight PB or saline, and were sacrificed after 72-144 h. Non-repopulated FRG-KO mice were used as additional control. Comprehensive proteomics datasets were generated by merging data obtained by targeted as well as non-targeted proteomics approaches. For the first time, a novel proteomics workflow was established to comparatively analyze the effects of PB on human and murine proteins within one sample. Analysis of merged proteome data sets and bioinformatics data mining revealed comparable responses in murine and human hepatocytes with respect to nuclear receptor activation and induction of xenobiotic metabolism. By contrast, activation of MYC, a key regulator of proliferation, was predicted only for mouse but not human hepatocytes. Analyses of 5-bromo-2'-deoxyuridine incorporation confirmed this finding. In summary, this study for the first time presents a comprehensive proteomic analysis of CAR-dependent effects in human and mouse hepatocytes from humanized FRG-KO mice. The data support the hypothesis that PB does induce adaptive metabolic responses, but not hepatocellular proliferation in human hepatocytes in vivo.

Comparative Analysis of Transcriptional Responses to Genotoxic and Non-Genotoxic Agents in the Blood Cell Model TK6 and the Liver Model HepaRG.

Transcript signatures are a promising approach to identify and classify genotoxic and non-genotoxic compounds and are of interest as biomarkers or for future regulatory application. Not much data, however, is yet available about the concordance of transcriptional responses in different cell types or tissues. Here, we analyzed transcriptomic responses to selected genotoxic food contaminants in the human p53-competent lymphoblastoid cell line TK6 using RNA sequencing. Responses to treatment with five genotoxins, as well as with four non-genotoxic liver toxicants, were compared with previously published gene expression data from the human liver cell model HepaRG. A significant overlap of the transcriptomic changes upon genotoxic stress was detectable in TK6 cells, whereas the comparison with the HepaRG model revealed considerable differences, which was confirmed by bioinformatic data mining for cellular upstream regulators or pathways. Taken together, the study presents a transcriptomic signature for genotoxin exposure in the human TK6 blood cell model. The data demonstrate that responses in different cell models have considerable variations. Detection of a transcriptomic genotoxin signature in blood cells indicates that gene expression analyses of blood samples might be a valuable approach to also estimate responses to toxic exposure in target organs such as the liver.

Using single-plant-omics in the field to link maize genes to functions and phenotypes.

Most of our current knowledge on plant molecular biology is based on experiments in controlled laboratory environments. However, translating this knowledge from the laboratory to the field is often not straightforward, in part because field growth conditions are very different from laboratory conditions. Here, we test a new experimental design to unravel the molecular wiring of plants and study gene-phenotype relationships directly in the field. We molecularly profiled a set of individual maize plants of the same inbred background grown in the same field and used the resulting data to predict the phenotypes of individual plants and the function of maize genes. We show that the field transcriptomes of individual plants contain as much information on maize gene function as traditional laboratory-generated transcriptomes of pooled plant samples subject to controlled perturbations. Moreover, we show that field-generated transcriptome and metabolome data can be used to quantitatively predict individual plant phenotypes. Our results show that profiling individual plants in the field is a promising experimental design that could help narrow the lab-field gap.

Transcriptomics analysis of hepatotoxicity induced by the pesticides imazalil, thiacloprid and clothianidin alone or in binary mixtures in a 28-day study in female Wistar rats.

Co-occurrence of pesticide residues in food commodities raises a potential safety issue as their mixture effects on human health are largely unknown. In a previous study, we reported the toxicological effects (pathology and histopathology) of imazalil (IMZ), thiacloprid (THI), and clothianidin (CTD) alone and in binary mixtures in a 28-day oral gavage study in female Wistar rats. Five dose levels (up to 350 mg/kg body weight/day) ranging from a typical toxicological reference value to a clear effect dose were applied. In the present study, we undertook a transcriptomics analysis of rat livers by means of total RNA sequencing (RNA-Seq). Bioinformatic data analysis involving Ingenuity Pathway Analysis (IPA) was used to gain mechanistic information on hepatotoxicity-related pathways affected after treatment with the pesticides, alone and in mixtures. Our data show that 2986 genes were differentially regulated by CTD while IMZ and THI had effects on 194 and 225 genes, respectively. All three individual compounds shared a common subset of genes whose network is associated with xenobiotic metabolism and nuclear receptor activation. Similar networks were retrieved for the mixtures. Alterations in the expression of individual genes were in line with the assumption of dose addition. Our results bring new insight into the hepatotoxicity mechanisms of IMZ, THI, and CTD and their mixtures.

Pyrrolizidine alkaloid-induced transcriptomic changes in rat lungs in a 28-day subacute feeding study.

Pyrrolizidine alkaloids (PAs) are secondary plant metabolites synthesized by a wide range of plants as protection against herbivores. These toxins are found worldwide and pose a threat to human health. PAs induce acute effects like hepatic sinusoidal obstruction syndrome and pulmonary arterial hypertension. Moreover, chronic exposure to low doses can induce cancer and liver cirrhosis in laboratory animals. The mechanisms causing hepatotoxicity have been investigated previously. However, toxic effects in the lung are less well understood, and especially data on the correlation effects with individual chemical structures of different PAs are lacking. The present study focuses on the identification of gene expression changes in vivo in rat lungs after exposure to six structurally different PAs (echimidine, heliotrine, lasiocarpine, senecionine, senkirkine, and platyphylline). Rats were treated by gavage with daily doses of 3.3 mg PA/kg bodyweight for 28 days and transcriptional changes in the lung and kidney were investigated by whole-genome microarray analysis. The results were compared with recently published data on gene regulation in the liver. Using bioinformatics data mining, we identified inflammatory responses as a predominant feature in rat lungs. By comparison, in liver, early molecular consequences to PAs were characterized by alterations in cell-cycle regulation and DNA damage response. Our results provide, for the first time, information about early molecular effects in lung tissue after subacute exposure to PAs, and demonstrates tissue-specificity of PA-induced molecular effects.

Genome-Wide Approach to Identify Quantitative Trait Loci for Drought Tolerance in Tetraploid Potato (Solanum tuberosum L.).

Drought represents a major abiotic stress factor negatively affecting growth, yield and tuber quality of potatoes. Quantitative trait locus (QTL) analyses were performed in cultivated potatoes for drought tolerance index DRYM (deviation of relative starch yield from the experimental median), tuber starch content, tuber starch yield, tuber fresh weight, selected transcripts and metabolites under control and drought stress conditions. Eight genomic regions of major interest for drought tolerance were identified, three representing standalone DRYM QTL. Candidate genes, e.g., from signaling pathways for ethylene, abscisic acid and brassinosteroids, and genes encoding cell wall remodeling enzymes were identified within DRYM QTL. Co-localizations of DRYM QTL and QTL for tuber starch content, tuber starch yield and tuber fresh weight with underlying genes of the carbohydrate metabolism were observed. Overlaps of DRYM QTL with metabolite QTL for ribitol or galactinol may indicate trade-offs between starch and compatible solute biosynthesis. Expression QTL confirmed the drought stress relevance of selected transcripts by overlaps with DRYM QTL. Bulked segregant analyses combined with next-generation sequencing (BSAseq) were used to identify mutations in genes under the DRYM QTL on linkage group 3. Future analyses of identified genes for drought tolerance will give a better insight into drought tolerance in potatoes.

Polycyclic Aromatic Hydrocarbons Activate the Aryl Hydrocarbon Receptor and the Constitutive Androstane Receptor to Regulate Xenobiotic Metabolism in Human Liver Cells.

Polycyclic aromatic hydrocarbons (PAHs) are environmental pollutants produced by incomplete combustion of organic matter. They induce their own metabolism by upregulating xenobiotic-metabolizing enzymes such as cytochrome P450 monooxygenase 1A1 (CYP1A1) by activating the aryl hydrocarbon receptor (AHR). However, previous studies showed that individual PAHs may also interact with the constitutive androstane receptor (CAR). Here, we studied ten PAHs, different in carcinogenicity classification, for their potential to activate AHR- and CAR-dependent luciferase reporter genes in human liver cells. The majority of investigated PAHs activated AHR, while non-carcinogenic PAHs tended to activate CAR. We further characterized gene expression, protein abundancies and activities of the AHR targets CYP1A1 and 1A2, and the CAR target CYP2B6 in human HepaRG hepatoma cells. Enzyme induction patterns strongly resembled the profiles obtained at the receptor level, with AHR-activating PAHs inducing CYP1A1/1A2 and CAR-activating PAHs inducing CYP2B6. In summary, this study provides evidence that beside well-known activation of AHR, some PAHs also activate CAR, followed by subsequent expression of respective target genes. Furthermore, we found that an increased PAH ring number is associated with AHR activation as well as the induction of DNA double-strand breaks, whereas smaller PAHs activated CAR but showed no DNA-damaging potential.

Molecular signatures associated with increased freezing tolerance due to low temperature memory in Arabidopsis.

Alternating temperatures require fast and coordinated adaptation responses of plants. Cold acclimation has been extensively investigated and results in increased freezing tolerance in Arabidopsis thaliana. Here, we show that the two Arabidopsis accessions, Col-0 and N14, which differ in their freezing tolerance, showed memory of cold acclimation, that is, cold priming. Freezing tolerance was higher in plants exposed to cold priming at 4°C, a lag phase at 20°C, and a second triggering cold stress (4°C) than in plants that were only cold primed. To our knowledge, this is the first report on cold memory improving plant freezing tolerance. The triggering response was distinguishable from the priming response at the levels of gene expression (RNA-Seq), lipid (ultraperformance liquid chromatography-mass spectrometry), and metabolite composition (gas chromatography-mass spectrometry). Transcriptomic responses pointed to induced lipid, secondary metabolism, and stress in Col-0 and growth-related functions in N14. Specific accumulation of lipids included arabidopsides with possible functions as signalling molecules or precursors of jasmonic acid. Whereas cold-induced metabolites such as raffinose and its precursors were maintained in N14 during the lag phase, they were strongly accumulated in Col-0 after the cold trigger. This indicates genetic differences in the transcriptomic and metabolic patterns during cold memory.

Metabolite and transcript markers for the prediction of potato drought tolerance.

Potato (Solanum tuberosum L.) is one of the most important food crops worldwide. Current potato varieties are highly susceptible to drought stress. In view of global climate change, selection of cultivars with improved drought tolerance and high yield potential is of paramount importance. Drought tolerance breeding of potato is currently based on direct selection according to yield and phenotypic traits and requires multiple trials under drought conditions. Marker-assisted selection (MAS) is cheaper, faster and reduces classification errors caused by noncontrolled environmental effects. We analysed 31 potato cultivars grown under optimal and reduced water supply in six independent field trials. Drought tolerance was determined as tuber starch yield. Leaf samples from young plants were screened for preselected transcript and nontargeted metabolite abundance using qRT-PCR and GC-MS profiling, respectively. Transcript marker candidates were selected from a published RNA-Seq data set. A Random Forest machine learning approach extracted metabolite and transcript markers for drought tolerance prediction with low error rates of 6% and 9%, respectively. Moreover, by combining transcript and metabolite markers, the prediction error was reduced to 4.3%. Feature selection from Random Forest models allowed model minimization, yielding a minimal combination of only 20 metabolite and transcript markers that were successfully tested for their reproducibility in 16 independent agronomic field trials. We demonstrate that a minimum combination of transcript and metabolite markers sampled at early cultivation stages predicts potato yield stability under drought largely independent of seasonal and regional agronomic conditions.

The drought response of potato reference cultivars with contrasting tolerance.

Systems responses to drought stress of four potato reference cultivars with differential drought tolerance (Solanum tuberosum L.) were investigated by metabolome profiling and RNA sequencing. Systems analysis was based on independent field and greenhouse trials. Robust differential drought responses across all cultivars under both conditions comprised changes of proline, raffinose, galactinol, arabitol, arabinonic acid, chlorogenic acid and 102 transcript levels. The encoded genes contained a high proportion of heat shock proteins and proteins with signalling or regulatory functions, for example, a homolog of abscisic acid receptor PYL4. Constitutive differences of the tolerant compared with the sensitive cultivars included arbutin, octopamine, ribitol and 248 transcripts. The gene products of many of these transcripts were pathogen response related, such as receptor kinases, or regulatory proteins, for example, a homolog of the Arabidopsis FOUR LIPS MYB-regulator of stomatal cell proliferation. Functional enrichment analyses imply heat stress as a major acclimation component of potato leaves to long-term drought stress. Enhanced heat stress during drought can be caused by loss of transpiration cooling. This effect and CO2 limitation are the main consequences of drought-induced or abscisic acid-induced stomatal closure. Constitutive differences in metabolite and transcript levels between tolerant and sensitive cultivars indicate interactions of drought tolerance and pathogen resistance in potato.

Whole-genome sequencing of tetraploid potato varieties reveals different strategies for drought tolerance.

Climate changes leading to increasingly longer seasonal drought periods in large parts of the world increase the necessity for breeding drought-tolerant crops. Cultivated potato (Solanum tuberosum), the third most important vegetable crop worldwide, is regarded as drought-sensitive due to its shallow root architecture. Two German tetraploid potato cultivars differing in drought tolerance and their F1-progeny were evaluated under various drought scenarios. Bulked segregant analyses were combined with whole-genome sequencing (BSA-Seq) using contrasting bulks of drought-tolerant and drought-sensitive F1-clones. Applying QTLseqr, 15 QTLs comprising 588,983 single nucleotide polymorphisms (SNPs) in 2325 genes associated with drought stress tolerance were identified. SeqSNP analyses in an association panel of 34 mostly starch potato varieties using 1-8 SNPs for each of 188 selected genes narrowed the number of candidate genes down to 10. In addition, ent-kaurene synthase B was the only gene present under QTL 10. Eight of the identified genes (StABP1, StBRI1, StKS, StLEA, StPKSP1, StPKSP2, StYAB5, and StZOG1) address plant development, the other three genes (StFATA, StHGD and StSYP) contribute to plant protection under drought stress. Allelic variation in these genes might be explored in future breeding for drought-tolerant potato varieties.

Corrigendum: Proximal hyperspectral imaging detects diurnal and drought-induced changes in maize physiology.

[This corrects the article DOI: 10.3389/fpls.2021.640914.].

Filling Analytical Gaps in Allergen Detection─Real-Time PCR for the Detection of Commercially Relevant Cephalopods and Gastropods in Food.

Mollusks belong to the group of shellfish, which are considered to be among the elicitors of severe food allergies worldwide. In recent years, numerous PCR detection methods have been developed for other shellfish such as crustaceans. However, cephalopods and gastropods were not considered in the development of these shellfish detection systems. In this study, we have developed highly specific real-time PCR methods for the comprehensive detection of all commercially relevant cephalopod species and the gastropod families Helicidae, Buccinidae, and Muricidae in food matrices. In total, we cross-tested over 100 animal and plant species to show the specificity of our systems. The limit of detection (LOD12) was set at 1 pg of cephalopod and gastropod DNA or 10 ppm (mg/kg) spiked in a vegetarian food product. The robustness of the protocol was confirmed by testing multiple parameters while cooking and autoclaving of samples ensured the practical applicability of the systems.

Monitoring of drought stress and transpiration rate using proximal thermal and hyperspectral imaging in an indoor automated plant phenotyping platform.

BACKGROUND: Thermography is a popular tool to assess plant water-use behavior, as plant temperature is influenced by transpiration rate, and is commonly used in field experiments to detect plant water deficit. Its application in indoor automated phenotyping platforms is still limited and mainly focuses on differences in plant temperature between genotypes or treatments, instead of estimating stomatal conductance or transpiration rate. In this study, the transferability of commonly used thermography analysis protocols from the field to greenhouse phenotyping platforms was evaluated. In addition, the added value of combining thermal infrared (TIR) with hyperspectral imaging to monitor drought effects on plant transpiration rate (E) was evaluated. RESULTS: The sensitivity of commonly used TIR indices to detect drought-induced and genotypic differences in water status was investigated in eight maize inbred lines in the automated phenotyping platform PHENOVISION. Indices that normalized plant temperature for vapor pressure deficit and/or air temperature at the time of imaging were most sensitive to drought and could detect genotypic differences in the plants' water-use behavior. However, these indices were not strongly correlated to stomatal conductance and E. The canopy temperature depression index, the crop water stress index and the simplified stomatal conductance index were more suitable to monitor these traits, and were consequently used to develop empirical E prediction models by combining them with hyperspectral indices and/or environmental variables. Different modeling strategies were evaluated, including single index-based, machine learning and mechanistic models. Model comparison showed that combining multiple TIR indices in a random forest model can improve E prediction accuracy, and that the contribution of the hyperspectral data is limited when multiple indices are used. However, the empirical models trained on one genotype were not transferable to all eight inbred lines. CONCLUSION: Overall, this study demonstrates that existing TIR indices can be used to monitor drought stress and develop E prediction models in an indoor setup, as long as the indices normalize plant temperature for ambient air temperature or relative humidity.

Identification of microRNAs Implicated in Modulating Senecionine-Induced Liver Toxicity in HepaRG Cells.

1,2-unsaturated Pyrrolizidine Alkaloids (PAs) are secondary plant metabolites that occur as food contaminants. Upon consumption, they can cause severe liver damage. PAs have been shown to induce apoptosis, to have cytotoxic and genotoxic effects, and to impair bile acid homeostasis in the human hepatoma cell line HepaRG. The major mode of action of PAs is DNA- and protein-adduct formation. Beyond that, nuclear receptor activation has only been observed for one receptor and two PAs, yielding the possibility that other cellular mediators are involved in PA-mediated toxicity. Here, the mode of action of Senecionine (Sc), a prominent and ubiquitous representative of hepatotoxic PAs, was investigated by analyzing 7 hepatic microRNAs (miRNAs) in HepaRG cells. Ultimately, 11 target genes that were predicted with Ingenuity Pathway Analysis software (IPA) were found to be significantly downregulated, while their assigned miRNAs showed significant upregulation of gene expression. According to IPA, these targets are positively correlated with apoptosis and cellular death and are involved in diseases such as hepatocellular carcinoma. Subsequent antagomiR-inhibition analysis revealed a significant correlation between PA-induced miRNA-4434 induction and P21-Activated Kinase-1 (PAK1) downregulation. PAK1 downregulation is usually associated with cell cycle arrest, suggesting a new function of Sc-mediated toxicity in human liver cells.