Anthracycline Toxicity: Light at the End of the Tunnel?

Anthracycline-induced cardiotoxicity (AIC) is a serious and common side effect of anthracycline therapy. Identification of genes and genetic variants associated with AIC risk has clinical potential as a cardiotoxicity predictive tool and to allow the development of personalized therapies. In this review, we provide an overview of the function of known AIC genes identified by association studies and categorize them based on their mechanistic implication in AIC. We also discuss the importance of functional validation of AIC-associated variants in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) to advance the implementation of genetic predictive biomarkers. Finally, we review how patient-specific hiPSC-CMs can be used to identify novel patient-relevant functional targets and for the discovery of cardioprotectant drugs to prevent AIC. Implementation of functional validation and use of hiPSC-CMs for drug discovery will identify the next generation of highly effective and personalized cardioprotectants and accelerate the inclusion of approved AIC biomarkers into clinical practice.

The burden of cardiovascular disease and risk for subsequent major adverse cardiovascular events in survivors of childhood cancer: a prospective, longitudinal analysis from the St Jude Lifetime Cohort Study.

BACKGROUND: The effect of the increasing lifetime burden of non-major cardiovascular conditions on risk for a subsequent major adverse cardiovascular event among survivors of childhood cancer has not been assessed. We aimed to characterise the prevalence of major adverse cardiovascular events and their association with the cumulative burden of non-major adverse cardiovascular events in childhood cancer survivors. METHODS: This is a longitudinal cohort study with participant data obtained from an ongoing cohort study at St Jude Children's Research Hospital: the St Jude Lifetime Cohort Study (SJLIFE). Prospective clinical follow-up was of 5-year survivors of childhood cancer who were diagnosed when aged younger than 25 years from 1962 to 2012. Age-frequency, sex-frequency, and race-frequency matched community-control participants completed a similar one-time clinical assessment. 22 cardiovascular events were graded using a St Jude Children's Research Hospital-modified version of the National Cancer Institute Common Terminology Criteria for Adverse Events (version 4.03). Cumulative incidence and burden of the primary outcome of major adverse cardiovascular events (cardiomyopathy, myocardial infarction, stroke, and other cardiovascular-related mortality) were estimated. Rate ratios (RR) of the association of major adverse cardiovascular events with 22 non-major adverse cardiovascular events were estimated using multivariable piecewise-exponential regression adjusting for attained age, age at diagnosis, sex, race and ethnicity, treatment era, diagnosis of diabetes, and exposure to cardiotoxic cancer therapies. The St Jude Lifetime Cohort study is registered with ClinicalTrials.gov, NCT00760656, and is ongoing. FINDINGS: 9602 5-year survivors of childhood cancer, and 737 community controls were included in the longitudinal follow-up (from Sept 13, 2007, to Dec 17, 2021). The median follow-up was 20·3 years (IQR 12·0-31·4) from the date of primary cancer diagnosis (4311 [44.9%] were females). By the age of 50 years (analysis stopped at age 50 years due to the low number of participants older than that age), the cumulative incidence of major adverse cardiovascular events among survivors was 17·7% (95% CI 15·9-19·5) compared with 0·9% (0·0-2·1) in the community controls. The cumulative burden of major adverse cardiovascular events in survivors was 0·26 (95% CI 0·23-0·29) events per survivor compared with 0·009 (0·000-0·021) events per community control participant. Increasing cumulative burden of grade 1-4 non-major adverse cardiovascular events was associated with an increased future risk of major adverse cardiovascular events (one condition: RR 4·3, 95% CI 3·1-6·0; p<0·0001; two conditions: 6·6, 4·6-9·5; p<0·0001; and three conditions: 7·7, 5·1-11·4; p<0·0001). Increased risk for major adverse cardiovascular events was observed with specific subclinical conditions (eg, grade 1 arrhythmias [RR 1·5, 95% CI 1·2-2·0; p=0·0017]), grade 2 left ventricular systolic dysfunction (2·2, 1·6-3·1; p<0·0001), grade 2 valvular disorders (2·2, 1·2-4·0; p=0·013), but not grade 1 hypercholesterolaemia, grade 1-2 hypertriglyceridaemia, or grade 1-2 vascular stenosis. INTERPRETATION: Among an ageing cohort of survivors of childhood cancer, the accumulation of non-major adverse cardiovascular events, including subclinical conditions, increased the risk of major adverse cardiovascular events and should be the focus of interventions for early detection and prevention of major adverse cardiovascular events. FUNDING: The US National Cancer Institute and the American Lebanese Syrian Associated Charities.

