Unraveling the molecular mechanisms of Ace2-mediated post-COVID-19 cognitive dysfunction through systems genetics approach.

The dysregulation of Angiotensin-converting enzyme 2 (ACE2) in central nervous system is believed associates with COVID-19 induced cognitive dysfunction. However, the detailed mechanism remains largely unknown. In this study, we performed a comprehensive system genetics analysis on hippocampal ACE2 based on BXD mice panel. Expression quantitative trait loci (eQTLs) mapping showed that Ace2 was strongly trans-regulated, and the elevation of Ace2 expression level was significantly correlated with impaired cognitive functions. Further Gene co-expression analysis showed that Ace2 may be correlated with the membrane proteins in Calcium signaling pathway. Further, qRT-PCR confirmed that SARS-CoV-2 spike S1 protein upregulated ACE2 expression together with eight membrane proteins in Calcium Signaling pathway. Moreover, such elevation can be attenuated by recombinant ACE2. Collectively, our findings revealed a potential mechanism of Ace2 in cognitive dysfunction, which could be beneficial for COVID-19-induced cognitive dysfunction prevention and potential treatment.

Identification of ethanol analgesia quantitative trait loci and candidate genes in BXD recombinant inbred mouse lines.

Alcohol consumption produces acute analgesic effects, and people experiencing pain conditions may drink alcohol to alleviate discomfort. However, tolerance to the analgesic properties of alcohol could prompt escalating consumption and dependence. Both nociception and alcohol-induced analgesia are under significant genetic control. Understanding the genetic architecture of these processes could inform better treatment options for people with pain conditions. This study aims to identify quantitative trait loci (QTL) driving variation in ethanol-induced analgesia across BXD recombinant inbred mouse lines. Male and female mice from 62 BXD strains received ethanol or saline oral gavage for five days and were tested for hot plate (HP) latency at baseline, Day 1, and Day 5. QTL mapping of HP phenotypes identified a significant provisional QTL on chromosome 17 for Day 1 HP latency in mice receiving ethanol. An additional highly suggestive QTL was present on chromosome 9 for the difference in pre- and post-ethanol thermal nociception. Candidate genes within QTL support intervals were provisionally identified using HP phenotypic correlations to transcriptomic database, expression QTL analysis, and other bioinformatics inquiries. The combined behavioral and bioinformatic analyses yielded strong ethanol analgesia candidate genes, specifically Myo6. Thus, the results of this genetic study of ethanol-induced analgesia in BXD mouse strains may contribute significantly to our understanding of the molecular basis for individual variation in the analgesic response to acute ethanol.

Novel Pathogenic Mutation of P209L in TRPC6 Gene Causes Adult Focal Segmental Glomerulosclerosis.

Focal segmental glomerulosclerosis (FSGS) is a leading kidney disease, clinically associated with proteinuria and progressive renal failure. The occurrence of this disease is partly related to gene mutations. We describe a single affected family member who presented with FSGS. We used high-throughput sequencing, sanger sequencing to identify the pathogenic mutations, and a systems genetics analysis in the BXD mice was conducted to explore the genetic regulatory mechanisms of pathogenic genes in the development of FSGS. We identified high urinary protein (++++) and creatinine levels (149 µmol/L) in a 29-year-old male diagnosed with a 5-year history of grade 2 hypertension. Histopathology of the kidney biopsy showed stromal hyperplasia at the glomerular segmental sclerosis and endothelial cell vacuolation degeneration. Whole-exome sequencing followed by Sanger sequencing revealed a heterozygous missense mutation (c.643C > T) in exon 2 of TRPC6, leading to the substitution of arginine with tryptophan at position 215 (p.Arg215Trp). Systems genetics analysis of the 53 BXD mice kidney transcriptomes identified Pygm as the upstream regulator of Trpc6. Those two genes are jointly involved in the regulation of FSGS mainly via Wnt and Hippo signaling pathways. We present a novel variant in the TRPC6 gene that causes FSGS. Moreover, our data suggested TRPC6 works with PYGM, as well as Wnt and Hippo signaling pathways to regulate renal function, which could guide future clinical prevention and targeted treatment for FSGS outcomes.

