Genome-wide association study identifies novel genetic variants contributing to variation in blood metabolite levels

Metabolites are small molecules involved in cellular metabolism, which can be detected in biological samples using metabolomic techniques. Here we present the results of genome-wide association and meta-analyses for variation in the blood serum levels of 129 metabolites as measured by the Biocrates metabolomic platform. In a discovery sample of 7,478 individuals of European descent, we find 4,068 genome- and metabolome-wide significant (Z-test, P<1.09 × 10−9) associations between single-nucleotide polymorphisms (SNPs) and metabolites, involving 59 independent SNPs and 85 metabolites. Five of the fifty-nine independent SNPs are new for serum metabolite levels, and were followed-up for replication in an independent sample (N=1,182). The novel SNPs are located in or near genes encoding metabolite transporter proteins or enzymes (SLC22A16, ARG1, AGPS and ACSL1) that have demonstrated biomedical or pharmaceutical importance. The further characterization of genetic influences on metabolic phenotypes is important for progress in biological and medical research.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated approximately 2,000, approximately 3,700 and approximately 9,500 SNPs explained approximately 21%, approximately 24% and approximately 29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/beta-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Estimation and partitioning of polygenic variation captured by common SNPs for Alzheimer's disease, multiple sclerosis and endometriosis

Common diseases such as endometriosis (ED), Alzheimer's disease (AD) and multiple sclerosis (MS) account for a significant proportion of the health care burden in many countries. Genome-wide association studies (GWASs) for these diseases have identified a number of individual genetic variants contributing to the risk of those diseases. However, the effect size for most variants is small and collectively the known variants explain only a small proportion of the estimated heritability. We used a linear mixed model to fit all single nucleotide polymorphisms (SNPs) simultaneously, and estimated genetic variances on the liability scale using SNPs from GWASs in unrelated individuals for these three diseases. For each of the three diseases, case and control samples were not all genotyped in the same laboratory. We demonstrate that a careful analysis can obtain robust estimates, but also that insufficient quality control (QC) of SNPs can lead to spurious results and that too stringent QC is likely to remove real genetic signals. Our estimates show that common SNPs on commercially available genotyping chips capture significant variation contributing to liability for all three diseases. The estimated proportion of total variation tagged by all SNPs was 0.26 (SE 0.04) for ED, 0.24 (SE 0.03) for AD and 0.30 (SE 0.03) for MS. Further, we partitioned the genetic variance explained into five categories by a minor allele frequency (MAF), by chromosomes and gene annotation. We provide strong evidence that a substantial proportion of variation in liability is explained by common SNPs, and thereby give insights into the genetic architecture of the diseases.

Variation in BMPR1B, TGFRB1 and BMPR2 and control of dizygotic twinning

Genes in the TGF9 signaling pathway play important roles in the regulation of ovarian follicle growth and ovulation rate. Mutations in three genes in this pathway, growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and the bone morphogenetic protein receptor B 1 (BMPRB1), influence dizygotic (DZ) twinning rates in sheep. To date, only variants in GDF9 and BMP15, but not their receptors transforming growth factor ss receptor 1 (TGFBR1), bone morphogenetic protein receptor 2 (BMPR2) and BMPR1B, have been investigated with respect to their roles in human DZ twinning. We screened for rare and novel variants in TGFBR1, BMPR2 and BMPR1B in mothers of dizygotic twins (MODZT) from twin-dense families, and assessed association between genotyped and imputed variants and DZ twinning in another large sample of MODZT. Three novel variants were found: a deep intronic variant in BMPR2, and one intronic and one non-synonymous exonic variant in BMPRB1 which would result in the replacement of glutamine by glutamic acid at amino acid position 294 (p.Gln294Glu). None of these variants were predicted to have major impacts on gene function. However, the p.Gln294Glu variant changes the same amino acid as a sheep BMPR1B functional variant and may have functional consequences. Six BMPR1B variants were marginally associated with DZ twinning in the larger case-control sample, but these were no longer significant once multiple testing was taken into account. Our results suggest that variation in the TGF9 signaling pathway type II receptors has limited effects on DZ twinning rates in humans.