Concomitant Wilms tumor and autosomal dominant polycystic kidney disease.

BACKGROUND: Concomitant Wilms tumor (WT) and autosomal dominant polycystic kidney disease (ADPKD) is exceedingly rare, presenting a diagnostic and technical challenge to pediatric surgical oncologists. The simultaneous workup and management of these disease processes are incompletely described. PROCEDURE: We performed a retrospective analysis of patients treated at our institution with concomitant diagnoses of WT and ADPKD. We also review the literature on the underlying biology and management principles of these conditions. RESULTS: We present three diverse cases of concomitant unilateral WT and ADPKD who underwent nephrectomy. One patient had preoperative imaging consistent with ADPKD with confirmatory testing postoperatively, one was found to have contralateral renal cysts intraoperatively with confirmatory imaging post nephrectomy, and one was diagnosed in childhood post nephrectomy. All patients are alive at last follow-up, and the patient with longest follow-up has progressed to end-stage kidney failure requiring transplantation and dialysis in adulthood. All patients underwent germline testing and were found to have no cancer predisposition syndrome or pathogenic or likely pathogenic variants for WT. CONCLUSION: Concomitant inheritance of ADPKD and development of WT are extremely rare, and manifestations of ADPKD may not present until late childhood or adulthood. ADPKD is not a known predisposing condition for WT. When ADPKD diagnosis is made by family history, imaging, and/or genetic testing before WT diagnosis and treatment, the need for extensive preoperative characterization of cystic kidney lesions in children and increased risk of post-nephrectomy kidney failure warrant further discussion of surgical approach and perioperative management strategies.

Keratinocyte carcinomas in survivors of childhood cancer: A report from the childhood cancer survivor study.

BACKGROUND: Childhood cancer survivors (CCS) are at increased risk for keratinocyte carcinomas (KC) however, the long-term incidence of single and multiple KC is not well established. OBJECTIVE: Identify risk factors and quantify KC cumulative incidence and multiple-incidence burden in CCS. METHODS: KC were identified among Childhood Cancer Survivor Study participants, a cohort of 5-year cancer survivors diagnosed <21 years of age between 1970 and 1999 in North America. Cumulative incidence was estimated and multivariable models assessed relative rates of KC associated with survivor and treatment characteristics. RESULTS: Among 25,658 participants, 1446 developed 5363 KC (93.5% basal cell carcinoma, 6.7% squamous cell carcinoma; mean age 37.0 years (range 7.3-67.4), mean latency 25.7 years; 95.3% White and 88.4% with radiotherapy). Mean lesion count was 3.7 with 26.1% experiencing ≥4. Radiotherapy imparted a 4.5-fold increase in the rate of any KC and 9.4-fold increase in the rate of ≥4 KC. Allogeneic and autologous hematopoietic cell transplant were associated with a 3.4- and 2.3-fold increased rate of KC, respectively. LIMITATIONS: Participant self-reporting of some data including race without skin phototype and past medical history may have impacted analysis. CONCLUSIONS: The burden of KC in CCS remains high, but predictable risk factors should guide screening.

Development and validation of age-specific risk prediction models for primary ovarian insufficiency in long-term survivors of childhood cancer: a report from the Childhood Cancer Survivor Study and St Jude Lifetime Cohort.

