In recent years pleiotropy, the phenomenon of one genetic locus influencing several traits, has become a widely researched field in human genetics. With the increasing availability of genome-wide ...association study summary statistics, as well as the establishment of deeply phenotyped sample collections, it is now possible to systematically assess the genetic overlap between multiple traits and diseases. In addition to increasing power to detect associated variants, multi-trait methods can also aid our understanding of how different disorders are aetiologically linked by highlighting relevant biological pathways. A plethora of available tools to perform such analyses exists, each with their own advantages and limitations. In this review, we outline some of the currently available methods to conduct multi-trait analyses. First, we briefly introduce the concept of pleiotropy and outline the current landscape of pleiotropy research in human genetics; second, we describe analytical considerations and analysis methods; finally, we discuss future directions for the field.
Head and neck squamous cell carcinoma (HNSCC) remain a substantial burden to global health. Cell-free circulating tumour DNA (ctDNA) is an emerging biomarker but has not been studied sufficiently in ...HNSCC.
We conducted a single-centre prospective cohort study to investigate ctDNA in patients with p16-negative HNSCC who received curative-intent primary surgical treatment. Whole-exome sequencing was performed on formalin-fixed paraffin-embedded (FFPE) tumour tissue. We utilised RaDaR
, a highly sensitive personalised assay using deep sequencing for tumour-specific variants, to analyse serial pre- and post-operative plasma samples for evidence of minimal residual disease and recurrence.
In 17 patients analysed, personalised panels were designed to detect 34 to 52 somatic variants. Data show ctDNA detection in baseline samples taken prior to surgery in 17 of 17 patients. In post-surgery samples, ctDNA could be detected at levels as low as 0.0006% variant allele frequency. In all cases with clinical recurrence to date, ctDNA was detected prior to progression, with lead times ranging from 108 to 253 days.
This study illustrates the potential of ctDNA as a biomarker for detecting minimal residual disease and recurrence in HNSCC and demonstrates the feasibility of personalised ctDNA assays for the detection of disease prior to clinical recurrence.
The epidemiologic link between schizophrenia (SCZ) and type 2 diabetes (T2D) remains poorly understood. Here, we investigate the presence and extent of a shared genetic background between SCZ and T2D ...using genome-wide approaches. We performed a genome-wide association study (GWAS) and polygenic risk score analysis in a Greek sample collection (GOMAP) comprising three patient groups: SCZ only (n = 924), T2D only (n = 822), comorbid SCZ and T2D (n = 505); samples from two separate Greek cohorts were used as population-based controls (n = 1,125). We used genome-wide summary statistics from two large-scale GWAS of SCZ and T2D from the PGC and DIAGRAM consortia, respectively, to perform genetic overlap analyses, including a regional colocalisation test. We show for the first time that patients with comorbid SCZ and T2D have a higher genetic predisposition to both disorders compared to controls. We identify five genomic regions with evidence of colocalising SCZ and T2D signals, three of which contain known loci for both diseases. We also observe a significant excess of shared association signals between SCZ and T2D at nine out of ten investigated p value thresholds. Finally, we identify 29 genes associated with both T2D and SCZ, several of which have been implicated in biological processes relevant to these disorders. Together our results demonstrate that the observed comorbidity between SCZ and T2D is at least in part due to shared genetic mechanisms.
High-altitude adaptation in Tibetans is influenced by introgression of a 32.7-kb haplotype from the Denisovans, an extinct branch of archaic humans, lying within the endothelial PAS domain protein 1 ...(
EPAS1
), and has also been reported in Sherpa. We genotyped 19 variants in this genomic region in 1507 Eurasian individuals, including 1188 from Bhutan and Nepal residing at altitudes between 86 and 4550 m above sea level. Derived alleles for five SNPs characterizing the core Denisovan haplotype (AGGAA) were present at high frequency not only in Tibetans and Sherpa, but also among many populations from the Himalayas, showing a significant correlation with altitude (Spearman’s correlation coefficient = 0.75,
p
value 3.9 × 10
−11
). Seven East- and South-Asian 1000 Genomes Project individuals shared the Denisovan haplotype extending beyond the 32-kb region, enabling us to refine the haplotype structure and identify a candidate regulatory variant (rs370299814) that might be interacting in an additive manner with the derived G allele of rs150877473, the variant previously associated with high-altitude adaptation in Tibetans. Denisovan-derived alleles were also observed at frequencies of 3–14 % in the 1000 Genomes Project African samples. The closest African haplotype is, however, separated from the Asian high-altitude haplotype by 22 mutations whereas only three mutations, including rs150877473, separate the Asians from the Denisovan, consistent with distant shared ancestry for African and Asian haplotypes and Denisovan adaptive introgression.
Recent genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) associated with glucose levels. We tested the hypothesis here whether the cumulative effect of ...glucose raising SNPs, assessed via a score, is associated with glucose levels. A total of 1,434 participants of Greek descent from the THISEAS study and 1,160 participants form the GOMAP study were included in this analysis. We developed a genetic risk score (GRS), based on the known glucose-raising loci, in order to investigate the cumulative effect of known glucose loci on glucose levels. In the THISEAS study, the GRS score was significantly associated with increased glucose levels (mmol/L) (β ± SE: 0.024 ± 0.004, P = 8.27e-07). The effect of the genetic risk score was also significant in the GOMAP study (β ± SE: 0.011 ± 0.005, P = 0.031). In the meta-analysis of the two studies both scores were significantly associated with higher glucose levels GRS: β ± SE: 0.019 ± 0.003, P = 1.41e-09. Also, variants at the SLC30A8, PROX1, MTNR1B, ADRA2A, G6PC2, LPIN3 loci indicated nominal evidence for association with glucose levels (p < 0.05). We replicate associations of the established glucose raising variants in the Greek population and confirm directional consistency of effects (binomial sign test p = 6.96e-05). We also demonstrate that the cumulative effect of the established glucose loci yielded a significant association with increasing glucose levels.
Schizophrenia (SCZ) is associated with increased risk of type 2 diabetes (T2D). The potential diabetogenic effect of concomitant application of psychotropic treatment classes in patients with SCZ has ...not yet been evaluated. The overarching goal of the Genetic Overlap between Metabolic and Psychiatric disease (GOMAP) study is to assess the effect of pharmacological, anthropometric, lifestyle and clinical measurements, helping elucidate the mechanisms underlying the aetiology of T2D.
The GOMAP case-control study (Genetic Overlap between Metabolic and Psychiatric disease) includes hospitalized patients with SCZ, some of whom have T2D. We enrolled 1653 patients with SCZ; 611 with T2D and 1042 patients without T2D. This is the first study of SCZ and T2D comorbidity at this scale in the Greek population. We retrieved detailed information on first- and second-generation antipsychotics (FGA, SGA), antidepressants and mood stabilizers, applied as monotherapy, 2-drug combination, or as 3- or more drug combination. We assessed the effects of psychotropic medication, body mass index, duration of schizophrenia, number of hospitalizations and physical activity on risk of T2D. Using logistic regression, we calculated crude and adjusted odds ratios (OR) to identify associations between demographic factors and the psychiatric medications.
Patients with SCZ on a combination of at least three different classes of psychiatric drugs had a higher risk of T2D OR 1.81 (95% CI 1.22-2.69); p = 0.003 compared to FGA alone therapy, after adjustment for age, BMI, sex, duration of SCZ and number of hospitalizations. We did not find evidence for an association of SGA use or the combination of drugs belonging to two different classes of psychiatric medications with increased risk of T2D 1.27 (0.84-1.93), p = 0.259 and 0.98 (0.71-1.35), p = 0.885, respectively compared to FGA use.
We find an increased risk of T2D in patients with SCZ who take a combination of at least three different psychotropic medication classes compared to patients whose medication consists only of one or two classes of drugs.
Osteoarthritis is a common complex disease imposing a large public-health burden. Here, we performed a genome-wide association study for osteoarthritis, using data across 16.5 million variants from ...the UK Biobank resource. After performing replication and meta-analysis in up to 30,727 cases and 297,191 controls, we identified nine new osteoarthritis loci, in all of which the most likely causal variant was noncoding. For three loci, we detected association with biologically relevant radiographic endophenotypes, and in five signals we identified genes that were differentially expressed in degraded compared with intact articular cartilage from patients with osteoarthritis. We established causal effects on osteoarthritis for higher body mass index but not for triglyceride levels or genetic predisposition to type 2 diabetes.
Osteoarthritis (OA) is a complex disease, but its genetic aetiology remains poorly characterised. To identify novel susceptibility loci for OA, we carried out a genome-wide association study (GWAS) ...in individuals from the largest UK-based OA collections to date.
We carried out a discovery GWAS in 5414 OA individuals with knee and/or hip total joint replacement (TJR) and 9939 population-based controls. We followed-up prioritised variants in OA subjects from the interim release of the UK Biobank resource (up to 12 658 cases and 50 898 controls) and our lead finding in operated OA subjects from the full release of UK Biobank (17 894 cases and 89 470 controls). We investigated its functional implications in methylation, gene expression and proteomics data in primary chondrocytes from 12 pairs of intact and degraded cartilage samples from patients undergoing TJR.
We detect a genome-wide significant association at rs10116772 with TJR (P=3.7×10
; for allele A: OR (95% CI) 0.97 (0.96 to 0.98)), an intronic variant in
, which is expressed in cartilage. Variants in strong correlation with rs10116772 have been associated with elevated plasma glucose levels and diabetes.
We identify a novel susceptibility locus for OA that has been previously implicated in diabetes and glycaemic traits.
Cohort 1 of the phase 1B NABUCCO trial showed high pathological complete response (pCR) rates with preoperative ipilimumab plus nivolumab in stage III urothelial cancer (UC). In cohort 2, the aim was ...dose adjustment to optimize responses. Additionally, we report secondary endpoints, including efficacy and tolerability, in cohort 2 and the association of presurgical absence of circulating tumor DNA (ctDNA) in urine and plasma with clinical outcome in both cohorts. Thirty patients received two cycles of either ipilimumab 3 mg kg
plus nivolumab 1 mg kg
(cohort 2A) or ipilimumab 1 mg kg
plus nivolumab 3 mg kg
(cohort 2B), both followed by nivolumab 3 mg kg
. We observed a pCR in six (43%) patients in cohort 2A and a pCR in one (7%) patient in cohort 2B. Absence of urinary ctDNA correlated with pCR in the bladder (ypT0Nx) but not with progression-free survival (PFS). Absence of plasma ctDNA correlated with pCR (odds ratio: 45.0; 95% confidence interval (CI): 4.9-416.5) and PFS (hazard ratio: 10.4; 95% CI: 2.9-37.5). Our data suggest that high-dose ipilimumab plus nivolumab is required in stage III UC and that absence of ctDNA in plasma can predict PFS. ClinicalTrials.gov registration: NCT03387761 .
Osteoarthritis (OA) is a common complex disease with high public health burden and no curative therapy. High bone mineral density (BMD) is associated with an increased risk of developing OA, ...suggesting a shared underlying biology. Here, we performed the first systematic overlap analysis of OA and BMD on a genome wide scale. We used summary statistics from the GEFOS consortium for lumbar spine (n = 31,800) and femoral neck (n = 32,961) BMD, and from the arcOGEN consortium for three OA phenotypes (hip, ncases=3,498; knee, ncases=3,266; hip and/or knee, ncases=7,410; ncontrols=11,009). Performing LD score regression we found a significant genetic correlation between the combined OA phenotype (hip and/or knee) and lumbar spine BMD (rg=0.18, P = 2.23 × 10-2), which may be driven by the presence of spinal osteophytes. We identified 143 variants with evidence for cross-phenotype association which we took forward for replication in independent large-scale OA datasets, and subsequent meta-analysis with arcOGEN for a total sample size of up to 23,425 cases and 236,814 controls. We found robustly replicating evidence for association with OA at rs12901071 (OR 1.08 95% CI 1.05-1.11, Pmeta=3.12 × 10-10), an intronic variant in the SMAD3 gene, which is known to play a role in bone remodeling and cartilage maintenance. We were able to confirm expression of SMAD3 in intact and degraded cartilage of the knee and hip. Our findings provide the first systematic evaluation of pleiotropy between OA and BMD, highlight genes with biological relevance to both traits, and establish a robust new OA genetic risk locus at SMAD3.
PDF
