A comparison of methods for detecting DNA methylation from long-read sequencing of human genomes.

BACKGROUND: Long-read sequencing can enable the detection of base modifications, such as CpG methylation, in single molecules of DNA. The most commonly used methods for long-read sequencing are nanopore developed by Oxford Nanopore Technologies (ONT) and single molecule real-time (SMRT) sequencing developed by Pacific Bioscience (PacBio). In this study, we systematically compare the performance of CpG methylation detection from long-read sequencing. RESULTS: We demonstrate that CpG methylation detection from 7179 nanopore-sequenced DNA samples is highly accurate and consistent with 132 oxidative bisulfite-sequenced (oxBS) samples, isolated from the same blood draws. We introduce quality filters for CpGs that further enhance the accuracy of CpG methylation detection from nanopore-sequenced DNA, while removing at most 30% of CpGs. We evaluate the per-site performance of CpG methylation detection across different genomic features and CpG methylation rates and demonstrate how the latest R10.4 flowcell chemistry and base-calling algorithms improve methylation detection from nanopore sequencing. Additionally, we show how the methylation detection of 50 SMRT-sequenced genomes compares to nanopore sequencing and oxBS. CONCLUSIONS: This study provides the first systematic comparison of CpG methylation detection tools for long-read sequencing methods. We compare two commonly used computational methods for the detection of CpG methylation in a large number of nanopore genomes, including samples sequenced using the latest R10.4 nanopore flowcell chemistry and 50 SMRT sequenced samples. We provide insights into the strengths and limitations of each sequencing method as well as recommendations for standardization and evaluation of tools designed for genome-scale modified base detection using long-read sequencing.

Predicting the presence of coronary plaques featuring high-risk characteristics using polygenic risk scores and targeted proteomics in patients with suspected coronary artery disease.

BACKGROUND: The presence of coronary plaques with high-risk characteristics is strongly associated with adverse cardiac events beyond the identification of coronary stenosis. Testing by coronary computed tomography angiography (CCTA) enables the identification of high-risk plaques (HRP). Referral for CCTA is presently based on pre-test probability estimates including clinical risk factors (CRFs); however, proteomics and/or genetic information could potentially improve patient selection for CCTA and, hence, identification of HRP. We aimed to (1) identify proteomic and genetic features associated with HRP presence and (2) investigate the effect of combining CRFs, proteomics, and genetics to predict HRP presence. METHODS: Consecutive chest pain patients (n = 1462) undergoing CCTA to diagnose obstructive coronary artery disease (CAD) were included. Coronary plaques were assessed using a semi-automatic plaque analysis tool. Measurements of 368 circulating proteins were obtained with targeted Olink panels, and DNA genotyping was performed in all patients. Imputed genetic variants were used to compute a multi-trait multi-ancestry genome-wide polygenic score (GPSMult). HRP presence was defined as plaques with two or more high-risk characteristics (low attenuation, spotty calcification, positive remodeling, and napkin ring sign). Prediction of HRP presence was performed using the glmnet algorithm with repeated fivefold cross-validation, using CRFs, proteomics, and GPSMult as input features. RESULTS: HRPs were detected in 165 (11%) patients, and 15 input features were associated with HRP presence. Prediction of HRP presence based on CRFs yielded a mean area under the receiver operating curve (AUC) ± standard error of 73.2 ± 0.1, versus 69.0 ± 0.1 for proteomics and 60.1 ± 0.1 for GPSMult. Combining CRFs with GPSMult increased prediction accuracy (AUC 74.8 ± 0.1 (P = 0.004)), while the inclusion of proteomics provided no significant improvement to either the CRF (AUC 73.2 ± 0.1, P = 1.00) or the CRF + GPSMult (AUC 74.6 ± 0.1, P = 1.00) models, respectively. CONCLUSIONS: In patients with suspected CAD, incorporating genetic data with either clinical or proteomic data improves the prediction of high-risk plaque presence. TRIAL REGISTRATION: https://clinicaltrials.gov/ct2/show/NCT02264717 (September 2014).

Obesity Variants in the GIPR Gene Are not Associated With Risk of Fracture or Bone Mineral Density.

CONTEXT: It is not clear if antagonizing the GIP (glucose-dependent insulinotropic polypeptide) receptor (GIPR) for treatment of obesity is likely to increase the risk of fractures, or to lower bone mineral density (BMD) beyond what is expected with rapid weight loss. OBJECTIVE: The objective of this study was to investigate the risk of fracture and BMD of sequence variants in GIPR that reduce the activity of the GIP receptor and have been associated with reduced body mass index (BMI). METHODS: We analyzed the association of 3 missense variants in GIPR, a common variant, rs1800437 (p.Glu354Gln), and 2 rare variants, rs139215588 (p.Arg190Gln) and rs143430880 (p.Glu288Gly), as well as a burden of predicted loss-of-function (LoF) variants with risk of fracture and with BMD in a large meta-analysis of up to 1.2 million participants. We analyzed associations with fractures at different skeletal sites in the general population: any fractures, hip fractures, vertebral fractures and forearm fractures, and specifically nonvertebral and osteoporotic fractures in postmenopausal women. We also evaluated associations with BMD at the lumbar spine, femoral neck, and total body measured with dual-energy x-ray absorptiometry (DXA), and with BMD estimated from heel ultrasound (eBMD). RESULTS: None of the 3 missense variants in GIPR was significantly associated with increased risk of fractures or with lower BMD. Burden of LoF variants in GIPR was not associated with fractures or with BMD measured with clinically validated DXA, but was associated with eBMD. CONCLUSION: Missense variants in GIPR, or burden of LoF variants in the gene, are not associated with risk of fractures or with lower BMD.

