Genetic dissection of the pluripotent proteome through multi-omics data integration

Algorithms and tools for data-driven omics integration to achieve multilayer biological insights: a narrative review

Systems biology is a holistic approach to biological sciences that combines experimental and computational strategies, aimed at integrating information from different scales of biological processes to unravel pathophysiological mechanisms and behaviours. In this scenario, high-throughput technologies have been playing a major role in providing huge amounts of omics data, whose integration would offer unprecedented possibilities in gaining insights on diseases and identifying potential biomarkers. In the present review, we focus on strategies that have been applied in literature to integrate genomics, transcriptomics, proteomics, and metabolomics in the year range 2018-2024. Integration approaches were divided into three main categories: statistical-based approaches, multivariate methods, and machine learning/artificial intelligence techniques. Among them, statistical approaches (mainly based on correlation) were the ones with a slightly higher prevalence, followed by multivariate approaches, and machine learning techniques. Integrating multiple biological layers has shown great potential in uncovering molecular mechanisms, identifying putative biomarkers, and aid classification, most of the time resulting in better performances when compared to single omics analyses. However, significant challenges remain. The high-throughput nature of omics platforms introduces issues such as variable data quality, missing values, collinearity, and dimensionality. These challenges further increase when combining multiple omics datasets, as the complexity and heterogeneity of the data increase with integration. We report different strategies that have been found in literature to cope with these challenges, but some open issues still remain and should be addressed to disclose the full potential of omics integration.

Natural genetic variation in wild-derived mice controls host survival and transcriptional responses during endotoxic shock

Innate immune cells sense microbial danger signals, resulting in dynamic transcriptional reprogramming and rapid inflammatory responses. If not properly regulated, such responses can be detrimental to the host, as is seen in septic shock. A better understanding of the genetic regulation of responses during endotoxemia could provide potential therapeutic insights. However, the majority of animal model studies have been performed using classic inbred laboratory strains of mice, capturing limited genetic diversity. Here, we compared classic inbred C57BL/6 (B6) mice with wild-derived and genetically divergent PWD/PhJ (PWD) mice using in vivo and in vitro models of endotoxemia. Compared with B6 mice, PWD mice were markedly resistant to bacterial lipopolysaccharide (LPS)-induced endotoxic shock. Using LPS stimulation of bone marrow derived dendritic cells (BMDC) and RNA sequencing, we demonstrate that B6 and PWD BMDCs exhibit partially overlapping yet highly divergent transcriptional responses, with B6 skewed toward stereotypical proinflammatory pathway activation, and PWD engaging regulatory or developmental pathways. To dissect genetic regulation of inflammatory responses by allelic variants, we used BMDCs from a sub-consomic strain carrying a ∼50 Mb PWD-derived portion of chromosome 11 on the B6 background. This identified a subset of cis-regulated and a large number of trans-regulated genes. Bioinformatic analyses identified candidate trans regulators encoded in the chromosome 11 locus as transcription factors Irf1, Ncor1, and Srebf1. Our results demonstrate that natural genetic variation controls host survival and transcriptional reprogramming during endotoxemia, suggesting possibilities for prediction of sepsis risk and/or personalized therapeutic interventions.

MicroRNA Gets a Mighty Award

Recent advancements in microRNAs (miRNAs) research have revealed their key roles in both normal physiological processes and pathological conditions, leading to potential applications in diagnostics and therapeutics. However, the path forward is fraught with several scientific and technical challenges. This review article briefly explores the milestones of the discovery, biogenesis, functions, and application for clinical diagnostic and therapeutic strategies of miRNAs. The potential challenges and future directions are also discussed to fully harness their capabilities.

Cross cell-type systems genetics reveals the influence of eQTL at multiple points in the developmental trajectory of mouse neural progenitor cells

