Cis-expression Quantitative Trait Loci Research Articles

Integration of genetic variants and gene network for drug repurposing in colorectal cancer

Even though many genetic risk loci for human diseases have been identified and comprehensively cataloged, strategies to guide clinical research by integrating the extensive results of genetic studies and biological resources are still limited. Moreover, integrative analyses that provide novel insights into disease biology are expected to be especially useful for drug discovery. Herein, we used text mining of genetic studies on colorectal cancer (CRC) and assigned biological annotations to identified risk genes in order to discover novel drug targets and potential drugs for repurposing. Risk genes for CRC were obtained from PubMed text mining, and for each gene, six functional and bioinformatic annotations were analyzed. The annotations include missense mutations, cis-expression quantitative trait loci (cis-eQTL), molecular pathway analyses, protein-protein interactions (PPIs), a genetic overlap with knockout mouse phenotypes, and primary immunodeficiency (PID). We then prioritized the biological risk candidate genes according to a scoring system of the six functional annotations. Each functional annotation was assigned one point, and those genes with a score ≥2 were designated “biological CRC risk genes”. Using this method, we revealed 82 biological CRC risk genes, which were mapped to 128 genes in an expanded PPI network. Further utilizing DrugBank and the Therapeutic Target Database, we found 21 genes in our list that are targeted by 166 candidate drugs. Based on data from ClinicalTrials.gov and literature review, we found four known target genes with six drugs for clinical treatment in CRC, and three target genes with nine drugs supported by previous preclinical results in CRC. Additionally, 12 genes are targeted by 32 drugs approved for other indications, which can possibly be repurposed for CRC treatment. Finally, analysis from Connectivity Map (CMap) showed that 18 drugs have a high potential for CRC.

Evaluation of relationship between DNA methyltransferase 3 β gene and the risk of hip osteoarthritis: A case-control study based on a Han Chinese population.

Osteoarthritis (OA) is a chronic degenerative joint disease. Early studies have indicated that genetic and environmental factors contribute to the risk of OA. However, the etiology of OA remains unknown. Our study aimed to evaluate the association of DNMT3B gene with the risk of hip OA in Han Chinese individuals. A total of 2070 subjects were recruited into the study, including 658 patients with hip OA and 1412 healthy controls. Twelve tag single nucleotide polymorphisms (SNPs) were selected and genotyped in our samples. Genetic associations between DNMT3B gene and the risk of hip OA were examined at both the single marker and haplotype levels. Cis-expression quantitative trait loci signals that achieve genome-wide significance of targeted SNPs from multiple types of human tissues were extracted from the GTEx database. Significant signals were identified for SNP rs2424905 in 4 genetic models. The T allele was significantly associated with an increased risk of hip OA (odds ratio=1.53; 95% CI=1.28-1.83). The T allele was also significantly associated with higher Kellgren-Lawrence grade in the patients with hip OA (χ2 =32.70, P=1.37×10-6 ). Moreover, SNP rs2424905 was significantly associated with the gene expression level of multiple genes, including DNMT3B, C20orf203, COMMD7, EFCAB8, MAPRE1, and RP5-1085F17.3, from several types of human tissues. Our results indicated that rs2424905 of DNMT3B gene contributed to the risk of hip OA and its clinical severity in a Han Chinese population. These findings suggested that rs2424905 of DNMT3B could be a promising genetic marker to assess susceptibility to hip OA in Han Chinese populations.

Genome-wide transcriptome analysis identifies novel dysregulated genes implicated in Alzheimer's pathology.

Abnormal gene expression patterns may contribute to the onset and progression of late-onset Alzheimer's disease (LOAD). We performed transcriptome-wide meta-analysis (N = 1440) of blood-based microarray gene expression profiles as well as neuroimaging and cerebrospinal fluid (CSF) endophenotype analysis. We identified and replicated five genes (CREB5, CD46, TMBIM6, IRAK3, and RPAIN) as significantly dysregulated in LOAD. The most significantly altered gene, CREB5, was also associated with brain atrophy and increased amyloid beta (Aβ) accumulation, especially in the entorhinal cortex region. cis-expression quantitative trait loci mapping analysis of CREB5 detected five significant associations (P < 5 × 10-8 ), where rs56388170 (most significant) was also significantly associated with global cortical Aβ deposition measured by [18 F]Florbetapir positron emission tomography and CSF Aβ1-42 . RNA from peripheral blood indicated a differential gene expression pattern in LOAD. Genes identified have been implicated in biological processes relevant to Alzheimer's disease. CREB, in particular, plays a key role in nervous system development, cell survival, plasticity, and learning and memory.

