Genetic alterations in microsatellite-unstable colorectal cancer

Colorectal cancer (CRC) is the third most common cancer in Finland. Of all CRC tumors, 15% display microsatellite-instability (MSI) caused by defective cellular mismatch repair. Cells displaying MSI accumulate a high number of mutations genome-wide, especially in short repeat areas, microsatellites. When targeting genes essential for cell growth or death, MSI can promote tumorigenesis. In non-coding areas, microsatellite mutations are generally considered as passenger events. Since the discovery of MSI and its linkage to cancer, more that 200 genes have been investigated for a role in MSI tumorigenesis. Although various criteria have been suggested for MSI target gene identification, the challenge has been to distinguish driver mutations from passenger mutations. This study aimed to clarify these key issues in the research field of MSI cancer. Prior to this, background mutation rate in MSI cancer has not been studied in a large-scale. We investigated the background mutation rate in MSI CRC by analyzing the spectrum of microsatellite mutations in non-coding areas. First, semenogelin I was studied for a possible role in MSI carcinogenesis. The intronic T9 repeat of semenogelin I was frequently mutated but no evidence for selection during tumorigenesis was obtained. Second, a sequencing approach was utilized to evaluate the general background mutation rate in MSI CRC. Both intronic and intergenic repeats harbored extremely high mutation rates of ≤ 87% and intergenic repeats were more unstable than the intronic repeats. As mutation rates of presumably neutral microsatellites can be high in MSI CRC in the absence of apparent selection pressure, high mutation frequency alone is not sufficient evidence for identification of driver MSI target genes. Next, an unbiased approach was designed to identify the mutatome of MSI CRC. By combining expression array data and a database search we identified novel genes possibly related to MSI CRC carcinogenesis. One of the genes was studied further. In the functional analysis this gene was observed to cause an abnormal cancer-prone cellular phenotype, possibly through altered responses to DNA damage. In our recent study, smooth muscle myosin heavy chain 11 (MYH11) was identified as a novel MSI CRC gene. Additionally, MYH11 has a well established role in acute myeloid leukemia (AML) through an oncogenic fusion protein CBFB-MYH11. We investigated further the role of MYH11 in AML by sequencing. Three novel missense variants of MYH11 were identified. None of the variants were present in the population-based control material. One of the identified variants, V71A, lies in the N-terminal SH3-like domain of MYH11 of unknown function. The other two variants, K1059E and R1792Q are located in the coil-coiled myosin rod essential for the regulation and filament formation of MYH11. The variant K1059E lies in the close proximity of the K1044N that has been functionally assessed in our earlier work of CRC and has been reported to cause total loss of MYH11 protein regulation. As the functional significance of the three novel variants examined in this work remains unknown, future studies should clarify the further role of MYH11 in AML leukaemogenesis and in other malignancies.

Rapid microsatellite marker development for African mahogany ( Khaya senegalensis, Meliaceae) using next-generation sequencing and assessment of its intra-specific genetic diversity.

Khaya senegalensis (African mahogany or dry-zone mahogany) is a high-value hardwood timber species with great potential for forest plantations in northern Australia. The species is distributed across the sub-Saharan belt from Senegal to Sudan and Uganda. Because of heavy exploitation and constraints on natural regeneration and sustainable planting, it is now classified as a vulnerable species. Here, we describe the development of microsatellite markers for K. senegalensis using next-generation sequencing to assess its intra-specific diversity across its natural range, which is a key for successful breeding programs and effective conservation management of the species. Next-generation sequencing yielded 93943 sequences with an average read length of 234bp. The assembled sequences contained 1030 simple sequence repeats, with primers designed for 522 microsatellite loci. Twenty-one microsatellite loci were tested with 11 showing reliable amplification and polymorphism in K. senegalensis. The 11 novel microsatellites, together with one previously published, were used to assess 73 accessions belonging to the Australian K. senegalensis domestication program, sampled from across the natural range of the species. STRUCTURE analysis shows two major clusters, one comprising mainly accessions from west Africa (Senegal to Benin) and the second based in the far eastern limits of the range in Sudan and Uganda. Higher levels of genetic diversity were found in material from western Africa. This suggests that new seed collections from this region may yield more diverse genotypes than those originating from Sudan and Uganda in eastern Africa.

