Genetic and Functional Evidence Implicating DLL1 as the Gene That Influences Susceptibility to Visceral Leishmaniasis at Chromosome 6q27

Background. Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. Methods. Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. Results. Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 x 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 x 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 x 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. Conclusions. DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

Fine mapping and replication of QTL in outbred chicken advanced intercross lines

Background: Linkage mapping is used to identify genomic regions affecting the expression of complex traits. However, when experimental crosses such as F2 populations or backcrosses are used to map regions containing a Quantitative Trait Locus (QTL), the size of the regions identified remains quite large, i.e. 10 or more Mb. Thus, other experimental strategies are needed to refine the QTL locations. Advanced Intercross Lines (AIL) are produced by repeated intercrossing of F2 animals and successive generations, which decrease linkage disequilibrium in a controlled manner. Although this approach is seen as promising, both to replicate QTL analyses and fine-map QTL, only a few AIL datasets, all originating from inbred founders, have been reported in the literature. Methods: We have produced a nine-generation AIL pedigree (n = 1529) from two outbred chicken lines divergently selected for body weight at eight weeks of age. All animals were weighed at eight weeks of age and genotyped for SNP located in nine genomic regions where significant or suggestive QTL had previously been detected in the F2 population. In parallel, we have developed a novel strategy to analyse the data that uses both genotype and pedigree information of all AIL individuals to replicate the detection of and fine-map QTL affecting juvenile body weight. Results: Five of the nine QTL detected with the original F2 population were confirmed and fine-mapped with the AIL, while for the remaining four, only suggestive evidence of their existence was obtained. All original QTL were confirmed as a single locus, except for one, which split into two linked QTL. Conclusions: Our results indicate that many of the QTL, which are genome-wide significant or suggestive in the analyses of large intercross populations, are true effects that can be replicated and fine-mapped using AIL. Key factors for success are the use of large populations and powerful statistical tools. Moreover, we believe that the statistical methods we have developed to efficiently study outbred AIL populations will increase the number of organisms for which in-depth complex traits can be analyzed.

Novel Statistical Methods in Quantitative Genetics : Modeling Genetic Variance for Quantitative Trait Loci Mapping and Genomic Evaluation

This thesis develops and evaluates statistical methods for different types of genetic analyses, including quantitative trait loci (QTL) analysis, genome-wide association study (GWAS), and genomic evaluation. The main contribution of the thesis is to provide novel insights in modeling genetic variance, especially via random effects models. In variance component QTL analysis, a full likelihood model accounting for uncertainty in the identity-by-descent (IBD) matrix was developed. It was found to be able to correctly adjust the bias in genetic variance component estimation and gain power in QTL mapping in terms of precision.  Double hierarchical generalized linear models, and a non-iterative simplified version, were implemented and applied to fit data of an entire genome. These whole genome models were shown to have good performance in both QTL mapping and genomic prediction. A re-analysis of a publicly available GWAS data set identified significant loci in Arabidopsis that control phenotypic variance instead of mean, which validated the idea of variance-controlling genes.  The works in the thesis are accompanied by R packages available online, including a general statistical tool for fitting random effects models (hglm), an efficient generalized ridge regression for high-dimensional data (bigRR), a double-layer mixed model for genomic data analysis (iQTL), a stochastic IBD matrix calculator (MCIBD), a computational interface for QTL mapping (qtl.outbred), and a GWAS analysis tool for mapping variance-controlling loci (vGWAS).

Modelling dominance in a flexible intercross analysis

The aim of this paper is to develop a flexible model for analysis of quantitative trait loci (QTL) in outbred line crosses, which includes both additive and dominance effects. Our flexible intercross analysis (FIA) model accounts for QTL that are not fixed within founder lines and is based on the variance component framework. Genome scans with FIA are performed using a score statistic, which does not require variance component estimation. RESULTS: Simulations of a pedigree with 800 F2 individuals showed that the power of FIA including both additive and dominance effects was almost 50% for a QTL with equal allele frequencies in both lines with complete dominance and a moderate effect, whereas the power of a traditional regression model was equal to the chosen significance value of 5%. The power of FIA without dominance effects included in the model was close to those obtained for FIA with dominance for all simulated cases except for QTL with overdominant effects. A genome-wide linkage analysis of experimental data from an F2 intercross between Red Jungle Fowl and White Leghorn was performed with both additive and dominance effects included in FIA. The score values for chicken body weight at 200 days of age were similar to those obtained in FIA analysis without dominance. CONCLUSION: We have extended FIA to include QTL dominance effects. The power of FIA was superior, or similar, to standard regression methods for QTL effects with dominance. The difference in power for FIA with or without dominance is expected to be small as long as the QTL effects are not overdominant. We suggest that FIA with only additive effects should be the standard model to be used, especially since it is more computationally efficient.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in upto 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 x 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for similar to 2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous systemin obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

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 ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼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/β-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.2014

Suppression by ABA of salicylic acid and lignin accumulation and the expression of multiple genes, in Arabidopsis infected with Pseudomonas syringae pv. tomato

