Invertebrate predation and trophic cascades in a pelagic food web : The multiple roles of Chaoborus flavicans (Meigen) in a clay-turbid lake

In lake ecosystems, both fish and invertebrate predators have dramatic effects on their prey communities. Fish predation selects large cladocerans while invertebrate predators prefer prey of smaller size. Since invertebrate predators are the preferred food items for fish, their occurrence at high densities is often connected with the absence or low number of fish. It is generally believed that invertebrate predators can play a significant role only if the density of planktivorous fish is low. However, in eutrophic clay-turbid Lake Hiidenvesi (southern Finland), a dense population of predatory Chaoborus flavicans larvae coexists with an abundant fish population. The population covers the stratifying area of the lake and attains a maximum population density of 23000 ind. m-2. This thesis aims to clarify the effects of Chaoborus flavicans on the zooplankton community and the environmental factors facilitating the coexistence of fish and invertebrate predators. In the stratifying area of Lake Hiidenvesi, the seasonal succession of cladocerans was exceptional. The spring biomass peak of cladocerans was missing and the highest biomass occurred in midsummer. In early summer, the consumption rate by chaoborids clearly exceeded the production rate of cladocerans and each year the biomass peak of cladocerans coincided with the minimum chaoborid density. In contrast, consumption by fish was very low and each study year cladocerans attained maximum biomass simultaneously with the highest consumption by smelt (Osmerus eperlanus). The results indicated that Chaoborus flavicans was the main predator of cladocerans in the stratifying area of Lake Hiidenvesi. The clay turbidity strongly contributed to the coexistence of chaoborids and smelt at high densities. Turbidity exceeding 30 NTU combined with light intensity below 0.1 μE m-2 s-1provides an efficient daytime refuge for chaoborids, but turbidity alone is not an adequate refuge unless combined with low light intensity. In the non-stratifying shallow basins of Lake Hiidenvesi, light intensity exceeds this level during summer days at the bottom of the lake, preventing Chaoborus forming a dense population in the shallow parts of the lake. Chaoborus can be successful particularly in deep, clay-turbid lakes where they can remain high in the water column close to their epilimnetic prey. Suspended clay alters the trophic interactions by weakening the link between fish and Chaoborus, which in turn strengthens the effect of Chaoborus predation on crustacean zooplankton. Since food web management largely relies on manipulations of fish stocks and the cascading effects of such actions, the validity of the method in deep clay-turbid lakes may be questioned.

Structure and function of GDNF receptor alpha splice variants

The growth factors of the glial cell line-derived neurotrophic factor (GDNF) family consisting of GDNF, neurturin (NRTN), artemin (ARTN) and persephin (PSPN), are involved in the development, differentiation and maintenance of many types of neurons. They also have important functions outside the nervous system in the development of kidney, testis and thyroid gland. Each of these GFLs preferentially binds to one of the glycosylphosphatidylinositol (GPI)-anchored GDNF family receptors α (GFRα). GDNF binds to GFRα1, NRTN to GFRα2, ARTN to GFRα3 and PSPN to GFRα4. The GFLs in the complex with their cognate GFRα receptors all bind to and signal through the receptor tyrosine kinase RET. Alternative splicing of the mouse GFRα4 gene yields three splice isoforms. These had been described as putative GPI-anchored, transmembrane and soluble forms. My goal was to characterise the function of the different forms of mouse GFRα4. I firstly found that the putative GPI-anchored GFRα4 (GFRα4-GPI) is glycosylated, membrane-bound, GPI-anchored and interacts with PSPN and RET. We also showed that mouse GFRα4-GPI mediates PSPN-induced phosphorylation of RET, promotes PSPN-dependent neuronal differentiation of the rat pheochromocytoma cell line PC6-3 and PSPN-dependent survival of cerebellar granule neurons (CGN). However, although this receptor can mediate PSPN-signalling and activate RET, GFRα4-GPI does not recruit RET into lipid rafts. The recruitment of RET into lipid rafts has previously been thought to be a crucial event for GDNF- and GFL-mediated signalling via RET. I secondly demonstrated that the putative transmembrane GFRα4 (GFRα4-TM) is indeed a real transmembrane GFRα4 protein. Although it has a weak binding capacity for PSPN, it can not mediate PSPN-dependent phosphorylation of RET, neuronal differentiation or survival. These data show that GFRα4-TM is inactive as a receptor for PSPN. Surprisingly, GFRα4-TM can negatively regulate PSPN-mediated signalling via GFRα4-GPI. GFRα4-TM interacts with GFRα4-GPI and blocks PSPN-induced phosphorylation of RET, neuronal differentiation as well as survival. Taken together, our data show that GFRα4-TM may act as a dominant negative inhibitor of PSPN-mediated signaling. The most exciting part of my work was the finding that the putative soluble GFRα4 (GFRα4-sol) can form homodimers and function as an agonist of the RET receptor. In the absence of PSPN, GFRα4-sol can promote the phosphorylation of RET, trigger the activation of the PI-3K/AKT pathway, induce neuronal differentiation and support the survival of CGN. Our findings are in line with a recent publication showing the GFRα4-sol might contribute to the inherited cancer syndrome multiple endocrine neoplasia type 2. Our data provide an explanation to how GFRα4-sol may cause or modify the disease. Mammalian GFRα4 receptors all lack the first Cys-rich domain which is present in other GFRα receptors. In the final part of my work I have studied the function of this particular domain. I created a truncated GFRα1 construct lacking the first Cys-rich domain. Using binding assays in both cellular and cell-free systems, phosphorylation assays with RET, as well as neurite outgrowth assays, we found that the first Cys-rich domain contributes to an optimal function of GFRα1, by stabilizing the interaction between GDNF and GFRα1.

