Systems genetics approaches to understand complex traits

From molecular mechanisms of prostate cancer to translational applications: based on multi-omics fusion analysis and intelligent medicine

Prostate cancer is the most common cancer in men worldwide and has a high mortality rate. The complex and heterogeneous development of prostate cancer has become a core obstacle in the treatment of prostate cancer. Simultaneously, the issues of overtreatment in early-stage diagnosis, oligometastasis and dormant tumor recognition, as well as personalized drug utilization, are also specific concerns that require attention in the clinical management of prostate cancer. Some typical genetic mutations have been proved to be associated with prostate cancer's initiation and progression. However, single-omic studies usually are not able to explain the causal relationship between molecular alterations and clinical phenotypes. Exploration from a systems genetics perspective is also lacking in this field, that is, the impact of gene network, the environmental factors, and even lifestyle behaviors on disease progression. At the meantime, current trend emphasizes the utilization of artificial intelligence (AI) and machine learning techniques to process extensive multidimensional data, including multi-omics. These technologies unveil the potential patterns, correlations, and insights related to diseases, thereby aiding the interpretable clinical decision making and applications, namely intelligent medicine. Therefore, there is a pressing need to integrate multidimensional data for identification of molecular subtypes, prediction of cancer progression and aggressiveness, along with perosonalized treatment performing. In this review, we systematically elaborated the landscape from molecular mechanism discovery of prostate cancer to clinical translational applications. We discussed the molecular profiles and clinical manifestations of prostate cancer heterogeneity, the identification of different states of prostate cancer, as well as corresponding precision medicine practices. Taking multi-omics fusion, systems genetics, and intelligence medicine as the main perspectives, the current research results and knowledge-driven research path of prostate cancer were summarized.

Systems genetics analysis of human body fat distribution genes identifies adipocyte processes

Excess abdominal fat is a sexually dimorphic risk factor for cardio-metabolic disease and is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Whereas this trait is highly heritable, few causal genes are known. We aimed to identify novel drivers of WHRadjBMI using systems genetics. We used two independent cohorts of adipose tissue gene expression and constructed sex- and depot-specific Bayesian networks to model gene-gene interactions from 8,492 genes. Using key driver analysis, we identified genes that, in silico and putatively in vitro, regulate many others. 51-119 key drivers in each network were replicated in both cohorts. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We overexpressed or down-regulated seven key driver genes in human subcutaneous pre-adipocytes. Key driver genes ANAPC2 and RSPO1 inhibited adipogenesis, whereas PSME3 increased adipogenesis. RSPO1 increased Wnt signaling activity. In differentiated adipocytes, MIGA1 and UBR1 down-regulation led to mitochondrial dysfunction. These five genes regulate adipocyte function, and we hypothesize that they regulate fat distribution.

Mitochondrial disorders leading to Alzheimer's disease-perspectives of diagnosis and treatment

Alzheimer's disease (AD) is a neurodegenerative disorder and the most common cause of dementia globally. The pathogenesis of AD remains still unclear. The three main features of AD are extracellular deposits of amyloid beta (Aβ) plaque, accumulation of abnormal formation hyper-phosphorylated tau protein, and neuronal loss. Mitochondrial impairment plays an important role in the pathogenesis of AD. There are problems with decreased activity of multiple complexes, disturbed mitochondrial fusion, and fission or formation of reactive oxygen species (ROS). Moreover, mitochondrial transport is impaired in AD. Mouse models in many research show disruptions in anterograde and retrograde transport. Both mitochondrial transportation and network impairment have a huge impact on synapse loss and, as a result, cognitive impairment. One of the very serious problems in AD is also disruption of insulin signaling which impairs mitochondrial Aβ removal.Discovering precise mechanisms leading to AD enables us to find new treatment possibilities. Recent studies indicate the positive influence of metformin or antioxidants such as MitoQ, SS-31, SkQ, MitoApo, MitoTEMPO, and MitoVitE on mitochondrial functioning and hence prevent cognitive decline. Impairments in mitochondrial fission may be treated with mitochondrial division inhibitor-1 or ceramide.

Single-Cell Gene-Regulatory Networks of Advanced Symptomatic Atherosclerosis

BACKGROUND

While our understanding of the single-cell gene expression patterns underlying the transformation of vascular cell types during the progression of atherosclerosis is rapidly improving, the clinical and pathophysiological relevance of these changes remains poorly understood.

METHODS

Single-cell RNA sequencing data generated with SmartSeq2 (≈8000 genes/cell) in 16 588 single cells isolated during atherosclerosis progression in Ldlr-/-Apob100/100 mice with human-like plasma lipoproteins and from humans with asymptomatic and symptomatic carotid plaques was clustered into multiple subtypes. For clinical and pathophysiological context, the advanced-stage and symptomatic subtype clusters were integrated with 135 tissue-specific (atherosclerotic aortic wall, mammary artery, liver, skeletal muscle, and visceral and subcutaneous, fat) gene-regulatory networks (GRNs) inferred from 600 coronary artery disease patients in the STARNET (Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task) study.

RESULTS

Advanced stages of atherosclerosis progression and symptomatic carotid plaques were largely characterized by 3 smooth muscle cells (SMCs), and 3 macrophage subtype clusters with extracellular matrix organization/osteogenic (SMC), and M1-type proinflammatory/Trem2-high lipid-associated (macrophage) phenotypes. Integrative analysis of these 6 clusters with STARNET revealed significant enrichments of 3 arterial wall GRNs: GRN33 (macrophage), GRN39 (SMC), and GRN122 (macrophage) with major contributions to coronary artery disease heritability and strong associations with clinical scores of coronary atherosclerosis severity. The presence and pathophysiological relevance of GRN39 were verified in 5 independent RNAseq data sets obtained from the human coronary and aortic artery, and primary SMCs and by targeting its top-key drivers, FRZB and ALCAM in cultured human coronary artery SMCs.

CONCLUSIONS

By identifying and integrating the most gene-rich single-cell subclusters of atherosclerosis to date with a coronary artery disease framework of GRNs, GRN39 was identified and independently validated as being critical for the transformation of contractile SMCs into an osteogenic phenotype promoting advanced, symptomatic atherosclerosis.

On knowing a gene: A distributional hypothesis of gene function

As words can have multiple meanings that depend on sentence context, genes can have various functions that depend on the surrounding biological system. This pleiotropic nature of gene function is limited by ontologies, which annotate gene functions without considering biological contexts. We contend that the gene function problem in genetics may be informed by recent technological leaps in natural language processing, in which representations of word semantics can be automatically learned from diverse language contexts. In contrast to efforts to model semantics as "is-a" relationships in the 1990s, modern distributional semantics represents words as vectors in a learned semantic space and fuels current advances in transformer-based models such as large language models and generative pre-trained transformers. A similar shift in thinking of gene functions as distributions over cellular contexts may enable a similar breakthrough in data-driven learning from large biological datasets to inform gene function.

Macrophage transplantation rescues RNASET2-deficient leukodystrophy by replacing deficient microglia in a zebrafish model

RNASET2-deficient leukodystrophy is a rare infantile white matter disorder mimicking a viral infection and resulting in severe psychomotor impairments. Despite its severity, there is little understanding of cellular mechanisms of pathogenesis and no treatments. Recent research using the rnaset2 mutant zebrafish model has suggested that microglia may be the drivers of the neuropathology, due to their failure to digest apoptotic debris during neurodevelopment. Therefore, we developed a strategy for microglial replacement through transplantation of adult whole kidney marrow-derived macrophages into embryonic hosts. Using live imaging, we revealed that transplant-derived macrophages can engraft within host brains and express microglia-specific markers, suggesting the adoption of a microglial phenotype. Tissue-clearing strategies revealed the persistence of transplanted cells in host brains beyond embryonic stages. We demonstrated that transplanted cells clear apoptotic cells within the brain, as well as rescue overactivation of the antiviral response otherwise seen in mutant larvae. RNA sequencing at the point of peak transplant-derived cell engraftment confirms that transplantation can reduce the brain-wide immune response and particularly, the antiviral response, in rnaset2-deficient brains. Crucially, this reduction in neuroinflammation resulted in behavioral rescue-restoring rnaset2 mutant motor activity to wild-type (WT) levels in embryonic and juvenile stages. Together, these findings demonstrate the role of microglia as the cellular drivers of neuropathology in rnaset2 mutants and that macrophage transplantation is a viable strategy for microglial replacement in the zebrafish. Therefore, microglia-targeted interventions may have therapeutic benefits in RNASET2-deficient leukodystrophy.

