Systems Genetics Analyses Reveals Key Genes Related to Behavioral Traits in the Striatum of CFW Mice.

The striatum plays a central role in directing many complex behaviors ranging from motor control to action choice and reward learning. In our study, we used 55 male CFW mice with rapid decay linkage disequilibrium to systematically mine the striatum-related behavioral functional genes by analyzing their striatal transcriptomes and 79 measured behavioral phenotypic data. By constructing a gene coexpression network, we clustered the genes into 13 modules, with most of them being positively correlated with motor traits. Based on functional annotations as well as Fisher's exact and hypergeometric distribution tests, brown and magenta modules were identified as core modules. They were significantly enriched for striatal-related functional genes. Subsequent Mendelian randomization analysis verified the causal relationship between the core modules and dyskinesia. Through the intramodular gene connectivity analysis, Adcy5 and Kcnma1 were identified as brown and magenta module hub genes, respectively. Knock outs of both Adcy5 and Kcnma1 lead to motor dysfunction in mice, and KCNMA1 acts as a risk gene for schizophrenia and smoking addiction in humans. We also evaluated the cellular composition of each module and identified oligodendrocytes in the striatum to have a positive role in motor regulation.

Identification of hub genes and potential molecular mechanisms related to drug sensitivity in acute myeloid leukemia based on machine learning.

Background: Acute myeloid leukemia (AML) is the most common form of leukemia among adults and is characterized by uncontrolled proliferation and clonal expansion of hematopoietic cells. There has been a significant improvement in the treatment of younger patients, however, prognosis in the elderly AML patients remains poor. Methods: We used computational methods and machine learning (ML) techniques to identify and explore the differential high-risk genes (DHRGs) in AML. The DHRGs were explored through multiple in silico approaches including genomic and functional analysis, survival analysis, immune infiltration, miRNA co-expression and stemness features analyses to reveal their prognostic importance in AML. Furthermore, using different ML algorithms, prognostic models were constructed and validated using the DHRGs. At the end molecular docking studies were performed to identify potential drug candidates targeting the selected DHRGs. Results: We identified a total of 80 DHRGs by comparing the differentially expressed genes derived between AML patients and normal controls and high-risk AML genes identified by Cox regression. Genetic and epigenetic alteration analyses of the DHRGs revealed a significant association of their copy number variations and methylation status with overall survival (OS) of AML patients. Out of the 137 models constructed using different ML algorithms, the combination of Ridge and plsRcox maintained the highest mean C-index and was used to build the final model. When AML patients were classified into low- and high-risk groups based on DHRGs, the low-risk group had significantly longer OS in the AML training and validation cohorts. Furthermore, immune infiltration, miRNA coexpression, stemness feature and hallmark pathway analyses revealed significant differences in the prognosis of the low- and high-risk AML groups. Drug sensitivity and molecular docking studies revealed top 5 drugs, including carboplatin and austocystin-D that may significantly affect the DHRGs in AML. Conclusion: The findings from the current study identified a set of high-risk genes that may be used as prognostic and therapeutic markers for AML patients. In addition, significant use of the ML algorithms in constructing and validating the prognostic models in AML was demonstrated. Although our study used extensive bioinformatics and machine learning methods to identify the hub genes in AML, their experimental validations using knock-out/-in methods would strengthen our findings.

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

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

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

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

5-methylcytosine RNA modification regulators-based patterns and features of immune microenvironment in acute myeloid leukemia.

