Rare genetic variants underlie outlying levels of DNA methylation and gene-expression.

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on DNA methylation (DNAm) as effect sizes are expected to be larger for molecular traits compared with complex traits. Here, we investigate DNAm in healthy ageing populations-the Lothian Birth Cohorts of 1921 and 1936-and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic single nucleotide polymorphisms (SNPs) within 1 kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels toward the tails of the population distribution. This study demonstrates the role of rare SNPs in the phenotypic variation of DNAm and the effect of extreme levels of DNAm on gene expression.

Germline variants are associated with increased primary melanoma tumor thickness at diagnosis.

Germline genetic variants have been identified, which predispose individuals and families to develop melanoma. Tumor thickness is the strongest predictor of outcome for clinically localized primary melanoma patients. We sought to determine whether there is a heritable genetic contribution to variation in tumor thickness. If confirmed, this will justify the search for specific genetic variants influencing tumor thickness. To address this, we estimated the proportion of variation in tumor thickness attributable to genome-wide genetic variation (variant-based heritability) using unrelated patients with measured primary cutaneous melanoma thickness. As a secondary analysis, we conducted a genome-wide association study (GWAS) of tumor thickness. The analyses utilized 10 604 individuals with primary cutaneous melanoma drawn from nine GWAS datasets from eight cohorts recruited from the general population, primary care and melanoma treatment centers. Following quality control and filtering to unrelated individuals with study phenotypes, 8125 patients were used in the primary analysis to test whether tumor thickness is heritable. An expanded set of 8505 individuals (47.6% female) were analyzed for the secondary GWAS meta-analysis. Analyses were adjusted for participant age, sex, cohort and ancestry. We found that 26.6% (SE 11.9%, P = 0.0128) of variation in tumor thickness is attributable to genome-wide genetic variation. While requiring replication, a chromosome 11 locus was associated (P < 5 × 10-8) with tumor thickness. Our work indicates that sufficiently large datasets will enable the discovery of genetic variants associated with greater tumor thickness, and this will lead to the identification of host biological processes influencing melanoma growth and invasion.

Genomic characterisation of the overlap of endometriosis with 76 comorbidities identifies pleiotropic and causal mechanisms underlying disease risk.

Comorbid conditions can be driven by underlying pleiotropic and causal mechanisms that can provide insights into shared molecular and biological processes contributing to disease risk. Endometriosis is a chronic condition affecting one in nine women of reproductive age and poses many challenges including lengthy diagnostic delays and limited treatment efficacy owing to poor understanding of disease aetiology. To shed light on the underlying biological mechanisms and to identify potential risk factors, we examine the epidemiological and genomic relationship between endometriosis and its comorbidities. In the UK Biobank 292 ICD10 codes were epidemiologically correlated with endometriosis diagnosis, including gynaecological, immune, infection, pain, psychiatric, cancer, gastrointestinal, urinary, bone and cardiovascular traits. A subset of the identified comorbidities (n = 76) underwent follow-up genetic analysis. Whilst Mendelian randomisation suggested causality was not responsible for most comorbid relationships, 22 traits were genetically correlated with endometriosis, including pain, gynaecological and gastrointestinal traits, suggestive of a shared genetic background. Pleiotropic genetic variants and genes were identified using gene-based and colocalisation analysis. Shared genetic risk factors and potential target genes suggest a diverse collection of biological systems are involved in these comorbid relationships including coagulation factors, development of the female reproductive tract and cell proliferation. These findings highlight the diversity of traits with epidemiological and genomic overlap with endometriosis and implicate a key role for pleiotropy in the comorbid relationships.

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.

Polygenic risk score phenome-wide association study reveals an association between endometriosis and testosterone.

