Joint mouse-human phenome-wide association to test gene function and disease risk

Systems genetics and bioinformatics analyses using ESR1-correlated genes identify potential candidates underlying female bone development

Estrogen receptor α (ESR1) is involved in E2 signaling and plays a major role in postmenopausal bone loss. However, the molecular network underlying ESR1 has not been explored. We used systems genetics and bioinformatics to identify important genes associated with Esr1 in postmenopausal bone loss. We identified ~2300 Esr1-coexpressed genes in female BXD bone femur, functional analysis of which revealed 'osteoblast signaling' as the most enriched pathway. PPI network led to the identification of 25 'female bone candidates'. The gene-regulatory analysis revealed RUNX2 as a key TF. ANKRD1 and RUNX2 were significantly different between osteoporosis patients and healthy controls. Sp7, Col1a1 and Pth1r correlated with multiple femur bone phenotypes in BXD mice. miR-3121-3p targeted Csf1, Ankrd1, Sp7 and Runx2. β-estradiol treatment markedly increased the expression of these candidates in mouse osteoblast. Our study revealed that Esr1-correlated genes Ankrd1, Runx2, Csf1 and Sp7 may play important roles in female bone development.

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.

A mouse model of the 3-hit effects of stress: Genotype controls the effects of life adversities in females

Helplessness is a dysfunctional coping response to stressors associated with different psychiatric conditions. The present study tested the hypothesis that early and adult adversities cumulate to produce helplessness depending on the genotype (3-hit hypothesis of psychopathology). To this aim, we evaluated whether Chronic Unpredictable Stress (CUS) differently affected coping and mesoaccumbens dopamine (DA) responses to stress challenge by adult mice of the C57BL/6J (B6) and DBA/2J (D2) inbred strains depending on early life experience (Repeated Cross Fostering, RCF). Three weeks of CUS increased the helplessness expressed in the Forced Swimming Test (FST) and the Tail Suspension Test by RCF-exposed female mice of the D2 strain. Moreover, female D2 mice with both RCF and CUS experiences showed inhibition of the stress-induced extracellular DA outflow in the Nucleus Accumbens, as measured by in vivo microdialysis, during and after FST. RCF-exposed B6 mice, instead, showed reduced helplessness and increased mesoaccumbens DA release. The present results support genotype-dependent additive effects of early experiences and adult adversities on behavioral and neural responses to stress by female mice. To our knowledge, this is the first report of a 3-hit effect in an animal model. Finally, the comparative analyses of behavioral and neural phenotypes expressed by B6 and D2 mice suggest some translationally relevant hypotheses of genetic risk factors for psychiatric disorders.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Disentangling archaic introgression and genomic signatures of selection at human immunity genes

Pathogens and infectious diseases have imposed exceptionally strong selective pressure on ancient and modern human genomes and contributed to the current variation in many genes. There is evidence that modern humans acquired immune variants through interbreeding with ancient hominins, but the impact of such variants on human traits is not fully understood. The main objectives of this research were to infer the genetic signatures of positive selection that may be involved in adaptation to infectious diseases and to investigate the function of Neanderthal alleles identified within a set of 50 Lithuanian genomes. Introgressed regions were identified using the machine learning tool ArchIE. Recent positive selection signatures were analysed using iHS. We detected high-scoring signals of positive selection at innate immunity genes (EMB, PARP8, HLAC, and CDSN) and evaluated their interactions with the structural proteins of pathogens. Interactions with human immunodeficiency virus (HIV) 1 and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) were identified. Overall, genomic regions introgressed from Neanderthals were shown to be enriched in genes related to immunity, keratinocyte differentiation, and sensory perception.

Polyglucosan body density in the aged mouse hippocampus is controlled by a novel modifier locus on chromosome 1

Aging can be associated with the accumulation of hypobranched glycogen molecules (polyglucosan bodies, PGBs), particularly in astrocytes of the hippocampus. While PGBs have a detrimental effect on cognition in diseases such as adult polyglucosan body disease and Lafora disease, the underlying mechanism and clinical relevance of age-related PGB accumulation remains unknown. Here, we have investigated the genetic basis and functional impact of age-related PGB accumulation in 32 fully sequenced BXD-type strains of mice which exhibit a 400-fold variation in PGB burden in 16-18 month old females. We mapped a major locus controlling PGB density in the hippocampus to chromosome 1 at 72-75 Mb (linkage of 4.9 -logP), which we defined as the Pgb1 locus. To identify potentially causal gene variants within Pgb1, we generated extensive hippocampal transcriptome datasets and identified two strong candidate genes for which mRNA correlates with PGB density-Smarcal1 and Usp37. In addition, both Smarcal1 and Usp37 contain non-synonymous allele variations likely to impact protein function. A phenome-wide association analysis highlighted a trans-regulatory effect of the Pgb1 locus on expression of Hp1bp3, a gene known to play a role in age-related changes in learning and memory. To investigate the potential impact of PGBs on cognition, we performed conditioned fear memory testing on strains displaying varying degrees of PGB burden, and a phenome-wide association scan of ~12,000 traits. Importantly, we did not find any evidence suggesting a negative impact of PGB burden on cognitive capacity. Taken together, we have identified a major modifier locus controlling PGB burden in the hippocampus and shed light on the genetic architecture and clinical relevance of this strikingly heterogeneous hippocampal phenotype.

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.

Tandem Mass Tag-Based Proteomic Profiling Identifies Biomarkers in Drainage Fluid for Early Detection of Anastomotic Leakage after Rectal Cancer Resection

Anastomotic leakage (AL), one of the most severe complications in rectal surgery, is often diagnosed late because of the low specificity of the clinical symptoms and limitations of current clinical investigations. Identification of patients with early AL remains challenging. Here, we explored the protein expression profiles of AL patients to provide potential biomarkers to identify AL in patients who undergo surgery for rectal cancer. We screened differentially expressed proteins (DEPs) in drainage fluid from AL and non-AL patients using a tandem mass tag method. A total of 248 DEPs, including 98 upregulated and 150 downregulated proteins, were identified between AL and non-AL groups. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses suggested that DEPs were enriched in neutrophil degranulation, bacterial infection, proteolysis, hemostasis, and complement and coagulation cascades. The results of enzyme-linked immunosorbent assay validated that the expression of the top three upregulated DEPs, AMY2A, RETN, and CELA3A, was significantly increased in the drainage fluid of AL patients, compared with that of non-AL patients (AMY2A, P = 0.001; RETN, P < 0.0001; and CELA3A, P = 0.023). Thus, our findings provide several potential biomarkers for the early diagnosis of AL after rectal cancer resection.

Adaptor protein complex 2 in the orbitofrontal cortex predicts alcohol use disorder

