A versatile gene-based test for genome-wide association studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

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.

Variants in EMX2 and PTEN do not contribute to risk of endometriosis

Endometriosis has a genetic component, and significant linkage has been found to a region on chromosome 10q. Two candidate genes, EMX2 and PTEN, implicated in both endometriosis and endometrial cancer, lie on chromosome 10q. We hypothesized that variation in EMX2 and/or PTEN could contribute to the risk of endometriosis and may account for some of the linkage signal on 10q. We genotyped single nucleotide polymorphisms (SNPs) in a case-control design to evaluate association between endometriosis and common variations in these two genes. The genotyping and statistical analysis were based on samples collected from Australian volunteers. The cases were 768 unrelated women with surgically confirmed endometriosis selected from affected sister pair (ASP) families participating in the Australian Genes behind Endometriosis Study. The controls were 768 female participants in twin studies who, based on screening questions, did not have a diagnosis of endometriosis. Genotypes of 22 SNPs in the EMX2 gene and 15 SNPs in the PTEN gene were the main outcome measures. Statistical analysis provided measures of linkage disequilibrium and association. Permutation testing showed no globally significant association between any SNPs or haplotypes and endometriosis for either gene. It is unlikely that the EMX2 or PTEN gene variants investigated contribute to risk for initiation and/or development of endometriosis.

ssSNPer: identifying statistically similar SNPs to aid interpretation of genetic association studies

ssSNPer is a novel user-friendly web interface that provides easy determination of the number and location of untested HapMap SNPs, in the region surrounding a tested HapMap SNP, which are statistically similar and would thus produce comparable and perhaps more significant association results. Identification of ssSNPs can have crucial implications for the interpretation of the initial association results and the design of follow-up studies. AVAILABILITY: http://fraser.qimr.edu.au/general/daleN/ssSNPer/

Association study of the dystrobrevin-binding gene with schizophrenia in Australian and Indian samples

Numerous studies have reported association between variants in the dystrobrevin binding protein 1 (dysbindin) gene (DTNBP1) and schizophrenia. However, the pattern of results is complex and to date, no specific risk marker or haplotype has been consistently identified. The number of single nucleotide polymorphisms (SNPs) tested in these studies has ranged from 5 to 20. We attempted to replicate previous findings by testing 16 SNPs in samples of 41 Australian pedigrees, 194 Australian cases and 180 controls, and 197 Indian pedigrees. No globally significant evidence for association was observed in any sample, despite power calculations indicating sufficient power to replicate several previous findings. Possible explanations for our results include sample differences in background linkage disequilibrium and/or risk allele effect size, the presence of multiple risk alleles upon different haplotypes, or the presence of a single risk allele upon multiple haplotypes. Some previous associations may also represent false positives. Examination of Caucasian HapMap phase II genotype data spanning the DTNBP1 region indicates upwards of 40 SNPs are required to satisfactorily assess all nonredundant variation within DTNBP1 and its potential regulatory regions for association with schizophrenia. More comprehensive studies in multiple samples will be required to determine whether specific DTNBP1 variants function as risk factors for schizophrenia.

Separate and interacting effects within the catechol-O-methyltransferase (COMT) are associated with schizophrenia

Several lines of evidence have implicated the catechol-O-methyltransferase (COMT) gene as a candidate for schizophrenia (SZ) susceptibility, not only because it encodes a key dopamine catabolic enzyme but also because it maps to the velocardiofacial syndrome region of chromosome 22q11 which has long been associated with SZ predisposition. The interest in COMT as a candidate SZ risk factor has led to numerous case-control and family-based studies, with the majority placing emphasis on examining a functional Val/Met polymorphism within this enzyme. Unfortunately, these studies have continually produced conflicting results. To assess the genetic contribution of other COMT variants to SZ susceptibility, we investigated three single-nucleotide polymorphisms (SNPs) (rs737865, rs4633, rs165599) in addition to the Val/Met variant (rs4680) in a highly selected sample of Australian Caucasian families containing 107 patients with SZ. The Val/Met and rs4633 variants showed nominally significant associations with SZ (P<0.05), although neither of the individual SNPs remained significant after adjusting for multiple testing (most significant P=0.1174). However, haplotype analyses showed strong evidence of an association; the most significant being the three-marker haplotype rs737865-rs4680-rs165599 (global P=0.0022), which spans more than 26 kb. Importantly, conditional analyses indicated the presence of two separate and interacting effects within this haplotype, irrespective of gender. In addition, our results indicate the Val/Met polymorphism is not disease-causing and is simply in strong linkage disequilibrium with a causative effect, which interacts with another as yet unidentified variant approximately 20 kb away. These results may help explain the inconsistent results reported on the Val/Met polymorphism and have important implications for future investigations into the role of COMT in SZ susceptibility.

