Molecular basis of colorectal cancer predisposition

Colorectal cancer (CRC) is one of the most frequent malignancies in Western countries. Inherited factors have been suggested to be involved in 35% of CRCs. The hereditary CRC syndromes explain only ~6% of all CRCs, indicating that a large proportion of the inherited susceptibility is still unexplained. Much of the remaining genetic predisposition for CRC is probably due to undiscovered low-penetrance variations. This study was conducted to identify germline and somatic changes that contribute to CRC predisposition and tumorigenesis. MLH1 and MSH2, that underlie Hereditary non-polyposis colorectal cancer (HNPCC) are considered to be tumor suppressor genes; the first hit is inherited in the germline and somatic inactivation of the wild type allele is required for tumor initiation. In a recent study, frequent loss of the mutant allele in HNPCC tumors was detected and a new model, arguing against the two-hit hypothesis, was proposed for somatic HNPCC tumorigenesis. We tested this hypothesis by conducting LOH analysis on 25 colorectal HNPCC tumors with a known germline mutation in the MLH1 or MSH2 genes. LOH was detected in 56% of the tumors. All the losses targeted the wild type allele supporting the classical two-hit model for HNPCC tumorigenesis. The variants 3020insC, R702W and G908R in NOD2 predispose to Crohn s disease. Contribution of NOD2 to CRC predisposition has been examined in several case-control series, with conflicting results. We have previously shown that 3020insC does not predispose to CRC in Finnish CRC patients. To expand our previous study the variants R702W and G908R were genotyped in a population-based series of 1042 Finnish CRC patients and 508 healthy controls. Association analyses did not show significant evidence for association of the variants with CRC. Single nucleotide polymorphism (SNP) rs6983267 at chromosome 8q24 was the first CRC susceptibility variant identified through genome-wide association studies. To characterize the role of rs6983267 in CRC predisposition in the Finnish population, we genotyped the SNP in the case-control material of 1042 cases and 1012 controls and showed that G allele of rs6983267 is associated with the increased risk of CRC (OR 1.22; P=0.0018). Examination of allelic imbalance in the tumors heterozygous for rs6983267 revealed that copy number increase affected 22% of the tumors and interestingly, it favored the G allele. By utilizing a computer algorithm, Enhancer Element Locator (EEL), an evolutionary conserved regulatory motif containing rs6983267 was identified. The SNP affected the binding site of TCF4, a transcription factor that mediates Wnt signaling in cells, and has proven to be crucial in colorectal neoplasia. The preferential binding of TCF4 to the risk allele G was showed in vitro and in vivo. The element drove lacZ marker gene expression in mouse embryos in a pattern that is consistent with genes regulated by the Wnt signaling pathway. These results suggest that rs6983267 at 8q24 exerts its effect in CRC predisposition by regulating gene expression. The most obvious target gene for the enhancer element is MYC, residing ~335 kb downstream, however further studies are required to establish the transcriptional target(s) of the predicted enhancer element.

Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and approximately 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 x 10(-)(8)), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

Microarrays in Molecular Profiling of Ewing Sarcoma Family of Tumors and CDKN2A Aberrations