Hepatokine ITIH3 protects against hepatic steatosis by downregulating mitochondrial bioenergetics and de novo lipogenesis

Recent studies demonstrate that liver secretory proteins, also known as hepatokines, regulate normal development, obesity, and simple steatosis to non-alcoholic steatohepatitis (NASH) progression. Using a panel of ∼100 diverse inbred strains of mice and a cohort of bariatric surgery patients, we found that one such hepatokine, inter-trypsin inhibitor heavy chain 3 (ITIH3), was progressively lower in severe non-alcoholic fatty liver disease (NAFLD) disease states highlighting an inverse relationship between Itih3/ITIH3 expression and NAFLD severity. Follow-up animal and cell culture models demonstrated that hepatic ITIH3 overexpression lowered liver triglyceride and lipid droplet accumulation, respectively. Conversely, ITIH3 knockdown in mice increased the liver triglyceride in two independent NAFLD models. Mechanistically, ITIH3 reduced mitochondrial respiration and this, in turn, reduced liver triglycerides, via downregulated de novo lipogenesis. This was accompanied by increased STAT1 signaling and Stat3 expression, both of which are known to protect against NAFLD/NASH. Our findings indicate hepatokine ITIH3 as a potential biomarker and/or treatment for NAFLD.

Systems Approaches to Cell Culture-Derived Extracellular Vesicles for Acute Kidney Injury Therapy: Prospects and Challenges

Acute kidney injury (AKI) is a heterogeneous syndrome, comprising diverse etiologies of kidney insults that result in high mortality and morbidity if not well managed. Although great efforts have been made to investigate underlying pathogenic mechanisms of AKI, there are limited therapeutic strategies available. Extracellular vesicles (EV) are membrane-bound vesicles secreted by various cell types, which can serve as cell-free therapy through transfer of bioactive molecules. In this review, we first overview the AKI syndrome and EV biology, with a particular focus on the technical aspects and therapeutic application of cell culture-derived EVs. Second, we illustrate how multi-omic approaches to EV miRNA, protein, and genomic cargo analysis can yield new insights into their mechanisms of action and address unresolved questions in the field. We then summarize major experimental evidence regarding the therapeutic potential of EVs in AKI, which we subdivide into stem cell and non-stem cell-derived EVs. Finally, we highlight the challenges and opportunities related to the clinical translation of animal studies into human patients.

Informing immunotherapy with multi-omics driven machine learning

Progress in sequencing technologies and clinical experiments has revolutionized immunotherapy on solid and hematologic malignancies. However, the benefits of immunotherapy are limited to specific patient subsets, posing challenges for broader application. To improve its effectiveness, identifying biomarkers that can predict patient response is crucial. Machine learning (ML) play a pivotal role in harnessing multi-omic cancer datasets and unlocking new insights into immunotherapy. This review provides an overview of cutting-edge ML models applied in omics data for immunotherapy analysis, including immunotherapy response prediction and immunotherapy-relevant tumor microenvironment identification. We elucidate how ML leverages diverse data types to identify significant biomarkers, enhance our understanding of immunotherapy mechanisms, and optimize decision-making process. Additionally, we discuss current limitations and challenges of ML in this rapidly evolving field. Finally, we outline future directions aimed at overcoming these barriers and improving the efficiency of ML in immunotherapy research.

Combinatorial in vivo genome editing identifies widespread epistasis during lung tumorigenesis

Lung adenocarcinoma, the most common subtype of lung cancer, is genomically complex, with tumors containing tens to hundreds of non-synonymous mutations. However, little is understood about how genes interact with each other to enable tumorigenesis in vivo , largely due to a lack of methods for investigating genetic interactions in a high-throughput and multiplexed manner. Here, we employed a novel platform to generate tumors with all pairwise inactivation of ten tumor suppressor genes within an autochthonous mouse model of oncogenic KRAS-driven lung cancer. By quantifying the fitness of tumors with every single and double mutant genotype, we show that most tumor suppressor genetic interactions exhibited negative epistasis, with diminishing returns on tumor fitness. In contrast, Apc inactivation showed positive epistasis with the inactivation of several other genes, including dramatically synergistic effects on tumor fitness in combination with Lkb1 or Nf1 inactivation. This approach has the potential to expand the scope of genetic interactions that may be functionally characterized in vivo , which could lead to a better understanding of how complex tumor genotypes impact each step of carcinogenesis.

Changes in environmental exposures over decades may influence the genetic architecture of severe spermatogenic failure

STUDY QUESTION

Do the genetic determinants of idiopathic severe spermatogenic failure (SPGF) differ between generations?

SUMMARY ANSWER

Our data support that the genetic component of idiopathic SPGF is impacted by dynamic changes in environmental exposures over decades.

WHAT IS KNOWN ALREADY

The idiopathic form of SPGF has a multifactorial etiology wherein an interaction between genetic, epigenetic, and environmental factors leads to the disease onset and progression. At the genetic level, genome-wide association studies (GWASs) allow the analysis of millions of genetic variants across the genome in a hypothesis-free manner, as a valuable tool for identifying susceptibility risk loci. However, little is known about the specific role of non-genetic factors and their influence on the genetic determinants in this type of conditions.

STUDY DESIGN, SIZE, DURATION

Case-control genetic association analyses were performed including a total of 912 SPGF cases and 1360 unaffected controls.

PARTICIPANTS/MATERIALS, SETTING, METHODS

All participants had European ancestry (Iberian and German). SPGF cases were diagnosed during the last decade either with idiopathic non-obstructive azoospermia (n = 547) or with idiopathic non-obstructive oligozoospermia (n = 365). Case-control genetic association analyses were performed by logistic regression models considering the generation as a covariate and by in silico functional characterization of the susceptibility genomic regions.

MAIN RESULTS AND THE ROLE OF CHANCE

This analysis revealed 13 novel genetic association signals with SPGF, with eight of them being independent. The observed associations were mostly explained by the interaction between each lead variant and the age-group. Additionally, we established links between these loci and diverse non-genetic factors, such as toxic or dietary habits, respiratory disorders, and autoimmune diseases, which might potentially influence the genetic architecture of idiopathic SPGF.

LARGE SCALE DATA

GWAS data are available from the authors upon reasonable request.

LIMITATIONS, REASONS FOR CAUTION

Additional independent studies involving large cohorts in ethnically diverse populations are warranted to confirm our findings.

WIDER IMPLICATIONS OF THE FINDINGS

Overall, this study proposes an innovative strategy to achieve a more precise understanding of conditions such as SPGF by considering the interactions between a variable exposome through different generations and genetic predisposition to complex diseases.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the "Plan Andaluz de Investigación, Desarrollo e Innovación (PAIDI 2020)" (ref. PY20_00212, P20_00583), the Spanish Ministry of Economy and Competitiveness through the Spanish National Plan for Scientific and Technical Research and Innovation (ref. PID2020-120157RB-I00 funded by MCIN/ AEI/10.13039/501100011033), and the 'Proyectos I+D+i del Programa Operativo FEDER 2020' (ref. B-CTS-584-UGR20). ToxOmics-Centre for Toxicogenomics and Human Health, Genetics, Oncology and Human Toxicology, is also partially supported by the Portuguese Foundation for Science and Technology (Projects: UIDB/00009/2020; UIDP/00009/2020). The authors declare no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

The genetic and dietary landscape of the muscle insulin signalling network

Metabolic disease is caused by a combination of genetic and environmental factors, yet few studies have examined how these factors influence signal transduction, a key mediator of metabolism. Using mass spectrometry-based phosphoproteomics, we quantified 23,126 phosphosites in skeletal muscle of five genetically distinct mouse strains in two dietary environments, with and without acute in vivo insulin stimulation. Almost half of the insulin-regulated phosphoproteome was modified by genetic background on an ordinary diet, and high-fat high-sugar feeding affected insulin signalling in a strain-dependent manner. Our data revealed coregulated subnetworks within the insulin signalling pathway, expanding our understanding of the pathway's organisation. Furthermore, associating diverse signalling responses with insulin-stimulated glucose uptake uncovered regulators of muscle insulin responsiveness, including the regulatory phosphosite S469 on Pfkfb2, a key activator of glycolysis. Finally, we confirmed the role of glycolysis in modulating insulin action in insulin resistance. Our results underscore the significance of genetics in shaping global signalling responses and their adaptability to environmental changes, emphasising the utility of studying biological diversity with phosphoproteomics to discover key regulatory mechanisms of complex traits.