Nuclear Import Mechanisms of STAT and NF-kB Transcription Factors

The eukaryotic cell nucleoplasm is separated from the cytoplasm by the nuclear envelope. This compartmentation of eukaryotic cells requires that all nuclear proteins must be transported from the cytoplasm into the nucleus. Transport of macromolecules between the nucleus and the cytoplasm occurs through nuclear pore complexes (NPCs). Proteins to be targeted into the nucleus by the classical nuclear import system contain nuclear localization signals (NLSs), which are recognized by importin alpha, the NLS receptor. Importin alpha binds to importin beta, which docks the importin-cargo complex on the cytoplasmic side of the NPC and mediates the movement of the complex into the nucleus. Presently six human importin alpha isoforms have been identified. Transcription factors are among the most important regulators of gene expression in eukaryotic organisms. Transcription factors bind to specific DNA sequences on target genes and modulate the activity of the target gene. Many transcription factors, including signal transducers and activators of transcription (STAT) and nuclear factor kB (NF-kB), reside in the cytoplasm in an inactive form, and upon activation they are rapidly transported into the nucleus. In the nucleus STATs and NF-kB regulate the activity of genes whose products are critical in controlling numerous cellular and organismal processes, such as inflammatory and immune responses, cell growth, differentiation and survival. The aim of this study was to investigate the nuclear import mechanisms of STAT and NF-kB transcription factors. This work shows that STAT1 homodimers and STAT1/STAT2 heterodimers bind specifically and directly to importin alpha5 molecule via unconventional dimer-specific NLSs. Importin alpha molecules have two regions, which have been shown to directly interact with the amino acids in the NLS of the cargo molecule. The Arm repeats 2-4 comprise the N-terminal NLS binding site and Arm repeats 7-8 the C-terminal NLS binding site. In this work it is shown that the binding site for STAT1 homodimers and STAT1/STAT2 heterodimers is composed of Arm repeats 8 and 9 of importin alpha5 molecule. This work demonstrates that all NF-kB proteins are transported into the nucleus by importin alpha molecules. In addition, NLS was identified in RelB protein. The interactions between NF-kB proteins and importin alpha molecules were found to be directly mediated by the NLSs of NF-kB proteins. Moreover, we found that p50 binds to the N-terminal and p65 to the C-terminal NLS binding site of importin alpha3. The results from this thesis work identify previously uncharacterized mechanisms in nuclear import of STAT and NF-kB. These findings provide new insights into the molecular mechanisms regulating the signalling cascades of these important transcription factors from the cytoplasm into the nucleus to the target genes.

Lipid-Containing Icosahedral dsDNA Bacteriophages: Entry, Exit and Structure

