Search

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

The original version of this Article contained errors in the spelling of the authors Fan Liu and M. Arfan Ikram, which were incorrectly given as Fan Lui and Arfan M. Ikram. In addition, the original version of this Article also contained errors in the author affiliations which are detailed in the associated Publisher Correction.

The total number of acquired melanocytic nevi on the skin is strongly correlated with melanoma risk. Here we report a meta-analysis of 11 nevus GWAS from Australia, Netherlands, UK, and USA comprising 52,506 individuals. We confirm known loci including MTAP, PLA2G6, and IRF4, and detect novel SNPs in KITLG and a region of 9q32. In a bivariate analysis combining the nevus results with a recent melanoma GWAS meta-analysis (12,874 cases, 23,203 controls), SNPs near GPRC5A, CYP1B1, PPARGC1B, HDAC4, FAM208B, DOCK8, and SYNE2 reached global significance, and other loci, including MIR146A and OBFC1, reached a suggestive level. Overall, we conclude that most nevus genes affect melanoma risk (KITLG an exception), while many melanoma risk loci do not alter nevus count. For example, variants in TERC and OBFC1 affect both traits, but other telomere length maintenance genes seem to affect melanoma risk only. Our findings implicate multiple pathways in nevogenesis.

Hair color is one of the most recognizable visual traits in European populations and is under strong genetic control. Here we report the results of a genome-wide association study meta-analysis of almost 300,000 participants of European descent. We identified 123 autosomal and one X-chromosome loci significantly associated with hair color; all but 13 are novel. Collectively, single-nucleotide polymorphisms associated with hair color within these loci explain 34.6% of red hair, 24.8% of blond hair, and 26.1% of black hair heritability in the study populations. These results confirm the polygenic nature of complex phenotypes and improve our understanding of melanin pigment metabolism in humans.

Asthma, hay fever (or allergic rhinitis) and eczema (or atopic dermatitis) often coexist in the same individuals, partly because of a shared genetic origin. To identify shared risk variants, we performed a genome-wide association study (GWAS; n = 360,838) of a broad allergic disease phenotype that considers the presence of any one of these three diseases. We identified 136 independent risk variants (P < 3 × 10-8), including 73 not previously reported, which implicate 132 nearby genes in allergic disease pathophysiology. Disease-specific effects were detected for only six variants, confirming that most represent shared risk factors. Tissue-specific heritability and biological process enrichment analyses suggest that shared risk variants influence lymphocyte-mediated immunity. Six target genes provide an opportunity for drug repositioning, while for 36 genes CpG methylation was found to influence transcription independently of genetic effects. Asthma, hay fever and eczema partly coexist because they share many genetic risk variants that dysregulate the expression of immune-related genes.

© 2016 The Author(s). 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.

Eczema often precedes the development of asthma in a disease course called the a € atopic marcha €. To unravel the genes underlying this characteristic pattern of allergic disease, we conduct a multi-stage genome-wide association study on infantile eczema followed by childhood asthma in 12 populations including 2,428 cases and 17,034 controls. Here we report two novel loci specific for the combined eczema plus asthma phenotype, which are associated with allergic disease for the first time; rs9357733 located in EFHC1 on chromosome 6p12.3 (OR 1.27; P=2.1 × 10 a'8) and rs993226 between TMTC2 and SLC6A15 on chromosome 12q21.3 (OR 1.58; P=5.3 × 10 a'9). Additional susceptibility loci identified at genome-wide significance are FLG (1q21.3), IL4/KIF3A (5q31.1), AP5B1/OVOL1 (11q13.1), C11orf30/LRRC32 (11q13.5) and IKZF3 (17q21). We show that predominantly eczema loci increase the risk for the atopic march. Our findings suggest that eczema may play an important role in the development of asthma after eczema.

Thirteen common susceptibility loci have been reproducibly associated with cutaneous malignant melanoma (CMM). We report the results of an international 2-stage meta-analysis of CMM genome-wide association studies (GWAS). This meta-analysis combines 11 GWAS (5 previously unpublished) and a further three stage 2 data sets, totaling 15,990 CMM cases and 26,409 controls. Five loci not previously associated with CMM risk reached genome-wide significance (P < 5 × 10(-8)), as did 2 previously reported but unreplicated loci and all 13 established loci. Newly associated SNPs fall within putative melanocyte regulatory elements, and bioinformatic and expression quantitative trait locus (eQTL) data highlight candidate genes in the associated regions, including one involved in telomere biology.

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.

