Bias and mean squared error in Mendelian randomization with invalid instrumental variables.

Mendelian randomization (MR) is a statistical method that utilizes genetic variants as instrumental variables (IVs) to investigate causal relationships between risk factors and outcomes. Although MR has gained popularity in recent years due to its ability to analyze summary statistics from genome-wide association studies (GWAS), it requires a substantial number of single nucleotide polymorphisms (SNPs) as IVs to ensure sufficient power for detecting causal effects. Unfortunately, the complex genetic heritability of many traits can lead to the use of invalid IVs that affect both the risk factor and the outcome directly or through an unobserved confounder. This can result in biased and imprecise estimates, as reflected by a larger mean squared error (MSE). In this study, we focus on the widely used two-stage least squares (2SLS) method and derive formulas for its bias and MSE when estimating causal effects using invalid IVs. Using those formulas, we identify conditions under which the 2SLS estimate is unbiased and reveal how the independent or correlated pleiotropic effects influence the accuracy and precision of the 2SLS estimate. We validate these formulas through extensive simulation studies and demonstrate the application of those formulas in an MR study to evaluate the causal effect of the waist-to-hip ratio on various sleeping patterns. Our results can aid in designing future MR studies and serve as benchmarks for assessing more sophisticated MR methods.

Cohort-based smoothing methods for age-specific contact rates.

The use of social contact rates is widespread in infectious disease modeling since it has been shown that they are key driving forces of important epidemiological parameters. Quantification of contact patterns is crucial to parameterize dynamic transmission models and to provide insights on the (basic) reproduction number. Information on social interactions can be obtained from population-based contact surveys, such as the European Commission project POLYMOD. Estimation of age-specific contact rates from these studies is often done using a piecewise constant approach or bivariate smoothing techniques. For the latter, typically, smoothness is introduced in the dimensions of the respondent's and contact's age (i.e., the rows and columns of the social contact matrix). We propose a smoothing constrained approach-taking into account the reciprocal nature of contacts-introducing smoothness over the diagonal (including all subdiagonals) of the social contact matrix. This modeling approach is justified assuming that when people age their contact behavior changes smoothly. We call this smoothing from a cohort perspective. Two approaches that allow for smoothing over social contact matrix diagonals are proposed, namely (i) reordering of the diagonal components of the contact matrix and (ii) reordering of the penalty matrix ensuring smoothness over the contact matrix diagonals. Parameter estimation is done in the likelihood framework by using constrained penalized iterative reweighted least squares. A simulation study underlines the benefits of cohort-based smoothing. Finally, the proposed methods are illustrated on the Belgian POLYMOD data of 2006. Code to reproduce the results of the article can be downloaded on this GitHub repository https://github.com/oswaldogressani/Cohort_smoothing.

Exploring Phylogenetic Signal in Multivariate Phenotypes by Maximizing Blomberg's K.

Due to the hierarchical structure of the tree of life, closely related species often resemble each other more than distantly related species; a pattern termed phylogenetic signal. Numerous univariate statistics have been proposed as measures of phylogenetic signal for single phenotypic traits, but the study of phylogenetic signal for multivariate data, as is common in modern biology, remains challenging. Here we introduce a new method to explore phylogenetic signal in multivariate phenotypes. Our approach decomposes the data into linear combinations with maximal (or minimal) phylogenetic signal, as measured by Blomberg's K. The loading vectors of these phylogenetic components or K-components can be biologically interpreted, and scatterplots of the scores can be used as a low-dimensional ordination of the data that maximally (or minimally) preserves phylogenetic signal. We present algebraic and statistical properties, along with two new summary statistics, KA and KG, of phylogenetic signal in multivariate data. Simulation studies showed that KA and KG have higher statistical power than the previously suggested statistic Kmult, especially if phylogenetic signal is low or concentrated in a few trait dimensions. In two empirical applications to vertebrate cranial shape (crocodyliforms and papionins), we found statistically significant phylogenetic signal concentrated in a few trait dimensions. The finding that phylogenetic signal can be highly variable across the dimensions of multivariate phenotypes has important implications for current maximum likelihood approaches to phylogenetic signal in multivariate data.

