Drafting the Kensington Pride Mango Genome

The mango industry in Australia is worth in excess of $150 million annually with the Kensington Pride (KP) cultivar capturing 60% of the domestic market. Valued by consumers for desirable taste and colour characteristics, KP has been used extensively as a parent in the Department of Agriculture and Fisheries’ (Queensland, Australia) mango breeding program with over 400 hybrid trees sharing KP as the male parent. In order to gain a better understanding of Australia’s most significant mango variety, Horticulture Innovation Australia had led an international collaboration between the Queensland Department of Agriculture and Fisheries (Australia), the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT, India) and the Beijing Genomics Institute (China) to sequence the KP genome. Preliminary de novo assembly of illumina short read sequence data suggests that the KP genome is highly heterozygous and has an estimated genome size of 407 Mb. As refinements and additional sequence data are added to the assembly, a more complete picture of the mango genome will be elucidated.

The Queensland Mango Genomics Initiative

Mango is an important industry for Queensland, Australia, with an annual value exceeding $80 million. The Kensington Pride cultivar, prized by consumers for desirable taste and colour characteristics, commands 60% of the domestic market though this market share has declined in recent years as new varieties, such as Calypso™, get established with consumers. In 2005, the Queensland Government's Department of Agriculture and Fisheries commenced the Mango Genomics Initiative. This project brought together multidisciplinary teams of breeders, pathologists, sensory scientists, flavour chemists and molecular biologists to develop a suite of tools and inter-related data sets to support the accelerated development of new commercial mango varieties. An overview of the Mango Genomics Initiative will be presented here culminating in the generation of a draft Kensington Pride mango genome sequence.

Temporal Variation in Physiological Biomarkers in Black Flying-Foxes (Pteropus alecto), Australia

Bats of the genus Pteropus (Pteropodidae) are recognised as the natural host of multiple emerging pathogenic viruses of animal and human health significance, including henipaviruses, lyssaviruses and ebolaviruses. Some studies have suggested that physiological and ecological factors may be associated with Hendra virus infection in flying-foxes in Australia; however, it is essential to understand the normal range and seasonal variability of physiological biomarkers before seeking physiological associations with infection status. We aimed to measure a suite of physiological biomarkers in P. alecto over time to identify any seasonal fluctuations and to examine possible associations with life-cycle and environmental stressors. We sampled 839 adult P. alecto in the Australian state of Queensland over a 12-month period. The adjusted population means of every assessed hematologic and biochemical parameter were within the reported reference range on every sampling occasion. However, within this range, we identified significant temporal variation in these parameters, in urinary parameters and body condition, which primarily reflected the normal annual life cycle. We found no evident effect of remarkable physiological demands or nutritional stress, and no indication of clinical disease driving any parameter values outside the normal species reference range. Our findings identify underlying temporal physiological changes at the population level that inform epidemiological studies and assessment of putative physiological risk factors driving Hendra virus infection in P. alecto. More broadly, the findings add to the knowledge of Pteropus populations in terms of their relative resistance and resilience to emerging infectious disease.

Genetics of T cell co-stimulatory receptors -CD28, CTLA4, ICOS and PDCD1 in immunity and transplantation