Acute myeloid leukemia (AML) is a highly heterogeneous malignant disease of the blood cell. The current therapies for AML are unsatisfactory and the molecular mechanisms underlying AML are unclear. 5-methylcytosine (m5C) is an important posttranscriptional modification of mRNA, and is involved in the regulation of mRNA stability, translation, and other aspects of RNA metabolism. However, based on our knowledge of published literature, the role of the m5C regulators has not been explored in AML till date. In this study, we clarified the expression and gene variants of m5C regulators in AML and found that most m5C regulators were differentially expressed and correlated with disease prognosis. We also found that the methylation status of certain m5C regulators (e.g., DNMT3A, DNMT3B) affects the survival of AML patients. Two m5C modification subtypes, and high- and low-risk subgroups identified based on the expression of m5C regulators showed significant differences in the prognosis as well as immune cell infiltration. In addition, most of the m5C regulators were found to be correlated with miRNA expression in AML, as well as IC50 values of many drugs. The miRNA and GSVA analysis were used to identify the different miRNAs and KEGG or hallmark pathways between high- and low-risk subgroups. We also built a prognostic model based on m5C regulators, which was validated by two GSE databases. To verify the reliability of our analysis and conclusions, qPCR was used to identify the expressions of m5C regulators between normal and AML. In summary, we comprehensively explored the molecular characteristics of m5C regulators and built a prognostic model in AML. We proposed new mechanistic insights into the role of m5C in multiple databases and clinical data, which may pave novel ways for the development of therapeutic strategies.

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.

Comparative Root Transcriptome Profiling and Gene Regulatory Network Analysis between Eastern and Western Carrot (Daucus carota L.) Cultivars Reveals Candidate Genes for Vascular Tissue Patterning.

Carrot (Daucus carota L.) is a highly consumed vegetable rich in carotenoids, known for their potent antioxidant, anti-inflammatory, and immune-protecting properties. While genetic and molecular studies have largely focused on wild and Western carrot cultivars (cvs), little is known about the evolutionary interactions between closely related Eastern and Western cvs. In this study, we conducted comparative transcriptome profiling of root tissues from Eastern (UHSBC-23-1) and Western (UHSBC-100) carrot cv. to better understand differentially expressed genes (DEGs) associated with storage root development and vascular cambium (VC) tissue patterning. Through reference-guided TopHat mapping, we achieved an average mapping rate of 73.87% and identified a total of 3544 DEGs (p < 0.05). Functional annotation and gene ontology classification revealed 97 functional categories, including 33 biological processes, 19 cellular components, 45 metabolic processes, and 26 KEGG pathways. Notably, Eastern cv. exhibited enrichment in cell wall, plant-pathogen interaction, and signal transduction terms, while Western cv. showed dominance in photosynthesis, metabolic process, and carbon metabolism terms. Moreover, constructed gene regulatory network (GRN) for both cvs. obtained orthologs with 1222 VC-responsive genes of Arabidopsis thaliana. In Western cv, GRN revealed VC-responsive gene clusters primarily associated with photosynthetic processes and carbon metabolism. In contrast, Eastern cv. exhibited a higher number of stress-responsive genes, and transcription factors (e.g., MYB15, WRKY46, AP2/ERF TF connected via signaling pathways with NAC036) were identified as master regulators of xylem vessel differentiation and secondary cell wall thickening. By elucidating the comparative transcriptome profiles of Eastern and Western cvs. for the first time, our study provides valuable insights into the differentially expressed genes involved in root development and VC tissue patterning. The identification of key regulatory genes and their roles in these processes represents a significant advancement in our understanding of the evolutionary relations and molecular mechanisms underlying secondary growth of carrot and regulation by vascular cambium.

Exploring Early-Stage Retinal Neurodegeneration in Murine Pigmentary Glaucoma: Insights From Gene Networks and miRNA Regulation Analyses.

