Evolution and detection of cyanobacterial hepatotoxin synthetase genes

Mass occurrences (blooms) of cyanobacteria are common in aquatic environments worldwide. These blooms are often toxic, due to the presence of hepatotoxins or neurotoxins. The most common cyanobacterial toxins are hepatotoxins: microcystins and nodularins. In freshwaters, the main producers of microcystins are Microcystis, Anabaena, and Planktothrix. Nodularins are produced by strains of Nodularia spumigena in brackish waters. Toxic and nontoxic strains of cyanobacteria co-occur and cannot be differentiated by conventional microscopy. Molecular biological methods based on microcystin and nodularin synthetase genes enable detection of potentially hepatotoxic cyanobacteria. In the present study, molecular detection methods for hepatotoxin-producing cyanobacteria were developed, based on microcystin synthetase gene E (mcyE) and the orthologous nodularin synthetase gene F (ndaF) sequences. General primers were designed to amplify the mcyE/ndaF gene region from microcystin-producing Anabaena, Microcystis, Planktothrix, and Nostoc, and nodularin-producing Nodularia strains. The sequences were used for phylogenetic analyses to study how cyanobacterial mcy genes have evolved. The results showed that mcy genes and microcystin are very old and were already present in the ancestor of many modern cyanobacterial genera. The results also suggested that the sporadic distribution of biosynthetic genes in modern cyanobacteria is caused by repeated gene losses in the more derived lineages of cyanobacteria and not by horizontal gene transfer. Phylogenetic analysis also proposed that nda genes evolved from mcy genes. The frequency and composition of the microcystin producers in 70 lakes in Finland were studied by conventional polymerase chain reaction (PCR). Potential microcystin producers were detected in 84% of the lakes, using general mcyE primers, and in 91% of the lakes with the three genus-specific mcyE primers. Potential microcystin-producing Microcystis were detected in 70%, Planktothrix in 63%, and Anabaena in 37% of the lakes. The presence and co-occurrence of potential microcystin producers were more frequent in eutrophic lakes, where the total phosphorus concentration was high. The PCR results could also be associated with various environmental factors by correlation and regression analyses. In these analyses, the total nitrogen concentration and pH were both associated with the presence of multiple microcystin-producing genera and partly explained the probability of occurrence of mcyE genes. In general, the results showed that higher nutrient concentrations increased the occurrence of potential microcystin producers and the risk for toxic bloom formation. Genus-specific probe pairs for microcystin-producing Anabaena, Microcystis, Planktothrix, and Nostoc, and nodularin-producing Nodularia were designed to be used in a DNA-chip assay. The DNA-chip can be used to simultaneously detect all these potential microcystin/nodularin producers in environmental water samples. The probe pairs detected the mcyE/ndaF genes specifically and sensitively when tested with cyanobacterial strains. In addition, potential microcystin/nodularin producers were identified in lake and Baltic Sea samples by the DNA-chip almost as sensitively as by quantitative real-time PCR (qPCR), which was used to validate the DNA-chip results. Further improvement of the DNA-chip assay was achieved by optimization of the PCR, the first step in the assay. Analysis of the mcy and nda gene clusters from various hepatotoxin-producing cyanobacteria was rewarding; it revealed that the genes were ancient. In addition, new methods detecting all the main producers of hepatotoxins could be developed. Interestingly, potential microcystin-producing cyanobacterial strains of Microcystis, Planktothrix, and Anabaena, co-occurred especially in eutrophic and hypertrophic lakes. Protecting waters from eutrophication and restoration of lakes may thus decrease the prevalence of toxic cyanobacteria and the frequency of toxic blooms.

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.

The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Polymorphic low penetrance genes and breast cancer : The role of genes involved in metabolism of xenobiotics, estrogens and reactive oxygen species

