Meta-analysis of tumor- and T cell-intrinsic mechanisms of sensitization to checkpoint inhibition.

Checkpoint inhibitors (CPIs) augment adaptive immunity. Systematic pan-tumor analyses may reveal the relative importance of tumor-cell-intrinsic and microenvironmental features underpinning CPI sensitization. Here, we collated whole-exome and transcriptomic data for >1,000 CPI-treated patients across seven tumor types, utilizing standardized bioinformatics workflows and clinical outcome criteria to validate multivariable predictors of CPI sensitization. Clonal tumor mutation burden (TMB) was the strongest predictor of CPI response, followed by total TMB and CXCL9 expression. Subclonal TMB, somatic copy alteration burden, and histocompatibility leukocyte antigen (HLA) evolutionary divergence failed to attain pan-cancer significance. Dinucleotide variants were identified as a source of immunogenic epitopes associated with radical amino acid substitutions and enhanced peptide hydrophobicity/immunogenicity. Copy-number analysis revealed two additional determinants of CPI outcome supported by prior functional evidence: 9q34 (TRAF2) loss associated with response and CCND1 amplification associated with resistance. Finally, single-cell RNA sequencing (RNA-seq) of clonal neoantigen-reactive CD8 tumor-infiltrating lymphocytes (TILs), combined with bulk RNA-seq analysis of CPI-responding tumors, identified CCR5 and CXCL13 as T-cell-intrinsic markers of CPI sensitivity.

Multivariate genetic architecture reveals testosterone-driven sexual antagonism in contemporary humans.

Sex difference (SD) is ubiquitous in humans despite shared genetic architecture (SGA) between the sexes. A univariate approach, i.e., studying SD in single traits by estimating genetic correlation, does not provide a complete biological overview, because traits are not independent and are genetically correlated. The multivariate genetic architecture between the sexes can be summarized by estimating the additive genetic (co)variance across shared traits, which, apart from the cross-trait and cross-sex covariances, also includes the cross-sex-cross-trait covariances, e.g., between height in males and weight in females. Using such a multivariate approach, we investigated SD in the genetic architecture of 12 anthropometric, fat depositional, and sex-hormonal phenotypes. We uncovered sexual antagonism (SA) in the cross-sex-cross-trait covariances in humans, most prominently between testosterone and the anthropometric traits - a trend similar to phenotypic correlations. 27% of such cross-sex-cross-trait covariances were of opposite sign, contributing to asymmetry in the SGA. Intriguingly, using multivariate evolutionary simulations, we observed that the SGA acts as a genetic constraint to the evolution of SD in humans only when selection is sexually antagonistic and not concordant. Remarkably, we found that the lifetime reproductive success in both the sexes shows a positive genetic correlation with anthropometric traits, but not with testosterone. Moreover, we demonstrated that genetic variance is depleted along multivariate trait combinations in both the sexes but in different directions, suggesting absolute genetic constraint to evolution. Our results indicate that testosterone drives SA in contemporary humans and emphasize the necessity and significance of using a multivariate framework in studying SD.

dCCA: detecting differential covariation patterns between two types of high-throughput omics data.

MOTIVATION: The advent of multimodal omics data has provided an unprecedented opportunity to systematically investigate underlying biological mechanisms from distinct yet complementary angles. However, the joint analysis of multi-omics data remains challenging because it requires modeling interactions between multiple sets of high-throughput variables. Furthermore, these interaction patterns may vary across different clinical groups, reflecting disease-related biological processes. RESULTS: We propose a novel approach called Differential Canonical Correlation Analysis (dCCA) to capture differential covariation patterns between two multivariate vectors across clinical groups. Unlike classical Canonical Correlation Analysis, which maximizes the correlation between two multivariate vectors, dCCA aims to maximally recover differentially expressed multivariate-to-multivariate covariation patterns between groups. We have developed computational algorithms and a toolkit to sparsely select paired subsets of variables from two sets of multivariate variables while maximizing the differential covariation. Extensive simulation analyses demonstrate the superior performance of dCCA in selecting variables of interest and recovering differential correlations. We applied dCCA to the Pan-Kidney cohort from the Cancer Genome Atlas Program database and identified differentially expressed covariations between noncoding RNAs and gene expressions. AVAILABILITY AND IMPLEMENTATION: The R package that implements dCCA is available at https://github.com/hwiyoungstat/dCCA.