Alcohol use disorder (AUD) is a life-threatening disease characterized by compulsive drinking, cognitive deficits, and social impairment that continue despite negative consequences. The inability of individuals with AUD to regulate drinking may involve functional deficits in cortical areas that normally balance actions that have aspects of both reward and risk. Among these, the orbitofrontal cortex (OFC) is critically involved in goal-directed behavior and is thought to maintain a representation of reward value that guides decision making. In the present study, we analyzed post-mortem OFC brain samples collected from age- and sex-matched control subjects and those with AUD using proteomics, bioinformatics, machine learning, and reverse genetics approaches. Of the 4,500+ total unique proteins identified in the proteomics screen, there were 47 proteins that differed significantly by sex that were enriched in processes regulating extracellular matrix and axonal structure. Gene ontology enrichment analysis revealed that proteins differentially expressed in AUD cases were involved in synaptic and mitochondrial function, as well as transmembrane transporter activity. Alcohol-sensitive OFC proteins also mapped to abnormal social behaviors and social interactions. Machine learning analysis of the post-mortem OFC proteome revealed dysregulation of presynaptic (e.g., AP2A1) and mitochondrial proteins that predicted the occurrence and severity of AUD. Using a reverse genetics approach to validate a target protein, we found that prefrontal Ap2a1 expression significantly correlated with voluntary alcohol drinking in male and female genetically diverse mouse strains. Moreover, recombinant inbred strains that inherited the C57BL/6J allele at the Ap2a1 interval consumed higher amounts of alcohol than those that inherited the DBA/2J allele. Together, these findings highlight the impact of excessive alcohol consumption on the human OFC proteome and identify important cross-species cortical mechanisms and proteins that control drinking in individuals with AUD.

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.

Adaptor protein complex 2 in the orbitofrontal cortex predicts alcohol use disorder

Alcohol use disorder (AUD) is a life-threatening disease characterized by compulsive drinking, cognitive deficits, and social impairment that continue despite negative consequences. The inability of individuals with AUD to regulate drinking may involve functional deficits in cortical areas that normally balance actions that have aspects of both reward and risk. Among these, the orbitofrontal cortex (OFC) is critically involved in goal-directed behavior and is thought to maintain a representation of reward value that guides decision making. In the present study, we analyzed post-mortem OFC brain samples collected from age- and sex-matched control subjects and those with AUD using proteomics, bioinformatics, machine learning, and reverse genetics approaches. Of the 4,500+ total unique proteins identified in the proteomics screen, there were 47 proteins that differed significantly by sex that were enriched in processes regulating extracellular matrix and axonal structure. Gene ontology enrichment analysis revealed that proteins differentially expressed in AUD cases were involved in synaptic and mitochondrial function, as well as transmembrane transporter activity. Alcohol-sensitive OFC proteins also mapped to abnormal social behaviors and social interactions. Machine learning analysis of the post-mortem OFC proteome revealed dysregulation of presynaptic (e.g., AP2A1) and mitochondrial proteins that predicted the occurrence and severity of AUD. Using a reverse genetics approach to validate a target protein, we found that prefrontal Ap2a1 expression significantly correlated with voluntary alcohol drinking in male and female genetically diverse mouse strains. Moreover, recombinant inbred strains that inherited the C57BL/6J allele at the Ap2a1 interval consumed higher amounts of alcohol than those that inherited the DBA/2J allele. Together, these findings highlight the impact of excessive alcohol consumption on the human OFC proteome and identify important cross-species cortical mechanisms and proteins that control drinking in individuals with AUD.

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.

Effects of ketamine-induced H3K9 hypoacetylation during pregnancy on cardiogenesis of mouse offspring

BACKGROUND

Prenatal exposure to adverse factors can cause congenital heart defects. Ketamine, a widely used anesthetic drug, produces several adverse reactions such as tachycardia, hypertension, and laryngospasm, especially in pediatric patients. This study aimed to detect the effects of ketamine exposure during pregnancy on the cardiogenesis of mouse offspring and the potential mechanisms.

METHODS

In this study, ketamine at an addictive dose (5 mg/kg) was administered to mice during early gestation to explore the epigenetic mechanism of its causing cardiac dysplasia. The cardiac morphology of the mouse offspring was observed through hematoxylin-eosin staining and transmission electron microscopy. The heart function of one-month-old neonates was detected by echocardiography. The expression of cardiomyogenesis-related genes was detected by western blot and RT-qPCR. The acetylation level of histone H3K9 at the Mlc2 promoter and its deacetylase level and activity were detected by CHIP-qPCR, RT-qPCR, and ELISA, respectively.

RESULTS

Our data revealed that ketamine exposure during pregnancy could cause cardiac enlargement, myocardial sarcomere disorganization, and decreased cardiac contractile function in mouse offspring. Moreover, ketamine reduced the expression of Myh6, Myh7, Mlc2, Mef2c, and cTnI. The histone H3K9 acetylation level at the Mlc2 promoter was down-regulated by increasing the histone deacetylase activity and HDAC3 level upon ketamine administration.

CONCLUSIONS

Our work indicates that H3K9 acetylation is a vital player in cardiac dysplasia in offspring caused by prenatal ketamine exposure and HDAC3 is a key regulatory factor.

Genetic susceptibility to diabetic kidney disease is linked to promoter variants of XOR

The lifetime risk of kidney disease in people with diabetes is 10-30%, implicating genetic predisposition in the cause of diabetic kidney disease (DKD). Here we identify an expression quantitative trait loci (QTLs) in the cis-acting regulatory region of the xanthine dehydrogenase, or xanthine oxidoreductase (Xor), a binding site for C/EBPβ, to be associated with diabetes-induced podocyte loss in DKD in male mice. We examine mouse inbred strains that are susceptible (DBA/2J) and resistant (C57BL/6J) to DKD, as well as a panel of recombinant inbred BXD mice, to map QTLs. We also uncover promoter XOR orthologue variants in humans associated with high risk of DKD. We introduced the risk variant into the 5'-regulatory region of XOR in DKD-resistant mice, which resulted in increased Xor activity associated with podocyte depletion, albuminuria, oxidative stress and damage restricted to the glomerular endothelium, which increase further with type 1 diabetes, high-fat diet and ageing. Therefore, differential regulation of Xor contributes to phenotypic consequences with diabetes and ageing.

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.

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

Genetic differences among mammalian hosts and Mycobacterium tuberculosis ( Mtb ) strains determine diverse 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 BXD family of mouse 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 for "endophenotypes", each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted QTL mapping of these bacterial fitness endophenotypes and identified 140 h ost- p athogen quantitative trait loci ( hp QTL). We identified 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.

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.