Background: The Ewing sarcoma family of tumors (ESFT) are rare but highly malignant neoplasms that occur mainly in bone or but also in soft tissue. ESFT affects patients typically in their second decade of life, whereby children and adolescents bear the heaviest incidence burden. Despite recent advances in the clinical management of ESFT patients, their prognosis and survival are still disappointingly poor, especially in cases with metastasis. No targeted therapy for ESFT patients is currently available. Moreover, based merely on current clinical and biological characteristics, accurate classification of ESFT patients often fails at the time of diagnosis. Therefore, there is a constant need for novel molecular biomarkers to be applied in tandem with conventional parameters to further intensify ESFT risk-stratification and treatment selection, and ultimately to develop novel targeted therapies. In this context, a greater understanding of the genetics and immune characteristics of ESFT is needed. Aims: This study sought to open novel insights into gene copy number changes and gene expression in ESFT and, further, to enlighten the role of inflammation in ESFT. For this purpose, microarrays were used to provide gene-level information on a genomewide scale. In addition, this study focused on screening of 9p21.3 deletion sizes and frequencies in ESFT and, in another pediatric cancer, acute lymphocytic leukemia (ALL), in order to define more exact criteria for highrisk patient selection and to provide data for developing a more reliable diagnostic method to detect CDKN2A deletions. Results: In study I, 20 novel ESFT-associated suppressor genes and oncogenes were pinpointed using combined array CGH and expression analysis. In addition, interesting chromosomal rearrangements were identified: (1) Duplication of derivative chromosome der(22)(11;22) was detected in three ESFT patients. This duplication included the EWSR1-FLI1 fusion gene leading to increase in its copy number; (2) Cryptic amplifications on chromosomes 20 and 22 were detected, suggesting a novel translocation between chromosomes 20 and 22, which most probably produces a fusion between EWSR1 and NFATC2. In study II, bioinformatic analysis of ESFT expression profiles showed that inflammatory gene activation is detectable in ESFT patient samples and that the activation is characterized by macrophage gene expression. Most interestingly, ESFT patient samples were shown to express certain inflammatory genes that were prognostically significant. High local expression of C5 and JAK1 at the tumor site was shown to associate with favorable clinical outcome, whereas high local expression of IL8 was shown to be detrimental. Studies III and IV showed that the smallest overlapping region of deletion in 9p21.3 includes CDKN2A in all cases and that the length of this region is 12.2 kb in both Ewing sarcoma and ALL. Furthermore, our results showed that the most widely used commercial CDKN2A FISH probe creates false negative results in the narrowest microdeletion cases (<190 kb). Therefore, more accurate methods should be developed for the detection of deletions in the CDKN2A locus. Conclusions: This study provides novel insights into the genetic changes involved in the biology of ESFT, in the interaction between ESFT cells and immune system, and in the inactivation of CDKN2A. Novel ESFT biomarker genes identified in this study serve as a useful resource for future studies and in developing novel therapeutic strategies to improve the survival of patients with ESFT.

White shark (Carcharodon carcharias) microsatellite genotypes (six loci) from regions along Australia's southern and eastern coastline.

Despite international protection of white sharks (Carcharodon carcharias), important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with six nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (FST = 0.142, p < 0.001), implying female natal philopatry. This concords with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited microsatellite markers (FST = 0.009, p <0.05), suggesting that males may also exhibit some degree of reproductive philopatry. Five sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal or migration resulting in breeding may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and below the threshold at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least one, possibly two orders of magnitude below our historical effective size estimates. Further population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail. Reference: Blower, D. C., Pandolfi, J. M., Gomez-Cabrera, M. del C., Bruce, B. D. & Ovenden, J. R. (In press - April 2012). Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes. Marine Ecology Progress Series.

Pleistocene isolation, secondary introgression and restricted contemporary gene flow in the pig-eye shark, Carcharhinus amboinensis across northern Australia

We examine the structure and phylogeography of the pig-eye shark (Carcharhinus amboinensis) common in shallow coastal environments in northern Australia using two types of genetic markers, two mitochondrial (control region and NADH hydrogenase 4) and two nuclear (microsatellite and Rag 1) DNA. Two populations were defined within northern Australia on the basis of mitochondrial DNA evidence, but this result was not supported by nuclear microsatellite or Rag 1 markers. One possibility for this structure might be sex-specific behaviours such as female philopatry, although we argue it is doubtful that sufficient time has elapsed for any potential signatures from this behaviour to be expressed in nuclear markers. It is more likely that the observed pattern represents ancient populations repeatedly isolated and connected during episodic sea level changes during the Pleistocene epoch, until current day with restricted contemporary gene flow maintaining population genetic structure. Our results show the need for an understanding of both the history and ecology of a species in order to interpret patterns in genetic structure.

Population Expansion and Genetic Structure in Carcharhinus brevipinna in the Southern Indo-Pacific

Background:Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa. Methods and Findings:Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and FST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717–0.03508, p values ≤ 0.0013; pairwise FST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia’s east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).Conclusions: The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

An improved real-time PCR assay for the detection of Old World screwworm flies