Leveraging the transcriptome to further our understanding of GWAS findings: eQTLs associated with genes related to LDL and LDL subclasses, in a cohort of African Americans

Background: GWAS discoveries often pose a significant challenge in terms of understanding their underlying mechanisms. Further research, such as an integration with expression quantitative trait locus (eQTL) analyses, are required to decipher the mechanisms connecting GWAS variants to phenotypes. An eQTL analysis was conducted on genes associated with low-density lipoprotein (LDL) cholesterol and its subclasses, with the aim of pinpointing genetic variants previously implicated in GWAS studies focused on lipid-related traits. Notably, the study cohort consisted of African Americans, a population characterized by a heightened prevalence of hypercholesterolemia. Methods: A comprehensive differential expression (DE) analysis was undertaken, with a dataset of 17,948 protein-coding mRNA transcripts extracted from the whole-blood transcriptomes of 416 samples to identify mRNA transcripts associated with LDL, with further granularity delineated between small LDL and large LDL subclasses. Subsequently, eQTL analysis was conducted with a subset of 242 samples for which whole-genome sequencing data were available to identify single-nucleotide polymorphisms (SNPs) associated with the LDL-related mRNA transcripts. Lastly, plausible functional connections were established between the identified eQTLs and genetic variants reported in the GWAS catalogue. Results: DE analysis revealed 1,048, 284, and 94 mRNA transcripts that exhibited differential expression in response to LDL, small LDL, and large LDL, respectively. The eQTL analysis identified a total of 9,950 significant SNP-mRNA associations involving 6,955 SNPs including a subset 101 SNPs previously documented in GWAS of LDL and LDL-related traits. Conclusion: Through comprehensive differential expression analysis, we identified numerous mRNA transcripts responsive to LDL, small LDL, and large LDL. Subsequent eQTL analysis revealed a rich landscape of eQTL-mRNA associations, including a subset of eQTL reported in GWAS studies of LDL and related traits. The study serves as a testament to the important role of integrative genomics in unraveling the enigmatic GWAS relationships between genetic variants and the complex fabric of human traits and diseases.

Integrative Metabolomics and Proteomics Allow the Global Intracellular Characterization of Bacillus subtilis Cells and Spores

Reliable and comprehensive multi-omics analysis is essential for researchers to understand and explore complex biological systems more completely. Bacillus subtilis (B. subtilis) is a model organism for Gram-positive spore-forming bacteria, and in-depth insight into the physiology and molecular basis of spore formation and germination in this organism requires advanced multilayer molecular data sets generated from the same sample. In this study, we evaluated two monophasic methods for polar and nonpolar compound extraction (acetonitrile/methanol/water; isopropanol/water, and 60% ethanol) and two biphasic methods (chloroform/methanol/water, and methyl tert-butyl ether/methanol/water) on coefficients of variation of analytes, identified metabolite composition, and the quality of proteomics profiles. The 60% EtOH protocol proved to be the easiest in sample processing and was more amenable to automation. Collectively, we annotated 505 and 484 metabolites and identified 1665 and 1562 proteins in B. subtilis vegetative cells and spores, respectively. We also show differences between vegetative cells and spores from a multi-omics perspective and demonstrate that an integrative multi-omics analysis can be implemented from one sample using the 60% EtOH protocol. The results obtained by the 60% EtOH protocol provide comprehensive insight into differences in the metabolic and protein makeup of B. subtilis vegetative cells and spores.

Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation

Skeletal muscle, the largest human organ by weight, is relevant to several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing the relevant cell types, regulatory elements, target genes, and causal variants. Here, we used genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing 456,880 nuclei. We identified 13 cell types that collectively represented 983,155 ATAC summits. We integrated genetic variation to discover 6,866 expression quantitative trait loci (eQTL) and 100,928 chromatin accessibility QTL (caQTL) (5% FDR) across the five most abundant cell types, cataloging caQTL peaks that atlas-level snATAC maps often miss. We identified 1,973 eGenes colocalized with caQTL and used mediation analyses to construct causal directional maps for chromatin accessibility and gene expression. 3,378 genome-wide association study (GWAS) signals across 43 relevant traits colocalized with sn-e/caQTL, 52% in a cell-specific manner. 77% of GWAS signals colocalized with caQTL and not eQTL, highlighting the critical importance of population-scale chromatin profiling for GWAS functional studies. GWAS-caQTL colocalization showed distinct cell-specific regulatory paradigms. For example, a C2CD4A/B T2D GWAS signal colocalized with caQTL in muscle fibers and multiple chromatin loop models nominated VPS13C, a glucose uptake gene. Sequence of the caQTL peak overlapping caSNP rs7163757 showed allelic regulatory activity differences in a human myocyte cell line massively parallel reporter assay. These results illuminate the genetic regulatory architecture of human skeletal muscle at high-resolution epigenomic, transcriptomic, and cell state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.

Effect of Genetic and Dietary Perturbation of Glycine Metabolism on Atherosclerosis in Humans and Mice

Objective

Epidemiological and genetic studies have reported inverse associations between circulating glycine levels and risk of coronary artery disease (CAD). However, these findings have not been consistently observed in all studies. We sought to evaluate the causal relationship between circulating glycine levels and atherosclerosis using large-scale genetic analyses in humans and dietary supplementation experiments in mice.

Methods

Serum glycine levels were evaluated for association with prevalent and incident CAD in the UK Biobank. A multi-ancestry genome-wide association study (GWAS) meta-analysis was carried out to identify genetic determinants for circulating glycine levels, which were then used to evaluate the causal relationship between glycine and risk of CAD by Mendelian randomization (MR). A glycine feeding study was carried out with atherosclerosis-prone apolipoprotein E deficient (ApoE-/-) mice to determine the effects of increased circulating glycine levels on amino acid metabolism, metabolic traits, and aortic lesion formation.

Results

Among 105,718 subjects from the UK Biobank, elevated serum glycine levels were associated with significantly reduced risk of prevalent CAD (Quintile 5 vs. Quintile 1 OR=0.76, 95% CI 0.67-0.87; P<0.0001) and incident CAD (Quintile 5 vs. Quintile 1 HR=0.70, 95% CI 0.65-0.77; P<0.0001) in models adjusted for age, sex, ethnicity, anti-hypertensive and lipid-lowering medications, blood pressure, kidney function, and diabetes. A meta-analysis of 13 GWAS datasets (total n=230,947) identified 61 loci for circulating glycine levels, of which 26 were novel. MR analyses provided modest evidence that genetically elevated glycine levels were causally associated with reduced systolic blood pressure and risk of type 2 diabetes, but did provide evidence for an association with risk of CAD. Furthermore, glycine-supplementation in ApoE-/- mice did not alter cardiometabolic traits, inflammatory biomarkers, or development of atherosclerotic lesions.

Conclusions

Circulating glycine levels were inversely associated with risk of prevalent and incident CAD in a large population-based cohort. While substantially expanding the genetic architecture of circulating glycine levels, MR analyses and in vivo feeding studies in humans and mice, respectively, did not provide evidence that the clinical association of this amino acid with CAD represents a causal relationship, despite being associated with two correlated risk factors.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Factors influencing physical therapy utilization after shoulder surgery: a retrospective review

Background

Postoperative physical therapy (PT) is a cornerstone to achieve optimal patient outcomes. Access to postoperative PT can be limited by insurance type, coverage, and cost. With copayments (CP) for PT as high as $75 per visit, PT can be costprohibitive for patients. The purpose of this study was to evaluate factors affecting PT utilization among patients that underwent shoulder surgery.

Methods

A retrospective analysis was performed of 80 shoulder surgery patients with postoperative PT sessions attended at a single institution from 2017 to 2019. Patients were divided based on insurance type: private insurance (PI), and Medicare with or without supplemental insurance (MI), and CP or no copayment. Demographics, CP, total, and postoperative number of PT sessions utilized was collected and analyzed.

Results

The cohort had 53 females and an average age of 62. There was no significant difference between PI and MI at baseline other than surgery performed (P = .03), older MI group (69 years vs. 56 years: P < .01), and more females in PI group (76% vs. 55%; P = .05). There was no significant difference in the number of PT sessions between groups. The PI group was more likely to have a CP (P < .01). The CP group more often had PI and significantly more total PT visits (P = .05), while the no copayment group more often had Medicare (P < .01). CP was not independently associated with a change in the number of PT visits or total PT visits.

