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.

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.

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

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

The Combination of Quantitative Proteomics and Systems Genetics Analysis Reveals that PTN Is Associated with Sleep-Loss-Induced Cognitive Impairment.

Sleep loss is associated with cognitive dysfunction. However, the detailed mechanisms remain unclear. In this study, we established a para-chlorophenylalanine (PCPA)-induced insomniac mouse model with impaired cognitive function. Mass-spectrometry-based proteomics showed that the expression of 164 proteins was significantly altered in the hippocampus of the PCPA mice. To identify critical regulators among the potential markers, a transcriptome-wide association screening was performed in the BXD mice panel. Among the candidates, the expression of pleiotrophin (Ptn) was significantly associated with cognitive functions, indicating that Ptn-mediates sleep-loss-induced cognitive impairment. Gene co-expression analysis further revealed the potential mechanism by which Ptn mediates insomnia-induced cognitive impairment via the MAPK signaling pathway; that is, the decreased secretion of Ptn induced by insomnia leads to reduced binding to Ptprz1 on the postsynaptic membrane with the activation of the MAPK pathway via Fos and Nr4a1, further leading to the apoptosis of neurons. In addition, Ptn is genetically trans-regulated in the mouse hippocampus and implicated in neurodegenerative diseases in human genome-wide association studies. Our study provides a novel biomarker for insomnia-induced cognitive impairment and a new strategy for seeking neurological biomarkers by the integration of proteomics and systems genetics.

Joint analysis of WES and RNA-Seq identify signature genes related to metastasis in prostate cancer.

Prostate cancer (PCa) has a certain degree of heritability, and metastasis occurs as cancer progresses. However, its underlying mechanism remains largely unknown. We sequenced four cases of cancer without metastasis, four metastatic cancer, and four benign hyperplasia tissues as controls. A total of 1839 damaging mutations were identified. Pathway analysis, gene clustering, and weighted gene co-expression network analysis were employed to find characteristics associated with metastasis. Chr19 had the most mutation density and 1p36 had the highest mutation frequency across the genome. These mutations occurred in 1630 genes, including the most frequently mutated genes TTN and PLEC, and dozens of metastasis-related genes, such as FOXA1, NCOA1, CD34, and BRCA2. Ras signalling and arachidonic acid metabolism were uniquely enriched in metastatic cancer. Gene programmes 10 and 11 showed the signatures indicating the occurrence of metastasis better. A module (135 genes) was specifically associated with metastasis. Of them, 67.41% reoccurred in program 10, with 26 genes further retained as the signature genes related to PCa metastasis, including AGR3, RAPH1, SOX14, DPEP1, and UBL4A. Our study provides new molecular perspectives on PCa metastasis. The signature genes and pathways could be served as potential therapeutic targets for metastasis or cancer progression.

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.

HOXD8 suppresses renal cell carcinoma growth by upregulating SHMT1 expression.

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

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.

Silencing DTX3L Inhibits the Progression of Cervical Carcinoma by Regulating PI3K/AKT/mTOR Signaling Pathway.

Cervical carcinoma (CC) is the second most prevalent gynecologic cancer in females across the world. To obtain a better understanding of the mechanisms underlying the development of CC, high-resolution label-free mass spectrometry was performed on CC and adjacent normal tissues from eight patients. A total of 2631 proteins were identified, and 46 significant differently expressed proteins (DEPs) were found between CC and normal tissues (p < 0.01, fold change >10 or <0.1). Ingenuity pathway analysis revealed that the majority of the proteins were involved in the regulation of eIF4 and p70S6K signaling and mTOR signaling. Among 46 DEPs, Integrinß6 (ITGB6), PPP1CB, TMPO, PTGES3 (P23) and DTX3L were significantly upregulated, while Desmin (DES) was significantly downregulated in CC tissues compared with the adjacent normal tissues. In in vivo and in vitro experiments, DTX3L knockdown suppressed CC cell proliferation, migration, invasion and xenograft tumorigenesis, and enhanced cell apoptosis. Combination of silencing DTX3L and cisplatin treatment induced higher apoptosis percentage compared to cisplatin treatment alone. Moreover, DTX3L silencing inhibited the PI3K/AKT/mTOR signal pathway. Thus, our results suggested DTX3L could regulate CC progression through the PI3K/AKT/mTOR signal pathway and is potentially a novel biomarker and therapeutic target for CC.

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.