Systems genetics of influenza A virus-infected mice identifies TRIM21 as a critical regulator of pulmonary innate immune response.

Tripartite motif 21 (TRIM21) is a cytosolic Fc receptor that targets antibody-bound, internalized pathogens for destruction. Apart from this intrinsic defense role, TRIM21 is implicated in autoimmune diseases, inflammation, and autophagy. Whether TRIM21 participates in host interactions with influenza A virus (IAV), however, is unknown. By computational modeling of body weight and lung transcriptome data from the BXD parents (C57BL/6 J (B6) and DBA/2 J (D2)) and 41 BXD mouse strains challenged by IAV, we reveal that a Trim21-associated gene network modulates the early host responses to IAV infection. Trim21 transcripts were significantly upregulated in infected mice of both B6 and D2 backgrounds. Its expression was significantly higher in infected D2 than in infected B6 early after infection and significantly correlated with body weight loss. We identified significant trans-eQTL on chromosome 14 that regulates Trim21 expression. Nr1d2 and Il3ra were among the strongest candidate genes. Pathway analysis found Trim21 to be involved in inflammation and immunity related pathways, such as inflammation signaling pathways (TNF, IL-17, and NF-κB), viral detection signaling pathways (NOD-like and RIG-I-like), influenza, and other respiratory viral infections. Knockdown of TRIM21 in human lung epithelial A549 cells significantly augmented IAV-induced expression of IFNB1, IFNL1, CCL5, CXCL10, and IFN-stimulated genes including DDX58 and IFIH1, among others. Our data suggest that a TRIM21-associated gene network is involved in several aspects of inflammation and viral detection mechanisms during IAV infection. We identify and validate TRIM21 as a critical regulator of innate immune responses to IAV in human lung epithelial cells.

Systems genetics analysis reveals the common genetic basis for pain sensitivity and cognitive function.

BACKGROUND: There is growing evidence of a strong correlation between pain sensitivity and cognitive function under both physiological and pathological conditions. However, the detailed mechanisms remain largely unknown. In the current study, we sought to explore candidate genes and common molecular mechanisms underlying pain sensitivity and cognitive function with a transcriptome-wide association study using recombinant inbred mice from the BXD family. METHODS: The pain sensitivity determined by Hargreaves' paw withdrawal test and cognition-related phenotypes were systematically analyzed in 60 strains of BXD mice and correlated with hippocampus transcriptomes, followed by quantitative trait locus (QTL) mapping and systems genetics analysis. RESULTS: The pain sensitivity showed significant variability across the BXD strains and co-varies with cognitive traits. Pain sensitivity correlated hippocampual genes showed a significant involvement in cognition-related pathways, including glutamatergic synapse, and PI3K-Akt signaling pathway. Moreover, QTL mapping identified a genomic region on chromosome 4, potentially regulating the variation of pain sensitivity. Integrative analysis of expression QTL mapping, correlation analysis, and Bayesian network modeling identified Ring finger protein 20 (Rnf20) as the best candidate. Further pathway analysis indicated that Rnf20 may regulate the expression of pain sensitivity and cognitive function through the PI3K-Akt signaling pathway, particularly through interactions with genes Ppp2r2b, Ppp2r5c, Col9a3, Met, Rps6, Tnc, and Kras. CONCLUSIONS: Our study demonstrated that pain sensitivity is associated with genetic background and Rnf20-mediated PI3K-Akt signaling may involve in the regulation of pain sensitivity and cognitive functions.

TEAD1, MYO7A and NDUFC2 are novel functional genes associated with glucose metabolism in BXD recombinant inbred population.