New SNPs for population genetic analysis reveal possible cryptic speciation of eastern Australian sea mullet (Mugil cephalus)

Sustainable management of sea mullet (Mugil cephalus) fisheries needs to account for recent observations of regional-scale differentiation. Population genetic analysis is sought to assess the situation of this ecologically and economically important fish species in eastern Australian waters. Here, we report (i) new population genetic markers [single nucleotide polymorphisms (SNPs) and potential microsatellites], (ii) first estimates of spatial genetic differentiation and (iii) prospective power tests for designing more comprehensive studies. Six DNA samples from three sampling regions (North Queensland, South Queensland and central New South Wales) on the eastern coast of Australia were used to prepare restriction site associated DNA (RAD) tag libraries from genomic DNA digested with EcoRI and MseI. A pooled sample of regional RAD tag libraries was sequenced using the Roche GS-FLX Titanium platform. A total of 172837 raw reads (17.4Mbp) were retrieved, 95500 of which were used to discover 1267 SNPs and 1417 microsatellites. A subset of 161 SNPs was validated based on 63 additional DNA samples genotyped using the Sequenom MassArray (iPLEX Gold chemistry). Altogether 92 SNPs (57%) were confirmed, with 40% of these marking fixed variants between northern and southern sampling regions. Our preliminary findings indicate a multispecies fishery stock of M. cephalus in eastern Australian waters, but suggest that strong genetic differentiation occurs north of major fishing grounds. Low potential differentiation within major fishing grounds (e.g. FST=0.0025) can be resolved with a likely power 67% by using standard sample sizes of 50 and validated subsets of available markers.

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.

Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

Natural variation and Biogeography of the melon fruit fly, Zeugodacus cucurbitae (Diptera : Tephritidae), in Southeast-Asia and the west-pacific

This thesis used multidisciplinary approaches which greatly enhance our understanding of population structure and can be particularly powerful tools for resolving variation of melon fly over geographic and temporal scales, and for determining invasive pathways. The results from this thesis reinforce the value of integrating multiple data sets to better understand and resolve natural variation within an important pest to determine whether there are cryptic species, discrete lineages or host races, and to identify dispersal pathways in an invasive pest. These results are instructive for regional biosecurity, trade and quarantine, and provide important background for future area-wide management programmes. The integrative methodology adopted in this thesis is applicable to a variety of other insect pests.

Reproduction and population structure in European aspen

