A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Identification of ALK in Thinness.

There is considerable inter-individual variability in susceptibility to weight gain despite an equally obesogenic environment in large parts of the world. Whereas many studies have focused on identifying the genetic susceptibility to obesity, we performed a GWAS on metabolically healthy thin individuals (lowest 6th percentile of the population-wide BMI spectrum) in a uniquely phenotyped Estonian cohort. We discovered anaplastic lymphoma kinase (ALK) as a candidate thinness gene. In Drosophila, RNAi mediated knockdown of Alk led to decreased triglyceride levels. In mice, genetic deletion of Alk resulted in thin animals with marked resistance to diet- and leptin-mutation-induced obesity. Mechanistically, we found that ALK expression in hypothalamic neurons controls energy expenditure via sympathetic control of adipose tissue lipolysis. Our genetic and mechanistic experiments identify ALK as a thinness gene, which is involved in the resistance to weight gain.

Genetic drivers and cellular selection of female mosaic X chromosome loss.

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Evaluation of polygenic scoring methods in five biobanks shows larger variation between biobanks than methods and finds benefits of ensemble learning.

Methods of estimating polygenic scores (PGSs) from genome-wide association studies are increasingly utilized. However, independent method evaluation is lacking, and method comparisons are often limited. Here, we evaluate polygenic scores derived via seven methods in five biobank studies (totaling about 1.2 million participants) across 16 diseases and quantitative traits, building on a reference-standardized framework. We conducted meta-analyses to quantify the effects of method choice, hyperparameter tuning, method ensembling, and the target biobank on PGS performance. We found that no single method consistently outperformed all others. PGS effect sizes were more variable between biobanks than between methods within biobanks when methods were well tuned. Differences between methods were largest for the two investigated autoimmune diseases, seropositive rheumatoid arthritis and type 1 diabetes. For most methods, cross-validation was more reliable for tuning hyperparameters than automatic tuning (without the use of target data). For a given target phenotype, elastic net models combining PGS across methods (ensemble PGS) tuned in the UK Biobank provided consistent, high, and cross-biobank transferable performance, increasing PGS effect sizes (ß coefficients) by a median of 5.0% relative to LDpred2 and MegaPRS (the two best-performing single methods when tuned with cross-validation). Our interactively browsable online-results and open-source workflow prspipe provide a rich resource and reference for the analysis of polygenic scoring methods across biobanks.

Genetic insights into the age-specific biological mechanisms governing human ovarian aging.

There is currently little evidence that the genetic basis of human phenotype varies significantly across the lifespan. However, time-to-event phenotypes are understudied and can be thought of as reflecting an underlying hazard, which is unlikely to be constant through life when values take a broad range. Here, we find that 74% of 245 genome-wide significant genetic associations with age at natural menopause (ANM) in the UK Biobank show a form of age-specific effect. Nineteen of these replicated discoveries are identified only by our modeling framework, which determines the time dependency of DNA-variant age-at-onset associations without a significant multiple-testing burden. Across the range of early to late menopause, we find evidence for significantly different underlying biological pathways, changes in the signs of genetic correlations of ANM to health indicators and outcomes, and differences in inferred causal relationships. We find that DNA damage response processes only act to shape ovarian reserve and depletion for women of early ANM. Genetically mediated delays in ANM were associated with increased relative risk of breast cancer and leiomyoma at all ages and with high cholesterol and heart failure for late-ANM women. These findings suggest that a better understanding of the age dependency of genetic risk factor relationships among health indicators and outcomes is achievable through appropriate statistical modeling of large-scale biobank data.

GWAS meta-analyses clarify the genetics of cervical phenotypes and inform risk stratification for cervical cancer.

