Inter-individual variation in ovarian reserve after gonadotoxic treatment in female childhood cancer survivors - a genome-wide association study: results from PanCareLIFE.

OBJECTIVE: We aimed to discover new variants associated with low ovarian reserve after gonadotoxic treatment among adult female childhood cancer survivors using a genome-wide association study approach. DESIGN: Genome-wide association study. SUBJECTS: A discovery cohort of adult female childhood cancer survivors, from the pan-European PanCareLIFE cohort (n=743; median age: 25.8 years), excluding those who received bilateral ovarian irradiation, bilateral oophorectomy, central nerve system or total body irradiation, or stem cell transplantation. Replication was attempted in the USA-based St. Jude Lifetime Cohort (n=391; median age: 31.3 years). EXPOSURE: Female childhood cancer survivors are at risk of therapy-related gonadal impairment. Alkylating agents are well-established risk factors, and the inter-individual variability in gonadotoxicity may be explained by genetic polymorphisms. Data were collected in real-life conditions and cyclophosphamide equivalent dose was used to quantify alkylation agent exposure. INTERVENTION: No intervention was performed. MAIN OUTCOME MEASURE: Anti-Müllerian hormone (AMH) levels served as a proxy for ovarian function and findings were combined in a meta-analysis. RESULTS: Three genome-wide significant (<5.0x10-8) and 16 genome-wide suggestive (<5.0x10-6) loci were associated with log-transformed AMH levels, adjusted for cyclophosphamide equivalent dose of alkylating agents, age at diagnosis, and age at study in the PanCareLIFE cohort. Based on effect allele frequency (EAF) (>0.01 if not genome-wide significant), p-value (<5.0×10-6), and biological relevance, 15 SNPs were selected for replication. None of the SNPs were statistically significantly associated with AMH levels. A meta-analysis indicated that rs78861946 was associated at borderline genome-wide statistical significance (Reference/effect allele: C/T; EAF: 0.04, Beta (SE): -0.484 (0.091), p-value= 9.39×10-8). CONCLUSION: This study found no genetic variants associated with a lower ovarian reserve after gonadotoxic treatment, as the findings of this GWAS were not statistically significant replicated in the replication cohort. Suggestive evidence for potential importance of one variant is briefly discussed, but the lack of statistical significance calls for larger cohort sizes. As the population of childhood cancer survivors is increasing, large-scale and systematic research is needed to identify genetic variants that could aid predictive risk models of gonadotoxicity and as well as fertility preservation options for childhood cancer survivors.

Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.

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

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Unraveling Dexamethasone-Induced Neurobehavioral and Sleep Problems in Children With ALL: Which Determinants Are Important?

PURPOSE: Dexamethasone, the preferred corticosteroid in most treatment protocols for pediatric acute lymphoblastic leukemia (ALL), can induce undesirable side effects. Neurobehavioral and sleep problems are frequently reported, but the interpatient variability is high. We therefore aimed to identify determinants for parent-reported dexamethasone-induced neurobehavioral and sleep problems in pediatric ALL. METHODS: Our prospective study included patients with medium-risk ALL and their parents during maintenance treatment. Patients were assessed before and after one 5-day dexamethasone course. Primary end points were parent-reported dexamethasone-induced neurobehavioral and sleep problems, measured with the Strengths and Difficulties Questionnaire and Sleep Disturbance Scale for Children, respectively. Analyzed determinants included patient and parent demographics, disease and treatment characteristics, parenting stress (Parenting Stress Index and Distress Thermometer for Parents), dexamethasone pharmacokinetics, and genetic variation (candidate single-nucleotide polymorphisms rs41423247 and rs4918). Statistically significant determinants identified in univariable logistic regression analyses were incorporated in a multivariable model. RESULTS: We included 105 patients: median age was 5.4 years (range, 3.0-18.8) and 61% were boys. Clinically relevant dexamethasone-induced neurobehavioral and sleep problems were reported by parents in 70 (67%) and 61 (59%) patients, respectively. In our multivariable regression models, we identified parenting stress as a significant determinant for parent-reported neurobehavioral (odds ratio [OR], 1.16; 95% CI, 1.07 to 1.26) and sleep problems (OR, 1.06; 95% CI, 1.02 to 1.10). Furthermore, parents who experienced more stress before start of a dexamethasone course reported more sleep problems in their child (OR, 1.16; 95% CI, 1.02 to 1.32). CONCLUSION: We identified parenting stress, and not dexamethasone pharmacokinetics, genetic variation, patient/parent demographics, or disease/treatment characteristics, as a significant determinant for parent-reported dexamethasone-induced neurobehavioral and sleep problems. Parenting stress may be a modifiable target to reduce these problems.