MANOCCA: a robust and computationally efficient test of covariance in high-dimension multivariate omics data.

Multivariate analysis is becoming central in studies investigating high-throughput molecular data, yet, some important features of these data are seldom explored. Here, we present MANOCCA (Multivariate Analysis of Conditional CovAriance), a powerful method to test for the effect of a predictor on the covariance matrix of a multivariate outcome. The proposed test is by construction orthogonal to tests based on the mean and variance and is able to capture effects that are missed by both approaches. We first compare the performances of MANOCCA with existing correlation-based methods and show that MANOCCA is the only test correctly calibrated in simulation mimicking omics data. We then investigate the impact of reducing the dimensionality of the data using principal component analysis when the sample size is smaller than the number of pairwise covariance terms analysed. We show that, in many realistic scenarios, the maximum power can be achieved with a limited number of components. Finally, we apply MANOCCA to 1000 healthy individuals from the Milieu Interieur cohort, to assess the effect of health, lifestyle and genetic factors on the covariance of two sets of phenotypes, blood biomarkers and flow cytometry-based immune phenotypes. Our analyses identify significant associations between multiple factors and the covariance of both omics data.

Early Transcriptional Divergence Marks Virus-Specific Primary Human CD8+ T Cells in Chronic versus Acute Infection.

Distinct molecular pathways govern the differentiation of CD8+ effector T cells into memory or exhausted T cells during acute and chronic viral infection, but these are not well studied in humans. Here, we employed an integrative systems immunology approach to identify transcriptional commonalities and differences between virus-specific CD8+ T cells from patients with persistent and spontaneously resolving hepatitis C virus (HCV) infection during the acute phase. We observed dysregulation of metabolic processes during early persistent infection that was linked to changes in expression of genes related to nucleosomal regulation of transcription, T cell differentiation, and the inflammatory response and correlated with subject age, sex, and the presence of HCV-specific CD4+ T cell populations. These early changes in HCV-specific CD8+ T cell transcription preceded the overt establishment of T cell exhaustion, making this signature a prime target in the search for the regulatory origins of T cell dysfunction in chronic viral infection.

An extensive microRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma.

By analyzing gene expression data in glioblastoma in combination with matched microRNA profiles, we have uncovered a posttranscriptional regulation layer of surprising magnitude, comprising more than 248,000 microRNA (miR)-mediated interactions. These include ∼7,000 genes whose transcripts act as miR "sponges" and 148 genes that act through alternative, nonsponge interactions. Biochemical analyses in cell lines confirmed that this network regulates established drivers of tumor initiation and subtype implementation, including PTEN, PDGFRA, RB1, VEGFA, STAT3, and RUNX1, suggesting that these interactions mediate crosstalk between canonical oncogenic pathways. siRNA silencing of 13 miR-mediated PTEN regulators, whose locus deletions are predictive of PTEN expression variability, was sufficient to downregulate PTEN in a 3'UTR-dependent manner and to increase tumor cell growth rates. Thus, miR-mediated interactions provide a mechanistic, experimentally validated rationale for the loss of PTEN expression in a large number of glioma samples with an intact PTEN locus.

Evolution of acoustic signals associated with cooperative parental behavior in a poison frog.