Common variation in the CYP17A1 and IFIT1 genes on chromosome 10 does not contribute to the risk of endometriosis

Endometriosis is a complex disease involving multiple susceptibility genes and environmental factors. Our previous studies on endometriosis identified a region of significant linkage on chromosome 10q. Two biological candidate genes (CYP17A1 and IFIT1) located on chromosome 10q, have previously been implicated in endometriosis and/or uterine function. We hypothesized that variation in CYP17A1 and/or IFIT1 could contribute to the risk of endometriosis and may account for some of the linkage signal on chromosome 10q. We genotyped 17 single nucleotide polymorphisms (SNPs) in the CYP17A1 and IFIT1 genes including SNP rs743572 previously associated with endometriosis in 768 endometriosis cases and 768 unrelated controls. We found no evidence for association between endometriosis and individual SNPs or SNP haplotypes in CYP17A1 and IFIT1. Common variation in these genes does not appear to be a major contributor to endometriosis susceptibility in our Australian sample.

Genetic variation in tumour necrosis factor and lymphotoxin is not associated with endometriosis in an Australian sample

BACKGROUND Tumour necrosis factor (TNF) is a pleiotropic cytokine with a wide range of immunoregulatory effects. Variation in the promoter region of TNF and the neighbouring lymphotoxin alpha (LTA) gene might be associated with endometriosis. METHODS We examined the association between endometriosis and common single-nucleotide polymorphisms (SNPs) or haplotypes in the TNF/LTA region in an Australian sample by analysing 26 SNPs in 958 endometriosis cases and 959 unrelated controls. We selected functional SNPs in the coding and the promoter region of the TNF gene and HapMap tagging SNPs and typed them on a Sequenom MassARRAY platform. A key SNP (rs1800630) in the promoter region typed in previous studies did not give reliable results. Therefore, we also examined a statistically identical (r(2) = 1) SNP (siSNP) (rs2844482), identified using the web based program ssSNPer. RESULTS Genotype completion rate was 99.5% for SNPs spanning a region of 15.5 kb across the TNF/LTA locus. There was no evidence for association between endometriosis and TNF/LTA SNPs or SNP haplotypes in our case-control study. CONCLUSIONS Our data suggest both TNF and LTA genes are not major susceptibility genes for endometriosis.

KRAS variation and risk of endometriosis

Endometriosis is a common gynaecological disease with symptoms of pelvic pain and infertility which affects 7-10% of women in their reproductive years. Activation of an oncogenic allele of Kirsten rat sarcoma viral oncogene homologue (KRAS) in the reproductive tract of mice resulted in the development of endometriosis. We hypothesized that variation in KRAS may influence risk of endometriosis in humans. Thirty tagSNPs spanning a region of 60.7 kb across the KRAS locus were genotyped using iPLEX chemistry on a MALDI-TOF MassARRAY platform in 959 endometriosis cases and 959 unrelated controls, and data were analysed for association with endometriosis. Genotypes were obtained for most individuals with a mean completion rate of 99.1%. We identified six haplotype blocks across the KRAS locus in our sample. There were no significant differences between cases and controls in the frequencies of individual single-nucleotide polymorphisms (SNPs) or haplotypes. We also developed a rapid method to screen for 11 common KRAS and BRAF mutations on the Sequenom MassARRAY system. The assay detected all mutations previously identified by direct sequencing in a panel of positive controls. No germline variants for KRAS or BRAF were detected. Our results demonstrate that any risk of endometriosis in women because of common variation in KRAS must be very small.

Variation in the walking time to bus stop by the degree of transit captivity