Weighted gene co-expression network analysis identifies key hub genes and pathways in acute myeloid leukemia

Introduction: Acute myeloid leukemia (AML) is the most common type of leukemia in adults. However, there is a gap in understanding the molecular basis of the disease, partly because key genes associated with AML have not been extensively explored. In the current study, we aimed to identify genes that have strong association with AML based on a cross-species integrative approach. Methods: We used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify co-expressed gene modules significantly correlated with human AML, and further selected the genes exhibiting a significant difference in expression between AML and healthy mouse. Protein-protein interactions, transcription factors, gene function, genetic regulation, and coding sequence variants were integrated to identify key hub genes in AML. Results: The cross-species approach identified a total of 412 genes associated with both human and mouse AML. Enrichment analysis confirmed an association of these genes with hematopoietic and immune-related functions, phenotypes, processes, and pathways. Further, the integrated analysis approach identified a set of important module genes including Nfe2, Trim27, Mef2c, Ets1, Tal1, Foxo1, and Gata1 in AML. Six of these genes (except ETS1) showed significant differential expression between human AML and healthy samples in an independent microarray dataset. All of these genes are known to be involved in immune/hematopoietic functions, and in transcriptional regulation. In addition, Nfe2, Trim27, Mef2c, and Ets1 harbor coding sequence variants, whereas Nfe2 and Trim27 are cis-regulated, making them attractive candidates for validation. Furthermore, subtype-specific analysis of the hub genes in human AML indicated high expression of NFE2 across all the subtypes (M0 through M7) and enriched expression of ETS1, LEF1, GATA1, and TAL1 in M6 and M7 subtypes. A significant correlation between methylation status and expression level was observed for most of these genes in AML patients. Conclusion: Findings from the current study highlight the importance of our cross-species approach in the identification of multiple key candidate genes in AML, which can be further studied to explore their detailed role in leukemia/AML.

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.

Multi-omics analysis reveals expression complexity and functional diversity of mouse kinome

Protein kinases are a crucial component of signaling pathways involved in a wide range of cellular responses, including growth, proliferation, differentiation, and migration. Systematic investigation of protein kinases is critical to better understand phosphorylation-mediated signaling pathways and may provide insights into the development of potential therapeutic drug targets. Here we perform a systems-level analysis of the mouse kinome by analyzing multi-omics data. We used bulk and single-cell transcriptomic data from the C57BL/6J mouse strain to define tissue- and cell-type-specific expression of protein kinases, followed by investigating variations in sequence and expression between C57BL/6J and DBA/2J strains. We then profiled a deep brain phosphoproteome from C57BL/6J and DBA/2J strains as well as their reciprocal hybrids to infer the activity of the mouse kinome. Finally, we performed phenome-wide association analysis using the BXD recombinant inbred (RI) mice (a cross between C57BL/6J and DBA/2J strains) to identify any associations between variants in protein kinases and phenotypes. Collectively, our study provides a comprehensive analysis of the mouse kinome by investigating genetic sequence variation, tissue-specific expression patterns, and associations with downstream phenotypes.

UPRmt and coordinated UPRER in type 2 diabetes

The mitochondrial unfolded protein response (UPR^mt) is a molecular mechanism that maintains mitochondrial proteostasis under stress and is closely related to various metabolic diseases, such as type 2 diabetes (T2D). Similarly, the unfolded protein response of the endoplasmic reticulum (UPR^ER) is responsible for maintaining proteomic stability in the endoplasmic reticulum (ER). Since the mitochondria and endoplasmic reticulum are the primary centers of energy metabolism and protein synthesis in cells, respectively, a synergistic mechanism must exist between UPR^mt and UPR^ER to cooperatively resist stresses such as hyperglycemia in T2D. Increasing evidence suggests that the protein kinase RNA (PKR)-like endoplasmic reticulum kinase (PERK) signaling pathway is likely an important node for coordinating UPR^mt and UPR^ER. The PERK pathway is activated in both UPR^mt and UPR^ER, and its downstream molecules perform important functions. In this review, we discuss the mechanisms of UPR^mt, UPR^ER and their crosstalk in T2D.

Meiotic drive in house mice: mechanisms, consequences, and insights for human biology

Meiotic drive occurs when one allele at a heterozygous site cheats its way into a disproportionate share of functional gametes, violating Mendel's law of equal segregation. This genetic conflict typically imposes a fitness cost to individuals, often by disrupting the process of gametogenesis. The evolutionary impact of meiotic drive is substantial, and the phenomenon has been associated with infertility and reproductive isolation in a wide range of organisms. However, cases of meiotic drive in humans remain elusive, a finding that likely reflects the inherent challenges of detecting drive in our species rather than unique features of human genome biology. Here, we make the case that house mice (Mus musculus) present a powerful model system to investigate the mechanisms and consequences of meiotic drive and facilitate translational inferences about the scope and potential mechanisms of drive in humans. We first detail how different house mouse resources have been harnessed to identify cases of meiotic drive and the underlying mechanisms utilized to override Mendel's rules of inheritance. We then summarize the current state of knowledge of meiotic drive in the mouse genome. We profile known mechanisms leading to transmission bias at several established drive elements. We discuss how a detailed understanding of meiotic drive in mice can steer the search for drive elements in our own species. Lastly, we conclude with a prospective look into how new technologies and molecular tools can help resolve lingering mysteries about the prevalence and mechanisms of selfish DNA transmission in mammals.

Genetic characteristics involving the PD-1/PD-L1/L2 and CD73/A2aR axes and the immunosuppressive microenvironment in DLBCL

Background

Targeting the PD-1/PD-L1/L2 (programmed cell death protein 1/programmed cell death ligand 1/ligand 2) pathway combined with other immunosuppressive signalings, such as CD73/A2aR (A2a adenosine receptor) adenosine signaling, has emerged as a promising strategy for cancer treatment. The genetic characteristics of these immune checkpoints need to be further investigated in diffuse large B-cell lymphoma (DLBCL).

Methods

We performed whole-exome sequencing/targeted deep sequencing to investigate the genetic characteristics of PD-1/PD-L1/L2 and CD73/A2aR. The immunosuppressive effect of these two pathways on the tumor microenvironment was evaluated via RNA sequencing. Single-cell RNA sequencing was further applied to investigate the dysfunctional CD8+ T cells. In addition, multiplex immunofluorescence staining was used to quantitatively assess the expression of dysfunctional CD8⁺ T cells in DLBCL.

Results