The emergence of complex social interactions is predicted to be an important selective force in the diversification of communication systems. Parental care presents a key social context in which to study the evolution of novel signals, as care often requires communication and behavioral coordination between parents and is an evolutionary stepping-stone toward increasingly complex social systems. Anuran amphibians (frogs and toads) are a classic model of acoustic communication and the vocal repertoires of many species have been characterized in the contexts of advertisement, courtship, and aggression, yet quantitative descriptions of calls elicited in the context of parental care are lacking. The biparental poison frog, Ranitomeya imitator, exhibits a remarkable parenting behavior in which females, cued by the calls of their male partners, feed tadpoles unfertilized eggs. Here, we characterized and compared calls across three social contexts, for the first time including a parental care context. We found that egg-feeding calls share some properties with both advertisement and courtship calls but also had unique properties. Multivariate analysis revealed high classification success for advertisement and courtship calls but misclassified nearly half of egg feeding calls as either advertisement or courtship calls. Egg feeding and courtship calls both contained less identity information than advertisement calls, as expected for signals used in close-range communication where uncertainty about identity is low and additional signal modalities may be used. Taken together, egg-feeding calls likely borrowed and recombined elements of both ancestral call types to solicit a novel, context-dependent parenting response.

Comprehensive assessment of drought resistance and recovery in kiwifruit genotypes using multivariate analysis.

As global climate change persists, ongoing warming exposes plants, including kiwifruit, to repeated cycles of drought stress and rewatering, necessitating the identification of drought-resistant genotypes for breeding purposes. To better understand the physiological mechanisms underlying drought resistance and recovery in kiwifruit, moderate (40-45% field capacity) and severe (25-30% field capacity) drought stresses were applied, followed by rewatering (80-85% field capacity) to eight kiwifruit rootstocks in this study. We then conducted a multivariate analysis of 20 indices for the assessment of drought resistance and recovery capabilities. Additionally, we identified four principal components, each playing a vital role in coping with diverse water conditions. Three optimal indicator groups were pinpointed, enhancing precision in kiwifruit drought resistance and recovery assessment and simplifying the evaluation system. Finally, MX-1 and HW were identified as representative rootstocks for future research on kiwifruit's responses to moderate and severe drought stresses. This study not only enhances our understanding of the response mechanisms of kiwifruit rootstocks to progressive drought stress and recovery but also provides theoretical guidance for reliable screening of drought-adaptive kiwifruit genotypes.

Covariate-guided Bayesian mixture of spline experts for the analysis of multivariate high-density longitudinal data.

With rapid development of techniques to measure brain activity and structure, statistical methods for analyzing modern brain-imaging data play an important role in the advancement of science. Imaging data that measure brain function are usually multivariate high-density longitudinal data and are heterogeneous across both imaging sources and subjects, which lead to various statistical and computational challenges. In this article, we propose a group-based method to cluster a collection of multivariate high-density longitudinal data via a Bayesian mixture of smoothing splines. Our method assumes each multivariate high-density longitudinal trajectory is a mixture of multiple components with different mixing weights. Time-independent covariates are assumed to be associated with the mixture components and are incorporated via logistic weights of a mixture-of-experts model. We formulate this approach under a fully Bayesian framework using Gibbs sampling where the number of components is selected based on a deviance information criterion. The proposed method is compared to existing methods via simulation studies and is applied to a study on functional near-infrared spectroscopy, which aims to understand infant emotional reactivity and recovery from stress. The results reveal distinct patterns of brain activity, as well as associations between these patterns and selected covariates.

Decoding context memories for threat in large-scale neural networks.

Humans are often tasked with determining the degree to which a given situation poses threat. Salient cues present during prior events help bring online memories for context, which plays an informative role in this process. However, it is relatively unknown whether and how individuals use features of the environment to retrieve context memories for threat, enabling accurate inferences about the current level of danger/threat (i.e. retrieve appropriate memory) when there is a degree of ambiguity surrounding the present context. We leveraged computational neuroscience approaches (i.e. independent component analysis and multivariate pattern analyses) to decode large-scale neural network activity patterns engaged during learning and inferring threat context during a novel functional magnetic resonance imaging task. Here, we report that individuals accurately infer threat contexts under ambiguous conditions through neural reinstatement of large-scale network activity patterns (specifically striatum, salience, and frontoparietal networks) that track the signal value of environmental cues, which, in turn, allows reinstatement of a mental representation, primarily within a ventral visual network, of the previously learned threat context. These results provide novel insight into distinct, but overlapping, neural mechanisms by which individuals may utilize prior learning to effectively make decisions about ambiguous threat-related contexts as they navigate the environment.