Conclusions

The utilization of PT after shoulder surgery was found to not be influenced by insurance type or CP as determined by the number of PT sessions attended. Further investigations are necessary to better understand the relationship between CP and different insurance types and develop effective strategies to increase access to PT for postoperative shoulder patients.

Genetic risk assessment based on association and prediction studies

The genetic basis of phenotypic emergence provides valuable information for assessing individual risk. While association studies have been pivotal in identifying genetic risk factors within a population, complementing it with insights derived from predictions studies that assess individual-level risk offers a more comprehensive approach to understanding phenotypic expression. In this study, we established personalized risk assessment models using single-nucleotide polymorphism (SNP) data from 200 Korean patients, of which 100 experienced hepatitis B surface antigen (HBsAg) seroclearance and 100 patients demonstrated high levels of HBsAg. The risk assessment models determined the predictive power of the following: (1) genome-wide association study (GWAS)-identified candidate biomarkers considered significant in a reference study and (2) machine learning (ML)-identified candidate biomarkers with the highest feature importance scores obtained by using random forest (RF). While utilizing all features yielded 64% model accuracy, using relevant biomarkers achieved higher model accuracies: 82% for 52 GWAS-identified candidate biomarkers, 71% for three GWAS-identified biomarkers, and 80% for 150 ML-identified candidate biomarkers. Findings highlight that the joint contributions of relevant biomarkers significantly influence phenotypic emergence. On the other hand, combining ML-identified candidate biomarkers into the pool of GWAS-identified candidate biomarkers resulted in the improved predictive accuracy of 90%, demonstrating the capability of ML as an auxiliary analysis to GWAS. Furthermore, some of the ML-identified candidate biomarkers were found to be linked with hepatocellular carcinoma (HCC), reinforcing previous claims that HCC can still occur despite the absence of HBsAg.

X-ray cross-complementing family: the bridge linking DNA damage repair and cancer

Genomic instability is a common hallmark of human tumours. As a carrier of genetic information, DNA is constantly threatened by various damaging factors that, if not repaired in time, can affect the transmission of genetic information and lead to cellular carcinogenesis. In response to these threats, cells have evolved a range of DNA damage response mechanisms, including DNA damage repair, to maintain genomic stability. The X-ray repair cross-complementary gene family (XRCC) comprises an important class of DNA damage repair genes that encode proteins that play important roles in DNA single-strand breakage and DNA base damage repair. The dysfunction of the XRCC gene family is associated with the development of various tumours. In the context of tumours, mutations in XRCC and its aberrant expression, result in abnormal DNA damage repair, thus contributing to the malignant progression of tumour cells. In this review, we summarise the significant roles played by XRCC in diverse tumour types. In addition, we discuss the correlation between the XRCC family members and tumour therapeutic sensitivity.

Systems genetics analysis of human body fat distribution genes identifies Wnt signaling and mitochondrial activity in adipocytes

BACKGROUND

Excess fat in the abdomen is a sexually dimorphic risk factor for cardio-metabolic disease. The relative storage between abdominal and lower-body subcutaneous adipose tissue depots is approximated by the waist-to-hip ratio adjusted for body mass index (WHRadjBMI). Genome-wide association studies (GWAS) identified 346 loci near 495 genes associated with WHRadjBMI. Most of these genes have unknown roles in fat distribution, but many are expressed and putatively act in adipose tissue. We aimed to identify novel sex- and depot-specific drivers of WHRadjBMI using a systems genetics approach.

METHODS

We used two independent cohorts of adipose tissue gene expression with 362 - 444 males and 147 - 219 females, primarily of European ancestry. We constructed sex- and depot- specific Bayesian networks to model the gene-gene interactions from 8,492 adipose tissue genes. Key driver analysis identified genes that, in silico and putatively in vitro, regulate many others, including the 495 WHRadjBMI GWAS genes. Key driver gene function was determined by perturbing their expression in human subcutaneous pre-adipocytes using lenti-virus or siRNA.

RESULTS

51 - 119 key drivers in each network were replicated in both cohorts. We used single-cell expression data to select replicated key drivers expressed in adipocyte precursors and mature adipocytes, prioritized genes which have not been previously studied in adipose tissue, and used public human and mouse data to nominate 53 novel key driver genes (10 - 21 from each network) that may regulate fat distribution by altering adipocyte function. In other cell types, 23 of these genes are found in crucial adipocyte pathways: Wnt signaling or mitochondrial function. We selected seven genes whose expression is highly correlated with WHRadjBMI to further study their effects on adipogenesis/Wnt signaling (ANAPC2, PSME3, RSPO1, TYRO3) or mitochondrial function (C1QTNF3, MIGA1, PSME3, UBR1).Adipogenesis was inhibited in cells overexpressing ANAPC2 and RSPO1 compared to controls. RSPO1 results are consistent with a positive correlation between gene expression in the subcutaneous depot and WHRadjBMI, therefore lower relative storage in the subcutaneous depot. RSPO1 inhibited adipogenesis by increasing β-catenin activation and Wnt-related transcription, thus repressing PPARG and CEBPA. PSME3 overexpression led to more adipogenesis than controls. In differentiated adipocytes, MIGA1 and UBR1 downregulation led to mitochondrial dysfunction, with lower oxygen consumption than controls; MIGA1 knockdown also lowered UCP1 expression.

SUMMARY

ANAPC2, MIGA1, PSME3, RSPO1, and UBR1 affect adipocyte function and may drive body fat distribution.

The Combination of Quantitative Proteomics and Systems Genetics Analysis Reveals that PTN Is Associated with Sleep-Loss-Induced Cognitive Impairment

Sleep loss is associated with cognitive dysfunction. However, the detailed mechanisms remain unclear. In this study, we established a para-chlorophenylalanine (PCPA)-induced insomniac mouse model with impaired cognitive function. Mass-spectrometry-based proteomics showed that the expression of 164 proteins was significantly altered in the hippocampus of the PCPA mice. To identify critical regulators among the potential markers, a transcriptome-wide association screening was performed in the BXD mice panel. Among the candidates, the expression of pleiotrophin (Ptn) was significantly associated with cognitive functions, indicating that Ptn-mediates sleep-loss-induced cognitive impairment. Gene co-expression analysis further revealed the potential mechanism by which Ptn mediates insomnia-induced cognitive impairment via the MAPK signaling pathway; that is, the decreased secretion of Ptn induced by insomnia leads to reduced binding to Ptprz1 on the postsynaptic membrane with the activation of the MAPK pathway via Fos and Nr4a1, further leading to the apoptosis of neurons. In addition, Ptn is genetically trans-regulated in the mouse hippocampus and implicated in neurodegenerative diseases in human genome-wide association studies. Our study provides a novel biomarker for insomnia-induced cognitive impairment and a new strategy for seeking neurological biomarkers by the integration of proteomics and systems genetics.

Classic genetic and hormonal switches during fetal sex development and beyond

Critical genetic and hormonal switches characterize fetal sex development in humans. They are decisive for gonadal sex determination and subsequent differentiation of the genital and somatic sex phenotype. Only at the first glace these switches seem to behave like the dual 0 and 1 system in computer sciences and lead invariably to either typically male or female phenotypes. More recent data indicate that this model is insufficient. In addition, in case of distinct mutations, many of these switches may act variably, causing a functional continuum of alterations of gene functions and -dosages, enzymatic activities, sex hormone levels, and sex hormone sensitivity, giving rise to a broad clinical spectrum of biological differences of sex development (DSD) and potentially diversity of genital and somatic sex phenotypes. The gonadal anlage is initially a bipotential organ that can develop either into a testis or an ovary. Sex-determining region Y (SRY) is the most important upstream switch of gonadal sex determination inducing SOX9 further downstream, leading to testicular Sertoli cell differentiation and the repression of ovarian pathways. If SRY is absent (virtually "switched off"), e. g., in 46,XX females, RSPO1, WNT4, FOXL2, and other factors repress the male pathway and promote ovarian development. Testosterone and its more potent derivative, dihydrotestosterone (DHT) as well as AMH, are the most important upstream hormonal switches in phenotypic sex differentiation. Masculinization of the genitalia, i. e., external genital midline fusion forming the scrotum, growth of the genital tubercle, and Wolffian duct development, occurs in response to testosterone synthesized by steroidogenic cells in the testis. Müllerian ducts will not develop into a uterus and fallopian tubes in males due to Anti-Müllerian-Hormone (AMH) produced by the Sertoli cells. The functionality of these two hormone-dependent switches is ensured by their corresponding receptors, the intracellular androgen receptor (AR) and the transmembrane AMH type II receptor. The absence of high testosterone and high AMH is crucial for anatomically female genital development during fetal life. Recent technological advances, including single-cell and spatial transcriptomics, will likely shed more light on the nature of these molecular switches.

