A model of macadamia with application to pruning in orchards

A self-organising model of macadamia, expressed using L-Systems, was used to explore aspects of canopy management. A small set of parameters control the basic architecture of the model, with a high degree of self-organisation occurring to determine the fate and growth of buds. Light was sensed at the leaf level and used to represent vigour and accumulated basipetally. Buds also sensed light so as to provide demand in the subsequent redistribution of the vigour. Empirical relationships were derived from a set of 24 completely digitised trees after conversion to multiscale tree graphs (MTG) and analysis with the OpenAlea software library. The ability to write MTG files was embedded within the model so that various tree statistics could be exported for each run of the model. To explore the parameter space a series of runs was completed using a high-throughput computing platform. When combined with MTG generation and analysis with OpenAlea it provided a convenient way in which thousands of simulations could be explored. We allowed the model trees to develop using self-organisation and simulated cultural practices such as hedging, topping, removal of the leader and limb removal within a small representation of an orchard. The model provides insight into the impact of these practices on potential for growth and the light distribution within the canopy and to the orchard floor by coupling the model with a path-tracing program to simulate the light environment. The lessons learnt from this will be applied to other evergreen, tropical fruit and nut trees.

Latent variable models in statistical genetics

Understanding the complexities that are involved in the genetics of multifactorial diseases is still a monumental task. In addition to environmental factors that can influence the risk of disease, there is also a number of other complicating factors. Genetic variants associated with age of disease onset may be different from those variants associated with overall risk of disease, and variants may be located in positions that are not consistent with the traditional protein coding genetic paradigm. Latent Variable Models are well suited for the analysis of genetic data. A latent variable is one that we do not directly observe, but which is believed to exist or is included for computational or analytic convenience in a model. This thesis presents a mixture of methodological developments utilising latent variables, and results from case studies in genetic epidemiology and comparative genomics. Epidemiological studies have identified a number of environmental risk factors for appendicitis, but the disease aetiology of this oft thought useless vestige remains largely a mystery. The effects of smoking on other gastrointestinal disorders are well documented, and in light of this, the thesis investigates the association between smoking and appendicitis through the use of latent variables. By utilising data from a large Australian twin study questionnaire as both cohort and case-control, evidence is found for the association between tobacco smoking and appendicitis. Twin and family studies have also found evidence for the role of heredity in the risk of appendicitis. Results from previous studies are extended here to estimate the heritability of age-at-onset and account for the eect of smoking. This thesis presents a novel approach for performing a genome-wide variance components linkage analysis on transformed residuals from a Cox regression. This method finds evidence for a dierent subset of genes responsible for variation in age at onset than those associated with overall risk of appendicitis. Motivated by increasing evidence of functional activity in regions of the genome once thought of as evolutionary graveyards, this thesis develops a generalisation to the Bayesian multiple changepoint model on aligned DNA sequences for more than two species. This sensitive technique is applied to evaluating the distributions of evolutionary rates, with the finding that they are much more complex than previously apparent. We show strong evidence for at least 9 well-resolved evolutionary rate classes in an alignment of four Drosophila species and at least 7 classes in an alignment of four mammals, including human. A pattern of enrichment and depletion of genic regions in the profiled segments suggests they are functionally significant, and most likely consist of various functional classes. Furthermore, a method of incorporating alignment characteristics representative of function such as GC content and type of mutation into the segmentation model is developed within this thesis. Evidence of fine-structured segmental variation is presented.

Predicting loss and fragmentation of habitat of the vulnerable subtropical rainforest tree Macadamia integrifolia with models developed from compiled ecological data