Multiple sclerosis risk variants regulate gene expression in innate and adaptive immune cells.

At least 200 single-nucleotide polymorphisms (SNPs) are associated with multiple sclerosis (MS) risk. A key function that could mediate SNP-encoded MS risk is their regulatory effects on gene expression. We performed microarrays using RNA extracted from purified immune cell types from 73 untreated MS cases and 97 healthy controls and then performed Cis expression quantitative trait loci mapping studies using additive linear models. We describe MS risk expression quantitative trait loci associations for 129 distinct genes. By extending these models to include an interaction term between genotype and phenotype, we identify MS risk SNPs with opposing effects on gene expression in cases compared with controls, namely, rs2256814 MYT1 in CD4 cells (q = 0.05) and rs12087340 RF00136 in monocyte cells (q = 0.04). The rs703842 SNP was also associated with a differential effect size on the expression of the METTL21B gene in CD8 cells of MS cases relative to controls (q = 0.03). Our study provides a detailed map of MS risk loci that function by regulating gene expression in cell types relevant to MS.

Multifaceted genome-wide study identifies novel regulatory loci in SLC22A11 and ZNF45 for body mass index in Indians.

Obesity, a risk factor for multiple diseases (e.g. diabetes, hypertension, cancers) originates through complex interactions between genes and prevailing environment (food habit and lifestyle) that varies across populations. Indians exhibit a unique obesity phenotype with high abdominal adiposity for a given body weight compared to matched white populations suggesting presence of population-specific genetic and environmental factors influencing obesity. However, Indian population-specific genetic contributors for obesity have not been explored yet. Therefore, to identify potential genetic contributors, we performed a two-staged genome-wide association study (GWAS) for body mass index (BMI), a common measure to evaluate obesity in 5973 Indian adults and the lead findings were further replicated in 1286 Indian adolescents. Our study revealed novel association of variants-rs6913677 in BAI3 gene (p = 1.08 × 10-8) and rs2078267 in SLC22A11 gene (p = 4.62 × 10-8) at GWAS significance, and of rs8100011 in ZNF45 gene (p = 1.04 × 10-7) with near GWAS significance. As genetic loci may dictate the phenotype through modulation of epigenetic processes, we overlapped genetic data of identified signals with their DNA methylation patterns in 236 Indian individuals and performed methylation quantitative trait loci (meth-QTL) analysis. Further, functional roles of discovered variants and underlying genes were speculated using publicly available gene regulatory databases (ENCODE, JASPAR, GeneHancer, GTEx). The identified variants in BAI3 and SLC22A11 genes were found to dictate methylation patterns at unique CpGs harboring critical cis-regulatory elements. Further, BAI3, SLC22A11 and ZNF45 variants were located in repressive chromatin, active enhancer, and active chromatin regions, respectively, in human subcutaneous adipose tissue in ENCODE database. Additionally, these genomic regions represented potential binding sites for key transcription factors implicated in obesity and/or metabolic disorders. Interestingly, GTEx portal identify rs8100011 as a robust cis-expression quantitative trait locus (cis-eQTL) in subcutaneous adipose tissue (p = 1.6 × 10-7), and ZNF45 gene expression in skeletal muscle of Indian subjects showed an inverse correlation with BMI indicating its possible role in obesity. In conclusion, our study discovered 3 novel population-specific functional genetic variants (rs6913677, rs2078267, rs8100011) in 2 novel (SLC22A11 and ZNF45) and 1 earlier reported gene (BAI3) for BMI in Indians. Our study decodes key genomic loci underlying obesity phenotype in Indians that may serve as prospective drug targets in future.