BACKGROUND: Endometriosis affects 1 in 9 women, yet it is poorly understood with long diagnostic delays, invasive diagnoses, and poor treatment outcomes. Characterised by the presence of endometrial-like tissue outside of the uterus, its main symptoms are pain and infertility. Endometriosis often co-occurs with other conditions, which may provide insights into the origins of endometriosis. METHODS: Here a polygenic risk score phenome-wide association study of endometriosis was conducted in the UK Biobank to investigate the pleiotropic effects of a genetic liability to endometriosis. The relationship between the polygenic risk score for endometriosis and health conditions, blood and urine biomarkers and reproductive factors were investigated separately in females, males and females without an endometriosis diagnosis. The relationship between endometriosis and the blood and urine biomarkers was further investigated using genetic correlation and Mendelian randomisation approaches to identify causal relationships. RESULTS: Multiple health conditions, blood and urine biomarkers and reproductive factors were associated with genetic liability to endometriosis in each group, indicating many endometriosis comorbidities are not dependent on the physical manifestation of endometriosis. Differences in the associated traits between males and females highlighted the importance of sex-specific pathways in the overlap of endometriosis with many other traits. Notably, an association of genetic liability to endometriosis with lower testosterone levels was identified. Follow-up analysis utilising Mendelian randomisation approaches suggested lower testosterone may be causal for both endometriosis and clear cell ovarian cancer. CONCLUSIONS: This study highlights the diversity of the pleiotropic effects of genetic risk to endometriosis irrespective of a diagnosis of endometriosis. A key finding was the identification of a causal effect of the genetic liability to lower testosterone on endometriosis using Mendelian randomisation.

Gene expression profiles separate endometriosis lesion subtypes and indicate a sensitivity of endometrioma to estrogen suppressive treatments through elevated ESR2 expression.

BACKGROUND: Endometriosis is a common, gynaecological disease characterised by the presence of endometrial-like cells growing outside the uterus. Lesions appear at multiple locations, present with variation in appearance, size and depth of invasion. Despite hormones being the recommended first-line treatment, their efficacy, success and side effects vary widely amongst study populations. Current, hormonal medication for endometriosis is designed to suppress systemic oestrogen. Whether these hormones can influence the lesions themselves is not yet clear. Evidence of hormone receptor expression in endometriotic lesions and their ability to respond is conflicting. A variation in their expression, activation of transcriptional co-regulators and the potential to respond may contribute to their variation in patient outcomes. Identifying patients who would benefit from hormonal treatments remain an important goal in endometriosis research. METHODS: Using gene expression data from endometriosis lesions including endometrioma (OMA, n = 28), superficial peritoneal lesions (SUP, n = 72) and deeply infiltrating lesions (DIE, n = 78), we performed principal component analysis, differential gene expression and gene correlation analyses to assess the impact of menstrual stage, lesion subtype and hormonal treatment on the gene expression. RESULTS: The gene expression profiles did not vary based on menstrual stage, but could distinguish lesion subtypes with OMA significantly differentiating from both SUP and DIE. Additionally, the effect of oestrogen suppression medication altered the gene expression profile in OMA, while such effect was not observed in SUP or DIE. Analysis of the target receptors for hormonal medication indicated ESR2 was differentially expressed in OMA and that genes that correlated with ESR2 varied significantly between medicated and non-medicated OMA samples. CONCLUSIONS: Our results demonstrate of the different lesion types OMA present with strongest response to hormonal treatment directly through ESR2. The data suggests that there may be the potential to target treatment options to individual patients based on pre-surgical diagnoses.

GWAS of Dizygotic Twinning in an Enlarged Australian Sample of Mothers of DZ Twins.

Female fertility is a complex trait with age-specific changes in spontaneous dizygotic (DZ) twinning and fertility. To elucidate factors regulating female fertility and infertility, we conducted a genome-wide association study (GWAS) on mothers of spontaneous DZ twins (MoDZT) versus controls (3273 cases, 24,009 controls). This is a follow-up study to the Australia/New Zealand (ANZ) component of that previously reported (Mbarek et al., 2016), with a sample size almost twice that of the entire discovery sample meta-analysed in the previous article (and five times the ANZ contribution to that), resulting from newly available additional genotyping and representing a significant increase in power. We compare analyses with and without male controls and show unequivocally that it is better to include male controls who have been screened for recent family history, than to use only female controls. Results from the SNP based GWAS identified four genomewide significant signals, including one novel region, ZFPM1 (Zinc Finger Protein, FOG Family Member 1), on chromosome 16. Previous signals near FSHB (Follicle Stimulating Hormone beta subunit) and SMAD3 (SMAD Family Member 3) were also replicated (Mbarek et al., 2016). We also ran the GWAS with a dominance model that identified a further locus ADRB2 on chr 5. These results have been contributed to the International Twinning Genetics Consortium for inclusion in the next GWAS meta-analysis (Mbarek et al., in press).

New concepts on the etiology of endometriosis.