BACKGROUND: Female survivors of childhood cancer are at risk for primary ovarian insufficiency (POI), defined as the cessation of gonadal function before the age of 40 years. We aimed to develop and validate models to predict age-specific POI risk among long-term survivors of childhood cancer. METHODS: To develop models to predict age-specific POI risk for the ages of 21-40 years, we used data from the Childhood Cancer Survivor Study (CCSS). Female survivors aged 18 years or older at their latest follow-up, with self-reported menstrual history information and free of subsequent malignant neoplasms within 5 years of diagnosis, were included. We evaluated models that used algorithms based on statistical or machine learning to consider all predictors, including cancer treatments. Cross-validated prediction performance metrics (eg, area under the receiver operating characteristic curve [AUROC]) were compared to select the best-performing models. For external validation of the models, we used data from 5-year survivors in the St Jude Lifetime Cohort (SJLIFE) with ovarian status clinically ascertained using hormone measurements (menopause defined by follicle stimulating hormone >30 mIU/mL and oestradiol <17 pg/mL) and medical chart or questionnaire review. We also evaluated an SJLIFE-based polygenic risk score for POI among 1985 CCSS survivors with genotype data available. FINDINGS: 7891 female CCSS survivors (922 with POI) were included in the development of the POI risk prediction model, and 1349 female SJLIFE survivors (101 with POI) were included in the validation study. Median follow-up from cancer diagnosis was 23·7 years (IQR 18·3-30·0) in CCSS and 15·1 years (10·4-22·9) in SJLIFE. Between the ages of 21 and 40 years, POI prevalence increased from 7·9% (95% CI 7·3-8·5) to 18·6% (17·3-20·0) in CCSS and 7·3% (5·8-8·9) to 14·9% (11·6-19·1) in SJLIFE. Age-specific logistic regression models considering ovarian radiation dosimetry or prescribed pelvic and abdominal radiation dose, along with individual chemotherapy predictors, performed well in CCSS. In the SJLIFE validation, the prescribed radiation dose model performed well (AUROC 0·88-0·95), as did a simpler model that considered any exposures to pelvic or abdominal radiotherapy or alkylators (0·82-0·90). Addition of the polygenic risk predictor significantly improved the average positive predictive value (from 0·76 [95% CI 0·63-0·89] to 0·87 [0·80-0·94]; p=0·029) among CCSS survivors treated with ovarian radiation and chemotherapy. INTERPRETATION: POI risk prediction models using treatment information showed robust prediction performance in adult survivors of childhood cancer. FUNDING: Canadian Institutes of Health Research, US National Cancer Institute.

A novel transcription factor combination for direct reprogramming to a spontaneously contracting human cardiomyocyte-like state.

The reprogramming of somatic cells to a spontaneously contracting cardiomyocyte-like state using defined transcription factors has proven successful in mouse fibroblasts. However, this process has been less successful in human cells, thus limiting the potential clinical applicability of this technology in regenerative medicine. We hypothesized that this issue is due to a lack of cross-species concordance between the required transcription factor combinations for mouse and human cells. To address this issue, we identified novel transcription factor candidates to induce cell conversion between human fibroblasts and cardiomyocytes, using the network-based algorithm Mogrify. We developed an automated, high-throughput method for screening transcription factor, small molecule, and growth factor combinations, utilizing acoustic liquid handling and high-content kinetic imaging cytometry. Using this high-throughput platform, we screened the effect of 4960 unique transcription factor combinations on direct conversion of 24 patient-specific primary human cardiac fibroblast samples to cardiomyocytes. Our screen revealed the combination of MYOCD, SMAD6, and TBX20 (MST) as the most successful direct reprogramming combination, which consistently produced up to 40% TNNT2+ cells in just 25 days. Addition of FGF2 and XAV939 to the MST cocktail resulted in reprogrammed cells with spontaneous contraction and cardiomyocyte-like calcium transients. Gene expression profiling of the reprogrammed cells also revealed the expression of cardiomyocyte associated genes. Together, these findings indicate that cardiac direct reprogramming in human cells can be achieved at similar levels to those attained in mouse fibroblasts. This progress represents a step forward towards the clinical application of the cardiac direct reprogramming approach.

Generalizability of "GWAS Hits" in Clinical Populations: Lessons from Childhood Cancer Survivors.

With mounting interest in translating genome-wide association study (GWAS) hits from large meta-analyses (meta-GWAS) in diverse clinical settings, evaluating their generalizability in target populations is crucial. Here, we consider long-term survivors of childhood cancers from the St. Jude Lifetime Cohort Study, and we show the limited generalizability of 1,376 robust SNP associations reported in the general population across 12 complex anthropometric and cardiometabolic phenotypes (n = 2,231; observed-to-expected replication ratio = 0.70, p = 6.2 × 10-8). An examination of five comparable phenotypes in a second independent cohort of survivors from the Childhood Cancer Survivor Study corroborated the overall limited generalizability of meta-GWAS hits to survivors (n = 4,212; observed-to-expected replication ratio = 0.55, p = 5.6 × 10-15). Finally, in direct comparisons of survivor samples against independent equivalently powered general population samples from the UK Biobank, we consistently observed lower meta-GWAS hit replication rates and poorer polygenic risk score predictive performance in survivor samples for multiple phenotypes. As a possible explanation, we found that meta-GWAS hits were less likely to be replicated in survivors who had been exposed to cancer therapies that are associated with phenotype risk. Examination of complementary DNA methylation data in a subset of survivors revealed that treatment-related methylation patterns at genomic sites linked to meta-GWAS hits may disrupt established genetic signals in survivors.

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.