GENEHUNTER: your 'one-stop shop' for statistical genetic analysis?

The past decade has brought a proliferation of statistical genetic (linkage) analysis techniques, incorporating new methodology and/or improvement of existing methodology in gene mapping, specifically targeted towards the localization of genes underlying complex disorders. Most of these techniques have been implemented in user-friendly programs and made freely available to the genetics community. Although certain packages may be more 'popular' than others, a common question asked by genetic researchers is 'which program is best for me?'. To help researchers answer this question, the following software review aims to summarize the main advantages and disadvantages of the popular GENEHUNTER package.

Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

Quantitative trait loci detection and benefits from marker-assisted selection in dairy cattle

Knowing the chromosomal areas or actual genes affecting the traits under selection would add more information to be used in the selection decisions which would potentially lead to higher genetic response. The first objective of this study was to map quantitative trait loci (QTL) affecting economically important traits in the Finnish Ayrshire population. The second objective was to investigate the effects of using QTL information in marker-assisted selection (MAS) on the genetic response and the linkage disequilibrium between the different parts of the genome. Whole genome scans were carried out on a grand-daughter design with 12 half-sib families and a total of 493 sons. Twelve different traits were studied: milk yield, protein yield, protein content, fat yield, fat content, somatic cell score (SCS), mastitis treatments, other veterinary treatments, days open, fertility treatments, non-return rate, and calf mortality. The average spacing of the typed markers was 20 cM with 2 to 14 markers per chromosome. Associations between markers and traits were analyzed with multiple marker regression. Significance was determined by permutation and genome-wise P-values obtained by Bonferroni correction. The benefits from MAS were investigated by simulation: a conventional progeny testing scheme was compared to a scheme where QTL information was used within families to select among full-sibs in the male path. Two QTL on different chromosomes were modelled. The effects of different starting frequencies of the favourable alleles and different size of the QTL effects were evaluated. A large number of QTL, 48 in total, were detected at 5% or higher chromosome-wise significance. QTL for milk production were found on 8 chromosomes, for SCS on 6, for mastitis treatments on 1, for other veterinary treatments on 5, for days open on 7, for fertility treatments on 7, for calf mortality on 6, and for non-return rate on 2 chromosomes. In the simulation study the total genetic response was faster with MAS than with conventional selection and the advantage of MAS persisted over the studied generations. The rate of response and the difference between the selection schemes reflected clearly the changes in allele frequencies of the favourable QTL. The disequilibrium between the polygenes and QTL was always negative and it was larger with larger QTL size. The disequilibrium between the two QTL was larger with QTL of large effect and it was somewhat larger with MAS for scenarios with starting frequencies below 0.5 for QTL of moderate size and below 0.3 for large QTL. In conclusion, several QTL affecting economically important traits of dairy cattle were detected. Further studies are needed to verify these QTL, check their presence in the present breeding population, look for pleiotropy and fine map the most interesting QTL regions. The results of the simulation studies show that using MAS together with embryo transfer to pre-select young bulls within families is a useful approach to increase the genetic merit of the AI-bulls compared to conventional selection.

Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

NEESTIMATOR v2: re-implementation of software for the estimation of contemporary effective population size (N-e) from genetic data

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.