Fluorescence in situ hybridization and spectral imaging of coral-associated bacterial communities

Microbial communities play important roles in the functioning of coral reef communities. However, extensive autofluorescence of coral tissues and endosymbionts limits the application of standard fluorescence in situ hybridization (FISH) techniques for the identification of the coral-associated bacterial communities. This study overcomes these limitations by combining FISH and spectral imaging.

Exploring hydrocarbonoclastic bacterial communities in the estuarine surface microlayer

Bacteria that degrade polycyclic aromatic hydrocarbons (PAHs) in the estuarine surface microlayer (SML) of the Ria de Aveiro, Portugal—which is chronically polluted with oil hydrocarbons (OH)—were isolated and characterized; Pseudomonas was dominant among the PAH-degrading bacteria. Screening for PAH dioxygenase genes detected almost identical nahAc genes (encoding the alpha subunits of naphthalene dioxygenase) in 2 phylogenetically distinct isolates: Pseudomonas sp. and an unknown species of the family Enterobacteriaceae; this suggested that horizontal transfer of nah genes might be involved in PAH degradation in the SML. We also investigated the effect of PAH contamination on the spatial variability of the bacterioneuston along a gradient of pollution in the estuarine system of the Ria de Aveiro. Culture-independent techniques—fluorescence in situ hy - bridization (FISH) and denaturing-gradient gel electrophoresis (DGGE)—revealed a similar structure among the bacterioneuston communities along the estuary. In contrast, we detected differences in the relative abundance and diversity of organisms of the Gammaproteobacteria, including those of the genus Pseudomonas (which belongs to the Gammaproteobacteria). This is the first insight into the hydrocarbonoclastic bacterial communities in the SML of an estuarine area polluted with hydrocarbons. Our findings highlight the importance of SML-adapted hydrocarbonoclastic bacterioneuston as a potential source of new PAH-degrading bacteria (including new pseudomonads) with potential use in the bioremediation of hydrocarbon-polluted ecosystems.

Origin and ecological selection of core and food-specific bacterial communities associated with meat and seafood spoilage

The microbial spoilage of meat and seafood products with short shelf lives is responsible for a significant amount of food waste. Food spoilage is a very heterogeneous process, involving the growth of various, poorly characterized bacterial communities. In this study, we conducted 16S ribosomal RNA gene pyrosequencing on 160 samples of fresh and spoiled foods to comparatively explore the bacterial communities associated with four meat products and four seafood products that are among the most consumed food items in Europe. We show that fresh products are contaminated in part by a microbiota similar to that found on the skin and in the gut of animals. However, this animal-derived microbiota was less prevalent and less abundant than a core microbiota, psychrotrophic in nature, mainly originated from the environment (water reservoirs). We clearly show that this core community found on meat and seafood products is the main reservoir of spoilage bacteria. We also show that storage conditions exert strong selective pressure on the initial microbiota: alpha diversity in fresh samples was 189 +/- 58 operational taxonomic units (OTUs) but dropped to 27 +/- 12 OTUs in spoiled samples. The OTU assemblage associated with spoilage was shaped by low storage temperatures, packaging and the nutritional value of the food matrix itself. These factors presumably act in tandem without any hierarchical pattern. Most notably, we were also able to identify putative new clades of dominant, previously undescribed bacteria occurring on spoiled seafood, a finding that emphasizes the importance of using culture-independent methods when studying food microbiota.

High diversity of skin-associated bacterial communities of marine fishes is promoted by their high variability among body parts, individuals and species

Animal-associated microbiotas form complex communities, which are suspected to play crucial functions for their host fitness. However, the biodiversity of these communities, including their differences between host species and individuals, has been scarcely studied, especially in case of skin-associated communities. In addition, the intraindividual variability (i.e. between body parts) has never been assessed to date. The objective of this study was to characterize skin bacterial communities of two teleostean fish species, namely the European seabass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata), using a high-throughput DNA sequencing method. In order to focus on intrinsic factors of host-associated bacterial community variability, individuals of the two species were raised in controlled conditions. Bacterial diversity was assessed using a set of four complementary indices, describing the taxonomic and phylogenetic facets of biodiversity and their respective composition (based on presence/absence data) and structure (based on species relative abundances) components. Variability of bacterial diversity was quantified at the interspecific, interindividual and intraindividual scales. We demonstrated that fish surfaces host highly diverse bacterial communities, whose composition was very different from that of surrounding bacterioplankton. This high total biodiversity of skin-associated communities was supported by the important variability, between host species, individuals and the different body parts (dorsal, anal, pectoral and caudal fins).