Clinical and genetic risk factors for radiation-associated ototoxicity: A report from the Childhood Cancer Survivor Study and the St. Jude Lifetime Cohort.

BACKGROUND: Cranial radiation therapy (CRT) is associated with ototoxicity, which manifests as hearing loss and tinnitus. The authors sought to identify clinical determinants and genetic risk factors for ototoxicity among adult survivors of pediatric cancer treated with CRT. METHODS: Logistic regression evaluated associations of tinnitus (n = 1991) and hearing loss (n = 2198) with nongenetic risk factors and comorbidities among CRT-treated survivors in the Childhood Cancer Survivor Study. Genome-wide association studies (GWASs) of CRT-related tinnitus and hearing loss were also performed. RESULTS: Males were more likely to report CRT-related tinnitus (9.4% vs 5.4%; P = 5.1 × 10-4 ) and hearing loss (14.0% vs 10.7%; P = .02) than females. Survivors with tinnitus or hearing loss were more likely to experience persistent dizziness or vertigo (tinnitus: P < 2 × 10-16 ; hearing loss: P = 6.4 × 10-9 ), take antidepressants (tinnitus: P = .02; hearing loss: P = .01), and report poorer overall health (tinnitus: P = 1.5 × 10-6 ; hearing loss: P = 1.7 × 10-6 ) in comparison with controls. GWAS of CRT-related tinnitus revealed a genome-wide significant signal in chromosome 1 led by rs203248 (P = 1.5 × 10-9 ), whereas GWAS of CRT-related hearing loss identified rs332013 (P = 5.8 × 10-7 ) in chromosome 8 and rs67522722 (P = 7.8 × 10-7 ) in chromosome 6 as nearly genome-wide significant. A replication analysis identified rs67522722, intronic to ATXN1, as being significantly associated with CRT-related hearing loss (P = .03) and de novo hearing loss (P = 3.6 × 10-4 ). CONCLUSIONS: CRT-associated ototoxicity was associated with sex, several neuro-otological symptoms, increased antidepressant use, and poorer self-reported health. GWAS of CRT-related hearing loss identified rs67522722, which was supported in an independent cohort of survivors. LAY SUMMARY: Hearing loss and subjective tinnitus (the perception of noise or ringing in the ear) are long-term side effects of cancer treatment and are common in children treated with radiation to the brain. These toxicities can affect childhood development and potentially contribute to serious learning and behavioral difficulties. This study's data indicate that males are at greater risk for hearing loss and tinnitus than females after radiation therapy to the brain. Those who develop these toxicities are more likely to use antidepressants and report poorer overall health. Health care providers can improve the management of survivors by informing patients and/or their parents of these risks.

Genetic variation in the body mass index of adult survivors of childhood acute lymphoblastic leukemia: A report from the Childhood Cancer Survivor Study and the St. Jude Lifetime Cohort.

BACKGROUND: Treatment characteristics such as cranial radiation therapy (CRT) do not fully explain adiposity risk in childhood acute lymphoblastic leukemia (ALL) survivors. This study was aimed at characterizing genetic variation related to adult body mass index (BMI) among survivors of childhood ALL. METHODS: Genetic associations of BMI among 1458 adult survivors of childhood ALL (median time from diagnosis, 20 years) were analyzed by multiple approaches. A 2-stage genome-wide association study in the Childhood Cancer Survivor Study (CCSS) and the St. Jude Lifetime Cohort Study (SJLIFE) was performed. BMI was a highly polygenic trait in the general population. Within the known loci, the BMI percent variance explained was estimated, and additive interactions (chi-square test) with CRT in the CCSS were evaluated. The role of DNA methylation in CRT interaction was further evaluated in a subsample of ALL survivors. RESULTS: In a meta-analysis of the CCSS and SJLIFE, 2 novel loci associated with adult BMI among survivors of childhood ALL (LINC00856 rs575792008 and EMR1 rs62123082; PMeta < 5E-8) were identified. It was estimated that the more than 700 known loci explained 6.2% of the variation in adult BMI in childhood ALL survivors. Within the known loci, significant main effects for 23 loci and statistical interactions with CRT at 9 loci (P < 7.0E-5) were further identified. At 2 CRT-interacting loci, DNA methylation patterns may have differed by age. CONCLUSIONS: Adult survivors of childhood ALL have genetic heritability for BMI similar to that observed in the general population. This study provides evidence that treatment with CRT can modify the effect of genetic variants on adult BMI in childhood ALL survivors.

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.