Risk Assessment for Hip and Knee Osteoarthritis Using Polygenic Risk Scores.

OBJECTIVE: Polygenic risk scores (PRS) allow risk stratification using common single-nucleotide polymorphisms (SNPs), and clinical applications are currently explored for several diseases. This study was undertaken to assess the risk of hip and knee osteoarthritis (OA) using PRS. METHODS: We analyzed 12,732 individuals from a population-based cohort from the Rotterdam Study (n = 11,496), a clinical cohort (Cohort Hip and Cohort Knee [CHECK] study; n = 908), and a high-risk cohort of overweight women (Prevention of Knee OA in Overweight Females [PROOF] study; n = 328), for the association of the PRS with prevalence/incidence of radiographic OA, of clinical OA, and of total hip replacement (THR) or total knee replacement (TKR). The hip PRS and knee PRS contained 44 and 24 independent SNPs, respectively, and were derived from a recent genome-wide association study meta-analysis. Standardized PRS (with Z transformation) were used in all analyses. RESULTS: We found a stronger association of the PRS for clinically defined OA compared to radiographic OA phenotypes, and we observed the highest PRS risk stratification for TKR/THR. The odds ratio (OR) per SD was 1.3 for incident THR (95% confidence interval [95% CI] 1.1-1.5) and 1.6 (95% CI 1.3-1.9) for incident TKR in the Rotterdam Study. The knee PRS was associated with incident clinical knee OA in the CHECK study (OR 1.3 [95% CI 1.1-1.5]), but not for the PROOF study (OR 1.2 [95% CI 0.8-1.7]). The OR for OA increased gradually across the PRS distribution, up to 2.1 (95% CI 1.4-3.2) for individuals with the 10% highest PRS compared to the middle 50% of the PRS distribution. CONCLUSION: Our findings validated the association of PRS across OA definitions. Since OA is becoming frequent and primary prevention is not commonly applicable, PRS-based risk assessment could play a role in OA prevention. However, the utility of PRS is dependent on the setting. Further studies are needed to test the integration of genetic risk assessment in diverse health care settings.

Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

A bioinformatics pipeline for estimating mitochondrial DNA copy number and heteroplasmy levels from whole genome sequencing data.

Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease.

Association between an AMH promoter polymorphism and serum AMH levels in PCOS patients.