Genome-wide association studies (GWAS) have successfully identified associations for cervical cancer, but the underlying mechanisms of cervical biology and pathology remain uncharacterised. Our GWAS meta-analyses fill this gap, as we characterise the genetic architecture of cervical phenotypes, including cervical ectropion, cervicitis, cervical dysplasia, as well as up to 9229 cases and 490 304 controls for cervical cancer from diverse ancestries. Leveraging the latest computational methods and gene expression data, we refine the association signals for cervical cancer and propose potential causal variants and genes at each locus. We prioritise PAX8/PAX8-AS1, LINC00339, CDC42, CLPTM1L, HLA-DRB1 and GSDMB as the most likely candidate genes for cervical cancer signals, providing insights into cervical cancer pathogenesis and supporting the involvement of reproductive tract development, immune response and cellular proliferation/apoptosis. We construct a genetic risk score (GRS) that is associated with cervical cancer [hazard ratios (HR) = 3.1 (1.7-5.6) for the top 15% vs lowest 15% of individuals], and with other HPV- and immune-system-related diagnoses in a phenome-wide association study analysis. Our results propose valuable leads for further functional studies and present a GRS for cervical cancer that allows additional risk stratification and could potentially be used to personalise the conventional screening strategies for groups more susceptible to cervical cancer.

The individual and global impact of copy-number variants on complex human traits.

The impact of copy-number variations (CNVs) on complex human traits remains understudied. We called CNVs in 331,522 UK Biobank participants and performed genome-wide association studies (GWASs) between the copy number of CNV-proxy probes and 57 continuous traits, revealing 131 signals spanning 47 phenotypes. Our analysis recapitulated well-known associations (e.g., 1q21 and height), revealed the pleiotropy of recurrent CNVs (e.g., 26 and 16 traits for 16p11.2-BP4-BP5 and 22q11.21, respectively), and suggested gene functionalities (e.g., MARF1 in female reproduction). Forty-eight CNV signals (38%) overlapped with single-nucleotide polymorphism (SNP)-GWASs signals for the same trait. For instance, deletion of PDZK1, which encodes a urate transporter scaffold protein, decreased serum urate levels, while deletion of RHD, which encodes the Rhesus blood group D antigen, associated with hematological traits. Other signals overlapped Mendelian disorder regions, suggesting variable expressivity and broad impact of these loci, as illustrated by signals mapping to Rotor syndrome (SLCO1B1/3), renal cysts and diabetes syndrome (HNF1B), or Charcot-Marie-Tooth (PMP22) loci. Total CNV burden negatively impacted 35 traits, leading to increased adiposity, liver/kidney damage, and decreased intelligence and physical capacity. Thirty traits remained burden associated after correcting for CNV-GWAS signals, pointing to a polygenic CNV architecture. The burden negatively correlated with socio-economic indicators, parental lifespan, and age (survivorship proxy), suggesting a contribution to decreased longevity. Together, our results showcase how studying CNVs can expand biological insights, emphasizing the critical role of this mutational class in shaping human traits and arguing in favor of a continuum between Mendelian and complex diseases.

Improving GWAS discovery and genomic prediction accuracy in biobank data.

Genetically informed, deep-phenotyped biobanks are an important research resource and it is imperative that the most powerful, versatile, and efficient analysis approaches are used. Here, we apply our recently developed Bayesian grouped mixture of regressions model (GMRM) in the UK and Estonian Biobanks and obtain the highest genomic prediction accuracy reported to date across 21 heritable traits. When compared to other approaches, GMRM accuracy was greater than annotation prediction models run in the LDAK or LDPred-funct software by 15% (SE 7%) and 14% (SE 2%), respectively, and was 18% (SE 3%) greater than a baseline BayesR model without single-nucleotide polymorphism (SNP) markers grouped into minor allele frequency-linkage disequilibrium (MAF-LD) annotation categories. For height, the prediction accuracy R2 was 47% in a UK Biobank holdout sample, which was 76% of the estimated [Formula: see text]. We then extend our GMRM prediction model to provide mixed-linear model association (MLMA) SNP marker estimates for genome-wide association (GWAS) discovery, which increased the independent loci detected to 16,162 in unrelated UK Biobank individuals, compared to 10,550 from BoltLMM and 10,095 from Regenie, a 62 and 65% increase, respectively. The average [Formula: see text] value of the leading markers increased by 15.24 (SE 0.41) for every 1% increase in prediction accuracy gained over a baseline BayesR model across the traits. Thus, we show that modeling genetic associations accounting for MAF and LD differences among SNP markers, and incorporating prior knowledge of genomic function, is important for both genomic prediction and discovery in large-scale individual-level studies.

Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry.

BACKGROUND: Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes. RESULTS: The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata. CONCLUSIONS: Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.

