Educational attainment: A genome wide association study in 9538 Australians

BACKGROUND Correlations between Educational Attainment (EA) and measures of cognitive performance are as high as 0.8. This makes EA an attractive alternative phenotype for studies wishing to map genes affecting cognition due to the ease of collecting EA data compared to other cognitive phenotypes such as IQ. METHODOLOGY In an Australian family sample of 9538 individuals we performed a genome-wide association scan (GWAS) using the imputed genotypes of approximately 2.4 million single nucleotide polymorphisms (SNP) for a 6-point scale measure of EA. Top hits were checked for replication in an independent sample of 968 individuals. A gene-based test of association was then applied to the GWAS results. Additionally we performed prediction analyses using the GWAS results from our discovery sample to assess the percentage of EA and full scale IQ variance explained by the predicted scores. RESULTS The best SNP fell short of having a genome-wide significant p-value (p = 9.77x10(-7)). In our independent replication sample six SNPs among the top 50 hits pruned for linkage disequilibrium (r(2)<0.8) had a p-value<0.05 but only one of these SNPs survived correction for multiple testing--rs7106258 (p = 9.7*10(-4)) located in an intergenic region of chromosome 11q14.1. The gene based test results were non-significant and our prediction analyses show that the predicted scores explained little variance in EA in our replication sample. CONCLUSION While we have identified a polymorphism chromosome 11q14.1 associated with EA, further replication is warranted. Overall, the absence of genome-wide significant p-values in our large discovery sample confirmed the high polygenic architecture of EA. Only the assembly of large samples or meta-analytic efforts will be able to assess the implication of common DNA polymorphisms in the etiology of EA.

Modern diagnosis of Epstein-Barr virus infections and post-transplant lymphoproliferative disease

Infection by Epstein-Barr virus (EBV) occurs in approximately 95% of the world s population. EBV was the first human virus implicated in oncogenesis. Characteristic for EBV primary infection are detectable IgM and IgG antibodies against viral capsid antigen (VCA). During convalescence the VCA IgM disappears while the VCA IgG persists for life. Reactivations of EBV occur both among immunocompromised and immunocompetent individuals. In serological diagnosis, measurement of avidity of VCA IgG separates primary from secondary infections. However, in serodiagnosis of mononucleosis it is quite common to encounter, paradoxically, VCA IgM together with high-avidity VCA IgG, indicating past immunity. We determined the etiology of this phenomenon and found that, among patients with cytomegalovirus (CMV) primary infection a large proportion (23%) showed antibody profiles of EBV reactivation. In contrast, EBV primary infection did not appear to induce immunoreactivation of CMV. EBV-associated post-transplant lymphoproliferative disease (PTLD) is a life threatening complication of allogeneic stem cell or solid organ transplantation. PTLD may present with a diverse spectrum of clinical symptoms and signs. Due to rapidity of PTLD progression especially after stem cell transplantation, the diagnosis must be obtained quickly. Pending timely detection, the evolution of the fatal disease may be halted by reduction of immunosuppression. A promising new PTLD treatment (also in Finland) is based on anti-CD-20 monoclonal antibodies. Diagnosis of PTLD has been demanding because of immunosuppression, blood transfusions and the latent nature of the virus. We set up in 1999 to our knowledge first in Finland for any microbial pathogen a real-time quantitative PCR (qPCR) for detection of EBV DNA in blood serum/plasma. In addition, we set up an in situ hybridisation assay for EBV RNA in tissue sections. In collaboration with a group of haematologists at Helsinki University Central Hospital we retrospectively determined the incidence of PTLD among 257 allogenic stem cell transplantations (SCT) performed during 1994-1999. Post-mortem analysis revealed 18 cases of PTLD. From a subset of PTLD cases (12/18) and a series of corresponding controls (36), consecutive samples of serum were studied by the new EBV-qPCR. All the PTLD patients were positive for EBV-DNA with progressively rising copy numbers. In most PTLD patients EBV DNA became detectable within 70 days of SCT. Of note, the appearance of EBV DNA preceded the PTLD symptoms (fever, lymphadenopathy, atypical lymphocytes). Among the SCT controls, EBV DNA occurred only sporadically, and the EBV-DNA levels remained relatively low. We concluded that EBV qPCR is a highly sensitive (100%) and specific (96%) new diagnostic approach. We also looked for and found risk factors for the development of PTLD. Together with a liver transplantation group at the Transplantation and Liver Surgery Clinic we wanted to clarify how often and how severely do EBV infections occur after liver transplantation. We studied by the EBV qPCR 1284 plasma samples obtained from 105 adult liver transplant recipients. EBV DNA was detected in 14 patients (13%) during the first 12 months. The peak viral loads of 13 asymptomatic patients were relatively low (<6600/ml), and EBV DNA subsided quickly from circulation. Fatal PTLD was diagnosed in one patient. Finally, we wanted to determine the number and clinical significance of EBV infections of various types occurring among a large, retrospective, nonselected cohort of allogenic SCT recipients. We analysed by EBV qPCR 5479 serum samples of 406 SCT recipients obtained during 1988-1999. EBV DNA was seen in 57 (14%) patients, of whom 22 (5%) showed progressively rising and ultimately high levels of EBV DNA (median 54 million /ml). Among the SCT survivors, EBV DNA was transiently detectable in 19 (5%) asymptomatic patients. Thereby, low-level EBV-DNA positivity in serum occurs relatively often after SCT and may subside without specific treatment. However, high molecular copy numbers (>50 000) are diagnostic for life-threatening EBV infection. We furthermore developed a mathematical algorithm for the prediction of development of life-threatening EBV infection.