In this thesis three icosahedral lipid-containing double-stranded (ds) deoxyribonucleic acid (DNA) bacteriophages have been studied: PRD1, Bam35 and P23-77. The work focuses on the entry, exit and structure of the viruses. PRD1 is the type member of the Tectiviridae family, infecting a variety of Gram-negative bacteria. The PRD1 receptor binding complex, consisting of the penton protein P31, the spike protein P5 and the receptor binding protein P2 recognizes a specific receptor on the host surface. In this study we found that the transmembrane protein P16 has an important stabilization function as the fourth member of the receptor binding complex and protein P16 may have a role in the formation of a tubular membrane structure, which is needed in the ejection of the genome into the cell. Phage Bam35 (Tectiviridae), which infects Gram-positive hosts, has been earlier found to resemble PRD1 in morphology and genome organization The uncharacterized early and late events in the Bam35 life cycle were studied by electrochemical methods. Physiological changes in the beginning of the infection were found to be similar in both lysogenic and nonlysogenic cell lines, Bam35 inducing a temporal decrease of membrane voltage and K+ efflux. At the end of the infection cycle physiological changes were observed only in the nonlysogenic cell line. The strong K+ efflux 40 min after infection and the induced premature cell lysis propose that Bam35 has a similar holin-endolysin lysis system to that of PRD1. Thermophilic icosahedral dsDNA Thermus phages P23-65H, P23-72 and P23-77 have been proposed to belong to the Tectiviridae family. In this study these phages were compared to each other. Analysis of structural protein patterns and stability revealed these phages to be very similar but not identical. The most stable of the studied viruses, P23-77, was further analyzed in more detail. Cryo-electron microscopy and three-dimensional image reconstruction was used to determine the structure of virus to 14 Å resolution. Results of thin layer chromatography for neutral lipids together with analysis of the three dimensional reconstruction of P23-77 virus particle revealed the presence of an internal lipid membrane. The overall capsid architecture of P23-77 is similar to PRD1 and Bam35, but most closely it resembles the structure of the capsid of archaeal virus SH1. This complicates the classification of dsDNA, internal lipid-containing icosahedral viruses.

Molecular genetics of RECQL4 syndromes

RAPADILINO syndrome is an autosomally resessively inherited condition that belongs to a group of rare syndromes more common in Finland than in other parts of the world. RAPADILINO is characterized by pre- and postnatal growth retardation, radial ray defects, diarrhoea of unknown aetiology during chilhood, a facial resemblance with other patients and normal intelligence. In Finland, 15 patients with this condition have been found which compares with only five patients in other parts of the world. We found RECQL4 gene mutations in RAPADILINO patients and proved this syndrome to be allelic with a subgroup of Rothmund-Thomson syndrome (RTS). Later we found RECQL4 mutations in patients with Baller-Gerold syndrome (BGS). These three syndromes share clinical findings and differential diagnostics rely on poikiloderma and craniosynostosis not seen in RAPADILINO syndrome. We found five different mutations in the Finnish RAPADILINO patients. The g.2545delT mutation is the founder mutation in the Finnish population as all the patients are either homozygotes or compound heterozygotes for it. This mutation leads to the inframe skipping of exon seven from mRNA. The protein encoded by this mutant mRNA lacks the nuclear retention signal and thus leads to the mislocalization of the mutant protein. The genotype-phenotype correlation is not straightforward but it seems that RAPADILINO could be due to alteration in protein function and truncating mutations in both alleles are more common among RTS patients. RTS patients with RECQL4 mutations have an elevated risk for osteosarcoma, but their risk to develop other types of malignancies is not increased.Two Finnish RAPADILINO patients have been diagnosed with osteosarcoma, but in addition to this we have found an excess of lymphoma cases among the Finnish RAPADILINO patients. This difference between cancer types could be due to different mutations found in these syndromes. The mutation screening of the patients will help to differentiate patients who have RECQL4 mutations and thus the elevated cancer risk. Patients will benefit from the follow up since early detection of malignancies is important for the treatment.

Expression, Enzymatic Activities and Subcellular Localization of Hepatitis E Virus and Semliki Forest Virus Replicase Proteins