AIM: The liver is an important metabolic organ that governs glucolipid metabolism, and its dysfunction may cause non-alcoholic fatty liver disease, type 2 diabetes mellitus, dyslipidaemia, etc. We aimed to systematic investigate the key factors related to hepatic glucose metabolism, which may be beneficial for understanding the underlying pathogenic mechanisms for obesity and diabetes mellitus. MATERIALS AND METHODS: Oral glucose tolerance test (OGTT) phenotypes and liver transcriptomes of BXD mice under chow and high-fat diet conditions were collected from GeneNetwork. QTL mapping was conducted to pinpoint genomic regions associated with glucose homeostasis. Candidate genes were further nominated using a multi-criteria approach and validated to confirm their functional relevance in vitro. RESULTS: Our results demonstrated that plasma glucose levels in OGTT were significantly affected by both diet and genetic background, with six genetic regulating loci were mapped on chromosomes 1, 4, and 7. Moreover, TEAD1, MYO7A and NDUFC2 were identified as the candidate genes. Functionally, siRNA-mediated TEAD1, MYO7A and NDUFC2 knockdown significantly decreased the glucose uptake and inhibited the transcription of genes related to insulin and glucose metabolism pathways. CONCLUSIONS: Our study contributes novel insights to the understanding of hepatic glucose metabolism, demonstrating the impact of TEAD1, MYO7A and NDUFC2 on mitochondrial function in the liver and their regulatory role in maintaining in glucose homeostasis.

When diet meets genetics.

Gene expression profiling of a diverse mouse population helps to decipher how a fat-rich diet contributes to inflammatory bowel disease.

Genetic and dietary modulators of the inflammatory response in the gastrointestinal tract of the BXD mouse genetic reference population.

Inflammatory gut disorders, including inflammatory bowel disease (IBD), can be impacted by dietary, environmental, and genetic factors. While the incidence of IBD is increasing worldwide, we still lack a complete understanding of the gene-by-environment interactions underlying inflammation and IBD. Here, we profiled the colon transcriptome of 52 BXD mouse strains fed with a chow or high-fat diet (HFD) and identified a subset of BXD strains that exhibit an IBD-like transcriptome signature on HFD, indicating that an interplay of genetics and diet can significantly affect intestinal inflammation. Using gene co-expression analyses, we identified modules that are enriched for IBD-dysregulated genes and found that these IBD-related modules share cis-regulatory elements that are responsive to the STAT2, SMAD3, and REL transcription factors. We used module quantitative trait locus analyses to identify genetic loci associated with the expression of these modules. Through a prioritization scheme involving systems genetics in the mouse and integration with external human datasets, we identified Muc4 and Epha6 as the top candidates mediating differences in HFD-driven intestinal inflammation. This work provides insights into the contribution of genetics and diet to IBD risk and identifies two candidate genes, MUC4 and EPHA6, that may mediate IBD susceptibility in humans.

HOXD8 suppresses renal cell carcinoma growth by upregulating SHMT1 expression.

Amplification of amino acids synthesis is reported to promote tumorigenesis. The serine/glycine biosynthesis pathway is a reversible conversion of serine and glycine catalyzed by cytoplasmic serine hydroxymethyltransferase (SHMT)1 and mitochondrial SHMT2; however, the role of SHTM1 in renal cell carcinoma (RCC) is still unclear. We found that low SHMT1 expression is correlated with poor survival of RCC patients. The in vitro study showed that overexpression of SHMT1 suppressed RCC proliferation and migration. In the mouse tumor model, SHMT1 significantly retarded RCC tumor growth. Furthermore, by gene network analysis, we found several SHMT1-related genes, among which homeobox D8 (HOXD8) was identified as the SHMT1 regulator. Knockdown of HOXD8 decreased SHMT1 expression, resulting in faster RCC growth, and rescued the SHMT1 overexpression-induced cell migration defects. Additionally, ChIP assay found the binding site of HOXD8 to SHMT1 promoter was at the -456~-254 bp region. Taken together, SHMT1 functions as a tumor suppressor in RCC. The transcription factor HOXD8 can promote SHMT1 expression and suppress RCC cell proliferation and migration, which provides new mechanisms of SHMT1 in RCC tumor growth and might be used as a potential therapeutic target candidate for clinical treatment.