Estimated number of adult survivors of childhood cancer in United States with cancer-predisposing germline variants.

PURPOSE: To estimate the absolute number of adult survivors of childhood cancer in the U.S. population who carry a pathogenic or likely pathogenic variant in a cancer predisposition gene. METHODS: Using the Surveillance, Epidemiology, and End Results (SEER) Program, we estimated the number of childhood cancer survivors on December 31, 2016 for each childhood cancer diagnosis, multiplied this by the proportion of carriers of pathogenic/likely pathogenic variants in the St. Jude Lifetime Cohort (SJLIFE) study, and projected the resulting number onto the U.S. RESULTS: Based on genome sequence data, 11.8% of 2450 SJLIFE participants carry a pathogenic/likely pathogenic variant in one of 156 cancer predisposition genes. Given this information, we estimate that 21 800 adult survivors of childhood cancer in the United States carry a pathogenic/likely pathogenic variant in one of these genes. The highest estimated absolute number of variant carriers are among survivors of central nervous system tumors (n = 4300), particularly astrocytoma (n = 1800) and other gliomas (n = 1700), acute lymphoblastic leukemia (n = 4300), and retinoblastoma (n = 3500). The most frequently mutated genes are RB1 (n = 3000), NF1 (n = 2300), and BRCA2 (n = 800). CONCLUSION: Given the increasing number of childhood cancer survivors in the United States, clinicians should counsel survivors regarding their potential genetic risk, consider referral for genetic counseling and testing, and, as appropriate, implement syndrome-specific cancer surveillance or risk-reducing measures.

Endometriosis risk alleles at 1p36.12 act through inverse regulation of CDC42 and LINC00339.

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

The genetic regulation of transcription in human endometrial tissue.

Study question: Do genetic effects regulate gene expression in human endometrium? Summary answer: This study demonstrated strong genetic effects on endometrial gene expression and some evidence for genetic regulation of gene expression in a menstrual cycle stage-specific manner. What is known already: Genetic effects on expression levels for many genes are tissue specific. Endometrial gene expression varies across menstrual cycle stages and between individuals, but there are limited data on genetic control of expression in endometrium. Study design, size, duration: We analysed genome-wide genotype and gene expression data to map cis expression quantitative trait loci (eQTL) in endometrium. Participants/materials, setting, methods: We recruited 123 women of European ancestry. DNA samples from blood were genotyped on Illumina HumanCoreExome chips. Total RNA was extracted from endometrial tissues. Whole-transcriptome profiles were characterized using Illumina Human HT-12 v4.0 Expression Beadchips. We performed eQTL mapping with ~8 000 000 genotyped and imputed single nucleotide polymorphisms (SNPs) and 12 329 genes. Main results and the role of chance: We identified a total of 18 595 cis SNP-probe associations at a study-wide level of significance (P < 1 × 10-7), which correspond to independent eQTLs for 198 unique genes. The eQTLs with the largest effect in endometrial tissue were rs4902335 for CHURC1 (P = 1.05 × 10-32) and rs147253019 for ZP3 (P = 8.22 × 10-30). We further performed a context-specific eQTL analysis to investigate if genetic effects on gene expression regulation act in a menstrual cycle-specific manner. Interestingly, five cis-eQTLs were identified with a significant stage-by-genotype interaction. The strongest stage interaction was the eQTL for C10ORF33 (PYROXD2) with SNP rs2296438 (P = 2.0 × 10-4), where we observe a 2-fold difference in the average expression levels of heterozygous samples depending on the stage of the menstrual cycle. Large scale data: The summary eQTL results are publicly available to browse or download. Limitations, reasons for caution: A limitation of the present study was the relatively modest sample size. It was not powered to identify trans-eQTLs and larger sample sizes will also be needed to provide better power to detect cis-eQTLs and cycle stage-specific effects, given the substantial changes in expression across the menstrual cycle for many genes. Wider implications of the findings: Identification of endometrial eQTLs provides a platform for better understanding genetic effects on endometriosis risk and other endometrial-related pathologies. Study funding/competing interest(s): Funding for this work was provided by NHMRC Project Grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321 and APP1083405. There are no competing interests.

