A Penalized Multi-trait Mixed Model for Association Mapping in Pedigree-based GWAS

In genome-wide association studies (GWAS), penalization is an important approach for identifying genetic markers associated with trait while mixed model is successful in accounting for a complicated dependence structure among samples. Therefore, penalized linear mixed model is a tool that combines the advantages of penalization approach and linear mixed model. In this study, a GWAS with multiple highly correlated traits is analyzed. For GWAS with multiple quantitative traits that are highly correlated, the analysis using traits marginally inevitably lose some essential information among multiple traits. We propose a penalized-MTMM, a penalized multivariate linear mixed model that allows both the within-trait and between-trait variance components simultaneously for multiple traits. The proposed penalized-MTMM estimates variance components using an AI-REML method and conducts variable selection and point estimation simultaneously using group MCP and sparse group MCP. Best linear unbiased predictor (BLUP) is used to find predictive values and the Pearson's correlations between predictive values and their corresponding observations are used to evaluate prediction performance. Both prediction and selection performance of the proposed approach and its comparison with the uni-trait penalized-LMM are evaluated through simulation studies. We apply the proposed approach to a GWAS data from Genetic Analysis Workshop (GAW) 18

A mendelian randomization study investigates the causal relationship between immune cell phenotypes and cerebral aneurysm

Background: Cerebral aneurysms (CAs) are a significant cerebrovascular ailment with a multifaceted etiology influenced by various factors including heredity and environment. This study aimed to explore the possible link between different types of immune cells and the occurrence of CAs.Methods: We analyzed the connection between 731 immune cell signatures and the risk of CAs by using publicly available genetic data. The analysis included four immune features, specifically median brightness levels (MBL), proportionate cell (PC), definite cell (DC), and morphological attributes (MA). Mendelian randomization (MR) analysis was conducted using the instrumental variables (IVs) derived from the genetic variation linked to CAs.Results: After multiple test adjustment based on the FDR method, the inverse variance weighted (IVW) method revealed that 3 immune cell phenotypes were linked to the risk of CAs. These included CD45 on HLA DR+NK (odds ratio (OR), 1.116; 95% confidence interval (CI), 1.001–1.244; p = 0.0489), CX3CR1 on CD14− CD16− (OR, 0.973; 95% CI, 0.948–0.999; p = 0.0447). An immune cell phenotype CD16− CD56 on NK was found to have a significant association with the risk of CAs in reverse MR study (OR, 0.950; 95% CI, 0.911–0.990; p = 0.0156).Conclusion: Our investigation has yielded findings that support a substantial genetic link between immune cells and CAs, thereby suggesting possible implications for future clinical interventions

Identification and validation of an immune-related gene prognostic signature for clear cell renal carcinoma

Clear Cell Renal Carcinoma (ccRCC) accounts for nearly 80% of renal carcinoma cases, and immunotherapy plays an important role in ccRCC therapy. However, the responses to immunotherapy and overall survival for ccRCC patients are still hard to predict. Here, we constructed an immune-related predictive signature using 19 genes based on TCGA datasets. We also analyzed its relationships between disease prognosis, infiltrating immune cells, immune subtypes, mutation load, immune dysfunction, immune escape, etc. We found that our signature can distinguish immune characteristics and predict immunotherapeutic response for ccRCC patients with better prognostic prediction value than other immune scores. The expression levels of prognostic genes were determined by RT-qPCR assay. This signature may help to predict overall survival and guide the treatment for patients with ccRCC

Microbial carbon use efficiency promotes global soil carbon storage

Soils store more carbon than other terrestrial ecosystems$^{1,2}$. How soil organic carbon (SOC) forms and persists remains uncertain$^{1,3}$, which makes it challenging to understand how it will respond to climatic change$^{3,4}$. It has been suggested that soil microorganisms play an important role in SOC formation, preservation and loss$^{5–7}$. Although microorganisms affect the accumulation and loss of soil organic matter through many pathways$^{4,6,8–11}$, microbial carbon use efficiency (CUE) is an integrative metric that can capture the balance of these processes$^{12,13}$. Although CUE has the potential to act as a predictor of variation in SOC storage, the role of CUE in SOC persistence remains unresolved$^{7,14,15}$. Here we examine the relationship between CUE and the preservation of SOC, and interactions with climate, vegetation and edaphic properties, using a combination of global-scale datasets, a microbial-process explicit model, data assimilation, deep learning and meta-analysis. We find that CUE is at least four times as important as other evaluated factors, such as carbon input, decomposition or vertical transport, in determining SOC storage and its spatial variation across the globe. In addition, CUE shows a positive correlation with SOC content. Our findings point to microbial CUE as a major determinant of global SOC storage. Understanding the microbial processes underlying CUE and their environmental dependence may help the prediction of SOC feedback to a changing climate

Microbial carbon use efficiency promotes global soil carbon storage