GWAS of butyrylcholinesterase activity identifies four novel loci, independent effects within BCHE and secondary associations with metabolic risk factors

Serum butyrylcholinesterase (BCHE) activity is associated with obesity, blood pressure and biomarkers of cardiovascular and diabetes risk. We have conducted a genome-wide association scan to discover genetic variants affecting BCHE activity, and to clarify whether the associations between BCHE activity and cardiometabolic risk factors are caused by variation in BCHE or whether BCHE variation is secondary to the metabolic abnormalities. We measured serum BCHE in adolescents and adults from three cohorts of Australian twin and family studies. The genotypes from approximately 2.4 million single-nucleotide polymorphisms (SNPs) were available in 8791 participants with BCHE measurements. We detected significant associations with BCHE activity at three independent groups of SNPs at the BCHE locus (P = 5.8 x 10(-262), 7.8 x 10(-47), 2.9 x 10(-12)) and at four other loci: RNPEP (P = 9.4 x 10(-16)), RAPH1-ABI2 (P = 4.1 x 10(-18)), UGT1A1 (P = 4.0 x 10(-8)) and an intergenic region on chromosome 8 (P = 1.4 x 10(-8)). These loci affecting BCHE activity were not associated with metabolic risk factors. On the other hand, SNPs in genes previously associated with metabolic risk had effects on BCHE activity more often than can be explained by chance. In particular, SNPs within FTO and GCKR were associated with BCHE activity, but their effects were partly mediated by body mass index and triglycerides, respectively. We conclude that variation in BCHE activity is due to multiple variants across the spectrum from uncommon/large effect to common/small effect, and partly results from (rather than causes) metabolic abnormalities.

Insights into the molecular genetics of celiac disease : Applying family- and population-based strategies