Abscisic acid (ABA) has been implicated in determining the outcome of interactions between many plants and their pathogens. We had previously shown that increased concentrations of ABA within leaves of Arabidopsis induced susceptibility towards an avirulent strain of Pseudomonas syringae pathovar (pv.) tomato. We now show that ABA induces susceptibility via suppression of the accumulation of components crucial for a resistance response. Lignin and salicylic acid concentrations in leaves were increased during a resistant interaction but reduced when plants were treated with ABA. The reduction in lignin and salicylic acid production was independent of the development of the hypersensitive response (HR), indicating that, in this host-pathogen system, HR is not required for resistance. Genome-wide gene expression analysis using microarrays showed that treatment with ABA suppressed the expression of many defence-related genes, including those important for phenylpropanoid biosynthesis and those encoding resistance-related proteins. Together, these results show that resistance induction in Arabidopsis to an avirulent strain of P. syringae pv. tomato is regulated by ABA.

Disease fingerprinting with cDNA microarrays reveals distinct gene expression profiles in lethal type 1 and type 2 cytokine-mediated inflammatory reactions

Development of polarized immune responses controls resistance and susceptibility to many microorganisms. However, studies of several infectious, allergic, and autoimmune diseases have shown that chronic type-1 and type-2 cytokine responses can also cause significant morbidity and mortality if left unchecked. We used mouse cDNA microarrays to molecularly phenotype the gene expression patterns that characterize two disparate but equally lethal forms of liver pathology that develop in Schistosoma mansoni infected mice polarized for type-1 and type-2 cytokine responses. Hierarchical clustering analysis identified at least three groups of genes associated with a polarized type-2 response and two linked with an extreme type-1 cytokine phenotype. Predictions about liver fibrosis,  apoptosis, and granulocyte recruitment and activation generated by the microarray studies were confirmed later by traditional biological assays. The data show that cDNA microarrays are useful not only for determining  coordinated gene expression profiles but are also highly effective for molecularly “fingerprinting” diseased tissues. Moreover, they illustrate the potential of genome-wide approaches for generating comprehensive views on the molecular and biochemical mechanisms regulating infectious  disease pathogenesis.

Investigation of plasmodiophora brassicae (clubroot disease) in vegetable brassica using arabidopsis thaliana as a model system

Clubroot, caused by Plasmodiophora brassicae, is the most devastating soil-borne disease of vegetable brassicas. It occurs all over the world and is responsible for crop losses of up to 10% every year. In Australia, the disease is being managed effectively with chemicals and cultural practices, but ideally control can be improved in the long term by the introduction of resistant cultivars. The life cycle ofP. brassicae and mode of action of plant resistance has not been fully elucidated because of the technical difficulties of working with an obligate, soil-borne plant pathogen. However, Arabidopsis thaliana, which is a host ofP. brassicae, has great potential as a model system for studying the life cycle, the infection process and development of resistance. We have developed a sand-liquid-culture system for growing Arabidopsis that allows easy observation of all life stages and, most importantly, the primary plasmodial stages within the root hair. The method was first optimised for observations of the lifecycle of the pathogen in a susceptible Arabidopsis ecotype (Col-3) where all stages of the lifecycle have now been observed and characterised. Further screening of Arabidopsis ecotypes for disease resistance has utilised one of the most virulent Australian pathotypes of brassica (ECD number 16/19/31). To date, Arabidopsis ecotype Ta-0 has shown a level of tolerance to the disease even though the roots get infected. It has been reported earlier that resistance toP. brassicae in Arabidopsis is due to one or a small number of genes. To examine changes in gene expression during the early, critical stages of infection, RNA was extracted from the susceptible and resistant ecotypes at two time points, 4 days and 17 days after inoculation. Microarray analysis will be used to investigate genome wide changes in gene expression during infection but also to identify candidate genes that may confer resistance to Australian isolates of the pathogen.

Abscisic acid regulation of plant defence responses during pathogen attack

