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 individual's 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.

GWAS of DNA methylation variation within imprinting control regions suggests parent-of-origin association.

Imprinting control regions (ICRs) play a fundamental role in establishing and maintaining the non-random monoallelic expression of certain genes, via common regulatory elements such as non-coding RNAs and differentially methylated regions (DMRs) of DNA. We recently surveyed DNA methylation levels within four ICRs (H19-ICR, IGF2-DMR, KvDMR, and NESPAS-ICR) in whole-blood genomic DNA from 128 monozygotic (MZ) and 128 dizygotic (DZ) human twin pairs. Our analyses revealed high individual variation and intra-domain covariation in methylation levels across CpGs and emphasized the interaction between epigenetic variation and the underlying genetic sequence in a parent-of-origin fashion. Here, we extend our analysis to conduct two genome-wide screenings of single nucleotide polymorphisms (SNPs) underlying either intra-domain covariation or parent-of-origin-dependent association with methylation status at individual CpG sites located within ICRs. Although genome-wide significance was not surpassed due to sample size limitations, the most significantly associated SNPs found through multiple-trait genome-wide association (MQFAM) included the previously described rs10732516, which is located in the vicinity of the H19-ICR. Similarly, we identified an association between rs965808 and methylation status within the NESPAS-ICR. This SNP is positioned within an intronic region of the overlapping genes GNAS and GNAS-AS1, which are imprinted genes regulated by the NESPAS-ICR. Sixteen other SNPs located in regions apart from the analyzed regions displayed suggestive association with intra-domain methylation. Additionally, we identified 13 SNPs displaying parent-of-origin association with individual methylation sites through family-based association testing. In this exploratory study, we show the value and feasibility of using alternative GWAS approaches in the study of the interaction between epigenetic state and genetic sequence within imprinting regulatory domains. Despite the relatively small sample size, we identified a number of SNPs displaying suggestive association either in a domain-wide or in a parent-of-origin fashion. Nevertheless, these associations will require future experimental validation or replication in larger and independent samples.

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.

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.

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.

Two corpora lutea seen at 6-13 weeks' gestation infers dizygosity among spontaneous same-sexed dichorionic twins.

Current ultrasound techniques can accurately determine the chorionicity of twins, but not zygosity. We previously proposed that the zygosity of spontaneously conceived twins can be determined at early ultrasound, where 2 corpora lutea infers dizygosity, and 1 implies monozygosity. Here we did a case series, comparing zygosity predicted using this method with definitive DNA genotyping of twins after birth. We retrospectively identified 14 ultrasound reports of spontaneous twin pregnancies at 6(+0 days) to 13+6 weeks' gestation, where both ovaries were seen and the number of corpora lutea documented. We visited all twin pairs, obtained buccal smears, and determined zygosity by genotyping 9 independent microsatellite markers. All 8 cases where 2 corpora lutea were seen were dizygotic pregnancies. One further case where 3 corpora lutea were seen was also dizygotic. All 3 sets of monozygotic twins had 1 corpus luteum. There were 2 cases incorrectly assigned, where 1 corpus luteum was seen in dizygotic pregnancies. We conclude if 2 corpora lutea are seen at a first trimester ultrasound of spontaneously conceived dichorionic twins, they appear to be almost certainly dizygotic. However, if 1 corpus luteum is seen in dichorionic twins, zygosity cannot be determined with certainty since it is either monozygotic, or dizygotic where a second corpus luteum has been missed.

Genetic linkage findings for DSM-IV nicotine withdrawal in two populations.

