Mapping quantitative trait loci in a wild population using linkage and linkage disequilibrium analyses

Historical information can be used, in addition to pedigree, traits and genotypes, to map quantitative trait locus (QTL) in general populations via maximum likelihood estimation of variance components. This analysis is known as linkage disequilibrium (LD) and linkage mapping, because it exploits both linkage in families and LD at the population level. The search for QTL in the wild population of Soay sheep on St. Kilda is a proof of principle. We analysed the data from a previous study and confirmed some of the QTLs reported. The most striking result was the confirmation of a QTL affecting birth weight that had been reported using association tests but not when using linkage-based analyses. Copyright © Cambridge University Press 2010.

Evolutionary genetics and human assisted movement of a globally invasive pest (Russian wheat aphid : Diuraphis noxia)

This PhD study has examined the population genetics of the Russian wheat aphid (RWA, Diuraphis noxia), one of the world’s most invasive agricultural pests, throughout its native and introduced global range. Firstly, this study investigated the geographic distribution of genetic diversity within and among RWA populations in western China. Analysis of mitochondrial data from 18 sites provided evidence for the long-term existence and expansion of RWAs in western China. The results refute the hypothesis that RWA is an exotic species only present in China since 1975. The estimated date of RWA expansion throughout western China coincides with the debut of wheat domestication and cultivation practices in western Asia in the Holocene. It is concluded that western China represents the limit of the far eastern native range of this species. Analysis of microsatellite data indicated high contemporary gene flow among northern populations in western China, while clear geographic isolation between northern and southern populations was identified across the Tianshan mountain range and extensive desert regions. Secondly, this study analyzed the worldwide pathway of invasion using both microsatellite and endosymbiont genetic data. Individual RWAs were obtained from native populations in Central Asia and the Middle East and invasive populations in Africa and the Americas. Results indicated two pathways of RWA invasion from 1) Syria in the Middle East to North Africa and 2) Turkey to South Africa, Mexico and then North and South America. Very little clone diversity was identified among invasive populations suggesting that a limited founder event occurred together with predominantly asexual reproduction and rapid population expansion. The most likely explanation for the rapid spread (within two years) from South Africa to the New World is by human movement, probably as a result of the transfer of wheat breeding material. Furthermore, the mitochondrial data revealed the presence of a universal haplotype and it is proposed that this haplotype is representative of a wheat associated super-clone that has gained dominance worldwide as a result of the widespread planting of domesticated wheat. Finally, this study examined salivary gland gene diversity to determine whether a functional basis for RWA invasiveness could be identified. Peroxidase DNA sequence data were obtained for a selection of worldwide RWA samples. Results demonstrated that most native populations were polymorphic while invasive populations were monomorphic, supporting previous conclusions relating to demographic founder effects in invasive populations. Purifying selection most likely explains the existence of a universal allele present in Middle Eastern populations, while balancing selection was evident in East Asian populations. Selection acting on the peroxidase gene may provide an allele-dependent advantage linked to the successful establishment of RWAs on wheat, and ultimately their invasion potential. In conclusion, this study is the most comprehensive molecular genetic investigation of RWA population genetics undertaken to date and provides significant insights into the source and pathway of global invasion and the potential existence of a wheat-adapted genotype that has colonised major wheat growing countries worldwide except for Australia. This research has major biosecurity implications for Australia’s grain industry.

Complete mitochondrial genome sequencing reveals novel haplotypes in a Polynesian population

