SCAN database: facilitating integrative analyses of cytosine modification and expression QTL

Endophenotype 2.0: updated definitions and criteria for endophenotypes of psychiatric disorders, incorporating new technologies and findings

Recent genetic studies have linked numerous loci to psychiatric disorders. However, the biological pathways that connect these genetic associations to psychiatric disorders' specific pathophysiological processes are largely unclear. Endophenotypes, first defined over five decades ago, are heritable traits, independent of disease state that are associated with a disease, encompassing a broad range of neurophysiological, biochemical, endocrinological, neuroanatomical, cognitive, and neuropsychological characteristics. Considering the advancements in genetics and genomics over recent decades, we propose a revised definition of endophenotypes as 'genetically influenced phenotypes linked to disease or treatment characteristics and their related events.' We also updated endophenotype criteria to include (1) reliable measurement, (2) association with the disease or its related events, and (3) genetic mediation. 'Genetic mediation' is introduced to differentiate between causality and pleiotropic effects and allows non-linear relationships. Furthermore, this updated Endophenotype 2.0 framework expands to encompass genetically regulated responses to disease-related factors, including environmental risks, illness progression, treatment responses, and resilience phenotypes, which may be state-dependent. This broadened definition paves the way for developing new endophenotypes crucial for genetic analyses in psychiatric disorders. Integrating genetics, genomics, and diverse endophenotypes into multi-dimensional mechanistic models is vital for advancing our understanding of psychiatric disorders. Crucially, elucidating the biological underpinnings of endophenotypes will enhance our grasp of psychiatric genetics, thereby improving disease risk prediction and treatment approaches.

Genetic Screening of Targeted Region on the Chromosome 22q11.2 in Patients with Microtia and Congenital Heart Defect

Microtia is a congenital malformation characterized by a small, abnormally shaped auricle (pinna) ranging in severity. Congenital heart defect (CHD) is one of the comorbid anomalies with microtia. However, the genetic basis of the co-existence of microtia and CHD remains unclear. Copy number variations (CNVs) of 22q11.2 contribute significantly to microtia and CHD, respectively, thus suggesting a possible shared genetic cause embedded in this genomic region. In this study, 19 sporadic patients with microtia and CHD, as well as a nuclear family, were enrolled for genetic screening of single nucleotide variations (SNVs) and CNVs in 22q11.2 by target capture sequencing. We detected a total of 105 potential deleterious variations, which were enriched in ear- or heart-development-related genes, including TBX1 and DGCR8. The gene burden analysis also suggested that these genes carry more deleterious mutations in the patients, as well as several other genes associated with cardiac development, such as CLTCL1. Additionally, a microduplication harboring SUSD2 was validated in an independent cohort. This study provides new insights into the underlying mechanisms for the comorbidity of microtia and CHD focusing on chromosome 22q11.2, and suggests that a combination of genetic variations, including SNVs and CNVs, may play a crucial role instead of single gene mutation.

The Type 2 Diabetes Knowledge Portal: An open access genetic resource dedicated to type 2 diabetes and related traits

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results.

Bevacizumab-induced hypertension and proteinuria: a genome-wide study of more than 1000 patients

BACKGROUND

Hypertension and proteinuria are common bevacizumab-induced toxicities. No validated biomarkers are available for identifying patients at risk of these toxicities.

METHODS

A genome-wide association study (GWAS) meta-analysis was performed in 1039 bevacizumab-treated patients of European ancestry in four clinical trials (CALGB 40502, 40503, 80303, 90401). Grade ≥2 hypertension and proteinuria were recorded (CTCAE v.3.0). Single-nucleotide polymorphism (SNP)-toxicity associations were determined using a cause-specific Cox model adjusting for age and sex.

RESULTS