Allelic variation in transcription factor PtoWRKY68 contributes to drought tolerance in Populus

Drought stress limits woody species productivity and influences tree distribution. However, dissecting the molecular mechanisms that underpin drought responses in forest trees can be challenging due to trait complexity. Here, using a panel of 300 Chinese white poplar (Populus tomentosa) accessions collected from different geographical climatic regions in China, we performed a genome-wide association study (GWAS) on seven drought-related traits and identified PtoWRKY68 as a candidate gene involved in the response to drought stress. A 12-bp insertion and/or deletion and three nonsynonymous variants in the PtoWRKY68 coding sequence categorized natural populations of P. tomentosa into two haplotype groups, PtoWRKY68hap1 and PtoWRKY68hap2. The allelic variation in these two PtoWRKY68 haplotypes conferred differential transcriptional regulatory activities and binding to the promoters of downstream abscisic acid (ABA) efflux and signaling genes. Overexpression of PtoWRKY68hap1 and PtoWRKY68hap2 in Arabidopsis (Arabidopsis thaliana) ameliorated the drought tolerance of two transgenic lines and increased ABA content by 42.7% and 14.3% compared to wild-type plants, respectively. Notably, PtoWRKY68hap1 (associated with drought tolerance) is ubiquitous in accessions in water-deficient environments, whereas the drought-sensitive allele PtoWRKY68hap2 is widely distributed in well-watered regions, consistent with the trends in local precipitation, suggesting that these alleles correspond to geographical adaptation in Populus. Moreover, quantitative trait loci analysis and an electrophoretic mobility shift assay showed that SHORT VEGETATIVE PHASE (PtoSVP.3) positively regulates the expression of PtoWRKY68 under drought stress. We propose a drought tolerance regulatory module in which PtoWRKY68 modulates ABA signaling and accumulation, providing insight into the genetic basis of drought tolerance in trees. Our findings will facilitate molecular breeding to improve the drought tolerance of forest trees.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

Quantifying common and distinct information in single-cell multimodal data with Tilted Canonical Correlation Analysis

Multimodal single-cell technologies profile multiple modalities for each cell simultaneously, enabling a more thorough characterization of cell populations. Existing dimension-reduction methods for multimodal data capture the "union of information," producing a lower-dimensional embedding that combines the information across modalities. While these tools are useful, we focus on a fundamentally different task of separating and quantifying the information among cells that is shared between the two modalities as well as unique to only one modality. Hence, we develop Tilted Canonical Correlation Analysis (Tilted-CCA), a method that decomposes a paired multimodal dataset into three lower-dimensional embeddings-one embedding captures the "intersection of information," representing the geometric relations among the cells that is common to both modalities, while the remaining two embeddings capture the "distinct information for a modality," representing the modality-specific geometric relations. We analyze single-cell multimodal datasets sequencing RNA along surface antibodies (i.e., CITE-seq) as well as RNA alongside chromatin accessibility (i.e., 10x) for blood cells and developing neurons via Tilted-CCA. These analyses show that Tilted-CCA enables meaningful visualization and quantification of the cross-modal information. Finally, Tilted-CCA's framework allows us to perform two specific downstream analyses. First, for single-cell datasets that simultaneously profile transcriptome and surface antibody markers, we show that Tilted-CCA helps design the target antibody panel to complement the transcriptome best. Second, for developmental single-cell datasets that simultaneously profile transcriptome and chromatin accessibility, we show that Tilted-CCA helps identify development-informative genes and distinguish between transient versus terminal cell types.

Intermediate Molecular Phenotypes to Identify Genetic Markers of Anthracycline-Induced Cardiotoxicity Risk

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.

Crosstalk between ferroptosis and chondrocytes in osteoarthritis: a systematic review of in vivo and in vitro studies

Purpose

Recent scientific reports have revealed a close association between ferroptosis and the occurrence and development of osteoarthritis (OA). Nevertheless, the precise mechanisms by which ferroptosis influences OA and how to hobble OA progression by inhibiting chondrocyte ferroptosis have not yet been fully elucidated. This study aims to conduct a comprehensive systematic review (SR) to address these gaps.

Methods

Following the guidelines of the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) 2020, we conducted a comprehensive search of the Embase, Ovid, ProQuest, PubMed, Scopus, the Cochrane Library, and Web of Science databases to identify relevant studies that investigate the association between ferroptosis and chondrocytes in OA. Our search included studies published from the inception of these databases until January 31st, 2023. Only studies that met the predetermined quality criteria were included in this SR.

Results

In this comprehensive SR, a total of 21 studies that met the specified criteria were considered suitable and included in the current updated synthesis. The mechanisms underlying chondrocyte ferroptosis and its association with OA progression involve various biological phenomena, including mitochondrial dysfunction, dysregulated iron metabolism, oxidative stress, and crucial signaling pathways.

Conclusion

Ferroptosis in chondrocytes has opened an entirely new chapter for the investigation of OA, and targeted regulation of it is springing up as an attractive and promising therapeutic tactic for OA.

Systematic review registration

https://inplasy.com/inplasy-2023-3-0044/, identifier INPLASY202330044.

Predicting mechanisms of action at genetic loci associated with discordant effects on type 2 diabetes and abdominal fat accumulation

Obesity is a major risk factor for cardiovascular disease, stroke, and type 2 diabetes (T2D). Excessive accumulation of fat in the abdomen further increases T2D risk. Abdominal obesity is measured by calculating the ratio of waist-to-hip circumference adjusted for the body-mass index (WHRadjBMI), a trait with a significant genetic inheritance. Genetic loci associated with WHRadjBMI identified in genome-wide association studies are predicted to act through adipose tissues, but many of the exact molecular mechanisms underlying fat distribution and its consequences for T2D risk are poorly understood. Further, mechanisms that uncouple the genetic inheritance of abdominal obesity from T2D risk have not yet been described. Here we utilize multi-omic data to predict mechanisms of action at loci associated with discordant effects on abdominal obesity and T2D risk. We find six genetic signals in five loci associated with protection from T2D but also with increased abdominal obesity. We predict the tissues of action at these discordant loci and the likely effector Genes (eGenes) at three discordant loci, from which we predict significant involvement of adipose biology. We then evaluate the relationship between adipose gene expression of eGenes with adipogenesis, obesity, and diabetic physiological phenotypes. By integrating these analyses with prior literature, we propose models that resolve the discordant associations at two of the five loci. While experimental validation is required to validate predictions, these hypotheses provide potential mechanisms underlying T2D risk stratification within abdominal obesity.

Genetic architecture of heart mitochondrial proteome influencing cardiac hypertrophy

Mitochondria play an important role in both normal heart function and disease etiology. We report analysis of common genetic variations contributing to mitochondrial and heart functions using an integrative proteomics approach in a panel of inbred mouse strains called the Hybrid Mouse Diversity Panel (HMDP). We performed a whole heart proteome study in the HMDP (72 strains, n=2-3 mice) and retrieved 848 mitochondrial proteins (quantified in ≥50 strains). High-resolution association mapping on their relative abundance levels revealed three trans-acting genetic loci on chromosomes (chr) 7, 13 and 17 that regulate distinct classes of mitochondrial proteins as well as cardiac hypertrophy. DAVID enrichment analyses of genes regulated by each of the loci revealed that the chr13 locus was highly enriched for complex-I proteins (24 proteins, P=2.2E-61), the chr17 locus for mitochondrial ribonucleoprotein complex (17 proteins, P=3.1E-25) and the chr7 locus for ubiquinone biosynthesis (3 proteins, P=6.9E-05). Follow-up high resolution regional mapping identified NDUFS4, LRPPRC and COQ7 as the candidate genes for chr13, chr17 and chr7 loci, respectively, and both experimental and statistical analyses supported their causal roles. Furthermore, a large cohort of Diversity Outbred mice was used to corroborate Lrpprc gene as a driver of mitochondrial DNA (mtDNA)-encoded gene regulation, and to show that the chr17 locus is specific to heart. Variations in all three loci were associated with heart mass in at least one of two independent heart stress models, namely, isoproterenol-induced heart failure and diet-induced obesity. These findings suggest that common variations in certain mitochondrial proteins can act in trans to influence tissue-specific mitochondrial functions and contribute to heart hypertrophy, elucidating mechanisms that may underlie genetic susceptibility to heart failure in human populations.