Variant in the synaptonemal complex protein SYCE2 associates with pregnancy loss through effect on recombination.

Two-thirds of all human conceptions are lost, in most cases before clinical detection. The lack of detailed understanding of the causes of pregnancy losses constrains focused counseling for future pregnancies. We have previously shown that a missense variant in synaptonemal complex central element protein 2 (SYCE2), in a key residue for the assembly of the synaptonemal complex backbone, associates with recombination traits. Here we show that it also increases risk of pregnancy loss in a genome-wide association analysis on 114,761 women with reported pregnancy loss. We further show that the variant associates with more random placement of crossovers and lower recombination rate in longer chromosomes but higher in the shorter ones. These results support the hypothesis that some pregnancy losses are due to failures in recombination. They further demonstrate that variants with a substantial effect on the quality of recombination can be maintained in the population.

Loss-of-function variants in ITSN1 confer high risk of Parkinson's disease.

Parkinson's disease (PD) is a debilitating neurodegenerative disorder and its rising global incidence highlights the need for the identification of modifiable risk factors. In a gene-based burden test of rare variants (8647 PD cases and 777,693 controls) we discovered a novel association between loss-of-function variants in ITSN1 and PD. This association was further supported with burden data from the Neurodegenerative Disease Knowledge Portal and the Accelerating Medicines Partnership Parkinson's Disease Knowledge Platform. Our findings show that Rho GTPases and disruptions in synaptic vesicle transport may be involved in the pathogenesis of PD, pointing to the possibility of novel therapeutic approaches.

GWAS meta-analysis reveals key risk loci in essential tremor pathogenesis.

Essential tremor (ET) is a prevalent neurological disorder with a largely unknown underlying biology. In this genome-wide association study meta-analysis, comprising 16,480 ET cases and 1,936,173 controls from seven datasets, we identify 12 sequence variants at 11 loci. Evaluating mRNA expression, splicing, plasma protein levels, and coding effects, we highlight seven putative causal genes at these loci, including CA3 and CPLX1. CA3 encodes Carbonic Anhydrase III and carbonic anhydrase inhibitors have been shown to decrease tremors. CPLX1, encoding Complexin-1, regulates neurotransmitter release. Through gene-set enrichment analysis, we identify a significant association with specific cell types, including dopaminergic and GABAergic neurons, as well as biological processes like Rho GTPase signaling. Genetic correlation analyses reveals a positive association between ET and Parkinson's disease, depression, and anxiety-related phenotypes. This research uncovers risk loci, enhancing our knowledge of the complex genetics of this common but poorly understood disorder, and highlights CA3 and CPLX1 as potential therapeutic targets.

Genome-wide association meta-analysis identifies five loci associated with postpartum hemorrhage.

Bleeding in early pregnancy and postpartum hemorrhage (PPH) bear substantial risks, with the former closely associated with pregnancy loss and the latter being the foremost cause of maternal death, underscoring the severe impact on maternal-fetal health. We identified five genetic loci linked to PPH in a meta-analysis. Functional annotation analysis indicated candidate genes HAND2, TBX3 and RAP2C/FRMD7 at three loci and showed that at each locus, associated variants were located within binding sites for progesterone receptors. There were strong genetic correlations with birth weight, gestational duration and uterine fibroids. Bleeding in early pregnancy yielded no genome-wide association signals but showed strong genetic correlation with various human traits, suggesting a potentially complex, polygenic etiology. Our results suggest that PPH is related to progesterone signaling dysregulation, whereas early bleeding is a complex trait associated with underlying health and possibly socioeconomic status and may include genetic factors that have not yet been identified.

The correlation between CpG methylation and gene expression is driven by sequence variants.

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Start codon variant in LAG3 is associated with decreased LAG-3 expression and increased risk of autoimmune thyroid disease.

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10-16) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10-3) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns.

X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements.

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

Biallelic variants in POPDC2 cause a novel autosomal recessive syndrome presenting with cardiac conduction defects and variable hypertrophic cardiomyopathy.

POPDC2 encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss of Popdc2 in mice results in sinus pauses and bradycardia and morpholino knockdown of popdc2 in zebrafish results in atrioventricular (AV) block. We identified bi-allelic variants in POPDC2 in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identified POPDC2 variants are predicted to diminish the ability of POPDC2 to bind cAMP. In in vitro electrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1-POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants in POPDC1 and POPDC2 and presence of sinus node disease in POPDC2, but not in POPDC1 related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence for POPDC2 as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functional POPDC2 display sinus and AV node dysfunction.

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.