STUDY QUESTION: Do polymorphisms in the anti-Müllerian hormone (AMH) promoter have an effect on AMH levels in patients with polycystic ovary syndrome (PCOS)? SUMMARY ANSWER: We have identified a novel AMH promoter polymorphism rs10406324 that is associated with lower serum AMH levels and is suggested to play a role in the mechanism of regulation of AMH gene expression in women. WHAT IS KNOWN ALREADY: Follicle number is positively correlated with serum AMH levels, reflected by elevated AMH levels in women with PCOS. In addition, it is suggested that AMH production per follicle is higher in women with PCOS than in normo-ovulatory women, implying an altered regulation of AMH in PCOS. STUDY DESIGN, SIZE, DURATION: A discovery cohort of 655 PCOS women of Northern European ancestry and both an internal and external validation PCOS cohort (n = 458 and n = 321, respectively) were included in this study. Summary-level data of an AMH genome-wide association study meta-analysis including 7049 normo-ovulatory women was included as a control cohort. A genetic approach was taken through association analysis and in silico analysis of the associated variants in the AMH promoter. In vitro analysis was performed to investigate the functional mechanisms. PARTICIPANTS/MATERIALS, SETTING, METHODS: All common two-allelic single-nucleotide polymorphisms (SNPs) in the region Chr19:2 245 353-2 250 827 bp (Build 37) were selected for the analysis. Linear regression analyses were performed to determine the association between SNPs in the AMH promoter region and serum AMH levels. For the in silico analysis, the webtools 'HaploReg' v4.1 for ENCODE prediction weight matrices and 'atSNP' were used. In vitro analysis was performed using KK1 cells, a mouse granulosa cell line and COV434 cells, a human granulosa tumor cell line. Cells were transfected with the reference or the variant human AMH promoter reporter construct together with several transcription factors (TFs). Dual-Glo® Luciferase Assay was performed to measure the luciferase activity. MAIN RESULTS AND THE ROLE OF CHANCE: Polymorphism rs10406324 was significantly associated with serum AMH levels in all three PCOS cohorts. Carriers of the minor allele G had significantly lower log-transformed serum AMH levels compared to non-carriers (P = 8.58 × 10-8, P = 1.35 × 10-3 and P = 1.24 × 10-3, respectively). This result was validated in a subsequent meta-analysis (P = 3.24 × 10-12). Interestingly, rs10406324 was not associated with follicle count, nor with other clinical traits. Also, in normo-ovulatory women, the minor allele of this variant was associated with lower serum AMH levels (P = 1.04 × 10-5). These findings suggest that polymorphism rs10406324 plays a role in the regulation of AMH expression, irrespective of clinical background. In silico analysis suggested a decreased binding affinity of the TFs steroidogenenic factor 1, estrogen-related receptor alpha and glucocorticoid receptor to the minor allele G variant, however in vitro analysis did not show a difference in promoter activity between the A and G allele. LIMITATIONS, REASONS FOR CAUTION: Functional analyses were performed in a mouse and a human granulosa cell line using an AMH promoter reporter construct. This may have limited assessment of the impact of the polymorphism on higher order chromatin structures. Human granulosa cells generated from induced pluripotent stem cells, combined with gene editing, may provide a method to elucidate the exact mechanism behind the decrease in serum AMH levels in carriers of the -210 G allele. We acknowledge that the lack of follicle number in the external validation and the control cohort is a limitation of the paper. Although we observed that the association between rs10406324 and AMH levels was independent of follicle number in our discovery and internal validation PCOS cohorts, we cannot fully rule out that the observed effects on serum AMH levels are, in part, caused by differences in follicle number. WIDER IMPLICATIONS OF THE FINDINGS: These results suggest that variations in serum AMH levels are not only caused by differences in follicle number but also by genetic factors. Therefore, the genetic context should be taken into consideration when assessing serum AMH levels in women. This may have clinical consequences when serum AMH levels are used as a marker for the polycystic ovarian morphology phenotype. STUDY FUNDING/COMPETING INTEREST(S): No external funding was used. J.S.E.L. has received consultancy fees from the following companies: Ferring, Roche Diagnostics and Ansh Labs and has received travel reimbursement from Ferring. J.A.V. has received royalties from AMH assays, paid to the institute/lab with no personal financial gain. The other authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

89Zr-3,2-HOPO-Mesothelin Antibody PET Imaging Reflects Tumor Uptake of Mesothelin-Targeted 227Th-Conjugate Therapy in Mice.