Measures of transit accessibility are important in evaluating transit services, planning for future services and investment on land use development. Existing tools measure transit accessibility using averaged walking distance or walking time to public transit. Although the mode captivity may have significant implications on one’s willingness to walk to use public transit, this has not been addressed in the literature to date. Failed to distinguish transit captive users may lead to overestimated ridership and spatial coverage of transit services. The aim of this research is to integrate the concept of transit captivity into the analysis of walking access to public transit. The conventional way of defining “captive” and “choice” transit users showed no significant difference in their walking times according to a preliminary analysis. A cluster analysis technique is used to further divide “choice” users by three main factors, namely age group, labour force status and personal income. After eliminating “true captive” users, defined as those without driver’s licence or without a car in respective household, “non-true captive” users were classified into a total of eight groups having similar socio-economic characteristics. The analysis revealed significant differences in the walking times and patterns by their level of captivity to public transit. This paper challenges the rule-of-thumb of 400m walking distance to bus stops. In average, people’s willingness to walk dropped drastically at 268m and continued to drop constantly until it reached the mark of 670m, where there was another drastic drop of 17%, which left with only 10% of the total bus riders willing to walk 670m or more. This research found that mothers working part time were the ones with lowest transit captivity and thus most sensitive to the walking time, followed by high-income earners and the elderly. The level of captivity increases when public transit users earned lesser income, such as students and students working part time.

Investigating variation in bus ridership rates within localised areas of a sub-tropical city

This thesis makes a significant contribution to knowledge and understanding of 'Human Travel Behaviour' in relation to transportation research. It holds some important merits that have not been proposed before. It develops a new, comprehensive and meaningful relationship that includes bus transit ridership change due to weather variables, seasonality and transit quality of service within a single daily ridership rate estimation model. The research incorporated both temporal and spatial influences on ridership within a modelling structure, named as the Nested Model Structure. It provides a complete picture of ridership variation across the sub-tropical city of Brisbane, Australia.

Diurnal variation of anterior scleral and conjunctival thickness

Purpose To examine whether anterior scleral and conjunctival thickness undergoes significant diurnal variation over a 24-hour period. Methods Nineteen healthy young adults (mean age 22 ± 2 years) with minimal refractive error (mean spherical equivalent refraction -0.08 ± 0.39 D), had measures of anterior scleral and conjunctival thickness collected using anterior segment optical coherence tomography (AS-OCT) at seven measurement sessions over a 24-hour period. The thickness of the temporal anterior sclera and conjunctiva were determined at 6 locations (each separated by 0.5 mm) at varying distances from the scleral spur for each subject at each measurement session. Results Both the anterior sclera and conjunctiva were found to undergo significant diurnal variations in thickness over a 24-hour period (both p <0.01). The sclera and conjunctiva exhibited a similar pattern of diurnal change, with a small magnitude thinning observed close to midday, and a larger magnitude thickening observed in the early morning immediately after waking. The amplitude of diurnal thickness change was larger in the conjunctiva (mean amplitude 69 ± 29 μm) compared to the sclera (21 ± 8 μm). The conjunctiva exhibited its smallest magnitude of change at the scleral spur location (mean amplitude 56 ± 17 μm) whereas the sclera exhibited its largest magnitude of change at this location (52 ± 21 μm). Conclusions This study provides the first evidence of diurnal variations occurring in the thickness of the anterior sclera and conjunctiva. Studies requiring precise measures of these anatomical layers should therefore take time of day into consideration. The majority of the observed changes occurred in the early morning immediately after waking and were of larger magnitude in the conjunctiva compared to the sclera. Thickness changes at other times of the day were of smaller magnitude and generally not statistically significant.

HR-MAS NMR spectroscopy shows regional metabolite variation in bovine articular cartilage