The high risk of metabolic disease traits in Polynesians may be partly explained by elevated prevalence of genetic variants involved in energy metabolism. The genetics of Polynesian populations has been shaped by island hoping migration events which have possibly favoured thrifty genes. The aim of this study was to sequence the mitochondrial genome in a group of Maoris in an effort to characterise genome variation in this Polynesian population for use in future disease association studies. We sequenced the complete mitochondrial genomes of 20 non-admixed Maori subjects using Affymetrix technology. DNA diversity analyses showed the Maori group exhibited reduced mitochondrial genome diversity compared to other worldwide populations, which is consistent with historical bottleneck and founder effects. Global phylogenetic analysis positioned these Maori subjects specifically within mitochondrial haplogroup - B4a1a1. Interestingly, we identified several novel variants that collectively form new and unique Maori motifs – B4a1a1c, B4a1a1a3 and B4a1a1a5. Compared to ancestral populations we observed an increased frequency of non-synonymous coding variants of several mitochondrial genes in the Maori group, which may be a result of positive selection and/or genetic drift effects. In conclusion, this study reports the first complete mitochondrial genome sequence data for a Maori population. Overall, these new data reveal novel mitochondrial genome signatures in this Polynesian population and enhance the phylogenetic picture of maternal ancestry in Oceania. The increased frequency of several mitochondrial coding variants makes them good candidates for future studies aimed at assessment of metabolic disease risk in Polynesian populations.

Genome-wide association study reveals three susceptibility loci for common migraine in the general population

Migraine is a common, heterogeneous and heritable neurological disorder. Its pathophysiology is incompletely understood, and its genetic influences at the population level are unknown. In a population-based genome-wide analysis including 5,122 migraineurs and 18,108 non-migraineurs, rs2651899 (1p36.32, PRDM16), rs10166942 (2q37.1, TRPM8) and rs11172113 (12q13.3, LRP1) were among the top seven associations (P < 5 × 10(-6)) with migraine. These SNPs were significant in a meta-analysis among three replication cohorts and met genome-wide significance in a meta-analysis combining the discovery and replication cohorts (rs2651899, odds ratio (OR) = 1.11, P = 3.8 × 10(-9); rs10166942, OR = 0.85, P = 5.5 × 10(-12); and rs11172113, OR = 0.90, P = 4.3 × 10(-9)). The associations at rs2651899 and rs10166942 were specific for migraine compared with non-migraine headache. None of the three SNP associations was preferential for migraine with aura or without aura, nor were any associations specific for migraine features. TRPM8 has been the focus of neuropathic pain models, whereas LRP1 modulates neuronal glutamate signaling, plausibly linking both genes to migraine pathophysiology.

Investigation of the 1758G>C and 2880A>G variants within the NCOA3 gene in a breast cancer affected Australian population

NCOA3 is a known low to moderate-risk breast cancer susceptibility gene, amplified in 5–10% and over expressed in about 60% of breast tumours. Additionally, this over expression is associated with Tamoxifen resistance and poor prognosis. Previously, two variants of NCOA3, 1758G > C and 2880A > G have been associated with breast cancer in two independent populations. Here we assessed the influence of the two NCOA3 variants on breast cancer risk by genotyping an Australian case–control study population. 172 cases and 178 controls were successfully genotyped for the 1758G > C variant and 186 cases and 182 controls were successfully genotyped for the 2880A > G variant using high-resolution melt analysis (HRM). The genotypes of the 1758G > C variant were validated by sequencing. χ2 tests were performed to determine if significant differences exist in the genotype and allele frequencies between the cases and controls. χ2 analysis returned no statistically significant difference (p > 0.05) for genotype frequencies between cases and controls for 1758G > C (χ2 = 0.97, p = 0.6158) or 2880A > G (χ2 = 2.09, p = 0.3516). Similarly, no statistical difference was observed for allele frequencies for 1758G > C (χ2 = 0.07, p = 0.7867) or 2880A > G (χ2 = 0.04, p = 0.8365). Haplotype analysis of the two SNPs also showed no difference between the cases and the controls (p = 0.9585). Our findings in an Australian Caucasian population composed of breast cancer sufferers and an age matched control population did not support the findings of previous studies demonstrating that these markers play a significant role in breast cancer susceptibility. Here, no significant difference was detected between breast cancer patients and healthy matched controls by either the genotype or allele frequencies for the investigated variants (all p ≥ 0.05). While an association of the two variants and breast cancer was not detected in our case–control study population, exploring these variants in a larger population of the same kind may obtain results in concordance with previous studies. Given the importance of NCOA3 and its involvement in biological processes involved in breast cancer and the possible implications variants of the gene could have on the response to Tamoxifen therapy, NCOA3 remains a candidate for further investigations.