Celiac disease, or gluten intolerance, is triggered by dietary glutens in genetically susceptible individuals and it affects approximately 1% of the Caucasian population. The best known genetic risk factors for celiac disease are HLA DQ2 and DQ8 heterodimers, which are necessary for the development of the disease. However, they alone are not sufficient for disease induction, other risk factors are required. This thesis investigated genetic factors for celiac disease, concentrating on susceptibility loci on chromosomes 5q31-q33, 19p13 and 2q12 previously reported in genome-wide linkage and association studies. In addition, a novel genotyping method for the detection of HLA DQ2 and DQ8 coding haplotypes was validated. This study was conducted using Finnish and Hungarian family materials, and Finnish, Hungarian and Italian case-control materials. Genetic linkage and association were analysed in these materials using candidate gene and fine-mapping approaches. The results confirmed linkage to celiac disease on the chromosomal regions 5q31-q33 and 19p13. Fine-mapping on chromosome 5q31-q33 revealed several modest associations in the region, and highlighted the need for further investigations to locate the causal risk variants. The MYO9B gene on chromosome 19p13 showed evidence for linkage and association particularly with dermatitis herpetiformis, the skin manifestation of celiac disease. This implies a potential difference in the genetic background of the intestinal and skin forms of the disease, although studies on larger samplesets are required. The IL18RAP locus on chromosome 2q12, shown to be associated with celiac disease in a previous genome-wide association study and a subsequent follow-up, showed association in the Hungarian population in this study. The expression of IL18RAP was further investigated in small intestinal tissue and in peripheral blood mononuclear cells. The results showed that IL18RAP is expressed in the relevant tissues. Two putative isoforms of IL18RAP were detected by Western blot analysis, and the results suggested that the ratios and total levels of these isoforms may contribute to the aetiology of celiac disease. A novel genotyping method for celiac disease-associated HLA haplotypes was also validated in this thesis. The method utilises single-nucleotide polymorphisms tagging these HLA haplotypes with high sensitivity and specificity. Our results suggest that this method is transferable between populations, and it is suitable for large-scale analysis. In conclusion, this doctorate study provides an insight into the roles of the 5q31-q33, MYO9B, IL18RAP and HLA loci in the susceptibility to celiac disease in the Finnish, Hungarian and Italian populations, highlighting the need for further studies at these genetic loci and examination of the function of the candidate genes.

Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and approximately 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 x 10(-)(8)), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

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.

A simple and fast two-locus quality control test to detect false positives due to batch effects in genome-wide association studies

The impact of erroneous genotypes having passed standard quality control (QC) can be severe in genome-wide association studies, genotype imputation, and estimation of heritability and prediction of genetic risk based on single nucleotide polymorphisms (SNP). To detect such genotyping errors, a simple two-locus QC method, based on the difference in test statistic of association between single SNPs and pairs of SNPs, was developed and applied. The proposed approach could detect many problematic SNPs with statistical significance even when standard single SNP QC analyses fail to detect them in real data. Depending on the data set used, the number of erroneous SNPs that were not filtered out by standard single SNP QC but detected by the proposed approach varied from a few hundred to thousands. Using simulated data, it was shown that the proposed method was powerful and performed better than other tested existing methods. The power of the proposed approach to detect erroneous genotypes was approximately 80% for a 3% error rate per SNP. This novel QC approach is easy to implement and computationally efficient, and can lead to a better quality of genotypes for subsequent genotype-phenotype investigations.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

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.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Integrating sorghum whole genome sequence information with a compendium of sorghum QTL studies reveals uneven distribution of QTL and of gene-rich regions with significant implications for crop improvement.

A comprehensive analysis was conducted using 48 sorghum QTL studies published from 1995 to 2010 to make information from historical sorghum QTL experiments available in a form that could be more readily used by sorghum researchers and plant breeders. In total, 771 QTL relating to 161 unique traits from 44 studies were projected onto a sorghum consensus map. Confidence intervals (CI) of QTL were estimated so that valid comparisons could be made between studies. The method accounted for the number of lines used and the phenotypic variation explained by individual QTL from each study. In addition, estimated centimorgan (cM) locations were calculated for the predicted sorghum gene models identified in Phytozome (JGI GeneModels SBI v1.4) and compared with QTL distribution genome-wide, both on genetic linkage (cM) and physical (base-pair/bp) map scales. QTL and genes were distributed unevenly across the genome. Heterochromatic enrichment for QTL was observed, with approximately 22% of QTL either entirely or partially located in the heterochromatic regions. Heterochromatic gene enrichment was also observed based on their predicted cM locations on the sorghum consensus map, due to suppressed recombination in heterochromatic regions, in contrast to the euchromatic gene enrichment observed on the physical, sequence-based map. The finding of high gene density in recombination-poor regions, coupled with the association with increased QTL density, has implications for the development of more efficient breeding systems in sorghum to better exploit heterosis. The projected QTL information described, combined with the physical locations of sorghum sequence-based markers and predicted gene models, provides sorghum researchers with a useful resource for more detailed analysis of traits and development of efficient marker-assisted breeding strategies.