Mechanical stress is an important external factor effecting the development and maintenance of articular cartilage. The metabolite profile of diseased cartilage has been well studied but there is limited information about the variation in metabolite profile of healthy cartilage. With the importance of load in maintaining healthy cartilage, regional differences in metabolite profile associated with differences in load may provide information on how load contributes to the maintenance of healthy cartilage. HR-MAS NMR spectroscopy allows the assessment of tissue samples without modification and was used for assessing the difference in metabolic profile between the load bearing and non-load bearing regions of the bovine articular cartilage. In this preliminary study, we examined cartilage from tibia and femur of four knee joints. Sixteen pairs of 1D-NOESY spectra were acquired. Principle component analysis (PCA) identified chemical shifts responsible for variance. SBASE (AMIX) and the Human Metabolome Database were used in conjunction with previous reported cartilage data for identifying metabolites associated with the PCA results. The major contributors to load-related differences in metabolite profile were N-acetyl groups, lactate and phosphocholine peaks. Integrals of these regions were further analysed using a Student's t-test. In load bearing cartilage regions. N-acetyl groups and phosphocholine were found at significantly higher concentration (p < 0.05 and p < 0.005, respectively) in both femur and tibia, while lactate was reduced in load bearing cartilage (p < 0.005). The results of this pilot HR-MAS NMR study demonstrate its ability to provide useful metabolite information for healthy cartilage.

Natural variation and Biogeography of the melon fruit fly, Zeugodacus cucurbitae (Diptera : Tephritidae), in Southeast-Asia and the west-pacific

This thesis used multidisciplinary approaches which greatly enhance our understanding of population structure and can be particularly powerful tools for resolving variation of melon fly over geographic and temporal scales, and for determining invasive pathways. The results from this thesis reinforce the value of integrating multiple data sets to better understand and resolve natural variation within an important pest to determine whether there are cryptic species, discrete lineages or host races, and to identify dispersal pathways in an invasive pest. These results are instructive for regional biosecurity, trade and quarantine, and provide important background for future area-wide management programmes. The integrative methodology adopted in this thesis is applicable to a variety of other insect pests.

Regional variation in social isolation amongst older Australians

Regional studies globally have a strong focus on understanding the causes of variation in the economic performance and wellbeing of regions and this emphasis acknowledges that the strength of the local or regional economy plays a determinant role in shaping quality of life. Regional research has been less active in considering spatial variation in other factors that are critical to individual and societal wellbeing. For example, the regional studies community has been absent from the debate on the social determinants of health and how these influences vary spatially. This paper considers the results of a cross sectional survey of Australians aged 65 and over that focussed on social connections and wellbeing. It examines regional variations in the incidence of social isolation within the older population. It finds that while the incidence of self-reported social isolation amongst older persons is broadly consistent with earlier studies, it demonstrates a spatial patterning that is unexpected. The paper considers community-building activities in addressing the impacts of social isolation, including the role of urban design, and suggests that there is a need to supplement the national overview presented there through more detailed studies focussed on individual localities.

Chromosomal rearrangements, phenotypic variation and modularity: A case study from a contact zone between house mouse Robertsonian races in Central Italy

The Western European house mouse, Mus musculus domesticus, is well-known for the high frequency of Robertsonian fusions that have rapidly produced more than 50 karyotipic races, making it an ideal model for studying the mechanisms of chromosomal speciation. The mouse mandible is one of the traits studied most intensively to investigate the effect of Robertsonian fusions on phenotypic variation within and between populations. This complex bone structure has also been widely used to study the level of integration between different morphogenetic units. Here, with the aim of testing the effect of different karyotypic assets on the morphology of the mouse mandible and on its level of modularity, we performed morphometric analyses of mice from a contact area between two highly metacentric races in Central Italy. We found no difference in size, while the mandible shape was found to be different between the two Robertsonian races, even after accounting for the genetic relationships among individuals and geographic proximity. Our results support the existence of two modules that indicate a certain degree of evolutionary independence, but no difference in the strength of modularity between chromosomal races. Moreover, the ascending ramus showed more pronounced interpopulation/race phenotypic differences than the alveolar region, an effect that could be associated to their different polygenic architecture. This study suggests that chromosomal rearrangements play a role in the house mouse phenotypic divergence, and that the two modules of the mouse mandible are differentially affected by environmental factors and genetic makeup.