Unraveling the Host Genetic Background Effect on Internal Organ Weight Influenced by Obesity and Diabetes Using Collaborative Cross Mice

Type 2 diabetes mellitus (T2DM) is a severe chronic epidemic that results from the body's improper usage of the hormone insulin. Globally, 700 million people are expected to have received a diabetes diagnosis by 2045, according to the International Diabetes Federation (IDF). Cancer and macro- and microvascular illnesses are only a few immediate and long-term issues it could lead to. T2DM accelerates the effect of organ weights by triggering a hyperinflammatory response in the body's organs, inhibiting tissue repair and resolving inflammation. Understanding how genetic variation translates into different clinical presentations may highlight the mechanisms through which dietary elements may initiate or accelerate inflammatory disease processes and suggest potential disease-prevention techniques. To address the host genetic background effect on the organ weight by utilizing the newly developed mouse model, the Collaborative Cross mice (CC). The study was conducted on 207 genetically different CC mice from 8 CC lines of both sexes. The experiment started with 8-week-old mice for 12 weeks. During this period, one group maintained a standard chow diet (CHD), while the other group maintained a high-fat diet (HFD). In addition, body weight was recorded bi-weekly, and at the end of the study, a glucose tolerance test, as well as tissue collection (liver, spleen, heart), were conducted. Our study observed a strong effect of HFD on blood glucose clearance among different CC lines. The HFD decreased the blood glucose clearance displayed by the significant Area Under Curve (AUC) values in both populations. In addition, variation in body weight changes among the different CC lines in response to HFD. The female liver weight significantly increased compared to males in the overall population when exposed to HFD. Moreover, males showed higher heritability values than females on the same diet. Regardless of the dietary challenge, the liver weight in the overall male population correlated positively with the final body weight. The liver weight results revealed that three different CC lines perform well under classification models. The regression results also varied among organs. Accordingly, the differences among these lines correspond to the genetic variance, and we suspect that some genetic factors invoke different body responses to HFD. Further investigations, such as quantitative trait loci (QTL) analysis and genomic studies, could find these genetic elements. These findings would prove critical factors for developing personalized medicine, as they could indicate future body responses to numerous situations early, thus preventing the development of complex diseases.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

The EN-TEx resource of multi-tissue personal epigenomes & variant-impact models

Understanding how genetic variants impact molecular phenotypes is a key goal of functional genomics, currently hindered by reliance on a single haploid reference genome. Here, we present the EN-TEx resource of 1,635 open-access datasets from four donors (∼30 tissues × ∼15 assays). The datasets are mapped to matched, diploid genomes with long-read phasing and structural variants, instantiating a catalog of >1 million allele-specific loci. These loci exhibit coordinated activity along haplotypes and are less conserved than corresponding, non-allele-specific ones. Surprisingly, a deep-learning transformer model can predict the allele-specific activity based only on local nucleotide-sequence context, highlighting the importance of transcription-factor-binding motifs particularly sensitive to variants. Furthermore, combining EN-TEx with existing genome annotations reveals strong associations between allele-specific and GWAS loci. It also enables models for transferring known eQTLs to difficult-to-profile tissues (e.g., from skin to heart). Overall, EN-TEx provides rich data and generalizable models for more accurate personal functional genomics.

Identification of the regulatory mechanism of ACE2 in COVID-19-induced kidney damage with systems genetics approach

Studies showed that SARS-CoV-2 can directly target the kidney and induce renal damage. As the cell surface receptor for SARS-CoV-2 infection, the angiotensin-converting enzyme 2 (ACE2) plays a pivotal role for renal physiology and function. Thus, it is important to understand ACE2 through which pathway influences the pathogenesis of renal damage induced by COVID-19. In this study, we first performed an eQTL mapping for Ace2 in kidney tissues in 53 BXD mice strains. Results demonstrated that Ace2 is highly expressed and strongly controlled by a genetic locus on chromosome 16 in the kidney, with six genes (Dnase1, Vasn, Usp7, Abat, Mgrn1, and Rbfox1) dominated as the upstream modulator, as they are highly correlated with Ace2 expression. Gene co-expression analysis showed that Ace2 co-variates are significantly involved in the renin-angiotensin system (RAS) pathway which acts as a reno-protector. Importantly, we also found that Ace2 is positively correlated with Pdgf family members, particularly Pdgfc, which showed the most association among the 76 investigated growth factors. Mammalian Phenotype Ontology enrichment indicated that the cognate transcripts for both Ace2 and Pdgfc were mainly involved in regulating renal physiology and morphology. Among which, Cd44, Egfr, Met, Smad3, and Stat3 were identified as hub genes through protein-protein interaction analysis. Finally, in aligning with our systems genetics findings, we found ACE2, pdgf family members, and RAS genes decreased significantly in the CAKI-1 kidney cancer cells treated with S protein and receptor binding domain structural protein. Collectively, our data suggested that ACE2 work with RAS, PDGFC, as well as their cognate hub genes to regulate renal function, which could guide for future clinical prevention and targeted treatment for COVID-19-induced renal damage outcomes. KEY MESSAGES: • Ace2 is highly expressed and strongly controlled by a genetic locus on chromosome 16 in the kidney. • Ace2 co-variates are enriched in the RAS pathway. • Ace2 is strongly correlated with the growth factor Pdgfc. • Ace2 and Pdgfc co-expressed genes involved in the regulation of renal physiology and morphology. • SARS-CoV-2 spike glycoprotein induces down-regulation of Ace2, RAS, and Pdgfc.

Intermediate molecular phenotypes to identify genetic markers of anthracycline-induced cardiotoxicity risk

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex disease whose polygenic component is mainly unidentified. We propose that levels of intermediate molecular phenotypes in the myocardium associated with histopathological damage could explain CDA susceptibility; so that variants of genes encoding these intermediate molecular phenotypes could identify patients susceptible to this complication. A genetically heterogeneous cohort of mice generated by backcrossing (N = 165) was treated with doxorubicin and docetaxel. Cardiac histopathological damage was measured by fibrosis and cardiomyocyte size by an Ariol slide scanner. We determine intramyocardial levels of intermediate molecular phenotypes of CDA associated with histopathological damage and quantitative trait loci (ipQTLs) linked to them. These ipQTLs seem to contribute to the missing heritability of CDA because they improve the heritability explained by QTL directly linked to CDA (cda-QTLs) through genetic models. Genes encoding these molecular subphenotypes were evaluated as genetic markers of CDA in three cancer patient cohorts (N = 517) whose cardiac damage was quantified by echocardiography or Cardiac Magnetic Resonance. Many SNPs associated with CDA were found using genetic models. LASSO multivariate regression identified two risk score models, one for pediatric cancer patients and the other for women with breast cancer. Molecular intermediate phenotypes associated with heart damage can identify genetic markers of CDA risk, thereby allowing a more personalized patient management. A similar strategy could be applied to identify genetic markers of other complex trait diseases.

A Critical Perspective on the (Neuro)biological Foundations of Language and Linguistic Cognition

The biological foundations of language reflect assumptions about the way language and biology relate to one another, and with the rise of biological studies of language, we appear to have come closer to a deep understanding of linguistic cognition-the part of cognition constituted by language. This article argues that relations of neurobiological and genetic instantiation between linguistic cognition and the underlying biological substrate are ultimately irrelevant to understanding the higher-level structure and form of language. Linguistic patterns and those that make up the character of cognition constituted by language do not simply arise from the biological substrate because higher-level structures typically assume forms based on constraints that only emerge once these new levels are constructed. The goal is not to show how the mapping problem between linguistic cognition and neurobiology can be solved. Rather, the goal is to show the mapping problem ceases to exist once a different understanding of language-(neuro)biology relations is embraced. With this goal, this article first uncovers a number of logical and conceptual fallacies in strategies deployed in understanding language-(neuro)biology relations. After having shown these flaws, the article offers an alternative view of language-biology relations that shows how biological constraints shape language (nature and form), making it what it is.

