Megan J. Campbell

Common sequence variants on 20q11.22 confer melanoma susceptibility

We conducted a genome-wide association pooling study for cutaneous melanoma and performed validation in samples totaling 2,019 cases and 2,105 controls. Using pooling, we identified a new melanoma risk locus on chromosome 20 (rs910873 and rs1885120), with replication in two further samples (combined P < 1 × 10−15). The per allele odds ratio was 1.75 (1.53, 2.01), with evidence for stronger association in early-onset cases.

Common Variants in the Trichohyalin Gene Are Associated with Straight Hair in Europeans

Hair morphology is highly differentiated between populations and among people of European ancestry. Whereas hair morphology in East Asian populations has been studied extensively, relatively little is known about the genetics of this trait in Europeans. We performed a genome-wide association scan for hair morphology (straight, wavy, curly) in three Australian samples of European descent. All three samples showed evidence of association implicating the Trichohyalin gene (TCHH), which is expressed in the developing inner root sheath of the hair follicle, and explaining approximately 6% of variance (p=1.5x10(-31)). These variants are at their highest frequency in Northern Europeans, paralleling the distribution of the straight-hair EDAR variant in Asian populations.

Common variants in TMPRSS6 are associated with iron status and erythrocyte volume

We report a genome-wide association study to iron status. We identify an association of SNPs in TPMRSS6 to serum iron (rs855791, combined P = 1.5 × 10−20), transferrin saturation (combined P = 2.2 × 10−23) and erythrocyte mean cell volume (MCV, combined P = 1.1 × 10−10). We also find suggestive evidence of association with blood hemoglobin levels (combined P = 5.3 × 10−7). These findings demonstrate the involvement of TMPRSS6 in control of iron homeostasis and in normal erythropoiesis.

Increased DNA Methylation at the AXIN1 Gene in a Monozygotic Twin from a Pair Discordant for a Caudal Duplication Anomaly

The AXIN1 gene has been implicated in caudal duplication anomalies. Its coding region was sequenced in both members of a monozygotic (MZ) twin pair discordant for a caudal duplication anomaly, but no mutation was found. Using bisulfite sequencing, we examined methylation at the promoter region of the AXIN1 gene in these twins and in twin and age-matched singleton controls. Methylation of the promoter region in peripheral blood mononucleated cells was variable among individuals, including MZ pairs. In the MZ pair discordant for the caudal duplication, this region of the affected twin was significantly more methylated than that of the unaffected twin (P < .0001), which was significantly more methylated than those of the controls (P = .02). We have confirmed that this CpG island does function as a promoter in vitro and that its activity is inversely proportional to the extent of methylation. This finding raises the possibility that hypermethylation of the AXIN1 promoter, by mechanisms as yet undetermined, is associated with the malformation. This case may be paradigmatic for some cases of MZ discordance.

The perception of quinine taste intensity is associated with common genetic variants in a bitter receptor cluster on chromosome 12

The perceived taste intensities of quinine HCl, caffeine, sucrose octaacetate (SOA) and propylthiouracil (PROP) solutions were examined in 1457 twins and their siblings. Previous heritability modeling of these bitter stimuli indicated a common genetic factor for quinine, caffeine and SOA (22–28%), as well as separate specific genetic factors for PROP (72%) and quinine (15%). To identify the genes involved, we performed a genome-wide association study with the same sample as the modeling analysis, genotyped for approximately 610 000 single-nucleotide polymorphisms (SNPs). For caffeine and SOA, no SNP association reached a genome-wide statistical criterion. For PROP, the peak association was within TAS2R38 (rs713598, A49P, P = 1.6 × 10−104), which accounted for 45.9% of the trait variance. For quinine, the peak association was centered in a region that contains bitter receptor as well as salivary protein genes and explained 5.8% of the trait variance (TAS2R19, rs10772420, R299C, P = 1.8 × 10−15). We confirmed this association in a replication sample of twins of similar ancestry (P = 0.00001). The specific genetic factor for the perceived intensity of PROP was identified as the gene previously implicated in this trait (TAS2R38). For quinine, one or more bitter receptor or salivary proline-rich protein genes on chromosome 12 have alleles which affect its perception but tight linkage among very similar genes precludes the identification of a single causal genetic variant.

A Variant in LIN28B Is Associated with 2D:4D Finger-Length Ratio, a Putative Retrospective Biomarker of Prenatal Testosterone Exposure

The ratio of the lengths of an individuals second to fourth digit (2D:4D) is commonly used as a noninvasive retrospective biomarker for prenatal androgen exposure. In order to identify the genetic determinants of 2D:4D, we applied a genome-wide association approach to 1507 11-year-old children from the Avon Longitudinal Study of Parents and Children (ALSPAC) in whom 2D:4D ratio had been measured, as well as a sample of 1382 12- to 16-year-olds from the Brisbane Adolescent Twin Study. A meta-analysis of the two scans identified a single variant in the LIN28B gene that was strongly associated with 2D:4D (rs314277: p = 4.1 x 10(-8)) and was subsequently independently replicated in an additional 3659 children from the ALSPAC cohort (p = 1.53 x 10(-6)). The minor allele of the rs314277 variant has previously been linked to increased height and delayed age at menarche, but in our study it was associated with increased 2D:4D in the direction opposite to that of previous reports on the correlation between 2D:4D and age at menarche. Our findings call into question the validity of 2D:4D as a simplistic retrospective biomarker for prenatal testosterone exposure.