Expression Levels of the Tnni3k Gene in the Heart Are Highly Associated with Cardiac and Glucose Metabolism-Related Phenotypes and Functional Pathways.

BACKGROUND: Troponin-I interacting kinase encoded by the TNNI3K gene is expressed in nuclei and Z-discs of cardiomyocytes. Mutations in TNNI3K were identified in patients with cardiac conduction diseases, arrhythmias, and cardiomyopathy. METHODS: We performed cardiac gene expression, whole genome sequencing (WGS), and cardiac function analysis in 40 strains of BXD recombinant inbred mice derived from C57BL/6J (B6) and DBA/2J (D2) strains. Expression quantitative trait loci (eQTLs) mapping and gene enrichment analysis was performed, followed by validation of candidate Tnni3k-regulatory genes. RESULTS: WGS identified compound splicing and missense T659I Tnni3k variants in the D2 parent and some BXD strains (D allele) and these strains had significantly lower Tnni3k expression than those carrying wild-type Tnni3k (B allele). Expression levels of Tnni3k significantly correlated with multiple cardiac (heart rate, wall thickness, PR duration, and T amplitude) and metabolic (glucose levels and insulin resistance) phenotypes in BXDs. A significant cis-eQTL on chromosome 3 was identified for the regulation of Tnni3k expression. Furthermore, Tnni3k-correlated genes were primarily involved in cardiac and glucose metabolism-related functions and pathways. Genes Nodal, Gnas, Nfkb1, Bmpr2, Bmp7, Smad7, Acvr1b, Acvr2b, Chrd, Tgfb3, Irs1, and Ppp1cb were differentially expressed between the B and D alleles. CONCLUSIONS: Compound splicing and T659I Tnni3k variants reduce cardiac Tnni3k expression and Tnni3k levels are associated with cardiac and glucose metabolism-related phenotypes.

Effects of genetics and sex on adolescent behaviors following neonatal ethanol exposure in BXD recombinant inbred strains.

Introduction: Fetal alcohol spectrum disorders (FASD) are the leading preventable neurodevelopmental disorders and two hallmark symptoms of FASD are abnormal behavior, and cognitive and learning deficits. The severity of alcohol's teratogenic effects on the developing brain is influenced by genetics and sex. We previously identified recombinant inbred BXD mouse strains that show differential vulnerability to ethanol-induced cell death in the developing hippocampus, a brain region important in learning and memory. The present study aimed to test the hypothesis that strains with increased vulnerability to ethanol-induced cell death in the hippocampus have concomitant deficits in multiple hippocampal-related behaviors during adolescence. Methods: The current study evaluated the effects of developmental ethanol exposure on adolescent behavior in two BXD strains that show high cell death (BXD48a, BXD100), two that show low cell death (BXD60, BXD71), and the two parental strains (C57BL/6 J (B6), DBA/2 J (D2)). On postnatal day 7, male and female neonatal pups were treated with ethanol (5.0 g/kg) or saline given in two equal doses 2 h apart. Adolescent behavior was assessed across multiple behavioral paradigms including the elevated plus maze, open field, Y-maze, and T-maze. Results: Our results demonstrate that the effects of developmental ethanol exposure on adolescent behavioral responses are highly dependent on strain. The low cell death strains, BXD60 and BXD71, showed minimal effect of ethanol exposure on all behavioral measures but did present sex differences. The parental -B6 and D2-strains and high cell death strains, BXD48a and BXD100, showed ethanol-induced effects on activity-related or anxiety-like behaviors. Interestingly, the high cell death strains were the only strains that showed a significant effect of postnatal ethanol exposure on hippocampal-dependent spatial learning and memory behaviors. Discussion: Overall, we identified effects of ethanol exposure, strain, and/or sex on multiple behavioral measures. Interestingly, the strains that showed the most effects of postnatal ethanol exposure on adolescent behavior were the BXD strains that show high ethanol-induced cell death in the neonatal hippocampus, consistent with our hypothesis. Additionally, we found evidence for interactions among strain and sex, demonstrating that these factors have a complex effect on alcohol responses and that both are important considerations.