SP140 was identified as a novel translocation partner for PD-L1, and a new inversion was detected between PD-L1 and PD-L2, both leading to the upregulation of PD-L1 expression. CD73 genetic mutations did not increase mRNA and protein expression. Patients with genetically altered CD73 tended to have a better overall survival than patients with wild-type CD73. Both PD-1/PD-L1 and CD73/A2aR signaling mediated the immunosuppressive microenvironment in DLBCL. The numbers of CD8⁺ T cells with PD-1 and A2aR expression were positively correlated with the number of dysfunctional CD8⁺ T cells (R²=0.974, p=0.013). According to the grades of dysfunctional CD8⁺ T cells we defined, grade 1 dysfunctional CD8⁺ T cells, with either PD-1⁺ or A2aR⁺, were significantly associated with poorer survival than grade 0 dysfunctional CD8⁺ T cells, with both PD-1⁻ and A2aR⁻; and patients with grade 2 dysfunctional CD8⁺ T cells showed the worst clinical outcomes.

Conclusions

This study describes the additional genetic basis of PD-L1 overexpression and characterizes certain genetic alterations of CD73/A2aR in DLBCL. The degree of T-cell dysfunction is correlated with clinical outcomes. Strategies that reverse T-cell dysfunction by inhibiting PD-1/PD-L1/L2, particularly in combination with CD73/A2aR, may show potential as effective therapeutic options for DLBCL.

Selection of Cashmere Fineness Functional Genes by Translatomics

Cashmere fineness is an important index to evaluate cashmere quality. Liaoning Cashmere Goat (LCG) has a large cashmere production and long cashmere fiber, but its fineness is not ideal. Therefore, it is important to find genes involved in cashmere fineness that can be used in future endeavors aiming to improve this phenotype. With the continuous advancement of research, the regulation of cashmere fineness has made new developments through high-throughput sequencing and genome-wide association analysis. It has been found that translatomics can identify genes associated with phenotypic traits. Through translatomic analysis, the skin tissue of LCG sample groups differing in cashmere fineness was sequenced by Ribo-seq. With these data, we identified 529 differentially expressed genes between the sample groups among the 27197 expressed genes. From these, 343 genes were upregulated in the fine LCG group in relation to the coarse LCG group, and 186 were downregulated in the same relationship. Through GO enrichment analysis and KEGG enrichment analysis of differential genes, the biological functions and pathways of differential genes can be found. In the GO enrichment analysis, 491 genes were significantly enriched, and the functional region was mainly in the extracellular region. In the KEGG enrichment analysis, the enrichment of the human papillomavirus infection pathway was seen the most. We found that the COL6A5 gene may affect cashmere fineness.

NHR-80 senses the mitochondrial UPR to rewire citrate metabolism for lipid accumulation in Caenorhabditis elegans

Mitochondria are known as the powerhouse of the cell. Dysfunction of mitochondria homeostasis induces the mitochondrial unfolded protein response (UPR^mt), altering cellular metabolism. How cells sense the UPR^mt to rewire metabolism is largely unknown. Here, we show that inactivation of either the citric/tricarboxylic acid (TCA) cycle enzymes aco-2 or idha-1, which encode aconitase and isocitrate dehydrogenase respectively, leads to citrate accumulation. In Caenorhabditis elegans, both in vitro and in vivo, citrate accumulation consequently triggers the UPR^mt and also promotes lipid accumulation. The transcription factor DVE-1 binds to the promoter of the nuclear hormone receptor nhr-80 to transactivate its expression. NHR-80 then upregulates lipogenesis and lipid accumulation, shifting excess citrate for use in lipogenesis and for storage as triacylglycerol in lipid droplets. Inactivation of DVE-1 or NHR-80 fully abolishes the citrate-induced lipid accumulation. Therefore, our work uncovers a DVE-1-NHR-80-lipogenesis axis linking the transmission of the mitochondrial stress signal to lipid metabolism.

Genetic architecture of protein expression and its regulation in the mouse brain

BACKGROUND

Natural variation in protein expression is common in all organisms and contributes to phenotypic differences among individuals. While variation in gene expression at the transcript level has been extensively investigated, the genetic mechanisms underlying variation in protein expression have lagged considerably behind. Here we investigate genetic architecture of protein expression by profiling a deep mouse brain proteome of two inbred strains, C57BL/6 J (B6) and DBA/2 J (D2), and their reciprocal F1 hybrids using two-dimensional liquid chromatography coupled with tandem mass spectrometry (LC/LC-MS/MS) technology.

RESULTS

By comparing protein expression levels in the four mouse strains, we observed 329 statistically significant differentially expressed proteins between the two parental strains and characterized the genetic basis of protein expression. We further applied a proteogenomic approach to detect variant peptides and define protein allele-specific expression (pASE), identifying 33 variant peptides with cis-effects and 17 variant peptides showing trans-effects. Comparison of regulation at transcript and protein levels show a significant divergence.

CONCLUSIONS

The results provide a comprehensive analysis of genetic architecture of protein expression and the contribution of cis- and trans-acting regulatory differences to protein expression.

BXD Recombinant Inbred Mice as a Model to Study Neurotoxicity

BXD recombinant inbred (RI) lines represent a genetic reference population derived from a cross between C57BL/6J mice (B6) and DBA/2J mice (D2), which through meiotic recombination events possesses recombinant chromosomes containing B6 or D2 haplotype segments. The quantitative trait loci (QTLs) are the locations of segregating genetic polymorphisms and are fundamental to understanding genetic diversity in human disease susceptibility and severity. QTL mapping represents the typical approach for identifying naturally occurring polymorphisms that influence complex phenotypes. In this process, genotypic values at markers of known genomic locations are associated with phenotypic values measured in a segregating population. Indeed, BXD RI strains provide a powerful tool to study neurotoxicity induced by different substances. In this review, we describe the use of BXD RI lines to understand the underlying mechanisms of neurotoxicity in response to ethanol and cocaine, as well as metals and pesticide exposures.

Pursuit of precision medicine: Systems biology approaches in Alzheimer's disease mouse models

Alzheimer's disease (AD) is a complex disease that is mediated by numerous factors and manifests in various forms. A systems biology approach to studying AD involves analyses of various body systems, biological scales, environmental elements, and clinical outcomes to understand the genotype to phenotype relationship that potentially drives AD development. Currently, there are many research investigations probing how modifiable and nonmodifiable factors impact AD symptom presentation. This review specifically focuses on how imaging modalities can be integrated into systems biology approaches using model mouse populations to link brain level functional and structural changes to disease onset and progression. Combining imaging and omics data promotes the classification of AD into subtypes and paves the way for precision medicine solutions to prevent and treat AD.

Paraquat Toxicogenetics: Strain-Related Reduction of Tyrosine Hydroxylase Staining in Substantia Nigra in Mice