Genomic Basis for Individual Differences in Susceptibility to the Neurotoxic Effects of Diesel Exhaust.

Air pollution is a known environmental health hazard. A major source of air pollution includes diesel exhaust (DE). Initially, research on DE focused on respiratory morbidities; however, more recently, exposures to DE have been associated with neurological developmental disorders and neurodegeneration. In this study, we investigated the effects of sub-chronic inhalation exposure to DE on neuroinflammatory markers in two inbred mouse strains and both sexes, including whole transcriptome examination of the medial prefrontal cortex. We exposed aged male and female C57BL/6J (B6) and DBA/2J (D2) mice to DE, which was cooled and diluted with HEPA-filtered compressed air for 2 h per day, 5 days a week, for 4 weeks. Control animals were exposed to HEPA-filtered air on the same schedule as DE-exposed animals. The prefrontal cortex was harvested and analyzed for proinflammatory cytokine gene expression (Il1ß, Il6, Tnfα) and transcriptome-wide response by RNA-seq. We observed differential cytokine gene expression between strains and sexes in the DE-exposed vs. control-exposed groups for Il1ß, Tnfα, and Il6. For RNA-seq, we identified 150 differentially expressed genes between air and DE treatment related to natural killer cell-mediated cytotoxicity per Kyoto Encyclopedia of Genes and Genomes pathways. Overall, our data show differential strain-related effects of DE on neuroinflammation and neurotoxicity and demonstrate that B6 are more susceptible than D2 to gene expression changes due to DE exposures than D2. These results are important because B6 mice are often used as the default mouse model for DE studies and strain-related effects of DE neurotoxicity warrant expanded studies.

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.

Identifying modifier genes for hypertrophic cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) severity greatly varies among patients even with the same HCM gene mutations. This variation is largely regulated by modifier gene(s), which, however, remain largely unknown. The current study is aimed to identify modifier genes using BXD strains, a large murine genetic reference population (GRP) derived from crosses between C57BL/6 J (B6) and D2 DBA/2 J (D2) mice. D2 mice natualy carrythe genetic basis and phenotypes of HCM. METHODS: Myocardial hypertrophy, the major phenotype of HCM, was determined by cardiomyocyte size on cardiac sections in 30 BXD strains, and their parental B6 and D2 strains and morphometric analysis was performed. Quantitative Trait Locus (QTL) mapping for cardiomyocyte sizes was conducted with WebQTL in GeneNetwork. Correlation of cardiomyocyte size and cardiac gene expression in BXDs accessed from GeneNetwork were evaluated. QTL candidate genes associated with cardiomyocyte sizes were prioritized based on the score system. RESULTS: Cardiomyocyte size varied significantly among BXD strains. Interval mapping on cardiomyocyte size data showed a significant QTL on chromosome (Chr) 2 at 66- 73.5 Mb and a suggestive QTL on Chr 5 at 20.9-39.7 Mb. Further score system revealed a high QTL score for Xirp2 in Chr 2. Xirp2 encodes xin actin-binding repeat containing 2, which is highly expressed in cardiac tissue and associate with cardiomyopathy and heart failure. In Chr5 QTL, Nos3, encoding nitric oxide synthase 3, received the highest score, which is significantly correlated with cardiomyocyte size. CONCLUSION: These results indicate that Xirp2 and Nos3 serve as novel candidate modifier genes for myocardial hypertrophy in HCM. These candidate genes will be validated in our future studies.

