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.

Assessing the genetic risk of nodular melanoma using a candidate gene approach.

BACKGROUND: Nodular melanoma (NM) is a challenge to diagnose early due to its rapid growth and more atypical clinical presentation, making it the largest contributor to melanoma mortality. OBJECTIVES: Our study aim was to perform a rare-variant allele (RVA) analysis of whole-exome sequencing of patients with NM and non-NM (minor allele frequency ≤ 1% non-Finnish European) for a set of 500 candidate genes potentially implicated in melanoma. METHODS: This study recruited 131 participants with NM and 194 with non-NM from South-east Queensland and patients with NM from Victoria to perform a comparative analysis of possible genetic differences or similarities between the two melanoma cohorts. RESULTS: Phenotypic analysis revealed that a majority of patients diagnosed with NM were older males with a higher frequency of fair skin and red hair than is seen in the general population. The distribution of common melanoma polygenic risk scores was similar in patients with NM and non-NM, with over 28% in the highest quantile of scores. There was also a similar frequency of carriage of familial/high-penetrant melanoma gene and loss-of-function variants. We identified 39 genes by filtering 500 candidate genes based on the greatest frequency in NM compared with non-NM cases. The genes with RVAs of greatest frequency in NM included PTCH1, ARID2 and GHR. Rare variants in the SMO gene, which interacts with PTCH1 as ligand and receptor, were also identified, providing evidence that the Hedgehog pathway may contribute to NM risk. There was a cumulative effect in carrying multiple rare variants in the NM-associated genes. A 14.8-fold increased ratio for NM compared with non-NM was seen when two RVAs of the 39 genes were carried by a patient. CONCLUSIONS: This study highlights the importance of considering frequency of RVA to identify those at risk of NM in addition to known high penetrance genes.

Uncovering molecular mechanisms for amelanotic/hypopigmented primary cutaneous melanoma.

BACKGROUND: A portion of approximately 2-20% of cutaneous melanoma (CM) are diagnosed as amelanotic/hypopigmented melanoma (AHM) and represent a challenge for early diagnosis. OBJECTIVES: Since the degree to which somatic mutations and copy number aberrations (CNA) in genes associated with skin-lightening or albinism may contribute to the loss of tumour pigmentation in AHM samples has not yet been addressed, we have investigated loss of function mutations of key pigmentation genes in matched germline and AHM as well as pigmented melanoma (PM) tumour DNA samples. METHODS: An analysis of clinical and histopathological characteristics together with whole exome sequencing data of 34 fresh frozen primary CM, graded according to the amount of pigmentation present was performed. Together with germline and somatic variant analysis, 30 samples were previously analysed for CNA changes. This study focussed on germline and somatic variants in the coding region of 16 genes known to be associated with albinism/hypopigmentation or variation in human pigmentation in all samples. Chromosomal regions encompassing these 16 genes were examined for DNA copy loss or gain. RESULTS: The finding that red hair related MC1R and TYR R402Q loss of activity gene variant alleles and genotypes are associated with AHM was validated in this study. Germline AHM-related gene variants were enriched in 70% (n=7 of 10) of AHM patients vs 8.3% (n=2 of 24) of PM patients. This surprisingly high frequency of rare germline variants in AHM patients constitutes the "first hit" and confirms that AHM patients are more likely to be albinism allele carriers than patients with PM. Next, in CNA analysis of each tumour sample, 50% (n=4 of 8) AHM samples with a pigmentation gene variant had LOH in the region containing the corresponding gene, and 25% (=2 of 8) had loss-of-heterozygosity (LOH) in chromosomal regions of two AHM-related genes. CONCLUSIONS: This study proposes that the likely molecular mechanism for development of amelanogenesis in AHM is carriage of an albinism/hypopigmentation allele followed by LOH of the corresponding gene in the tumour.

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

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

The MC1R r allele does not increase melanoma risk in MITF E318K carriers.