Bacterial communities of the specialty crop phyllosphere: response to biological soil amendment use, rainfall, and insect visitation

Microorganisms in the plant rhizosphere, the zone under the influence of roots, and phyllosphere, the aboveground plant habitat, exert a strong influence on plant growth, health, and protection. Tomatoes and cucumbers are important players in produce safety, and the microbial life on their surfaces may contribute to their fitness as hosts for foodborne pathogens such as Salmonella enterica and Listeria monocytogenes. External factors such as agricultural inputs and environmental conditions likely also play a major role. However, the relative contributions of the various factors at play concerning the plant surface microbiome remain obscure, although this knowledge could be applied to crop protection from plant and human pathogens. Recent advances in genomic technology have made investigations into the diversity and structure of microbial communities possible in many systems and at multiple scales. Using Illumina sequencing to profile particular regions of the 16S rRNA gene, this study investigates the influences of climate and crop management practices on the field-grown tomato and cucumber microbiome. The first research chapter (Chapter 3) involved application of 4 different soil amendments to a tomato field and profiling of harvest-time phyllosphere and rhizosphere microbial communities. Factors such as water activity, soil texture, and field location influenced microbial community structure more than soil amendment use, indicating that field conditions may exert more influence on the tomato microbiome than certain agricultural inputs. In Chapter 4, the impact of rain on tomato and cucumber-associated microbial community structures was evaluated. Shifts in bacterial community composition and structure were recorded immediately following rain events, an effect which was partially reversed after 4 days and was strongest on cucumber fruit surfaces. Chapter 5 focused on the contribution of insect visitors to the tomato microbiota, finding that insects introduced diverse bacterial taxa to the blossom and green tomato fruit microbiome. This study advances our understanding of the factors that influence the microbiomes of tomato and cucumber. Farms are complex environments, and untangling the interactions between farming practices, the environment, and microbial diversity will help us develop a comprehensive understanding of how microbial life, including foodborne pathogens, may be influenced by agricultural conditions.

Using the Gini coefficient with BIOLOG substrate utilisation data to provide an alternative quantitative measure for comparing bacterial soil communities

A measure quantifying unequal use of carbon sources, the Gini coefficient (G), has been developed to allow comparisons of the observed functional diversity of bacterial soil communities. This approach was applied to the analysis of substrate utilisation data obtained from using BIOLOG microtiter plates in a study which compared decomposition processes in two contrasting plant substrates in two different soils. The relevance of applying the Gini coefficient as a measure of observed functional diversity, for soil bacterial communities is evaluated against the Shannon index (H) and average well colour development (AWCD), a measure of the total microbial activity. Correlation analysis and analysis of variance of the experimental data show that the Gini coefficient, the Shannon index and AWCD provided similar information when used in isolation. However, analyses based on the Gini coefficient and the Shannon index, when total activity on the microtiter plates was maintained constant (i.e. AWCD as a covariate), indicate that additional information about the distribution of carbon sources being utilised can be obtained. We demonstrate that the Lorenz curve and its measure of inequality, the Gini coefficient, provides not only comparable information to AWCD and the Shannon index but when used together with AWCD encompasses measures of total microbial activity and absorbance inequality across all the carbon sources. This information is especially relevant for comparing the observed functional diversity of soil microbial communities.

Vertical distribution of bacterial and archaeal communities along discrete layers of a deep-sea cold sediment sample at the East Pacific Rise (similar to 13 degrees N)

The community structure and vertical distribution of prokaryotes in a deep-sea (ca. 3,191 m) cold sediment sample (ca. 43 cm long) collected at the East Pacific Rise (EPR) similar to 13 degrees N were studied with 16SrDNA-based molecular analyses. Total community DNA was extracted from each of four discrete layers EPRDS-1, -2, -3 and -4 (from top to bottom) and 16S rDNA were amplified by PCR. Cluster analysis of DGGE profiles revealed that the bacterial communities shifted sharply between EPRDS-1 and EPRDS-2 in similarity coefficient at merely 49%. Twenty-three sequences retrieved from DGGE bands fell into 11 groups based on BLAST and bootstrap analysis. The dominant groups in the bacterial communities were Chloroflexi, Gamma proteobacteria, Actinobacterium and unidentified bacteria, with their corresponding percentages varying along discrete layers. Pairwise Fst (F-statistics) values between the archaeal clone libraries indicated that the archaeal communities changed distinctly between EPRDS-2 and EPRDS-3. Sequences from the archaeal libraries were divided to eight groups. Crenarchaea Marine Group I (MGI) was prevalent in EPRDS-1 at 83%, while Uncultured Crenarchaea group II B (UCII B) abounded in EPRDS-4 at 61%. Our results revealed that the vertically stratified distribution of prokaryotic communities might be in response to the geochemical settings and suggested that the sampling area was influenced by hydrothermalism. The copresence of members related to hydrothermalism and cold deep-sea environments in the microbial community indicated that the area might be a transitional region from hydrothermal vents to cold deep-sea sediments.

Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe

Soil microorganisms are critical to ecosystem functioning and the maintenance of soil fertility. However, despite global increases in the inputs of nitrogen (N) and phosphorus (P) to ecosystems due to human activities, we lack a predictive understanding of how microbial communities respond to elevated nutrient inputs across environmental gradients. Here we used high-throughput sequencing of marker genes to elucidate the responses of soil fungal, archaeal, and bacterial communities using an N and P addition experiment replicated at 25 globally distributed grassland sites. We also sequenced metagenomes from a subset of the sites to determine how the functional attributes of bacterial communities change in response to elevated nutrients. Despite strong compositional differences across sites, microbial communities shifted in a consistent manner with N or P additions, and the magnitude of these shifts was related to the magnitude of plant community responses to nutrient inputs. Mycorrhizal fungi and methanogenic archaea decreased in relative abundance with nutrient additions, as did the relative abundances of oligotrophic bacterial taxa. The metagenomic data provided additional evidence for this shift in bacterial life history strategies because nutrient additions decreased the average genome sizes of the bacterial community members and elicited changes in the relative abundances of representative functional genes. Our results suggest that elevated N and P inputs lead to predictable shifts in the taxonomic and functional traits of soil microbial communities, including increases in the relative abundances of faster-growing, copiotrophic bacterial taxa, with these shifts likely to impact belowground ecosystems worldwide.

Bacterial Community Structure in the Hindgut of Wild and Captive Dugongs (Dugong dugon)

Dugongs (Dugong dugon) are marine mammals that obtain nutrients through hindgut fermentation of seagrass, however, the microbes responsible have not been identified. This study used denaturing gradient gel electrophoresis (DGGE) and 454-pyrosequencing to profile hindgut bacterial communities in wild dugongs. Faecal samples obtained from 32 wild dugongs representing four size/maturity classes, and two captive dugongs fed on cos lettuce were screened using DGGE. Partial 16S rRNA gene profiles of hindgut bacteria from wild dugong calves and juveniles were grouped together and were different to those in subadults and adults. Marked differences between hindgut bacterial communities of wild and captive dugongs were also observed, except for a single captive whose profile resembled wild adults following an unsuccessful reintroduction to the wild. Pyrosequencing of hindgut communities in two wild dugongs confirmed the stability of bacterial populations, and Firmicutes (average 75.6% of Operational Taxonomic Units [OTUs]) and Bacteroidetes (19.9% of OTUs) dominated. Dominant genera were Roseburia, Clostridium, and Bacteroides. Hindgut microbial composition and diversity in wild dugongs is affected by ontogeny and probably diet. In captive dugongs, the absence of the dominant bacterial DNA bands identified in wild dugongs is probably dependent upon prevailing diet and other captive conditions such as the use of antibiotics. This study represents a first step in the characterisation of a novel microbial ecosystem-the marine hindgut of Sirenia.

Impact of the Human Bacterial Environment on Mycobacteriosis and Allergy