The European aspen (Populus tremula) is a keystone species for biodiversity in boreal forests. However, the future of aspen may be threatened, because large aspens have mostly been removed from managed forests, whereas regeneration and the long-term persistence of mature trees are subjects of concern in protected areas. Aspen is a pioneer tree, and it can reproduce both sexually by seed and asexually by root suckers. Through asexual reproduction aspen forms clones, groups of genetically identical trees (ramets). In my thesis, I have studied the structure of aspen populations in terms of number, size, clonal and demographic properties. Additionally, I have investigated the emergence and survival of seedlings as well as the seed quantity and quality in crosses between the European and hybrid aspen. To study the regeneration and population structure, mature aspens were recorded in old-growth and managed forests in eastern Finland based on a large-scale inventory (11 400 ha). In addition, small aspen trees were surveyed on sample plots. Clonal structure was investigated both by morphological characters and by DNA-based markers (microsatellites). Seedling emergence and survival was studied with two sowing experiments. With crosses between European and hybrid aspens we wanted to study whether elevated temperatures due to climate change would benefit the different crosses of European and hybrid aspen unequally and thus affect the gene flow between the two species. The average volumes of mature aspen were 5.3 m3/ha in continuous old-growth, and 0.8 m3/ha in managed forests. Results indicate also that large aspen trees in managed forests are a legacy of the past less intensively managed forest landscapes. Long-term persistence of aspen in protected areas can only be secured by restoration measures creating sufficiently large gaps for regeneration. More emphasis should be given to sparing aspens in thinnings and to retaining of mature aspens in regeneration cutting in managed forests. Aspen was found to be spatially aggregated in the landscape. This could be explained by site type, disturbance history and / or limitations in seed dispersal. Clonal structure does not explain the spatial aggregation, since average size of the clones was only 2.3 ramets, and most clones (70 %) consisted of just one ramet. The small size of the clones suggests that most of them are relatively young. Therefore, sexual reproduction may be more common than has previously been thought. Seedling emergence was most successful in mineral soil especially, when the site had been burned. Only few seedlings occurred on humus. Survival of the seedlings was low, and strongly dependent on moisture, but also on seedbed conditions. The seeds were found to maintain their germinability longer than has earlier been thought to be possible. Interspecific crosses produced more seeds with higher quality than intraspecific crosses. When temperature was elevated, germination of hybrid aspen seeds increased more than seeds from P. tremula x P. tremula crosses. These results suggest that hybrid aspen may have a significant genetic impact on the European aspen, and this effect may become strengthened by climate warming.

Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.

Population genetics of the invasive plant species Impatiens glandulifera in Southern Finland

Biological invasions affect biodiversity worldwide, and, consequently, the invaded ecosystems may suffer from significant losses in economic and cultural values. Impatiens glandulifera Royle (Balsaminaceae) is an invasive annual herb, native to the western Himalayas and introduced into Europe in the 19th century as a garden ornamental plant. The massive invasion of I. glandulifera is due to its high reproductive output, rapid growth and its ability to outcompete native species. In Finland, the first observations regarding the presence of I. glandulifera date from the year 1947, and today it is considered a serious problem in riparian habitats. The aim of this master’s thesis research is to reveal the population genetic structure of I. glandulifera in Finland and to find out whether there have been one or multiple invasions in Finland. The study focuses on investigating the origin of I. glandulifera in Southern Finland, by comparing plant samples from the Helsinki region with those from its native region and other regions of invasion. Samples from four populations in Helsinki and from the United Kingdom, Canada, India and Pakistan were collected and genotyped using 11 microsatellite markers. The genetic analyses were evaluated using the programs Arlequin and Structure. The results of the genetic analyses suggested that I. glandulifera has been introduced to Finland more than once. Multiple introductions are supported by the higher level of genetic diversity detected within and among Finnish populations than would be expected for a single introduction. Results of the Bayesian Structure analysis divided the four Finnish populations into four clusters. This geographical structure was further supported by pairwise Fst values among populations. The causes and potential consequences of such multiple introductions of I. glandulifera in Finland and further perspectives are discussed.

The many ways to delimit species: hairs, genes and surface chemistry

Species identification forms the basis for understanding the diversity of the living world, but it is also a prerequisite for understanding many evolutionary patterns and processes. The most promising approach for correctly delimiting and identifying species is to integrate many types of information in the same study. Our aim was to test how cuticular hydro- carbons, traditional morphometrics, genetic polymorphisms in nuclear markers (allozymes and DNA microsatellites) and DNA barcoding (partial mitochondrial COI gene) perform in delimiting species. As an example, we used two closely related Formica ants, F. fusca and F. lemani, sampled from a sympatric population in the northern part of their distribu- tion. Morphological characters vary and overlap in different parts of their distribution areas, but cuticular hydrocarbons include a strong taxonomic signal and our aim is to test the degree to which morphological and genetic data correspond to the chemical data. In the morphological analysis, species were best separated by the combined number of hairs on pro- notum and mesonotum, but individual workers overlapped in hair numbers, as previously noted by several authors. Nests of the two species were separated but not clustered according to species in a Principal Component Analysis made on nuclear genetic data. However, model-based Bayesian clustering resulted in perfect separation of the species and gave no indication of hybridization. Furthermore, F. lemani and F. fusca did not share any mitochondrial haplotypes, and the species were perfectly separated in a phylogenetic tree. We conclude that F. fusca and F. lemani are valid species that can be separated in our study area relatively well with all methods employed. However, the unusually small genetic differen- tiation in nuclear markers (FST = 0.12) shows that they are closely related, and occasional hybridization between F. fusca and F. lemani cannot be ruled out.