The main genetic determinant of soluble interleukin 6 receptor (sIL-6R) levels is the missense variant rs2228145 that maps to the cleavage site of IL-6R. For each Ala allele, sIL-6R serum levels increase by ≈ 20 ng ml(-1) and asthma risk by 1.09-fold. However, this variant does not explain the total heritability for sIL-6R levels. Additional independent variants in IL6R may therefore contribute to variation in sIL-6R levels and influence asthma risk. We imputed 471 variants in IL6R and tested these for association with sIL-6R serum levels in 360 individuals. An intronic variant (rs12083537) was associated with sIL-6R levels independently of rs4129267 (P=0.0005), a proxy single-nucleotide polymorphism for rs2228145. A significant and consistent association for rs12083537 was observed in a replication panel of 354 individuals (P=0.033). Each rs12083537:A allele increased sIL-6R serum levels by 2.4 ng ml(-1). Analysis of mRNA levels in two cohorts did not identify significant associations between rs12083537 and IL6R transcription levels. On the other hand, results from 16,705 asthmatics and 30,809 controls showed that the rs12083537:A allele increased asthma risk by 1.04-fold (P=0.0419). Genetic risk scores based on IL6R regulatory variants may prove useful in explaining variation in clinical response to tocilizumab, an anti-IL-6R monoclonal antibody.

RATIONALE: Asthma has substantial morbidity and mortality and a strong genetic component, but identification of genetic risk factors is limited by availability of suitable studies. OBJECTIVES: To test if population-based cohorts with self-reported physician-diagnosed asthma and genome-wide association (GWA) data could be used to validate known associations with asthma and identify novel associations. METHODS: The APCAT (Analysis in Population-based Cohorts of Asthma Traits) consortium consists of 1,716 individuals with asthma and 16,888 healthy controls from six European-descent population-based cohorts. We examined associations in APCAT of thirteen variants previously reported as genome-wide significant (P<5 x 10(-8)) and three variants reported as suggestive (P<5× 10(-7)). We also searched for novel associations in APCAT (Stage 1) and followed-up the most promising variants in 4,035 asthmatics and 11,251 healthy controls (Stage 2). Finally, we conducted the first genome-wide screen for interactions with smoking or hay fever. MAIN RESULTS: We observed association in the same direction for all thirteen previously reported variants and nominally replicated ten of them. One variant that was previously suggestive, rs11071559 in RORA, now reaches genome-wide significance when combined with our data (P = 2.4 × 10(-9)). We also identified two genome-wide significant associations: rs13408661 near IL1RL1/IL18R1 (P(Stage1+Stage2) = 1.1x10(-9)), which is correlated with a variant recently shown to be associated with asthma (rs3771180), and rs9268516 in the HLA region (P(Stage1+Stage2) = 1.1x10(-8)), which appears to be independent of previously reported associations in this locus. Finally, we found no strong evidence for gene-environment interactions with smoking or hay fever status. CONCLUSIONS: Population-based cohorts with simple asthma phenotypes represent a valuable and largely untapped resource for genetic studies of asthma.

We performed a genome-wide association study of melanoma in a discovery cohort of 2,168 Australian individuals with melanoma and 4,387 control individuals. In this discovery phase, we confirm several previously characterized melanoma-associated loci at MC1R, ASIP and MTAP-CDKN2A. We selected variants at nine loci for replication in three independent case-control studies (Europe: 2,804 subjects with melanoma, 7,618 control subjects; United States 1: 1,804 subjects with melanoma, 1,026 control subjects; United States 2: 585 subjects with melanoma, 6,500 control subjects). The combined meta-analysis of all case-control studies identified a new susceptibility locus at 1q21.3 (rs7412746, P = 9.0 × 10(-11), OR in combined replication cohorts of 0.89 (95% CI 0.85-0.95)). We also show evidence suggesting that melanoma associates with 1q42.12 (rs3219090, P = 9.3 × 10(-8)). The associated variants at the 1q21.3 locus span a region with ten genes, and plausible candidate genes for melanoma susceptibility include ARNT and SETDB1. Variants at the 1q21.3 locus do not seem to be associated with human pigmentation or measures of nevus density.

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.

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