Paraquat (PQ) is a putative risk factor for the development of sporadic Parkinson's disease. To model a possible genetic basis for individual differences in susceptibility to exposure to PQ, we recently examined the effects of paraquat on tyrosine hydroxylase (TH)-containing neurons in the substantia nigra pars compacta (SNc) of six members of the BXD family of mice (n = 2-6 per strain). We injected males with 5 mg/kg paraquat weekly three times. The density of TH+ neurons counted by immunocytochemistry at 200x in eight or more sections through the SNc is reduced in five of the six strains relative to control (N = 4 ± 2 mice per strain). TH+ loss ranged from 0 to 20% with an SEM of 1%. The heritability was estimated using standard ANOVA and jackknife resampling and is 0.37 ± 0.05 in untreated animals and 0.47 ± 0.04 in treated animals. These results demonstrate genetic modulation and GxE variation in susceptibility to PQ exposure and the loss of TH staining in the substantia nigra.

Mitochondrial cristae-remodeling protein OPA1 in POMC neurons couples Ca2+ homeostasis with adipose tissue lipolysis

Appropriate cristae remodeling is a determinant of mitochondrial function and bioenergetics and thus represents a crucial process for cellular metabolic adaptations. Here, we show that mitochondrial cristae architecture and expression of the master cristae-remodeling protein OPA1 in proopiomelanocortin (POMC) neurons, which are key metabolic sensors implicated in energy balance control, is affected by fluctuations in nutrient availability. Genetic inactivation of OPA1 in POMC neurons causes dramatic alterations in cristae topology, mitochondrial Ca2+ handling, reduction in alpha-melanocyte stimulating hormone (α-MSH) in target areas, hyperphagia, and attenuated white adipose tissue (WAT) lipolysis resulting in obesity. Pharmacological blockade of mitochondrial Ca2+ influx restores α-MSH and the lipolytic program, while improving the metabolic defects of mutant mice. Chemogenetic manipulation of POMC neurons confirms a role in lipolysis control. Our results unveil a novel axis that connects OPA1 in POMC neurons with mitochondrial cristae, Ca2+ homeostasis, and WAT lipolysis in the regulation of energy balance.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Mitochondrial function in development and disease

Mitochondria are organelles with vital functions in almost all eukaryotic cells. Often described as the cellular ‘powerhouses’ due to their essential role in aerobic oxidative phosphorylation, mitochondria perform many other essential functions beyond energy production. As signaling organelles, mitochondria communicate with the nucleus and other organelles to help maintain cellular homeostasis, allow cellular adaptation to diverse stresses, and help steer cell fate decisions during development. Mitochondria have taken center stage in the research of normal and pathological processes, including normal tissue homeostasis and metabolism, neurodegeneration, immunity and infectious diseases. The central role that mitochondria assume within cells is evidenced by the broad impact of mitochondrial diseases, caused by defects in either mitochondrial or nuclear genes encoding for mitochondrial proteins, on different organ systems. In this Review, we will provide the reader with a foundation of the mitochondrial ‘hardware’, the mitochondrion itself, with its specific dynamics, quality control mechanisms and cross-organelle communication, including its roles as a driver of an innate immune response, all with a focus on development, disease and aging. We will further discuss how mitochondrial DNA is inherited, how its mutation affects cell and organismal fitness, and current therapeutic approaches for mitochondrial diseases in both model organisms and humans.

Summary: Mitochondria have a plethora of functions beyond metabolism. This Review discusses the emerging and multifaceted roles of mitochondria in different model organisms and human disease biology.

Rare functional genetic variants in COL7A1, COL6A5, COL1A2 and COL5A2 frequently occur in Chiari Malformation Type 1

Chiari Malformation Type 1 (CM-1) is characterized by herniation of the cerebellar tonsils below the foramen magnum and the presence of headaches and other neurologic symptoms. Cranial bone constriction is suspected to be the most common biologic mechanism leading to CM-1. However, other mechanisms may also contribute, particularly in the presence of connective tissue disorders (CTDs), such as Ehlers Danlos Syndrome (EDS). Accumulating data suggest CM-1 with connective tissue disorders (CTD+) may have a different patho-mechanism and different genetic risk factors than CM-1 without CTDs (CTD-). To identify CM-1 genetic risk variants, we performed whole exome sequencing on a single large, multiplex family from Spain and targeted sequencing on a cohort of 186 unrelated adult, Caucasian females with CM-1. Targeted sequencing captured the coding regions of 21 CM-1 and EDS candidate genes, including two genes identified in the Spanish family. Using gene burden analysis, we compared the frequency of rare, functional variants detected in CM-1 cases versus publically available ethnically-matched controls from gnomAD. A secondary analysis compared the presence of rare variants in these genes between CTD+ and CTD- CM-1 cases. In the Spanish family, rare variants co-segregated with CM-1 in COL6A5, ADGRB3 and DST. A variant in COL7A1 was present in affected and unaffected family members. In the targeted sequencing analysis, rare variants in six genes (COL7A1, COL5A2, COL6A5, COL1A2, VEGFB, FLT1) were significantly more frequent in CM-1 cases compared to public controls. In total, 47% of CM-1 cases presented with rare variants in at least one of the four significant collagen genes and 10% of cases harbored variants in multiple significant collagen genes. Moreover, 26% of CM-1 cases presented with rare variants in the COL6A5 gene. We also identified two genes (COL7A1, COL3A1) for which the burden of rare variants differed significantly between CTD+ and CTD- CM-1 cases. A higher percentage of CTD+ patients had variants in COL7A1 compared to CTD+ patients, while CTD+ patients had fewer rare variants in COL3A1 than did CTD- patients. In summary, rare variants in several collagen genes are particularly frequent in CM-1 cases and those in COL6A5 co-segregated with CM-1 in a Spanish multiplex family. COL6A5 has been previously associated with musculoskeletal phenotypes, but this is the first association with CM-1. Our findings underscore the contribution of rare genetic variants in collagen genes to CM-1, and suggest that CM-1 in the presence and absence of CTD symptoms is driven by different genes.

Characterizing modifier genes of cardiac fibrosis phenotype in hypertrophic cardiomyopathy

BACKGROUND

Clinical phenotypes of hypertrophic cardiomyopathy (HCM) vary greatly even among patients with the same gene mutations. This variability is largely regulated by unidentified modifier loci. The purpose of the study is to identify modifier genes for cardiac fibrosis-a major phenotype of HCM-using the BXD family, a murine cohort.

METHODS

The relative severity of cardiac fibrosis was estimated by quantitation of cardiac collagen volume fraction (CCVF) across 66 members of the BXD family. Quantitative trait locus (QTL) mapping for cardiac fibrosis was done using GeneNetwork. Candidate modifier loci and genes associated with fibrosis were prioritized based on an explicit scoring system. Networks of correlation between fibrosis and cardiac transcriptomes were evaluated to generate causal models of disease susceptibility.