A Multivariate analysis on evoked components of Chinese semantic congruity: an OP-MEG study with EEG.

The application of wearable magnetoencephalography using optically-pumped magnetometers has drawn extensive attention in the field of neuroscience. Electroencephalogram system can cover the whole head and reflect the overall activity of a large number of neurons. The efficacy of optically-pumped magnetometer in detecting event-related components can be validated through electroencephalogram results. Multivariate pattern analysis is capable of tracking the evolution of neurocognitive processes over time. In this paper, we adopted a classical Chinese semantic congruity paradigm and separately collected electroencephalogram and optically-pumped magnetometer signals. Then, we verified the consistency of optically-pumped magnetometer and electroencephalogram in detecting N400 using mutual information index. Multivariate pattern analysis revealed the difference in decoding performance of these two modalities, which can be further validated by dynamic/stable coding analysis on the temporal generalization matrix. The results from searchlight analysis provided a neural basis for this dissimilarity at the magnetoencephalography source level and the electroencephalogram sensor level. This study opens a new avenue for investigating the brain's coding patterns using wearable magnetoencephalography and reveals the differences in sensitivity between the two modalities in reflecting neuron representation patterns.

Mass spectrometry imaging to explore molecular heterogeneity in cell culture.

Molecular analysis on the single-cell level represents a rapidly growing field in the life sciences. While bulk analysis from a pool of cells provides a general molecular profile, it is blind to heterogeneities between individual cells. This heterogeneity, however, is an inherent property of every cell population. Its analysis is fundamental to understanding the development, function, and role of specific cells of the same genotype that display different phenotypical properties. Single-cell mass spectrometry (MS) aims to provide broad molecular information for a significantly large number of cells to help decipher cellular heterogeneity using statistical analysis. Here, we present a sensitive approach to single-cell MS based on high-resolution MALDI-2-MS imaging in combination with MALDI-compatible staining and use of optical microscopy. Our approach allowed analyzing large amounts of unperturbed cells directly from the growth chamber. Confident coregistration of both modalities enabled a reliable compilation of single-cell mass spectra and a straightforward inclusion of optical as well as mass spectrometric features in the interpretation of data. The resulting multimodal datasets permit the use of various statistical methods like machine learning-driven classification and multivariate analysis based on molecular profile and establish a direct connection of MS data with microscopy information of individual cells. Displaying data in the form of histograms for individual signal intensities helps to investigate heterogeneous expression of specific lipids within the cell culture and to identify subpopulations intuitively. Ultimately, t-MALDI-2-MSI measurements at 2-µm pixel sizes deliver a glimpse of intracellular lipid distributions and reveal molecular profiles for subcellular domains.

Dengue virus infection modifies mosquito blood-feeding behavior to increase transmission to the host.

Mosquito blood-feeding behavior is a key determinant of the epidemiology of dengue viruses (DENV), the most-prevalent mosquito-borne viruses. However, despite its importance, how DENV infection influences mosquito blood-feeding and, consequently, transmission remains unclear. Here, we developed a high-resolution, video-based assay to observe the blood-feeding behavior of Aedes aegypti mosquitoes on mice. We then applied multivariate analysis on the high-throughput, unbiased data generated from the assay to ordinate behavioral parameters into complex behaviors. We showed that DENV infection increases mosquito attraction to the host and hinders its biting efficiency, the latter resulting in the infected mosquitoes biting more to reach similar blood repletion as uninfected mosquitoes. To examine how increased biting influences DENV transmission to the host, we established an in vivo transmission model with immuno-competent mice and demonstrated that successive short probes result in multiple transmissions. Finally, to determine how DENV-induced alterations of host-seeking and biting behaviors influence dengue epidemiology, we integrated the behavioral data within a mathematical model. We calculated that the number of infected hosts per infected mosquito, as determined by the reproduction rate, tripled when mosquito behavior was influenced by DENV infection. Taken together, this multidisciplinary study details how DENV infection modulates mosquito blood-feeding behavior to increase vector capacity, proportionally aggravating DENV epidemiology. By elucidating the contribution of mosquito behavioral alterations on DENV transmission to the host, these results will inform epidemiological modeling to tailor improved interventions against dengue.