Co-stimulatory signals are essential for the activation of naïve T cells and productive immune response. Naïve T cells receive first, antigen-specific signal through T cell receptor. Co-stimulatory receptors provide the second signal which can be either activating or inhibitory. The balance between signals determines the outcome of an immune response. CD28 is crucial for T cell activation; whereas cytotoxic T lymphocyte associated antigen 4 (CTLA4) mediates critical inhibitory signal. Inducible co-stimulator (ICOS) augments cytokine expression and plays role in immunoglobulin class switching. Programmed cell death 1 (PDCD1) acts as negative regulator of T cell proliferation and cytokine responses. The co-stimulatory receptor pathways are potentially involved in self-tolerance and thus, they provide a promising therapeutic strategy for autoimmune diseases and transplantation. The genes encoding CD28, CTLA4 and ICOS are located adjacently in the chromosome region 2q33. The PDCD1 gene maps further, to the region 2q37. CTLA4 and PDCD1 are associated with the risk of a few autoimmune diseases. There is strong linkage disequilibrium (LD) on the 2q33 region; the whole gene of CD28 exists in its own LD block but CTLA4 and the 5' part of ICOS are within a same LD block. The 3' part of ICOS and PDCD1 are in their own separate LD blocks. Extended haplotypes covering the 2q33 region can be identified. This study focuses on immune related conditions like coeliac disease (CD) which is a chronic inflammatory disease with autoimmune features. Immunoglobulin A deficiency (IgAD) belongs to the group of primary antibody deficiencies characterised by reduced levels of immunoglobulins. IgAD co-occurs often with coeliac disease. Renal transplantation is needed in the end stage kidney diseases. Transplantation causes strong immune response which is tried to suppress with drugs. All these conditions are multifactorial with complex genetic background and multiple environmental factors affecting the outcome. We have screened ICOS for polymorphisms by sequencing the exon regions. We detected 11 new variants and determined their frequencies in Finnish population. We have measured linkage disequilibrium on the 2q33 region in Finnish as well as other European populations and observed conserved haplotypes. We analysed genetic association and linkage of the co-stimulatory receptor gene region aiming to study if it is a common risk locus for immune diseases. The 2q33 region was replicated to be linked to coeliac disease in Finnish population and CTLA4-ICOS haplotypes were found to be associated with CD and IgAD being the first non-HLA risk locus common for CD and immunodeficiencies. We also showed association between ICOS and the outcome of kidney transplantation. Our results suggest new evidence for CTLA4-ICOS gene region to be involved in susceptibility of coeliac disease. The earlier published contradictory association results can be explained by involvement of both CTLA4 and ICOS in disease susceptibility. The pattern of variants acting together rather than a single polymorphism may confer the disease risk. These genes may predispose also to immunodeficiencies as well as decreased graft survival and delayed graft function. Consequently, the present study indicates that like the well established HLA locus, the co-stimulatory receptor genes predispose to variety of immune disorders.

Phylogeography and Population Genetics of the Moor frog (Rana arvalis Nilsson) in Northern Europe

The ongoing climate change along with increasing levels of pollutants, diseases, habitat loss and fragmentation constitute global threats to the persistence of many populations, species and ecosystems. However, for the long-term persistence of local populations, one of the biggest threats is the intrinsic loss of genetic variation. In order to adapt to changes in the environment, organisms must have a sufficient supply of heritable variation in traits important for their fitness. With a loss of genetic variation, the risk of extinction will increase. For conservational practices, one should therefore understand the processes that shape the genetic population structure and also the broader (historical) phylogenetic patterning of the species in focus. In this thesis, microsatellite markers were applied to study genetic diversity and population differentiation of the protected moor frog (Rana arvalis) in Fennoscandia from both historical (evolutionary) and applied (conservation) perspectives. The results demonstrate that R. arvalis populations are highly structured over rather short geographic distances. Moreover, the results suggest that R. arvalis recolonized Fennoscandia from two directions after the last ice age. This has had implications for the genetic structuring and population differentiation, especially in the northernmost parts where the two lineages have met. Compared to more southern populations, the genetic variation decreases and the interpopulation differentiation increases dramatically towards north. This could be an outcome of serial population bottlenecking along the recolonization route. Also, current isolation and small population sizes increase the effect of drift, thus reinforcing the observed pattern. The same pattern can also be seen in island populations. However, though R. arvalis on the island of Gotland has lost most of its neutral genetic variability, our results indicate that the levels of additive genetic variation have remained high. This conforms to the conjecture that though neutral markers are widely used in conservation purposes, they may be quite uninformative about the levels of genetic variation in ecologically important traits. Finally, the evolutionary impact of the typical amphibian mating behaviour on genetic diversity was investigated. Given the short time available for larval development, it is important that mating takes place as early as possible. The genetic data and earlier capture-recapture data suggest that R. arvalis gather at mating grounds they are familiar with. However, by forming leks in random to relatedness, and having multiple paternities in single clutches, the risk of inbreeding may be minimized in this otherwise highly philopatric species.

Molecular genetics of Meckel syndrome : Ciliary genes are defective in MKS