Investigation of two Wnt signalling pathway single nucleotide polymorphisms in a breast cancer-affected Australian population

In the mammary gland, Wnt signals are strongly implicated in initial development of the mammary rudiments and in the ductal branching and alveolar morphogenesis that occurs during pregnancy. Previously, we identified two Wnt signaling pathway-implicated genes, PPP3CA and MARK4, as having a role in more aggressive and potentially metastatic breast tumors. In this study, we examined two SNPs within PPP3CA and MARK4 in an Australian case-control study population for a potential role in human breast cancers. 182 cases and 180 controls were successfully genotyped for the PPP3CA SNP (rs2850328) and 182 cases and 177 controls were successfully genotyped for the MARK4 SNP (rs2395) using High Resolution Melt (HRM) analysis. Genotypes of randomly selected samples for both SNPs were validated by dye terminator sequencing. Chi-square tests were performed to determine any significant differences in the genotype and allele frequencies between the cases and controls. Chi-square analysis showed no statistically significant difference (p > .05) for genotype frequencies between cases and controls for rs2850328 (χ2 = 1.2, p = .5476) or rs2395 (χ2 = .3, p = .8608). Similarly, no statistical difference was observed for allele frequencies for rs2850328 (χ2 = .68, p = .4108) or rs2395 (χ2 = .02, p = .893). Even though an association of the polymorphisms rs2850328 and rs2395 and breast cancer was not detected in our case-control study population, other variants within the PPP3CA and MARK4 genes may still be associated with breast cancer, as both genes are implicated with processes involved in the disease as well as their mutual partaking in the Wnt signaling pathway.

Linkage mapping of CVD risk traits in the isolated Norfolk Island population

To understand the underlying genetic architecture of cardiovascular disease (CVD) risk traits, we undertook a genome-wide linkage scan to identify CVD quantitative trait loci (QTLs) in 377 individuals from the Norfolk Island population. The central aim of this research focused on the utilization of a genetically and geographically isolated population of individuals from Norfolk Island for the purposes of variance component linkage analysis to identify QTLs involved in CVD risk traits. Substantial evidence supports the involvement of traits such as systolic and diastolic blood pressures, high-density lipoprotein-cholesterol, low-density lipoprotein-cholesterol, body mass index and triglycerides as important risk factors for CVD pathogenesis. In addition to the environmental inXuences of poor diet, reduced physical activity, increasing age, cigarette smoking and alcohol consumption, many studies have illustrated a strong involvement of genetic components in the CVD phenotype through family and twin studies. We undertook a genome scan using 400 markers spaced approximately 10 cM in 600 individuals from Norfolk Island. Genotype data was analyzed using the variance components methods of SOLAR. Our results gave a peak LOD score of 2.01 localizing to chromosome 1p36 for systolic blood pressure and replicated previously implicated loci for other CVD relevant QTLs.

A population genomics overview of the neuronal nitric oxide synthase (nNOS) gene and its relationship to migraine susceptibility