The evolving role of TonEBP as an immunometabolic stress protein.

Tonicity-responsive enhancer-binding protein (TonEBP), which is also known as nuclear factor of activated T cells 5 (NFAT5), was discovered 20 years ago as a transcriptional regulator of the cellular response to hypertonic (hyperosmotic salinity) stress in the renal medulla. Numerous studies since then have revealed that TonEBP is a pleiotropic stress protein that is involved in a range of immunometabolic diseases. Some of the single-nucleotide polymorphisms (SNPs) in TONEBP introns are cis-expression quantitative trait loci that affect TONEBP transcription. These SNPs are associated with increased risk of type 2 diabetes mellitus, diabetic nephropathy, inflammation, high blood pressure and abnormal plasma osmolality, indicating that variation in TONEBP expression might contribute to these phenotypes. In addition, functional studies have shown that TonEBP is involved in the pathogenesis of rheumatoid arthritis, atherosclerosis, diabetic nephropathy, acute kidney injury, hyperlipidaemia and insulin resistance, autoimmune diseases (including type 1 diabetes mellitus and multiple sclerosis), salt-sensitive hypertension and hepatocellular carcinoma. These pathological activities of TonEBP are in contrast to the protective actions of TonEBP in response to hypertonicity, bacterial infection and DNA damage induced by genotoxins. An emerging theme is that TonEBP is a stress protein that mediates the cellular response to a range of pathological insults, including excess caloric intake, inflammation and oxidative stress.

Enhancer Domains Predict Gene Pathogenicity and Inform Gene Discovery in Complex Disease

Non-coding transcriptional regulatory elements are critical for controlling the spatiotemporal expression of genes. Here, we demonstrate that the sizes and number of enhancers linked to a gene reflect its disease pathogenicity. Moreover, genes with redundant enhancer domains are depleted of cis-acting genetic variants that disrupt gene expression, and they are buffered against the effects of disruptive non-coding mutations. Our results demonstrate that dosage-sensitive genes have evolved a robustness to the disruptive effects of genetic variation by expanding their regulatory domains. This solves a puzzle about why genes associated with human disease are depleted of cis-eQTLs (cis-expression quantitative trait loci), suggesting that this relationship might complicate gene identification in causal genome-wide association studies (GWASs) using eQTL information, and establishes a framework for identifying non-coding regulatory variation with phenotypic consequences.

Genomic dissection of 43 serum urate-associated loci provides multiple insights into molecular mechanisms of urate control.

High serum urate is a prerequisite for gout and associated with metabolic disease. Genome-wide association studies (GWAS) have reported dozens of loci associated with serum urate control; however, there has been little progress in understanding the molecular basis of the associated loci. Here, we employed trans-ancestral meta-analysis using data from European and East Asian populations to identify 10 new loci for serum urate levels. Genome-wide colocalization with cis-expression quantitative trait loci (eQTL) identified a further five new candidate loci. By cis- and trans-eQTL colocalization analysis, we identified 34 and 20 genes, respectively, where the causal eQTL variant has a high likelihood that it is shared with the serum urate-associated locus. One new locus identified was SLC22A9 that encodes organic anion transporter 7 (OAT7). We demonstrate that OAT7 is a very weak urate-butyrate exchanger. Newly implicated genes identified in the eQTL analysis include those encoding proteins that make up the dystrophin complex, a scaffold for signaling proteins and transporters at the cell membrane; MLXIP that, with the previously identified MLXIPL, is a transcription factor that may regulate serum urate via the pentose-phosphate pathway and MRPS7 and IDH2 that encode proteins necessary for mitochondrial function. Functional fine mapping identified six loci (RREB1, INHBC, HLF, UBE2Q2, SFMBT1 and HNF4G) with colocalized eQTL containing putative causal SNPs. This systematic analysis of serum urate GWAS loci identified candidate causal genes at 24 loci and a network of previously unidentified genes likely involved in control of serum urate levels, further illuminating the molecular mechanisms of urate control.

Integrating DNA sequencing and transcriptomic data for association analyses of low-frequency variants and lipid traits.