Multivariate pattern analysis: a method and software to reveal, quantify, and visualize predictive association patterns in multicollinear data.

BACKGROUND: Strongly multicollinear covariates, such as those typically represented in metabolomics applications, represent a challenge for multivariate regression analysis. These challenges are commonly circumvented by reducing the number of covariates to a subset of linearly independent variables, but this strategy may lead to loss of resolution and thus produce models with poorer interpretative potential. The aim of this work was to implement and illustrate a method, multivariate pattern analysis (MVPA), which can handle multivariate covariates without compromising resolution or model quality. RESULTS: MVPA has been implemented in an open-source R package of the same name, mvpa. To facilitate the usage and interpretation of complex association patterns, mvpa has also been integrated into an R shiny app, mvpaShiny, which can be accessed on www.mvpashiny.org . MVPA utilizes a general projection algorithm that embraces a diversity of possible models. The method handles multicollinear and even linear dependent covariates. MVPA separates the variance in the data into orthogonal parts within the frame of a single joint model: one part describing the relations between covariates, outcome, and explanatory variables and another part describing the "net" predictive association pattern between outcome and explanatory variables. These patterns are visualized and interpreted in variance plots and plots for pattern analysis and ranking according to variable importance. Adjustment for a linear dependent covariate is performed in three steps. First, partial least squares regression with repeated Monte Carlo resampling is used to determine the number of predictive PLS components for a model relating the covariate to the outcome. Second, postprocessing of this PLS model by target projection provided a single component expressing the predictive association pattern between the outcome and the covariate. Third, the outcome and the explanatory variables were adjusted for the covariate by using the target score in the projection algorithm to obtain "net" data. We illustrate the main features of MVPA by investigating the partial mediation of a linearly dependent metabolomics descriptor on the association pattern between a measure of insulin resistance and lifestyle-related factors. CONCLUSIONS: Our method and implementation in R extend the range of possible analyses and visualizations that can be performed for complex multivariate data structures. The R packages are available on github.com/liningtonlab/mvpa and github.com/liningtonlab/mvpaShiny.

Focusing on the Abnormal Events of NPC1, NPC2, and NPC1L1 in Pan-Cancer and Further Constructing LUAD and KICH Prediction Models.

Cancer's high incidence and death rate jeopardize human health and life, and it has become a global public health issue. Some members of NPCs have been studied in a few cancers, but comprehensive and prognostic analysis is lacking in most cancers. In this study, we used the Cancer Genome Atlas (TCGA) data genomics and transcriptome technology to examine the differential expression and prognosis of NPCs in 33 cancer samples, as well as to investigate NPCs mutations and their effect on patient prognosis and to evaluate the methylation level of NPCs in cancer. The linked mechanisms and medication resistance were subsequently investigated in order to investigate prospective tumor therapy approaches. The relationships between NPCs and immune infiltration, immune cells, immunological regulatory substances, and immune pathways were also investigated. Finally, the LUAD and KICH prognostic prediction models were built using univariate and multivariate COX regression analysis. Additionally, the mRNA and protein levels of NPCs were also identified.

Multivariate analysis of a missense variant in CREBRF reveals associations with measures of adiposity in people of Polynesian ancestries.