Nicotine withdrawal (NW) is both an important contributor to difficulty quitting cigarettes and because of mood-related withdrawal symptoms a problem of particular relevance to psychiatry. Twin-studies suggest that genetic factors influence NW (heritability = 45%). Only one previous linkage study has published findings on NW [Swan et al. (2006); Am J Med Genet Part B 141B:354-360; LOD = 2.7; Chr. 6 at 159 cM]. As part of an international consortium, genome-wide scans (using over 360 autosomal microsatellite markers) and telephone diagnostic interviews were conducted on 289 Australian (AUS) and 161 Finnish (FIN, combined (COMB) N = 450 families) families ascertained from twin registries through index-cases with a lifetime history of cigarette smoking. The statistical approach used an affected-sib-pair design (at least two adult full siblings reported a history of DSM-IV NW) and conducted the linkage analyses using MERLIN. Linkage signals with LOD scores >1.5 were found on two chromosomes: 6 (FIN: LOD = 1.93 at 75 cM) and 11 at two different locations (FIN: LOD = 3.55 at 17 cM, and AUS: LOD = 1.68 with a COMB: LOD = 2.30 at 123 cM). The multipoint LOD score of 3.55 on chromosome 11p15 in FIN met genomewide significance (P = 0.013 with 1,000 simulations). At least four strong candidate genes lie within or near this peak on chromosome 11: DRD4, TPH, TH, and CHRNA10. Other studies have reported that chromosome 11 may harbor genes associated with various aspects of smoking behavior. This study adds to that literature by highlighting evidence for NW.

Evaluation of multiple displacement amplification in a 5 cM STR genome-wide scan.

Multiple displacement amplification (MDA) has emerged as a promising new method of whole genome amplification (WGA) with the potential to generate virtually unlimited genome-equivalent DNA from only a small amount of seed DNA. To date, genome-wide high marker density assessments of MDA-DNA have focussed mainly upon suitability for single nucleotide polymorphism (SNP) genotyping applications. Suitability for short tandem repeat (STR) genotyping has not been investigated in great detail, despite their inherent instability during DNA replication, and the obvious challenge that this presents to WGA techniques. Here, we aimed to assess the applicability of MDA in STR genotyping by conducting a genome-wide scan of 768 STR markers for MDAs of 15 high quality genomic DNAs. We found that MDA genotyping call and accuracy rates were only marginally lower than for genomic DNA. Pooling of three replicate MDAs resulted in a small increase in both call rate and genotyping accuracy. We identified 34 STRs (4.4% of total markers) of which five essentially failed with MDA samples, and 29 of which showed elevated genotyping failures/discrepancies in the MDAs. We emphasise the importance of DNA and MDA quality checks, and the use of appropriate controls to identify problematic STR markers.

Research funding for telemedicine: an Australian perspective.

Winning research funding is one of the most difficult challenges faced by researchers, especially with falling success rates and shrinking budgets. Telemedicine researchers can find it especially hard to win funding as they are often researching small changes to the health system that whilst important for patient care are often not as competitive as proposals that promise to cure diseases. In a climate of both tight health funding and tight research funding, telemedicine researchers should emphasise the potential for their research to add value and lower costs in order to increase their chances of winning funding.

The high costs of getting ethical and site-specific approvals for multi-centre research.

BACKGROUND: Multi-centre studies generally cost more than single-centre studies because of larger sample sizes and the need for multiple ethical approvals. Multi-centre studies include clinical trials, clinical quality registries, observational studies and implementation studies. We examined the costs of two large Australian multi-centre studies in obtaining ethical and site-specific approvals. METHODS: We collected data on staff time spent on approvals and expressed the overall cost as a percent of the total budget. RESULTS: The total costs of gaining approval were 38 % of the budget for a study of 50 centres (mean cost AUD $6960 per site) and 2 % for a study of 11 centres (mean cost AUD $2300 per site). Seventy-five and 90 % of time was spent on repeated tasks, respectively, and many time-consuming tasks, such as reformatting documents, did nothing to improve the study design or participant safety. CONCLUSIONS: Improvements have been made to the ethical approval application system, but more gains could be made without increasing risks of harm to research participants. We propose that ethical review bodies and individual sites publish statistics on how long they take to process approvals which could then be nationally benchmarked.

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.

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.

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 x 10(-20)), transferrin saturation (combined P = 2.2 x 10(-23)) and erythrocyte mean cell volume (MCV, combined P = 1.1 x 10(-10)). We also find suggestive evidence of association with blood hemoglobin levels (combined P = 5.3 x 10(-7)). These findings demonstrate the involvement of TMPRSS6 in control of iron homeostasis and in normal erythropoiesis.

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 x 10(-15)). The per allele odds ratio was 1.75 (1.53, 2.01), with evidence for stronger association in early-onset cases.