My work describes two sectors of the human bacterial environment: 1. The sources of exposure to infectious non-tuberculous mycobacteria. 2. Bacteria in dust, reflecting the airborne bacterial exposure in environments protecting from or predisposing to allergic disorders. Non-tuberculous mycobacteria (NTM) transmit to humans and animals from the environment. Infection by NTM in Finland has increased during the past decade beyond that by Mycobacterium tuberculosis. Among the farm animals, porcine mycobacteriosis is the predominant NTM disease in Finland. Symptoms of mycobacteriosis are found in 0.34 % of slaughtered pigs. Soil and drinking water are suspected as sources for humans and bedding materials for pigs. To achieve quantitative data on the sources of human and porcine NTM exposure, methods for quantitation of environmental NTM are needed. We developed a quantitative real-time PCR method, utilizing primers targeted at the 16S rRNA gene of the genus of Mycobacterium. With this method, I found in Finnish sphagnum peat, sandy soils and mud high contents of mycobacterial DNA, 106 to 107 genome equivalents per gram. A similar result was obtained by a method based on the Mycobacterium-specific hybridization of 16S rRNA. Since rRNA is found mainly in live cells, this result shows that the DNA detected by qPCR mainly represented live mycobacteria. Next, I investigated the occurrence of environmental mycobacteria in the bedding materials obtained from 5 pig farms with high prevalence (>4 %) of mycobacteriosis. When I used for quantification the same qPCR methods as for the soils, I found that piggery samples contained non-mycobacterial DNA that was amplified in spite of several mismatches with the primers. I therefore improved the qPCR assay by designing Mycobacterium-specific detection probes. Using the probe qPCR assay, I found 105 to 107 genome equivalents of mycobacterial DNA in unused bedding materials and up to 1000 fold more in the bedding collected after use in the piggery. This result shows that there was a source of mycobacteria in the bedding materials purchased by the piggery and that mycobacteria increased in the bedding materials during use in the piggery. Allergic diseases have reached epidemic proportions in urbanized countries. At the same time, childhood in rural environment or simple living conditions appears to protect against allergic disorders. Exposure to immunoreactive microbial components in rural environments seems to prevent allergies. I searched for differences in the bacterial communities of two indoor dusts, an urban house dust shown to possess immunoreactivity of the TH2-type and a farm barn dust with TH1-activity. The immunoreactivities of the dusts were revealed by my collaborators, in vitro in human dendritic cells and in vivo in mouse. The dusts accumulated >10 years in the respiratory zone (>1.5 m above floor), thus reflecting the long-term content of airborne bacteria at the two sites. I investigated these dusts by cloning and sequencing of bacterial 16S rRNA genes from dust contained DNA. From the TH2-active urban house dust, I isolated 139 16S rRNA gene clones. The most prevalent genera among the clones were Corynebacterium (5 species, 34 clones), Streptococcus (8 species, 33 clones), Staphylococcus (5 species, 9 clones) and Finegoldia (1 species, 9 clones). Almost all of these species are known as colonizers of the human skin and oral cavity. Species of Corynebacterium and Streptococcus have been reported to contain anti-inflammatory lipoarabinomannans and immunmoreactive beta-glucans respectively. Streptococcus mitis, found in the urban house dust is known as an inducer of TH2 polarized immunity, characteristic of allergic disorders. I isolated 152 DNA clones from the TH1-active farm barn dust and found species quite different from those found from the urban house dust. Among others, I found DNA clones representing Bacillus licheniformis, Acinetobacter lwoffii and Lactobacillus each of which was recently reported to possess anti-allergy immunoreactivity. Moreover, the farm barn dust contained dramatically higher bacterial diversity than the urban house dust. Exposure to this dust thus stimulated the human dendritic cells by multiple microbial components. Such stimulation was reported to promote TH1 immunity. The biodiversity in dust may thus be connected to its immunoreactivity. Furthermore, the bacterial biomass in the farm barn dust consisted of live intact bacteria mainly. In the urban house dust only ~1 % of the biomass appeared as intact bacteria, as judged by microscoping. Fragmented microbes may possess bioactivity different from that of intact cells. This was recently shown for moulds. If this is also valid for bacteria, the different immunoreactivities of the two dusts may be explained by the intactness of dustborne bacteria. Based on these results, we offer three factors potentially contributing to the polarized immunoreactivities of the two dusts: (i) the species-composition, (ii) the biodiversity and (iii) the intactness of the dustborne bacterial biomass. The risk of childhood atopic diseases is 4-fold lower in the Russian compared with the Finnish Karelia. This difference across the country border is not explainable by different geo-climatic factors or genetic susceptibilities of the two populations. Instead, the explanation must be lifestyle-related. It has already been reported that the microbiological quality of drinking water differs on the two sides of the borders. In collaboration with allergists, I investigated dusts collected from homes in the Russian Karelia and in the Finnish Karelia. I found that bacterial 16S rRNA genes cloned from the Russian Karelian dusts (10 homes, 234 clones) predominantly represented Gram-positive taxa (the phyla Actinobacteria and Firmicutes, 67%). The Russian Karelian dusts contained nine-fold more of muramic acid (60 to 70 ng mg-1) than the Finnish Karelian dusts (3 to 11 ng mg-1). Among the DNA clones isolated from the Finnish side (n=231), Gram-negative taxa (40%) outnumbered the Gram-positives (34%). Out of the 465 DNA clones isolated from the Karelian dusts, 242 were assigned to cultured validly described bacterial species. In Russian Karelia, animal-associated species e.g. Staphylococcus and Macrococcus were numerous (27 clones, 14 unique species). This finding may connect to the difference in the prevalence of allergy, as childhood contacts with pets and farm animals have been connected with low allergy risk. Plant-associated bacteria and plant-borne 16S rRNA genes (chloroplast) were frequent among the DNA clones isolated from the Finnish Karelia, indicating components originating from plants. In conclusion, my work revealed three major differences between the bacterial communtites in the Russian and in the Finnish Karelian homes: (i) the high prevalence of Gram-positive bacteria on the Russian side and of Gram-negative bacteria on the Finnish side and (ii) the rich presence of animal-associated bacteria on the Russian side whereas (iii) plant-associated bacteria prevailed on the Finnish side. One or several of these factors may connect to the differences in the prevalence of allergy.

