Melanoma as a detection-dependent cancer: proxy measures of screening and surveillance associated with 3-year incidence of detection procedures for skin cancer in a national prospective study.

BACKGROUND: Increased levels of skin detection activities are suspected of driving recent rapid increases in melanoma incidence. While ecological studies report rising rates of skin biopsies in many industrialised populations, few studies have accessed individual-level clinical data to test the hypothesis that people exposed to skin screening and surveillance go on to experience higher levels of biopsies and excisions, and higher incidence of melanoma. OBJECTIVE: To measure the incidence rates of detection events and melanoma diagnoses in a national cohort during follow-up, stratified according to screening or surveillance activities in the index year. METHODS: We conducted a prospective cohort study of 10,674,200 men and women using linked data from Australia's universal health insurance scheme (2011-2016). Participants who underwent biopsies or who were diagnosed in skin cancer in 2011/2012 were excluded. We then defined participants as surveilled, screened, or unscreened for skin cancer based on medical claims occurring in 2013. We compared rates of biopsies, excisions for suspected melanoma, and excisions for histologically confirmed melanoma during follow-up (2014-2016). RESULTS: After adjusting for socio-demographic factors, screened people were significantly more likely than unscreened people to undergo skin biopsies (RR 2.59, 95% CI 2.56-2.63), excisions for suspected melanoma (RR 2.57, 95%CI 2.53-2.60), and excisions for confirmed melanoma (hazard ratio 3.32, 95% CI 3.24-3.40) during follow-up. Similarly elevated rates of subsequent diagnostic events were observed for the surveilled group. Importantly, rates of detection and rates of melanoma remained elevated in each of the subsequent years of follow-up among those who were screened or surveilled in the index year. CONCLUSIONS: People undergoing skin screening or surveillance subsequently experience higher rates of diagnostic scrutiny and higher rates of melanoma.

No evidence that retinol is protective for skin cancer.

With over 1.5 million new cases annually, skin cancers are the most commonly diagnosed group of cancers worldwide. Among these, melanoma and keratinocyte cancers (KC), comprising squamous cell carcinoma (SCC) and basal cell carcinoma (BCC), are predominant. Retinol, a vitamin A derivative, is essential in the regulation of growth and differentiation of epidermal cells. Moreover, retinol exhibits antioxidant properties, protecting the skin against ultra-violet (UV) radiation induced oxidative damage. Existing research on the impact of retinol on melanoma, SCC and BCC development shows mixed results. Several dietary intake studies have suggested that higher retinol levels reduce skin cancer risk, however, others have failed to find this association. We used two-sample Mendelian randomization (MR) to explore if there is a causal relationship between retinol and the risk of developing melanoma, SCC or BCC. Genetically predicted circulating retinol levels were obtained from a genome wide association study (GWAS) meta-analysis of the INTERVAL (N=11,132) and METSIM (N=6,136) cohorts. Melanoma (30,134 cases and 375,188 controls), SCC (10,557 cases and 537,850 controls) and BCC (36,479 cases and 540,185 controls) risks were derived from published GWAS meta-analyses. We conducted two MR approaches. In the first MR we used a single SNP (rs10882283) that is associated with the levels of Retinol Binding Protein 4 (RBP4) as an instrument variable (IV) for circulating retinol levels. In the second MR we used all independent genetic variants that were strongly associated (P < 5 × 10-8) with retinol levels as IVs. Odds ratios (OR) for skin cancer were calculated for a one standard deviation (SD) increase in genetically predicted retinol levels. The single IV approach revealed that retinol levels were not significantly associated with risk of melanoma (OR = 1.04 [95% confidence interval 0.83, 1.31], P = 0.72), SCC (OR = 1.15 [0.87, 1.51], P = 0.32) or BCC (OR = 1.06 [0.90, 1.23], P = 0.50). Similar null results were observed with the multiple IV approach for melanoma (OR = 1.03 [0.95, 1.11], P = 0.54), SCC (OR = 1.01 [0.91, 1.13], P = 0.83), and BCC (OR = 1.04 [0.96, 1.12], P = 0.38). In conclusion, we found no evidence that circulating retinol levels were causally associated with the development of melanoma, SCC and BCC.