Habitat models are widely used in ecology, however there are relatively few studies of rare species, primarily because of a paucity of survey records and lack of robust means of assessing accuracy of modelled spatial predictions. We investigated the potential of compiled ecological data in developing habitat models for Macadamia integrifolia, a vulnerable mid-stratum tree endemic to lowland subtropical rainforests of southeast Queensland, Australia. We compared performance of two binomial models—Classification and Regression Trees (CART) and Generalised Additive Models (GAM)—with Maximum Entropy (MAXENT) models developed from (i) presence records and available absence data and (ii) developed using presence records and background data. The GAM model was the best performer across the range of evaluation measures employed, however all models were assessed as potentially useful for informing in situ conservation of M. integrifolia, A significant loss in the amount of M. integrifolia habitat has occurred (p < 0.05), with only 37% of former habitat (pre-clearing) remaining in 2003. Remnant patches are significantly smaller, have larger edge-to-area ratios and are more isolated from each other compared to pre-clearing configurations (p < 0.05). Whilst the network of suitable habitat patches is still largely intact, there are numerous smaller patches that are more isolated in the contemporary landscape compared with their connectedness before clearing. These results suggest that in situ conservation of M. integrifolia may be best achieved through a landscape approach that considers the relative contribution of small remnant habitat fragments to the species as a whole, as facilitating connectivity among the entire network of habitat patches.

Bayesian methodology for genetics of complex diseases

Genetic research of complex diseases is a challenging, but exciting, area of research. The early development of the research was limited, however, until the completion of the Human Genome and HapMap projects, along with the reduction in the cost of genotyping, which paves the way for understanding the genetic composition of complex diseases. In this thesis, we focus on the statistical methods for two aspects of genetic research: phenotype definition for diseases with complex etiology and methods for identifying potentially associated Single Nucleotide Polymorphisms (SNPs) and SNP-SNP interactions. With regard to phenotype definition for diseases with complex etiology, we firstly investigated the effects of different statistical phenotyping approaches on the subsequent analysis. In light of the findings, and the difficulties in validating the estimated phenotype, we proposed two different methods for reconciling phenotypes of different models using Bayesian model averaging as a coherent mechanism for accounting for model uncertainty. In the second part of the thesis, the focus is turned to the methods for identifying associated SNPs and SNP interactions. We review the use of Bayesian logistic regression with variable selection for SNP identification and extended the model for detecting the interaction effects for population based case-control studies. In this part of study, we also develop a machine learning algorithm to cope with the large scale data analysis, namely modified Logic Regression with Genetic Program (MLR-GEP), which is then compared with the Bayesian model, Random Forests and other variants of logic regression.

The genetics of cutaneous melanoma

Although germline mutations in CDKN2A are present in approximately 25% of large multicase melanoma families, germline mutations are much rarer in the smaller melanoma families that make up most individuals reporting a family history of this disease. In addition, only three families worldwide have been reported with germline mutations in a gene other than CDKN2A (i.e., CDK4). Accordingly, current genomewide scans underway at the National Human Genome Research Institute hope to reveal linkage to one or more chromosomal regions, and ultimately lead to the identification of novel genes involved in melanoma predisposition. Both CDKN2A and PTEN have been identified as genes involved in sporadic melanoma development; however, mutations are more common in cell lines than uncultured tumors. A combination of cytogenetic, molecular, and functional studies suggests that additional genes involved in melanoma development are located to chromosomal regions 1p, 6q, 7p, 11q, and possibly also 9p and 10q. With the near completion of the human genome sequencing effort, combined with the advent of high throughput mutation analyses and new techniques including cDNA and tissue microarrays, the identification and characterization of additional genes involved in melanoma pathogenesis seem likely in the near future.

Quality of life and genetics in men with prostate cancer