System biology analysis reveals the role of voltage-dependent anion channel in mitochondrial dysfunction during non-alcoholic fatty liver disease progression into hepatocellular carcinoma.

Non-alcoholic fatty liver disease (NAFLD) is one of the most common causes of hepatocellular carcinoma (HCC), but the underlying mechanisms behind the correlation of NAFLD with HCC are unclear. We aimed to uncover the genes and potential mechanisms that drive this progression. This study uncovered the genes and potential mechanisms through a multiple 'omics integration approach. Quantitative proteomics combined with phenotype-association analysis was performed. To investigate the potential mechanisms, a comprehensive transcriptome/lipidome/phenome-wide association analysis was performed in genetic reference panel BXD mice strains. The quantitative proteomics combined with phenotype-association results showed that VDAC1 was significantly increased in tumor tissues and correlated with NAFLD-related traits. Gene co-expression network analysis indicated that VDAC1 is involved in mitochondria dysfunction in the tumorigenic/tumor progression. The association between VDAC1 and mitochondria dysfunction can be explained by the fact that VDAC1 was associated with mitochondria membrane lipids cardiolipin (CL) composition shift. VDAC1 was correlated with the suppression of mature specie CL(LLLL) and elevation level of nascent CL species. Such profiling shift was supported by the significant positive correlation between VDAC1 and PTPMT1, as well as negative correlation with CL remodeling enzyme Tafazzin (TAZ). This study confirmed that the expression of VADC1 was dysregulated in NAFLD-driven HCC and associated with NAFLD progression. The VDAC1-driven mitochondria dysfunction is associated with cardiolipin composition shift, which causes alteration of mitochondria membrane properties.

A systems genetics approach to revealing the Pdgfb molecular network of the retina.

Purpose: Platelet-derived growth factor (PDGF) signaling is well known to be involved in vascular retinopathies. Among the PDGF family, the subunit B (PDGFB) protein is considered a promising therapeutic target. This study aimed to identify the genes and potential pathways through which PDGFB affects retinal phenotypes by using a systems genetics approach. Methods: Gene expression data had been previously generated in a laboratory for the retinas of 75 C57BL/6J(B6) X DBA/2J (BXD) recombinant inbred (RI) strains. Using this data, the genetic correlation method was used to identify genes correlated to Pdgfb. A correlation between intraocular pressure (IOP) and Pdgfb was calculated based on the Pearson correlation coefficient. A gene set enrichment analysis and the STRING database were used to evaluate gene function and to construct protein-protein interaction (PPI) networks. Results: Pdgfb was a cis-regulated gene in the retina; its expression had a significant correlation with IOP (r = 0.305; p value = 0.012). The expression levels of 2,477 genes also had significant correlations with Pdgfb expressions (p<0.05), among which Atf4 was the most positively correlated (r = 0.628; p value = 1.29e-10). Thus, Atf4 was highly expressed in the retina and shared the transcription factor (TF) Hnf4a binding site with Pdgfb. Gene Ontology and a pathway analysis revealed that Pdgfb and its covariates were highly involved in mitogen-activated protein kinase (MAPK) and vascular endothelial growth factor (VEGF) pathways. A generated gene network indicated that Pdgfb was directly connected to and interacted with other genes with similar biologic functions. Conclusions: A systems genetics analysis revealed that Pdgfb had significant interactions with Atf4 and other genes in MAPK and VEGF pathways, through which Pdgfb was important in maintaining retina function. These findings provided basic information regarding the Pdgfb regulation mechanism and potential therapy for vascular retinopathies.

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.

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 Dissection of Hypertrophic Cardiomyopathy with Myocardial RNA-Seq.