A Genome-Wide Association Study to Identify Potential Germline Copy Number Variants for Sporadic Breast Cancer Susceptibility.

Breast cancer (BC) predisposition in populations arises from both genetic and nongenetic risk factors. Structural variations such as copy number variations (CNVs) are heritable determinants for disease susceptibility. The primary objectives of this study are (1) to identify CNVs associated with sporadic BC using a genome-wide association study (GWAS) design; (2) to utilize 2 distinct CNV calling algorithms to identify concordant CNVs as a strategy to reduce false positive associations in the hypothesis-generating GWAS discovery phase, and (3) to identify potential candidate CNVs for follow-up replication studies. We used Affymetrix SNP Array 6.0 data profiled on Caucasian subjects (422 cases/348 controls) to call CNVs using algorithms implemented in Nexus Copy Number and Partek Genomics Suite software. Nexus algorithm identified CNVs associated with BC (731 autosomal CNVs with >5% frequency in the total sample and Q < 0.05). Thirteen CNVs were identified when Partek algorithm-called CNVs were overlapped with Nexus-identified CNVs; these CNVs showed concordances for frequency, effect size, and direction. Coding genes present within BC-associated CNVs were known to play a role in disease etiology and prognosis. Long noncoding RNAs identified within CNVs showed tissue-specific expression, indicating potential functional relevance of the findings. The identified candidate CNVs warrant independent replication.

Endometrial vezatin and its association with endometriosis risk.