The mesothelin (MSLN)-targeted 227Th conjugate is a novel α-therapy developed to treat MSLN-overexpressing cancers. We radiolabeled the same antibody-chelator conjugate with 89Zr to evaluate whether PET imaging with 89Zr-MSLN matches 227Th-MSLN tumor uptake, biodistribution, and antitumor activity. Methods: Serial PET imaging with protein doses of 4, 20, or 40 µg of 89Zr-MSLN and 89Zr-control was performed up to 168 h after tracer injection in human tumor-bearing nude mice with high (HT29-MSLN) and low (BxPc3) MSLN expression. 89Zr-MSLN and 227Th-MSLN ex vivo tumor uptake and biodistribution were compared at 6 time points in HT29-MSLN and in medium-MSLN-expressing (OVCAR-3) tumor-bearing mice. 89Zr-MSLN PET imaging was performed before 227Th-MSLN treatment in HT29-MSLN and BxPc3 tumor-bearing mice. Results: 89Zr-MSLN PET imaging showed an SUVmean of 2.2 ± 0.5 in HT29-MSLN tumors. Ex vivo tumor uptake was 10.6% ± 2.4% injected dose per gram at 168 h. 89Zr-MSLN tumor uptake was higher than uptake of 89Zr-control (P = 0.0043). 89Zr-MSLN and 227Th-MSLN showed comparable tumor uptake and biodistribution in OVCAR-3 and HT29-MSLN tumor-bearing mice. Pretreatment SUVmean was 2.2 ± 0.2 in HT29-MSLN tumors, which decreased in volume on 227Th-MSLN treatment. BxPc3 tumors showed an SUVmean of 1.2 ± 0.3 and remained similar in size after 227Th-MSLN treatment. Conclusion: 89Zr-MSLN PET imaging reflected MSLN expression and matched 227Th-MSLN tumor uptake and biodistribution. Our data support the clinical exploration of 89Zr-MSLN PET imaging together with 227Th-MSLN therapy, both using the same antibody-chelator conjugate.

Toll-Like Receptor 1 Locus Re-examined in a Genome-Wide Association Study Update on Anti-Helicobacter pylori IgG Titers.

BACKGROUND & AIMS: A genome-wide significant association between anti-Helicobacter pylori (H pylori) IgG titers and Toll-like receptor (TLR1/6/10) locus on 4p14 was demonstrated for individuals of European ancestry, but not uniformly replicated. We re-investigated this association in an updated genome-wide association study (GWAS) meta-analysis for populations with low gastric cancer incidence, address potential causes of cohort heterogeneity, and explore functional implications of genetic variation at the TLR1/6/10 locus. METHODS: The dichotomous GWAS (25% individuals exhibiting highest anti-H pylori IgG titers vs remaining 75%) included discovery and replication sampls of, respectively, n = 15,685 and n = 9676, all of European ancestry. Longitudinal analysis of serologic data was performed on H pylori-eradicated subjects (n = 132) and patients under surveillance for premalignant gastric lesions (n = 107). TLR1/6/10 surface expression, TLR1 mRNA, and cytokine levels were measured in leukocyte subsets of healthy subjects (n = 26) genotyped for TLR1/6/10 variants. RESULTS: The association of the TLR1/6/10 locus with anti-H pylori IgG titers (rs12233670; ß = -0.267 ± SE 0.034; P = 4.42 × 10-15) presented with high heterogeneity and failed replication. Anti-H pylori IgG titers declined within 2-4 years after eradication treatment (P = 0.004), and decreased over time in patients with premalignant gastric lesions (P < 0.001). Variation at the TLR1/6/10 locus affected TLR1-mediated cytokine production and TLR1 surface expression on monocytes (P = 0.016) and neutrophils (P = 0.030), but not mRNA levels. CONCLUSIONS: The association between anti-H pylori IgG titers and TLR1/6/10 locus was not replicated across cohorts, possibly owing to dependency of anti-H pylori IgG titers on therapy, clearance, and antibody decay. H pylori-mediated immune cell activation is partly mediated via TLR1 signaling, which in turn is affected by genetic variation.

Effect of Genetic Variation in CYP450 on Gonadal Impairment in a European Cohort of Female Childhood Cancer Survivors, Based on a Candidate Gene Approach: Results from the PanCareLIFE Study.