The most significant SNP associated with hypertension with concordant effect in three out of the four studies (p-value <0.05 for each study) was rs6770663 (A > G) in KCNAB1, with the G allele increasing the risk of hypertension (p-value = 4.16 × 10-6). The effect of the G allele was replicated in ECOG-ACRIN E5103 in 582 patients (p-value = 0.005). The meta-analysis of all five studies for rs6770663 led to p-value = 7.73 × 10-8, close to genome-wide significance. The most significant SNP associated with proteinuria was rs339947 (C > A, between DNAH5 and TRIO), with the A allele increasing the risk of proteinuria (p-value = 1.58 × 10-7).

CONCLUSIONS

The results from the largest study of bevacizumab toxicity provide new markers of drug safety for further evaluations. SNP in KCNAB1 validated in an independent dataset provides evidence toward its clinical applicability to predict bevacizumab-induced hypertension. ClinicalTrials.gov Identifier: NCT00785291 (CALGB 40502); NCT00601900 (CALGB 40503); NCT00088894 (CALGB 80303) and NCT00110214 (CALGB 90401).

PIK3R5 genetic predictors of hypertension induced by VEGF-pathway inhibitors

No biomarkers are available to predict patients at risk of developing hypertension induced by VEGF-pathway inhibitors. This study aimed to identify predictive biomarkers of hypertension induced by these drugs using a discovery-replication approach. The discovery set included 140 sorafenib-treated patients (TARGET study) genotyped for 973 SNPs in 56 genes. The most statistically significant SNPs associated with grade ≥2 hypertension were tested for association with grade ≥2 hypertension in the replication set of a GWAS of 1039 bevacizumab-treated patients from four clinical trials (CALGB/Alliance). In the discovery set, rs444904 (G > A) in PIK3R5 was associated with an increased risk of sorafenib-induced hypertension (p = 0.006, OR = 3.88 95% CI 1.54-9.81). In the replication set, rs427554 (G > A) in PIK3R5 (in complete linkage disequilibrium with rs444904) was associated with an increased risk of bevacizumab-induced hypertension (p = 0.008, OR = 1.39, 95% CI 1.09-1.78). This study identified a predictive marker of drug-induced hypertension that should be evaluated for other VEGF-pathway inhibitors.ClinicalTrials.gov Identifier:NCT00073307 (TARGET).

Similar burden of pathogenic coding variants in exceptionally long-lived individuals and individuals without exceptional longevity

Centenarians (exceptionally long‐lived individuals—ELLI) are a unique segment of the population, exhibiting long human lifespan and healthspan, despite generally practicing similar lifestyle habits as their peers. We tested disease‐associated mutation burden in ELLI genomes by determining the burden of pathogenic variants reported in the ClinVar and HGMD databases using data from whole exome sequencing (WES) conducted in a cohort of ELLI, their offspring, and control individuals without antecedents of familial longevity (n = 1879), all descendent from the founder population of Ashkenazi Jews. The burden of pathogenic variants did not differ between the three groups. Additional analyses of variants subtypes and variant effect predictor (VEP) biotype frequencies did not reveal a decrease of pathogenic or loss‐of‐function (LoF) variants in ELLI and offspring compared to the control group. Case–control pathogenic variants enrichment analyses conducted in ELLI and controls also did not identify significant differences in any of the variants between the groups and polygenic risk scores failed to provide a predictive model. Interestingly, cancer and Alzheimer's disease‐associated variants were significantly depleted in ELLI compared to controls, suggesting slower accumulation of mutation. That said, polygenic risk score analysis failed to find any predictive variants among the functional variants tested. The high similarity in the burden of pathogenic variation between ELLI and individuals without familial longevity supports the notion that extension of lifespan and healthspan in ELLI is not a consequence of pathogenic variant depletion but rather a result of other genomic, epigenomic, or potentially nongenomic properties.

An initial genetic analysis of gemcitabine-induced high-grade neutropenia in pancreatic cancer patients in CALGB 80303 (Alliance)

OBJECTIVES

One of the standard of care regimens for advanced pancreatic cancer is gemcitabine-based chemotherapy. The efficacy of gemcitabine is limited by dose-limiting hematologic toxicities especially neutropenia. Uncovering the variability of these toxicities attributed to germline DNA variation is of great importance.

PATIENTS AND METHODS