Genome-Wide Association Studies across Environmental and Genetic Contexts Reveal Complex Genetic Architecture of Symbiotic Extended Phenotypes

A goal of modern biology is to develop the genotype-phenotype (G→P) map, a predictive understanding of how genomic information generates trait variation that forms the basis of both natural and managed communities. As microbiome research advances, however, it has become clear that many of these traits are symbiotic extended phenotypes, being governed by genetic variation encoded not only by the host's own genome, but also by the genomes of myriad cryptic symbionts. Building a reliable G→P map therefore requires accounting for the multitude of interacting genes and even genomes involved in symbiosis. Here, we use naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti paired with two genotypes of the host Medicago truncatula in four genome-wide association studies (GWAS) to determine the genomic architecture of a key symbiotic extended phenotype-partner quality, or the fitness benefit conferred to a host by a particular symbiont genotype, within and across environmental contexts and host genotypes. We define three novel categories of loci in rhizobium genomes that must be accounted for if we want to build a reliable G→P map of partner quality; namely, (i) loci whose identities depend on the environment, (ii) those that depend on the host genotype with which rhizobia interact, and (iii) universal loci that are likely important in all or most environments. IMPORTANCE Given the rapid rise of research on how microbiomes can be harnessed to improve host health, understanding the contribution of microbial genetic variation to host phenotypic variation is pressing, and will better enable us to predict the evolution of (and select more precisely for) symbiotic extended phenotypes that impact host health. We uncover extensive context-dependency in both the identity and functions of symbiont loci that control host growth, which makes predicting the genes and pathways important for determining symbiotic outcomes under different conditions more challenging. Despite this context-dependency, we also resolve a core set of universal loci that are likely important in all or most environments, and thus, serve as excellent targets both for genetic engineering and future coevolutionary studies of symbiosis.

Relative importance of genotype, gene expression, and DNA methylation on complex traits in perennial ryegrass

The growing demand for food and feed crops in the world because of growing population and more extreme weather events requires high-yielding and resilient crops. Many agriculturally important traits are polygenic, controlled by multiple regulatory layers, and with a strong interaction with the environment. In this study, 120 F₂ families of perennial ryegrass (Lolium perenne L.) were grown across a water gradient in a semifield facility with subsoil irrigation. Genomic (single-nucleotide polymorphism [SNP]), transcriptomic (gene expression [GE]), and DNA methylomic (MET) data were integrated with feed quality trait data collected from control and drought sections in the semifield facility, providing a treatment effect. Deep root length (DRL) below 110 cm was assessed with convolutional neural network image analysis. Bayesian prediction models were used to partition phenotypic variance into its components and evaluated the proportion of phenotypic variance in all traits captured by different regulatory layers (SNP, GE, and MET). The spatial effects and effects of SNP, GE, MET, the interaction between GE and MET (GE × MET) and GE × treatment (GE_Control and GE_Drought ) interaction were investigated. Gene expression explained a substantial part of the genetic and spatial variance for all the investigated phenotypes, whereas MET explained residual variance not accounted for by SNPs or GE. For DRL, MET also contributed to explaining spatial variance. The study provides a statistically elegant analytical paradigm that integrates genomic, transcriptomic, and MET information to understand the regulatory mechanisms of polygenic effects for complex traits.

Network reconstruction for trans acting genetic loci using multi-omics data and prior information

BACKGROUND

Molecular measurements of the genome, the transcriptome, and the epigenome, often termed multi-omics data, provide an in-depth view on biological systems and their integration is crucial for gaining insights in complex regulatory processes. These data can be used to explain disease related genetic variants by linking them to intermediate molecular traits (quantitative trait loci, QTL). Molecular networks regulating cellular processes leave footprints in QTL results as so-called trans-QTL hotspots. Reconstructing these networks is a complex endeavor and use of biological prior information can improve network inference. However, previous efforts were limited in the types of priors used or have only been applied to model systems. In this study, we reconstruct the regulatory networks underlying trans-QTL hotspots using human cohort data and data-driven prior information.

METHODS

We devised a new strategy to integrate QTL with human population scale multi-omics data. State-of-the art network inference methods including BDgraph and glasso were applied to these data. Comprehensive prior information to guide network inference was manually curated from large-scale biological databases. The inference approach was extensively benchmarked using simulated data and cross-cohort replication analyses. Best performing methods were subsequently applied to real-world human cohort data.

RESULTS

Our benchmarks showed that prior-based strategies outperform methods without prior information in simulated data and show better replication across datasets. Application of our approach to human cohort data highlighted two novel regulatory networks related to schizophrenia and lean body mass for which we generated novel functional hypotheses.

CONCLUSIONS

We demonstrate that existing biological knowledge can improve the integrative analysis of networks underlying trans associations and generate novel hypotheses about regulatory mechanisms.

Identification of Cardiovascular Disease-Related Genes Based on the Co-Expression Network Analysis of Genome-Wide Blood Transcriptome

Inference of co-expression network and identification of disease-related modules and gene sets can help us understand disease-related molecular pathophysiology. We aimed to identify a cardiovascular disease (CVD)-related transcriptomic signature, specifically, in peripheral blood tissue, based on differential expression (DE) and differential co-expression (DcoE) analyses. Publicly available blood sample datasets for coronary artery disease (CAD) and acute coronary syndrome (ACS) statuses were integrated to establish a co-expression network. A weighted gene co-expression network analysis was used to construct modules that include genes with highly correlated expression values. The DE criterion is a linear regression with module eigengenes for module-specific genes calculated from principal component analysis and disease status as the dependent and independent variables, respectively. The DcoE criterion is a paired t-test for intramodular connectivity between disease and matched control statuses. A total of 21 and 23 modules were established from CAD status- and ACS-related datasets, respectively, of which six modules per disease status (i.e., obstructive CAD and ACS) were selected based on the DE and DcoE criteria. For each module, gene–gene interactions with extremely high correlation coefficients were individually selected under the two conditions. Genes displaying a significant change in the number of edges (gene–gene interaction) were selected. A total of 6, 10, and 7 genes in each of the three modules were identified as potential CAD status-related genes, and 14 and 8 genes in each of the two modules were selected as ACS-related genes. Our study identified gene sets and genes that were dysregulated in CVD blood samples. These findings may contribute to the understanding of CVD pathophysiology.

Systematic multi-omics reveals the overactivation of T cell receptor signaling in immune system following bisphenol A exposure

Bisphenol A (BPA) is pervasive in the environment, and exposure to BPA may increase the incidence of noncommunicable diseases like autoimmune diseases and cancer. Although BPA causes immunological problems at the cellular level, no system-level research has been conducted on this. Hence, in this study, we aimed to gain a better understanding of the biological response to BPA exposure and its association with immunological disorders. For that, we explored the transcriptome and the proteomic modifications at the systems and cellular levels following BPA exposure. Our integrated multi-omics data showed the alteration of the T cell receptor (TCR) signaling pathway at both levels. The proportion of enlarged T cells increased with upregulation of CD69, a surface marker of early T cell activation, even though the number of T cells reduced after BPA exposure. Additionally, on BPA exposure, the levels of pLCK and pSRC increased in T cells, while that of pLAT decreased. Following BPA exposure, we investigated cytokine profiles and discovered that chitinase 3 Like 1 and matrix metalloproteinase 9 were enriched in T cells. These results indicated that T cells were hyperactivated by CD69 stimulation, and phosphorylation of SRC accelerated on BPA exposure. Hence, alteration in the TCR signaling pathway during development and differentiation due to BPA exposure could lead to insufficient and hasty activation of TCR signaling in T cells, which could modify cytokine profiles, leading to increased environmental susceptibility to chronic inflammation or diseases, increasing the chance of autoimmune diseases and cancer. This study enhances our understanding of the effects of environmental perturbations on immunosuppression at molecular, cellular, and systematic levels following pubertal BPA exposure, and may help develop better predictive, preventative, and therapeutic techniques.

The Future of Biomarkers in Veterinary Medicine: Emerging Approaches and Associated Challenges

Simple Summary

In this review we seek to outline the role of new technologies in biomarker discovery, particularly within the veterinary field and with an emphasis on ‘omics’, as well as to examine why many biomarkers-despite much excitement-have not yet made it to clinical practice. Further we emphasise the critical need for close collaboration between clinicians, researchers and funding bodies and the need to set clear goals for biomarker requirements and realistic application in the clinical setting, ensuring that biomarker type, method of detection and clinical utility are compatible, and adequate funding, time and sample size are available for all phases of development.