BACKGROUND: Female childhood cancer survivors (CCSs) carry a risk of therapy-related gonadal dysfunction. Alkylating agents (AA) are well-established risk factors, yet inter-individual variability in ovarian function is observed. Polymorphisms in CYP450 enzymes may explain this variability in AA-induced ovarian damage. We aimed to evaluate associations between previously identified genetic polymorphisms in CYP450 enzymes and AA-related ovarian function among adult CCSs. METHODS: Anti-Müllerian hormone (AMH) levels served as a proxy for ovarian function in a discovery cohort of adult female CCSs, from the pan-European PanCareLIFE cohort (n = 743; age (years): median 25.8, interquartile range (IQR) 22.1-30.6). Using two additive genetic models in linear and logistic regression, nine genetic variants in three CYP450 enzymes were analyzed in relation to cyclophosphamide equivalent dose (CED) score and their impact on AMH levels. The main model evaluated the effect of the variant on AMH and the interaction model evaluated the modifying effect of the variant on the impact of CED score on log-transformed AMH levels. Results were validated, and meta-analysis performed, using the USA-based St. Jude Lifetime Cohort (n = 391; age (years): median 31.3, IQR 26.6-37.4). RESULTS: CYP3A4*3 was significantly associated with AMH levels in the discovery and replication cohort. Meta-analysis revealed a significant main deleterious effect (Beta (95% CI): -0.706 (-1.11--0.298), p-value = 7 × 10-4) of CYP3A4*3 (rs4986910) on log-transformed AMH levels. CYP2B6*2 (rs8192709) showed a significant protective interaction effect (Beta (95% CI): 0.527 (0.126-0.928), p-value = 0.01) on log-transformed AMH levels in CCSs receiving more than 8000 mg/m2 CED. CONCLUSIONS: Female CCSs CYP3A4*3 carriers had significantly lower AMH levels, and CYP2B6*2 may have a protective effect on AMH levels. Identification of risk-contributing variants may improve individualized counselling regarding the treatment-related risk of infertility and fertility preservation options.

Genome-wide association studies identify 137 genetic loci for DNA methylation biomarkers of aging.

BACKGROUND: Biological aging estimators derived from DNA methylation data are heritable and correlate with morbidity and mortality. Consequently, identification of genetic and environmental contributors to the variation in these measures in populations has become a major goal in the field. RESULTS: Leveraging DNA methylation and SNP data from more than 40,000 individuals, we identify 137 genome-wide significant loci, of which 113 are novel, from genome-wide association study (GWAS) meta-analyses of four epigenetic clocks and epigenetic surrogate markers for granulocyte proportions and plasminogen activator inhibitor 1 levels, respectively. We find evidence for shared genetic loci associated with the Horvath clock and expression of transcripts encoding genes linked to lipid metabolism and immune function. Notably, these loci are independent of those reported to regulate DNA methylation levels at constituent clock CpGs. A polygenic score for GrimAge acceleration showed strong associations with adiposity-related traits, educational attainment, parental longevity, and C-reactive protein levels. CONCLUSION: This study illuminates the genetic architecture underlying epigenetic aging and its shared genetic contributions with lifestyle factors and longevity.

A comparison of genotyping arrays.

Array technology to genotype single-nucleotide variants (SNVs) is widely used in genome-wide association studies (GWAS), clinical diagnostics, and linkage studies. Arrays have undergone a tremendous growth in both number and content over recent years making a comprehensive comparison all the more important. We have compared 28 genotyping arrays on their overall content, genome-wide coverage, imputation quality, presence of known GWAS loci, mtDNA variants and clinically relevant genes (i.e., American College of Medical Genetics (ACMG) actionable genes, pharmacogenetic genes, human leukocyte antigen (HLA) genes and SNV density). Our comparison shows that genome-wide coverage is highly correlated with the number of SNVs on the array but does not correlate with imputation quality, which is the main determinant of GWAS usability. Average imputation quality for all tested arrays was similar for European and African populations, indicating that this is not a good criterion for choosing a genotyping array. Rather, the additional content on the array, such as pharmacogenetics or HLA variants, should be the deciding factor. As the research question of a study will in large part determine which class of genes are of interest, there is not just one perfect array for all different research questions. This study can thus help as a guideline to determine which array best suits a study's requirements.