Genome-wide screen identifies host loci that modulate Mycobacterium tuberculosis fitness in immunodivergent mice.

Genetic differences among mammalian hosts and among strains of Mycobacterium tuberculosis (Mtb) are well-established determinants of tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host-pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the highly diverse BXD family of strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb-resistant C57BL/6J (B6 or B) and Mtb-susceptible DBA/2J (D2 or D) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters of "endophenotypes," each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted quantitative trait loci (QTL) mapping of these bacterial fitness endophenotypes and identified 140 host-pathogen QTL (hpQTL). We located a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes: Rv0127 (mak), Rv0359 (rip2), Rv0955 (perM), and Rv3849 (espR). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

A new gene-scoring method for uncovering novel glaucoma-related genes using non-negative matrix factorization based on RNA-seq data.

Early diagnosis and treatment of glaucoma are challenging. The discovery of glaucoma biomarkers based on gene expression data could potentially provide new insights for early diagnosis, monitoring, and treatment options of glaucoma. Non-negative Matrix Factorization (NMF) has been widely used in numerous transcriptome data analyses in order to identify subtypes and biomarkers of different diseases; however, its application in glaucoma biomarker discovery has not been previously reported. Our study applied NMF to extract latent representations of RNA-seq data from BXD mouse strains and sorted the genes based on a novel gene scoring method. The enrichment ratio of the glaucoma-reference genes, extracted from multiple relevant resources, was compared using both the classical differentially expressed gene (DEG) analysis and NMF methods. The complete pipeline was validated using an independent RNA-seq dataset. Findings showed our NMF method significantly improved the enrichment detection of glaucoma genes. The application of NMF with the scoring method showed great promise in the identification of marker genes for glaucoma.

Gene network based analysis identifies a coexpression module involved in regulating plasma lipids with high-fat diet response.

Plasma lipids are modulated by gene variants and many environmental factors, including diet-associated weight gain. However, understanding how these factors jointly interact to influence molecular networks that regulate plasma lipid levels is limited. Here, we took advantage of the BXD recombinant inbred family of mice to query weight gain as an environmental stressor on plasma lipids. Coexpression networks were examined in both nonobese and obese livers, and a network was identified that specifically responded to the obesogenic diet. This obesity-associated module was significantly associated with plasma lipid levels and enriched with genes known to have functions related to inflammation and lipid homeostasis. We identified key drivers of the module, including Cidec, Cidea, Pparg, Cd36, and Apoa4. The Pparg emerged as a potential master regulator of the module as it can directly target 19 of the top 30 hub genes. Importantly, activation of this module is causally linked to lipid metabolism in humans, as illustrated by correlation analysis and inverse-variance weighed Mendelian randomization. Our findings provide novel insights into gene-by-environment interactions for plasma lipid metabolism that may ultimately contribute to new biomarkers, better diagnostics, and improved approaches to prevent or treat dyslipidemia in patients.

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.

Cardiac copper content and its relationship with heart physiology: Insights based on quantitative genetic and functional analyses using BXD family mice.