We hypothesized that variants in genes expressed as a consequence of interactions between ovarian cancer cells and the host micro-environment could contribute to cancer susceptibility. We therefore used a two-stage approach to evaluate common single nucleotide polymorphisms (SNPs) in 173 genes involved in stromal epithelial interactions in the Ovarian Cancer Association Consortium (OCAC). In the discovery stage, cases with epithelial ovarian cancer (n=675) and controls (n=1,162) were genotyped at 1,536 SNPs using an Illumina GoldenGate assay. Based on Positive Predictive Value estimates, three SNPs-PODXL rs1013368, ITGA6 rs13027811, and MMP3 rs522616-were selected for replication using TaqMan genotyping in up to 3,059 serous invasive cases and 8,905 controls from 16 OCAC case-control studies. An additional 18 SNPs with Pper-allele<0.05 in the discovery stage were selected for replication in a subset of five OCAC studies (n=1,233 serous invasive cases; n=3,364 controls). The discovery stage associations in PODXL, ITGA6, and MMP3 were attenuated in the larger replication set (adj. Pper-allele>or=0.5). However genotypes at TERT rs7726159 were associated with ovarian cancer risk in the smaller, five-study replication study (Pper-allele=0.03). Combined analysis of the discovery and replication sets for this TERT SNP showed an increased risk of serous ovarian cancer among non-Hispanic whites [adj. ORper-allele 1.14 (1.04-1.24) p=0.003]. Our study adds to the growing evidence that, like the 8q24 locus, the telomerase reverse transcriptase locus at 5p15.33, is a general cancer susceptibility locus.

Previous studies have successfully identified genetic variants in several genes associated with human iris (eye) color; however, they all used simplified categorical trait information. Here, we quantified continuous eye color variation into hue and saturation values using high-resolution digital full-eye photographs and conducted a genome-wide association study on 5,951 Dutch Europeans from the Rotterdam Study. Three new regions, 1q42.3, 17q25.3, and 21q22.13, were highlighted meeting the criterion for genome-wide statistically significant association. The latter two loci were replicated in 2,261 individuals from the UK and in 1,282 from Australia. The LYST gene at 1q42.3 and the DSCR9 gene at 21q22.13 serve as promising functional candidates. A model for predicting quantitative eye colors explained over 50% of trait variance in the Rotterdam Study. Over all our data exemplify that fine phenotyping is a useful strategy for finding genes involved in human complex traits.

Background: Linkage studies have implicated the 2q33, 9p21, 11q13 and 20q13 regions in the regulation of allergic disease. The aim of this study was to test genetic variants in candidate genes from these regions for association with specific asthma traits. Methods: Ninety-five single nucleotide polymorphisms (SNP) located in eight genes (CD28, CTLA4, ICOS, ADAM23, ADAMTSL1, MS4A2, CDH26 and HRH3) were genotyped in >5000 individuals from Australian (n = 1162), Dutch (n = 99) and Danish (n = 303) families. Traits tested included doctor-diagnosed asthma, atopy, airway obstruction, total serum immunoglobulin (Ig) E levels and eosinophilia. Association was tested using both multivariate and univariate methods, with gene-wide thresholds for significance determined through simulation. Gene-by-gene and gene-by-environment analyses were also performed. Results: There was no overall evidence for association with seven of the eight genes tested when considering all genetic variation assayed in each gene. The exception was MS4A2 on chromosome 11q13, which showed weak evidence for association with IgE (gene-wide P < 0.05, rs502581). There were no significant gene-by-gene or gene-by-environment interaction effects after accounting for the number of tests performed. Conclusions: The individual variants genotyped in the 2q33, 9p21 and 20q13 regions do not explain a large fraction of the variation in the quantitative traits tested or have a major impact on asthma or atopy risk. Our results are consistent with a weak effect of MS4A2 polymorphisms on the variation of total IgE levels. © 2009 John Wiley & Sons A/S.

BACKGROUND: Spontaneous dizygotic (DZ) twinning in humans is under genetic control. In sheep, heterozygous loss of function mutations in bone morphogenetic protein 15 (BMP15) increase ovulation and hence twinning rates. METHODS: To investigate the role of BMP15 in human twinning, we typed 14 common variants, 4 rare novel variants initially detected by sequencing 279 mothers of DZ twins (MODZT) and 17 variants previously associated with premature ovarian failure (POF) in 933 DZ twinning families. We also typed five additional POF associated GDF9 variants. RESULTS: There was some evidence for association between DZ twinning and a common intronic BMP15 variant (rs3897937), but this was not significant after correction for multiple testing. Three of the four novel variants (p.Pro174Ser, p.Ala311Thr and p.Arg392Thr) occurred in 1-5 MODZT but were not detected in 1512 controls. We also detected three POF associated mutations in both BMP15 and GDF9 at low frequencies in MODZT and controls. CONCLUSIONS: We conclude that neither rare nor common BMP15 variants play a significant role in the variation in human DZ twinning. © The Author 2008. Published by Oxford University Press on behalf of the European Society of Human Reproduction and Embryology. All rights reserved.