As family history has been established as a risk factor for prostate cancer, attempts have been made to isolate predisposing genetic variants that are related to hereditary prostate cancer. With many genetic variants still to be identified and investigated, it is not yet possible to fully understand the impact of genetic variants on prostate cancer development. The high survival rates among men with prostate cancer have meant that other issues, such as quality of life (QoL), have also become important. Through their effect on a person’s health, a range of inherited genetic variants may potentially influence QoL in men with prostate cancer, even prior to treatment. Until now, limited research has been conducted on the relationship between genetics and QoL. Thus, this study contributes to an emerging field by aiming to identify certain genetic variants related to the QoL found in men with prostate cancer. It is hoped that this study may lead to future research that will identify men who have an increased risk of a poor QoL following prostate cancer treatment, which will aid in developing treatments that are individually tailored to support them. Previous studies have established that genetic variants of Vascular Endothelial Growth Factor (VEGF) and Insulin-like Growth Factor 1 (IGF-1) may play a role in prostate cancer development. VEGF and IGF-1 have also been reported to be associated with QoL in people with ovarian cancer and colorectal cancer, respectively. This study completed a series of secondary analyses using two major data-sets (from 850 men newly diagnosed with prostate cancer, and approximately 550 men from the general Queensland population), in which genetic variants of VEGF and IGF-1 were investigated for associations with prostate cancer susceptibility and QoL. The first aim of this research was to investigate genetic variants in the VEGF and IGF-I gene for an association with the risk of prostate cancer. It was found that one IGF-1 genetic variant (rs35765) had a statistically significant association with prostate cancer (p = 0.04), and one VEGF genetic variant (rs2146323) had a statistically significant association with advanced prostate cancer (p = 0.02). The estimates suggest that carriers of the CA and AA genotype for rs35765 may have a reduced risk of developing prostate cancer (Odds Ratio (OR) = 0.72, 95% Confidence Interval (CI) = 0.55, 0.95, OR = 0.60, 95% CI = 0.26, 1.39, respectively). Meanwhile, carriers of the CA and AA genotype for rs2146323 may be at increased risk of advanced prostate cancer, which was determined by a Gleason score of above 7 (OR = 1.72, 95% CI = 1.12, 2.63, OR = 1.90, 95% CI = 1.08, 3.34, respectively). Utilising the widely used short-form health survey, the SF-36v2, the second aim of this study was to investigate the relationship between prostate cancer and QoL prior to treatment. Assessing QoL at this time-point was important as little research has been conducted to evaluate if prostate cancer affects QoL regardless of treatment. The analyses found that mean SF-36v2 scale scores related to physical health were higher by at least 0.3 Standard Deviations (SD) among men with prostate cancer than the general population comparison group. This difference was considered clinically significant (defined by group differences in mean SF-36v2 scores by at least 0.3 SD). These differences were also statistically significant (p<0.05). Mean QoL scale scores related to mental health were similar between men with prostate cancer and those from the general population comparison group. The third aim of this study was to investigate genetic variants in the VEGF and IGF-1 gene for an association with QoL in prostate cancer patients prior to their treatment. It was essential to evaluate these relationships prior to treatment, before the involvement of these genes was potentially interrupted by treatment. The analyses found that some genetic variants had a small clinically significant association (0.3 SD) to some QoL domains experienced by these men. However, most relationships were not statistically significant (p>0.05). Most of the associations found identified that a small sub-group of men with prostate cancer (approximately 2%) reported, on average, a slightly better QoL than the majority of the prostate cancer patients. The fourth aim of this research was to investigate whether associations between genetic variants in VEGF and IGF-1 and QoL were specific to men with prostate cancer, or were also applicable to the general male population. It was found that twenty out of one-hundred relationships between the genetic variants of VEGF and IGF-1 and QoL health-measures and scales examined differed between these groups. In the majority of the relationships involving VEGF SNPs that differed, a clinically significant difference (0.3 or more SD) between mean scores among the genotype groups in prostate cancer patients was found, while mean scores among men from the general-population comparison group were similar. For example, prostate cancer participants who carried at least one T allele (CT or TT genotype) for rs3024994 had a clinically significant higher (0.3 SD) mean QoL score in terms of the role-physical scale, than participants who carried the CC genotype. This was not seen among men from the general population sample, as the mean score was similar between genotype groups. The opposite was seen in regards to the IGF-1 SNPs examined. Overall, these relationships were not considered to directly impact on the clinical options for men with prostate cancer. As this study utilised secondary data from two separate studies, there are a number of important limitations that should be acknowledged including issues of multiple comparisons, power, and missing or unavailable data. It is recommended that this study be replicated as a better-designed study that takes greater consideration of the many factors involved in prostate cancer and QoL. Investigation into other genetic variants of VEGF or IGF-1 is also warranted, as is consideration of other genes and their relationship with QoL. Through identifying certain genetic variants that have a modest association to prostate cancer, this project adds to the knowledge surrounding VEGF and IGF-1 and their role in prostate cancer susceptibility. Importantly, this project has also introduced the potential role genetics plays in QoL, through investigating the relationships between genetic variants of VEGF and IGF-1 and QoL.