Meckel syndrome (MKS, MIM 249000) is a severe developmental disorder that leads to death already in utero or shortly after birth. MKS diagnosis can be established by a careful ultrasound examination already at 11-14 weeks of gestation. The main features of MKS are occipital meningoencephalocele, cystic kidney dysplasia and fibrotic changes of the liver. In addition, polydactyly is frequently reported in the cases. The aim of the study was to characterize the molecular and functional defects in MKS. In this study we were able to identify two major MKS mutations in Finnish population, which cover over 90% of the cases. The first mutation is a 29 bp intronic deletion in the MKS1 gene (c.1483-7_35del) that is found in 70% of the families and the second is a C>T substitution in the coding region of CC2D2A (c.1762C>T), that is found in 20% of the MKS families. Both of these mutations result in abnormal splicing. The discovery of the disease genes has revealed that MKS is caused by primary cilia dysfunction. MKS1 gene has a conserved B9 domain, and it is found in the predicted ciliary proteome. CC2D2A protein is also found in the predicted ciliary proteome and it has a Ca2+ binding domain. The number of genes behind MKS has increased rapidly in the past years and to date, mutations have been identified in five genes (MKS1, TMEM67/MKS3, CEP290/MKS4, RPGRIP1L/MKS5 and CC2D2A/MKS6). Identification of the disease genes mutations has also revealed that MKS is an allelic disorder with other syndromes with overlapping phenotypes. Disorders that are caused by primary cilia dysfunction are collectively known as ciliopathies. Sequence analysis of all the known MKS genes in Finnish and non-Finnish families available to us, where the mutation was still unknown, revealed mutations in 14 out of the 30 families included in the study. When we collected all the reported mutations in MKS genes in different syndromes we could see that there was clearly a genotype-syndrome correlation between the mutations and the syndromes, since the same pair of mutations has never been reported in different syndromes. The basic molecular events behind MKS will not only give us information of this syndrome, but also significant novel information on early fetal development in general.

Genomic approach to organ-specific growth control

Growth is a fundamental aspect of life cycle of all organisms. Body size varies highly in most animal groups, such as mammals. Moreover, growth of a multicellular organism is not uniform enlargement of size, but different body parts and organs grow to their characteristic sizes at different times. Currently very little is known about the molecular mechanisms governing this organ-specific growth. The genome sequencing projects have provided complete genomic DNA sequences of several species over the past decade. The amount of genomic sequence information, including sequence variants within species, is constantly increasing. Based on the universal genetic code, we can make sense of this sequence information as far as it codes proteins. However, less is known about the molecular mechanisms that control expression of genes, and about the variations in gene expression that underlie many pathological states in humans. This is caused in part by lack of information about the second genetic code that consists of the binding specificities of transcription factors and the combinatorial code by which transcription factor binding sites are assembled to form tissue-specific and/or ligand-regulated enhancer elements. This thesis presents a high-throughput assay for identification of transcription factor binding specificities, which were then used to measure the DNA binding profiles of transcription factors involved in growth control. We developed ‘enhancer element locator’, a computational tool, which can be used to predict functional enhancer elements. A genome-wide prediction of human and mouse enhancer elements generated a large database of enhancer elements. This database can be used to identify target genes of signaling pathways, and to predict activated transcription factors based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c- and N-Myc genes, which has implications to organ specific growth control and tumor type specificity of oncogenes. Furthermore, we were able to locate a variation in a single nucleotide, which carries a susceptibility to colorectal cancer, to an enhancer element and propose a mechanism by which this SNP might be involved in generation of colorectal cancer.

Molecular genetics of tooth agenesis