Transcriptome-wide association studies (TWAS) integrate genome-wide association studies (GWAS) and transcriptomic data to showcase their improved statistical power of identifying gene-trait associations while, importantly, offering further biological insights. TWAS have thus far focused on common variants as available from GWAS. Compared with common variants, the findings for or even applications to low-frequency variants are limited and their underlying role in regulating gene expression is less clear. To fill this gap, we extend TWAS to integrating whole genome sequencing data with transcriptomic data for low-frequency variants. Using the data from the Framingham Heart Study, we demonstrate that low-frequency variants play an important and universal role in predicting gene expression, which is not completely due to linkage disequilibrium with the nearby common variants. By including low-frequency variants, in addition to common variants, we increase the predictivity of gene expression for 79% of the examined genes. Incorporating this piece of functional genomic information, we perform association testing for five lipid traits in two UK10K whole genome sequencing cohorts, hypothesizing that cis-expression quantitative trait loci, including low-frequency variants, are more likely to be trait-associated. We discover that two genes, LDLR and TTC22, are genome-wide significantly associated with low-density lipoprotein cholesterol based on 3203 subjects and that the association signals are largely independent of common variants. We further demonstrate that a joint analysis of both common and low-frequency variants identifies association signals that would be missed by testing on either common variants or low-frequency variants alone.

Genome-wide association study in Chinese cohort reveals a new susceptibility locus within the 12q13.13 region for hypospadias

Hypospadias risk–associated gene variants have been reported in populations of European descent using genome-wide association studies (GWAS). There is little known at present about any possible hypospadias risk associations in Han Chinese populations. To systematically investigate hypospadias risk–associated gene variants in Chinese patients, we performed the first GWAS in a Han Chinese cohort consisting of 200 moderate-severe hypospadias cases and 1008 unaffected controls. Suggestive loci (p<1×10−4) were replicated in 118 cases and 383 controls, as well as in a second independent validation population of 137 cases and 190 controls. Regulatory were then conducted for the functional analyses of candidate variants. We identified rs11170516 with the risk allele G within the SP1/SP7 region that was independently associated with moderate-severe Hypospadias [SP1/SP7, rs11170516, Pcombine = 3.5×10-9, odds ratio (OR) = 1.96 (1.59-2.44)]. Results also suggested that rs11170516 is associated with the expression of SP1 as a cis-expression quantitative trait loci (cis-eQTL). We performed the first GWAS of moderate-severe hypospadias in a Han Chinese cohort, and identified one novel susceptibility cis-acting regulatory locus at 12q13.13, which may regulate a variety of hypospadias-related pathways by affecting proximal SP1 gene expression. This study complements known common hypospadias risk-associated variants and provides the possible role of cis-acting regulatory variant in causing hypospadias.

Genome-wide association study in Chinese cohort identifies one novel hypospadias risk associated locus at 12q13.13

BackgroundHypospadias risk–associated gene variants have been reported in populations of European descent using genome-wide association studies (GWASs). There is little known at present about any possible hypospadias risk associations in Han Chinese populations.MethodsTo systematically investigate hypospadias risk–associated gene variants in Chinese patients, we performed the first GWAS in a Han Chinese cohort consisting of 197 moderate-severe hypospadias cases and 933 unaffected controls. Suggestive loci (p < 1 × 10− 4) were replicated in 118 cases and 383 controls, as well as in a second independent validation population of 137 cases and 190 controls. Regulatory and protein-protein interactions (PPIs) were then conducted for the functional analyses of candidate variants.ResultsWe identified rs11170516 with the risk allele G within the SP1/SP7 region that was independently associated with moderate-severe hypospadias [SP1/SP7, rs11170516, Pcombine = 3.5 × 10− 9, odds ratio (OR) = 1.96 (1.59–2.44)]. Results also suggested that rs11170516 is associated with the expression of SP1 as a cis-expression quantitative trait locus (cis-eQTL). Protein SP1 could affect the risk of hypospadias via PPIs.ConclusionsWe performed the first GWAS of moderate-severe hypospadias in a Han Chinese cohort, and identified one novel susceptibility cis-acting regulatory locus at 12q13.13, which may regulate a variety of hypospadias-related pathways by affecting proximal SP1 gene expression and subsequent PPIs. This study complements known common hypospadias risk-associated variants and provides the possible role of cis-acting regulatory variant in causing hypospadias.