BACKGROUND: Population-wide screening for melanoma is not cost-effective, but genetic characterization could facilitate risk stratification and targeted screening. Common Melanocortin-1 receptor (MC1R) red hair colour (RHC) variants and Microphthalmia-associated transcription factor (MITF) E318K separately confer moderate melanoma susceptibility, but their interactive effects are relatively unexplored. OBJECTIVES: To evaluate whether MC1R genotypes differentially affect melanoma risk in MITF E318K+ vs. E318K- individuals. MATERIALS AND METHODS: Melanoma status (affected or unaffected) and genotype data (MC1R and MITF E318K) were collated from research cohorts (five Australian and two European). In addition, RHC genotypes from E318K+ individuals with and without melanoma were extracted from databases (The Cancer Genome Atlas and Medical Genome Research Bank, respectively). χ2 and logistic regression were used to evaluate RHC allele and genotype frequencies within E318K+/- cohorts depending on melanoma status. Replication analysis was conducted on 200 000 general-population exomes (UK Biobank). RESULTS: The cohort comprised 1165 MITF E318K- and 322 E318K+ individuals. In E318K- cases MC1R R and r alleles increased melanoma risk relative to wild type (wt), P < 0.001 for both. Similarly, each MC1R RHC genotype (R/R, R/r, R/wt, r/r and r/wt) increased melanoma risk relative to wt/wt (P < 0.001 for all). In E318K+ cases, R alleles increased melanoma risk relative to the wt allele [odds ratio (OR) 2.04 (95% confidence interval 1.67-2.49); P = 0.01], while the r allele risk was comparable with the wt allele [OR 0.78 (0.54-1.14) vs. 1.00, respectively]. E318K+ cases with the r/r genotype had a lower but not significant melanoma risk relative to wt/wt [OR 0.52 (0.20-1.38)]. Within the E318K+ cohort, R genotypes (R/R, R/r and R/wt) conferred a significantly higher risk compared with non-R genotypes (r/r, r/wt and wt/wt) (P < 0.001). UK Biobank data supported our findings that r did not increase melanoma risk in E318K+ individuals. CONCLUSIONS: RHC alleles/genotypes modify melanoma risk differently in MITF E318K- and E318K+ individuals. Specifically, although all RHC alleles increase risk relative to wt in E318K- individuals, only MC1R R increases melanoma risk in E318K+ individuals. Importantly, in the E318K+ cohort the MC1R r allele risk is comparable with wt. These findings could inform counselling and management for MITF E318K+ individuals.

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

Multiplex melanoma families are enriched for polygenic risk.

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

The Heritability of Twinning in Seven Large Historic Pedigrees.

It is widely recognized that dizygotic twinning (DZT) runs in families, but estimates of heritability from twin and family data are remarkably scarce and vary considerably. Here, we traced seven large, sometimes historical, multigeneration pedigrees from West Africans, fin de siècle French Jews, Canadians (two pedigrees), and the French royal family, in which twin births were recorded. We estimated heritability of twinning (of all types) as zygosity information was not available, diluting the true DZT heritability by a third or so. The estimates in the range 8-20% are remarkably consistent across time (8-19 generations) and ethnicities and also consistent with twin and family estimates.

Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes.

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4 Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

'Essential Tremor' Phenotype in FMR1 Premutation/Gray Zone Sibling Series: Exploring Possible Genetic Modifiers.