CALGB 80303 was a randomized study in advanced pancreatic cancer patients treated with gemcitabine with or without bevacizumab. The study protocol included genotyping of genes of gemcitabine disposition (CDA, DCTD, SLC29A1, SLC28A1, and SLC29A2), as well as a genome-wide analysis. The clinical phenotype was time to early high-grade neutropenia event accounting for progression or death or other treatment-terminating adverse events as competing for informative events. The inference was carried out on the basis of the association between genotype and cause-specific hazard of a neutropenic event.

RESULTS

The primary analyses were carried out on the basis of 294 genetically estimated European pancreatic cancer patients. For CDA rs2072671 (A>C), AC and CC patients had a lower risk of neutropenia than AA patients (P=0.01, hazard ratio: 0.61, 95% confidence interval: 0.41-0.89). For SLC28A1 rs3825876 (G>A), AA patients have a higher risk of neutropenia than GA and GG patients (P=0.02, hazard ratio: 1.51, 95% confidence interval: 1.06-2.16). CDA rs2072671 was associated with increased mRNA expression in whole blood in three studies (P=2.7e-14, 6.61e-62, and 9.70e-65). In the genome-wide analysis, variants in TGFB2 were among the top hits (lowest P=1.62e-06) but had no effect in luciferase assays.

CONCLUSION

This is the first genetic analysis of gemcitabine-induced neutropenia using a competing risk model in a prospective randomized clinical study has proposed a potentially novel mechanism of the protective effect of the CDA rs2072671 variant. Further confirmation is needed.

Cancer-specific expression quantitative loci are affected by expression dysregulation

Expression quantitative trait loci (eQTLs) have been touted as the missing piece that can bridge the gap between genetic variants and phenotypes. Over the past decade, we have witnessed a sharp rise of effort in the identification and application of eQTLs. The successful application of eQTLs relies heavily on their reproducibility. The current eQTL databases such as Genotype-Tissue Expression (GTEx) were populated primarily with eQTLs deriving from germline single nucleotide polymorphisms and normal tissue gene expression. The novel scenarios that employ eQTL models for prediction purposes often involve disease phenotypes characterized by altered gene expressions. To evaluate eQTL reproducibility across diverse data sources and the effect of disease-specific gene expression alteration on eQTL identification, we conducted an eQTL study using 5178 samples from The Cancer Genome Atlas (TCGA). We found that the reproducibility of eQTLs between normal and tumor tissues was low in terms of the number of shared eQTLs. However, among the shared eQTLs, the effect directions were generally concordant. This suggests that the source of the gene expression (normal or tumor tissue) has a strong effect on the detectable eQTLs and the effect direction of the eQTLs. Additional analyses demonstrated good directional concordance of eQTLs between GTEx and TCGA. Furthermore, we found that multi-tissue eQTLs may exert opposite effects across multiple tissue types. In summary, our results suggest that eQTL prediction models need to carefully address tissue and disease dependency of eQTLs. Tissue-disease-specific eQTL databases can afford more accurate prediction models for future studies.

Exploratory analysis of genetic variants influencing molecular traits in cerebral cortex of suicide completers

Genetic factors have been implicated in suicidal behavior. It has been suggested that one of the roles of genetic factors in suicide could be represented by the effect of genetic variants on gene expression regulation. Alteration in the expression of genes participating in multiple biological systems in the suicidal brain has been demonstrated, so it is imperative to identify genetic variants that could influence gene expression or its regulatory mechanisms. In this study, we integrated DNA methylation, gene expression, and genotype data from the prefrontal cortex of suicides to identify genetic variants that could be factors in the regulation of gene expression, generally called quantitative trait locus (xQTLs). We identify 6,224 methylation quantitative trait loci and 2,239 expression quantitative trait loci (eQTLs) in the prefrontal cortex of suicide completers. The xQTLs identified influence the expression of genes involved in neurodevelopment and cell organization. Two of the eQTLs identified (rs8065311 and rs1019238) were previously associated with cannabis dependence, highlighting a candidate genetic variant for the increased suicide risk in subjects with substance use disorders. Our findings suggest that genetic variants may regulate gene expression in the prefrontal cortex of suicides through the modulation of promoter and enhancer activity, and to a lesser extent, binding transcription factors.