Purpose: Glaucoma is a group of heterogeneous optic neuropathies characterized by the progressive degeneration of retinal ganglion cells. However, the underlying mechanisms have not been understood completely. We aimed to elucidate the genetic network associated with the development of pigmentary glaucoma with DBA/2J (D2) mouse model of glaucoma and corresponding genetic control D2-Gpnmb (D2G) mice carrying the wild type (WT) Gpnmb allele. Methods: Retinas isolated from 13 D2 and 12 D2G mice were subdivided into 2 age groups: pre-onset (1-6 months: samples were collected at approximately 1-2, 2-4, and 5-6 months) and post-onset (7-15 months: samples were collected at approximately 7-9, 10-12, and 13-15 months) glaucoma were compared. Differential gene expression (DEG) analysis and gene-set enrichment analyses were performed. To identify micro-RNAs (miRNAs) that target Gpnmb, miRNA expression levels were correlated with time point matched mRNA expression levels. A weighted gene co-expression network analysis (WGCNA) was performed using the reference BXD mouse population. Quantitative real-time PCR (qRT-PCR) was used to validate Gpnmb and miRNA expression levels. Results: A total of 314 and 86 DEGs were identified in the pre-onset and post-onset glaucoma groups, respectively. DEGs in the pre-onset glaucoma group were associated with the crystallin gene family, whereas those in the post-onset group were related to innate immune system response. Of 1329 miRNAs predicted to target Gpnmb, 3 miRNAs (miR-125a-3p, miR-3076-5p, and miR-214-5p) were selected. A total of 47 genes demonstrated overlapping with the identified DEGs between D2 and D2G, segregated into their time-relevant stages. Gpnmb was significantly downregulated, whereas 2 out of 3 miRNAs were significantly upregulated (P < 0.05) in D2 mice at both 3-and 10-month time points. Conclusions: These findings suggest distinct gene-sets involved in pre-and post-glaucoma in the D2 mouse. We identified three miRNAs regulating Gpnmb in the development of murine pigmentary glaucoma.

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.

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

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

Identification of Potential Biomarkers for Group I Pulmonary Hypertension Based on Machine Learning and Bioinformatics Analysis.

A number of processes and pathways have been reported in the development of Group I pulmonary hypertension (Group I PAH); however, novel biomarkers need to be identified for a better diagnosis and management. We employed a robust rank aggregation (RRA) algorithm to shortlist the key differentially expressed genes (DEGs) between Group I PAH patients and controls. An optimal diagnostic model was obtained by comparing seven machine learning algorithms and was verified in an independent dataset. The functional roles of key DEGs and biomarkers were analyzed using various in silico methods. Finally, the biomarkers and a set of key candidates were experimentally validated using patient samples and a cell line model. A total of 48 key DEGs with preferable diagnostic value were identified. A gradient boosting decision tree algorithm was utilized to build a diagnostic model with three biomarkers, PBRM1, CA1, and TXLNG. An immune-cell infiltration analysis revealed significant differences in the relative abundances of seven immune cells between controls and PAH patients and a correlation with the biomarkers. Experimental validation confirmed the upregulation of the three biomarkers in Group I PAH patients. In conclusion, machine learning and a bioinformatics analysis along with experimental techniques identified PBRM1, CA1, and TXLNG as potential biomarkers for Group I PAH.

Differential Expression of MicroRNAs and Predicted Drug Target in Amyotrophic Lateral Sclerosis.

ALS (Amyotrophic Lateral Sclerosis) is a rare type of neurodegenerative disease. It shows progressive degradation of motor neurons in the brain and spinal cord. At present, there is no treatment available that can completely cure ALS. The available treatments can only increase a patient's life span by a few months. Recently, microRNAs (miRNAs), a sub-class of small non-coding RNAs have been shown to play an essential role in the diagnosis, prognosis, and therapy of ALS. Our study focuses on analyzing differential miRNA profiles and predicting drug targets in ALS using bioinformatics and computational approach. The study identifies eight highly differentially expressed miRNAs in ALS patients, four of which are novel. We identified 42 hub genes for these eight highly expressed miRNAs with Amyloid Precursor Protein (APP) as a candidate gene among them for highly expressed down-regulated miRNA, hsa-miR-455-3p using protein-protein interaction network and Cytoscape analysis. A novel association has been found between hsa-miR-455-3p/APP/serotonergic pathway using KEGG pathway analysis. Also, molecular docking studies have revealed curcumin as a potential drug target that may be used for the treatment of ALS. Thus, the present study has identified four novel miRNA biomarkers: hsa-miR-3613-5p, hsa-miR-24, hsa-miR-3064-5p, and hsa-miR-4455. There is a formation of a novel axis, hsa-miR-455-3p/APP/serotonergic pathway, and curcumin is predicted as a potential drug target for ALS.

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis.

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

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.