The minor allele of rs373863828, a missense variant in CREB3 Regulatory Factor, is associated with several cardiometabolic phenotypes in Polynesian peoples. To better understand the variant, we tested the association of rs373863828 with a panel of correlated phenotypes (body mass index [BMI], weight, height, HDL cholesterol, triglycerides, and total cholesterol) using multivariate Bayesian association and network analyses in a Samoa cohort (n = 1632), Aotearoa New Zealand cohort (n = 1419), and combined cohort (n = 2976). An expanded set of phenotypes (adding estimated fat and fat-free mass, abdominal circumference, hip circumference, and abdominal-hip ratio) was tested in the Samoa cohort (n = 1496). In the Samoa cohort, we observed significant associations (log10 Bayes Factor [BF] ≥ 5.0) between rs373863828 and the overall phenotype panel (8.81), weight (8.30), and BMI (6.42). In the Aotearoa New Zealand cohort, we observed suggestive associations (1.5 < log10 BF < 5) between rs373863828 and the overall phenotype panel (4.60), weight (3.27), and BMI (1.80). In the combined cohort, we observed concordant signals with larger log10 BFs. In the Samoa-specific expanded phenotype analyses, we also observed significant associations between rs373863828 and fat mass (5.65), abdominal circumference (5.34), and hip circumference (5.09). Bayesian networks provided evidence for a direct association of rs373863828 with weight and indirect associations with height and BMI.

Multi-trait transcriptome-wide association studies with probabilistic Mendelian randomization.

A transcriptome-wide association study (TWAS) integrates data from genome-wide association studies and gene expression mapping studies for investigating the gene regulatory mechanisms underlying diseases. Existing TWAS methods are primarily univariate in nature, focusing on analyzing one outcome trait at a time. However, many complex traits are correlated with each other and share a common genetic basis. Consequently, analyzing multiple traits jointly through multivariate analysis can potentially improve the power of TWASs. Here, we develop a method, moPMR-Egger (multiple outcome probabilistic Mendelian randomization with Egger assumption), for analyzing multiple outcome traits in TWAS applications. moPMR-Egger examines one gene at a time, relies on its cis-SNPs that are in potential linkage disequilibrium with each other to serve as instrumental variables, and tests its causal effects on multiple traits jointly. A key feature of moPMR-Egger is its ability to test and control for potential horizontal pleiotropic effects from instruments, thus maximizing power while minimizing false associations for TWASs. In simulations, moPMR-Egger provides calibrated type I error control for both causal effects testing and horizontal pleiotropic effects testing and is more powerful than existing univariate TWAS approaches in detecting causal associations. We apply moPMR-Egger to analyze 11 traits from 5 trait categories in the UK Biobank. In the analysis, moPMR-Egger identified 13.15% more gene associations than univariate approaches across trait categories and revealed distinct regulatory mechanisms underlying systolic and diastolic blood pressures.

Clinical Features of Vaccine-Induced Immune Thrombocytopenia and Thrombosis.

BACKGROUND: Vaccine-induced immune thrombocytopenia and thrombosis (VITT) is a new syndrome associated with the ChAdOx1 nCoV-19 adenoviral vector vaccine against severe acute respiratory syndrome coronavirus 2. Data are lacking on the clinical features of and the prognostic criteria for this disorder. METHODS: We conducted a prospective cohort study involving patients with suspected VITT who presented to hospitals in the United Kingdom between March 22 and June 6, 2021. Data were collected with the use of an anonymized electronic form, and cases were identified as definite or probable VITT according to prespecified criteria. Baseline characteristics and clinicopathological features of the patients, risk factors, treatment, and markers of poor prognosis were determined. RESULTS: Among 294 patients who were evaluated, we identified 170 definite and 50 probable cases of VITT. All the patients had received the first dose of ChAdOx1 nCoV-19 vaccine and presented 5 to 48 days (median, 14) after vaccination. The age range was 18 to 79 years (median, 48), with no sex preponderance and no identifiable medical risk factors. Overall mortality was 22%. The odds of death increased by a factor of 2.7 (95% confidence interval [CI], 1.4 to 5.2) among patients with cerebral venous sinus thrombosis, by a factor of 1.7 (95% CI, 1.3 to 2.3) for every 50% decrease in the baseline platelet count, by a factor of 1.2 (95% CI, 1.0 to 1.3) for every increase of 10,000 fibrinogen-equivalent units in the baseline d-dimer level, and by a factor of 1.7 (95% CI, 1.1 to 2.5) for every 50% decrease in the baseline fibrinogen level. Multivariate analysis identified the baseline platelet count and the presence of intracranial hemorrhage as being independently associated with death; the observed mortality was 73% among patients with platelet counts below 30,000 per cubic millimeter and intracranial hemorrhage. CONCLUSIONS: The high mortality associated with VITT was highest among patients with a low platelet count and intracranial hemorrhage. Treatment remains uncertain, but identification of prognostic markers may help guide effective management. (Funded by the Oxford University Hospitals NHS Foundation Trust.).