Genetic variation leads to phenotypic variability in pluripotent stem cells that presents challenges for regenerative medicine. Although recent studies have investigated the impact of genetic variation on pluripotency maintenance and differentiation capacity, less is known about how genetic variants affecting the pluripotent state influence gene regulation in later stages of development. Here, we characterized expression of more than 12,000 genes in 127 donor-matched Diversity Outbred (DO) mouse embryonic stem cell (mESC) and neural progenitor cell (mNPC) lines. Quantitative trait locus (QTL) mapping identified 2,947 expression QTL (eQTL) unique to DO mNPCs and 1,113 eQTL observed in both mNPCs and mESCs with highly concordant allele effects. We mapped three eQTL hotspots on Chromosomes (Chrs) 1, 10, and 11 that were unique to mNPCs. Target genes of the Chr 1 hotspot were overrepresented for those involved in mRNA processing, DNA repair, chromatin organization, protein degradation, and cell cycle. Mediation analysis of the Chr 1 hotspot identified Rnf152 as the best candidate mediator expressed in mNPCs, while cross-cell type mediation using mESC gene expression along with partial correlation analysis strongly implicated genetic variant(s) affecting Pign expression in the mESC state as regulating the mNPC Chr 1 eQTL hotspot. Together these findings highlight that many local eQTL confer similar effects on gene expression in multiple cell states; distant eQTL in DO mNPCs are numerous and largely unique to that cell state, with many co-localizing to mNPC-specific hotspots; and mediation analysis across cell types suggests that expression of Pign early in development (mESCs) shapes the transcriptome of the more specialized mNPC state.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein coexpression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

An in vitro neurogenetics platform for precision disease modeling in the mouse

The power and scope of disease modeling can be markedly enhanced through the incorporation of broad genetic diversity. The introduction of pathogenic mutations into a single inbred mouse strain sometimes fails to mimic human disease. We describe a cross-species precision disease modeling platform that exploits mouse genetic diversity to bridge cell-based modeling with whole organism analysis. We developed a universal protocol that permitted robust and reproducible neural differentiation of genetically diverse human and mouse pluripotent stem cell lines and then carried out a proof-of-concept study of the neurodevelopmental gene DYRK1A. Results in vitro reliably predicted the effects of genetic background on Dyrk1a loss-of-function phenotypes in vivo. Transcriptomic comparison of responsive and unresponsive strains identified molecular pathways conferring sensitivity or resilience to Dyrk1a1A loss and highlighted differential messenger RNA isoform usage as an important determinant of response. This cross-species strategy provides a powerful tool in the functional analysis of candidate disease variants identified through human genetic studies.

Cell morphology QTL reveal gene by environment interactions in a genetically diverse cell population

Genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening as they offer power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.

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.

Species-wide quantitative transcriptomes and proteomes reveal distinct genetic control of gene expression variation in yeast

Gene expression varies between individuals and corresponds to a key step linking genotypes to phenotypes. However, our knowledge regarding the species-wide genetic control of protein abundance, including its dependency on transcript levels, is very limited. Here, we have determined quantitative proteomes of a large population of 942 diverse natural Saccharomyces cerevisiae yeast isolates. We found that mRNA and protein abundances are weakly correlated at the population gene level. While the protein co-expression network recapitulates major biological functions, differential expression patterns reveal proteomic signatures related to specific populations. Comprehensive genetic association analyses highlight that genetic variants associated with variation in protein (pQTL) and transcript (eQTL) levels poorly overlap (3.6%). Our results demonstrate that transcriptome and proteome are governed by distinct genetic bases, likely explained by protein turnover. It also highlights the importance of integrating these different levels of gene expression to better understand the genotype-phenotype relationship.

Highlights

At the level of individual genes, the abundance of transcripts and proteins is weakly correlated within a species ( ρ = 0.165). While the proteome is not imprinted by population structure, co-expression patterns recapitulate the cellular functional landscapeWild populations exhibit a higher abundance of respiration-related proteins compared to domesticated populationsLoci that influence protein abundance differ from those that impact transcript levels, likely because of protein turnover.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.