Endometriosis is a serious, chronic disorder where endometrial tissue grows outside the uterus, causing severe pelvic pain and infertility. It affects 11% of women. Endometriosis is a multifactorial disorder of unclear etiology, although retrograde menstruation plays a major role. It has a genetic component with over 40 genetic risk factors mapped, although their mechanism of action is still emerging. New evidence suggests a role for retrograde menstruation of endometrial stem/progenitor cells, now that identifying markers of these cells are available. Recent lineage tracing and tissue clearing microscopy and 3D reconstruction has provided new understanding of endometrial glandular structure, particularly the horizontal orientation and interconnection of basalis glands. New sequencing technologies, particularly whole genome DNA sequencing are revealing somatic mutations, including in cancer driver genes, in normal and eutopic endometrium of patients with endometriosis, as well as ectopic endometriotic lesions. Methylome sequencing is offering insight into the regulation of genes and the role of the environmental factors. Single cell RNA sequencing reveals the transcriptome of individual endometrial cells, shedding new light on the diversity and range of cellular subpopulations of the major cell types present in the endometrium and in endometriotic lesions. New endometrial epithelial organoid cultures replicating glandular epithelium are providing tractable models for studying endometriosis. Organoids derived from menstrual fluid offer a non-invasive source of endometrial tissue and a new avenue for testing drugs and developing personalized medicine for treating endometriosis. These new approaches are rapidly advancing our understanding of endometriosis etiology.

Multiplex melanoma families are enriched for polygenic risk.

Cancers, including cutaneous melanoma, can cluster in families. In addition to environmental etiological factors such as ultraviolet radiation, cutaneous melanoma has a strong genetic component. Genetic risks for cutaneous melanoma range from rare, high-penetrance mutations to common, low-penetrance variants. Known high-penetrance mutations account for only about half of all densely affected cutaneous melanoma families, and the causes of familial clustering in the remainder are unknown. We hypothesize that some clustering is due to the cumulative effect of a large number of variants of individually small effect. Common, low-penetrance genetic risk variants can be combined into polygenic risk scores. We used a polygenic risk score for cutaneous melanoma to compare families without known high-penetrance mutations with unrelated melanoma cases and melanoma-free controls. Family members had significantly higher mean polygenic load for cutaneous melanoma than unrelated cases or melanoma-free healthy controls (Bonferroni-corrected t-test P = 1.5 × 10-5 and 6.3 × 10-45, respectively). Whole genome sequencing of germline DNA from 51 members of 21 families with low polygenic risk for melanoma identified a CDKN2A p.G101W mutation in a single family but no other candidate high-penetrance melanoma susceptibility genes. This work provides further evidence that melanoma, like many other common complex disorders, can arise from the joint action of multiple predisposing factors, including rare high-penetrance mutations, as well as via a combination of large numbers of alleles of small effect.

Molecular Support for Heterogonesis Resulting in Sesquizygotic Twinning.

Sesquizygotic multiple pregnancy is an exceptional intermediate between monozygotic and dizygotic twinning. We report a monochorionic twin pregnancy with fetal sex discordance. Genotyping of amniotic fluid from each sac showed that the twins were maternally identical but chimerically shared 78% of their paternal genome, which makes them genetically in between monozygotic and dizygotic; they are sesquizygotic. We observed no evidence of sesquizygosis in 968 dizygotic twin pairs whom we screened by means of pangenome single-nucleotide polymorphism genotyping. Data from published repositories also show that sesquizygosis is a rare event. Detailed genotyping implicates chimerism arising at the juncture of zygotic division, termed heterogonesis, as the likely initial step in the causation of sesquizygosis.

Changing genetic architecture of body mass index from infancy to early adulthood: an individual based pooled analysis of 25 twin cohorts.

BACKGROUND: Body mass index (BMI) shows strong continuity over childhood and adolescence and high childhood BMI is the strongest predictor of adult obesity. Genetic factors strongly contribute to this continuity, but it is still poorly known how their contribution changes over childhood and adolescence. Thus, we used the genetic twin design to estimate the genetic correlations of BMI from infancy to adulthood and compared them to the genetic correlations of height. METHODS: We pooled individual level data from 25 longitudinal twin cohorts including 38,530 complete twin pairs and having 283,766 longitudinal height and weight measures. The data were analyzed using Cholesky decomposition offering genetic and environmental correlations of BMI and height between all age combinations from 1 to 19 years of age. RESULTS: The genetic correlations of BMI and height were stronger than the trait correlations. For BMI, we found that genetic correlations decreased as the age between the assessments increased, a trend that was especially visible from early to middle childhood. In contrast, for height, the genetic correlations were strong between all ages. Age-to-age correlations between environmental factors shared by co-twins were found for BMI in early childhood but disappeared altogether by middle childhood. For height, shared environmental correlations persisted from infancy to adulthood. CONCLUSIONS: Our results suggest that the genes affecting BMI change over childhood and adolescence leading to decreasing age-to-age genetic correlations. This change is especially visible from early to middle childhood indicating that new genetic factors start to affect BMI in middle childhood. Identifying mediating pathways of these genetic factors can open possibilities for interventions, especially for those children with high genetic predisposition to adult obesity.