Fragile X-associated tremor/ataxia syndrome (FXTAS) occurs in carriers of fragile X mental retardation 1 (FMR1) X-linked small CGG expansion (gray zone [GZ] and premutation [PM]) alleles, containing 41-200 repeats. Major features comprise kinetic tremor, gait ataxia, cognitive decline and cerebellar peduncular white matter lesions, but atypical/incomplete FXTAS may occur. We explored the possibility of polygenic effects modifying the FXTAS spectrum phenotypes. We used three motor scales and selected cognitive tests in a series of three males and three females from a single sibship carrying PM or GZ alleles (44 to 75 repeats). The molecular profiles from these siblings were determined by genomewide association study with single-nucleotide polymorphism (SNP) genotyping by Illumina Global Screening Array. Nonparametric linkage analysis was applied and Parkinson's disease (PD) polygenic risk scores (PRSs) were calculated for all the siblings, based on 107 known risk variants. All male and female siblings manifested similar kinetic tremor phenotypes. In contrast to FXTAS, they showed negligible gait ataxia, and few white matter lesions on MRI. Cognitive functioning was unaffected. Suggestive evidence of linkage to a broad region of the short arm of chromosome 10 was obtained, and median PD PRS for the sibship fell within the top 30% of a sample of over 500,000 UK and Australian controls. The genomewide study results are suggestive of modifying effects of genetic risk loci linked to PD, on the neurological phenotype of FMR1-CGG small expansion carriers, resulting in an oligosymptomatic kinetic tremor seen in FXTAS spectrum, but also consistent with essential tremor.

Noncoding Variations in the Gene Encoding Ceramide Synthase 6 are Associated with Type 2 Diabetes in a Large Indigenous Australian Pedigree.

Type 2 diabetes (T2D) is a chronic disease that disproportionately affects Indigenous Australians. We have previously reported the localization of a novel T2D locus by linkage analysis to chromosome 2q24 in a large admixed Indigenous Australian pedigree (Busfield et al. (2002). American Journal of Human Genetics, 70, 349-357). Here we describe fine mapping of this region in this pedigree, with the identification of SNPs showing strong association with T2D: rs3845724 (diabetes p = 7 × 10-4), rs4668106 (diabetes p = 9 × 10-4) and rs529002 (plasma glucose p = 3 × 10-4). These associations were successfully replicated in an independent collection of Indigenous Australian T2D cases and controls. These SNPs all lie within the gene encoding ceramide synthase 6 (CERS6) and thus may regulate ceramide synthesis.

Variation at DENND1B and Asthma on the Island of Tristan da Cunha.

A high prevalence of asthma has been documented among the inhabitants of Tristan da Cunha, an isolated island in the South Atlantic. The population derives from just 28 founders. We performed lung function testing, including methacholine inhalation challenge, allergen skin prick testing, and collected DNA from essentially all of the current island population (269 individuals), and genotyped a panel of 43 single-nucleotide polymorphisms (SNPs) reported as associated with asthma and atopy. We carried out a mixed-model association analysis using the known pedigree. There were 96 individuals diagnosed as asthmatic (36%), and heritability estimates were similar to those from nonisolated population samples (multifactorial threshold model, h2 = 48%). The first component from a genetic principal components analysis using the entire SNP panel was nonlinearly associated with asthma, with the maximum risk to those intermediate to reference (Human Genome Diversity Project) European and African samples means. The single most strongly associated SNP was rs2786098 (p = 5.5 × 10-5), known to regulate the gene DENND1B. This explained approximately one-third of the trait heritability, with an allelic odds ratio for the A allele of 2.6. Among A/A carriers, 10 out of 12 individuals were asthmatic. The rs2786098*A variant was initially reported to decrease the risk of childhood (atopic) asthma in European but slightly increase the risk in African-descended populations, and does significantly alter Th2 cell function. Despite an absence of overall association with this variant in recent asthma genome wide association studies meta-analyses, an effect may exist on the particular genetic background of the Tristan da Cunha population.

The Relationship Between Personality and Somatic and Psychological Distress: A Comparison of Chinese and Australian Adolescents.

The extent to which correlations between personality domains and physical and psychological health generalize cross-culturally is unclear. We compared the strength of associations between the personality domains and somatic and psychological distress in Chinese (N = 2069) and a genetically informative sample of Australian (N = 2936) adolescents. We also examined the genetic and environmental etiology between personality, somatic and psychological distress in an Australian sample of 390 monozygotic twins and 698 dizygotic twins. In both populations, personality was assessed using the Junior Eysenck Personality Questionnaire. Somatic and psychological distress was assessed using the Somatic and Psychological Health Report. We found significant cultural differences in the relationship between adolescents' personality traits and somatic and psychological distress. Extraversion was positively associated with somatic distress in the Chinese but not in Australian adolescents. In the Australian twins, genetic covariation between neuroticism and somatic and psychological distress was stronger compared to the genetic associations between either psychoticism or extraversion with psychological distress.