The landscape of autosomal-recessive pathogenic variants in European populations reveals phenotype-specific effects.

The number and distribution of recessive alleles in the population for various diseases are not known at genome-wide-scale. Based on 6,447 exome sequences of healthy, genetically unrelated Europeans of two distinct ancestries, we estimate that every individual is a carrier of at least 2 pathogenic variants in currently known autosomal-recessive (AR) genes and that 0.8%-1% of European couples are at risk of having a child affected with a severe AR genetic disorder. This risk is 16.5-fold higher for first cousins but is significantly more increased for skeletal disorders and intellectual disabilities due to their distinct genetic architecture.

Patterns of genetic connectedness between modern and medieval Estonian genomes reveal the origins of a major ancestry component of the Finnish population.

The Finnish population is a unique example of a genetic isolate affected by a recent founder event. Previous studies have suggested that the ancestors of Finnic-speaking Finns and Estonians reached the circum-Baltic region by the 1st millennium BC. However, high linguistic similarity points to a more recent split of their languages. To study genetic connectedness between Finns and Estonians directly, we first assessed the efficacy of imputation of low-coverage ancient genomes by sequencing a medieval Estonian genome to high depth (23×) and evaluated the performance of its down-sampled replicas. We find that ancient genomes imputed from >0.1× coverage can be reliably used in principal-component analyses without projection. By searching for long shared allele intervals (LSAIs; similar to identity-by-descent segments) in unphased data for >143,000 present-day Estonians, 99 Finns, and 14 imputed ancient genomes from Estonia, we find unexpectedly high levels of individual connectedness between Estonians and Finns for the last eight centuries in contrast to their clear differentiation by allele frequencies. High levels of sharing of these segments between Estonians and Finns predate the demographic expansion and late settlement process of Finland. One plausible source of this extensive sharing is the 8th-10th centuries AD migration event from North Estonia to Finland that has been proposed to explain uniquely shared linguistic features between the Finnish language and the northern dialect of Estonian and shared Christianity-related loanwords from Slavic. These results suggest that LSAI detection provides a computationally tractable way to detect fine-scale structure in large cohorts.

Circulating anti-Müllerian hormone levels in pre-menopausal women: novel genetic insights from a genome-wide association meta-analysis.