Background: Copper (Cu) is essential for the functioning of various enzymes involved in important cellular and physiological processes. Although critical for normal cardiac function, excessive accumulation, or deficiency of Cu in the myocardium is detrimental to the heart. Fluctuations in cardiac Cu content have been shown to cause cardiac pathologies and imbalance in systemic Cu metabolism. However, the genetic basis underlying cardiac Cu levels and their effects on heart traits remain to be understood. Representing the largest murine genetic reference population, BXD strains have been widely used to explore genotype-phenotype associations and identify quantitative trait loci (QTL) and candidate genes. Methods: Cardiac Cu concentration and heart function in BXD strains were measured, followed by QTL mapping. The candidate genes modulating Cu homeostasis in mice hearts were identified using a multi-criteria scoring/filtering approach. Results: Significant correlations were identified between cardiac Cu concentration and left ventricular (LV) internal diameter and volumes at end-diastole and end-systole, demonstrating that the BXDs with higher cardiac Cu levels have larger LV chamber. Conversely, cardiac Cu levels negatively correlated with LV posterior wall thickness, suggesting that lower Cu concentration in the heart is associated with LV hypertrophy. Genetic mapping identified six QTLs containing a total of 217 genes, which were further narrowed down to 21 genes that showed a significant association with cardiac Cu content in mice. Among those, Prex1 and Irx3 are the strongest candidates involved in cardiac Cu modulation. Conclusion: Cardiac Cu level is significantly correlated with heart chamber size and hypertrophy phenotypes in BXD mice, while being regulated by multiple genes in several QTLs. Prex1 and Irx3 may be involved in modulating Cu metabolism and its downstream effects and warrant further experimental and functional validations.

Mouse population genetics phenocopies heterogeneity of human Chd8 haploinsufficiency.

Preclinical models of neurodevelopmental disorders typically use single inbred mouse strains, which fail to capture the genetic diversity and symptom heterogeneity that is common clinically. We tested whether modeling genetic background diversity in mouse genetic reference panels would recapitulate population and individual differences in responses to a syndromic mutation in the high-confidence autism risk gene, CHD8. We measured clinically relevant phenotypes in >1,000 mice from 33 strains, including brain and body weights and cognition, activity, anxiety, and social behaviors, using 5 behavioral assays: cued fear conditioning, open field tests in dark and bright light, direct social interaction, and social dominance. Trait disruptions mimicked those seen clinically, with robust strain and sex differences. Some strains exhibited large effect-size trait disruptions, sometimes in opposite directions, and-remarkably-others expressed resilience. Therefore, systematically introducing genetic diversity into models of neurodevelopmental disorders provides a better framework for discovering individual differences in symptom etiologies.

Echocardiography phenotyping in murine genetic reference population of BXD strains reveals significant QTLs associated with cardiac function and morphology.

The genetic reference population of recombinant inbred BXD mice has been derived from crosses between C57BL/6J and DBA/2J strains. The DBA/2J parent exhibits cardiomyopathy phenotypes, whereas C57BL/6J has normal heart. BXD mice are sequenced for studying genetic interactions in cardiomyopathies. The study aimed to assess cardiomyopathy traits in BXDs and investigate the quantitative genetic architecture of those traits. Echocardiography, blood pressure, and cardiomyocyte size parameters obtained from 44 strains of BXD family (n > 5/sex) at 4-5 mo of age were associated with heart transcriptomes and expression quantitative trait loci (eQTL) mapping was performed. More than twofold variance in ejection fraction (EF%), fractional shortening (FS%), left ventricular volumes (LVVols), internal dimensions (LVIDs), mass (LVM), and posterior wall (LVPW) thickness was found among BXDs. In male BXDs, eQTL mapping identified Ndrg4 on chromosome 8 QTL to be positively correlated with LVVol and LVID and negatively associated with cardiomyocyte diameter. In female BXDs, significant QTLs were found on chromosomes 7 and 3 to be associated with LVPW and EF% and FS%, respectively, and Josd2, Dap3, and Tpm3 were predicted as strong candidate genes. Our study found variable cardiovascular traits among BXD strains and identified multiple associated QTLs, suggesting an influence of genetic background on expression of echocardiographic and cardiomyocyte diameter traits. Increased LVVol and reduced EF% and FS% represented dilated cardiomyopathy, whereas increased LV mass and wall thickness indicated hypertrophic cardiomyopathy traits. The BXD family is ideal for identifying candidate genes, causal and modifier, that influence cardiovascular phenotypes.NEW & NOTEWORTHY This study aimed to establish a cardiac phenotype-genotype correlation in murine genetic reference population of BXD RI strains by phenotyping the echocardiography, blood pressure, and cardiomyocyte diameter traits and associating each collected phenotype with genetic background. Our study identified several QTLs and candidate genes that have significant association with cardiac hypertrophy, ventricular dilation, and function including systolic hyperfunction and dysfunction.