Adipose Tissue Gene Expression Associations Reveal Hundreds of Candidate Genes for Cardiometabolic Traits

Genome-wide association studies (GWASs) have identified thousands of genetic loci associated with cardiometabolic traits including type 2 diabetes (T2D), lipid levels, body fat distribution, and adiposity, although most causal genes remain unknown. We used subcutaneous adipose tissue RNA-seq data from 434 Finnish men from the METSIM study to identify 9,687 primary and 2,785 secondary cis-expression quantitative trait loci (eQTL; <1 Mb from TSS, FDR < 1%). Compared to primary eQTL signals, secondary eQTL signals were located further from transcription start sites, had smaller effect sizes, and were less enriched in adipose tissue regulatory elements compared to primary signals. Among 2,843 cardiometabolic GWAS signals, 262 colocalized by LD and conditional analysis with 318 transcripts as primary and conditionally distinct secondary cis-eQTLs, including some across ancestries. Of cardiometabolic traits examined for adipose tissue eQTL colocalizations, waist-hip ratio (WHR) and circulating lipid traits had the highest percentage of colocalized eQTLs (15% and 14%, respectively). Among alleles associated with increased cardiometabolic GWAS risk, approximately half (53%) were associated with decreased gene expression level. Mediation analyses of colocalized genes and cardiometabolic traits within the 434 individuals provided further evidence that gene expression influences variant-trait associations. These results identify hundreds of candidate genes that may act in adipose tissue to influence cardiometabolic traits.

Smoking Modifies the Genetic Risk for Early-Onset Periodontitis

Periodontitis has low-prevalence, highly severe disease manifestations with an early onset and rapid progression. The diagnosis is based on severe destruction of the alveolar bone in adolescents and young adults. Genetic susceptibility variants and smoking are well-established risk factors, but their interactions in modifying disease susceptibility have not been studied. We aimed to identify genetic risk variants of early-onset periodontitis that unmask their effects on tobacco smoke exposure. To this end, we analyzed 79,780,573 common variants in 741 northwest Europeans diagnosed to have >30% bone loss at >2 teeth before 35 y of age, using imputed genotypes of the OmniExpress BeadChip. Never versus ever smokers were compared in a logistic regression analysis via a case-only approach. To explore the effect of tobacco smoke on the expression of the G×S-associated genes, cultures of primary gingival fibroblasts (n = 9) were exposed to cigarette smoke extract, and transcripts were quantified by reverse transcription polymerase chain reaction. We identified 16 loci for which our analysis suggested an association with G×S increased disease risk (P < 5 × 10−5). Nine loci had previously been reported to be associated with spirometric measures of pulmonary function by an earlier G×S genome-wide association study. Genome-wide significant cis expression quantitative trait loci were reported for G×S-associated single-nucleotide polymorphisms at ST8SIA1 and SOST, indicating a causal role of these genes in tobacco-related etiopathology. Notably, SOST is a negative regulator of bone growth, and ST8SIA1 has a role in tissue remodeling. Cigarette smoke extract significantly altered the expression of 2 associated genes: SSH1 (P = 5 × 10−07), which is required for NF-κB activation and innate immune responses to bacterial invasion, and ST8SIA1 (P = 0.0048). We conclude that the genetic predisposition to early-onset periodontitis is in part triggered by smoking and that tobacco smoke directly affects the expression of genes involved in bone homeostasis, tissue repair, and immune response.

Genetic determinants controlling maize rubisco activase gene expression and a comparison with rice counterparts

