Systems genetics in diversity outbred mice inform BMD GWAS and identify determinants of bone strength

Unlocking metabolic insights with mouse genetic diversity

As part of EMBO Journal’s 2024 metabolism methods series, this commentary revisits the impact of genetics on metabolic studies, enabling dissection of novel mechanisms and phenotypes.

From Genomics to Metabolomics: Molecular Insights into Osteoporosis for Enhanced Diagnostic and Therapeutic Approaches

Osteoporosis (OP) is a prevalent skeletal disorder characterized by decreased bone mineral density (BMD) and increased fracture risk. The advancements in omics technologies-genomics, transcriptomics, proteomics, and metabolomics-have provided significant insights into the molecular mechanisms driving OP. These technologies offer critical perspectives on genetic predispositions, gene expression regulation, protein signatures, and metabolic alterations, enabling the identification of novel biomarkers for diagnosis and therapeutic targets. This review underscores the potential of these multi-omics approaches to bridge the gap between basic research and clinical applications, paving the way for precision medicine in OP management. By integrating these technologies, researchers can contribute to improved diagnostics, preventative strategies, and treatments for patients suffering from OP and related conditions.

KRAS mutation promotes the colonization of Fusobacterium nucleatum in colorectal cancer by down-regulating SERTAD4

This study explores and verifies potential molecular targets through which KRAS mutations regulate the colonization of Fusobacterium nucleatum (FN) in colorectal cancer (CRC). This study combined multiple bioinformatics methods and biological assays. Through The Cancer Genome Atlas, Gene Expression Omnibus, Human Protein Atlas, immunohistochemistry, and co-culture assays, we further confirmed the differential expression of SERTAD4 in CRC. We delved deeper into examining how expression of SERTAD4 is linked with immune cell infiltration and the enrichment of potential pathways. Lastly, through bacterial phenotypic assays, we validated the function of SERTAD4. As a molecule associated with KRAS mutations and FN infection, the expression levels of SERTAD4 were downregulated in CRC. The diagnostic efficacy of SERTAD4 for CRC is not inferior to that of CEA. Low expression of SERTAD4 is associated with poorer overall survival in CRC. Correlation analysis found that increased expression of SERTAD4 is associated with various immune cell infiltrations and immune checkpoint genes. Finally, bacterial adhesion and invasion assays verify that SERTAD4 inhibits the adhesion and invasion abilities of FN in CRC. This study demonstrates that SERTAD4 exerts a protective role in CRC by inhibiting the colonization of FN.

Trajectory-centric framework TrajAtlas reveals multi-scale differentiation heterogeneity among cells, genes, and gene modules in osteogenesis

Osteoblasts, the key cells responsible for bone formation and the maintenance of skeletal integrity, originate from a diverse array of progenitor cells. However, the mechanisms underlying osteoblast differentiation from these multiple osteoprogenitors remain poorly understood. To address this knowledge gap, we developed a comprehensive framework to investigate osteoblast differentiation at multiple scales, encompassing cells, genes, and gene modules. We constructed a reference atlas focused on differentiation, which incorporates various osteoprogenitors and provides a seven-level cellular taxonomy. To reconstruct the differentiation process, we developed a model that identifies the transcription factors and pathways involved in differentiation from different osteoprogenitors. Acknowledging that covariates such as age and tissue type can influence differentiation, we created an algorithm to detect differentially expressed genes throughout the differentiation process. Additionally, we implemented methods to identify conserved pseudotemporal gene modules across multiple samples. Overall, our framework systematically addresses the heterogeneity observed during osteoblast differentiation from diverse sources, offering novel insights into the complexities of bone formation and serving as a valuable resource for understanding osteogenesis.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease

A major fraction of loci identified by genome-wide association studies (GWASs) mediate alternative splicing, but mechanistic interpretation is hindered by the technical limitations of short-read RNA sequencing (RNA-seq), which cannot directly link splicing events to full-length protein isoforms. Long-read RNA-seq represents a powerful tool to characterize transcript isoforms, and recently, infer protein isoform existence. Here, we present an approach that integrates information from GWASs, splicing quantitative trait loci (sQTLs), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes that colocalized with BMD associations (H4PP ≥ 0.75). We generated PacBio Iso-Seq data (N = ∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were unannotated. By casting the sQTLs onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense-mediated decay and 190 that potentially resulted in the expression of unannotated protein isoforms. Finally, we functionally validated colocalizing sQTLs in TPM2, in which siRNA-mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization but exhibited no effect upon knockdown of the entire gene. Our approach should be to generalize across diverse clinical traits and to provide insights into protein isoform activities modulated by GWAS loci.

Cell type-specific network analysis in Diversity Outbred mice identifies genes potentially responsible for human bone mineral density GWAS associations

Genome-wide association studies (GWASs) have identified many sources of genetic variation associated with bone mineral density (BMD), a clinical predictor of fracture risk and osteoporosis. Aside from the identification of causal genes, other difficult challenges to informing GWAS include characterizing the roles of predicted causal genes in disease and providing additional functional context, such as the cell type predictions or biological pathways in which causal genes operate. Leveraging single-cell transcriptomics (scRNA-seq) can assist in informing BMD GWAS by linking disease-associated variants to genes and providing a cell type context for which these causal genes drive disease. Here, we use large-scale scRNA-seq data from bone marrow-derived stromal cells cultured under osteogenic conditions (BMSC-OBs) from Diversity Outbred (DO) mice to generate cell type-specific networks and contextualize BMD GWAS-implicated genes. Using trajectories inferred from the scRNA-seq data, we identify networks enriched with genes that exhibit the most dynamic changes in expression across trajectories. We discover 21 network driver genes, which are likely to be causal for human BMD GWAS associations that colocalize with expression/splicing quantitative trait loci (eQTL/sQTL). These driver genes, including Fgfrl1 and Tpx2, along with their associated networks, are predicted to be novel regulators of BMD via their roles in the differentiation of mesenchymal lineage cells. In this work, we showcase the use of single-cell transcriptomics from mouse bone-relevant cells to inform human BMD GWAS and prioritize genetic targets with potential causal roles in the development of osteoporosis.

Multi-scale cortical bone traits vary in females and males from two mouse models of genetic diversity

Understanding the genetic basis of cortical bone traits can allow for the discovery of novel genes or biological pathways regulating bone health. Mice are the most widely used mammalian model for skeletal biology and allow for the quantification of traits that cannot easily be evaluated in humans, such as osteocyte lacunar morphology. The goal of our study was to investigate the effect of genetic diversity on multi-scale cortical bone traits of 3 long bones in skeletally-mature mice. We measured bone morphology, mechanical properties, material properties, lacunar morphology, and mineral composition of mouse bones from 2 populations of genetic diversity. Additionally, we compared how intrabone relationships varied in the 2 populations. Our first population of genetic diversity included 72 females and 72 males from the 8 inbred founder strains used to create the Diversity Outbred (DO) population. These 8 strains together span almost 90% of the genetic diversity found in mice (Mus musculus). Our second population of genetic diversity included 25 genetically unique, outbred females and 25 males from the DO population. We show that multi-scale cortical bone traits vary significantly with genetic background; heritability values range from 21% to 99% indicating genetic control of bone traits across length scales. We show for the first time that lacunar shape and number are highly heritable. Comparing the 2 populations of genetic diversity, we show that each DO mouse does not resemble a single inbred founder, but instead the outbred mice display hybrid phenotypes with the elimination of extreme values. Additionally, intrabone relationships (eg, ultimate force vs. cortical area) were mainly conserved in our 2 populations. Overall, this work supports future use of these genetically diverse populations to discover novel genes contributing to cortical bone traits, especially at the lacunar length scale.

Dexamethasone-induced muscle atrophy and bone loss in six genetically diverse collaborative cross founder strains demonstrates phenotypic variability by Rg3 treatment

Background

Osteosarcopenia is a common condition characterized by the loss of both bone and muscle mass, which can lead to an increased risk of fractures and disability in older adults. The study aimed to elucidate the response of various mouse strains to treatment with Rg3, one of the leading ginsenosides, on musculoskeletal traits and immune function, and their correlation.

Methods

Six Collaborative Cross (CC) founder strains induced muscle atrophy and bone loss with dexamethasone (15 mg/kg) treatment for 1 month, and half of the mice for each strain were orally administered Rg3 (20 mg/kg). Different responses were observed depending on genetic background and Rg3 treatment.

Results

Rg3 significantly increased grip strength, running performance, and expression of muscle and bone health-related genes in a two-way analysis of variance considering the genetic backgrounds and Rg3 treatment. Significant improvements in grip strength, running performance, bone area, and muscle mass, and the increased gene expression were observed in specific strains of PWK/PhJ. For traits related to muscle, bone, and immune functions, significant correlations between traits were confirmed following Rg3 administration compared with control mice. The phenotyping analysis was compiled into a public web resource called Rg3-OsteoSarco.

Conclusion

This highlights the complex interplay between genetic determinants, pathogenesis of muscle atrophy and bone loss, and phytochemical bioactivity and the need to move away from single inbred mouse models to improve their translatability to genetically diverse humans. Rg3-OsteoSarco highlights the use of CC founder strains as a valuable tool in the field of personalized nutrition.

Neuroprotective effects of acteoside in a glaucoma mouse model by targeting Serta domain-containing protein 4

AIM

To explore the therapeutic effect and main molecular mechanisms of acteoside in a glaucoma model in DBA/2J mice.

METHODS

Proteomics was used to compare the differentially expressed proteins of C57 and DBA/2J mice. After acteoside administration in DBA/2J mice, anterior segment observation, intraocular pressure (IOP) monitoring, electrophysiology examination, and hematoxylin and eosin staining were used to analyze any potential effects. Immunohistochemistry (IHC) assays were used to verify the proteomics results. Furthermore, retinal ganglion cell 5 (RGC5) cell proliferation was assessed with cell counting kit-8 (CCK-8) assays. Serta domain-containing protein 4 (Sertad4) mRNA and protein expression levels were measured by qRT-PCR and Western blot analysis, respectively.

RESULTS

Proteomics analysis suggested that Sertad4 was the most significantly differentially expressed protein. Compared with the saline group, the acteoside treatment group showed decreased IOP, improved N1-P1 wave amplitudes, thicker retina, and larger numbers of cells in the ganglion cell layer (GCL). The IHC results showed that Sertad4 expression levels in DBA/2J mice treated with acteoside were significantly lower than in the saline group. Acteoside treatment could improve RGC5 cell survival and reduce the Sertad4 mRNA and protein expression levels after glutamate injury.

CONCLUSION

Sertad4 is differentially expressed in DBA/2J mice. Acteoside can protect RGCs from damage, possibly through the downregulation of Sertad4, and has a potential use in glaucoma treatment.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Long noncoding RNA Malat1 protects against osteoporosis and bone metastasis

MALAT1, one of the few highly conserved nuclear long noncoding RNAs (lncRNAs), is abundantly expressed in normal tissues. Previously, targeted inactivation and genetic rescue experiments identified MALAT1 as a suppressor of breast cancer lung metastasis. On the other hand, Malat1-knockout mice are viable and develop normally. On a quest to discover the fundamental roles of MALAT1 in physiological and pathological processes, we find that this lncRNA is downregulated during osteoclastogenesis in humans and mice. Remarkably, Malat1 deficiency in mice promotes osteoporosis and bone metastasis of melanoma and mammary tumor cells, which can be rescued by genetic add-back of Malat1. Mechanistically, Malat1 binds to Tead3 protein, a macrophage-osteoclast-specific Tead family member, blocking Tead3 from binding and activating Nfatc1, a master regulator of osteoclastogenesis, which results in the inhibition of Nfatc1-mediated gene transcription and osteoclast differentiation. Notably, single-cell transcriptome analysis of clinical bone samples reveals that reduced MALAT1 expression in pre-osteoclasts and osteoclasts is associated with osteoporosis and metastatic bone lesions. Altogether, these findings identify Malat1 as a lncRNA that protects against osteoporosis and bone metastasis.

Identification of Key Osteoporosis Genes Through Comparative Analysis of Men's and Women's Osteoblast Transcriptomes

Osteoporosis disproportionately affects older women, yet gender differences in human osteoblasts remain unexplored. Identifying mechanisms and biomarkers of osteoporosis will enable the development of preventative and therapeutic approaches. Transcriptome data of 187 osteoblast samples from men and women were compared. Differentially expressed genes (DEGs) were identified, and weighted gene co-expression network analysis (WGCNA) was used to discover co-expressed modules. Enrichment analysis was performed to annotate DEGs. Preservation analysis determined whether modules and pathways were similar between genders. Blood methylation, transcriptome data, mouse phenotype data, and drug treatment data were utilized to identify key osteoporosis genes. We identified 1460 DEGs enriched in immune response, neurogenesis, and GWAS osteoporosis-related genes. WGCNA uncovered 8 modules associated with immune response, development, collagen metabolism, mitochondrion, and amino acid synthesis. Preservation analysis indicated modules and pathways were generally similar between genders. Incorporating GWAS and mouse phenotype data revealed 9 key genes, including GMDS, SMOC2, SASH1, MMP2, AHCYL1, ARRDC2, IGHMBP2, ATP6V1A, and CTSK. These genes were differentially methylated in patient blood and differentiated high and low bone mineral density patients in pre- and postmenopausal women. Denosumab treatment in postmenopausal women down-regulated 6 key genes, up-regulated T cell proportions, and down-regulated fibroblast proportion. qRT-PCR was used to confirm the genes in postmenopausal women. We identified 9 key osteoporosis genes by comparing the transcriptome of osteoblasts in women and men. Our findings' clinical implications were confirmed by multi-omics data and qRT-PCR, and our study provides novel biomarkers and therapeutic targets for osteoporosis diagnosis and treatment.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Multi-view information fusion using multi-view variational autoencoder to predict proximal femoral fracture load

Background

Hip fracture occurs when an applied force exceeds the force that the proximal femur can support (the fracture load or "strength") and can have devastating consequences with poor functional outcomes. Proximal femoral strengths for specific loading conditions can be computed by subject-specific finite element analysis (FEA) using quantitative computerized tomography (QCT) images. However, the radiation and availability of QCT limit its clinical usability. Alternative low-dose and widely available measurements, such as dual energy X-ray absorptiometry (DXA) and genetic factors, would be preferable for bone strength assessment. The aim of this paper is to design a deep learning-based model to predict proximal femoral strength using multi-view information fusion.

Results

We developed new models using multi-view variational autoencoder (MVAE) for feature representation learning and a product of expert (PoE) model for multi-view information fusion. We applied the proposed models to an in-house Louisiana Osteoporosis Study (LOS) cohort with 931 male subjects, including 345 African Americans and 586 Caucasians. We performed genome-wide association studies (GWAS) to select 256 genetic variants with the lowest p-values for each proximal femoral strength and integrated whole genome sequence (WGS) features and DXA-derived imaging features to predict proximal femoral strength. The best prediction model for fall fracture load was acquired by integrating WGS features and DXA-derived imaging features. The designed models achieved the mean absolute percentage error of 18.04%, 6.84% and 7.95% for predicting proximal femoral fracture loads using linear models of fall loading, nonlinear models of fall loading, and nonlinear models of stance loading, respectively.

Conclusion

The proposed models are capable of predicting proximal femoral strength using WGS features and DXA-derived imaging features. Though this tool is not a substitute for predicting FEA using QCT images, it would make improved assessment of hip fracture risk more widely available while avoiding the increased radiation exposure from QCT.

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.

Bone Trans-omics: Integrating Omics to Unveil Mechanistic Molecular Networks Regulating Bone Biology and Disease

PURPOSE OF REVIEW

Recent advancements in "omics" technologies and bioinformatics have afforded researchers new tools to study bone biology in an unbiased and holistic way. The purpose of this review is to highlight recent studies integrating multi-omics data gathered from multiple molecular layers (i.e.; trans-omics) to reveal new molecular mechanisms that regulate bone biology and underpin skeletal diseases.

RECENT FINDINGS

Bone biologists have traditionally relied on single-omics technologies (genomics, transcriptomics, proteomics, and metabolomics) to profile measureable differences (both qualitative and quantitative) of individual molecular layers for biological discovery and to investigate mechanisms of disease. Recently, literature has grown on the implementation of integrative multi-omics to study bone biology, which combines computational and informatics support to connect multiple layers of data derived from individual "omic" platforms. This emerging discipline termed "trans-omics" has enabled bone biologists to identify and construct detailed molecular networks, unveiling new pathways and unexpected interactions that have advanced our mechanistic understanding of bone biology and disease. While the era of trans-omics is poised to revolutionize our capacity to answer more complex and diverse questions pertinent to bone pathobiology, it also brings new challenges that are inherent when trying to connect "Big Data" sets. A concerted effort between bone biologists and interdisciplinary scientists will undoubtedly be needed to extract physiologically and clinically meaningful data from bone trans-omics in order to advance its implementation in the field.

Single-Cell Integration of BMD GWAS Results Prioritize Candidate Genes Influencing Age-Related Bone Loss

The regulation of bone mineral density (BMD) is highly influenced by genetics and age. Although genome-wide association studies (GWAS) for BMD have uncovered many genes through their proximity to associated variants (variant nearest-neighbor [VNN] genes), the cell-specific mechanisms of each VNN gene remain unclear. This is primarily due to the inability to prioritize these genes by cell type and age-related expression. Using age-related transcriptomics, we found that the expression of many VNN genes was upregulated in the bone and marrow from aged mice. Candidate genes from GWAS were investigated using single-cell RNA-sequencing (scRNA-seq) datasets to enrich for cell-specific expression signatures. VNN candidate genes are highly enriched in osteo-lineage cells, osteocytes, hypertrophic chondrocytes, and Lepr+ mesenchymal stem cells. These data were used to generate a "blueprint" for Cre-loxp mouse line selection for functional validation of candidate genes and further investigation of their role in BMD maintenance throughout aging. In VNN-gene-enriched cells, Sparc, encoding the extracellular matrix (ECM) protein osteonectin, was robustly expressed. This, along with expression of numerous other ECM genes, indicates that many VNN genes likely have roles in ECM deposition by osteoblasts. Overall, we provide data supporting streamlined translation of GWAS candidate genes to potential novel therapeutic targets for the treatment of osteoporosis. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

Identification of Genetic and Genomic Influences on Progressive Ethanol Consumption in Diversity Outbred Mice

Genetic factors play a significant role in the risk for development of alcohol use disorder (AUD). Using 3-bottle choice intermittent access ethanol (IEA), we have employed the Diversity Outbred (DO) mouse panel as a model of alcohol use disorder in a genetically diverse population. Through use of gene expression network analysis techniques, in combination with expression quantitative trait loci (eQTL) mapping, we have completed an extensive analysis of the influence of genetic background on gene expression changes in the prefrontal cortex (PFC). This approach revealed that, in DO mice, genes whose expression was significantly disrupted by intermittent ethanol in the PFC also tended to be those whose expression correlated to intake. This finding is in contrast to previous studies of both mice and nonhuman primates. Importantly, these analyses identified genes involved in myelination in the PFC as significantly disrupted by IEA, correlated to ethanol intake, and having significant eQTLs. Genes that code for canonical components of the myelin sheath, such as Mbp, also emerged as key drivers of the gene expression response to intermittent ethanol drinking. Several regulators of myelination were also key drivers of gene expression, and had significant QTLs, indicating that genetic background may play an important role in regulation of brain myelination. These findings underscore the importance of disruption of normal myelination in the PFC in response to prolonged ethanol exposure, that genetic variation plays an important role in this response, and that this interaction between genetics and myelin disruption in the presence of ethanol may underlie previously observed behavioral changes under intermittent access ethanol drinking such as escalation of consumption.

Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Genetics of evolved load resistance in the skeletons of unusually large mice from Gough Island

A primary function of the skeleton is to resist the loads imparted by body weight. Genetic analyses have identified genomic regions that contribute to differences in skeletal load resistance between laboratory strains of mice, but these studies are usually restricted to 1 or 2 bones and leave open the question of how load resistance evolves in natural populations. To address these challenges, we examined the genetics of bone structure using the largest wild house mice on record, which live on Gough Island (GI). We measured structural traits connected to load resistance in the femur, tibia, scapula, humerus, radius, ulna, and mandible of GI mice, a smaller-bodied reference strain from the mainland, and 760 of their F2s. GI mice have bone geometries indicative of greater load resistance abilities but show no increase in bone mineral density compared to the mainland strain. Across traits and bones, we identified a total of 153 quantitative trait loci (QTL) that span all but one of the autosomes. The breadth of QTL detection ranges from a single bone to all 7 bones. Additive effects of QTL are modest. QTL for bone structure show limited overlap with QTL for bone length and width and QTL for body weight mapped in the same cross, suggesting a distinct genetic architecture for load resistance. Our findings provide a rare genetic portrait of the evolution of load resistance in a natural population with extreme body size.

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.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease

A major fraction of loci identified by genome-wide association studies (GWASs) lead to alterations in alternative splicing, but interpretation of how such alterations impact proteins is hindered by the technical limitations of short-read RNA-seq, which cannot directly link splicing events to full-length transcript or protein isoforms. Long-read RNA-seq represents a powerful tool to define and quantify transcript isoforms, and recently, infer protein isoform existence. Here we present a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H 4 PP ≥ 0.75). We generated deep coverage PacBio long-read RNA-seq data (N=∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were novel. By casting the colocalized sQTLs directly onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Using these data, we created one of the first proteome-scale resources defining full-length isoforms impacted by colocalized sQTLs. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense mediated decay (NMD) and 190 that potentially resulted in the expression of new protein isoforms. Finally, we identified colocalizing sQTLs in TPM2 for splice junctions between two mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization. We expect our approach to be widely generalizable across diverse clinical traits and accelerate system-scale analyses of protein isoform activities modulated by GWAS loci.

Long noncoding RNA Malat1 inhibits Tead3-Nfatc1-mediated osteoclastogenesis to suppress osteoporosis and bone metastasis

MALAT1, one of the few highly conserved nuclear long noncoding RNAs (IncRNAs), is abundantly expressed in normal tissues. Previously, targeted inactivation and genetic rescue experiments identified MALAT1 as a suppressor of breast cancer lung metastasis. On the other hand, Malat1-knockout mice are viable and develop normally. On a quest to discover new roles of MALAT1 in physiological and pathological processes, we found that this lncRNA is downregulated during osteoclastogenesis in humans and mice. Notably, Malat1 deficiency in mice promotes osteoporosis and bone metastasis, which can be rescued by genetic add-back of Malat1. Mechanistically, Malat1 binds to Tead3 protein, a macrophage-osteoclast-specific Tead family member, blocking Tead3 from binding and activating Nfatc1, a master regulator of osteoclastogenesis, which results in the inhibition of Nfatc1-mediated gene transcription and osteoclast differentiation. Altogether, these findings identify Malat1 as a lncRNA that suppresses osteoporosis and bone metastasis.

Genome-wide placental DNA methylations in fetal overgrowth and associations with leptin, adiponectin and fetal growth factors

BACKGROUND

Fetal overgrowth "programs" an elevated risk of type 2 diabetes in adulthood. Epigenetic alterations may be a mechanism in programming the vulnerability. We sought to characterize genome-wide alterations in placental gene methylations in fetal overgrowth and the associations with metabolic health biomarkers including leptin, adiponectin and fetal growth factors.

RESULTS

Comparing genome-wide placental gene DNA methylations in large-for-gestational-age (LGA, an indicator of fetal overgrowth, n = 30) versus optimal-for-gestational-age (OGA, control, n = 30) infants using the Illumina Infinium Human Methylation-EPIC BeadChip, we identified 543 differential methylation positions (DMPs; 397 hypermethylated, 146 hypomethylated) at false discovery rate < 5% and absolute methylation difference > 0.05 after adjusting for placental cell-type heterogeneity, maternal age, pre-pregnancy BMI and HbA1c levels during pregnancy. Twenty-five DMPs annotated to 20 genes (QSOX1, FCHSD2, LOC101928162, ADGRB3, GCNT1, TAP1, MYO16, NAV1, ATP8A2, LBXCOR1, EN2, INCA1, CAMTA2, SORCS2, SLC4A4, RPA3, UMAD1,USP53, OR2L13 and NR3C2) could explain 80% of the birth weight variations. Pathway analyses did not detect any statistically significant pathways after correcting for multiple tests. We validated a newly discovered differentially (hyper-)methylated gene-visual system homeobox 1 (VSX1) in an independent pyrosequencing study sample (LGA 47, OGA 47). Our data confirmed a hypermethylated gene-cadherin 13 (CDH13) reported in a previous epigenome-wide association study. Adiponectin in cord blood was correlated with its gene methylation in the placenta, while leptin and fetal growth factors (insulin, IGF-1, IGF-2) were not.

CONCLUSIONS

Fetal overgrowth may be associated with a large number of altered placental gene methylations. Placental VSX1 and CDH13 genes are hypermethylated in fetal overgrowth. Placental ADIPOQ gene methylations and fetal circulating adiponectin levels were correlated, suggesting the contribution of placenta-originated adiponectin to cord blood adiponectin.

Intramembranous bone regeneration in diversity outbred mice is heritable

There are over one million cases of failed bone repair in the U.S. annually, resulting in substantial patient morbidity and societal costs. Multiple candidate genes affecting bone traits such as bone mineral density have been identified in human subjects and animal models using genome-wide association studies (GWAS). This approach for understanding the genetic factors affecting bone repair is impractical in human subjects but could be performed in a model organism if there is sufficient variability and heritability in the bone regeneration response. Diversity Outbred (DO) mice, which have significant genetic diversity and have been used to examine multiple intact bone traits, would be an excellent possibility. Thus, we sought to evaluate the phenotypic distribution of bone regeneration, sex effects and heritability of intramembranous bone regeneration on day 7 following femoral marrow ablation in 47 12-week old DO mice (23 males, 24 females). Compared to a previous study using 4 inbred mouse strains, we found similar levels of variability in the amount of regenerated bone (coefficient of variation of 86 % v. 88 %) with approximately the same degree of heritability (0.42 v. 0.49). There was a trend toward more bone regeneration in males than females. The amount of regenerated bone was either weakly or not correlated with bone mass at intact sites, suggesting that the genetic factors responsible for bone regeneration and intact bone phenotypes are at least partially independent. In conclusion, we demonstrate that DO mice exhibit variation and heritability of intramembranous bone regeneration that will be suitable for future GWAS.

Diet X Gene Interactions Control Femoral Bone Adaptation to Low Dietary Calcium

Genetics and dietary calcium (Ca) are each critical regulators of peak bone mass but it is unclear how genetics alters the physiologic response of bone to dietary Ca restriction (RCR). Here, we conducted genetic mapping in C57BL/6J × DBA/2J (BXD) recombinant inbred mouse lines to identify environmentally sensitive loci controlling whole-bone mass (bone mineral density [BMD], bone mineral content [BMC]), distal trabecular bone, and cortical bone midshaft of the femur. Mice were fed adequate (basal) or low Ca diets from 4-12 weeks of age. Femurs were then examined by dual-energy X-ray absorptiometry (DXA) and micro-computed tomography (μCT). Body size-corrected residuals were used for statistical analysis, genetic mapping, and to estimate narrow sense heritability (h2). Genetics had a strong impact on femoral traits (eg, bone volume fraction [BV/TV] basal Ca, h2 = 0.60) as well as their RCR (eg, BV/TV, h2 = 0.32). Quantitative trait locus (QTL) mapping identified up to six loci affecting each bone trait. A subset of loci was detected in both diet groups, providing replication of environmentally robust genetic effects. Several loci control multiple bone phenotypes suggesting the existence of genetic pleiotropy. QTL controlling the bone RCR did not overlap with basal diet QTL, demonstrating genetic independence of those traits. Candidate genes underlying select multi-trait loci were prioritized by protein coding effects or gene expression differences in bone cells. These include candidate alleles in Rictor (chromosome [chr] 15) and Egfl7 (chr 2) at loci affecting bone in the basal or low Ca groups and in Msr1 (chr 8), Apc, and Camk4 (chr 18) at loci affecting RCR. By carefully controlling dietary Ca and measuring traits in age-matched mice we identified novel genetic loci determining bone mass/microarchitecture of the distal femur as well as their physiologic adaptation to inadequate dietary Ca intake. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

QTLViewer: an interactive webtool for genetic analysis in the Collaborative Cross and Diversity Outbred mouse populations

The Collaborative Cross and the Diversity Outbred mouse populations are related multiparental populations, derived from the same 8 isogenic founder strains. They carry >50 M known genetic variants, which makes them ideal tools for mapping genetic loci that regulate phenotypes, including physiological and molecular traits. Mapping quantitative trait loci requires statistical and computational training, which can present a barrier to access for some researchers. The QTLViewer software allows users to graphically explore Collaborative Cross and Diversity Outbred quantitative trait locus mapping and related analyses performed through the R/qtl2 package. Additionally, the QTLViewer website serves as a repository for published Collaborative Cross and Diversity Outbred studies, increasing the accessibility of these genetic resources to the broader scientific community.

Pharmacoepigenetics of hypertension: genome-wide methylation analysis of responsiveness to four classes of antihypertensive drugs using a double-blind crossover study design

Essential hypertension remains the leading risk factor of global disease burden, but its treatment goals are often not met. We investigated whether DNA methylation is associated with antihypertensive responses to a diuretic, a beta-blocker, a calcium channel blocker or an angiotensin receptor antagonist. In addition, since we previously showed an SNP at the transcription start site (TSS) of the catecholamine biosynthesis-related ACY3 gene to associate with blood pressure (BP) response to beta-blockers, we specifically analysed the association of methylation sites close to the ACY3 TSS with BP responses to beta-blockers. We conducted an epigenome-wide association study between leukocyte DNA methylation and BP responses to antihypertensive monotherapies in two hypertensive Finnish cohorts: the GENRES (https://clinicaltrials.gov/ct2/show/NCT03276598; amlodipine 5 mg, bisoprolol 5 mg, hydrochlorothiazide 25 mg, or losartan 50 mg daily) and the LIFE-Fin studies (https://clinicaltrials.gov/ct2/show/NCT00338260; atenolol 50 mg or losartan 50 mg daily). The monotherapy groups consisted of approximately 200 individuals each. We identified 64 methylation sites to suggestively associate (P < 1E-5) with either systolic or diastolic BP responses to a particular study drug in GENRES. These associations did not replicate in LIFE-Fin . Three methylation sites close to the ACY3 TSS were associated with systolic BP responses to bisoprolol in GENRES but not genome-wide significantly (P < 0.05). No robust associations between DNA methylation and BP responses to four different antihypertensive drugs were identified. However, the findings on the methylation sites close to the ACY3 TSS may support the role of ACY3 genetic and epigenetic variation in BP response to bisoprolol.

A generic hidden Markov model for multiparent populations

A common step in the analysis of multiparent populations (MPPs) is genotype reconstruction: identifying the founder origin of haplotypes from dense marker data. This process often makes use of a probability model for the pattern of founder alleles along chromosomes, including the relative frequency of founder alleles and the probability of exchanges among them, which depend on a model for meiotic recombination and on the mating design for the population. While the precise experimental design used to generate the population may be used to derive a precise characterization of the model for exchanges among founder alleles, this can be tedious, particularly given the great variety of experimental designs that have been proposed. We describe an approximate model that can be applied for a variety of MPPs. We have implemented the approach in the R/qtl2 software, and we illustrate its use in applications to publicly available data on Diversity Outbred and Collaborative Cross mice.

Functional Validation of Osteoporosis Genetic Findings Using Small Fish Models

The advancement of human genomics has revolutionized our understanding of the genetic architecture of many skeletal diseases, including osteoporosis. However, interpreting results from human association studies remains a challenge, since index variants often reside in non-coding regions of the genome and do not possess an obvious regulatory function. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary, such as the one offered by animal models. These models enable us to identify causal mechanisms, clarify the underlying biology, and apply interventions. Over the past several decades, small teleost fishes, mostly zebrafish and medaka, have emerged as powerful systems for modeling the genetics of human diseases. Due to their amenability to genetic intervention and the highly conserved genetic and physiological features, fish have become indispensable for skeletal genomic studies. The goal of this review is to summarize the evidence supporting the utility of Zebrafish (Danio rerio) for accelerating our understanding of human skeletal genomics and outlining the remaining gaps in knowledge. We provide an overview of zebrafish skeletal morphophysiology and gene homology, shedding light on the advantages of human skeletal genomic exploration and validation. Knowledge of the biology underlying osteoporosis through animal models will lead to the translation into new, better and more effective therapeutic approaches.

Evolutionary Genetic Signatures of Selection on Bone-Related Variation within Human and Chimpanzee Populations

Bone strength and the incidence and severity of skeletal disorders vary significantly among human populations, due in part to underlying genetic differentiation. While clinical models predict that this variation is largely deleterious, natural population variation unrelated to disease can go unnoticed, altering our perception of how natural selection has shaped bone morphologies over deep and recent time periods. Here, we conduct the first comparative population-based genetic analysis of the main bone structural protein gene, collagen type I α 1 (COL1A1), in clinical and 1000 Genomes Project datasets in humans, and in natural populations of chimpanzees. Contrary to predictions from clinical studies, we reveal abundant COL1A1 amino acid variation, predicted to have little association with disease in the natural population. We also find signatures of positive selection associated with intron haplotype structure, linkage disequilibrium, and population differentiation in regions of known gene expression regulation in humans and chimpanzees. These results recall how recent and deep evolutionary regimes can be linked, in that bone morphology differences that developed among vertebrates over 450 million years of evolution are the result of positive selection on subtle type I collagen functional variation segregating within populations over time.

Predicting the targets of IRF8 and NFATc1 during osteoclast differentiation using the machine learning method framework cTAP

BACKGROUND

Interferon regulatory factor-8 (IRF8) and nuclear factor-activated T cells c1 (NFATc1) are two transcription factors that have an important role in osteoclast differentiation. Thanks to ChIP-seq technology, scientists can now estimate potential genome-wide target genes of IRF8 and NFATc1. However, finding target genes that are consistently up-regulated or down-regulated across different studies is hard because it requires analysis of a large number of high-throughput expression studies from a comparable context.

METHOD

We have developed a machine learning based method, called, Cohort-based TF target prediction system (cTAP) to overcome this problem. This method assumes that the pathway involving the transcription factors of interest is featured with multiple "functional groups" of marker genes pertaining to the concerned biological process. It uses two notions, Gene-Present Sufficiently (GP) and Gene-Absent Insufficiently (GA), in addition to log2 fold changes of differentially expressed genes for the prediction. Target prediction is made by applying multiple machine-learning models, which learn the patterns of GP and GA from log2 fold changes and four types of Z scores from the normalized cohort's gene expression data. The learned patterns are then associated with the putative transcription factor targets to identify genes that consistently exhibit Up/Down gene regulation patterns within the cohort. We applied this method to 11 publicly available GEO data sets related to osteoclastgenesis.

RESULT

Our experiment identified a small number of Up/Down IRF8 and NFATc1 target genes as relevant to osteoclast differentiation. The machine learning models using GP and GA produced NFATc1 and IRF8 target genes different than simply using a log2 fold change alone. Our literature survey revealed that all predicted target genes have known roles in bone remodeling, specifically related to the immune system and osteoclast formation and functions, suggesting confidence and validity in our method.

CONCLUSION

cTAP was motivated by recognizing that biologists tend to use Z score values present in data sets for the analysis. However, using cTAP effectively presupposes assembling a sizable cohort of gene expression data sets within a comparable context. As public gene expression data repositories grow, the need to use cohort-based analysis method like cTAP will become increasingly important.

Intermittent fasting and caloric restriction interact with genetics to shape physiological health in mice

Dietary interventions can dramatically affect physiological health and organismal lifespan. The degree to which organismal health is improved depends upon genotype and the severity of dietary intervention, but neither the effects of these factors, nor their interaction, have been quantified in an outbred population. Moreover, it is not well understood what physiological changes occur shortly after dietary change and how these may affect the health of an adult population. In this article, we investigated the effect of 6-month exposure of either caloric restriction (CR) or intermittent fasting (IF) on a broad range of physiological traits in 960 1-year old Diversity Outbred mice. We found CR and IF affected distinct aspects of physiology and neither the magnitude nor the direction (beneficial or detrimental) of effects were concordant with the severity of the intervention. In addition to the effects of diet, genetic variation significantly affected 31 of 36 traits (heritabilities ranged from 0.04 to 0.65). We observed significant covariation between many traits that was due to both diet and genetics and quantified these effects with phenotypic and genetic correlations. We genetically mapped 16 diet-independent and 2 diet-dependent significant quantitative trait loci, both of which were associated with cardiac physiology. Collectively, these results demonstrate the degree to which diet and genetics interact to shape the physiological health of adult mice following 6 months of dietary intervention.

The Osteocyte Transcriptome: Discovering Messages Buried Within Bone

PURPOSE OF THE REVIEW

Osteocytes are cells embedded within the bone matrix, but their function and specific patterns of gene expression remain only partially defined; this is beginning to change with recent studies using transcriptomics. This unbiased approach can generate large amounts of data and is now being used to identify novel genes and signalling pathways within osteocytes both at baseline conditions and in response to stimuli. This review outlines the methods used to isolate cell populations containing osteocytes, and key recent transcriptomic studies that used osteocyte-containing preparations from bone tissue.

RECENT FINDINGS

Three common methods are used to prepare samples to examine osteocyte gene expression: digestion followed by sorting, laser capture microscopy, and the isolation of cortical bone shafts. All these methods present challenges in interpreting the data generated. Genes previously not known to be expressed by osteocytes have been identified and variations in osteocyte gene expression have been reported with age, sex, anatomical location, mechanical loading, and defects in bone strength. A substantial proportion of newly identified transcripts in osteocytes remain functionally undefined but several have been cross-referenced with functional data. Future work and improved methods (e.g. scRNAseq) likely provide useful resources for the study of osteocytes and important new information on the identity and functions of this unique cell type within the skeleton.

Utilizing preclinical models of genetic diversity to improve translation of phytochemical activities from rodents to humans and inform personalized nutrition

Mouse models are an essential tool in different areas of research, including nutrition and phytochemical research. Traditional inbred mouse models have allowed the discovery of therapeutical targets and mechanisms of action and expanded our knowledge of health and disease. However, these models lack the genetic variability typically found in human populations, which hinders the translatability of the results found in mice to humans. The development of genetically diverse mouse models, such as the collaborative cross (CC) or the diversity outbred (DO) models, has been a useful tool to overcome this obstacle in many fields, such as cancer, immunology and toxicology. However, these tools have not yet been widely adopted in the field of phytochemical research. As demonstrated in other disciplines, use of CC and DO models has the potential to provide invaluable insights for translation of phytochemicals from rodents to humans, which are desperately needed given the challenges and numerous failed clinical trials in this field. These models may prove informative for personalized use of phytochemicals in humans, including: predicting interindividual variability in phytochemical bioavailability and efficacy, identifying genetic loci or genes governing response to phytochemicals, identifying phytochemical mechanisms of action and therapeutic targets, and understanding the impact of genetic variability on individual response to phytochemicals. Such insights would prove invaluable for personalized implementation of phytochemicals in humans. This review will focus on the current work performed with genetically diverse mouse populations, and the research opportunities and advantages that these models can offer to phytochemical research.

Novel insights into the molecular mechanisms underlying risk of colorectal cancer from smoking and red/processed meat carcinogens by modeling exposure in normal colon organoids

Tobacco smoke and red/processed meats are well-known risk factors for colorectal cancer (CRC). Most research has focused on studies of normal colon biopsies in epidemiologic studies or treatment of CRC cell lines in vitro. These studies are often constrained by challenges with accuracy of self-report data or, in the case of CRC cell lines, small sample sizes and lack of relationship to normal tissue at risk. In an attempt to address some of these limitations, we performed a 24-hour treatment of a representative carcinogens cocktail in 37 independent organoid lines derived from normal colon biopsies. Machine learning algorithms were applied to bulk RNA-sequencing and revealed cellular composition changes in colon organoids. We identified 738 differentially expressed genes in response to carcinogens exposure. Network analysis identified significantly different modules of co-expression, that included genes related to MSI-H tumor biology, and genes previously implicated in CRC through genome-wide association studies. Our study helps to better define the molecular effects of representative carcinogens from smoking and red/processed meat in normal colon epithelial cells and in the etiology of the MSI-H subtype of CRC, and suggests an overlap between molecular mechanisms involved in inherited and environmental CRC risk.