Viruses are submicroscopic, infectious agents that are obligate intracellular parasites. They adopt various types of strategies for their parasitic replication and proliferation in infected cells. The nucleic acid genome of a virus contains information that redirects molecular machinery of the cell to the replication and production of new virions. Viruses that replicate in the cytoplasm and are unable to use the nuclear transcription machinery of the host cell have developed their own transcription and capping systems. This thesis describes replication strategies of two distantly related viruses, hepatitis E virus (HEV) and Semliki Forest virus (SFV), which belong to the alphavirus-like superfamily of positive-strand RNA viruses. We have demonstrated that HEV and SFV share a unique cap formation pathway specific for alphavirus-like superfamily. The capping enzyme first acts as a methyltransferase, catalyzing the transfer of a methyl group from S-adenosylmethionine to GTP to yield m7GTP. It then transfers the methylated guanosine to the end of viral mRNA. Both reactions are virus-specific and differ from those described for the host cell. Therefore, these capping reactions offer attractive targets for the development of antiviral drugs. Additionally, it has been shown that replication of SFV and HEV takes place in association with cellular membranes. The origin of these membranes and the intracellular localization of the components of the replication complex were studied by modern microscopy techniques. It was demonstrated that SFV replicates in cytoplasmic membranes that are derived from endosomes and lysosomes. According to our studies, site for HEV replication seems to be the intermediate compartment which mediates the traffic between endoplasmic reticulum and the Golgi complex. As a result of this work, a unique mechanism of cap formation for hepatitis E virus replicase has been characterized. It represents a novel target for the development of specific inhibitors against viral replication.

Functional characterization of MutL homologue mismatch repair proteins and their variants

Mismatch repair (MMR) mechanisms repair DNA damage occurring during replication and recombination. To date, five human MMR genes, MSH2, MHS6, MSH3, MLH1 and PMS2 are known to be involved in the MMR function. Human MMR proteins form 3 different heterodimers: MutSα (MSH2 and MSH6) and MutSβ (MSH2 and MSH3), which are needed for mismatch recognition and binding, and MutLα (MLH1 and PMS2), which is needed for mediating interactions between MutS homologues and other MMR proteins. The other two MutL homologues, MLH3 and PMS1, have been shown to heterodimerize with MLH1. However, the heterodimers MutLγ (MLH1and MLH3) and MutLβ (MLH1 and PMS1) are able to correct mismatches only with low or no efficiency, respectively. A deficient MMR mechanism is associated with the hereditary colorectal cancer syndrome called hereditary nonpolyposis colorectal cancer (HNPCC) or Lynch syndrome. HNPCC is the most common hereditary colorectal cancer syndrome and accounts for 2-5% of all colorectal cancer cases. HNPCC-associated mutations have been found in 5 MMR genes: MLH1, MSH2, MSH6, PMS2 and MLH3. Most of the mutations have been found in MLH1 and MSH2 (~90%) and are associated with typical HNPCC, while mutations in MSH6, PMS2 and MLH3 are mainly linked to putative HNPCC families lacking the characteristics of the syndrome. More data of MLH3 mutations are needed to assess the significance of its mutations in HNPCC. In this study, were functionally characterized 51 nontruncating mutations in the MLH1, MLH3 and MSH2 genes to address their pathogenic significance and mechanism of pathogenicity. Of the 36 MLH1 mutations, 22 were deficient in more than one assay, 2 variants were impaired only in one assay, and 12 variants behaved like the wild type protein, whereas all seven MLH3 mutants functioned like the wild type protein in the assays. To further clarify the role and relevance of MLH3 in MMR, we analyzed the subcellular localization of the native MutL homologue proteins. Our immunofluorescence analyses indicated that when all the three MutL homologues are natively expressed in human cells, endogenous MLH1 and PMS2 localize in the nucleus, whereas MLH3 stays in the cytoplasm. The coexpression of MLH3 with MLH1 results in its partial nuclear localization. Only one MSH2 mutation was pathogenic in the in vitro MMR assay. Our study on MLH1 mutations could clearly distinguish nontruncating alterations with severe functional defects from those not or only slightly impaired in protein function. However, our study on MLH3 mutations suggest that MLH3 mutations per se are not sufficient to trigger MMR deficiency and the continuous nuclear localization of MLH1 and PMS2 suggest that MutLα has a major activity in MMR in vivo. Together with our functional assays, this confirms that MutLγ is a less efficient MMR complex than MutLα.

Effects of large-scale human land use on Capercaillie (Tetrao urogallus L.) populations in Finland