Changes in the incidence of melanoma in Australia, 2006-2021, by age group and ancestry: a modelling study.

OBJECTIVES: To estimate the incidence of melanoma in Australia among people with ancestries associated with low, moderate, or high risk of melanoma, by sex and 5-year age group; to establish whether age-specific incidence rates by ancestry risk group have changed over time. STUDY DESIGN: Modelling study; United States (SEER database) melanoma incidence rates for representative ancestral populations and Australian census data (2006, 2011, 2016, 2021) used to estimate Australian melanoma incidence rates by ancestry-based risk. SETTING, PARTICIPANTS: Australia, 2006-2021. MAIN OUTCOME MEASURES: Age-specific invasive melanoma incidence rates, and average annual percentage change (AAPC) in age-specific melanoma rates, by ancestry-based risk group, sex, and 5-year age group. RESULTS: The proportion of people in Australia who reported high risk (European) ancestry declined from 85.3% in 2006 to 71.1% in 2021. The estimated age-standardised melanoma incidence rate was higher for people with high risk ancestry (2021: males, 82.2 [95% confidence interval {CI}, 80.5-83.8] cases per 100 000 population; females, 58.5 [95% CI, 57.0-59.9] cases per 100 000 population) than for all Australians (males, 67.8 [95% CI, 66.5-69.2] cases per 100 000 population; females, 45.4 [95% CI, 44.3-46.5] cases per 100 000 population). AAPCs were consistently positive for Australians aged 50 years or older, both overall and for people with high risk ancestry, but were statistically significant only for some age groups beyond 65 years. AAPCs were negative for people aged 34 years or younger, but were generally not statistically significant. CONCLUSIONS: Melanoma incidence has declined in some younger age groups in Australia, including among people with high risk ancestry. Social and behavioural changes over the same period that lead to lower levels of ultraviolet radiation exposure probably contributed to these changes.

Exploring the Germline Genetics of In Situ and Invasive Cutaneous Melanoma: A Genome-Wide Association Study Meta-Analysis.

Importance: It is unknown whether germline genetic factors influence in situ melanoma risk differently than invasive melanoma risk. Objective: To determine whether differences in risk of in situ melanoma and invasive melanoma are heritable. Design, Setting, and Participants: Three genome-wide association study meta-analyses were conducted of in situ melanoma vs controls, invasive melanoma vs controls, and in situ vs invasive melanoma (case-case) using 4 population-based genetic cohorts: the UK Biobank, the FinnGen cohort, the QSkin Sun and Health Study, and the Queensland Study of Melanoma: Environmental and Genetic Associations (Q-MEGA). Melanoma status was determined using International Statistical Classification of Diseases and Related Health Problems codes from cancer registry data. Data were collected from 1987 to 2022, and data were analyzed from September 2022 to June 2023. Exposure: In situ and invasive cutaneous melanoma. Main Outcomes and Measures: To test whether in situ and invasive melanoma have independent heritable components, genetic effect estimates were calculated for single-nucleotide variants (SNV; formerly single-nucleotide polymorphisms) throughout the genome for each melanoma. Then, SNV-based heritability was estimated, the genetic correlation between melanoma subtypes was assessed, and polygenic risk scores (PRS) were generated for in situ vs invasive status in Q-MEGA participants. Results: A total of 6 genome-wide significant loci associated with in situ melanoma and 18 loci with invasive melanoma were identified. A strong genetic correlation (genetic r = 0.96; 95% CI, 0.76-1.15) was observed between the 2 classifications. Notably, loci near IRF4, KLF4, and HULC had significantly larger effects for in situ melanoma compared with invasive melanoma, while MC1R had a significantly larger effect on invasive melanoma compared with in situ melanoma. Heritability estimates were consistent for both, with in situ melanoma heritability of 6.7% (95% CI, 4.1-9.3) and invasive melanoma heritability of 4.9% (95% CI, 2.8-7.2). Finally, a PRS, derived from comparing invasive melanoma with in situ melanoma genetic risk, was on average significantly higher in participants with invasive melanoma (odds ratio per 1-SD increase in PRS, 1.43; 95% CI, 1.16-1.77). Conclusions and Relevance: There is much shared genetic architecture between in situ melanoma and invasive melanoma. Despite indistinguishable heritability estimates between the melanoma classifications, PRS suggest germline genetics may influence whether a person gets in situ melanoma or invasive melanoma. PRS could potentially help stratify populations based on invasive melanoma risk, informing future screening programs without exacerbating the current burden of melanoma overdiagnosis.