RESULTS

CCVF levels varied greatly within this family. Interval mapping identified a significant CCVF-related QTL on chromosome (Chr) 2 in males, and a significant QTL on Chr 4 Mb in females. The scoring system highlighted two strong candidate genes in the Chr 2 locus-Nek6 and Nr6a1. Both genes are highly expressed in the heart. Cardiac Nek6 mRNA levels are significantly correlated with CCVF. Nipsnap3b and Fktn are lead candidate genes for the Chr 4 locus, and both are also highly expressed in heart. Cardiac Nipsnap3b gene expression correlates well with CCVF.

CONCLUSION

Our study demonstrated that candidate modifier genes of cardiac fibrosis phenotype in HCM are different in males and females. Nek6 and Nr6a1 are strong candidates in males, while Nipsnap3b and Fktn are top candidates in females.

Explaining Divergent Observations Regarding Osteocalcin/GPRC6A Endocrine Signaling

A new schema proposes that the bone-derived osteocalcin (Ocn) peptide hormone activates the G-protein-coupled receptor GPRC6A to directly regulate glucose and fat metabolism in liver, muscle, and fat, and to stimulate the release of metabolism-regulating hormones, including insulin, fibroblast growth factor 21, glucagon-like peptide 1, testosterone, and interleukin 6. Ocn/GPRC6A activation has also been implicated in cancer progression. GPRC6A is activated by cations, amino acids, and testosterone. The multiligand specificity, the regulation of energy metabolism in diverse tissues, and the coordinated release of metabolically active hormones make the GPRC6A endocrine networks unique. Recently, the significance of Ocn/GPRCA has been questioned. There is a lack of metabolic abnormalities in newly created genetically engineered Ocn- and Gprc6a-deficient mouse models. There are also paradoxical observations that GPRC6A may function as a tumor suppressor. In addition, discordant published studies have cast doubt on the function of the most prevalent uniquely human GPRC6A-KGKY polymorphism. Explanations for these divergent findings are elusive. We provide evidence that the metabolic susceptibility of genetically engineered Ocn- and Gprc6a-deficient mice is influenced by environmental challenges and genetic differences in mouse strains. In addition, the GPRC6A-KGKY polymorphism appears to be a gain-of-function variant. Finally, alternatively spliced isoforms of GPRC6A may alter ligand specificity and signaling that modulate oncogenic effects. Thus, genetic, post-translational and environmental factors likely account for the variable results regarding the functions of GPRC6A in animal models. Pending additional information, GPRC6A should remain a potential therapeutic target for regulating energy and fat metabolism, hormone production, and cancer progression.

A Bioinformatics Perspective on the Links Between Tetraspanin-Enriched Microdomains and Cardiovascular Pathophysiology

Background: Tetraspanins and integrins are integral membrane proteins. Tetraspanins interact with integrins to modulate the dynamics of adhesion, migration, proliferation, and signaling in the form of membrane domains called tetraspanin-enriched microdomains (TEMs). TEMs also contain other cell adhesion proteins like immunoglobulin superfamily (IgSF) proteins and claudins. Cardiovascular functions of these TEM proteins have emerged and remain to be further revealed. Objectives: The aims of this study are to explore the roles of these TEM proteins in the cardiovascular system using bioinformatics tools and databases and to highlight the TEM proteins that may functionally associate with cardiovascular physiology and pathology. Methods: For human samples, three databases-GTEx, NCBI-dbGaP, and NCBI-GEO-were used for the analyses. The dbGaP database was used for GWAS analysis to determine the association between target genes and human phenotypes. GEO is an NCBI public repository that archives genomics data. GTEx was used for the analyses of tissue-specific mRNA expression levels and eQTL. For murine samples, GeneNetwork was used to find gene-phenotype correlations and gene-gene correlations of expression levels in mice. The analysis of cardiovascular data was the focus of this study. Results: Some integrins and tetraspanins, such as ITGA8 and Cd151, are highly expressed in the human cardiovascular system. TEM components are associated with multiple cardiovascular pathophysiological events in humans. GWAS and GEO analyses showed that human Cd82 and ITGA9 are associated with blood pressure. Data from mice also suggest that various cardiovascular phenotypes are correlated with integrins and tetraspanins. For instance, Cd82 and ITGA9, again, have correlations with blood pressure in mice. Conclusion: ITGA9 is related to blood pressure in both species. KEGG analysis also linked ITGA9 to metabolism and MAPK signaling pathway. This work provides an example of using integrated bioinformatics approaches across different species to identify the connections of structurally and/or functionally related molecules to certain categories of diseases.

A platform for experimental precision medicine: The extended BXD mouse family

The challenge of precision medicine is to model complex interactions among DNA variants, phenotypes, development, environments, and treatments. We address this challenge by expanding the BXD family of mice to 140 fully isogenic strains, creating a uniquely powerful model for precision medicine. This family segregates for 6 million common DNA variants-a level that exceeds many human populations. Because each member can be replicated, heritable traits can be mapped with high power and precision. Current BXD phenomes are unsurpassed in coverage and include much omics data and thousands of quantitative traits. BXDs can be extended by a single-generation cross to as many as 19,460 isogenic F1 progeny, and this extended BXD family is an effective platform for testing causal modeling and for predictive validation. BXDs are a unique core resource for the field of experimental precision medicine.

Mitochondrial unfolded protein response: An emerging pathway in human diseases

Mitochondrial unfolded protein response (UPR^mt) is a mitochondria stress response, which the transcriptional activation programs of mitochondrial chaperone proteins and proteases are initiated to maintain proteostasis in mitochondria. Additionally, the activation of UPR^mt delays aging and extends lifespan by maintaining mitochondrial proteostasis. Growing evidences suggests that UPR^mt plays an important role in diverse human diseases, especially ageing-related diseases. Therefore, this review focuses on the role of UPR^mt in ageing and ageing-related neurodegenerative diseases such as Alzheimer's disease, Huntington's disease and Parkinson's disease. The activation of UPR^mt and the high expression of UPR^mt components contribute to longevity extension. The activation of UPR^mt may ameliorate Alzheimer's disease, Parkinson's disease and Huntington's disease. Besides, UPR^mt is also involved in the occurrence and development of cancers and heart diseases. UPR^mt contributes to the growth, invasive and metastasis of cancers. UPR^mt has paradoxical roles in heart diseases. UPR^mt not only protects against heart damage, but may sometimes aggravates the development of heart diseases. Considering the pleiotropic actions of UPR^mt system, targeting UPR^mt pathway may be a potent therapeutic avenue for neurodegenerative diseases, cancers and heart diseases.

The mitochondrial unfolded protein response and its diverse roles in cellular stress