Context: Genes from the ovarian bone morphogenetic signaling pathway (GDF9 and BMP15) are critical for normal human fertility. We previously identified a deletion mutation in GDF9 in sisters with spontaneous dizygotic (DZ) twins, but the prevalence of rare GDF9 variants in twinning families is unknown. Objective: The objective was to evaluate the frequency of rare variants in GDF9 in families with a history of DZ twinning. Design and Subjects: We recruited 3450 individuals from 915 DZ twinning families (1693 mothers of twins) and 1512 controls of Caucasian origin. One mother of DZ twins was selected from 279 of the 915 families, and a DNA sample was screened for rare variants in GDF9 using denaturant HPLC. Variants were confirmed by DNA sequencing and genotyped in the entire sample by matrix-assisted laser desorption ionization time of flight (MALDI-TOF) mass spectrometry. Results: We found two novel insertion/deletions (c.392-393insT, c.1268-1269delAA) and four missense alterations in the GDF9 sequence in mothers of twins. Two of the missense variants (c.307C>T, p.Pro103Ser and c.362C>T, p.Thr121Leu) were located in the proregion of GDF9 and two (c.1121C>T, p.Pro374Leu and c.1360C>T, p.Arg454Cys) in the mature protein region. For each variant, the frequencies were higher in cases compared with controls. The proportion of mothers of DZ twins carrying any variant (4.12%) was significantly higher (P < 0.0001) than the proportion of carriers in controls (2.29%). Conclusion: We describe new variants in the GDF9 gene that are significantly more common in mothers of DZ twins than controls, suggesting that rare GDF9 variants contribute to the likelihood of DZ twinning. Copyright © 2006 by The Endocrine Society.

BACKGROUND: The liability to asthma is influenced both by genetic and environmental factors. The objective of this study was to identify risk factors for asthma in young adult twin pairs during an 8-year period. METHODS: From the birth cohorts 1953-1982 of the Danish Twin Registry, 6,090 twin pairs who were initially unaffected with respect to asthma at a nationwide questionnaire-based study in 1994 participated in a similar follow-up study in 2002. Subjects were regarded incident asthma cases when responding affirmatively to the question 'Do you have, or have you ever had asthma'? in 2002. Pairs in which only one twin developed asthma -- discordant pairs -- were identified and conditional logistic regression was applied to detect effects of risk factors. RESULTS: A total of 126 monozygotic (MZ) and 273 dizygotic (DZ) discordant twin pairs were identified. In MZ twins hay fever (OR = 3.16, 95% CI: 1.29-7.73, P = 0.007) and exercise (OR for inactivity = 0.35, 95% CI: 0.13-0.91, P = 0.023) were significantly associated with asthma, whereas in DZ twins, hay fever (OR = 2.44, 95% CI: 1.44-4.13, P = 0.001), eczema (OR = 1.96, 95% CI: 1.02-3.78, P = 0.040), female sex (OR between males and females = 0.54, 95% CI: 0.36-0.80, P = 0.002), and increasing levels of body mass index (BMI; OR per unit = 1.11, 95% CI: 1.02-1.20, P = 0.009) were significant predictors of asthma. CONCLUSIONS: Hay fever, eczema, female sex, exercise and increasing levels of BMI were risk factors for asthma in young adults. The different risk profile observed in MZ twins compared with DZ twins may reflect an underlying genetic vulnerability shared between those risk factors and asthma.

After ingestion of a standardized dose of ethanol, alcohol concentrations were assessed, over 3.5 hours from blood (six readings) and breath (10 readings) in a sample of 412 MZ and DZ twins who took part in an Alcohol Challenge Twin Study (ACTS). Nearly all participants were subsequently genotyped on two polymorphic SNPs in the ADH1B and ADH1C loci known to affect in vitro ADH activity. In the DZ pairs, 14 microsatellite markers covering a 20.5 cM region on chromosome 4 that includes the ADH gene family were assessed, Variation in the timed series of autocorrelated blood and breath alcohol readings was studied using a bivariate simplex design. The contribution of a quantitative trait locus (QTL) or QTL's linked to the ADH region was estimated via a mixture of likelihoods weighted by identity-by-descent probabilities. The effects of allelic substitution at the ADH1B and ADH1C loci were estimated in the means part of the model simultaneously with the effects sex and age. There was a major contribution to variance in alcohol metabolism due to a QTL which accounted for about 64% of the additive genetic covariation common to both blood and breath alcohol readings at the first time point. No effects of the ADH1B*47His or ADH1C*349Ile alleles on in vivo metabolism were observed, although these have been shown to have major effects in vitro. This implies that there is a major determinant of variation for in vivo alcohol metabolism in the ADH region that is not accounted for by these polymorphisms. Earlier analyses of these data suggested that alcohol metabolism is related to drinking behavior and imply that this QTL may be protective against alcohol dependence. © 2005 Springer Science+Business Media, Inc.

