Robust Co-clustering to Discover Toxicogenomic Biomarkers and Their Regulatory Doses of Chemical Compounds Using Logistic Probabilistic Hidden Variable Model

Multivariate Analysis of Morpho-Physiological Traits Reveals Differential Drought Tolerance Potential of Bread Wheat Genotypes at the Seedling Stage

Drought is one of the foremost environmental stresses that can severely limit crop growth and productivity by disrupting various physiological processes. In this study, the drought tolerance potential of 127 diverse bread wheat genotypes was evaluated by imposing polyethylene glycol (PEG)-induced drought followed by multivariate analysis of several growth-related attributes. Results showed significant variations in the mean values of different morpho-physiological traits due to PEG-induced drought effects. Correlation analysis revealed that most of the studied traits were significantly correlated among them. The robust hierarchical co-clustering indicated that all the genotypes were clustered into four major groups, with cluster 4 (26 genotypes) being, in general, drought-tolerant followed by cluster 1 (19 genotypes) whereas, cluster 2 (55 genotypes) and 3 (27 genotypes) being drought-sensitive. Linear discriminant analysis (LDA) confirmed that around 90% of the genotypes were correctly assigned to clusters. Squared distance (D2) analysis indicated that the clusters differed significantly from each other. Principal component analysis (PCA) and genotype by trait biplot analysis showed that the first three components accounted for 71.6% of the total variation, with principal component (PC) 1 accounting for 35.4%, PC2 for 24.6% and PC3 for 11.6% of the total variation. Both PCA and LDA revealed that dry weights, tissue water content, cell membrane stability, leaf relative water content, root-shoot weight ratio and seedling vigor index played the most important discriminatory roles in explaining drought tolerance variations among 127 wheat genotypes. Our results conclude that the drought-tolerant and -sensitive wheat genotypes identified in this study would offer valuable genetic tools for further improvement of wheat productivity in arid and semi-arid regions during this time of unpredictable climate change.

An introduction to new robust linear and monotonic correlation coefficients

BACKGROUND

The most common measure of association between two continuous variables is the Pearson correlation (Maronna et al. in Safari an OMC. Robust statistics, 2019. https://login.proxy.bib.uottawa.ca/login?url=https://learning.oreilly.com/library/view/-/9781119214687/?ar&orpq&email=^u). When outliers are present, Pearson does not accurately measure association and robust measures are needed. This article introduces three new robust measures of correlation: Taba (T), TabWil (TW), and TabWil rank (TWR). The correlation estimators T and TW measure a linear association between two continuous or ordinal variables; whereas TWR measures a monotonic association. The robustness of these proposed measures in comparison with Pearson (P), Spearman (S), Quadrant (Q), Median (M), and Minimum Covariance Determinant (MCD) are examined through simulation. Taba distance is used to analyze genes, and statistical tests were used to identify those genes most significantly associated with Williams Syndrome (WS).

RESULTS

Based on the root mean square error (RMSE) and bias, the three proposed correlation measures are highly competitive when compared to classical measures such as P and S as well as robust measures such as Q, M, and MCD. Our findings indicate TBL2 was the most significant gene among patients diagnosed with WS and had the most significant reduction in gene expression level when compared with control (P value = 6.37E-05).

CONCLUSIONS

Overall, when the distribution is bivariate Log-Normal or bivariate Weibull, TWR performs best in terms of bias and T performs best with respect to RMSE. Under the Normal distribution, MCD performs well with respect to bias and RMSE; but TW, TWR, T, S, and P correlations were in close proximity. The identification of TBL2 may serve as a diagnostic tool for WS patients. A Taba R package has been developed and is available for use to perform all necessary computations for the proposed methods.

Assessment of Drugs Toxicity and Associated Biomarker Genes Using Hierarchical Clustering

Background and objectives: Assessment of drugs toxicity and associated biomarker genes is one of the most important tasks in the pre-clinical phase of drug development pipeline as well as in toxicogenomic studies. There are few statistical methods for the assessment of doses of drugs (DDs) toxicity and their associated biomarker genes. However, these methods consume more time for computation of the model parameters using the EM (expectation-maximization) based iterative approaches. To overcome this problem, in this paper, an attempt is made to propose an alternative approach based on hierarchical clustering (HC) for the same purpose. Methods and materials: There are several types of HC approaches whose performance depends on different similarity/distance measures. Therefore, we explored suitable combinations of distance measures and HC methods based on Japanese Toxicogenomics Project (TGP) datasets for better clustering/co-clustering between DDs and genes as well as to detect toxic DDs and their associated biomarker genes. Results: We observed that Word’s HC method with each of Euclidean, Manhattan, and Minkowski distance measures produces better clustering/co-clustering results. For an example, in the case of the glutathione metabolism pathway (GMP) dataset LOC100359539/Rrm2, Gpx6, RGD1562107, Gstm4, Gstm3, G6pd, Gsta5, Gclc, Mgst2, Gsr, Gpx2, Gclm, Gstp1, LOC100912604/Srm, Gstm4, Odc1, Gsr, Gss are the biomarker genes and Acetaminophen_Middle, Acetaminophen_High, Methapyrilene_High, Nitrofurazone_High, Nitrofurazone_Middle, Isoniazid_Middle, Isoniazid_High are their regulatory (associated) DDs explored by our proposed co-clustering algorithm based on the distance and HC method combination Euclidean: Word. Similarly, for the peroxisome proliferator-activated receptor signaling pathway (PPAR-SP) dataset Cpt1a, Cyp8b1, Cyp4a3, Ehhadh, Plin5, Plin2, Fabp3, Me1, Fabp5, LOC100910385, Cpt2, Acaa1a, Cyp4a1, LOC100365047, Cpt1a, LOC100365047, Angptl4, Aqp7, Cpt1c, Cpt1b, Me1 are the biomarker genes and Aspirin_Low, Aspirin_Middle, Aspirin_High, Benzbromarone_Middle, Benzbromarone_High, Clofibrate_Middle, Clofibrate_High, WY14643_Low, WY14643_High, WY14643_Middle, Gemfibrozil_Middle, Gemfibrozil_High are their regulatory DDs. Conclusions: Overall, the methods proposed in this article, co-cluster the genes and DDs as well as detect biomarker genes and their regulatory DDs simultaneously consuming less time compared to other mentioned methods. The results produced by the proposed methods have been validated by the available literature and functional annotation.

Robust identification of significant interactions between toxicogenomic biomarkers and their regulatory chemical compounds using logistic moving range chart

Identification of significant interactions between genes and chemical compounds/drugs is an important issue in toxicogenomic studies as well as in drug discovery and development. There are some online and offline computational tools for toxicogenomic data analysis to identify the biomarker genes and their regulatory chemical compounds/drugs. However, none of the researchers has considered yet the identification of significant interactions between genes and compounds. Therefore, in this paper, we have discussed two approaches namely moving range chart (MRC) and logistic moving range chart (LMRC) for the identification of significant up-regulatory (UpR) and down-regulatory (DnR) gene-compound interactions as well as toxicogenomic biomarkers and their regulatory chemical compounds/drugs. We have investigated the performance of both MRC and LMRC approaches using simulated datasets. Simulation results show that both approaches perform almost equally in absence of outliers. However, in presence of outliers, the LMRC shows much better performance than the MRC. In case of real life toxicogenomic data analysis, the proposed LMRC approach detected some important down-regulated biomarker genes those were not detected by other approaches. Therefore, in this paper, our proposal is to use LMRC for robust identification of significant interactions between genes and chemical compounds/drugs.