Possible modification of BRSK1 on the risk of alkylating chemotherapy-related reduced ovarian function.

STUDY QUESTION: Do genetic variations in the DNA damage response pathway modify the adverse effect of alkylating agents on ovarian function in female childhood cancer survivors (CCS)? SUMMARY ANSWER: Female CCS carrying a common BR serine/threonine kinase 1 (BRSK1) gene variant appear to be at 2.5-fold increased odds of reduced ovarian function after treatment with high doses of alkylating chemotherapy. WHAT IS KNOWN ALREADY: Female CCS show large inter-individual variability in the impact of DNA-damaging alkylating chemotherapy, given as treatment of childhood cancer, on adult ovarian function. Genetic variants in DNA repair genes affecting ovarian function might explain this variability. STUDY DESIGN, SIZE, DURATION: CCS for the discovery cohort were identified from the Dutch Childhood Oncology Group (DCOG) LATER VEVO-study, a multi-centre retrospective cohort study evaluating fertility, ovarian reserve and risk of premature menopause among adult female 5-year survivors of childhood cancer. Female 5-year CCS, diagnosed with cancer and treated with chemotherapy before the age of 25 years, and aged 18 years or older at time of study were enrolled in the current study. Results from the discovery Dutch DCOG-LATER VEVO cohort (n = 285) were validated in the pan-European PanCareLIFE (n = 465) and the USA-based St. Jude Lifetime Cohort (n = 391). PARTICIPANTS/MATERIALS, SETTING, METHODS: To evaluate ovarian function, anti-Müllerian hormone (AMH) levels were assessed in both the discovery cohort and the replication cohorts. Using additive genetic models in linear and logistic regression, five genetic variants involved in DNA damage response were analysed in relation to cyclophosphamide equivalent dose (CED) score and their impact on ovarian function. Results were then examined using fixed-effect meta-analysis. MAIN RESULTS AND THE ROLE OF CHANCE: Meta-analysis across the three independent cohorts showed a significant interaction effect (P = 3.0 × 10-4) between rs11668344 of BRSK1 (allele frequency = 0.34) among CCS treated with high-dose alkylating agents (CED score ≥8000 mg/m2), resulting in a 2.5-fold increased odds of a reduced ovarian function (lowest AMH tertile) for CCS carrying one G allele compared to CCS without this allele (odds ratio genotype AA: 2.01 vs AG: 5.00). LIMITATIONS, REASONS FOR CAUTION: While low AMH levels can also identify poor responders in assisted reproductive technology, it needs to be emphasized that AMH remains a surrogate marker of ovarian function. WIDER IMPLICATIONS OF THE FINDINGS: Further research, validating our findings and identifying additional risk-contributing genetic variants, may enable individualized counselling regarding treatment-related risks and necessity of fertility preservation procedures in girls with cancer. STUDY FUNDING/COMPETING INTEREST(S): This work was supported by the PanCareLIFE project that has received funding from the European Union's Seventh Framework Programme for research, technological development and demonstration under grant agreement no 602030. In addition, the DCOG-LATER VEVO study was funded by the Dutch Cancer Society (Grant no. VU 2006-3622) and by the Children Cancer Free Foundation (Project no. 20) and the St Jude Lifetime cohort study by NCI U01 CA195547. The authors declare no competing interests. TRIAL REGISTRATION NUMBER: N/A.

Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene.