Effects of cultivar on rodent damage in Australian macadamia orchards

The Black rat (Rattus rattus), a serious pest of Australian macadamia orchards has been estimated to cause up to 30% crop damage in Australian orchards. In recent years an increase in the number of commercially available cultivars has seen a change in orchard characteristics in Australia, primarily effecting fruiting and flowering patterns. This has been suggested to affect the feeding behaviour of rodents and in turn altered the damage process. In this study we compare the extent of damage in orchards containing one of three prevalent cultivars (A4/A16, A268 and HAES 344/741) and investigate the influence of these cultivars, particularly their distinctive fruiting traits, on rodent damage within the orchard. We demonstrate that the temporal pattern and extent of damage differs between cultivar types. Newer Australian macadamia cultivars tested in this study were found to be far more susceptible to rodent damage than the older Hawaiian developed cultivars, most likely due to an extended fruiting period and thinner shells. This has resulted in a more sustained period of crop damage than the patterns of crop damage observed in previous Australian studies. Crop damage caused by R. rattus is significantly higher in orchards that maintain high levels of canopy resources through the fruiting season and we postulate that this is due to the extended fruiting periods of the new cultivars used. The maintenance of canopy resource load in turn corresponds to high crop damage, in this study resulting in crop losses of up to 25%.

The spatial distribution of rodent damage in Australian macadamia (Macadamia integrifolia) orchards

The Black Rat (Rattus rattus), a global pest within the macadamia production industry, causes up to 30% crop damage in Australian orchards. During early stages of production in Australia, research demonstrated the importance of non crop adjacent habitats as significant in affecting the patterns of crop damage seen throughout orchards. Where once rodent damage was limited to the outside edges of orchard blocks, growers are now reporting finding crop damage throughout entire orchards. This study therefore aims to explore the spatial patterns of rodent distribution and damage now occurring in Australian macadamia orchards. We show that rodent damage and rodent distribution in these newer production regions differ from that shown in previous Australian research. Previous Australian research has shown damage patterns which were associated with the edges of orchard blocks however this study demonstrates a more widespread damage distribution. In the current study there is no relationship between rodent damage and the orchard edge. Arboreal rodent nests were identified within these newer orchard systems, suggesting rodents are residing within the tree component of the orchard system and not dependent on adjacent non-crop habitat for shelter. Results from this study confirm that rodents have modified their nesting and foraging behaviour in newer orchards systems in Australia. We suggest that this is a response of increased and prolonged availability of macadamia nuts in newer production regions enabling populations to be maintained throughout the year. Management strategies will require modification if control is to be achieved.

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.

Molecular genetics of migraine

Migraine is a common neurological disorder with a significantly heritable component. It is a complex disease and despite numerous molecular genetic studies, the exact pathogenesis causing the neurological disturbance remains poorly understood. Although several known molecular mechanisms have been associated with an increased risk for developing migraine, there remains significant scope for future studies. The majority of studies have investigated the most plausible candidate genes involved in common migraine pathogenesis utilising criteria that takes into account a combination of physiological functionality in conjunction with regions of genomic association. Thus, far genes involved in neurological, vascular or hormonal pathways have been identified and investigated on this basis. Genome-wide association studies (GWAS) studies have helped to identify novel regions that may be associated with migraine and have aided in providing the basis for further molecular investigations. However, further studies utilising sequencing technologies are required to characterise the genetic basis for migraine.