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.

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.

The TMEM43 S358L mutation affects cardiac, small intestine, and metabolic homeostasis in a knock-in mouse model.

The transmembrane protein 43 (TMEM43/LUMA) p.S358L mutation causes arrhythmogenic cardiomyopathy named as ARVC5, a fully penetrant disease with high risk of ventricular arrhythmias, sudden death, and heart failure. Male gender and vigorous exercise independently predicted deleterious outcome. Our systems genetics analysis revealed the importance of Tmem43 for cardiac and metabolic pathways associated with elevated lipid absorption from small intestine. This study sought to delineate gender-specific cardiac, intestinal, and metabolic phenotypes in vivo and investigate underlying pathophysiological mechanisms of S358L mutation. Serial echocardiography, surface electrocardiography (ECG), treadmill running, and body EchoMRI have been used in knock-in heterozygous (Tmem43WT/S358L), homozygous (Tmem43S358L), and wildtype (Tmem43WT) littermate mice. Electron microscopy, histology, immunohistochemistry, transcriptome, and protein analysis have been performed in cardiac and intestinal tissues. Systolic dysfunction was apparent in 3-mo-old Tmem43S358L and 6-mo-old Tmem43WT/S358L mutants. Both mutant lines displayed intolerance to acute stress at 6 mo of age, arrhythmias, fibro-fatty infiltration, and subcellular abnormalities in the myocardium. Microarray analysis found significantly differentially expressed genes between left ventricular (LV) and right ventricular (RV) myocardium. Mutants displayed diminished PPARG activities and significantly reduced TMEM43 and ß-catenin expression in the heart, whereas junctional plakoglobin (JUP) translocated into nuclei of mutant cardiomyocytes. Conversely, elongated villi, fatty infiltration, and overexpression of gut epithelial proliferation markers, ß-catenin and Ki-67, were evident in small intestine of mutants. We defined Tmem43 S358L-induced pathological effects on cardiac and intestinal homeostasis via distinctly disturbed WNT-ß-catenin and PPARG signaling thereby contributing to ARVC5 pathophysiology. Results suggest that cardiometabolic assessment in mutation carriers may be important for predictive and personalized care.NEW & NOTEWORTHY This manuscript describes the findings of our investigation of cardiac, small intestine, and metabolic features of Tmem43-S358L mouse model. By investigating interorgan pathologies, we uncovered multiple mechanisms of the S358L-induced disease, and these unique mechanisms likely appear to contribute to the disease pathogenesis. We hope our findings are important and novel and open new avenues in the hunting for additional diagnostic and therapeutic targets in subjects carrying TMEM43 mutation.

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.

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.

Acetylation of Checkpoint suppressor 1 enhances its stability and promotes the progression of triple-negative breast cancer.

Checkpoint suppressor 1 (CHES1), a transcriptional regulator, had been dysregulated in many types of malignancies including breast cancer, and its expression level is strongly associated with progression and prognosis of patients. However, the underlying regulatory mechanisms of CHES1 expression in the breast cancer and the effects of post-translational modifications (PTMs) on its functional performance remain to be fully investigated. Herein, we found that CHES1 had a high abundance in triple-negative breast cancer (TNBC) and its expression was tightly associated with malignant phenotype and poor outcomes of patients. Furthermore, we confirmed that CHES1 was an acetylated protein and its dynamic modification was mediated by p300 and HDAC1, and CHES1 acetylation enhanced its stability via decreasing its ubiquitination and degradation, which resulted in the high abundance of CHES1 in TNBC. RNA-seq and functional study revealed that CHES1 facilitated the activation of oncogenic genes and pathways leading to proliferation and metastasis of TNBC. Taken together, this research established a novel regulatory role of acetylation on the stability and activity of CHES1. The results demonstrate the significance of CHES1 acetylation and underlying mechanisms in the progression of TNBC, offering new potential candidate for molecular-targeted therapy in breast cancer.

Restrictive cardiomyopathy: from genetics and clinical overview to animal modeling.