A novel recurrent mutation in MITF predisposes to familial and sporadic melanoma.

So far, two genes associated with familial melanoma have been identified, accounting for a minority of genetic risk in families. Mutations in CDKN2A account for approximately 40% of familial cases, and predisposing mutations in CDK4 have been reported in a very small number of melanoma kindreds. Here we report the whole-genome sequencing of probands from several melanoma families, which we performed in order to identify other genes associated with familial melanoma. We identify one individual carrying a novel germline variant (coding DNA sequence c.G1075A; protein sequence p.E318K; rs149617956) in the melanoma-lineage-specific oncogene microphthalmia-associated transcription factor (MITF). Although the variant co-segregated with melanoma in some but not all cases in the family, linkage analysis of 31 families subsequently identified to carry the variant generated a log of odds (lod) score of 2.7 under a dominant model, indicating E318K as a possible intermediate risk variant. Consistent with this, the E318K variant was significantly associated with melanoma in a large Australian case-control sample. Likewise, it was similarly associated in an independent case-control sample from the United Kingdom. In the Australian sample, the variant allele was significantly over-represented in cases with a family history of melanoma, multiple primary melanomas, or both. The variant allele was also associated with increased naevus count and non-blue eye colour. Functional analysis of E318K showed that MITF encoded by the variant allele had impaired sumoylation and differentially regulated several MITF targets. These data indicate that MITF is a melanoma-predisposition gene and highlight the utility of whole-genome sequencing to identify novel rare variants associated with disease susceptibility.

Familial aggregation of albuminuria and arterial hypertension in an Aboriginal Australian community and the contribution of variants in ACE and TP53.

BACKGROUND: Aboriginal Australians are at high risk of cardiovascular, metabolic and renal diseases, resulting in a marked reduction in life expectancy when compared to the rest of the Australian population. This is partly due to recognized environmental and lifestyle risk factors, but a contribution of genetic susceptibility is also likely. METHODS: Using results from a comprehensive survey of one community (N = 1350 examined individuals), we have tested for familial aggregation of plasma glucose, arterial blood pressure, albuminuria (measured as urinary albumin to creatinine ratio, UACR) and estimated glomerular filtration rate (eGFR), and quantified the contribution of variation at four candidate genes (ACE; TP53; ENOS3; MTHFR). RESULTS: In the subsample of 357 individuals with complete genotype and phenotype data we showed that both UACR (h2 = 64%) and blood pressure (sBP h2 = 29%, dBP, h2 = 11%) were significantly heritable. The ACE insertion-deletion (P = 0.0009) and TP53 codon72 polymorphisms (P = 0.003) together contributed approximately 15% of the total heritability of UACR, with an effect of ACE genotype on BP also clearly evident. CONCLUSIONS: While the effects of the ACE insertion-deletion on risk of renal disease (especially in the setting of diabetes) are well recognized, this is only the second study to implicate p53 genotype as a risk factor for albuminuria - the other being an earlier study we performed in a different Aboriginal community (McDonald et al., J Am Soc Nephrol 13: 677-83, 2002). We conclude that there are significant genetic contributions to the high prevalence of chronic diseases observed in this population.

Genetics of skin color variation in Europeans: genome-wide association studies with functional follow-up.

In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations.

Genome-wide association analysis identifies 11 risk variants associated with the asthma with hay fever phenotype.