The ubiquitous chemical messenger molecule nitric oxide (NO) has been implicated in a diverse range of biological activities including neurotransmission, smooth muscle motility and mediation of nociception. Endogenous synthesis of NO by the neuronal isoform of the nitric oxide synthase gene family has an essential role within the central and peripheral nervous systems in addition to the autonomic innervation of cerebral blood vessels. To investigate the potential role of NO and more specifically the neuronal nitric oxide synthase (nNOS) gene in migraine susceptibility, we investigated two microsatellite repeat variants residing within the 5′ and 3′ regions of the nNOS gene. Population genomic evaluation of the two nNOS repeat variants indicated significant linkage disequilibrium between the two loci. Z-DNA conformational sequence structures within the 5′ region of the nNOS gene have the potential to enhance or repress gene promoter activity. We suggest that genetic analysis of this 5′ repeat variant is the more functional variant expressing gene wide information that could affect endogenous NO synthesis and potentially result in diseased states. However, no association with migraine (with or without aura) was seen in our extensive case-control cohort (n = 579 affected with matched controls), when both the 5′ and 3′ genetic variants were investigated.

An R2R3 MYB transcription factor determines red petal colour in an Actinidia (kiwifruit) hybrid population

Background Red colour in kiwifruit results from the presence of anthocyanin pigments. Their expression, however, is complex, and varies among genotypes, species, tissues and environments. An understanding of the biosynthesis, physiology and genetics of the anthocyanins involved, and the control of their expression in different tissues, is required. A complex, the MBW complex, consisting of R2R3-MYB and bHLH transcription factors together with a WD-repeat protein, activates anthocyanin 3-O-galactosyltransferase (F3GT1) to produce anthocyanins. We examined the expression and genetic control of anthocyanins in flowers of Actinidia hybrid families segregating for red and white petal colour. Results Four inter-related backcross families between Actinidia chinensis Planch. var. chinensis and Actinidia eriantha Benth. were identified that segregated 1:1 for red or white petal colour. Flower pigments consisted of five known anthocyanins (two delphinidin-based and three cyanidin-based) and three unknowns. Intensity and hue differed in red petals from pale pink to deep magenta, and while intensity of colour increased with total concentration of anthocyanin, no association was found between any particular anthocyanin data and hue. Real time qPCR demonstrated that an R2R3 MYB, MYB110a, was expressed at significant levels in red-petalled progeny, but not in individuals with white petals. A microsatellite marker was developed that identified alleles that segregated with red petal colour, but not with ovary, stamen filament, or fruit flesh colour in these families. The marker mapped to chromosome 10 in Actinidia. The white petal phenotype was complemented by syringing Agrobacterium tumefaciens carrying Actinidia 35S::MYB110a into the petal tissue. Red pigments developed in white petals both with, and without, co-transformation with Actinidia bHLH partners. MYB110a was shown to directly activate Actinidia F3GT1 in transient assays. Conclusions The transcription factor, MYB110a, regulates anthocyanin production in petals in this hybrid population, but not in other flower tissues or mature fruit. The identification of delphinidin-based anthocyanins in these flowers provides candidates for colour enhancement in novel fruits.

A population-specific HTR2B stop codon predisposes to severe impulsivity

Impulsivity, describing action without foresight, is an important feature of several psychiatric diseases, suicidality and violent behaviour. The complex origins of impulsivity hinder identification of the genes influencing it and the diseases with which it is associated. Here we perform exon-focused sequencing of impulsive individuals in a founder population, targeting fourteen genes belonging to the serotonin and dopamine domain. A stop codon in HTR2B was identified that is common (minor allele frequency > 1%) but exclusive to Finnish people. Expression of the gene in the human brain was assessed, as well as the molecular functionality of the stop codon, which was associated with psychiatric diseases marked by impulsivity in both population and family-based analyses. Knockout of Htr2b increased impulsive behaviours in mice, indicative of predictive validity. Our study shows the potential for identifying and tracing effects of rare alleles in complex behavioural phenotypes using founder populations, and indicates a role for HTR2B in impulsivity.

Glutathione transferase Ti and M1 genotype polymorphism in the normal population of Shanghai