INTRODUCTION: The androgen receptor (AR) gene exon 1 CAG repeat polymorphism encodes a string of 9-32 glutamines. Women with germline BRCA1 mutations who carry at least one AR allele with 28 or more repeats have been reported to have an earlier age at onset of breast cancer. METHODS: A total of 604 living female Australian and British BRCA1 and/or BRCA2 mutation carriers from 376 families were genotyped for the AR CAG repeat polymorphism. The association between AR genotype and disease risk was assessed using Cox regression. AR genotype was analyzed as a dichotomous covariate using cut-points previously reported to be associated with increased risk among BRCA1 mutation carriers, and as a continuous variable considering smaller allele, larger allele and average allele size. RESULTS: There was no evidence that the AR CAG repeat polymorphism modified disease risk in the 376 BRCA1 or 219 BRCA2 mutation carriers screened successfully. The rate ratio associated with possession of at least one allele with 28 or more CAG repeats was 0.74 (95% confidence interval 0.42-1.29; P = 0.3) for BRCA1 carriers, and 1.12 (95% confidence interval 0.55-2.25; P = 0.8) for BRCA2 carriers. CONCLUSION: The AR exon 1 CAG repeat polymorphism does not appear to have an effect on breast cancer risk in BRCA1 or BRCA2 mutation carriers.

A loss of function mutation in growth differentiation factor 9 (GDF9) in sheep causes increased ovulation rate and infertility in a dosage-sensitive manner. Spontaneous dizygotic (DZ) twinning in the human is under genetic control and women with a history of DZ twinning have an increased incidence of multiple follicle growth and multiple ovulation. We sequenced the GDF9 coding region in DNA samples from 20 women with DZ twins and identified a four-base pair deletion in GDF9 in two sisters with twins from one family. We screened a further 429 families and did not find the loss of function mutation in any other families. We genotyped eight single nucleotide polymorphisms across the GDF9 locus in 379 families with two sisters who have both given birth to spontaneous DZ twins (1527 individuals) and 226 triad families with mothers of twins and their parents (723 individuals). Using case control analysis and the transmission disequilibrium test we found no evidence for association between common variants in GDF9 and twinning in the families. We conclude that rare mutations in GDF9 may influence twinning, but twinning frequency is not associated with common variation in GDF9.

Background: Folate metabolism is critical to embryonic development, influencing neural tube defects (NTD) and recurrent early pregnancy loss. Polymorphisms in 5,10-methylenetetrahydrofolate reductase (MTHFR) have been associated with dizygotic (DZ) twinning through pregnancy loss. Methods: The C677T and A1298C polymorphisms in MTHFR were genotyped in 258 Australasian families (1016 individuals) and 118 Dutch families (462 individuals) of mothers of DZ twins and a population sample of 462 adolescent twin families (1861 individuals). Haplotypes were constructed from the alleles, and transmission of the MTHFR haplotypes to mothers of DZ twins and from parents to twins in the adolescent twin families analysed. Results: The C677T and A1298C were common in all three populations (frequencies > 0.29). There was strong linkage disequilibrium (D′ = 1) between the variants, showing that specific combinations of alleles (haplotypes) were transmitted together. Three haplotypes accounted for nearly all the variation. There was no evidence of any association between MTHFR genotype and twinning in mothers of twins, or of the loss of specific MTHFR genotypes during twin pregnancies. Conclusions: It is concluded that variation in twinning frequency is not associated with MTHFR genotype.

The tendency to dizygotic (DZ) twinning is inherited in both humans and sheep, and a fecundity gene in sheep (FecB) maps to sheep chromosome 6, syntenic with human 4q21-25. Our aim was to see whether a gene predisposing to human DZ twinning mapped to this region. DNA was collected from 169 pairs and 17 sets of 3 sisters (trios) from Australia and New Zealand who had each had spontaneous DZ twins, mostly before the age of 35, and from a replication sample of 111 families (92 affected sister pairs) from The Netherlands. Exclusion mapping was carried out after typing 26 markers on chromosome 4, of which 8 spanned the region likely to contain the human homologue of the sheep FecB gene. We used nonparametric affected sib pair methods for linkage analysis [ASPEX 2.2, Hinds and Risch, 1999]. Complete exclusion of linkage (lod < -2) of a gene conferring a relative risk for sibs as low as 1.5 (λ