STUDY QUESTION: Do endometriosis risk-associated single nucleotide polymorphisms (SNPs) found at the 12q22 locus have effects on vezatin ( ITALIC! VEZT) expression? SUMMARY ANSWER: The original genome-wide association study (GWAS) SNP (rs10859871), and other newly identified association signals, demonstrate strong evidence for ITALIC! cis-expression quantitative trait loci (eQTL) effects on ITALIC! VEZT expression. WHAT IS KNOWN ALREADY: GWAS have identified several disease-risk loci (SNPs) associated with endometriosis. The SNP rs10859871 is located within the ITALIC! VEZT gene. ITALIC! VEZT expression is altered in the endometrium of endometriosis patients and is an excellent candidate for having a causal role in endometriosis. Most of the SNPs identified from GWAS are not located within the coding region of the genome. However, they are likely to have an effect on the regulation of gene expression. Genetic variants that affect levels of gene expression are called expression quantitative trait loci (eQTL). STUDY DESIGN, SIZE, DURATION: Samples for genotyping and ITALIC! VEZT variant screening were drawn from women recruited for genetic studies in Australia/New Zealand and women undergoing surgery in a tertiary care centre. Coding variants for ITALIC! VEZT were screened in blood from 100 unrelated individuals (endometriosis-dense families) from the QIMR Berghofer Medical Research Institute dataset. SNPs at the 12q22 locus were imputed and reanalysed for their association with endometriosis. Reanalysis of endometriosis risk-association was performed on a final combined Australian dataset of 2594 cases and 4496 controls. Gene expression was performed on 136 endometrial samples. eQTL analysis in whole blood was performed on 862 individuals from the Brisbane Systems Genetics Study. Endometrial tissue-specific eQTL analysis was performed on 122 samples (eutopic endometrium) collected following laparoscopic surgery. VEZT protein expression studies employed ITALIC! n = 56 (western blotting) and ITALIC! n = 42 (immunohistochemistry) endometrial samples. PARTICIPANTS/MATERIALS, SETTING, METHODS: The women recruited for this study provided blood and/or endometrial tissue samples in a hospital setting. Genomic DNA was screened for common and coding variants. SNPs of interest in the 12q22 region were genotyped using Agena MassARRAY technology or Taqman SNP genotyping assay. Gene expression profiles from RNA extracted from blood and endometrial tissue samples were generated using Illumina whole-genome expression chips (Human HT-12 v4.0). Whole protein extracted from endometrium was used for VEZT western blots, and paraffin sections of endometrium were employed for VEZT immunohistochemistry semi-quantitative analysis. MAIN RESULTS AND THE ROLE OF CHANCE: A total of 11 coding variants of ITALIC! VEZT (including one novel variant) were identified from an endometriosis-dense cohort. Polymorphic coding and imputed SNPs were combined with previous GWAS data to reanalyse the endometriosis risk association of the 12q22 region. The disease association signal at 12q22 was due to coding variants in ITALIC! VEZT or ITALIC! FGD6 (FYVE, RhoGEF and PH domain-containing 6) and SNPs with the strongest signals were either intronic or intergenic. We found strong evidence for ITALIC! VEZT cis-eQTLs with the sentinel SNP (rs10859871) in blood and endometrium, where the endometriosis risk allele (C) was associated with an increase in ITALIC! VEZT expression. We could not demonstrate this genotype-specific effect on VEZT protein expression in endometrium. However, we did observe a menstrual cycle stage specific increase in VEZT protein expression in endometrial glands, specific to the secretory phase ( ITALIC! P = 2.0 × 10(-4)). LIMITATIONS, REASONS FOR CAUTION: In comparison to the blood sample datasets, the study numbers of endometrial tissues were substantially reduced. Protein studies failed to complement RNA results, also likely a reflection of the low study numbers in these experiments. ITALIC! In silico prediction tools used in this investigation are typically based on cell lines different to our tissues of interest, thus any functional annotations drawn from these approaches should be considered carefully. Therefore, functional studies on VEZT and related pathway components are still warranted to unequivocally implicate a causal role for VEZT in endometriosis pathophysiology. WIDER IMPLICATIONS OF THE FINDINGS: GWAS have proven to be very valuable tools for deciphering complex diseases. Endometriosis is a text-book example of a complex disease, involving genetic, lifestyle and environmental influences. Our focused investigation of the 12q22 region validates an association with increased endometriosis risk. Endometriosis risk SNPs (including rs10859871) located within this locus demonstrated evidence for ITALIC! cis-eQTLs on ITALIC! VEZT expression. By examining women who possess an enhanced genetic risk of developing endometriosis, we have identified an effect on ITALIC! VEZT expression and therefore a potential gene/gene pathway in endometriosis disease establishment and development. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants GNT1012245, GNT1026033, GNT1049472 and GNT1046880. G.W.M. is supported by the NHMRC Fellowship scheme (GNT1078399). S.J.H.-C. is supported by the J.N. Peters Bequest Fellowship. The authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Genetic burden associated with varying degrees of disease severity in endometriosis.

Endometriosis is primarily characterized by the presence of tissue resembling endometrium outside the uterine cavity and is usually diagnosed by laparoscopy. The most commonly used classification of disease, the revised American Fertility Society (rAFS) system to grade endometriosis into different stages based on disease severity (I to IV), has been questioned as it does not correlate well with underlying symptoms, posing issues in diagnosis and choice of treatment. Using two independent European genome-wide association (GWA) datasets and top-level classification of the endometriosis cases based on rAFS [minimal or mild (Stage A) and moderate-to-severe (Stage B) disease], we previously showed that Stage B endometriosis has greater contribution of common genetic variation to its aetiology than Stage A disease. Herein, we extend our previous analysis to four endometriosis stages [minimal (Stage I), mild (Stage II), moderate (Stage III) and severe (Stage IV) disease] based on the rAFS classification system and compared the genetic burden across stages. Our results indicate that genetic burden increases from minimal to severe endometriosis. For the minimal disease, genetic factors may contribute to a lesser extent than other disease categories. Mild and moderate endometriosis appeared genetically similar, making it difficult to tease them apart. Consistent with our previous reports, moderate and severe endometriosis showed greater genetic burden than minimal or mild disease. Overall, our results provide new insights into the genetic architecture of endometriosis and further investigation in larger samples may help to understand better the aetiology of varying degrees of endometriosis, enabling improved diagnostic and treatment modalities.

Functional evaluation of genetic variants associated with endometriosis near GREB1.

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.

Association between endometriosis and the interleukin 1A (IL1A) locus.