Restrictive cardiomyopathy (RCM), a potentially devastating heart muscle disorder, is characterized by diastolic dysfunction due to abnormal muscle relaxation and myocardial stiffness resulting in restrictive filling of the ventricles. Diastolic dysfunction is often accompanied by left atrial or bi-atrial enlargement and normal ventricular size and systolic function. RCM is the rarest form of cardiomyopathy, accounting for 2-5% of pediatric cardiomyopathy cases, however, survival rates have been reported to be 82%, 80%, and 68% at 1-, 2-, and 5-years after diagnosis, respectively. RCM can be idiopathic, familial, or secondary to a systemic disorder, such as amyloidosis, sarcoidosis, and hereditary hemochromatosis. Approximately 30% of cases are familial RCM, and the genes that have been linked to RCM are cTnT, cTnI, MyBP-C, MYH7, MYL2, MYL3, DES, MYPN, TTN, BAG3, DCBLD2, LNMA, and FLNC. Increased Ca2+ sensitivity, sarcomere disruption, and protein aggregates are some of the few mechanisms of pathogenesis that have been revealed by studies utilizing cell lines and animal models. Additional exploration into the pathogenesis of RCM is necessary to create novel therapeutic strategies to reverse restrictive cardiomyopathic phenotypes.

Systems genetics analysis defines importance of TMEM43/LUMA for cardiac- and metabolic-related pathways.

Broad cellular functions and diseases including muscular dystrophy, arrhythmogenic right ventricular cardiomyopathy (ARVC5) and cancer are associated with transmembrane protein43 (TMEM43/LUMA). The study aimed to investigate biological roles of TMEM43 through genetic regulation, gene pathways and gene networks, candidate interacting genes, and up- or downstream regulators. Cardiac transcriptomes from 40 strains of recombinant inbred BXD mice and two parental strains representing murine genetic reference population (GRP) were applied for genetic correlation, functional enrichment, and coexpression network analysis using systems genetics approach. The results were validated in a newly created knock-in Tmem43-S358L mutation mouse model (Tmem43S358L) that displayed signs of cardiac dysfunction, resembling ARVC5 phenotype seen in humans. We found high Tmem43 levels among BXDs with broad variability in expression. Expression of Tmem43 highly negatively correlated with heart mass and heart rate among BXDs, whereas levels of Tmem43 highly positively correlated with plasma high-density lipoproteins (HDL). Through finding differentially expressed genes (DEGs) between Tmem43S358L mutant and wild-type (Tmem43WT) lines, 18 pathways (out of 42 found in BXDs GRP) that are involved in ARVC, hypertrophic cardiomyopathy, dilated cardiomyopathy, nonalcoholic fatty liver disease, Alzheimer's disease, Parkinson's disease, and Huntington's disease were verified. We further constructed Tmem43-mediated gene network, in which Ctnna1, Adcy6, Gnas, Ndufs6, and Uqcrc2 were significantly altered in Tmem43S358L mice versus Tmem43WT controls. Our study defined the importance of Tmem43 for cardiac- and metabolism-related pathways, suggesting that cardiovascular disease-relevant risk factors may also increase risk of metabolic and neurodegenerative diseases via TMEM43-mediated pathways.

Combining whole exome sequencing with in silico analysis and clinical data to identify candidate variants in pediatric left ventricular noncompaction.

BACKGROUND: Understanding the overall variant burden in pediatric patients with left ventricular noncompaction (LVNC) has clinical implications. Whole exome sequencing (WES) allows detection of coding variants in both candidate cardiomyopathy genes and those included on commercial panels. Other lines of evidence, including in silico analysis, are necessary to reduce the overwhelming number of variants to those most likely having a phenotypic impact. METHODS: Five families, including five pediatric probands with LVNC, 5 other affected, and 10 unaffected family members, had WES performed, followed by bioinformatics filtering and Sanger sequencing. Review of the HGMD, variant classification by ACMG guidelines, and clinical information were used to further refine complex genotypes. RESULTS: One nonsense and eleven missense variants were identified. In Family 1, affected siblings carried digenic heterozygous variants: E1350K-MYH7 and A276V-ANKRD1. The proband also carried heterozygous W143X-NRG1. Four affected members of Family 2 carried K184Q-MYH7 while unaffected members did not. In Family 3, homozygous A161T-MYH7 and heterozygous P4935T-OBSCN variants were identified in the proband with the latter being absent in his unaffected brother. In Family 4, proband's father and half-sibling have mild hypertrabeculation and carry T3796I-PLEC. The proband, carrying T3796I-PLEC and V2878A-OBSCN, demonstrated higher trabeculation burden. The proband in Family 5 carried four variants, R3247W-PLEC, C92Y-ERG, T1233M-NCOR2, and E54K-HIST1H4B. Application of ACMG criteria and clinical data revealed that W143X-NRG1, P4935T-OBSCN, and V2878A-OBSCN likely have no phenotypic role. CONCLUSIONS: We report nine variants, including novel T3796I-PLEC and biallelic A161T-MYH7, likely contributing to phenotypes ranging from asymptomatic hypertrabeculation to severe LVNC with heart failure.