Congenital missing of teeth, tooth agenesis or hypodontia, is one of the most common developmental anomalies in man. The common forms in which one or a few teeth are absent, may cause occlusal or cosmetic harm, while severe forms which are relatively rare always require clinical attention to support and maintain the dental function. Observation of tooth agenesis is also important for diagnosis of malformation syndromes. Some external factors may cause developmental defects and agenesis in dentition. However, the role of inheritance in the etiology of tooth agenesis is well established by twin and family studies. Studies on familial tooth agenesis as well as mouse null mutants have also identified several genetic factors. However, these explain syndromic or rare dominant forms of tooth agenesis, whereas the genes and defects responsible for the majority of cases of tooth agenesis, especially the common and less severe forms, are largely unknown. In this study it was shown, that a dominant nonsense mutation in PAX9 was responsible for severe tooth agenesis (oligodontia) in a Finnish family. In a study of tooth agenesis associated with Wolf-Hirschhorn syndrome, it was shown that severe tooth agenesis was present if the causative deletion in 4p spanned the MSX1 locus. It was concluded that severe tooth agenesis was caused by haploinsufficiency of these transcription factors. A summary of the phenotypes associated with known defects in MSX1 and PAX9 showed that, despite similarities, they were significantly different, suggesting that the genes, in addition to known interactions, also have independent roles during the development of human dentition. The original aim of this work was to identify gene defects that underlie the common incisor and premolar hypodontia. After excluding several candidate genes, a genome-wide search was conducted in seven Finnish families in which this phenotype was inherited in an autosomal dominant manner. A promising locus for second premolar agenesis was identified in chromosome 18 in one family and this finding was supported by results from other families. The results also implied the existence of other loci both for second premolar agenesis and for incisor agenesis. On the other hand the results did not lend support for comprehensive involvement of the most obvious candidate genes in the etiology of incisor and premolar hypodontia. Rather, they suggest remarkable genetic heterogeneity of tooth agenesis. The available evidence suggests that quantitative defects during tooth development predispose to a failure to overcome a developmental threshold and to agenesis. The results of the study increase the understanding of the etiology and heredity of tooth agenesis. Further studies may lead to identification of novel genes that affect the development of teeth.

Mu in vitro Transposition Technology in Functional Genetics and Genomics: Applications on Mouse and Bacteriophages

Transposons are mobile elements of genetic material that are able to move in the genomes of their host organisms using a special form of recombination called transposition. Bacteriophage Mu was the first transposon for which a cell-free in vitro transposition reaction was developed. Subsequently, the reaction has been refined and the minimal Mu in vitro reaction is useful in the generation of comprehensive libraries of mutant DNA molecules that can be used in a variety of applications. To date, the functional genetics applications of Mu in vitro technology have been subjected to either plasmids or genomic regions and entire genomes of viruses cloned on specific vectors. This study expands the use of Mu in vitro transposition in functional genetics and genomics by describing novel methods applicable to the targeted transgenesis of mouse and the whole-genome analysis of bacteriophages. The methods described here are rapid, efficient, and easily applicable to a wide variety of organisms, demonstrating the potential of the Mu transposition technology in the functional analysis of genes and genomes. First, an easy-to-use, rapid strategy to generate construct for the targeted mutagenesis of mouse genes was developed. To test the strategy, a gene encoding a neuronal K+/Cl- cotransporter was mutagenised. After a highly efficient transpositional mutagenesis, the gene fragments mutagenised were cloned into a vector backbone and transferred into bacterial cells. These constructs were screened with PCR using an effective 3D matrix system. In addition to traditional knock-out constructs, the method developed yields hypomorphic alleles that lead into reduced expression of the target gene in transgenic mice and have since been used in a follow-up study. Moreover, a scheme is devised to rapidly produce conditional alleles from the constructs produced. Next, an efficient strategy for the whole-genome analysis of bacteriophages was developed based on the transpositional mutagenesis of uncloned, infective virus genomes and their subsequent transfer into susceptible host cells. Mutant viruses able to produce viable progeny were collected and their transposon integration sites determined to map genomic regions nonessential to the viral life cycle. This method, applied here to three very different bacteriophages, PRD1, ΦYeO3 12, and PM2, does not require the target genome to be cloned and is directly applicable to all DNA and RNA viruses that have infective genomes. The method developed yielded valuable novel information on the three bacteriophages studied and whole-genome data can be complemented with concomitant studies on individual genes. Moreover, end-modified transposons constructed for this study can be used to manipulate genomes devoid of suitable restriction sites.

The genetic associations of acute anterior uveitis and their overlap with the genetics of ankylosing spondylitis