Glutathione transferases are known to be important enzymes in the metabolism of xenobiotics. In humans genetic polymorphisms have been reported for the hGSTM1 and hGSTT1 genes leading to individual differences in susceptibility towards toxic effects, such as cancer. This study describes the distribution of the two polymorphisms of hGSTT1 and hGSTM1 in the normal Chinese population of Shanghai. Out of 219 healthy individuals having been genotyped for GSTTI and GSTMI, 108 (49%) were identified to be homozygously deficient for the GSTT1 gene and 107 (49%) for the GSTM1 gene.

Conservation genetics of the water mouse, Xeromys myoides Thomas, 1889

The water mouse, Xeromys myoides, is currently recognised as a vulnerable species in Australia, inhabiting a small number of distinct and isolated coastal regions of Queensland and the Northern Territory. An examination of the evolutionary history and contemporary influences shaping the genetic structure of this species is required to make informed conservation management decisions. Here, we report the first analysis undertaken on the phylogeography and population genetics of the water mouse across its mainland Australian distribution. Genetic diversity was assessed at two mitochondrial DNA (Cytochrome b, 1000 bp; D-loop, 400 bp) and eight microsatellite DNA loci. Very low genetic diversity was found, indicating that water mice underwent a recent expansion throughout their Australian range and constitute a single evolutionarily significant unit. Microsatellite analyses revealed that the highest genetic diversity was found in the Mackay region of central Queensland; population substructure was also identified, suggesting that local populations may be isolated in this region. Conversely, genetic diversity in the Coomera region of south-east Queensland was very low and the population in this region has experienced a significant genetic bottleneck. These results have significant implications for future management, particularly in terms of augmenting populations through translocations or reintroducing water mice in areas where they have gone extinct.

Association of heparan sulfate proteoglycans SDC1 and SDC4 polymorphisms with breast cancer in an Australian Caucasian population

Breast cancer is a common disease in both developing and developed countries with early identification and treatment improving prognosis and survival. Heparan sulfate proteoglycans (HSPGs) are key components of the extracellular matrix (ECM) that mediate cell adhesion, motility, proliferation, invasion and cell signalling. Members of the syndecan family of HSPGs have been identified to be involved in breast cancer progression through their varied interactions with a number of growth factors, ligands and receptors. Specifically, high expression levels of syndecan-1 (SDC1) have been demonstrated in more invasive breast tumours while elevated syndecan-4 (SDC4) levels have been identified to correspond with improved prognosis. With genetic changes in the syndecans and their association with breast cancers plausible, we examined two single nucleotide polymorphisms in SDC1 (rs1131351) and SDC4 (rs67068737) within an Australian Caucasian breast cancer case/control population. No association was found with SDC4 and breast cancer in our population. However, a significant association between SDC1 and breast cancer was identified in both our case/control population and in a replication cohort. When both populations were combined for analysis, this association became more significant (genotype, p = 0.0003; allele, p = 0.0001). This data suggests an increased risk of developing breast cancer associated with the presence of the C allele of the SDC1 rs1131351 single nucleotide polymorphism (SNP) and may provide a marker toward early breast cancer detection.

Brain fiber architecture, genetics, and intelligence: a high angular resolution diffusion imaging (HARDI) study

We developed an analysis pipeline enabling population studies of HARDI data, and applied it to map genetic influences on fiber architecture in 90 twin subjects. We applied tensor-driven 3D fluid registration to HARDI, resampling the spherical fiber orientation distribution functions (ODFs) in appropriate Riemannian manifolds, after ODF regularization and sharpening. Fitting structural equation models (SEM) from quantitative genetics, we evaluated genetic influences on the Jensen-Shannon divergence (JSD), a novel measure of fiber spatial coherence, and on the generalized fiber anisotropy (GFA) a measure of fiber integrity. With random-effects regression, we mapped regions where diffusion profiles were highly correlated with subjects' intelligence quotient (IQ). Fiber complexity was predominantly under genetic control, and higher in more highly anisotropic regions; the proportion of genetic versus environmental control varied spatially. Our methods show promise for discovering genes affecting fiber connectivity in the brain.