BACKGROUND: To date, no genome-wide association study (GWAS) has considered the combined phenotype of asthma with hay fever. Previous analyses of family data from the Tasmanian Longitudinal Health Study provide evidence that this phenotype has a stronger genetic cause than asthma without hay fever. OBJECTIVE: We sought to perform a GWAS of asthma with hay fever to identify variants associated with having both diseases. METHODS: We performed a meta-analysis of GWASs comparing persons with both physician-diagnosed asthma and hay fever (n = 6,685) with persons with neither disease (n = 14,091). RESULTS: At genome-wide significance, we identified 11 independent variants associated with the risk of having asthma with hay fever, including 2 associations reaching this level of significance with allergic disease for the first time: ZBTB10 (rs7009110; odds ratio [OR], 1.14; P = 4 × 10(-9)) and CLEC16A (rs62026376; OR, 1.17; P = 1 × 10(-8)). The rs62026376:C allele associated with increased asthma with hay fever risk has been found to be associated also with decreased expression of the nearby DEXI gene in monocytes. The 11 variants were associated with the risk of asthma and hay fever separately, but the estimated associations with the individual phenotypes were weaker than with the combined asthma with hay fever phenotype. A variant near LRRC32 was a stronger risk factor for hay fever than for asthma, whereas the reverse was observed for variants in/near GSDMA and TSLP. Single nucleotide polymorphisms with suggestive evidence for association with asthma with hay fever risk included rs41295115 near IL2RA (OR, 1.28; P = 5 × 10(-7)) and rs76043829 in TNS1 (OR, 1.23; P = 2 × 10(-6)). CONCLUSION: By focusing on the combined phenotype of asthma with hay fever, variants associated with the risk of allergic disease can be identified with greater efficiency.

Identification of a melanoma susceptibility locus and somatic mutation in TET2.

Although genetic studies have reported a number of loci associated with melanoma risk, the complex genetic architecture of the disease is not yet fully understood. We sought to identify common genetic variants associated with melanoma risk in a genome-wide association study (GWAS) of 2298 cases and 6654 controls. Thirteen of 15 known loci were replicated with nominal significance. A total of 69 single-nucleotide polymorphisms (SNPs) were selected for in silico replication in two independent melanoma GWAS datasets (a total of 5149 cases and 12 795 controls). Seven novel loci were nominally significantly associated with melanoma risk. These seven SNPs were further genotyped in 234 melanoma cases and 238 controls. The SNP rs4698934 was nominally significantly associated with melanoma risk. The combined odds ratio per T allele = 1.18; 95% confidence interval (1.10-1.25); combined P = 7.70 × 10(-) (7). This SNP is located in the intron of the TET2 gene on chromosome 4q24. In addition, a novel somatic mutation of TET2 was identified by next-generation sequencing in 1 of 22 sporadic melanoma cases. TET2 encodes a member of TET family enzymes that oxidizes 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). It is a putative epigenetic biomarker of melanoma as we previously reported, with observation of reduced TET2 transcriptional expression. This study is the first to implicate TET2 genetic variation and mutation in melanoma.

GWAS findings for human iris patterns: associations with variants in genes that influence normal neuronal pattern development.

Human iris patterns are highly variable. The origins of this variation are of interest in the study of iris-related eye diseases and forensics, as well as from an embryological developmental perspective, with regard to their possible relationship to fundamental processes of neurodevelopment. We have performed genome-wide association scans on four iris characteristics (crypt frequency, furrow contractions, presence of peripupillary pigmented ring, and number of nevi) in three Australian samples of European descent. Both the discovery (n = 2121) and replication (n = 499 and 73) samples showed evidence for association between (1) crypt frequency and variants in the axonal guidance gene SEMA3A (p = 6.6 × 10(-11)), (2) furrow contractions and variants within the cytoskeleton gene TRAF3IP1 (p = 2.3 × 10(-12)), and (3) the pigmented ring and variants in the well-known pigmentation gene SLC24A4 (p = 7.6 × 10(-21)). These replicated findings individually accounted for around 1.5%-3% of the variance for these iris characteristics. Because both SEMA3A and TRAFIP1 are implicated in pathways that control neurogenesis, neural migration, and synaptogenesis, we also examined the evidence of enhancement among such genes, finding enrichment for crypts and furrows. These findings suggest that genes involved in normal neuronal pattern development may also influence tissue structures in the human iris.

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.