The Capercaillie (Tetrao urogallus L.) is often used as a focal species for landscape ecological studies: the minimum size for its lekking area is 300 ha, and the annual home range for an individual may cover 30 80 km2. In Finland, Capercaillie populations have decreased by approximately 40 85%, with the declines likely to have started in the 1940s. Although the declines have partly stabilized from the 1990s onwards, it is obvious that the negative population trend was at least partly caused by changes in human land use. The aim of this thesis was to study the connections between human land use and Capercaillie populations in Finland, using several spatial and temporal scales. First, the effect of forest age structure on Capercaillie population trends was studied in 18 forestry board districts in Finland, during 1965 1988. Second, the abundances of Capercaillie and Moose (Alces alces L.) were compared in terms of several land-use variables on a scale of 50 × 50 km grids and in five regions in Finland. Third, the effects of forest cover and fine-grain forest fragmentation on Capercaillie lekking area persistence were studied in three study locations in Finland, on 1000 and 3000 m spatial scales surrounding the leks. The analyses considering lekking areas were performed with two definitions for forest: > 60 and > 152 m3ha 1 of timber volume. The results show that patterns and processes at large spatial scales strongly influence Capercaillie in Finland. In particular, in southwestern and eastern Finland, high forest cover and low human impact were found to be beneficial for this species. Forest cover (> 60 m3ha 1 of timber) surrounding the lekking sites positively affected lekking area persistence only at the larger landscape scale (3000 m radius). The effects of older forest classes were hard to assess due to scarcity of older forests in several study areas. Young and middle-aged forest classes were common in the vicinity of areas with high Capercaillie abundances especially in northern Finland. The increase in the amount of younger forest classes did not provide a good explanation for Capercaillie population decline in 1965 1988. In addition, there was no significant connection between mature forests (> 152 m3ha 1 of timber) and lekking area persistence in Finland. It seems that in present-day Finnish landscapes, area covered with old forest is either too scarce to efficiently explain the abundance of Capercaillie and the persistence of the lekking areas, or the effect of forest age is only important when considering smaller spatial scales than the ones studied in this thesis. In conclusion, larger spatial scales should be considered for assessing the future Capercaillie management. According to the proposed multi-level planning, the first priority should be to secure the large, regional-scale forest cover, and the second priority should be to maintain fine-grained, heterogeneous structure within the separate forest patches. A management unit covering hundreds of hectares, or even tens or hundreds of square kilometers, should be covered, which requires regional-level land-use planning and co-operation between forest owners.

The viral coat protein is regulated by HSP70 and HSP40 in Potato virus A infection

Plus-stranded (plus) RNA viruses multiply within a cellular environment as tightly integrated units and rely on the genetic information carried within their genomes for multiplication and, hence, persistence. The minimal genomes of plus RNA viruses are unable to encode the molecular machineries that are required for virus multiplication. This sets requisites for the virus, which must form compatible interactions with host components during multiplication to successfully utilize primary metabolites as building blocks or metabolic energy, and to divert the protein synthesis machinery for production of viral proteins. In fact, the emerging picture of a virus-infected cell displays tight integration with the virus, from simple host and virus protein interactions through to major changes in the physiological state of the host cell. This study set out to develop a method for the identification of host components, mainly host proteins, that interact with proteins of Potato virus A (PVA; Potyvirus) during infection. This goal was approached by developing affinity-tag based methods for the purification of viral proteins complexed with associated host proteins from infected plants. Using this method, host membrane-associated viral ribonucleoprotein (RNP) complexes were obtained, and several host and viral proteins could be identified as components of these complexes. One of the host proteins identified using this strategy was a member of the heat shock protein 70 (HSP70) family, and this protein was chosen for further analysis. To enable the analysis of viral gene expression, a second method was developed based on Agrobacterium-mediated virus genome delivery into plant cells, and detection of virally expressed Renilla luciferase (RLUC) as a quantitative measure of viral gene expression. Using this method, it was observed that down-regulation of HSP70 caused a PVA coat protein (CP)-mediated defect associated with replication. Further experimentation suggested that CP can inhibit viral gene expression and that a distinct translational activity coupled to replication, referred to as replication-associated translation (RAT), exists. Unlike translation of replication-deficient viral RNA, RAT was dependent on HSP70 and its co-chaperone CPIP. HSP70 and CPIP together regulated CP turnover by promoting its modification by ubiquitin. Based on these results, an HSP70 and CPIP-driven mechanism that functions to regulate CP during viral RNA replication and/or translation is proposed, possibly to prevent premature particle assembly caused by CP association with viral RNA.