The Old World screwworm (OWS) fly, Chrysomya bezziana, is a serious pest of livestock, wildlife and humans in tropical Africa, parts of the Middle East, the Indian subcontinent, south-east Asia and Papua New Guinea. Although to date Australia remains free of OWS flies, an incursion would have serious economic and animal welfare implications. For these reasons Australia has an OWS fly preparedness plan including OWS fly surveillance with fly traps. The recent development of an improved OWS fly trap and synthetic attractant and a specific and sensitive real-time PCR molecular assay for the detection of OWS flies in trap catches has improved Australia's OWS fly surveillance capabilities. Because all Australian trap samples gave negative results in the PCR assay, it was deemed necessary to include a positive control mechanism to ensure that fly DNA was being successfully extracted and amplified and to guard against false negative results. A new non-competitive internal amplification control (IAC) has been developed that can be used in conjunction with the OWS fly PCR assay in a multiplexed single-tube reaction. The multiplexed assay provides an indicator of the performance of DNA extraction and amplification without greatly increasing labour or reagent costs. The fly IAC targets a region of the ribosomal 16S mitochondrial DNA which is conserved across at least six genera of commonly trapped flies. Compared to the OWS fly assay alone, the multiplexed OWS fly and fly IAC assay displayed no loss in sensitivity or specificity for OWS fly detection. The multiplexed OWS fly and fly IAC assay provides greater confidence for trap catch samples returning negative OWS fly results. © 2014 International Atomic Energy Agency.

Plasmopara sphagneticolae sp. nov. (Peronosporales) on Sphagneticola (Asteraceae) in Australia

A specimen of downy mildew on leaves of Sphagneticola trilobata found in northern Queensland was identified by a systematic approach as a novel species of Plasmopara. A new species, Plasmopara sphagneticolae, is proposed for this specimen, which differs from other species of Plasmopara by morphology, host range, and sequence data from nuclear-ribosomal DNA and mitochondrial DNA. Plasmopara sphagneticolae, together with P. halstedii, are downy mildews found on host species in the tribe Heliantheae (Asteraceae). Plasmopara halstedii causes downy mildew on Helianthus annuus, and is not present on sunflower in Australia. Phylogenetic analysis of the large subunit region of ribosomal DNA showed that P. sphagneticolae was sister to P. halstedii on sunflower.

Population genetics of Australian white sharks reveals fine-scale spatial structure, transoceanic dispersal events and low effective population sizes

Despite international protection of white sharks Carcharodon carcharias, important conservation parameters such as abundance, population structure and genetic diversity are largely unknown. The tissue of 97 predominately juvenile white sharks sampled from spatially distant eastern and southwestern Australian coastlines was sequenced for the mitochondrial DNA (mtDNA) control region and genotyped with 6 nuclear-encoded microsatellite loci. MtDNA population structure was found between the eastern and southwestern coasts (F-ST = 0.142, p < 0.0001), implying female reproductive philopatry. This concurs with recent satellite and acoustic tracking findings which suggest the sustained presence of discrete east coast nursery areas. Furthermore, population subdivision was found between the same regions with biparentally inherited micro satellite markers (F-ST = 0.009, p < 0.05), suggesting that males may also exhibit some degree of reproductive philopatry; 5 sharks captured along the east coast had mtDNA haplotypes that resembled western Indian Ocean sharks more closely than Australian/New Zealand sharks, suggesting that transoceanic dispersal, or migration resulting in breeding, may occur sporadically. Our most robust estimate of contemporary genetic effective population size was low and close to thresholds at which adaptive potential may be lost. For a variety of reasons, these contemporary estimates were at least 1, possibly 2, orders of magnitude below our historical effective size estimates. Population decline could expose these genetically isolated populations to detrimental genetic effects. Regional Australian white shark conservation management units should be implemented until genetic population structure, size and diversity can be investigated in more detail.

Integrating different approaches in the definition of biological stocks: A northern Australian multi-jurisdictional fisheries example using grey mackerel, Scomberomorus semifasciatus

A holistic approach to stock structure studies utilises multiple different techniques on the same individuals sampled from selected populations and combines results across spatial and temporal scales to produce a weight of evidence conclusion. It is the most powerful and reliable source of information to use in formulating resource management and monitoring plans. Few examples of the use of a holistic approach in stock structure studies exist, although more recently this is changing. Using such an approach makes integration of results from each technique challenging. An integrated stock definition (ISD) approach for holistic stock structure studies was developed in this study to aid in the appropriate interpretation of stock structure results to guide the determination of fishery management units. The ISD approach is applied herein to a study of the northern Australian endemic grey mackerel, Scomberomorus semifasciatus (Scombridae). Analyses of genetic (mitochondrial DNA and microsatellites), parasite, otolith stable isotope, and growth data are synthesised to determine the stock structure of S. semifasciatus across northern Australia. Integration of the results from all techniques identified at least six S. semifasciatus stocks for management purposes. Further, the use of the ISD approach provided a simple basis for integrating multiple techniques and for their interpretation. The use of this holistic approach was a powerful tool in providing greater certainty about the appropriate management units for S. semifasciatus. Future stock structure studies investigating spatial management questions in the fisheries context should adopt a holistic approach and apply the ISD approach for a more accurate definition of biological stocks to improve fisheries management.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.