Various endogenous and exogenous factors have been reported to increase the risk of breast cancer. Many of those are related to prolonged lifetime exposure to estrogens. Furthermore, a positive family history of breast cancer and certain benign breast diseases are known to increase the risk of breast cancer. The role of lifestyle factors, such as use of alcohol and smoking has been an area of intensive study. Alcohol has been found to increase the risk of breast cancer, whereas the role of smoking has remained obscure. A multitude of enzymes are involved in the metabolism of estrogens and xenobiotics including the carcinogens found in tobacco smoke. Many of the metabolic enzymes exhibit genetic polymorphisms that can lead to inter-individual differences in their abilities to modify hazardous substrates. Therefore, in presence of a given chemical exposure, one subgroup of women may be more susceptible to breast carcinogenesis, since they carry unfavourable forms of the polymorphic genes involved in the metabolism of the chemical. In this work, polymorphic genes encoding for cytochrome P450 (CYP) 1A1 and 1B1, N-acetyl transferase 2 (NAT2), sulfotransferase 1A1 (SULT1A1), manganese superoxide dismutase (MnSOD) and vitamin D receptor (VDR) were investigated in relation to breast cancer susceptibility in a Finnish population. CYP1A1, CYP1B1 and SULT1A1 are involved in the metabolism of both estrogens and xenobiotics, whereas NAT2 is involved only in the latter. MnSOD is an antioxidant enzyme protecting cells from oxidative damage. VDR, in turn, mediates the effects of the active form of vitamin D (1,25(OH)2D3, calcitriol) on maintenance of calcium homeostasis and it has anti-proliferative effects in many cancer cells. A 1.3-fold (95% CIs 1.01-1.73) increased risk of breast cancer was seen among women who carried the NAT2 slow acetylator genotype and a 1.5-fold (95% CI 1.1-2.0) risk was found in women with a MnSOD variant A allele containing genotypes compared to women with the NAT2 rapid acetylator genotype or to those with the MnSOD VV genotype, respectively. Instead, women with the VDR a allele containing genotypes were found to be at a decreased risk for breast cancer (OR 0.73; 95% CI 0.54-0.98) compared to women with the AA genotype. No significant overall associations were found between SULT1A1 or CYP genotypes and breast cancer risk, whereas a combination of the CYP1B1 432Val allele containing genotypes with the NAT2 slow acetylator genotypes posed a 1.5-fold (95% CI 1.03-2.24) increased risk. Moreover, NAT2 slow acetylator genotype was found to be confined to women with an advanced stage of breast cancer (stages III and IV). Further evidence for the association of xenobiotic metabolising genes with breast cancer risk was found when active smoking was taken into account. Women who smoked less than 10 cigarettes/day and carried at least one CYP1B1 432Val variant allele, were at 3.1-fold (95% CI 1.32-7.12) risk of breast cancer compared to women who smoked the same amount but did not carry the variant allele. Furthermore, the risk was significantly increased with increasing number of the CYP1B1 432Val alleles (p for trend 0.005). In addition, women who smoked less than 5 pack-years and carried the NAT2 slow acetylator genotype were at a 2.6-fold (95% CI 1.01-6.48) increased risk of breast cancer compared to women who smoked the same amount but carried the NAT2 rapid acetylator genotype. Furthermore, the combination of the CYP1B1 432Val allele and the NAT2 slow acetylator genotype increased the risk of breast cancer by 2.5-fold (95% CI 1.11-5.45) among ever smokers. Instead, the MnSOD A allele was found to be a risk factor among postmenopausal long-term smokers (>15 years of smoking) (OR 5.1; 95% CI 1.4-18.4) or among postmenopausal women who had smoked more than 10 cigarettes/day (OR 5.5; 95% CI 1.3-23.4) compared to women who had similar smoking habits but carried the MnSOD V/V genotype. Similarly, within subgroups of postmenopausal women who were using oral contraceptives, hormone replacement therapy or alcohol, women carrying the MnSOD A allele genotypes seemed to be at increased risk of breast cancer compared to women with the MnSOD V/V genotype. A positive family history of breast cancer and high parity were shown to be inversely associated with breast cancer risk among women carrying the VDR ApaI a allele or among premenopausal women carrying the SULT1A1*2 allele, respectively.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.