Genetic architecture of subcortical brain structures in 38,851 individuals

Subcortical brain structures are integral to motion, consciousness, emotions and learning. We identified common genetic variation related to the volumes of the nucleus accumbens, amygdala, brainstem, caudate nucleus, globus pallidus, putamen and thalamus, using genome-wide association analyses in almost 40,000 individuals from CHARGE, ENIGMA and UK Biobank. We show that variability in subcortical volumes is heritable, and identify 48 significantly associated loci (40 novel at the time of analysis). Annotation of these loci by utilizing gene expression, methylation and neuropathological data identified 199 genes putatively implicated in neurodevelopment, synaptic signaling, axonal transport, apoptosis, inflammation/infection and susceptibility to neurological disorders. This set of genes is significantly enriched for Drosophila orthologs associated with neurodevelopmental phenotypes, suggesting evolutionarily conserved mechanisms. Our findings uncover novel biology and potential drug targets underlying brain development and disease.

PancanQTL: systematic identification of cis-eQTLs and trans-eQTLs in 33 cancer types

Expression quantitative trait locus (eQTL) analysis, which links variations in gene expression to genotypes, is essential to understanding gene regulation and to interpreting disease-associated loci. Currently identified eQTLs are mainly in samples of blood and other normal tissues. However, no database comprehensively provides eQTLs in large number of cancer samples. Using the genotype and expression data of 9196 tumor samples in 33 cancer types from The Cancer Genome Atlas (TCGA), we identified 5 606 570 eQTL-gene pairs in the cis-eQTL analysis and 231 210 eQTL-gene pairs in the trans-eQTL analysis. We further performed survival analysis and identified 22 212 eQTLs associated with patient overall survival. Furthermore, we linked the eQTLs to genome-wide association studies (GWAS) data and identified 337 131 eQTLs that overlap with existing GWAS loci. We developed PancanQTL, a user-friendly database (http://bioinfo.life.hust.edu.cn/PancanQTL/), to store cis-eQTLs, trans-eQTLs, survival-associated eQTLs and GWAS-related eQTLs to enable searching, browsing and downloading. PancanQTL could help the research community understand the effects of inherited variants in tumorigenesis and development.

The framework for population epigenetic study

At present, understanding of DNA methylation at the population level is still limited. Here, we first extended the classical framework of population genetics, such as single nucleotide polymorphism allele frequency, linkage disequilibrium (LD), LD block and haplotype, to epigenetics. Then, as an example, we compared the DNA methylation disequilibrium (MD) maps between HapMap CEU (Caucasian residents of European ancestry from Utah) population and YRI (Yoruba people from Ibadan) population (lymphoblastoid cell lines). We analyzed the differences and similarities between CEU and YRI from the following aspects: SMP (single methylation polymorphism) allele frequency, SMP allele association, MD, MD block and methylation haplotype (meplotype) frequency. The results showed that CEU and YRI had similar distribution of SMP allele frequency, and shared many MD block region. We believe that the framework of population genetics can be used in the population epigenetics. The population epigenetic framework also has potential prospects in the study of complex diseases, such as epigenome-wide association study.

Genome-wide association study of maternal genetic effects and parent-of-origin effects on food allergy

Previous genetic studies of food allergy (FA) have mainly focused on inherited genotypic effects. The role of parental genotypic effects remains largely unexplored. Leveraging existing genome-wide association study (GWAS) data generated from the Chicago Food Allergy Study, we examined maternal genotypic and parent-of-origin (PO) effects using multinomial likelihood ratio tests in 588 complete and incomplete Caucasian FA trios. We identified 1 single nucleotide polymorphism with significant (P < 5×10) maternal effect on any FA (rs4235235), which is located in a noncoding RNA (LOC101927947) with unknown function. We also identified 3 suggestive (P < 5×10) loci with maternal genetic effects: 1 for any FA (rs976078, in a gene desert region on 13q31.1) and 2 for egg allergy (rs1343795 and rs4572450, in the ZNF652 gene, where genetic variants have been associated with atopic dermatitis). Three suggestive loci with PO effect were observed: 1 for peanut allergy (rs4896888 in the ADGB gene) and 2 for any FA in boys only (rs1036504 and rs2917750 in the IQCE gene). Findings from this family-based GWAS of FA provided some preliminary evidence on maternal genotypic or PO effects on FA. Additional family-based studies are needed to confirm our findings and gain new insight into maternal and paternal genetic contribution to FA.

Genetic-epigenetic interactions in cis: a major focus in the post-GWAS era

Studies on genetic-epigenetic interactions, including the mapping of methylation quantitative trait loci (mQTLs) and haplotype-dependent allele-specific DNA methylation (hap-ASM), have become a major focus in the post-genome-wide-association-study (GWAS) era. Such maps can nominate regulatory sequence variants that underlie GWAS signals for common diseases, ranging from neuropsychiatric disorders to cancers. Conversely, mQTLs need to be filtered out when searching for non-genetic effects in epigenome-wide association studies (EWAS). Sequence variants in CCCTC-binding factor (CTCF) and transcription factor binding sites have been mechanistically linked to mQTLs and hap-ASM. Identifying these sites can point to disease-associated transcriptional pathways, with implications for targeted treatment and prevention.

Epigenetic contribution of the myosin light chain kinase gene to the risk for acute respiratory distress syndrome

Acute respiratory distress syndrome (ARDS) is a devastating clinical syndrome with a considerable case fatality rate (∼30%-40%). Health disparities exist with African descent (AD) subjects exhibiting greater mortality than European descent (ED) individuals. Myosin light chain kinase is encoded by MYLK, whose genetic variants are implicated in ARDS pathogenesis and may influence ARDS mortality. As baseline population-specific epigenetic changes, that is, cytosine modifications, have been observed between AD and ED individuals, epigenetic variations in MYLK may provide insights into ARDS disparities. We compared methylation levels of MYLK cytosine-guanine dinucleotides (CpGs) between ARDS patients and intensive care unit (ICU) controls overall and by ethnicity in a nested case-control study of 39 ARDS cases and 75 non-ARDS ICU controls. Two MYLK CpG sites (cg03892735 and cg23344121) were differentially modified between ARDS subjects and controls (P < 0.05; q < 0.25) in a logistic regression model, where no effect modification by ethnicity or age was found. One CpG site was associated with ARDS in patients aged <58 years, cg19611163 (intron 19, 20). Two CpG sites were associated with ARDS in EDs only, gene body CpG (cg01894985, intron 2, 3) and CpG (cg16212219, intron 31, 32), with higher modification levels exhibited in ARDS subjects than controls. Cis-acting modified cytosine quantitative trait loci (mQTL) were identified using linear regression between local genetic variants and modification levels for 2 ARDS-associated CpGs (cg23344121 and cg16212219). In summary, these ARDS-associated MYLK CpGs with effect modification by ethnicity and local mQTL suggest that MYLK epigenetic variation and local genetic background may contribute to health disparities observed in ARDS.

Chronic Periodontitis Genome-wide Association Study in the Hispanic Community Health Study / Study of Latinos

Chronic periodontitis (CP) has a genetic component, particularly its severe forms. Evidence from genome-wide association studies (GWASs) has highlighted several potential novel loci. Here, the authors report the first GWAS of CP among a large community-based sample of Hispanics/Latinos. The authors interrogated a quantitative trait of CP (mean interproximal clinical attachment level determined by full-mouth periodontal examinations) among 10,935 adult participants (mean age: 45 y, range: 18 to 76 y) from the Hispanic Community Health Study / Study of Latinos. Genotyping was done with a custom Illumina Omni2.5M array, and imputation to approximately 20 million single-nucleotide polymorphisms was based on the 1000 Genomes Project phase 1 reference panel. Analyses were based on linear mixed models adjusting for sex, age, study design features, ancestry, and kinship and employed a conventional P < 5 × 10^-8 statistical significance threshold. The authors identified a genome-wide significant association signal in the 1q42.2 locus ( TSNAX-DISC1 noncoding RNA, lead single-nucleotide polymorphism: rs149133391, minor allele [C] frequency = 0.01, P = 7.9 × 10^-9) and 4 more loci with suggestive evidence of association ( P < 5 × 10^-6): 1q22 (rs13373934), 5p15.33 (rs186066047), 6p22.3 (rs10456847), and 11p15.1 (rs75715012). We tested these loci for replication in independent samples of European-American ( n = 4,402) and African-American ( n = 908) participants of the Atherosclerosis Risk in Communities study. There was no replication among the European Americans; however, the TSNAX-DISC1 locus replicated in the African-American sample (rs149133391, minor allele frequency = 0.02, P = 9.1 × 10^-3), while the 1q22 locus was directionally concordant and nominally significant (rs13373934, P = 4.0 × 10^-2). This discovery GWAS of interproximal clinical attachment level-a measure of lifetime periodontal tissue destruction-was conducted in a large, community-based sample of Hispanic/Latinos. It identified a genome-wide significant locus that was independently replicated in an African-American population. Identifying this genetic marker offers direction for interrogation in subsequent genomic and experimental studies of CP.

Linking short tandem repeat polymorphisms with cytosine modifications in human lymphoblastoid cell lines

Inter-individual variation in cytosine modifications has been linked to complex traits in humans. Cytosine modification variation is partially controlled by single nucleotide polymorphisms (SNPs), known as modified cytosine quantitative trait loci (mQTL). However, little is known about the role of short tandem repeat polymorphisms (STRPs), a class of structural genetic variants, in regulating cytosine modifications. Utilizing the published data on the International HapMap Project lymphoblastoid cell lines (LCLs), we assessed the relationships between 721 STRPs and the modification levels of 283,540 autosomal CpG sites. Our findings suggest that, in contrast to the predominant cis-acting mode for SNP-based mQTL, STRPs are associated with cytosine modification levels in both cis-acting (local) and trans-acting (distant) modes. In local scans within the ±1 Mb windows of target CpGs, 21, 9, and 21 cis-acting STRP-based mQTL were detected in CEU (Caucasian residents from Utah, USA), YRI (Yoruba people from Ibadan, Nigeria), and the combined samples, respectively. In contrast, 139,420, 76,817, and 121,866 trans-acting STRP-based mQTL were identified in CEU, YRI, and the combined samples, respectively. A substantial proportion of CpG sites detected with local STRP-based mQTL were not associated with SNP-based mQTL, suggesting that STRPs represent an independent class of mQTL. Functionally, genetic variants neighboring CpG-associated STRPs are enriched with genome-wide association study (GWAS) loci for a variety of complex traits and diseases, including cancers, based on the National Human Genome Research Institute (NHGRI) GWAS Catalog. Therefore, elucidating these STRP-based mQTL in addition to SNP-based mQTL can provide novel insights into the genetic architectures of complex traits.

The Human Phenotype Ontology: Semantic Unification of Common and Rare Disease

The Human Phenotype Ontology (HPO) is widely used in the rare disease community for differential diagnostics, phenotype-driven analysis of next-generation sequence-variation data, and translational research, but a comparable resource has not been available for common disease. Here, we have developed a concept-recognition procedure that analyzes the frequencies of HPO disease annotations as identified in over five million PubMed abstracts by employing an iterative procedure to optimize precision and recall of the identified terms. We derived disease models for 3,145 common human diseases comprising a total of 132,006 HPO annotations. The HPO now comprises over 250,000 phenotypic annotations for over 10,000 rare and common diseases and can be used for examining the phenotypic overlap among common diseases that share risk alleles, as well as between Mendelian diseases and common diseases linked by genomic location. The annotations, as well as the HPO itself, are freely available.