Hypertrophic cardiomyopathy (HCM) is an inherited disorder of the myocardium, and pathogenic mutations in the sarcomere genes myosin heavy chain 7 (MYH7) and myosin-binding protein C (MYBPC3) explain 60%-70% of observed clinical cases. The heterogeneity of phenotypes observed in HCM patients, however, suggests that novel causative genes or genetic modifiers likely exist. Here, we systemically evaluated RNA-seq data from 28 HCM patients and nine healthy controls with pathogenic variant identification, differential expression analysis, and gene co-expression and protein-protein interaction network analyses. We identified 43 potential pathogenic variants in 19 genes in 24 HCM patients. Genes with more than one variant included the following: MYBPC3, TTN, MYH7, PSEN2, and LDB3. A total of 2538 protein-coding genes, six microRNAs (miRNAs), and 1617 long noncoding RNAs (lncRNAs) were identified differentially expressed between the groups, including several well-characterized cardiomyopathy-related genes (ANKRD1, FHL2, TGFB3, miR-30d, and miR-154). Gene enrichment analysis revealed that those genes are significantly involved in heart development and physiology. Furthermore, we highlighted four subnetworks: mtDNA-subnetwork, DSP-subnetwork, MYH7-subnetwork, and MYBPC3-subnetwork, which could play significant roles in the progression of HCM. Our findings further illustrate that HCM is a complex disease, which results from mutations in multiple protein-coding genes, modulation by non-coding RNAs and perturbations in gene networks.

Investigating a downstream gene of Gpnmb using the systems genetics method.

Purpose: Glaucoma is characterized by optic nerve damage and retinal ganglion cell loss. The glycoprotein neuromedin B-associated (Gpnmb) gene is well-known to be involved in the glaucoma disease process. The purpose of this study is to identify a downstream gene through which Gpnmb affects the glaucoma phenotypes using a systems genetics approach. Methods: Retinal gene expression data for the BXD recombinant inbred (RI) strains (n=75) have previously been generated in our laboratory for a glaucoma study, and these data were used for genetic and bioinformatics analysis. Expression quantitative trait locus (eQTL) mapping and genetic correlation methods were used to identify a gene downstream of Gpnmb. Gene-set enrichment analysis was used to evaluate gene function and to construct coexpression networks. Results: The level of Gpnmb expression is associated with a highly statistically significant cis-eQTL. Stanniocalcin 1 (Stc1) has a significant trans-eQTL mapping to the Gpnmb locus. The expression of Gpnmb and Stc1 is highly correlated in the retina and other tissues, as well as with glaucoma-related phenotypes. Gene Ontology and pathway analysis showed that Stc1 and its covariates are highly associated with apoptosis, oxidative stress, and mitochondrial activity. A generated gene network indicated that Gpnmb and Stc1 are directly connected to and interact with other genes with similar biologic functions. Conclusions: These results suggest that Stc1 may be a downstream candidate of Gpnmb, and that both genes interact with other genes in a network to develop glaucoma pathogenesis through mechanisms such as apoptosis and oxidative stress.

A multiplex SNP genotyping by allele-specificspecific PCR based on stem-loop and universal fluorescent primers of Chr1daxin mice.

Single nucleotide polymorphisms (SNPs) are one of the most common markers in mammals. Rapid, accurate, and multiplex typing of SNPs is critical for subsequent biological and genetic research. In this study, we have developed a novel method for multiplex genotyping SNPs in mice. The method involves allele-specific PCR amplification of genomic DNA with two stem-loop primers accompanied by two different universal fluorescent primers. Blue and green fluorescent signals were conveniently detected on a DNA sequencer. We verified four SNPs of 65 mice based on the novel method, and it is well suited for multiplex genotyping as it requires only one reaction per sample in a single tube with multiplex PCR. The use of universal fluorescent primers greatly reduces the cost of designing different fluorescent probes for each SNP. Therefore, this method can be applied to many biological and genetic studies, such as multiple candidate gene testing, genome-wide association study, pharmacogenetics, and medical diagnostics.