STUDY QUESTION: Can a genome-wide association study (GWAS) meta-analysis, including a large sample of young premenopausal women from a founder population from Northern Finland, identify novel genetic variants for circulating anti-Müllerian hormone (AMH) levels and provide insights into single-nucleotide polymorphism enrichment in different biological pathways and tissues involved in AMH regulation? SUMMARY ANSWER: The meta-analysis identified a total of six loci associated with AMH levels at P < 5 × 10-8, three of which were novel in or near CHEK2, BMP4, and EIF4EBP1, as well as highlighted significant enrichment in renal system vasculature morphogenesis, and the pituitary gland as the top associated tissue in tissue enrichment analysis. WHAT IS KNOWN ALREADY: AMH is expressed by preantral and small antral stage ovarian follicles in women, and variation in age-specific circulating AMH levels has been associated with several health conditions. However, the biological mechanisms underlying the association between health conditions and AMH levels are not yet fully understood. Previous GWAS have identified loci associated with AMH levels in pre-menopausal women, in or near MCM8, AMH, TEX41, and CDCA7. STUDY DESIGN, SIZE, DURATION: We performed a GWAS meta-analysis for circulating AMH level measurements in 9668 pre-menopausal women. PARTICIPANTS/MATERIALS, SETTING, METHODS: We performed a GWAS meta-analysis in which we combined 2619 AMH measurements (at age 31 years) from a prospective founder population cohort (Northern Finland Birth Cohort 1966, NFBC1966) with a previous GWAS meta-analysis that included 7049 pre-menopausal women (age range 15-48 years) (N = 9668). NFBC1966 AMH measurements were quantified using an automated assay. We annotated the genetic variants, combined different data layers to prioritize potential candidate genes, described significant pathways and tissues enriched by the GWAS signals, identified plausible regulatory roles using colocalization analysis, and leveraged publicly available summary statistics to assess genetic and phenotypic correlations with multiple traits. MAIN RESULTS AND THE ROLE OF CHANCE: Three novel genome-wide significant loci were identified. One of these is in complete linkage disequilibrium with c.1100delC in CHEK2, which is found to be 4-fold enriched in the Finnish population compared to other European populations. We propose a plausible regulatory effect of some of the GWAS variants linked to AMH, as they colocalize with GWAS signals associated with gene expression levels of BMP4, TEX41, and EIFBP41. Gene set analysis highlighted significant enrichment in renal system vasculature morphogenesis, and tissue enrichment analysis ranked the pituitary gland as the top association. LARGE SCALE DATA: The GWAS meta-analysis summary statistics are available for download from the GWAS Catalogue with accession number GCST90428625. LIMITATIONS, REASONS FOR CAUTION: This study only included women of European ancestry and the lack of sufficiently sized relevant tissue data in gene expression datasets hinders the assessment of potential regulatory effects in reproductive tissues. WIDER IMPLICATIONS OF THE FINDINGS: Our results highlight the increased power of founder populations and larger sample sizes to boost the discovery of novel trait-associated variants underlying variation in AMH levels, which aided the characterization of GWAS signals enrichment in different biological pathways and plausible genetic regulatory effects linked with AMH level variation for the first time. STUDY FUNDING/COMPETING INTEREST(S): This work has received funding from the European Union's Horizon 2020 Research and Innovation Programme under the MATER Marie Sklodowska-Curie Grant Agreement No. 813707 and Oulu University Scholarship Foundation and Paulon Säätiö Foundation. (N.P.-G.), Academy of Finland, Sigrid Jusélius Foundation, Novo Nordisk, University of Oulu, Roche Diagnostics (T.T.P.). This work was supported by the Estonian Research Council Grant 1911 (R.M.). J.R. was supported by the European Union's Horizon 2020 Research and Innovation Program under Grant Agreements No. 874739 (LongITools), 824989 (EUCAN-Connect), 848158 (EarlyCause), and 733206 (LifeCycle). U.V. was supported by the Estonian Research Council grant PRG (PRG1291). The NFBC1966 received financial support from University of Oulu Grant No. 24000692, Oulu University Hospital Grant No. 24301140, and ERDF European Regional Development Fund Grant No. 539/2010 A31592. T.T.P. has received grants from Roche, Perkin Elmer, and honoraria for scientific presentations from Gedeon Richter, Exeltis, Astellas, Roche, Stragen, Astra Zeneca, Merck, MSD, Ferring, Duodecim, and Ajaton Terveys. For all other authors, there are no competing interests.

Genome-wide Study Identifies Association between HLA-B∗55:01 and Self-Reported Penicillin Allergy.

Hypersensitivity reactions to drugs are often unpredictable and can be life threatening, underscoring a need for understanding their underlying mechanisms and risk factors. The extent to which germline genetic variation influences the risk of commonly reported drug allergies such as penicillin allergy remains largely unknown. We extracted data from the electronic health records of more than 600,000 participants from the UK, Estonian, and Vanderbilt University Medical Center's BioVU biobanks to study the role of genetic variation in the occurrence of self-reported penicillin hypersensitivity reactions. We used imputed SNP to HLA typing data from these cohorts to further fine map the human leukocyte antigen (HLA) association and replicated our results in 23andMe's research cohort involving a total of 1.12 million individuals. Genome-wide meta-analysis of penicillin allergy revealed two loci, including one located in the HLA region on chromosome 6. This signal was further fine-mapped to the HLA-B∗55:01 allele (OR 1.41 95% CI 1.33-1.49, p value 2.04 × 10-31) and confirmed by independent replication in 23andMe's research cohort (OR 1.30 95% CI 1.25-1.34, p value 1.00 × 10-47). The lead SNP was also associated with lower lymphocyte counts and in silico follow-up suggests a potential effect on T-lymphocytes at HLA-B∗55:01. We also observed a significant hit in PTPN22 and the GWAS results correlated with the genetics of rheumatoid arthritis and psoriasis. We present robust evidence for the role of an allele of the major histocompatibility complex (MHC) I gene HLA-B in the occurrence of penicillin allergy.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Genome-wide association study meta-analysis supports association between MUC1 and ectopic pregnancy.