Acute anterior uveitis (AAU) involves inflammation of the iris and ciliary body of the eye. It occurs both in isolation and as a complication of ankylosing spondylitis (AS). It is strongly associated with HLA-B*27, but previous studies have suggested that further genetic factors may confer additional risk. We sought to investigate this using the Illumina Exomechip microarray, to compare 1504 cases with AS and AAU, 1805 with AS but no AAU and 21 133 healthy controls. We also used a heterogeneity test to test the differences in effect size between AS with AAU and AS without AAU. In the analysis comparing AS+AAU+ cases versus controls, HLA-B*27 and HLA-A*02:01 were significantly associated with the presence of AAU (P<10−300 and P=6 × 10−8, respectively). Secondary independent association with PSORS1C3 (P=4.7 × 10−5) and TAP2 (P=1.1 × 10−5) were observed in the major histocompatibility complex. There was a new suggestive association with a low-frequency variant at zinc-finger protein 154 in the AS without AAU versus control analysis (zinc-finger protein 154 (ZNF154), P=2.2 × 10−6). Heterogeneity testing showed that rs30187 in ERAP1 has a larger effect on AAU compared with that in AS alone. These findings also suggest that variants in ERAP1 have a differential impact on the risk of AAU when compared with AS, and hence the genetic risk for AAU differs from AS.

Genetics of ankylosing spondylitis - Insights into pathogenesis

Ankylosing spondylitis (AS), an immune-mediated arthritis, is the prototypic member of a group of conditions known as spondyloarthropathies that also includes reactive arthritis, psoriatic arthritis and enteropathic arthritis. Patients with these conditions share a clinical predisposition for spinal and pelvic joint dysfunction, as well as genetic associations, notably with HLA-B*27. Spondyloarthropathies are characterized by histopathological inflammation in entheses (regions of high mechanical stress where tendons and ligaments insert into bone) and in the subchondral bone marrow, and by abnormal osteoproliferation at involved sites. The association of AS with HLA-B*27, first described >40 years ago, led to hope that the cause of the disease would be rapidly established. However, even though many theories have been advanced to explain how HLA-B*27 is involved in AS, no consensus about the answers to this question has been reached, and no successful treatments have yet been developed that target HLA-B27 or its functional pathways. Over the past decade, rapid progress has been made in discovering further genetic associations with AS that have shed new light on the aetiopathogenesis of the disease. Some of these discoveries have driven translational ideas, such as the repurposing of therapeutics targeting the cytokines IL-12 and IL-23 and other factors downstream of this pathway. AS provides an excellent example of how hypothesis-free research can lead to major advances in understanding pathogenesis and to the development of innovative therapeutic strategies.

GDNF Receptors: Veterans and Novices

Neurotrophic factors play essential role in the development and functioning of the nervous system and other organs. Glial cell line-Derived Neurotrophic Factor (GDNF) family ligands (GFLs) are of particular interest because they promote the survival of dopaminergic neurons in vitro, in Parkinson s disease animal models and in patients. GDNF is also a potent survival factor for the central motoneurons and thus is considered as a potential lead for the treatment of amyotrophic lateral sclerosis. The survival promoting receptor complex for GFLs consists of a ligand-specific co-receptor, GFRα and a signal transducing module, receptor tyrosine kinase RET. At least GDNF and persephin, a GFL, have established functions outside central nervous system. GDNF is crucial for enteric nervous system and kidney development as well as for spermatogenesis. Persephin controls calcitonin secretion. Communication between cells often occurs in the extracellular matrix (ECM), a meshwork, which is secreted and deposited by the cells and is mainly composed of fibrillar proteins and polymerized sugars. We evaluated the relationship between GFLs and extracellular matrix components and demonstrated that three GFLs - GDNF, neurturin and artemin bind heparan sulfates with nanomolar affinities. The fourth member of the family - persephin binds these polysaccharides thousand times less tightly. GDNF, neurturin and artemin also bind with high affinity to heparan sulfate proteoglycan (HSPG) isolated from the nervous system, syndecan-3. GDNF signals through HSPGs, evoking Src family kinase activation. This signaling induces cell spreading, hippocampal neurite outgrowth in vitro and cellular migration. Specifically, GDNF signaling through syndecan-3 is important for embryonic cortical neuron migration. Syndecan-3-deficient mice, similarly to mice lacking GDNF, have less GABAergic neurons in their cortex, as compared to the wild-type mice. This fact provides indirect evidence that GDNF interaction with syndecan-3 is important for cortical brain development. Noteworthy, in non-neuronal tissues GFLs may signal via other syndecans. We also present the structural model for a GDNF co-receptor, GFRα1. The X-ray structure of the GFRα1 domain 3 was solved with 1.8 Å resolution, revealing a new protein fold. Later we also solved the structure of the truncated GFRα1 in the complex with GDNF and this model was confirmed by site-directed mutagenesis. In summary, our work contributed to the structural characterization of GFRα-based receptor complex and revealed a new receptor for GDNF, neurturin and artemin the HSPG syndecan-3. This information is critically important for the development of GFRα/RET agonists for the treatment of neurodegenerative diseases.