Funding Information: We thank H. Yang, M. Schrumpf, T. Wutzler, R. Zheng and H. Ma for their comments and suggestions on this study. This work was supported by the National Natural Science Foundation of China (42125503) and the National Key Research and Development Program of China (2020YFA0608000, 2020YFA0607900 and 2021YFC3101600). F.T. was financially supported by China Scholarship Council during his visit at Food and Agricultural Organization of the United Nations (201906210489) and the Max-Planck Institute for Biogeochemistry (202006210289). The contributions of Y.L. were supported through US National Science Foundation DEB 1655499 and 2242034, subcontract CW39470 from Oak Ridge National Laboratory (ORNL) to Cornell University, DOE De-SC0023514, and the USDA National Institute of Food and Agriculture. S.M. has received funding from the ERC under the European Union’s H2020 Research and Innovation Programme (101001608). The contributions of U.M. were supported through a US Department of Energy grant to the Sandia National Laboratories, which is a multi-mission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International, Inc., for the US Department of Energy’s National Nuclear Security Administration under contract DE-NA-0003525. We thank the WoSIS database ( https://www.isric.org/explore/wosis ) for providing the publicly available global-scale SOC database used in this study. Publisher Copyright: © 2023, The Author(s).Peer reviewedPublisher PD

Transient dynamics of terrestrial carbon storage : mathematical foundation and its applications

Terrestrial ecosystems have absorbed roughly 30 % of anthropogenic CO2 emissions over the past decades, but it is unclear whether this carbon (C) sink will endure into the future. Despite extensive modeling and experimental and observational studies, what fundamentally determines transient dynamics of terrestrial C storage under global change is still not very clear. Here we develop a new framework for understanding transient dynamics of terrestrial C storage through mathematical analysis and numerical experiments. Our analysis indicates that the ultimate force driving ecosystem C storage change is the C storage capacity, which is jointly determined by ecosystem C input (e.g., net primary production, NPP) and residence time. Since both C input and residence time vary with time, the C storage capacity is time-dependent and acts as a moving attractor that actual C storage chases. The rate of change in C storage is proportional to the C storage potential, which is the difference between the current storage and the storage capacity. The C storage capacity represents instantaneous responses of the land C cycle to external forcing, whereas the C storage potential represents the internal capability of the land C cycle to influence the C change trajectory in the next time step. The influence happens through redistribution of net C pool changes in a network of pools with different residence times. Moreover, this and our other studies have demonstrated that one matrix equation can replicate simulations of most land C cycle models (i.e., physical emulators). As a result, simulation outputs of those models can be placed into a three-dimensional (3-D) parameter space to measure their differences. The latter can be decomposed into traceable components to track the origins of model uncertainty. In addition, the physical emulators make data assimilation computationally feasible so that both C flux- and pool-related datasets can be used to better constrain model predictions of land C sequestration. Overall, this new mathematical framework offers new approaches to understanding, evaluating, diagnosing, and improving land C cycle models.This work was partially done through the working group, Nonautonomous Systems and Terrestrial Carbon Cycle, at the National Institute for Mathematical and Biological Synthesis, an institute sponsored by the National Science Foundation, the US Departmernt of Homeland Security, and the US Department of Agriculture through NSF award no. EF-0832858, with additional support from the University of Tennessee, Knoxville, Research in Yiqi Luo EcoLab was financially supported by US Department of Energy grants DE-SC0008270, DE-SC0014085, and US National Science Foundation (NSF) grants EF 1137293 and OIA-1301789.Ye

The immune enhancing effect of antimicrobial peptide LLv on broilers chickens

ABSTRACT: To evaluate the effect and its mechanism of heat-resistant antimicrobial peptide LLv on broilers, three hundred 1-day-old healthy AA+ female broilers were allocated into 5 groups with 6 replicates in each group and 10 birds in each replicate. Birds were given a basal diet, an antibiotic diet (10.2 mg/kg chlortetracycline hydrochloride), and the basal diet supplemented with 10, 50, and 100 mg/kg LLv for 42 d, respectively. Compared with the group which birds were fed an antibiotic-free basal diet (control group), supplementing 100 mg/kg LLv increased 21-day IgA, IgM, IL-4, AIV-Ab, IFN-γ levels and 42-day IgA, IgM, IL-4, AIV-Ab levels and reduced 42-day IL-1 levels in serum (P 0.05). Compared with antibiotic group, the 10 mg/kg LLv reduced 21-day sIgA content and the 50 mg/kg LLv reduced 42-d the expression rate of sIgA secretory cells in jejunal mucosa (P < 0.05). Compared with control group, the 100 mg/kg LLv increased the expression of TCR, IL-15, CD28, BAFF, CD86, CD83, MHC-II, and CD40 genes in jejunal mucosa at 21 d and 42 d (P < 0.05). Compared with antibiotic group, the 100 mg/kg LLv increased the expression of 21-day BAFF, CD40, MHC-II, CD83 genes and the expression of 42-day BAFF, TCR, IL-15, CD40, CD83 genes in jejunal mucosa (P < 0.05). The results showed that the addition of LLv to the ration had a promotional effect on the immune function of broiler chickens