Mitochondrial function is centrally involved in many cellular processes, such as energy production, metabolism of nucleotides, amino acids, and lipids, calcium buffering, and regulation of cell death. Multiple mechanisms are engaged under conditions of mitochondrial dysfunction to restore cellular and, subsequently, systemic functions. The mitochondrial unfolded protein response is a homeostatic mechanism that has attracted a lot of interest recently and has been described in several organisms, including humans. The mitochondrial unfolded protein response serves as a first-line-of-defence mechanism against stress to restore mitochondrial proteostasis and functions. Here, we discuss the canonical mechanisms via which the mitochondrial unfolded protein response is activated under stress and examine recent evidence that links the response with other processes that promote survival and the recovery of the mitochondrial network (i.e. the integrated stress response and mitophagy).

Genetic Dissection of the Regulatory Mechanisms of Ace2 in the Infected Mouse Lung

Acute lung injury (ALI) is an important cause of morbidity and mortality after viral infections, including influenza A virus H1N1, SARS-CoV, MERS-CoV, and SARS-CoV-2. The angiotensin I converting enzyme 2 (ACE2) is a key host membrane-bound protein that modulates ALI induced by viral infection, pulmonary acid aspiration, and sepsis. However, the contributions of ACE2 sequence variants to individual differences in disease risk and severity after viral infection are not understood. In this study, we quantified H1N1 influenza-infected lung transcriptomes across a family of 41 BXD recombinant inbred strains of mice and both parents—C57BL/6J and DBA/2J. In response to infection Ace2 mRNA levels decreased significantly for both parental strains and the expression levels was associated with disease severity (body weight loss) and viral load (expression levels of viral NA segment) across the BXD family members. Pulmonary RNA-seq for 43 lines was analyzed using weighted gene co-expression network analysis (WGCNA) and Bayesian network approaches. Ace2 not only participated in virus-induced ALI by interacting with TNF, MAPK, and NOTCH signaling pathways, but was also linked with high confidence to gene products that have important functions in the pulmonary epithelium, including Rnf128, Muc5b, and Tmprss2. Comparable sets of transcripts were also highlighted in parallel studies of human SARS-CoV-infected primary human airway epithelial cells. Using conventional mapping methods, we determined that weight loss at two and three days after viral infection maps to chromosome X—the location of Ace2. This finding motivated the hierarchical Bayesian network analysis, which defined molecular endophenotypes of lung infection linked to Ace2 expression and to a key disease outcome. Core members of this Bayesian network include Ace2, Atf4, Csf2, Cxcl2, Lif, Maml3, Muc5b, Reg3g, Ripk3, and Traf3. Collectively, these findings define a causally-rooted Ace2 modulatory network relevant to host response to viral infection and identify potential therapeutic targets for virus-induced respiratory diseases, including those caused by influenza and coronaviruses.

The Genetic Dissection of Ace2 Expression Variation in the Heart of Murine Genetic Reference Population

Background: A high inflammatory and cytokine burden that induces vascular inflammation, myocarditis, cardiac arrhythmias, and myocardial injury is associated with a lethal outcome in COVID-19. The SARS-CoV-2 virus utilizes the ACE2 receptor for cell entry in a similar way to SARS-CoV. This study investigates the regulation, gene network, and associated pathways of ACE2 that may be involved in inflammatory and cardiovascular complications of COVID-19.

Methods: Cardiovascular traits were determined in the one of the largest mouse genetic reference populations: BXD recombinant inbred strains using blood pressure, electrocardiography, and echocardiography measurements. Expression quantitative trait locus (eQTL) mapping, genetic correlation, and functional enrichment analysis were used to identify Ace2 regulation, gene pathway, and co-expression networks.

Results: A wide range of variation was found in expression of Ace2 among the BXD strains. Levels of Ace2 expression are negatively correlated with cardiovascular traits, including systolic and diastolic blood pressure and P wave duration and amplitude. Ace2 co-expressed genes are significantly involved in cardiac- and inflammatory-related pathways. The eQTL mapping revealed that Cyld is a candidate upstream regulator for Ace2. Moreover, the protein–protein interaction (PPI) network analysis inferred several potential key regulators (Cul3, Atf2, Vcp, Jun, Ppp1cc, Npm1, Mapk8, Set, Dlg1, Mapk14, and Hspa1b) for Ace2 co-expressed genes in the heart.

Conclusions:Ace2 is associated with blood pressure, atrial morphology, and sinoatrial conduction in BXD mice. Ace2 co-varies with Atf2, Cyld, Jun, Mapk8, and Mapk14 and is enriched in the RAS, TGFβ, TNFα, and p38α signaling pathways, involved in inflammation and cardiac damage. We suggest that all these novel Ace2-associated genes and pathways may be targeted for preventive, diagnostic, and therapeutic purposes in cardiovascular damage in patients with systemic inflammation, including COVID-19 patients.

The extracellular matrix phenome across species

Extracellular matrices are essential for cellular and organismal function. Recent genome-wide and phenome-wide association studies started to reveal a broad spectrum of phenotypes associated with genetic variants. However, the phenome or spectrum of all phenotypes associated with genetic variants in extracellular matrix genes is unknown. Here, we analyzed over two million recorded genotype-to-phenotype relationships across multiple species to define their extracellular matrix phenomes. By using the previously defined matrisomes of humans, mice, zebrafish, Drosophila, and C. elegans, we found that the extracellular matrix phenome comprises of 3-10% of the entire phenome. Collagens (COL1A1, COL2A1) and fibrillin (FBN1) are each associated with >150 distinct phenotypes in humans, whereas collagen COL4A1, Wnt- and sonic hedgehog (shh) signaling are predominantly associated in other species. We determined the phenotypic fingerprints of matrisome genes and found that MSTN, CTSD, LAMB2, HSPG2, and COL11A2 and their corresponding orthologues have the most phenotypes across species. Out of the 42,551 unique matrisome genotype-to-phenotype relationships across the five species with defined matrisomes, we have constructed interaction networks to identify the underlying molecular components connecting with orthologues phenotypes and with novel phenotypes. Thus, our networks provide a framework to predict unassessed phenotypes and their potential underlying molecular interactions. These frameworks inform on matrisome genotype-to-phenotype relationships and potentially provide a sophisticated choice of biological model system to study human phenotypes and diseases.

Genome-wide transcriptome architecture in a mouse model of Gulf War Illness