BackgroundRubisco activase (RCA) regulates the activity of Rubisco and is a key enzyme of photosynthesis. RCA expression was widely reported to affect plant photosynthesis and crop yield, but the molecular basis of natural variation in RCA expression in a wide range of maize materials has not been fully elucidated.ResultsIn this study, correlation analysis in approximately 200 maize inbred lines revealed a significantly positive correlation between the expression of maize RCA gene ZmRCAβ and grain yield. A genome-wide association study revealed both cis-expression quantitative trait loci (cis-eQTLs) and trans-eQTLs underlying the expression of ZmRCAβ, with the latter playing a more important role. Further allele mining and genetic transformation analysis showed that a 2-bp insertion and a 14-bp insertion in the promoter of ZmRCAβ conferred increased gene expression. Because rice is reported to have higher RCA gene expression than does maize, we subsequently compared the genetic factors underlying RCA gene expression between maize and rice. The promoter activity of the rice RCA gene was shown to be stronger than that of the maize RCA gene, suggesting that replacing the maize RCA gene promoter with that of the rice RCA gene would improve the expression of RCA in maize.ConclusionOur results revealed two DNA polymorphisms regulating maize RCA gene ZmRCAβ expression, and the RCA gene promoter activity of rice was stronger than that of maize. This work increased understanding of the genetic mechanism that underlies RCA gene expression and identify new targets for both genetic engineering and selection for maize yield improvement.

Identifying Putative Susceptibility Genes and Evaluating Their Associations with Somatic Mutations in Human Cancers

Genome-wide association studies (GWASs) have identified hundreds of genetic risk variants for human cancers. However, target genes for the majority of risk loci remain largely unexplored. It is also unclear whether GWAS risk-loci-associated genes contribute to mutational signatures and tumor mutational burden (TMB) in cancer tissues. We systematically conducted cis-expression quantitative trait loci (cis-eQTL) analyses for 294 GWAS-identified variants for six major types of cancer-colorectal, lung, ovary, prostate, pancreas, and melanoma-by using transcriptome data from the Genotype-Tissue Expression (GTEx) Project, the Cancer Genome Atlas (TCGA), and other public data sources. By using integrative analysis strategies, we identified 270 candidate target genes, including 99 with previously unreported associations, for six cancer types. By analyzing functional genomic data, our results indicate that 180 genes (66.7% of 270) had evidence of cis-regulation by putative functional variants via proximal promoter or distal enhancer-promoter interactions. Together with our previously reported associations for breast cancer risk, our results show that 24 genes are shared by at least two cancer types, including four genes for both breast and ovarian cancer. By integrating mutation data from TCGA, we found that expression levels of 33 and 66 putative susceptibility genes were associated with specific mutational signatures and TMB of cancer-driver genes, respectively, at a Bonferroni-corrected p < 0.05. Together, these findings provide further insight into our understanding of how genetic risk variants might contribute to carcinogenesis through the regulation of susceptibility genes that are related to the biogenesis of somatic mutations.

Abstract 1595: Functional genomic analysis of 18q21.1 locus identifies potential functional variants and genes for colorectal cancer

Abstract Colorectal cancer (CRC) is one of the most commonly diagnosed cancers in the world. Genome-wide association studies have identified approximately 50 genetic risk loci for CRC. We have previously identified SNP rs7229639, located in the second intronic region of the SMAD7 gene at locus 18q21.1, in relation to CRC risk in Asians. In this study, we performed functional genomic analysis to identify potential functional variants and target genes for this association. Functional annotation analysis using epigenetic data from ENCODE and Roadmap projects showed functional evidence of promoter and/or enhancer activities for four highly correlated potential functional SNPs, including rs7229639, rs77544449, rs112540150 and rs57668004. Using in vitro luciferase reporter assays in colon cancer HCT116 and PKO cells, we showed that the risk alleles of rs7229639 and rs77544449 increase promoter or enhancer activities, compared to reference alleles. The risk alleles of both SNPs also affect the ability of interacting with nuclear proteins in the electrophoretic mobility shift assays, compared to reference alleles. To identify candidate target genes for this risk locus, we performed cis-expression quantitative trait loci (eQTL) analysis using transcriptome data from the Genotype-Tissue Expression (GTEx) project, The Cancer Genome Atlas (TCGA) and Colonmics. We identified putative target gene LIPG, located 500 kb downstream of rs7229639, in all three datasets. No associations were observed for SMAD7, or eight other genes in the region, via eQTL analysis. To investigate the biological functions of LIPG gene in colorectal cancer, we knock-downed LIPG gene expression using siRNA in HCT116 and PKO cells, which resulted in a significant decrease in cell proliferation, colony formation, migration and invasion in both HCT116 and PKO cells. These results suggest that LIPG gene plays a significant role in colorectal tumorigenesis via disruption of cell behaviors. Taken together, results from our study suggest that the association between the rs7229639/18q21.1 locus and CRC risk may be mediated through SNPs that regulate expression of the LIPG gene. Citation Format: Jifeng Wang, Xingyi Guo, Shimian Qu, Joshua A. Bauer, Jirong Long, Wei Zheng, Qiuyin Cai. Functional genomic analysis of 18q21.1 locus identifies potential functional variants and genes for colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1595.