Ace2 and Tmprss2 Expressions Are Regulated by Dhx32 and Influence the Gastrointestinal Symptoms Caused by SARS-CoV-2.

Studies showed that the gastrointestinal (GI) tract is one of the most important pathways for SARS-CoV-2 infection and coronavirus disease 2019 (COVID-19). As SARS-CoV-2 cellular entry depends on the ACE2 receptor and TMPRSS2 priming of the spike protein, it is important to understand the molecular mechanisms through which these two proteins and their cognate transcripts interact and influence the pathogenesis of COVID-19. In this study, we quantified the expression, associations, genetic modulators, and molecular pathways for Tmprss2 and Ace2 mRNA expressions in GI tissues using a systems genetics approach and the expanded family of highly diverse BXD mouse strains. The results showed that both Tmprss2 and Ace2 are highly expressed in GI tissues with significant covariation. We identified a significant expression quantitative trait locus on chromosome 7 that controls the expression of both Tmprss2 and Ace2. Dhx32 was found to be the strongest candidate in this interval. Co-expression network analysis demonstrated that both Tmprss2 and Ace2 were located at the same module that is significantly associated with other GI-related traits. Protein-protein interaction analysis indicated that hub genes in this module are linked to circadian rhythms. Collectively, our data suggested that genes with circadian rhythms of expression may have an impact on COVID-19 disease, with implications related to the timing and treatment of COVID-19.

Deficiency in nebulin repeats of sarcomeric nebulette is detrimental for cardiomyocyte tolerance to exercise and biomechanical stress.

The actin-binding sarcomeric nebulette (NEBL) protein provides efficient contractile flexibility via interaction with desmin intermediate filaments. NEBL gene mutations affecting the nebulin repeat (NR) domain are known to induce cardiomyopathy. The study aimed to explore the roles of NEBL in exercise and biomechanical stress response. We ablated exon3 encoding the first NR of Nebl and created global Neblex3-/ex3- knockout mice. Cardiac function, structure, and transcriptome were assessed before and after a 4-wk treadmill regimen. A Nebl-based exercise signaling network was constructed using systems genetics methods. H9C2 and neonatal rat cardiomyocytes (NRCs) expressing wild-type or mutant NEBL underwent cyclic mechanical strain. Neblex3-/ex3- mice demonstrated diastolic dysfunction with preserved systolic function at 6 mo of age. After treadmill running, 4-mo-old Neblex3-/ex3- mice developed concentric cardiac hypertrophy and left ventricular dilation compared with running Nebl+/+ and sedentary Neblex3-/ex3- mice. Disturbance of sarcomeric Z-disks and thin filaments architecture and disruption of intercalated disks and mitochondria were found in exercised Neblex3-/ex3- mice. A Nebl-based exercise signaling network included Csrp3, Des, Fbox32, Jup, Myh6, and Myh7. Disturbed expression of TM1, DES, JUP, ß-catenin, MLP, α-actinin2, and vinculin proteins was demonstrated. In H9C2 cells, NEBL was recruited into focal adhesions at 24-h poststrain and redistributed along with F-actin at 72-h poststrain, suggesting time-dependent redistribution of NEBL in response to strain. NEBL mutations cause desmin disorganization in NRCs upon stretch. We conclude that Nebl's NR ablation causes disturbed sarcomere, Z-disks, and desmin organization, and prevents NEBL redistribution to focal adhesions in cardiomyocytes, weakening cardiac tolerance to biomechanical stress.NEW & NOTEWORTHY We demonstrate that ablation of first nebulin-repeats of sarcomeric nebulette (Nebl) causes diastolic dysfunction in Neblex3-/ex3- mice. Exercise-induced development of diastolic dysfunction, cardiac hypertrophy and ventricular dilation in knockouts. This was associated with sarcomere disturbance, intercalated disks disruption, and mitochondrial distortion upon stress and altered expression of genes involved in Nebl-based stress network. We demonstrate that G202R and A592 mutations alter actin and desmin expression causing disorganization of desmin filaments upon cyclic strain.