Identifying vital nodes for yeast network by dynamic network entropy.

BACKGROUND: The progress of the cell cycle of yeast involves the regulatory relationships between genes and the interactions proteins. However, it is still obscure which type of protein plays a decisive role in regulation and how to identify the vital nodes in the regulatory network. To elucidate the sensitive node or gene in the progression of yeast, here, we select 8 crucial regulatory factors from the yeast cell cycle to decipher a specific network and propose a simple mixed K2 algorithm to identify effectively the sensitive nodes and genes in the evolution of yeast. RESULTS: Considering the multivariate of cell cycle data, we first utilize the K2 algorithm limited to the stationary interval for the time series segmentation to measure the scores for refining the specific network. After that, we employ the network entropy to effectively screen the obtained specific network, and simulate the gene expression data by a normal distribution approximation and the screened specific network by the partial least squares method. We can conclude that the robustness of the specific network screened by network entropy is better than that of the specific network with the determined relationship by comparing the obtained specific network with the determined relationship. Finally, we can determine that the node CDH1 has the highest score in the specific network through a sensitivity score calculated by network entropy implying the gene CDH1 is the most sensitive regulatory factor. CONCLUSIONS: It is clearly of great potential value to reconstruct and visualize gene regulatory networks according to gene databases for life activities. Here, we present an available algorithm to achieve the network reconstruction by measuring the network entropy and identifying the vital nodes in the specific nodes. The results indicate that inhibiting or enhancing the expression of CDH1 can maximize the inhibition or enhancement of the yeast cell cycle. Although our algorithm is simple, it is also the first step in deciphering the profound mystery of gene regulation.

Inference of causal metabolite networks in the presence of invalid instrumental variables with GWAS summary data.

We propose structural equation models (SEMs) as a general framework to infer causal networks for metabolites and other complex traits. Traditionally SEMs are used only for individual-level data under the assumption that all instrumental variables (IVs) are valid. To overcome these limitations, we propose both one- and two-sample approaches for causal network inference based on SEMs that can: (1) perform causal analysis and discover causal relationships among multiple traits; (2) account for the possible presence of some invalid IVs; (3) allow for data analysis using only genome-wide association studies (GWAS) summary statistics when individual-level data are not available; (4) consider the possibility of bidirectional relationships between traits. Our method employs a simple stepwise selection to identify invalid IVs, thus avoiding false positives while possibly increasing true discoveries based on two-stage least squares (2SLS). We use both real GWAS data and simulated data to demonstrate the superior performance of our method over the standard 2SLS/SEMs. For real data analysis, our proposed approach is applied to a human blood metabolite GWAS summary data set to uncover putative causal relationships among the metabolites; we also identify some metabolites (putative) causal to Alzheimer's disease (AD), which, along with the inferred causal metabolite network, suggest some possible pathways of metabolites involved in AD.

Overcoming the impacts of two-step batch effect correction on gene expression estimation and inference.