Phylogeography unplugged: comparative surveys in the genomic era

In March 2012, the authors met at the National Evolutionary Synthesis Center (NESCent) in Durham, North Carolina, USA, to discuss approaches and cooperative ventures in Indo-Pacific phylogeography. The group emerged with a series of findings: (1) Marine population structure is complex, but single locus mtDNA studies continue to provide powerful first assessment of phylogeographic patterns. (2) These patterns gain greater significance/power when resolved in a diversity of taxa. New analytical tools are emerging to address these analyses with multi-taxon approaches. (3) Genome-wide analyses are warranted if selection is indicated by surveys of standard markers. Such indicators can include discordance between genetic loci, or between genetic loci and morphology. Phylogeographic information provides a valuable context for studies of selection and adaptation. (4) Phylogeographic inferences are greatly enhanced by an understanding of the biology and ecology of study organisms. (5) Thorough, range-wide sampling of taxa is the foundation for robust phylogeographic inference. (6) Congruent geographic and taxonomic sampling by the Indo-Pacific community of scientists would facilitate better comparative analyses. The group concluded that at this stage of technology and software development, judicious rather than wholesale application of genomics appears to be the most robust course for marine phylogeographic studies. Therefore, our group intends to affirm the value of traditional (''unplugged'') approaches, such as those based on mtDNA sequencing and microsatellites, along with essential field studies, in an era with increasing emphasis on genomic approaches.

Rapid divergence of ecotypes of an invasive plant

Invasive species demonstrate rapid evolution within a very short period of time allowing one to understand the underlying mechanism(s). Lantana camara, a highly invasive plant of the tropics and subtropics, has expanded its range and successfully established itself almost throughout India. In order to uncover the processes governing the invasion dynamics, 218 individuals from various locations across India were characterized with six microsatellites. By integrating genetic data with niche modelling, we examined the effect of drift and environmental selection on genetic divergence. We found multiple genetic clusters that were non-randomly distributed across space. Spatial autocorrelation revealed a strong fine-scale structure, i.e. isolation by distance. In addition, we obtained evidence of inhibitory effects of selection on gene flow, i.e. isolation by environmental distance. Perhaps, local adaptation in response to selection is offsetting gene flow and causing the populations to diverge. Niche models suggested that temperature and precipitation play a major role in the observed spatial distribution of this plant. Based on a non-random distribution of clusters, unequal gene flow among them and different bioclimatic niche requirements, we concluded that the emergence of ecotypes represented by two genetic clusters is underway. They may be locally adapted to specific climatic conditions, and perhaps at the very early stages of ecological divergence.

Multiple SNP Markers Reveal Fine-Scale Population and Deep Phylogeographic Structure in European Anchovy (Engraulis encrasicolus L.)

10 p.

Aquaculture Asia, Vol.14, No.2, pp.1-40, April-June 2009