BACKGROUND: Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology. METHODS: We used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees. RESULTS: We found two novel genome-wide significant loci for OA in the thumb joints. We identified WNT9A as a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints: TGFa, RUNX2, COL27A1, ASTN2, IL11 and GDF5 loci. CONCLUSIONS: We identified a robust novel genetic locus for hand OA on chromosome 1, of which WNT9A is the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes.

Short-Term, Combined Fasting and Exercise Improves Body Composition in Healthy Males.

Fasting enhances the beneficial metabolic outcomes of exercise; however, it is unknown whether body composition is favorably modified on the short term. A baseline-follow-up study was carried out to assess the effect of an established protocol involving short-term combined exercise with fasting on body composition. One hundred seven recreationally exercising males underwent a 10-day intervention across 15 fitness centers in the Netherlands involving a 3-day gradual decrease of food intake, a 3-day period with extremely low caloric intake, and a gradual 4-day increase to initial caloric intake, with daily 30-min submaximal cycling. Using dual-energy X-ray absorptiometry analysis, all subjects substantially lost total body mass (-3.9 ± 1.9 kg; p < .001) and fat mass (-3.3 ± 1.3 kg; p < .001). Average lean mass was lost (-0.6 ± 1.5 kg; p < .001), but lean mass as a percentage of total body mass was not reduced. The authors observed a loss of -3.9 ± 1.9% android fat over total fat mass (p < .001), a loss of -2.2 ± 1.9% gynoid over total fat mass (p < .001), and reduced android/gynoid ratios (-0.05 ± 0.1; p < .001). Analyzing 15 preselected single-nucleotide polymorphisms in 13 metabolism-related genes revealed trending associations for thyroid state-related single-nucleotide polymorphisms rs225014 (deiodinase 2) and rs35767 (insulin-like growth factor1), and rs1053049 (PPARD). In conclusion, a short period of combined fasting and exercise leads to a substantial loss of body and fat mass without a loss of lean mass as a percentage of total mass.

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity.

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.

Usefulness of current candidate genetic markers to identify childhood cancer patients at risk for platinum-induced ototoxicity: Results of the European PanCareLIFE cohort study.

BACKGROUND: Irreversible sensorineural hearing loss is a common side effect of platinum treatment with the potential to significantly impair the neurocognitive, social and educational development of childhood cancer survivors. Genetic association studies suggest a genetic predisposition for cisplatin-induced ototoxicity. Among other candidate genes, thiopurine methyltransferase (TPMT) is considered a critical gene for susceptibility to cisplatin-induced hearing loss in a pharmacogenetic guideline. The aim of this cross-sectional cohort study was to confirm the genetic associations in a large pan-European population and to evaluate the diagnostic accuracy of the genetic markers. METHODS: Eligibility criteria required patients to be aged less than 19 years at the start of chemotherapy, which had to include cisplatin and/or carboplatin. Patients were assigned to three phenotype categories: no, minor and clinically relevant hearing loss. Fourteen variants in eleven candidate genes (ABCC3, OTOS, TPMT, SLC22A2, NFE2L2, SLC16A5, LRP2, GSTP1, SOD2, WFS1 and ACYP2) were investigated. Multinomial logistic regression was performed to model the relationship between genetic predictors and platinum ototoxicity, adjusting for clinical risk factors. Additionally, measures of the diagnostic accuracy of the genetic markers were determined. RESULTS: 900 patients were included in this study. In the multinomial logistic regression, significant unique contributions were found from SLC22A2 rs316019, the age at the start of platinum treatment, cranial radiation and the interaction term [platinum compound]∗[cumulative dose of cisplatin]. The predictive performance of the genetic markers was poor compared with the clinical risk factors. CONCLUSIONS: PanCareLIFE is the largest study of cisplatin-induced ototoxicity to date and confirmed a role for the polyspecific organic cation transporter SLC22A2. However, the predictive value of the current genetic candidate markers for clinical use is negligible, which puts the value of clinical factors for risk assessment of cisplatin-induced ototoxicity back into the foreground.

Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.