Nonignorable technical variation is commonly observed across data from multiple experimental runs, platforms, or studies. These so-called batch effects can lead to difficulty in merging data from multiple sources, as they can severely bias the outcome of the analysis. Many groups have developed approaches for removing batch effects from data, usually by accommodating batch variables into the analysis (one-step correction) or by preprocessing the data prior to the formal or final analysis (two-step correction). One-step correction is often desirable due it its simplicity, but its flexibility is limited and it can be difficult to include batch variables uniformly when an analysis has multiple stages. Two-step correction allows for richer models of batch mean and variance. However, prior investigation has indicated that two-step correction can lead to incorrect statistical inference in downstream analysis. Generally speaking, two-step approaches introduce a correlation structure in the corrected data, which, if ignored, may lead to either exaggerated or diminished significance in downstream applications such as differential expression analysis. Here, we provide more intuitive and more formal evaluations of the impacts of two-step batch correction compared to existing literature. We demonstrate that the undesired impacts of two-step correction (exaggerated or diminished significance) depend on both the nature of the study design and the batch effects. We also provide strategies for overcoming these negative impacts in downstream analyses using the estimated correlation matrix of the corrected data. We compare the results of our proposed workflow with the results from other published one-step and two-step methods and show that our methods lead to more consistent false discovery controls and power of detection across a variety of batch effect scenarios. Software for our method is available through GitHub (https://github.com/jtleek/sva-devel) and will be available in future versions of the $\texttt{sva}$ R package in the Bioconductor project (https://bioconductor.org/packages/release/bioc/html/sva.html).

MuSiC2: cell-type deconvolution for multi-condition bulk RNA-seq data.

Cell-type composition of intact bulk tissues can vary across samples. Deciphering cell-type composition and its changes during disease progression is an important step toward understanding disease pathogenesis. To infer cell-type composition, existing cell-type deconvolution methods for bulk RNA sequencing (RNA-seq) data often require matched single-cell RNA-seq (scRNA-seq) data, generated from samples with similar clinical conditions, as reference. However, due to the difficulty of obtaining scRNA-seq data in diseased samples, only limited scRNA-seq data in matched disease conditions are available. Using scRNA-seq reference to deconvolve bulk RNA-seq data from samples with different disease conditions may lead to a biased estimation of cell-type proportions. To overcome this limitation, we propose an iterative estimation procedure, MuSiC2, which is an extension of MuSiC, to perform deconvolution analysis of bulk RNA-seq data generated from samples with multiple clinical conditions where at least one condition is different from that of the scRNA-seq reference. Extensive benchmark evaluations indicated that MuSiC2 improved the accuracy of cell-type proportion estimates of bulk RNA-seq samples under different conditions as compared with the traditional MuSiC deconvolution. MuSiC2 was applied to two bulk RNA-seq datasets for deconvolution analysis, including one from human pancreatic islets and the other from human retina. We show that MuSiC2 improves current deconvolution methods and provides more accurate cell-type proportion estimates when the bulk and single-cell reference differ in clinical conditions. We believe the condition-specific cell-type composition estimates from MuSiC2 will facilitate the downstream analysis and help identify cellular targets of human diseases.

Predicting miRNA-disease associations based on graph attention networks and dual Laplacian regularized least squares.

Increasing biomedical evidence has proved that the dysregulation of miRNAs is associated with human complex diseases. Identification of disease-related miRNAs is of great importance for disease prevention, diagnosis and remedy. To reduce the time and cost of biomedical experiments, there is a strong incentive to develop efficient computational methods to infer potential miRNA-disease associations. Although many computational approaches have been proposed to address this issue, the prediction accuracy needs to be further improved. In this study, we present a computational framework MKGAT to predict possible associations between miRNAs and diseases through graph attention networks (GATs) using dual Laplacian regularized least squares. We use GATs to learn embeddings of miRNAs and diseases on each layer from initial input features of known miRNA-disease associations, intra-miRNA similarities and intra-disease similarities. We then calculate kernel matrices of miRNAs and diseases based on Gaussian interaction profile (GIP) with the learned embeddings. We further fuse the kernel matrices of each layer and initial similarities with attention mechanism. Dual Laplacian regularized least squares are finally applied for new miRNA-disease association predictions with the fused miRNA and disease kernels. Compared with six state-of-the-art methods by 5-fold cross-validations, our method MKGAT receives the highest AUROC value of 0.9627 and AUPR value of 0.7372. We use MKGAT to predict related miRNAs for three cancers and discover that all the top 50 predicted results in the three diseases are confirmed by existing databases. The excellent performance indicates that MKGAT would be a useful computational tool for revealing disease-related miRNAs.

Improved motion correction in brain MRI using 3D radial trajectory and projection moment analysis.

PURPOSE: To develop a generalized rigid body motion correction method in 3D radial brain MRI to deal with continuous motion pattern through projection moment analysis. METHODS: An assumption was made that the multichannel coil moves with the head, which was achieved by using a flexible head coil. A two-step motion correction scheme was proposed to directly extract the motion parameters from the acquired k-space data using the analysis of center-of-mass with high noise robustness, which were used for retrospective motion correction. A recursive least-squares model was introduced to recursively estimate the motion parameters for every single spoke, which used the smoothness of motion and resulted in high temporal resolution and low computational cost. Five volunteers were scanned at 3 T using a 3D radial multidimensional golden-means trajectory with instructed motion patterns. The performance was tested through both simulation and in vivo experiments. Quantitative image quality metrics were calculated for comparison. RESULTS: The proposed method showed good accuracy and precision in both translation and rotation estimation. A better result was achieved using the proposed two-step correction compared to traditional one-step correction without significantly increasing computation time. Retrospective correction showed substantial improvements in image quality among all scans, even for stationary scans. CONCLUSIONS: The proposed method provides an easy, robust, and time-efficient tool for motion correction in brain MRI, which may benefit clinical diagnosis of uncooperative patients as well as scientific MRI researches.

An improved method to measure transfer functions using MRI.

PURPOSE: A previously published method for MRI-based transfer function assessment makes use of the so-called transceive phase assumption (TPA). This limits its applicability to shorter leads and/or lower field strengths. A new method is presented where the background electric field is determined from both B 1 + $$ {\mathrm{B}}_1^{+} $$ - and B 1 - $$ {\mathrm{B}}_1^{-} $$ -field distributions, avoiding the TPA and making it more generally applicable. THEORY AND METHODS: These B 1 $$ {\mathrm{B}}_1 $$ -distributions are determined from a spoiled gradient echo multiflip angle acquisition. From the separated B 1 $$ {\mathrm{B}}_1 $$ -components the background electrical field and the induced current are computed. Further improvement is achieved by recasting the B 1 $$ {\mathrm{B}}_1 $$ -field model as a "magnitude squared least squares" problem. The proposed reconstruction method is used to determine transfer functions of various copper wire lengths up to 40 cm inside an elliptical ASTM phantom. The method is first tested on EM-simulated data and subsequently phantom and bench measurements are used to determine transfer functions experimentally. RESULTS: In silica reconstructions demonstrate the validity of the proposed B 1 $$ {\mathrm{B}}_1 $$ -field model resulting in highly accurate reconstructed B 1 $$ {\mathrm{B}}_1 $$ -fields, currents, incident electric fields and transfer functions. The experimental results show slight deviations in the field model, however, resulting transfer functions are accurately determined with high similarity to simulations and comparable to bench measurements. CONCLUSION: A more generally applicable method for MRI-based transfer function assessment is presented. The proposed method circumvents phase assumptions making it applicable for longer objects and/or higher field strengths. Additional improvements are implemented in the B 1 $$ {\mathrm{B}}_1 $$ -mapping method and the solution algorithm.

Impacts of litter microbial community on litter decomposition in the absence of soil microorganisms.

What is the effect of phyllosphere microorganisms on litter decomposition in the absence of colonization by soil microorganisms? Here, we simulated the litter standing decomposition stage in the field to study the differences in the composition and structure of the phyllosphere microbial community after the mixed decomposition of Populus × canadensis and Pinus sylvestris var. mongolica litter. After 15 months of mixed decomposition, we discovered that litters that were not in contact with soil had an antagonistic effect (the actual decomposition rate was 18.18%, which is lower than the expected decomposition rate) and the difference between the litters themselves resulted in a negative response to litter decomposition. In addition, there was no significant difference in bacterial and fungal community diversity after litter decomposition. The litter bacterial community was negatively responsive to litter properties and positively responsive to the fungal community. Importantly, we found that bacterial communities had a greater impact on litter decomposition than fungi. This study has enriched our understanding of the decomposition of litter itself and provided a theoretical basis for further exploring the "additive and non-additive effects" of litter decomposition and the mechanism of microbial drive. IMPORTANCE: The study of litter decomposition mechanism plays an important role in the material circulation of the global ecosystem. However, previous studies have often looked at contact with soil as the starting point for decomposition. But actually, standing litter is very common in forest ecosystems. Therefore, we used field simulation experiments to simulate the decomposition of litters without contact with soil for 15 months, to explore the combined and non-added benefits of the decomposition of mixed litters, and to study the influence of microbial community composition on the decomposition rate while comparing the differences of microbial communities.

Unveiling the transferability of PLSR models for leaf trait estimation: lessons from a comprehensive analysis with a novel global dataset.

Leaf traits are essential for understanding many physiological and ecological processes. Partial least squares regression (PLSR) models with leaf spectroscopy are widely applied for trait estimation, but their transferability across space, time, and plant functional types (PFTs) remains unclear. We compiled a novel dataset of paired leaf traits and spectra, with 47 393 records for > 700 species and eight PFTs at 101 globally distributed locations across multiple seasons. Using this dataset, we conducted an unprecedented comprehensive analysis to assess the transferability of PLSR models in estimating leaf traits. While PLSR models demonstrate commendable performance in predicting chlorophyll content, carotenoid, leaf water, and leaf mass per area prediction within their training data space, their efficacy diminishes when extrapolating to new contexts. Specifically, extrapolating to locations, seasons, and PFTs beyond the training data leads to reduced R2 (0.12-0.49, 0.15-0.42, and 0.25-0.56) and increased NRMSE (3.58-18.24%, 6.27-11.55%, and 7.0-33.12%) compared with nonspatial random cross-validation. The results underscore the importance of incorporating greater spectral diversity in model training to boost its transferability. These findings highlight potential errors in estimating leaf traits across large spatial domains, diverse PFTs, and time due to biased validation schemes, and provide guidance for future field sampling strategies and remote sensing applications.

Postprandial Plasma and Whole Blood Amino Acids Are Largely Indicative of Dietary Amino Acids in Adult Dogs Consuming Diets with Increasing Whole Pulse Ingredient Inclusion.

BACKGROUND: Pulse ingredients often replace grains in grain-free dog diets owing to their high-protein content. However, research to ascertain the benefit of this modification is limited. OBJECTIVES: This study aimed to correlate food compounds in 1 corn-inclusive control diet and 3 grain-free diets with increasing inclusions of whole pulses (≤45%; Pulse15, Pulse30, and Pulse45), formulated to meet similar macronutrient and micronutrient targets with postprandial amino acids (AAs) in healthy dogs >20 wk. METHODS: Diets were analyzed for biochemical compounds using tandem mass spectrometry. Twenty-eight outdoor-housed, healthy, adult Siberian Huskies were allocated to diet, and meal responses were analyzed at baseline and weeks 2, 4, 8, 16, and 20 with samples collected at fasted and 15, 30, 60, 90, 120, and 180 min after meal presentation. Blood AAs were analyzed by ultra performance liquid chromatography and differences across week, treatment, and time postmeal were analyzed in SAS Studio. Partial least squares regression was performed in SAS Studio using biochemical compounds in the diet as predictor variables and blood AAs as response variables. RESULTS: In plasma, Pulse45 had ∼32% greater postprandial Asn than Pulse15, and the control diet had ∼34% greater postprandial Leu and ∼35% greater Pro than Pulse15 (P < 0.05). In whole blood, Pulse30 had ∼23% greater postprandial Lys than the control diet, and the control diet had ∼21% greater postprandial Met and ∼18% greater Pro than Pulse45 and Pulse30, respectively (P < 0.05). Several phospholipids were correlated with postprandial AAs. Compounds in the urea cycle and glycine and serine metabolism were more enriched (P < 0.05) in plasma and whole blood, respectively. CONCLUSIONS: In macronutrient-balanced and micronutrient-balanced canine diets that differ in their inclusion of corn-derived compared with pulse-derived ingredients, postprandial changes in circulating AAs are largely indicative of the dietary AAs. This helps further our understanding of AA metabolism in healthy dogs fed grain-free diets.

Development and Validation of SCACOMS, a Composite Scale for Assessing Disease Progression and Treatment Effects in Spinocerebellar Ataxia.

Spinocerebellar ataxias (SCA) are rare inherited neurodegenerative disorders characterized by a progressive impairment of gait, balance, limb coordination, and speech. There is currently no composite scale that includes multiple aspects of the SCA experience to assess disease progression and treatment effects. Applying the method of partial least squares (PLS) regression, we developed the Spinocerebellar Ataxia Composite Scale (SCACOMS) from two SCA natural history datasets (NCT01060371, NCT02440763). PLS regression selected items based on their ability to detect clinical decline, with optimized weights based on the item's degree of progression. Following model validation, SCACOMS was leveraged to examine disease progression and treatment effects in a 48-week SCA clinical trial cohort (NCT03701399). Items from the Clinical Global Impression-Global Improvement Scale (CGI-I), the Friedreich Ataxia Rating Scale (FARS) - functional stage, and the Modified Functional Scale for the Assessment and Rating of Ataxia (f-SARA) were objectively selected with weightings based on their sensitivity to clinical decline. The resulting SCACOMS exhibited improved sensitivity to disease progression and greater treatment effects (compared to the original scales from which they were derived) in a 48-week clinical trial of a novel therapeutic agent. The trial analyses also provided a SCACOMS-derived estimate of the temporal delay in SCA disease progression. SCACOMS is a useful composite measure, effectively capturing disease progression and highlighting treatment effects in patients with SCA. SCACOMS will be a powerful tool in future studies given its sensitivity to clinical decline and ability to detect a meaningful clinical impact of disease-modifying treatments.

Hybrid 2D/3D-quantitative structure-activity relationship studies on the bioactivities and molecular mechanism of antibacterial peptides.

Antimicrobial peptide (AMP) is the polypeptide, which protects the organism avoiding attack from pathogenic bacteria. Studies have shown that there were some antimicrobial peptides with molecular action mechanism involved in crossing the cell membrane without inducing severe membrane collapse, then interacting with cytoplasmic target-nucleic acid, and exerting antibacterial activity by interfacing the transmission of genetic information of pathogenic microorganisms. However, the relationship between the antibacterial activities and peptide structures was still unclear. Therefore, in the present work, a series of AMPs with a sequence of 20 amino acids was extracted from DBAASP database, then, quantitative structure-activity relationship (QSAR) methods were conducted on these peptides. In addition, novel antimicrobial peptides with  stronger antimicrobial activities were designed according to the information originated from the constructed models. Hence, the outcome of this study would lay a solid foundation for the in-silico design and exploration of novel antibacterial peptides with improved activity activities.

In Situ Eutectic Formation in a Polymeric Matrix via Hot-Melt Reactive Extrusion and the Use of Partial Least Squares Regression Modeling for Reaction Yield Determination.

There has been a significant volume of work investigating the design and synthesis of new crystalline multicomponent systems via examining complementary functional groups that can reliably interact through the formation of noncovalent bonds, such as hydrogen bonds (H-bonds). Crystalline multicomponent molecular adducts formed using this approach, such as cocrystals, salts, and eutectics, have emerged as drug product intermediates that can lead to effective drug property modifications. Recent advancement in the production for these multicomponent molecular adducts has moved from batch techniques that rely upon intensive solvent use to those that are solvent-free, continuous, and industry-ready, such as reactive extrusion. In this study, a novel eutectic system was found when processing albendazole and maleic acid at a 1:2 molar ratio and successfully prepared using mechanochemical methods including liquid-assisted grinding and hot-melt reactive extrusion. The produced eutectic was characterized to exhibit a 100 °C reduction in melting temperature and enhanced dissolution performance (>12-fold increase at 2 h point), when compared to the native drug compound. To remove handling of the eutectic as a formulation intermediate, an end-to-end continuous-manufacturing-ready process enables feeding of the raw parent reagents in their respective natural forms along with a chosen polymeric excipient, Eudragit EPO. The formation of the eutectic was confirmed to have taken place in situ in the presence of the polymer, with the reaction yield determined using a multivariate calibration model constructed by combining spectroscopic analysis with partial least-squares regression modeling. The ternary extrudates exhibited a dissolution profile similar to that of the 1:2 prepared eutectic, suggesting a physical distribution (or suspension) of the in situ synthesized eutectic contents within the polymeric matrix.

Looking beyond the mean: quantile regression for comparative physiologists.

Statistical analyses that physiologists use to test hypotheses predominantly centre on means, but the tail ends of the response distribution can behave quite differently and underpin important scientific phenomena. We demonstrate that quantile regression (QR) offers a way to bypass some limitations of least squares regression (LSR) by building a picture of independent variable effects across the whole distribution of a dependent variable. We used LSR and QR with simulated and real datasets. With simulated data, LSR showed no change in the mean response but missed significant effects in the tails of the distribution found using QR. With real data, LSR showed a significant change in the mean response but missed a lack of response in the upper quantiles which was biologically revealing. Together, this highlights that QR can help to ask and answer more questions about variation in nature.

Statistical inference with large-scale trait imputation.

Recently a nonparametric method called LS-imputation has been proposed for large-scale trait imputation based on a GWAS summary dataset and a large set of genotyped individuals. The imputed trait values, along with the genotypes, can be treated as an individual-level dataset for downstream genetic analyses, including those that cannot be done with GWAS summary data. However, since the covariance matrix of the imputed trait values is often too large to calculate, the current method imposes a working assumption that the imputed trait values are identically and independently distributed, which is incorrect in truth. Here we propose a "divide and conquer/combine" strategy to estimate and account for the covariance matrix of the imputed trait values via batches, thus relaxing the incorrect working assumption. Applications of the methods to the UK Biobank data for marginal association analysis showed some improvement by the new method in some cases, but overall the original method performed well, which was explained by nearly constant variances of and mostly weak correlations among imputed trait values.

Multiblock partial least squares and rank aggregation: Applications to detection of bacteriophages associated with antimicrobial resistance in the presence of potential confounding factors.

Urban environments, characterized by bustling mass transit systems and high population density, host a complex web of microorganisms that impact microbial interactions. These urban microbiomes, influenced by diverse demographics and constant human movement, are vital for understanding microbial dynamics. We explore urban metagenomics, utilizing an extensive dataset from the Metagenomics & Metadesign of Subways & Urban Biomes (MetaSUB) consortium, and investigate antimicrobial resistance (AMR) patterns. In this pioneering research, we delve into the role of bacteriophages, or "phages"-viruses that prey on bacteria and can facilitate the exchange of antibiotic resistance genes (ARGs) through mechanisms like horizontal gene transfer (HGT). Despite their potential significance, existing literature lacks a consensus on their significance in ARG dissemination. We argue that they are an important consideration. We uncover that environmental variables, such as those on climate, demographics, and landscape, can obscure phage-resistome relationships. We adjust for these potential confounders and clarify these relationships across specific and overall antibiotic classes with precision, identifying several key phages. Leveraging machine learning tools and validating findings through clinical literature, we uncover novel associations, adding valuable insights to our comprehension of AMR development.

Machine Learning Demonstrates Dominance of Physical Characteristics over Particle Composition in Coal Dust Toxicity.

Mine dust has been linked to the development of pneumoconiotic diseases such as silicosis and coal workers' pneumoconiosis. Currently, it is understood that the physicochemical and mineralogical characteristics drive the toxic nature of dust particles; however, it remains unclear which parameter(s) account for the differential toxicity of coal dust. This study aims to address this issue by demonstrating the use of the partial least squares regression (PLSR) machine learning approach to compare the influence of D50 sub 10 µm coal particle characteristics against markers of cellular damage. The resulting analysis of 72 particle characteristics against cytotoxicity and lipid peroxidation reflects the power of PLSR as a tool to elucidate complex particle-cell relationships. By comparing the relative influence of each characteristic within the model, the results reflect that physical characteristics such as shape and particle roughness may have a greater impact on cytotoxicity and lipid peroxidation than composition-based parameters. These results present the first multivariate assessment of a broad-spectrum data set of coal dust characteristics using latent structures to assess the relative influence of particle characteristics on cellular damage.