Estimating the magnitude and healthcare costs of melanoma in situ and thin invasive melanoma overdiagnosis in Australia.

BACKGROUND: Research suggests a high proportion of melanoma in situ (MIS) may be overdiagnosed, potentially contributing to overtreatment, patient harm and inflated costs for individuals and healthcare systems. However, Australia-wide estimates of the magnitude of melanoma overdiagnosis are potentially outdated and there has been no estimation of the cost to the healthcare system. OBJECTIVE: To estimate the magnitude and cost of overdiagnosed MIS and thin invasive melanomas in Australia. METHODS: Using two different methods for calculating lifetime risk, we used routinely collected, national-level data to estimate overdiagnosed MIS and thin invasive melanomas (stage IA) in Australia in 2017 and 2021, separately for men and women. We multiplied the number of overdiagnosed melanomas by the estimated annual cost of a MIS or thin invasive melanoma to quantify the financial burden of melanoma overdiagnosis to the Australian healthcare system in the year following diagnosis. RESULTS: We estimated that between 67-70% of MIS were overdiagnosed in 2017, rising to 71-76% in 2021, contributing to between 19,829 (95%CI: 19,553-20,105) and 20,811 (95%CI: 20,528-21,094) overdiagnosed MIS. In 2021, the estimated costs in Australia ranged between $17.7 million (95%CI: $17.4-17.9 million) and $18.6 million (95%CI: $18.3-18.8 million). We estimated that 22-29% of thin invasive melanomas were overdiagnosed in 2017, rising to 28-34% in 2021, contributing to between 2,831 (95%CI: 2,726-2,935) and 3,168 (95%CI: 3,058-3,279) overdiagnosed thin invasive melanomas. In 2021, the estimated costs from thin invasive melanoma overdiagnoses ranged between $2.5 million (95%CI: $2.4-2.6 million) and $2.8 million (95%CI: $2.7-2.9 million). CONCLUSIONS: Melanoma overdiagnosis is a growing clinical and public health problem in Australia, producing significant economic costs in the year following overdiagnosis. Limiting melanoma overdiagnosis may prevent unnecessary healthcare resource use and improve financial sustainability within the Australian healthcare system.

Phenotypic and genotypic risk factors for invasive melanoma by sex and body site.

BACKGROUND: Cutaneous melanoma incidence varies consistently across body sites between men and women, but the underlying causes of the differences remain unclear. To date, no prospective studies have examined risk factors for melanoma separately for men and women according to body site. METHODS: We compared the association between constitutional, genetic and environmental risk factors for invasive melanoma on different body sites separately for men and women in a population-based prospective cohort study of 17,774 men and 21,070 women aged between 40 and 69 years and residents of Queensland, Australia at baseline in 2011. Participants were followed until December 2021.We examined risk factors including hair colour, tanning ability, naevus density, and proxies for high cumulative sun exposure, all self-reported at baseline. We also examined polygenic risk score (PRS) derived from summary statistics from a melanoma genome-wide association study meta-analysis. RESULTS: During a median 10.4 years of follow-up, 455 men and 331 women developed an incident invasive melanoma; the mean age at diagnosis was lower in women than in men (62.6 vs. 65.0, respectively). The most common body site was the trunk in men (45.1%), and the upper (36.8%) and lower limbs (27.4%) in women. High naevus density and proxy measures of high cumulative sun exposure were similarly associated with melanoma at all sites in men and women. In both sexes, high genetic risk was associated with melanoma on all body sites except the head and neck. We observed differences between men and women in the association between PRS and melanoma of the trunk (highest vs. lowest tertile of PRS: HR 2.78, 95% CI 1.64-4.69 for men; 1.55, 95% CI 0.63-3.80 for women), and non-significant but large differences for the lower limbs (HR 5.25, 95% CI 1.80-15.27 for men; 1.75, 95% CI 0.88-3.47 for women). CONCLUSIONS: While there are a number of potential explanations for these findings, this raises the possibility that genetic factors other than those related to pigmentation and naevus phenotypes may play a role in the predilection for melanoma to arise on different sites between the sexes.