STUDY QUESTION: Can we identify genetic variants associated with ectopic pregnancy by undertaking the first genome-wide association study (GWAS) leveraging two large-scale biobank initiatives? SUMMARY ANSWER: We identified two novel genome-wide significant associations with ectopic pregnancy, highlighting MUC1 (mucin 1) as the most plausible affected gene. WHAT IS KNOWN ALREADY: Ectopic pregnancy is an important cause of maternal morbidity and mortality worldwide. Despite being a common early pregnancy complication, the genetic predisposition to this condition remains understudied and no large scale genetic studies have been performed so far. STUDY DESIGN, SIZE, DURATION: A GWAS meta-analysis including 7070 women with ectopic pregnancy and 248 810 controls from Estonian Biobank and the FinnGen study. PARTICIPANTS/MATERIALS, SETTING, METHODS: We identified ectopic pregnancy cases from national registers by ICD (International Classification of Disease) codes (ICD-10 O00), and all remaining women were considered controls. We carried out standard GWAS meta-analysis and additionally annotated GWAS signals, analysed co-localization with quantitative trait loci, estimated genetic correlations and identified associated phenotypes to characterize the genetic signals, as well as to analyse the genetic and phenotypic relationships with the condition. MAIN RESULTS AND THE ROLE OF CHANCE: We identified two genome-wide significant loci on chromosomes 1 (rs4971091, P = 5.32×10-9) and 10 (rs11598956, P = 2.41×10-8) potentially associated with ectopic pregnancy. Follow-up analyses propose MUC1, which codes for an epithelial glycoprotein with an important role in barrier function, as the most likely candidate gene for the association on chromosome 1. We also characterize the phenotypic and genetic correlations with other phenotypes, identifying a genetic correlation with smoking and diseases of the (genito)urinary and gastrointestinal system, and phenotypic correlations with various reproductive health diagnoses, reflecting the previously known epidemiological associations. LARGE SCALE DATA: The GWAS meta-analysis summary statistics are available from the GWAS Catalogue (GCST90272883). LIMITATIONS, REASONS FOR CAUTION: The main limitation is that the findings are based on European-based ancestry populations, with limited data on other populations, and we only captured maternal genomes. Additionally, further larger meta-analysis or independent studies are needed to validate these findings. WIDER IMPLICATIONS OF THE FINDINGS: This study encourages the use of large-scale genetic datasets to unravel genetic factors linked to ectopic pregnancy, which is difficult to study in experimental settings. Increased sample size might bring additional genetic factors associating with ectopic pregnancy and inform its heritability. Altogether, our results provide more insight into the biology of ectopic pregnancy and, accordingly, the biological processes governing embryo implantation. STUDY FUNDING/COMPETING INTEREST(S): N.P.G. was supported by MATER Marie Sklodowska-Curie which received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No. 813707. This study was funded by European Union through the European Regional Development Fund Project No. 2014-2020.4.01.15-0012 GENTRANSMED. Computations were performed in the High-Performance Computing Center of University of Tartu. The authors declare no competing interests.

Contribution of schizophrenia polygenic burden to longitudinal phenotypic variance in 22q11.2 deletion syndrome.

While the recurrent 22q11.2 deletion is one of the strongest genetic risk factors for schizophrenia (SCZ), variability of its associated neuropsychiatric endophenotypes reflects its incomplete penetrance for psychosis development. To assess whether this phenotypic variability is linked to common variants associated with SCZ, we studied the association between SCZ polygenic risk score (PRS) and longitudinally acquired phenotypic information of the Swiss 22q11.2DS cohort (n = 97, 50% females, mean age 17.7 yr, mean visit interval 3.8 yr). The SCZ PRS with the best predictive performance was ascertained in the Estonian Biobank (n = 201,146) with LDpred. The infinitesimal SCZ PRS model showed the strongest capacity in discriminating SCZ cases from controls with one SD difference in SCZ PRS corresponding to an odds ratio (OR) of 1.73 (95% CI 1.57-1.90, P = 1.47 × 10-29). In 22q11.2 patients, random-effects ordinal regression modelling using longitudinal data showed SCZ PRS to have the strongest effect on social anhedonia (OR = 2.09, P = 0.0002), and occupational functioning (OR = 1.82, P = 0.0003) within the negative symptoms course, and dysphoric mood (OR = 2.00, P = 0.002) and stress intolerance (OR = 1.76, P = 0.0002) within the general symptoms course. Genetic liability for SCZ was additionally associated with full scale cognitive decline (ß = -0.25, P = 0.02) and with longitudinal volumetric reduction of the right and left hippocampi (ß = -0.28, P = 0.005; ß = -0.23, P = 0.02, respectively). Our results indicate that the polygenic contribution to SCZ acts upon the threshold-lowering first hit (i.e., the deletion). It modifies the endophenotypes of 22q11.2DS and augments the derailment of developmental trajectories of negative and general symptoms, cognition, and hippocampal volume.