A co-evolutionary relationship exists between Endoraecium (Pucciniales) and its Acacia hosts in Australia

Endoraecium is a genus of rust fungi that infects several species of Acacia in Australia, South-East Asia and Hawaii. This study investigated the systematics of Endoraecium from 55 specimens in Australia based on a combined morphological and molecular approach. Phylogenetic analyses were conducted on partitioned datasets of loci from ribosomal and mitochondrial DNA. The recovered molecular phylogeny supported a recently published taxonomy based on morphology and host range that divided Endoraecium digitatum into five species. Spore morphology is synapomorphic and there is evidence Endoraecium co-evolved with its Acacia hosts. The broad host ranges of E. digitatum, E. parvum, E. phyllodiorum and E. violae-faustiae are revised in light of this study, and nine new species of Endoraecium are described from Australia based on host taxonomy, morphology and phylogenetic concordance.

Computational Identification of Recessive Mutations in Cancers using High Throughput SNP-arrays

This thesis presents a highly sensitive genome wide search method for recessive mutations. The method is suitable for distantly related samples that are divided into phenotype positives and negatives. High throughput genotype arrays are used to identify and compare homozygous regions between the cohorts. The method is demonstrated by comparing colorectal cancer patients against unaffected references. The objective is to find homozygous regions and alleles that are more common in cancer patients. We have designed and implemented software tools to automate the data analysis from genotypes to lists of candidate genes and to their properties. The programs have been designed in respect to a pipeline architecture that allows their integration to other programs such as biological databases and copy number analysis tools. The integration of the tools is crucial as the genome wide analysis of the cohort differences produces many candidate regions not related to the studied phenotype. CohortComparator is a genotype comparison tool that detects homozygous regions and compares their loci and allele constitutions between two sets of samples. The data is visualised in chromosome specific graphs illustrating the homozygous regions and alleles of each sample. The genomic regions that may harbour recessive mutations are emphasised with different colours and a scoring scheme is given for these regions. The detection of homozygous regions, cohort comparisons and result annotations are all subjected to presumptions many of which have been parameterized in our programs. The effect of these parameters and the suitable scope of the methods have been evaluated. Samples with different resolutions can be balanced with the genotype estimates of their haplotypes and they can be used within the same study.

The life cycle and genetic structure of the red alga Furcellaria lumbricalis on a salinity gradient

The life cycle and genetic diversity of the red alga Furcellaria lumbricalis (Hudson) Lamouroux were investigated in 15 populations in northern Europe. The occurrence of different life cycle phases and seasonality of reproduction were studied in four brackish populations in the northern Baltic Sea. Furthermore, a new method, based on genome screening with ISSR markers combined with a restriction-ligation method, was developed to discover microsatellite markers for population genetic analyses. The mitochondrial DNA cox2-3 spacer sequence and four microsatellite markers were used to examine the genetic diversity and differentiation of red algal populations in northern Europe. In addition, clonality and small-scale genetic structure of one Irish and four Baltic Sea populations were studied with microsatellite markers. It was discovered that at the low salinities of the northern Baltic Sea, only tetrasporophytes and males were present in the populations of F. lumbricalis and that winter was the main season for tetrasporangial production. Furthermore, the population occurring at the lowest salinity (3.6 practical salinity units, psu) did not produce spores. The size of the tetraspores was smaller in the Baltic Sea populations than that in the Irish population, and there were more deformed spores in the Baltic Sea populations than in the Irish populations. Studies with microsatellite markers indicated that clonality is a common phenomenon in the Baltic Sea populations of F. lumbricalis, although the proportion of clonal individuals varied among populations. Some genetic divergence occurred within locations both in Ireland and in the northern Baltic Sea. Even though no carpogonia were detected in the field samples during the study, the microsatellite data indicated that sexual reproduction occurs at least occasionally in the northern Baltic Sea. The genetic diversity of F. lumbricalis was highest in Brittany, France. Since no variation was discovered in the mtDNA cox2-3 spacer sequence, which is generally regarded as an informative phylogeographic marker in red algae, it can be assumed that the studied populations probably share the same origin.