Does Sex Matter? Temporal Analyses of Melanoma Trends among Men and Women Suggest Etiologic Heterogeneity.

The incidence and distribution of cutaneous melanoma differ between the sexes, but it is unclear whether these differences have been constant through time or across generations. We compared incidence trends by age, sex, and anatomic site by analyzing long-term melanoma data (1982-2018) in 3 populations residing at high-, moderate-, and low-ambient sun exposure: Queensland, Australia; United States White; and Scotland. We fit age-period-cohort models and compared trends in the male-to-female incidence rate ratio by site and sex. In men, melanoma incidence was always highest on the trunk; in women, incidence was historically highest on limbs, but there have been recent increases in truncal melanoma among females in all populations. The incidence rate ratio showed excess melanoma on the lower limb in females in most age groups in all populations. In contrast, there was a male excess of melanoma on the trunk (from about age 25 years) and head/neck (from about age 40 years), which increased with age. Birth cohort analyses identified turning points in incidence from high to low incidence among recent birth cohorts, which differed by population and site. Changing exposure to UVR is implicated, possibly superimposed upon innate differences between the sexes in site-specific susceptibility.

Comprehensive Sex-Stratified Genetic Analysis of 28 Blood Biomarkers and Depression Reveals a Significant Association between Depression and Low Levels of Total Protein in Females.

Introduction: Major depression (MD) is more common amongst women than men, and MD episodes have been associated with fluctuations in reproductive hormones amongst women. To investigate biological underpinnings of heterogeneity in MD, the associations between depression, stratified by sex and including perinatal depression (PND), and blood biomarkers, using UK Biobank (UKB) data, were evaluated, and extended to include the association of depression with biomarker polygenic scores (PGS), generated as proxy for each biomarker. Method: Using female (N = 39,761) and male (N = 38,821) UKB participants, lifetime MD and PND were tested for association with 28 blood biomarkers. A GWAS was conducted for each biomarker and genetic correlations with depression subgroups were estimated. Using independent data from the Australian Genetics of Depression Study, PGS were constructed for each biomarker, and tested for association with depression status (n [female cases/controls] = 9,006/6,442; n [male cases/controls] = 3,106/6,222). Regions of significant local genetic correlation between depression subgroups and biomarkers highlighted by the PGS analysis were identified. Results: Depression in females was significantly associated with levels of twelve biomarkers, including total protein (OR = 0.90, CI = [0.86, 0.94], p = 3.9 × 10-6) and vitamin D (OR = 0.94, CI = [0.90, 0.97], p = 2.6 × 10-4), and PND with five biomarker levels, also including total protein (OR = 0.88, CI = [0.81, 0.96], p = 4.7 × 10-3). Depression in males was significantly associated with levels of eleven biomarkers. In the independent Australian Genetics of Depression Study, PGS analysis found significant associations for female depression and PND with total protein (female depression: OR = 0.93, CI = [0.88, 0.98], p = 3.6 × 10-3; PND: OR = 0.91, CI = [0.86, 0.96], p = 1.1 × 10-3), as well as with vitamin D (female depression: OR = 0.93, CI = [0.89, 0.97], p = 2.0 × 10-3; PND: OR = 0.92, CI = [0.87, 0.97], p = 1.4 × 10-3). The male depression sample did not report any significant results, and the point estimate of total protein (OR = 0.98, CI = [0.92-1.04], p = 4.7 × 10-1) did not indicate any association. Local genetic correlation analysis highlighted significant genetic correlation between PND and total protein, located in 5q13.3 (rG = 0.68, CI = [0.33, 1.0], p = 3.6 × 10-4). Discussion and Conclusion: Multiple lines of evidence from genetic analysis highlight an association between total serum protein levels and depression in females. Further research involving prospective measurement of total protein and depressive symptoms is warranted.