EWI2 and its relatives in Tetraspanin-enriched membrane domains regulate malignancy.

Experimental studies on immunoglobulin superfamily (IgSF) member EWI2 reveal that it suppresses a variety of solid malignant tumors including brain, lung, skin, and prostate cancers in animal models and inhibits tumor cell movement and growth in vitro. While EWI2 appears to support myeloid leukemia in mouse models and maintain leukemia stem cells. Bioinformatics analyses suggest that EWI2 gene expression is downregulated in glioblastoma but upregulated in melanoma, pancreatic cancer, and liver cancer. The mechanism of action for EWI2 is linked to its inhibition of growth factor receptors and cell adhesion proteins through its associated tetraspanin-enriched membrane domains (TEMDs), by altering the cell surface clustering and endolysosome trafficking/turnover of these transmembrane proteins. Recent studies also show that EWI2 modulates the nuclear translocation of ERK and TFEB to change the activities of these gene expression regulators. For EWI2 relatives including FPRP, IgSF3, and CD101, although their roles in malignant diseases are not fully clear and remain to be determined experimentally, FPRP and IgSF3 likely promote the progression of solid malignant tumors while CD101 seems to modulate immune cells of tumor microenvironment. Distinctive from other tumor regulators, the impacts of EWI subfamily members on solid malignant tumors are likely to be context dependent. In other words, the effect of a given EWI subfamily member on a tumor probably depends on the molecular network and composition of TEMDs in that tumor. Collectively, EWI2 and its relatives are emerged as important regulators of malignant diseases with promising potentials to become anti-cancer therapeutics and cancer therapy targets.

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.

A systems genetics approach delineates the role of Bcl2 in leukemia pathogenesis.

Leukemia is a group of malignancies of the blood forming tissues, and is characterized by the uncontrolled proliferation of blood cells. In the United States, it accounts for approximately 3.5% and 4% of all cancer-related incidences and mortalities, respectively. The current study aimed to explore the role of Bcl2 and associated genes in leukemia pathogenesis using a systems genetics approach. The transcriptome data from BXD Recombinant Inbred (RI) mice was analyzed to identify the expression of Bcl2 in myeloid cells. eQTL mapping was performed to select the potential chromosomal region and subsequently identify the candidate gene modulating the expression of Bcl2. Furthermore, gene enrichment and protein-protein interaction (PPI) analyses of the Bcl2-coexpressed genes were performed to demonstrate the role of Bcl2 in leukemia pathogenesis. The Bcl2-coexpressed genes were found to be enriched in various hematopoietic system related functions, and multiple pathways related to signaling, immune response, and cancer. The PPI network analysis demonstrated direct interaction of hematopoietic function related genes, such as Bag3, Bak1, Bcl2l11, Bmf, Mapk9, Myc, Ppp2r5c, and Ppp3ca with Bcl2. The eQTL mapping identified a 4.5 Mb genomic region on chromosome 11, potentially regulating the expression of Bcl2. A multi-criteria filtering process identified Top2a, among the genes located in the mapped locus, as the best candidate upstream regulator for Bcl2 expression variation. Hence, the current study provides better insights into the role of Bcl2 in leukemia pathogenesis and demonstrates the significance of our approach in gaining new knowledge on leukemia. Furthermore, our findings from the PPI network analysis and eQTL mapping provide supporting evidence of leukemia-associated genes, which can be further explored for their functional importance in leukemia. DATA AVAILABILITY: The myeloid cell transcriptomic data of the BXD mice used in this study can be accessed through our GeneNetwork (http://www.genenetwork.org) with the accession number of GN144.

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.