Evolutionary Genetics in the WIld : from Populations to Individuals

Predicting evolutionary outcomes and reconstructing past evolutionary transitions are among the main goals of evolutionary biology. Ultimately, understanding the mechanisms of evolutionary change will also provide answers to the timely question of whether and how organisms will adapt to changing environmental conditions. In this thesis, I have investigated the relative roles of natural selection, random genetic drift and genetic correlations in the evolution of complex traits at different levels of organisation from populations to individuals. I have shown that natural selection has been the driving force behind body shape divergence of marine and freshwater threespine stickleback (Gasterosteus aculeatus) populations, while genetic drift may have played a significant role in the more fine scale divergence among isolated freshwater populations. These results are concurrent with the patterns that have emerged in the published studies comparing the relative importance of natural selection and genetic drift as explanations for population divergence in different traits and taxa. I have also shown that body shape and armour divergence among threespine stickleback populations is likely to be biased by the patterns of genetic variation and covariation. Body shape and armour variation along the most likely direction of evolution the direction of maximum genetic variance reflects the general patterns of variation observed wild populations across the distribution range of the threespine stickleback. Conversely, it appears that genetic correlations between the sexes have not imposed significant constraints on the evolution of sexual dimorphism in threespine stickleback body shape and armour. I have demonstrated that the patterns of evolution seen in the wild can be experimentally recreated to tease out the effects of different selection agents in detail. In addition, I have shown how important it is to take into account the correlative nature of traits, when making interpretations about the effects of natural selection on individual traits. Overall, this thesis provides a demonstration of how considering the relative roles of different mechanism of evolutionary change at different levels of organisation can aid in an emergence of a comprehensive picture of how adaptive divergence in wild populations occurs.

The Emerging Phenotype: Ecological Genetics of Color, Form and Sex in the Common Frog

One of the main aims of evolutionary biology is to explain why organisms vary phenotypically as they do. Proximately, this variation arises from genetic differences and from environmental influences, the latter of which is referred to as phenotypic plasticity. Phenotypic plasticity is thus a central concept in evolutionary biology, and understanding its relative importance in causing the phenotypic variation and differentiation is important, for instance in anticipating the consequences of human induced environmental changes. The aim of this thesis was to study geographic variation and local adaptation, as well as sex ratios and environmental sex reversal, in the common frog (Rana temporaria). These themes cover three different aspects of phenotypic plasticity, which emerges as the central concept for the thesis. The first two chapters address geographic variation and local adaptation in two potentially thermally adaptive traits, namely the degree of melanism and the relative leg length. The results show that although there is an increasing latitudinal trend in the degree of melanism in wild populations across Scandinavian Peninsula, this cline has no direct genetic basis and is thus environmentally induced. The second chapter demonstrates that although there is no linear, latitudinally ordered phenotypic trend in relative leg length that would be expected under Allen s rule an ecogeographical rule linking extremity length to climatic conditions there seems to be such a trend at the genetic level, hidden under environmental effects. The first two chapters thus view phenotypic plasticity through its ecological role and evolution, and demonstrate that it can both give rise to phenotypic variation and hide evolutionary patterns in studies that focus solely on phenotypes. The last three chapters relate to phenotypic plasticity through its ecological and evolutionary role in sex determination, and consequent effects on population sex ratio, genetic recombination and the evolution of sex chromosomes. The results show that while sex ratios are strongly female biased and there is evidence of environmental sex reversals, these reversals are unlikely to have caused the sex ratio skew, at least directly. The results demonstrate that environmental sex reversal can have an effect on the evolution of sex chromosomes, as the recombination patterns between them seem to be controlled by phenotypic, rather than genetic, sex. This potentially allows Y chromosomes to recombine, lending support for the recent hypothesis suggesting that sex-reversal may play an important role on the rejuvenation of Y chromosomes.