Abstract

New biomarkers promise to transform veterinary practice through rapid diagnosis of diseases, effective monitoring of animal health and improved welfare and production efficiency. However, the road from biomarker discovery to translation is not always straightforward. This review focuses on molecular biomarkers under development in the veterinary field, introduces the emerging technological approaches transforming this space and the role of ‘omics platforms in novel biomarker discovery. The vast majority of veterinary biomarkers are at preliminary stages of development and not yet ready to be deployed into clinical translation. Hence, we examine the major challenges encountered in the process of biomarker development from discovery, through validation and translation to clinical practice, including the hurdles specific to veterinary practice and to each of the ‘omics platforms–transcriptomics, proteomics, lipidomics and metabolomics. Finally, recommendations are made for the planning and execution of biomarker studies with a view to assisting the success of novel biomarkers in reaching their full potential.

Genome-wide association mapping of ethanol sensitivity in the Diversity Outbred mouse population

BACKGROUND

A strong predictor for the development of alcohol use disorder (AUD) is altered sensitivity to the intoxicating effects of alcohol. Individual differences in the initial sensitivity to alcohol are controlled in part by genetic factors. Mice offer a powerful tool to elucidate the genetic basis of behavioral and physiological traits relevant to AUD, but conventional experimental crosses have only been able to identify large chromosomal regions rather than specific genes. Genetically diverse, highly recombinant mouse populations make it possible to observe a wider range of phenotypic variation, offer greater mapping precision, and thus increase the potential for efficient gene identification.

METHODS

We have taken advantage of the Diversity Outbred (DO) mouse population to identify and precisely map quantitative trait loci (QTL) associated with ethanol sensitivity. We phenotyped 798 male J:DO mice for three measures of ethanol sensitivity: ataxia, hypothermia, and loss of the righting response. We used high-density MegaMUGA and GigaMUGA to obtain genotypes ranging from 77,808 to 143,259 SNPs. We also performed RNA sequencing in striatum to map expression QTLs and identify gene expression-trait correlations. We then applied a systems genetic strategy to identify narrow QTLs and construct the network of correlations that exists between DNA sequence, gene expression values, and ethanol-related phenotypes to prioritize our list of positional candidate genes.

RESULTS

We observed large amounts of phenotypic variation with the DO population and identified suggestive and significant QTLs associated with ethanol sensitivity on chromosomes 1, 2, and 16. The implicated regions were narrow (4.5-6.9 Mb in size) and each QTL explained ~4-5% of the variance.

CONCLUSIONS

Our results can be used to identify alleles that contribute to AUD in humans, elucidate causative biological mechanisms, or assist in the development of novel therapeutic interventions.

The impact of genomic variation on protein phosphorylation states and regulatory networks

Genomic variation impacts on cellular networks by affecting the abundance (e.g., protein levels) and the functional states (e.g., protein phosphorylation) of their components. Previous work has focused on the former, while in this context, the functional states of proteins have largely remained neglected. Here, we generated high‐quality transcriptome, proteome, and phosphoproteome data for a panel of 112 genomically well‐defined yeast strains. Genetic effects on transcripts were generally transmitted to the protein layer, but specific gene groups, such as ribosomal proteins, showed diverging effects on protein levels compared with RNA levels. Phosphorylation states proved crucial to unravel genetic effects on signaling networks. Correspondingly, genetic variants that cause phosphorylation changes were mostly different from those causing abundance changes in the respective proteins. Underscoring their relevance for cell physiology, phosphorylation traits were more strongly correlated with cell physiological traits such as chemical compound resistance or cell morphology, compared with transcript or protein abundance. This study demonstrates how molecular networks mediate the effects of genomic variants to cellular traits and highlights the particular importance of protein phosphorylation.

Powerful and robust inference of complex phenotypes' causal genes with dependent expression quantitative loci by a median-based Mendelian randomization

Isolating the causal genes from numerous genetic association signals in genome-wide association studies (GWASs) of complex phenotypes remains an open and challenging question. In the present study, we proposed a statistical approach, the effective-median-based Mendelian randomization (MR) framework, for inferring the causal genes of complex phenotypes with the GWAS summary statistics (named EMIC). The effective-median method solved the high false-positive issue in the existing MR methods due to either correlation among instrumental variables or noises in approximated linkage disequilibrium (LD). EMIC can further perform a pleiotropy fine-mapping analysis to remove possible false-positive estimates. With the usage of multiple cis-expression quantitative trait loci (eQTLs), EMIC was also more powerful than the alternative methods for the causal gene inference in the simulated datasets. Furthermore, EMIC rediscovered many known causal genes of complex phenotypes (schizophrenia, bipolar disorder, and total cholesterol) and reported many new and promising candidate causal genes. In sum, this study provided an efficient solution to discriminate the candidate causal genes from vast amounts of GWAS signals with eQTLs. EMIC has been implemented in our integrative software platform KGGSEE.

A Systematic Review of Tissue and Single Cell Transcriptome/Proteome Studies of the Brain in Multiple Sclerosis

Multiple sclerosis (MS) is an inflammatory demyelinating and degenerative disease of the central nervous system (CNS). Although inflammatory responses are efficiently treated, therapies for progression are scarce and suboptimal, and biomarkers to predict the disease course are insufficient. Cure or preventive measures for MS require knowledge of core pathological events at the site of the tissue damage. Novelties in systems biology have emerged and paved the way for a more fine-grained understanding of key pathological pathways within the CNS, but they have also raised questions still without answers. Here, we systemically review the power of tissue and single-cell/nucleus CNS omics and discuss major gaps of integration into the clinical practice. Systemic search identified 49 transcriptome and 11 proteome studies of the CNS from 1997 till October 2021. Pioneering molecular discoveries indicate that MS affects the whole brain and all resident cell types. Despite inconsistency of results, studies imply increase in transcripts/proteins of semaphorins, heat shock proteins, myelin proteins, apolipoproteins and HLAs. Different lesions are characterized by distinct astrocytic and microglial polarization, altered oligodendrogenesis, and changes in specific neuronal subtypes. In all white matter lesion types, CXCL12, SCD, CD163 are highly expressed, and STAT6- and TGFβ-signaling are increased. In the grey matter lesions, TNF-signaling seems to drive cell death, and especially CUX2-expressing neurons may be susceptible to neurodegeneration. The vast heterogeneity at both cellular and lesional levels may underlie the clinical heterogeneity of MS, and it may be more complex than the current disease phenotyping in the clinical practice. Systems biology has not solved the mystery of MS, but it has discovered multiple molecules and networks potentially contributing to the pathogenesis. However, these results are mostly descriptive; focused functional studies of the molecular changes may open up for a better interpretation. Guidelines for acceptable quality or awareness of results from low quality data, and standardized computational and biological pipelines may help to overcome limited tissue availability and the “snap shot” problem of omics. These may help in identifying core pathological events and point in directions for focus in clinical prevention.

Modeling common and rare genetic risk factors of neuropsychiatric disorders in human induced pluripotent stem cells

Recent genome-wide association studies (GWAS) and whole-exome sequencing of neuropsychiatric disorders, especially schizophrenia, have identified a plethora of common and rare disease risk variants/genes. Translating the mounting human genetic discoveries into novel disease biology and more tailored clinical treatments is tied to our ability to causally connect genetic risk variants to molecular and cellular phenotypes. When combined with the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/CRISPR-associated (Cas) nuclease-mediated genome editing system, human induced pluripotent stem cell (hiPSC)-derived neural cultures (both 2D and 3D organoids) provide a promising tractable cellular model for bridging the gap between genetic findings and disease biology. In this review, we first conceptualize the advances in understanding the disease polygenicity and convergence from the past decade of iPSC modeling of different types of genetic risk factors of neuropsychiatric disorders. We then discuss the major cell types and cellular phenotypes that are most relevant to neuropsychiatric disorders in iPSC modeling. Finally, we critically review the limitations of iPSC modeling of neuropsychiatric disorders and outline the need for implementing and developing novel methods to scale up the number of iPSC lines and disease risk variants in a systematic manner. Sufficiently scaled-up iPSC modeling and a better functional interpretation of genetic risk variants, in combination with cutting-edge CRISPR/Cas9 gene editing and single-cell multi-omics methods, will enable the field to identify the specific and convergent molecular and cellular phenotypes in precision for neuropsychiatric disorders.