Plasmids and aromatic degradation in Sphingomonas for bioremediation : Aromatic ring cleavage genes in soil and rhizosphere

Microbial degradation pathways play a key role in the detoxification and the mineralization of polyaromatic hydrocarbons (PAHs), which are widespread pollutants in soil and constituents of petroleum hydrocarbons. In microbiology the aromatic degradation pathways are traditionally studied from single bacterial strains with capacity to degrade certain pollutant. In soil the degradation of aromatics is performed by a diverse community of micro-organisms. The aim of this thesis was to study biodegradation on different levels starting from a versatile aromatic degrader Sphingobium sp. HV3 and its megaplasmid, extending to revelation of diversity of key catabolic enzymes in the environment and finally studying birch rhizoremediation in PAH-polluted soil. To understand biodegradation of aromatics on bacterial species level, the aromatic degradation capacity of Sphingobium sp. HV3 and the role of the plasmid pSKY4, was studied. Toluene, m-xylene, biphenyl, fluorene, phenanthrene were detected as carbon and energy sources of the HV3 strain. Tn5 transposon mutagenesis linked the degradation capacity of toluene, m-xylene, biphenyl and naphthalene to the pSKY4 plasmid and qPCR expression analysis showed that plasmid extradiol dioxygenases genes (bphC and xylE) are inducted by phenanthrene, m-xylene and biphenyl whereas the 2,4-dichlorophenoxyacetic acid herbicide induced the chlorocatechol 1,2-dioxygenase gene (tfdC) from the ortho-pathway. A method to study upper meta-pathway extradiol dioxygenase gene diversity in soil was developed. The extradiol dioxygenases catalyse cleavage of the aromatic ring between a hydroxylated carbon and an adjacent non-hydroxylated carbon (meta-cleavage). A high diversity of extradiol dioxygenases were detected from polluted soils. The detected extradiol dioxygenases showed sequence similarity to known catabolic genes of Alpha-, Beta-, and Gammaproteobacteria. Five groups of extradiol dioxygenases contained sequences with no close homologues in the database, representing novel genes. In rhizoremediation experiment with birch (Betula pendula) treatment specific changes of extradiol dioxygenase communities were shown. PAH pollution changed the bulk soil extradiol dioxygenase community structure and birch rhizosphere contained a more diverse extradiol dioxygenase community than the bulk soil showing a rhizosphere effect. The degradation of pyrene in soil was enhanced with birch seedlings compared to soil without birch. The complete 280,923 kb nucleotide sequence of pSKY4 plasmid was determined. The open reading frames of pSKY4 were divided into putative conjugative transfer, aromatic degradation, replication/maintaining and transposition/integration function-encoding proteins. Aromatic degradation orfs shared high similarity to corresponding genes in pNL1, a plasmid from the deep subsurface strain Novosphingobium aromaticivorans F199. The plasmid backbones were considerably more divergent with lower similarity, which suggests that the aromatic pathway has functioned as a plasmid independent mobile genetic element. The functional diversity of microbial communities in soil is still largely unknown. Several novel clusters of extradiol dioxygenases representing catabolic bacteria, whose function, biodegradation pathways and phylogenetic position is not known were amplified with single primer pair from polluted soils. These extradiol dioxygenase communities were shown to change upon PAH pollution, which indicates that their hosts function in PAH biodegradation in soil. Although the degradation pathways of specific bacterial species are substantially better depicted than pathways in situ, the evolution of degradation pathways for the xenobiotic compounds is largely unknown. The pSKY4 plasmid contains aromatic degradation genes in putative mobile genetic element causing flexibility/instability to the pathway. The localisation of the aromatic biodegradation pathway in mobile genetic elements suggests that gene transfer and rearrangements are a competetive advantage for Sphingomonas bacteria in the environment.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.