Genetic studies on recurrent miscarriage

Recurrent miscarriage (RM) is defined as three consecutive pregnancy failures and is estimated to affect ~1% of couples trying to conceive. The cause of RM remains unknown in approximately 50% of cases. In this study, it was hypothesized that some of the underlying factors yet to be discovered are genetic. The aim was to search for mutations in genes AMN, EPCR, TM, and p53 known to cause miscarriage in mouse models and thereby find new genetic causes for unexplained miscarriages in humans. In addition, the mitochondrial genome was studied because mitochondria are involved in processes important in early development. Furthermore, sex chromosome characteristics suggested to underlie miscarriage were also studied. A total of 40 couples and 8 women with unexplained RM were collected for this study and screened for mutations in the candidate genes. Six interesting exonic or potential splice site disrupting variations were detected. However, their phenotypic effects cannot be determined without further investigations. Additionally, an association between the C11992A polymorphism of the p53 gene and RM was detected. The results indicate that women carrying the C/A or A/A genotype have a two-fold higher risk for RM than women with a C/C genotype. This strengthens the results of previous studies reporting that p53 sequence variations may cause miscarriage. The role of variation C11992A in embryonic development is, however, difficult to predict without further studies When screening the mitochondrial genome a heteroplasmic mtDNA variation was found in an unexpected high number of women, as heteroplasmic variations are reported to be rare. One novel variation and 18 previously reported polymorphisms were detected in the mitochondrial genome. Although the detected variations are likely to be neutral polymorphisms, a role in the aetiology of miscarriage cannot be excluded as some mtDNA variations may be pathogenic only when a threshold is reached. Recent publications have reported skewed X chromosome inactivation and Y chromosome microdeletions to be associated with RM. Therefore, these sex chromosome abnormalities in the context of RM were investigated. No associations between skewed X chromosome inactivation or Y chromosome microdeletions and RM in the Finnish patients were detected. Data on ancestral birthplaces of the patients were collected to study any possible geographic clustering, which would indicate a common predisposing factor. The results showed clustering of the birthplaces in eastern Finland in a subset of patients. This suggests a possibility of an enriched susceptibility gene which may contribute to RM.

Economic modelling of grazing systems in the Fitzroy and Burdekin catchments

The Great Barrier Reef (GBR) is the largest reef system in the world; it covers an area of approximately 2,225,000 km² in the northern Queensland continental shelf. There are approximately 750 reefs that exist within 40 km of the Queensland coast. Recent research has identified that poor water quality is having negative impacts on the GBR (Haynes et al. 2007). The Fitzroy Basin covers 143,000 km² and is the largest catchment draining into the GBR as well as being one of the largest catchments in Australia (Karfs et al. 2009). The Burdekin Catchment is the second largest catchment entering into the GBR and covers 133,432 km².The prime determinant for the changes in water quality entering into the GBR have been attributed to grazing, with beef production the largest single land use industry comprising 90% of the land area (Karfs et al. 2009). Extensive beef production contributes over $1 billion dollars to the national economy annually and employs over 9000 people, many in rural communities (Gordon 2007). ‘Economic modelling of grazing systems in the Fitzroy and Burdekin catchments’ was a joint project with the Fitzroy Basin Association and the Queensland Department of Employment Economic Development and Innovation. The project was formed under the federally funded Caring For Our Country and the Reef Rescue programs. The project objectives were as follows; * Quantifying the costs of over-utilising available pasture and the resulting sediment leaving a representative farm for four of the major land systems in the Burdekin or Fitzroy catchments and identifying economically optimal pasture utilisation rates * Estimating the cost of reducing pasture utilisation rates below the determined optimal * Using this information, guide the selection of appropriate tools to achieve reduced utilisation rates e.g. extension process versus incentive payments or a combination of both * Model the biophysical and economic impacts of altering grazing systems to restore land condition e.g. from C condition to B condition for four land systems in the Burdekin or Fitzroy catchments.