Investigation of betaine as a novel psychotherapeutic for schizophrenia.

BackgroundBetaine is known to act against various biological stresses and its levels were reported to be decreased in schizophrenia patients. We aimed to test the role of betaine in schizophrenia pathophysiology, and to evaluate its potential as a novel psychotherapeutic.MethodsUsing Chdh (a gene for betaine synthesis)-deficient mice and betaine-supplemented inbred mice, we assessed the role of betaine in psychiatric pathophysiology, and its potential as a novel psychotherapeutic, by leveraging metabolomics, behavioral-, transcriptomics and DNA methylation analyses.FindingsThe Chdh-deficient mice revealed remnants of psychiatric behaviors along with schizophrenia-related molecular perturbations in the brain. Betaine supplementation elicited genetic background-dependent improvement in cognitive performance, and suppressed methamphetamine (MAP)-induced behavioral sensitization. Furthermore, betaine rectified the altered antioxidative and proinflammatory responses induced by MAP and in vitro phencyclidine (PCP) treatments. Betaine also showed a prophylactic effect on behavioral abnormality induced by PCP. Notably, betaine levels were decreased in the postmortem brains from schizophrenia, and a coexisting elevated carbonyl stress, a form of oxidative stress, demarcated a subset of schizophrenia with “betaine deficit-oxidative stress pathology”. We revealed the decrease of betaine levels in glyoxylase 1 (GLO1)-deficient hiPSCs, which shows elevated carbonyl stress, and the efficacy of betaine in alleviating it, thus supporting a causal link between betaine and oxidative stress conditions. Furthermore, a CHDH variant, rs35518479, was identified as a cis-expression quantitative trait locus (QTL) for CHDH expression in postmortem brains from schizophrenia, allowing genotype-based stratification of schizophrenia patients for betaine efficacy.InterpretationThe present study revealed the role of betaine in psychiatric pathophysiology and underscores the potential benefit of betaine in a subset of schizophrenia.FundThis study was supported by the Strategic Research Program for Brain Sciences from AMED (Japan Agency for Medical Research and Development) under Grant Numbers JP18dm0107083 and JP19dm0107083 (TY), JP18dm0107129 (MM), JP18dm0107086 (YK), JP18dm0107107 (HY), JP18dm0107104 (AK) and JP19dm0107119 (KH), by the Grant-in-Aid for Scientific Research on Innovative Areas from the MEXT under Grant Numbers JP18H05435 (TY), JP18H05433 (AH.-T), JP18H05428 (AH.-T and TY), and JP16H06277 (HY), and by JSPS KAKENHI under Grant Number JP17H01574 (TY). In addition, this study was supported by the Collaborative Research Project of Brain Research Institute, Niigata University under Grant Numbers 2018–2809 (YK) and RIKEN Epigenetics Presidential Fund (100214–201801063606-340120) (TY).

Tissue-specific sex differences in human gene expression.