Co-Inheritance of Variation in All-Cause Mortality and Biochemical Risk Factors.

Biomarkers may be useful endophenotypes for genetic studies if they share genetic sources of variation with the outcome, for example, with all-cause mortality. Australian adult study participants who had reported their parental survival information were included in the study: 14,169 participants had polygenic risk scores (PRS) from genotyping and up to 13,365 had biomarker results. We assessed associations between participants' biomarker results and parental survival, and between biomarker results and eight parental survival PRS at varying p-value cut-offs. Survival in parents was associated with participants' serum bilirubin, C-reactive protein, HDL cholesterol, triglycerides and uric acid, and with LDL cholesterol for participants' fathers but not for their mothers. PRS for all-cause mortality were associated with liver function tests (alkaline phosphatase, butyrylcholinesterase, gamma-glutamyl transferase), metabolic tests (LDL and HDL cholesterol, triglycerides, uric acid), and acute-phase reactants (C-reactive protein, globulins). Association between offspring biomarker results and parental survival demonstrates the existence of familial effects common to both, while associations between biomarker results and PRS for mortality favor at least a partial genetic cause of this covariation. Identification of genetic loci affecting mortality-associated biomarkers offers a route to the identification of additional loci affecting mortality.

Genetic analysis of endometriosis and depression identifies shared loci and implicates causal links with gastric mucosa abnormality.

Evidence from observational studies indicates that endometriosis and depression often co-occur. However, conflicting evidence exists, and the etiology as well as biological mechanisms underlying their comorbidity remain unknown. Utilizing genome-wide association study (GWAS) data, we comprehensively assessed the relationship between endometriosis and depression. Single nucleotide polymorphism effect concordance analysis (SECA) found a significant genetic overlap between endometriosis and depression (PFsig-permuted = 9.99 × 10-4). Linkage disequilibrium score regression (LDSC) analysis estimated a positive and highly significant genetic correlation between the two traits (rG = 0.27, P = 8.85 × 10-27). A meta-analysis of endometriosis and depression GWAS (sample size = 709,111), identified 20 independent genome-wide significant loci (P < 5 × 10-8), of which eight are novel. Mendelian randomization analysis (MR) suggests a causal effect of depression on endometriosis. Combining gene-based association results across endometriosis and depression GWAS, we identified 22 genes with a genome-wide significant Fisher's combined P value (FCPgene < 2.75 × 10-6). Genes with a nominal gene-based association (Pgene < 0.05) were significantly enriched across endometriosis and depression (Pbinomial-test = 2.90 × 10-4). Also, genes overlapping the two traits at Pgene < 0.1 (Pbinomial-test = 1.31 × 10-5) were significantly enriched for the biological pathways 'cell-cell adhesion', 'inositol phosphate metabolism', 'Hippo-Merlin signaling dysregulation' and 'gastric mucosa abnormality'. These results reveal a shared genetic etiology for endometriosis and depression. Indeed, additional analyses found evidence of a causal association between each of endometriosis and depression and at least one abnormal condition of gastric mucosa. Our study confirms the comorbidity of endometriosis and depression, implicates links with gastric mucosa abnormalities in their causal pathways and reveals potential therapeutic targets for further investigation.

Genetic analyses of gynecological disease identify genetic relationships between uterine fibroids and endometrial cancer, and a novel endometrial cancer genetic risk region at the WNT4 1p36.12 locus.