Gulf War Illness (GWI) is thought to be a chronic neuroimmune disorder caused by in-theater exposure during the 1990-1991 Gulf War. There is a consensus that the illness is caused by exposure to insecticides and nerve agent toxicants. However, the heterogeneity in both development of disease and clinical outcomes strongly suggests a genetic contribution. Here, we modeled GWI in 30 BXD recombinant inbred mouse strains with a combined treatment of corticosterone (CORT) and diisopropyl fluorophosphate (DFP). We quantified transcriptomes from 409 prefrontal cortex samples. Compared to the untreated and DFP treated controls, the combined treatment significantly activated pathways such as cytokine-cytokine receptor interaction and TNF signaling pathway. Protein-protein interaction analysis defined 6 subnetworks for CORT + DFP, with the key regulators being Cxcl1, Il6, Ccnb1, Tnf, Agt, and Itgam. We also identified 21 differentially expressed genes having significant QTLs related to CORT + DFP, but without evidence for untreated and DFP treated controls, suggesting regions of the genome specifically involved in the response to CORT + DFP. We identified Adamts9 as a potential contributor to response to CORT + DFP and found links to symptoms of GWI. Furthermore, we observed a significant effect of CORT + DFP treatment on the relative proportion of myelinating oligodendrocytes, with a QTL on Chromosome 5. We highlight three candidates, Magi2, Sema3c, and Gnai1, based on their high expression in the brain and oligodendrocyte. In summary, our results show significant genetic effects of the CORT + DFP treatment, which mirrors gene and protein expression changes seen in GWI sufferers, providing insight into the disease and a testbed for future interventions.

The shaping of immunological responses through natural selection after the Roma Diaspora

The Roma people are the largest transnational ethnic minority in Europe and can be considered the last human migration of South Asian origin into the continent. They left Northwest India approximately 1,000 years ago, reaching the Balkan Peninsula around the twelfth century and Romania in the fourteenth century. Here, we analyze whole-genome sequencing data of 40 Roma and 40 non-Roma individuals from Romania. We performed a genome-wide scan of selection comparing Roma, their local host population, and a Northwestern Indian population, to identify the selective pressures faced by the Roma mainly after they settled in Europe. We identify under recent selection several pathways implicated in immune responses, among them cellular metabolism pathways known to be rewired after immune stimulation. We validated the interaction between PIK3-mTOR-HIF-1α and cytokine response influenced by bacterial and fungal infections. Our results point to a significant role of these pathways for host defense against the most prevalent pathogens in Europe during the last millennium.

Sex and Strain Variation in Initial Sensitivity and Rapid Tolerance to Δ9-Tetrahydrocannabinol

Background and Objectives: For cannabis and other drugs of abuse, initial response and/or tolerance to drug effects can predict later dependence and problematic use. Our objective is to identify sex and genetic (strain) differences in initial response and rapid tolerance to Δ9–tetrahydrocannabinol (THC), the main psychoactive ingredient in cannabis, between highly genetically divergent inbred mouse strains—C57BL/6J (B6) and DBA/2J (D2).

Experimental Approach: Sex and strain responses relative to baseline were quantified following daily exposure (i.p.) to 10 mg/kg THC or vehicle (VEH) over the course of 5 days. Dependent measures included hypothermia (decreased body temperature) and ataxia (decreased spontaneous activity in the open field), and antinociception (increase in tail withdrawal latency to a thermal stimulus). Initial sensitivity to THC was defined as the difference in response between baseline and day 1. Rapid tolerance to THC was defined as the difference in response between days 1 and 2.

Results: B6 exhibited greater THC-induced motor activity suppression and initial sensitivity to ataxia relative to the D2 strain. Females demonstrated greater levels of THC-induced hypothermia and initial sensitivity relative to males. Higher levels of THC-induced antinociception and initial sensitivity were observed for D2 relative to B6. Rapid tolerance to THC was observed for hypothermia and antinociception. Much less tolerance was observed for THC-induced ataxia. D2 exhibited rapid tolerance to THC-induced hypothermia and antinociception at time points associated with peak THC initial response. Likewise, at the peak initial THC response time point, females demonstrated greater levels of rapid tolerance to hypothermic effects relative to males.

Conclusions: Both sex and genetic factors drive variation in initial response and rapid tolerance to the ataxic, antinociceptive, and hypothermic effects of THC. As these traits directly result from THC activation of the cannabinoid receptor 1, gene variants between B6 and D2 in cannabinoid signaling pathways are likely to mediate strain differences in response to THC.

An Age-Related Hearing Protection Locus on Chromosome 16 of BXD Strain Mice

Inbred mouse models are widely used to study age-related hearing loss (AHL). Many genes associated with AHL have been mapped in a variety of strains. However, little is known about gene variants that have the converse function—protective genes that confer strong resistance to hearing loss. Previously, we reported that C57BL/6J (B6) and DBA/2J (D2) strains share a common hearing loss allele in Cdh23. The cadherin 23 (Cdh23) gene is a key contributor to early-onset hearing loss in humans. In this study, we tested hearing across a large family of 54 BXD strains generated from B6 to D2 crosses. Five of 54 strains maintain the normal threshold (20 dB SPL) even at 2 years old—an age at which both parental strains are essentially deaf. Further analyses revealed an age-related hearing protection (ahp) locus on chromosome 16 (Chr 16) at 57~76 Mb with a maximum LOD of 5.7. A small number of BXD strains at 2 years with good hearing correspond roughly to the percentage of humans who have good hearing at 90 years old. Further studies to define candidate genes in the ahp locus and related molecular mechanisms involved in age-related resilience or resistance to AHL are warranted.

Genetic Architecture Associated With Familial Short Stature

CONTEXT

Human height is an inheritable, polygenic trait under complex and multilocus genetic regulation. Familial short stature (FSS; also called genetic short stature) is the most common type of short stature and is insufficiently known.

OBJECTIVE

To investigate the FSS genetic profile and develop a polygenic risk predisposition score for FSS risk prediction.

DESIGN AND SETTING

The FSS participant group of Han Chinese ancestry was diagnosed by pediatric endocrinologists in Taiwan.

PATIENTS AND INTERVENTIONS

The genetic profiles of 1163 participants with FSS were identified by using a bootstrapping subsampling and genome-wide association studies (GWAS) method.

MAIN OUTCOME MEASURES

Genetic profile, polygenic risk predisposition score for risk prediction.

RESULTS

Ten novel genetic single nucleotide polymorphisms (SNPs) and 9 reported GWAS human height-related SNPs were identified for FSS risk. These 10 novel SNPs served as a polygenic risk predisposition score for FSS risk prediction (area under the curve: 0.940 in the testing group). This FSS polygenic risk predisposition score was also associated with the height reduction regression tendency in the general population.

CONCLUSION

A polygenic risk predisposition score composed of 10 genetic SNPs is useful for FSS risk prediction and the height reduction tendency. Thus, it might contribute to FSS risk in the Han Chinese population from Taiwan.