Enhancing economic input to the CQSS2 Project report. Commissioned by the Fitzroy Basin Association.

The Fitzroy Basin is the second largest catchment area in Australia covering 143,00 km² and is the largest catchment for the Great Barrier Reef lagoon (Karfs et al., 2009). The Great Barrier Reef is the largest reef system in the world; it covers an area of approximately 225,000 km² in the northern Queensland continental shelf. There are approximately 750 reefs that exist within 40 km of the Queensland Coast (Haynes et al., 2007). The prime determinant for the changes in water quality have been attributed to grazing, with beef production the largest single land use industry comprising 90% of the land area (Karfs et al., 2009). In response to the depletion of water quality in the reef, in 2003 a Reef Water Quality plan was developed by the Australian and Queensland governments. The plan targets as a priority sediment contributions from grazing cattle in high risk catchments (The State of Queensland and Commonwealth of Australia, 2003). The economic incentive strategy designed includes analysing the costs and benefits of best management practice that will lead to improved water quality (The State of Queensland and Commonwealth of Australia, 2003).

Dengue virus infection : Diagnostics and molecular epidemiology

Dengue is a mosquito-borne viral disease caused by the four dengue virus serotypes (DENV-1-4) and is currently considered as the most important arthropod-borne viral disease in the world. Nearly half of the human population lives in risk areas, and 50-100 million infections occur yearly according to World Health Organization. The disease can vary from a mild febrile disease to severe haemorrhagic fever and shock. A secondary infection with heterologous serotype increases the risk for severe disease outcome. During the last three decades the impact of dengue has dramatically increased in the endemic areas including the tropics and subtropics of the world. The current situation with massive epidemics of severe disease forms has been associated with socio-ecological changes that have increased the transmission and enabled the co-circulation of different serotypes. Consequently, an increase of dengue has also been observed in travelers visiting these areas. Currently approximately 30 cases are diagnosed yearly in Finnish travelers. In travelers dengue is rarely a life-threatening disease, however in the current study, a fatality was documented in a young Finnish patient who experienced a prolonged primary dengue infection. To improve particularly early laboratory diagnostics, a novel real-time RT-PCR method was developed for the detection of DENV-1-4 RNA based on TaqMan chemistry. The method was shown to be sensitive and specific for detecting DENV RNA and suitable for diagnostic use. The newly developed real-time RT-PCR was compared to other available early diagnostic methods including IgM and NS1 antigen detection using a panel of selected patient samples. The results suggest that the best diagnostic rates are achieved by a combination of IgM with RNA or NS1 detection. The dengue virus strains studied here included the first DENV strains isolated from serum samples of Finnish travelers collected in 2000-2005. The results of sequence analysis demonstrated that the 11 isolates included all four DENV serotypes and presented a global sample of DENV strains from different geographical areas including Asia, Africa and South America. In the present study sequence analysis was also carried out for a collection of 23 novel DENV-2 isolates from Venezuelan patients collected in 1999-2005. The Venezuelan DENV-2 exclusively represented the American-Asian genotype, suggesting that no foreign DENV-2 lineages have recently been introduced to the country. The results also suggest that the DENV-2 viruses detected earlier from Venezuela have been maintained in the area where they have evolved into several lineages. This is in contrast to the pattern observed in some other dengue endemic areas, where introductions of novel virus types and lineages are frequently detected.