Phospholipids of lipid-containing bacteriophages and their transbilayer distribution

In this study we used electro-spray ionization mass-spectrometry to determine phospholipid class and molecular species compositions in bacteriophages PM2, PRD1, Bam35 and phi6 as well as their hosts. To obtain compositional data of the individual leaflets, phospholipid transbilayer distribution in the viral membranes was studied. We found that 1) the membranes of all studied bacteriophage are enriched in PG as compared to the host membranes, 2) molecular species compositions in the phage and host membranes are similar, and 3) phospholipids in the viral membranes are distributed asymmetrically with phosphatidylglycerol enriched in the outer leaflet and phosphatidylethanolamine in the inner one (except Bam35). Alternative models for selective incorporation of phospholipids to phages and for the origins of the asymmetric phospholipid transbilayer distribution are discussed. Notably, the present data are also useful when constructing high resolution structural models of bacteriophages, since diffraction methods cannot provide a detailed structure of the membrane due to high motility of the lipids and lack of symmetric organization of membrane proteins.

Export of organic matter, sulphate and base cations from boreal headwater catchments downstream to the coast: impacts of land use and climate

The terrestrial export of dissolved organic matter (DOM) is associated with climate, vegetation and land use, and thus is under the influence of climatic variability and human interference with terrestrial ecosystems, their soils and hydrological cycles. The present study provides an assessment of spatial variation of DOM concentrations and export, and interactions between DOM, catchment characteristics, land use and climatic factors in boreal catchments. The influence of catchment characteristics, land use and climatic drivers on the concentrations and export of total organic carbon (TOC), total organic nitrogen (TON) and dissolved organic phosphorus (DOP) was estimated using stream water quality, forest inventory and climatic data from 42 Finnish pristine forested headwater catchments, and water quality monitoring, GIS land use, forest inventory and climatic data from the 36 main Finnish rivers (and their sub-catchments) flowing to the Baltic Sea. Moreover, the export of DOM in relation to land use along a European climatic gradient was studied using river water quality and land use data from four European areas. Additionally, the role of organic and minerogenic acidity in controlling pH levels in Finnish rivers and pristine streams was studied by measuring organic anion, sulphate (SO4) and base cation (Ca, Mg, K and Na) concentrations. In all study catchments, TOC was a major fraction of DOM, with much lower proportions of TON and DOP. Moreover, most of TOC and TON was in a dissolved form. The correlation between TOC and TON concentrations was strong and TOC concentrations explained 78% of the variation in TON concentrations in pristine headwater streams. In a subgroup of 20 headwater catchments with similar climatic conditions and low N deposition in eastern Finland, the proportion of peatlands in the catchment and the proportion of Norway spruce (Picea abies Karsten) of the tree stand had the strongest correlation with the TOC and TON concentrations and export. In Finnish river basins, TOC export increased with the increasing proportion of peatland in the catchment, whereas TON export increased with increasing extent of agricultural land. The highest DOP concentrations and export were recorded in river basins with a high extent of agricultural land and urban areas, reflecting the influence of human impact on DOP loads. However, the most important predictor for TOC, TON and DOP export in Finnish rivers was the proportion of upstream lakes in the catchment. The higher the upstream lake percentage, the lower the export indicating organic matter retention in lakes. Molar TOC:TON ratio decreased from headwater catchments covered by forests and peatlands to the large river basins with mixed land use, emphasising the effect of the land use gradient on the stoichiometry of rivers. This study also demonstrated that the land use of the catchments is related to both organic and minerogenic acidity in rivers and pristine headwater streams. Organic anion dominated in rivers and streams situated in northern Finland, reflecting the higher extent of peatlands in these areas, whereas SO4 dominated in southern Finland and on western coastal areas, where the extent of fertile areas, agricultural land, urban areas, acid sulphate soils, and sulphate deposition is highest. High TOC concentrations decreased pH values in the stream and river water, whereas no correlation between SO4 concentrations and pH was observed. This underlines the importance of organic acids in controlling pH levels in Finnish pristine headwater streams and main rivers. High SO4 concentrations were associated with high base cation concentrations and fertile areas, which buffered the effects of SO4 on pH.