Nontargeted Volatile Metabolite Screening and Microbial Contamination Detection in Fermentation Processes by Headspace GC-IMS.

Gas chromatography combined with ion mobility spectrometry (GC-IMS) is a powerful separation and detection technique for volatile organic compounds (VOC). This combination is characterized by exceptionally low detection limits in the low ppbv range, high 2-dimensional selectivity, and robust operation. These qualities make it an ideal tool for nontarget screening approaches. Fermentation broths contain a substantial number of VOC, either from the medium or produced by microbial metabolism, that are currently not regularly measured for process monitoring. In this study, Escherichia coli, Saccharomyces cerevisiae, Levilactobacillus brevis, and Pseudomonas fluorescens were exemplarily used as model organisms and cultivated, and the headspace was analyzed by GC-IMS. Additionally, mixed cultures for every combination of two of the microorganisms were also characterized. Multivariate data analysis of the GC-IMS data revealed that it is possible to differentiate between the microorganisms using PLS-DA with a prediction accuracy of 0.92. The mixed cultures could be separated from the pure cultures with accuracies between 0.87 and 1.00 depending on the organism. GC-IMS data correlate with the optical density and can be used to follow and model growth curves. The root mean squared errors ranged between 10 and 20% of the maximum value, depending on the organism. Peak identification with reference compounds did not reveal specific marker compounds, rather the pattern was found to be responsible for the model performance.

Multivariate investigation of Moringa oleifera morpho-physiological and biochemical traits under various water regimes.

BACKGROUND: The climatic changes crossing the world menace the green life through limitation of water availability. The goal of this study was to determine whether Moringa oleifera Lam. trees cultivated under Tunisian arid climate, retain their tolerance ability to tolerate accentuated environmental stress factors such as drought and salinity. For this reason, the seeds of M. oleifera tree planted in Bouhedma Park (Tunisian arid area), were collected, germinated, and grown in the research area at the National Institute of Research in Rural Engineering, Waters and Forests (INRGREF) of Tunis (Tunisia). The three years aged trees were exposed to four water-holding capacities (25, 50, 75, and 100%) for 60 days to realise this work. RESULTS: Growth change was traduced by the reduction of several biometric parameters and fluorescence (Fv/Fm) under severe water restriction (25 and 50%). Whereas roots presented miraculous development in length face to the decrease of water availability (25 and 50%) in their rhizospheres. The sensitivity to drought-induced membrane damage (Malondialdehyde (MDA) content) and reactive oxygen species (ROS) liberation (hydrogen peroxide (H2O2) content) was highly correlated with ROS antiradical scavenging (ferric reducing antioxidant power (FRAP) and (2, 2'-diphenyl-1-picrylhydrazyle (DPPH)), phenolic components and osmolytes accumulation. The drought stress tolerance of M. oleifera trees was associated with a dramatic stimulation of superoxide dismutase (SOD), catalase (CAT), glutathione reductase (GR), ascorbate peroxidase (APX), and glutathione peroxidase (GPX) activities. CONCLUSION: Based on the several strategies adopted, integrated M. oleifera can grow under drought stress as accentuated adverse environmental condition imposed by climate change.