Impact of the Genome on the Epigenome Is Manifested in DNA Methylation Patterns of Imprinted Regions in Monozygotic and Dizygotic Twins

One of the best studied read-outs of epigenetic change is the differential expression of imprinted genes, controlled by differential methylation of imprinted control regions (ICRs). To address the impact of genotype on the epigenome, we performed a detailed study in 128 pairs of monozygotic (MZ) and 128 pairs of dizygotic (DZ) twins, interrogating the DNA methylation status of the ICRs of IGF2, H19, KCNQ1, GNAS and the non-imprinted gene RUNX1. While we found a similar overall pattern of methylation between MZ and DZ twins, we also observed a high degree of variability in individual CpG methylation levels, notably at the H19/IGF2 loci. A degree of methylation plasticity independent of the genome sequence was observed, with both local and regional CpG methylation changes, discordant between MZ and DZ individual pairs. However, concordant gains or losses of methylation, within individual twin pairs were more common in MZ than DZ twin pairs, indicating that de novo and/or maintenance methylation is influenced by the underlying DNA sequence. Specifically, for the first time we showed that the rs10732516 [A] polymorphism, located in a critical CTCF binding site in the H19 ICR locus, is strongly associated with increased hypermethylation of specific CpG sites in the maternal H19 allele. Together, our results highlight the impact of the genome on the epigenome and demonstrate that while DNA methylation states are tightly maintained between genetically identical and related individuals, there remains considerable epigenetic variation that may contribute to disease susceptibility.

Quantitative Trait Loci for CD4:CD8 Lymphocyte Ratio Are Associated with Risk of Type 1 Diabetes and HIV-1 Immune Control

Abnormal expansion or depletion of particular lymphocyte subsets is associated with clinical manifestations such as HIV progression to AIDS and autoimmune disease. We sought to identify genetic predictors of lymphocyte levels and reasoned that these may play a role in immune-related diseases. We tested 2.3 million variants for association with five lymphocyte subsets, measured in 2538 individuals from the general population, including CD4+ T cells, CD8+ T cells, CD56+ natural killer (NK) cells, and the derived measure CD4:CD8 ratio. We identified two regions of strong association. The first was located in the major histocompatibility complex (MHC), with multiple SNPs strongly associated with CD4:CD8 ratio (rs2524054, p = 2.1 x 10(-28)). The second region was centered within a cluster of genes from the Schlafen family and was associated with NK cell levels (rs1838149, p = 6.1 x 10(-14)). The MHC association with CD4:CD8 replicated convincingly (p = 1.4 x 10(-9)) in an independent panel of 988 individuals. Conditional analyses indicate that there are two major independent quantitative trait loci (QTL) in the MHC region that regulate CD4:CD8 ratio: one is located in the class I cluster and influences CD8 levels, whereas the second is located in the class II cluster and regulates CD4 levels. Jointly, both QTL explained 8% of the variance in CD4:CD8 ratio. The class I variants are also strongly associated with durable host control of HIV, and class II variants are associated with type-1 diabetes, suggesting that genetic variation at the MHC may predispose one to immune-related diseases partly through disregulation of T cell homeostasis.

Sequence variants in three loci influence monocyte counts and erythrocyte volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p=8.9x10(-14)). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p=1.2x10(-9), 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p=5.3x10(-9), 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values=0.001, 9.9x10(-5), and 7x10(-5), respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.

Estimation of the Rate of SNP Genotyping Errors From DNA Extracted From Different Tissues

High density single nucleotide polymorphism (SNP) genotyping panels provide an alternative to microsatellite markers for genome scans. However, genotype errors have a major impact on power to detect linkage or association and are difficult to detect for SNPs. We estimated error rates with the Affymetrix GeneChip SNP platform in samples from a family with a mixed set of monozygotic (MZ) and dizygotic (DZ) triplets using lymphocyte, buccal DNA and samples from whole genome amplification using the multiple displacement amplification (MDA) technique. The average call rate from 58,960 SNPs for five genomic samples was 99.48%. Comparison of results for the MZ twins showed only three discordant genotypes (concordance rate 99.995%). The mean concordance rate for comparisons of samples from lymphocyte and buccal DNA was 99.97%. Mendelian inconsistencies were identified in 46 SNPs with errors in one or more family members, a rate of 0.022%. Observed genotype concordance rates between parents, between parents and children, and among siblings were consistent with previously reported allele frequencies and Hardy-Weinberg equilibrium. Using the MDA technique, results for two samples had equivalent high accuracy to results with genomic samples. However, the SNP call rate for the remaining seven samples varied from 72.5% to 99.5%, with an average of 86.11%. Quality of the DNA sample following the MDA reaction appears to be the critical factor in SNP call rate for MDA samples. Our results demonstrate highly accurate and reproducible genotyping for the Affymetrix GeneChip Human Mapping Set in lymphocyte and buccal DNA samples.