STUDY QUESTION: Are single-nucleotide polymorphisms (SNPs) at the interleukin 1A (IL1A) gene locus associated with endometriosis risk? SUMMARY ANSWER: We found evidence for strong association between IL1A SNPs and endometriosis risk. WHAT IS KNOWN ALREADY: Genetic factors contribute substantially to the complex aetiology of endometriosis and the disease has an estimated heritability of ∼51%. We, and others, have conducted genome-wide association (GWA) studies for endometriosis, which identified a total of nine independent risk loci. Recently, two small Japanese studies reported eight SNPs (rs6542095, rs11677416, rs3783550, rs3783525, rs3783553, rs2856836, rs1304037 and rs17561) at the IL1A gene locus as suggestively associated with endometriosis risk. There is also evidence of a link between inflammation and endometriosis. STUDY DESIGN, SIZE, DURATION: We sought to further investigate the eight IL1A SNPs for association with endometriosis using an independent sample of 3908 endometriosis cases and 8568 controls of European and Japanese ancestry. The study was conducted between October 2013 and July 2014. PARTICIPANTS/MATERIALS, SETTING, METHODS: By leveraging GWA data from our previous multi-ethnic GWA meta-analysis for endometriosis, we imputed variants in the IL1A region, using a recent 1000 Genomes reference panel. After combining summary statistics for the eight SNPs from our European and Japanese imputed data with the published results, a fixed-effect meta-analysis was performed. An additional meta-analysis restricted to endometriosis cases with moderate-to-severe (revised American Fertility Society stage 3 or 4) disease versus controls was also performed. MAIN RESULTS AND THE ROLE OF CHANCE: All eight IL1A SNPs successfully replicated at P < 0.014 in the European imputed data with concordant direction and similar size to the effects reported in the original Japanese studies. Of these, three SNPs (rs6542095, rs3783550 and rs3783525) also showed association with endometriosis at a nominal P < 0.05 in our independent Japanese sample. Fixed-effect meta-analysis of the eight SNPs for moderate-to-severe endometriosis produced a genome-wide significant association for rs6542095 (odds ratio = 1.21; 95% confidence interval = 1.13-1.29; P = 3.43 × 10(-8)). LIMITATIONS, REASONS FOR CAUTION: The meta-analysis for moderate-to-severe endometriosis included results of moderate-to-severe endometriosis cases from our European data sets and all endometriosis cases from the Japanese data sets, as disease stage information was not available for endometriosis cases in the Japanese data sets. WIDER IMPLICATIONS OF THE FINDINGS: SNP rs6542095 is located ∼2.3 kb downstream of the IL1A gene and ∼6.9 kb upstream of cytoskeleton-associated protein 2-like (CKAP2L) gene. The IL1A gene encodes the IL1a protein, a member of the interleukin 1 cytokine family which is involved in various immune responses and inflammatory processes. These results provide important replication in an independent Japanese sample and, for the first time, association of the IL1A locus in endometriosis patients of European ancestry. SNPs within the IL1A locus may regulate other genes, but if IL1A is the target, our results provide supporting evidence for a link between inflammatory responses and the pathogenesis of endometriosis. STUDY FUNDING/COMPETING INTERESTS: The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study.

Independent Replication and Meta-Analysis for Endometriosis Risk Loci.

Endometriosis is a complex disease that affects 6-10% of women in their reproductive years and 20-50% of women with infertility. Genome-wide and candidate-gene association studies for endometriosis have identified 10 independent risk loci, and of these, nine (rs7521902, rs13394619, rs4141819, rs6542095, rs1519761, rs7739264, rs12700667, rs1537377, and rs10859871) are polymorphic in European populations. Here we investigate the replication of nine SNP loci in 998 laparoscopically and histologically confirmed endometriosis cases and 783 disease-free controls from Belgium. SNPs rs7521902, rs13394619, and rs6542095 show nominally significant (p < .05) associations with endometriosis, while the directions of effect for seven SNPs are consistent with the original reports. Association of rs6542095 at the IL1A locus with 'All' (p = .066) and 'Grade_B' (p = .01) endometriosis is noteworthy because this is the first successful replication in an independent population. Meta-analysis with the published results yields genome-wide significant evidence for rs7521902, rs13394619, rs6542095, rs12700667, rs7739264, and rs1537377. Notably, three coding variants in GREB1 (near rs13394619) and CDKN2B-AS1 (near rs1537377) also showed nominally significant associations with endometriosis. Overall, this study provides important replication in a uniquely characterized independent population, and indicates that the majority of the original genome-wide association findings are not due to chance alone.