Genetic predisposition to childhood obesity does not influence the risk of developing skin cancer in adulthood.

The relationship between body mass index (BMI) and melanoma and other skin cancers remains unclear. The objective of this study was to employ the Mendelian randomization (MR) approach to evaluate the effects of genetically predicted childhood adiposity on the risk of developing skin cancer later in life. Two-sample MR analyses were conducted using summary data from genome-wide association study (GWAS) meta-analyses of childhood BMI, melanoma, cutaneous squamous cell carcinoma (cSCC), and basal cell carcinoma (BCC). We used the inverse-variance-weighted (IVW) methods to obtain a pooled estimate across all genetic variants for childhood BMI. We performed multiple sensitivity analyses to evaluate the potential influence of various assumptions on our findings. We found no evidence that genetically predicted childhood BMI was associated with risks of developing melanoma, cSCC, or BCC in adulthood (OR, 95% CI: melanoma: 1.02 (0.93-1.13), cSCC 0.94 (0.79-1.11), BCC 0.97 (0.84-1.12)). Our findings do not support the conclusions from observational studies that childhood BMI is associated with increased risks of melanoma, cSCC, or BCC in adulthood. Intervening on childhood adiposity will not reduce the risk of common skin cancers later in life.

Environmental plastics in the context of UV radiation, climate change, and the Montreal Protocol.

There are close links between solar UV radiation, climate change, and plastic pollution. UV-driven weathering is a key process leading to the degradation of plastics in the environment but also the formation of potentially harmful plastic fragments such as micro- and nanoplastic particles. Estimates of the environmental persistence of plastic pollution, and the formation of fragments, will need to take in account plastic dispersal around the globe, as well as projected UV radiation levels and climate change factors.

Age group-specific changes in keratinocyte cancer treatment rates in Australia, 2012-2021: a retrospective cohort study based on MBS claims data.

OBJECTIVES: To examine recent changes in the numbers of Medicare-subsidised keratinocyte cancer excisions, particularly for younger people exposed to primary prevention campaigns since the early 1980s. STUDY DESIGN: Retrospective cohort study; analysis of administrative data. SETTING, PARTICIPANTS: Analysis of Medicare Benefits Schedule (MBS) claims data for procedures related to the diagnosis and treatment of keratinocyte cancer in Australia, 2012-2021. MAIN OUTCOME MEASURES: Age-standardised rates for MBS-subsidised claims for first surgical squamous cell carcinoma (SCC) and basal cell carcinoma (BCC) excisions, Mohs surgery, surgical excisions of benign lesions, skin biopsies, and cryotherapy or serial curettage of premalignant and malignant lesions, overall, and by sex, state/territory, and age group; average annual percentage change in rate for time intervals determined by joinpoint regression. RESULTS: In men, the age-standardised rate of BCC/SCC excisions increased by 1.9% (95% confidence interval [CI], 1.4-2.4%) per year during 2012-2019 (from 2931 to 3371 per 100 000 men) and then declined by 3.8% (95% CI, 0.5-7.0%) per year during 2019-2021 (to 3152 per 100 000). In women, the age-standardised rate increased by 2.2% (95% CI, 1.7-2.8%) per year during 2012-2019 (from 1798 to 2093 per 100 000 women); the decline to 1967 excisions per 100 000 women in 2021 was not statistically significant. BCC/SCC excision rates declined for men under 55 years of age (by 1.0-3.4% per year) and women under 45 years of age (by 1.7-2.3% per year). Age-standardised biopsy rates increased during 2012-2021 in all age groups (by 2.8-6.9% per year). CONCLUSIONS: Rates of MBS-subsidised treatment for keratinocyte cancers increased during 2012-2019, but BCC/SCC treatment rates declined among younger Australians, who have probably been exposed to less sunlight than earlier generations because of public health interventions and population-wide lifestyle changes related to technology use.

Genetic Analysis of Perceived Youthfulness Reveals Differences in How Men's and Women's Age Is Assessed.