Variation in the hindgut microbial communities of the Florida manatee, Trichechus manatus latirostris over winter in Crystal River, Florida

The Florida manatee, Trichechus manatus latirostris, is a hindgut-fermenting herbivore. In winter, manatees migrate to warm water overwintering sites where they undergo dietary shifts and may suffer from cold-induced stress. Given these seasonally induced changes in diet, the present study aimed to examine variation in the hindgut bacterial communities of wild manatees overwintering at Crystal River, west Florida. Faeces were sampled from 36 manatees of known sex and body size in early winter when manatees were newly arrived and then in mid-winter and late winter when diet had probably changed and environmental stress may have increased. Concentrations of faecal cortisol metabolite, an indicator of a stress response, were measured by enzyme immunoassay. Using 454-pyrosequencing, 2027 bacterial operational taxonomic units were identified in manatee faeces following amplicon pyrosequencing of the 16S rRNA gene V3/V4 region. Classified sequences were assigned to eight previously described bacterial phyla; only 0.36% of sequences could not be classified to phylum level. Five core phyla were identified in all samples. The majority (96.8%) of sequences were classified as Firmicutes (77.3 ± 11.1% of total sequences) or Bacteroidetes (19.5 ± 10.6%). Alpha-diversity measures trended towards higher diversity of hindgut microbiota in manatees in mid-winter compared to early and late winter. Beta-diversity measures, analysed through permanova, also indicated significant differences in bacterial communities based on the season.

Bayesian statistical analysis of bacterial diversity

Bacteria play an important role in many ecological systems. The molecular characterization of bacteria using either cultivation-dependent or cultivation-independent methods reveals the large scale of bacterial diversity in natural communities, and the vastness of subpopulations within a species or genus. Understanding how bacterial diversity varies across different environments and also within populations should provide insights into many important questions of bacterial evolution and population dynamics. This thesis presents novel statistical methods for analyzing bacterial diversity using widely employed molecular fingerprinting techniques. The first objective of this thesis was to develop Bayesian clustering models to identify bacterial population structures. Bacterial isolates were identified using multilous sequence typing (MLST), and Bayesian clustering models were used to explore the evolutionary relationships among isolates. Our method involves the inference of genetic population structures via an unsupervised clustering framework where the dependence between loci is represented using graphical models. The population dynamics that generate such a population stratification were investigated using a stochastic model, in which homologous recombination between subpopulations can be quantified within a gene flow network. The second part of the thesis focuses on cluster analysis of community compositional data produced by two different cultivation-independent analyses: terminal restriction fragment length polymorphism (T-RFLP) analysis, and fatty acid methyl ester (FAME) analysis. The cluster analysis aims to group bacterial communities that are similar in composition, which is an important step for understanding the overall influences of environmental and ecological perturbations on bacterial diversity. A common feature of T-RFLP and FAME data is zero-inflation, which indicates that the observation of a zero value is much more frequent than would be expected, for example, from a Poisson distribution in the discrete case, or a Gaussian distribution in the continuous case. We provided two strategies for modeling zero-inflation in the clustering framework, which were validated by both synthetic and empirical complex data sets. We show in the thesis that our model that takes into account dependencies between loci in MLST data can produce better clustering results than those methods which assume independent loci. Furthermore, computer algorithms that are efficient in analyzing large scale data were adopted for meeting the increasing computational need. Our method that detects homologous recombination in subpopulations may provide a theoretical criterion for defining bacterial species. The clustering of bacterial community data include T-RFLP and FAME provides an initial effort for discovering the evolutionary dynamics that structure and maintain bacterial diversity in the natural environment.