Estimation of recurrent atherosclerotic cardiovascular event risk in patients with established cardiovascular disease: the updated SMART2 algorithm.

AIMS: The 10-year risk of recurrent atherosclerotic cardiovascular disease (ASCVD) events in patients with established ASCVD can be estimated with the Secondary Manifestations of ARTerial disease (SMART) risk score, and may help refine clinical management. To broaden generalizability across regions, we updated the existing tool (SMART2 risk score) and recalibrated it with regional incidence rates and assessed its performance in external populations. METHODS AND RESULTS: Individuals with coronary artery disease, cerebrovascular disease, peripheral artery disease, or abdominal aortic aneurysms were included from the Utrecht Cardiovascular Cohort-SMART cohort [n = 8355; 1706 ASCVD events during a median follow-up of 8.2 years (interquartile range 4.2-12.5)] to derive a 10-year risk prediction model for recurrent ASCVD events (non-fatal myocardial infarction, non-fatal stroke, or cardiovascular mortality) using a Fine and Gray competing risk-adjusted model. The model was recalibrated to four regions across Europe, and to Asia (excluding Japan), Japan, Australia, North America, and Latin America using contemporary cohort data from each target region. External validation used data from seven cohorts [Clinical Practice Research Datalink, SWEDEHEART, the international REduction of Atherothrombosis for Continued Health (REACH) Registry, Estonian Biobank, Spanish Biomarkers in Acute Coronary Syndrome and Biomarkers in Acute Myocardial Infarction (BACS/BAMI), the Norwegian COgnitive Impairment After STroke, and Bialystok PLUS/Polaspire] and included 369 044 individuals with established ASCVD of whom 62 807 experienced an ASCVD event. C-statistics ranged from 0.605 [95% confidence interval (CI) 0.547-0.664] in BACS/BAMI to 0.772 (95% CI 0.659-0.886) in REACH Europe high-risk region. The clinical utility of the model was demonstrated across a range of clinically relevant treatment thresholds for intensified treatment options. CONCLUSION: The SMART2 risk score provides an updated, validated tool for the prediction of recurrent ASCVD events in patients with established ASCVD across European and non-European populations. The use of this tool could allow for a more personalized approach to secondary prevention based upon quantitative rather than qualitative estimates of residual risk.

Uniting biobank resources reveals novel genetic pathways modulating susceptibility for atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is a common chronic inflammatory skin disease with high heritability. Previous genome-wide association studies have identified several loci predisposing to AD. These findings explain approximately 30% of the variance in AD susceptibility, suggesting that further work is required to fully understand the genetic underpinnings. OBJECTIVE: We sought to gain additional understanding of the genetic contribution to AD risk by using biobank resources. METHODS: We completed a genome-wide meta-analysis of AD in 796,661 individuals (Ncases = 22,474) from the FinnGen study, the Estonian Biobank, and the UK Biobank. We further performed downstream in silico analyses to characterize the risk variants at the novel loci. RESULTS: We report 30 loci associating with AD (P < 5 × 10-8), 5 of which are novel. In 2 of the novel loci, we identified missense mutations with deleterious predictions in desmocollin 1 and serpin family B member 7, genes encoding proteins crucial to epidermal strength and integrity. CONCLUSIONS: These findings elucidate novel genetic pathways involved in AD pathophysiology. The likely involvement of desmocollin 1 and serpin family B member 7 in AD pathogenesis may offer opportunities for the development of novel treatment strategies for AD in the future.