Sustainable aquaculture GLOBALG.A.P. standard in Thai shrimp farms: Mission (im)possible? p.4 Leepaisomboon, T., Chuchird, N., Limsuwan, C., Steenbruggen, E.R., and Mungkung, R. The Victorian trout industry & the bushfires p.6 Mosig, J. Small-scale aquaculture in the Ayeyarwady Delta, Myanmar p.10 Edwards, P. The history, status, and future prospects of monosex tilapia culture in Thailand p.18 Belton, B., Turongruang, D., Bhujel, R. and Little, D.C. Mangroves of Nakhon Si Thammarat Province in southern Thailand: Species diversity, community structure and current status p.20 By Amarasinghe, M.D., Dulyapurk, V., Taparhudee, W., Yoonpundh, R. and Jumnongsong, S. Research and farming techniques Induced breeding of pacu (Piaractus brachypomus) in captivity with pituitary extract p.23 Chattarjee, N.R. and Mazumdar, B. Aquatic animal health Fumonisins - mycotoxins of increasing importance in fish! p. 24 Griessler, K. and Encarnação, P. Genetics and biodiversity Microsatellite DNA markers, a fisheries perspective. Part 1: The nature of microsatellites p.27 Sekar, M., Suresh, E., Kumar, N.S., Nayak, S.K., Balakrishna, C. Asia-Pacific Marine Finfi sh Aquaculture Magazine Formulated feed for tiger grouper grow-out p.30 Rachmansyah, Usman, Palinggi, N.N. and Williams, K. NACA Newsletter 36

Estrutura populacional e filogeografia de Panulirus argus (Latreille, 1804)

Panulirus argus (Latreille, 1804) é uma das principais espécies de lagosta no Atlântico, sendo um dos maiores recursos pesqueiros do Atlântico Ocidental, onde apresenta um alto valor comercial. A forte explotação da espécie resulta em uma grande pressão sobre suas populações. Recentemente, foi descoberto que sob o binômio P. argus estão contidas duas espécies crípticas que ocorrem em alopatria, uma na região do Caribe e outra na costa brasileira. Esta tese tem como objetivo estudar como se estruturam geneticamente as populações dessas duas espécies, com o propósito de fornecer mais informações para a determinação de estoques e um correto manejo das espécies, e analisar os processos históricos evolutivos que moldaram suas histórias demográficas. Para tal, foram estudados dois marcadores mitocondriais (região controle e o gene da Citocromo Oxidase I) e loci de microssatélites de indivíduos de 7 regiões do Caribe (Florida, Bahamas, Turks e Caicos, Porto Rico, Cuba, Colômbia e Venezuela) e 11 estados do Brasil (Pará, Maranhão, Piauí, Ceará, Rio Grande do Norte, Pernambuco, Alagoas, Bahia, Espírito Santo, Rio de Janeiro e São Paulo). Dentro de cada espécie foram observadas duas linhagens mitocondriais diferentes, que co-ocorriam, de maneira homogênea, ao longo de suas distribuições. Hipotetiso que essas linhagens foram formadas a partir de um evento de vicariância com contato secundário ou como consequência de um efeito gargalo seguido de expansão. As duas linhagens são evidentes nas sequências da região controle mitocondrial, mas no gene da COI foram evidentes apenas em P. cf. argus do Caribe. As linhagens do Brasil se separaram há aproximadamente 233 - 288 mil anos e cada uma sofreu expansão em tempos diferentes, a primeira se expandiu há 100 mil anos e a segunda linhagem há 50 mil anos. As linhagens do Caribe se separaram cerca de 1 milhão de anos atrás e possuem o mesmo tempo de expansão, 50 mil anos. Os microssatélites não revelaram subdivisão populacional para nenhuma das duas espécies, porém os marcadores, juntos, sugeriram um fluxo gênico diferenciado entre localidades expostas a diferentes correntes marítimas. Considerando que essas lagostas são intensamente explotadas, é importante ser cuidadoso no momento de definir estoques pesqueiros. Para a espécie do Brasil, dois estoques pesqueiros foram sugeridos, o primeiro do Pará à Bahia e o segundo do sul da Bahia a São Paulo. Para a espécie do Caribe, foi mantida e reforçada a hipótese de quatro estoques sugerida pela FAO (Norte, Sul, Centro-Norte e Centro-Sul).