Despite extensive sex differences in human complex traits and disease, the male and female genomes differ only in the sex chromosomes. This implies that most sex-differentiated traits are the result of differences in the expression of genes that are common to both sexes. While sex differences in gene expression have been observed in a range of different tissues, the biological mechanisms for tissue-specific sex differences (TSSDs) in gene expression are not well understood. A total of 30 640 autosomal and 1021 X-linked transcripts were tested for heterogeneity in sex difference effect sizes in n = 617 individuals across 40 tissue types in Genotype-Tissue Expression (GTEx). This identified 65 autosomal and 66 X-linked TSSD transcripts (corresponding to unique genes) at a stringent significance threshold. Results for X-linked TSSD transcripts showed mainly concordant direction of sex differences across tissues and replicate previous findings. Autosomal TSSD transcripts had mainly discordant direction of sex differences across tissues. The top cis-expression quantitative trait loci (eQTLs) across tissues for autosomal TSSD transcripts are located a similar distance away from the nearest androgen and estrogen binding motifs and the nearest enhancer, as compared to cis-eQTLs for transcripts with stable sex differences in gene expression across tissue types. Enhancer regions that overlap top cis-eQTLs for TSSD transcripts, however, were found to be more dispersed across tissues. These observations suggest that androgen and estrogen regulatory elements in a cis region may play a common role in sex differences in gene expression, but TSSD in gene expression may additionally be due to causal variants located in tissue-specific enhancer regions.

Genome-Wide Association between the 2q33.1 Locus and Intracranial Aneurysm Susceptibility: An Updated Meta-Analysis Including 18,019 Individuals.

Previous genome-wide association studies did not show a consistent association between the BOLL gene (rs700651, 2q33.1) and intracranial aneurysm (IA) susceptibility. We aimed to perform an updated meta-analysis for the potential IA-susceptibility locus in large-scale multi-ethnic populations. We conducted a systematic review of studies identified by an electronic search from January 1990 to March 2019. The overall estimates of the “G” allele of rs700651, indicating IA susceptibility, were calculated under the fixed- and random-effect models using the inverse-variance method. Subsequent in silico function and cis-expression quantitative trait loci (cis-eQTL) analyses were performed to evaluate biological functions and genotype-specific expressions in human tissues. We included 4513 IA patients and 13,506 controls from five studies with seven independent populations: three European-ancestry, three Japanese, and one Korean population. The overall result showed a genome-wide significance threshold between rs700651 and IA susceptibility after controlling for study heterogeneity (OR = 1.213, 95% CI: 1.135–1.296). Subsequent cis-eQTL analysis showed significant genome-wide expressions in three human tissues, i.e., testis (p = 8.04 × 10−15 for ANKRD44), tibial nerves (p = 3.18 × 10−10 for SF3B1), and thyroid glands (p = 4.61 × 10−9 for SF3B1). The rs700651 common variant of the 2q33.1 region may be involved in genetic mechanisms that increase the risk of IA and may play crucial roles in regulatory functions.

Progressive optic nerve changes in cavitary optic disc anomaly: integration of copy number alteration and cis-expression quantitative trait loci to assess disease etiology

BackgroundWe performed clinical and genetic characterization of a family with cavitary optic disc anomaly (CODA), an autosomal dominant condition that causes vision loss due to adult-onset maculopathy in the majority of cases. CODA is characterized by a variably excavated optic nerve appearance such as morning glory, optic pit, atypical coloboma, and severe optic nerve cupping.MethodsFour affected and fourteen unaffected family members of a multi-generation pedigree were phenotyped by visual acuity, intraocular pressure, dilated fundus examination, fundus photography, and optical coherence tomography. Genetic analysis was performed by breakpoint polymerase chain reaction (PCR), long range PCR, and direct Sanger sequencing. The functional relevance of the copy number alteration region was assessed by in silico analysis.ResultsWe found progressive optic nerve cupping in three affected members of a family with CODA. In one individual, an optic pit developed over time from a normal optic nerve. By two independent methods, we detected a previously described intergenic triplication that segregated with disease in all adults of the family. The copy number alteration was also detected in five children with normal optic nerves. eQTL analysis demonstrated that this CNA region regulates expression of up to 4 genes in cis.ConclusionsMorning glory, optic pit and atypical coloboma are currently considered congenital anomalies of the optic nerve, but our data indicate that in CODA, the excavated optic nerve appearance may develop after birth and into adulthood. In silico analysis of the CNA, may explain why vairable expressivity is observed in CODA.