Skin aging is a natural process that occurs over time but can be accelerated by sun exposure. Measuring skin age in a large population can provide insight into the extent of skin damage from sun exposure and skin cancer risk. Understanding the genetics of skin aging, within and across sexes (males and females), could improve our understanding of the genetic drivers of both skin aging and skin cancer. We used UK Biobank data to examine the genetic overlap between perceived youthfulness and traits relevant to actinic photoaging. Our GWAS identified 22 genome-wide significant loci for women and 43 for men. The genetic correlation (rg) between perceived youthfulness in men and women was significantly less than unity (rg = 0.75, 95% confidence interval = 0.69-0.80), suggesting a gene-by-sex interaction. In women, perceived youthfulness was modestly correlated with keratinocyte cancer (rg = -0.19) and skin tanning (rg = 0.18). In men, perceived youthfulness was correlated with male-pattern baldness (rg = -0.23). This suggests that the genetic architecture of perceived youthfulness may differ between sexes, with genes influencing skin tanning and skin cancer susceptibility driving the difference in women, whereas genes influencing male-pattern baldness and other puberty-related traits drive the difference in men. We recommend that future genetic analysis of skin aging include a sex-stratified component.

Plastics in the environment in the context of UV radiation, climate change and the Montreal Protocol: UNEP Environmental Effects Assessment Panel, Update 2023.

This Assessment Update by the Environmental Effects Assessment Panel (EEAP) of the United Nations Environment Programme (UNEP) considers the interactive effects of solar UV radiation, global warming, and other weathering factors on plastics. The Assessment illustrates the significance of solar UV radiation in decreasing the durability of plastic materials, degradation of plastic debris, formation of micro- and nanoplastic particles and accompanying leaching of potential toxic compounds. Micro- and nanoplastics have been found in all ecosystems, the atmosphere, and in humans. While the potential biological risks are not yet well-established, the widespread and increasing occurrence of plastic pollution is reason for continuing research and monitoring. Plastic debris persists after its intended life in soils, water bodies and the atmosphere as well as in living organisms. To counteract accumulation of plastics in the environment, the lifetime of novel plastics or plastic alternatives should better match the functional life of products, with eventual breakdown releasing harmless substances to the environment.

International patterns and trends in the incidence of melanoma and cutaneous squamous cell carcinoma, 1989-2020.

BACKGROUND: Cutaneous squamous cell carcinoma (cSCC) and melanoma have different associations with sun exposure. OBJECTIVES: To compare trends in the incidence rates of cSCC and melanoma, to provide insight into changing patterns of exposure to ultraviolet radiation (UVR). METHODS: We compared trends in the incidence of cSCC and melanoma in seven susceptible populations residing at mid-to-high latitudes: Finland, Norway, Sweden, Denmark, Scotland, the Netherlands and Tasmania (Australia). We fitted Joinpoint models to describe trends in age-standardized incidence rates for melanoma and cSCC and calculated the average annual percentage rate of change for the period 1989-2020 (1989-2018 for Tasmania). We calculated the incident rate ratio (IRR) as the ratio of the age-standardized rates (European Standard Population) for cSCC to melanoma and conducted age-period-cohort modelling to compare age, period and cohort effects. RESULTS: The ratio of cSCC-to-melanoma incidence increased with proximity to the equator and over time. In the most recent time period, the incidence of cSCC was higher than the incidence of melanoma for men and women in all seven populations. While the ratio of cSCC-to-melanoma incidence was higher for men vs. women, in most countries the cSCC-to-melanoma IRR increased over time to a greater extent in women than in men. Melanoma incidence was higher among younger people and cSCC incidence was higher among older people; the age at which the incidence of cSCC overtook the incidence of melanoma was progressively younger with proximity to the equator. CONCLUSIONS: Despite concerted international efforts to preserve the ozone layer over the past four decades resulting in significant reductions in surface ultraviolet B at mid-latitudes, the incidence of skin cancer, particularly cSCC, continues to rise in those regions. Our findings are consistent with a stronger association with age-associated cumulative sun exposure for cSCC vs. melanoma and suggest that women are currently receiving greater UV radiation exposure than in the past.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

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).