The plant hormone, abscisic acid (ABA), has previously been shown to have an impact on the resistance or susceptibility of plants to pathogens. In this thesis, it was shown that ABA had a regulatory effect on an extensive array of plant defence responses in three different plant and pathogen interaction combinations as well as following the application of an abiotic elicitor. In unique studies using ABA deficient mutants of Arabidopsis, exogenous ABA addition or ABA biosynthesis inhibitor application and simulated drought stress, ABA was shown to have a profound effect on the outcome of interactions between plants and pathogens of differing lifestyles and from different kingdoms. The systems used included a model plant and an important agricultural species: Arabidopsis thaliana (Arabidopsis) and Peronospora parasitica (a biotrophic Oomycete pathogen), Arabidopsis and Pseudomonas syringae pathovar tomato (a biotrophic bacterial pathogen) and an unrelated plant species, soybean (Glycine max) and Phytophthora sojae (a hemibiotrophic Oomycete pathogen), Generally, a higher than basal endogenous ABA concentration within plant tissues at the time of avirulent pathogen inoculation, caused an interaction shift towards what phenotypically resembled susceptibility. Conversely, a lower than basal endogenous ABA concentration in plants inoculated with a virulent pathogen caused a shift towards resistance. An extensive suppressive effect of ABA on defence responses was revealed by a range of techniques that included histochemical, biochemical and molecular approaches. A universal effect of ABA on suppression or induction of the phenylpropanoid pathway via regulation of the key entry point gene, phenylalanine ammonia-lyase (PAL), when stimulated by biotic or abiotic elicitors was shown. ABA also influenced a wide variety of other defence-related components such as: the development of a hypersensitive response (HR), the accumulation of the reactive oxyden species, hydrogen peroxide and the cell wall strengthening compounds lignin and callose, accumulation of SA and the phytoalexin, glyceollin and the transcription of the SA-dependent pathogenesis- related gene (PR-1). The near genome-wide microarray gene expression analysis of an ABA induced susceptible interaction also revealed an yet unprecedented insight into the great diversity of defence responses that were influenced by ABA that included: disease resistance like proteins, antimicrobial proteins as well as phenylpropanoid and tryptophan pathway enzymes. Subtle differences were found in the number and type of defence responses that were regulated by ABA in each type of plant and pathogen interaction that was studied. This thesis has clearly identified in plant/pathogen interactions previously unknown and important roles for ABA in the regulation of many defence responses.

Familial and genomic analyses of postural changes in systolic and diastolic blood pressure

The physiological adaptation to the erect posture involves integrated neural and cardiovascular responses that might be determined by genetic factors. We examined the familial- and individual-specific components of variance for postural changes in systolic and diastolic blood pressure in 767 volunteer nuclear adult families from the Victorian Family Heart Study. In 274 adult sibling pairs, we made a genome-wide scan using 400 markers for quantitative trait loci linked with the postural changes in systolic and diastolic pressures. Overall, systolic pressure did not change on standing, but there was considerable variation in this phenotype (SD=8.1 mm Hg). Familial analyses revealed that 25% of the variance of change in systolic pressure was attributable to genetic factors. In contrast, diastolic pressure increased by 6.3 mm Hg (SD=7.0 mm Hg) on standing and there was no evidence of contributory genetic factors. Multipoint quantitative genome linkage mapping suggested evidence (Z=3.2) of linkage of the postural change in systolic pressure to chromosome 12 but found no genome-wide evidence of linkage for the change in diastolic pressure. These findings suggest that genetic factors determine whether systolic pressure decreases or increases when one stands, possibly as the result of unidentified alleles on chromosome 12. The genetics of postural changes in systolic blood pressure might reflect the general buffering function of the baroreflex; thereby, the predisposition to sudden decreases or increases in systolic pressure might cause postural hypotension or vessel wall disruption, respectively.

Defence mechanisms of a resistant monocot model to phytophthora cinnamomi

The plant pathogen Phytophthora cinnamomi causes devastating disease in natural and agricultural systems worldwide. While some plants can survive, little is known about the underlying mechanisms of resistance. This research used histochemical and genome-wide analysis to identify key cellular and molecular defence mechanisms within the resistant plant Zea mays.

Cordycepin interferes with 3' end formation in yeast independently of its potential to terminate RNA chain elongation

Cordycepin (3′ deoxyadenosine) is a biologically active compound that, when incorporated during RNA synthesis in vitro, provokes chain termination due to the absence of a 3′ hydroxyl moiety. We were interested in the effects mediated by this drug in vivo and analyzed its impact on RNA metabolism of yeast. Our results support the view that cordycepin-triphosphate (CoTP) is the toxic component that is limiting cell growth through inhibition of RNA synthesis. Unexpectedly, cordycepin treatment modulated 3′ end heterogeneity of ACT1 and ASC1 mRNAs and rapidly induced extended transcripts derived from CYH2 and NEL025c loci. Moreover, cordycepin ameliorated the growth defects of poly(A) polymerase mutants and the pap1-1 mutation neutralized the effects of the drug on gene expression. Our observations are consistent with an epistatic relationship between poly(A) polymerase function and cordycepin action and suggest that a major mode of cordycepin activity reduces 3′ end formation efficiency independently of its potential to terminate RNA chain elongation. Finally, chemical-genetic profiling revealed genome-wide pathways linked to cordycepin activity and identified novel genes involved in poly(A) homeostasis.

Toward forward genetic screens in malaria-causing parasites using the piggyBac transposon

The ability to analyze gene function in malaria-causing Plasmodium parasites has received a boost with a recent paper in BMC Genomics that describes a genome-wide mutagenesis system in the rodent malaria species Plasmodium berghei using the transposon piggyBac. This advance holds promise for identifying and validating new targets for intervention against malaria. But further improvements are still needed for the full power of genome-wide molecular genetic screens to be utilized in this organism.

Disease gene prediction database

This database includes gene predictions for disease phenotypes based on published Genome-Wide Association Data. May be used to choose primers for phenotype-specific resquencing of patient DNA.
For each prediction for following data is listed: phenotype, predicted gene, significant SNP, datasource, datasource reference.