Microbial communities in Pb contaminated boreal forest soil

Lead contamination in the environment is of particular concern, as it is a known toxin. Until recently, however, much less attention has been given to the local contamination caused by activities at shooting ranges compared to large-scale industrial contamination. In Finland, more than 500 tons of Pb is produced each year for shotgun ammunition. The contaminant threatens various organisms, ground water and the health of human populations. However, the forest at shooting ranges usually shows no visible sign of stress compared to nearby clean environments. The aboveground biota normally reflects the belowground ecosystem. Thus, the soil microbial communities appear to bear strong resistance to contamination, despite the influence of lead. The studies forming this thesis investigated a shooting range site at Hälvälä in Southern Finland, which is heavily contaminated by lead pellets. Previously it was experimentally shown that the growth of grasses and degradation of litter are retarded. Measurements of acute toxicity of the contaminated soil or soil extracts gave conflicting results, as enchytraeid worms used as toxicity reporters were strongly affected, while reporter bacteria showed no or very minor decreases in viability. Measurements using sensitive inducible luminescent reporter bacteria suggested that the bioavailability of lead in the soil is indeed low, and this notion was supported by the very low water extractability of the lead. Nevertheless, the frequency of lead-resistant cultivable bacteria was elevated based on the isolation of cultivable strains. The bacterial and fungal diversity in heavily lead contaminated shooting sectors were compared with those of pristine sections of the shooting range area. The bacterial 16S rRNA gene and fungal ITS rRNA gene were amplified, cloned and sequenced using total DNA extracted from the soil humus layer as the template. Altogether, 917 sequenced bacterial clones and 649 sequenced fungal clones revealed a high soil microbial diversity. No effect of lead contamination was found on bacterial richness or diversity, while fungal richness and diversity significantly differed between lead contaminated and clean control areas. However, even in the case of fungi, genera that were deemed sensitive were not totally absent from the contaminated area: only their relative frequency was significantly reduced. Some operational taxonomic units (OTUs) assigned to Basidiomycota were clearly affected, and were much rarer in the lead contaminated areas. The studies of this thesis surveyed EcM sporocarps, analyzed morphotyped EcM root tips by direct sequencing, and 454-pyrosequenced fungal communities in in-growth bags. A total of 32 EcM fungi that formed conspicuous sporocarps, 27 EcM fungal OTUs from 294 root tips, and 116 EcM fungal OTUs from a total of 8 194 ITS2 454 sequences were recorded. The ordination analyses by non-parametric multidimensional scaling (NMS) indicated that Pb enrichment induced a shift in the EcM community composition. This was visible as indicative trends in the sporocarp and root tip datasets, but explicitly clear in the communities observed in the in-growth bags. The compositional shift in the EcM community was mainly attributable to an increase in the frequencies of OTUs assigned to the genus Thelephora, and to a decrease in the OTUs assigned to Pseudotomentella, Suillus and Tylospora in Pb-contaminated areas when compared to the control. The enrichment of Thelephora in contaminated areas was also observed when examining the total fungal communities in soil using DNA cloning and sequencing technology. While the compositional shifts are clear, their functional consequences for the dominant trees or soil ecosystem remain undetermined. The results indicate that at the Hälvälä shooting range, lead influences the fungal communities but not the bacterial communities. The forest ecosystem shows apparent functional redundancy, since no significant effects were seen on forest trees. Recently, by means of 454 pyrosequencing , the amount of sequences in a single analysis run can be up to one million. It has been applied in microbial ecology studies to characterize microbial communities. The handling of sequence data with traditional programs is becoming difficult and exceedingly time consuming, and novel tools are needed to handle the vast amounts of data being generated. The field of microbial ecology has recently benefited from the availability of a number of tools for describing and comparing microbial communities using robust statistical methods. However, although these programs provide methods for rapid calculation, it has become necessary to make them more amenable to larger datasets and numbers of samples from pyrosequencing. As part of this thesis, a new program was developed, MuSSA (Multi-Sample Sequence Analyser), to handle sequence data from novel high-throughput sequencing approaches in microbial community analyses. The greatest advantage of the program is that large volumes of sequence data can be manipulated, and general OTU series with a frequency value can be calculated among a large number of samples.