Endometriosis, polycystic ovary syndrome (PCOS) and uterine fibroids have been proposed as endometrial cancer risk factors; however, disentangling their relationships with endometrial cancer is complicated due to shared risk factors and comorbidities. Using genome-wide association study (GWAS) data, we explored the relationships between these non-cancerous gynecological diseases and endometrial cancer risk by assessing genetic correlation, causal relationships and shared risk loci. We found significant genetic correlation between endometrial cancer and PCOS, and uterine fibroids. Adjustment for genetically predicted body mass index (a risk factor for PCOS, uterine fibroids and endometrial cancer) substantially attenuated the genetic correlation between endometrial cancer and PCOS but did not affect the correlation with uterine fibroids. Mendelian randomization analyses suggested a causal relationship between only uterine fibroids and endometrial cancer. Gene-based analyses revealed risk regions shared between endometrial cancer and endometriosis, and uterine fibroids. Multi-trait GWAS analysis of endometrial cancer and the genetically correlated gynecological diseases identified a novel genome-wide significant endometrial cancer risk locus at 1p36.12, which replicated in an independent endometrial cancer dataset. Interrogation of functional genomic data at 1p36.12 revealed biologically relevant genes, including WNT4 which is necessary for the development of the female reproductive system. In summary, our study provides genetic evidence for a causal relationship between uterine fibroids and endometrial cancer. It further provides evidence that the comorbidity of endometrial cancer, PCOS and uterine fibroids may partly be due to shared genetic architecture. Notably, this shared architecture has revealed a novel genome-wide risk locus for endometrial cancer.

Educational attainment of same-sex and opposite-sex dizygotic twins: An individual-level pooled study of 19 twin cohorts.

Comparing twins from same- and opposite-sex pairs can provide information on potential sex differences in a variety of outcomes, including socioeconomic-related outcomes such as educational attainment. It has been suggested that this design can be applied to examine the putative role of intrauterine exposure to testosterone for educational attainment, but the evidence is still disputed. Thus, we established an international database of twin data from 11 countries with 88,290 individual dizygotic twins born over 100 years and tested for differences between twins from same- and opposite-sex dizygotic pairs in educational attainment. Effect sizes with 95% confidence intervals (CI) were estimated by linear regression models after adjusting for birth year and twin study cohort. In contrast to the hypothesis, no difference was found in women (ß = -0.05 educational years, 95% CI -0.11, 0.02). However, men with a same-sex co-twin were slightly more educated than men having an opposite-sex co-twin (ß = 0.14 educational years, 95% CI 0.07, 0.21). No consistent differences in effect sizes were found between individual twin study cohorts representing Europe, the USA, and Australia or over the cohorts born during the 20th century, during which period the sex differences in education reversed favoring women in the latest birth cohorts. Further, no interaction was found with maternal or paternal education. Our results contradict the hypothesis that there would be differences in the intrauterine testosterone levels between same-sex and opposite-sex female twins affecting education. Our findings in men may point to social dynamics within same-sex twin pairs that may benefit men in their educational careers.

Genetic Regulation of Physiological Reproductive Lifespan and Female Fertility.

There is substantial genetic variation for common traits associated with reproductive lifespan and for common diseases influencing female fertility. Progress in high-throughput sequencing and genome-wide association studies (GWAS) have transformed our understanding of common genetic risk factors for complex traits and diseases influencing reproductive lifespan and fertility. The data emerging from GWAS demonstrate the utility of genetics to explain epidemiological observations, revealing shared biological pathways linking puberty timing, fertility, reproductive ageing and health outcomes. The observations also identify unique genetic risk factors specific to different reproductive diseases impacting on female fertility. Sequencing in patients with primary ovarian insufficiency (POI) have identified mutations in a large number of genes while GWAS have revealed shared genetic risk factors for POI and ovarian ageing. Studies on age at menopause implicate DNA damage/repair genes with implications for follicle health and ageing. In addition to the discovery of individual genes and pathways, the increasingly powerful studies on common genetic risk factors help interpret the underlying relationships and direction of causation in the regulation of reproductive lifespan, fertility and related traits.

Genome-wide survey of parent-of-origin effects on DNA methylation identifies candidate imprinted loci in humans.

Genomic imprinting is an epigenetic mechanism leading to parent-of-origin silencing of alleles. So far, the precise number of imprinted regions in humans is uncertain. In this study, we leveraged genome-wide DNA methylation in whole blood measured longitudinally at three time points (birth, childhood and adolescence) and genome-wide association studies (GWAS) data in 740 mother-child duos from the Avon Longitudinal Study of parents and children to identify candidate imprinted loci. We reasoned that cis-meQTLs at genomic regions that were imprinted would show strong evidence of parent-of-origin associations with DNA methylation, enabling the detection of imprinted regions. Using this approach, we identified genome-wide significant cis-meQTLs that exhibited parent-of-origin effects (POEs) at 82 loci, 34 novel and 48 regions previously implicated in imprinting (3.7-10

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero.

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (ß = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.