Exploring the Regulation and Function of Rpl3l in the Development of Early-Onset Dilated Cardiomyopathy and Congestive Heart Failure Using Systems Genetics Approach.

BACKGROUND: Cardiomyopathies, diseases affecting the myocardium, are common causes of congestive heart failure (CHF) and sudden cardiac death. Recently, biallelic variants in ribosomal protein L3-like (RPL3L) have been reported to be associated with severe neonatal dilated cardiomyopathy (DCM) and CHF. This study employs a systems genetics approach to gain understanding of the regulatory mechanisms underlying the role of RPL3L in DCM. METHODS: Genetic correlation, expression quantitative trait loci (eQTL) mapping, differential expression analysis and comparative functional analysis were performed using cardiac gene expression data from the patients and murine genetic reference populations (GRPs) of BXD mice (recombinant inbred strains from a cross of C57BL/6J and DBA/2J mice). Additionally, immune infiltration analysis was performed to understand the relationship between DCM, immune cells and RPL3L expression. RESULTS: Systems genetics analysis identified high expression of Rpl3l mRNA, which ranged from 11.31 to 12.16 across murine GRPs of BXD mice, with an ~1.8-fold difference. Pathways such as "diabetic cardiomyopathy", "focal adhesion", "oxidative phosphorylation" and "DCM" were significantly associated with Rpl3l. eQTL mapping suggested Myl4 (Chr 11) and Sdha (Chr 13) as the upstream regulators of Rpl3l. The mRNA expression of Rpl3l, Myl4 and Sdha was significantly correlated with multiple echocardiography traits in BXD mice. Immune infiltration analysis revealed a significant association of RPL3L and SDHA with seven immune cells (CD4, CD8-naive T cell, CD8 T cell, macrophages, cytotoxic T cell, gamma delta T cell and exhausted T cell) that were also differentially infiltrated between heart samples obtained from DCM patients and normal individuals. CONCLUSIONS: RPL3L is highly expressed in the heart tissue of humans and mice. Expression of Rpl3l and its upstream regulators, Myl4 and Sdha, correlate with multiple cardiac function traits in murine GRPs of BXD mice, while RPL3L and SDHA correlate with immune cell infiltration in DCM patient hearts, suggesting important roles for RPL3L in DCM and CHF pathogenesis via immune inflammation, necessitating experimental validations of Myl4 and Sdha in Rpl3l regulation.

Effects of Genetics and Sex on Acute Gene Expression Changes in the Hippocampus Following Neonatal Ethanol Exposure in BXD Recombinant Inbred Mouse Strains.

Fetal alcohol spectrum disorders (FASD) are prevalent neurodevelopmental disorders. Genetics have been shown to have a role in the severity of alcohol's teratogenic effects on the developing brain. We previously identified recombinant inbred BXD mouse strains that show high (HCD) or low cell death (LCD) in the hippocampus following ethanol exposure. The present study aimed to identify gene networks that influence this susceptibility. On postnatal day 7 (3rd-trimester-equivalent), male and female neonates were treated with ethanol (5.0 g/kg) or saline, and hippocampi were collected 7hrs later. Using the Affymetrix microarray platform, ethanol-induced gene expression changes were identified in all strains with divergent expression sets found between sexes. Genes, such as Bcl2l11, Jun, and Tgfb3, showed significant strain-by-treatment interactions and were involved in many apoptosis pathways. Comparison of HCD versus LCD showed twice as many ethanol-induced genes changes in the HCD. Interestingly, these changes were regulated in the same direction suggesting (1) more perturbed effects in HCD compared to LCD and (2) limited gene expression changes that confer resistance to ethanol-induced cell death in LCD. These results demonstrate that genetic background and sex are important factors that affect differential cell death pathways after alcohol exposure during development that could have long-term consequences.