Multi-element fingerprinting and high throughput sequencing identify multiple elements that affect fungal communities in Quercus macrocarpa foliage

Diverse fungal mutualists, pathogens and saprobes colonize plant leaves. These fungi face a complex environment, in which stochastic dispersal interplays with abiotic and biotic filters. However, identification of the specific factors that drive the community assembly seems unattainable. We mined two broad data sets and identified chemical elements, to which dominant molecular operational taxonomic units (OTUs) in the foliage of a native tree respond most extremely. While many associations could be identified, potential complicating issues emerged. Those were related to unevenly distributed OTU frequency data, a large number of potentially explanatory variables and the disproportionate effects of outlier observations.Key words: community assembly, environmental filter, fungi, heavy metal enrichment, nutrient enrichment, oak, Quercus macrocarpaHyperdiverse fungal communities inhabit the foliage of most plants^1,2 and these fungal communities have been reported for virtually every plant that has been examined.³ Baas-Becking hypothesis states that environment selects microbial communities from the abundant and possibly globally distributed propagule pools.⁴ Although the foliage-associated communities—like other microbial communities—are suspected to be sensitive to environmental drivers, determination of the mechanisms that control the assembly of these foliar communities has remained difficult and elusive. Some of the proposed mechanisms include distance limitations to propagule dispersal,^5–7 volume limitations to propagule loads,⁷ or limitations set by the environmental conditions either on the scale of the site of fungal colonization⁸ or more broadly on a landscape level.^6,9 The forces that may control the fungal community assembly are overlaid by additional biotic controls that include compatibilities between the fungi and host species^10,11 or genotypes^6,12 and the competitive or facilitative interactions among the component fungal genotypes.^6,10–13 Although a variety of potential controls for the foliage-associated fungal communities have been speculated, very little consensus exists on the relative importance of the different drivers. For example, while macronutrient and heavy metal enrichment may have an influence on the composition fungal communities¹⁴ and populations,¹⁵ relative importance of various chemical elements in the foliage remains yet to be investigated.To evaluate the use of multi-element fingerprinting data produced by Inductively Coupled Plasma Mass Spectrometry (ICP-MS) in combination with high throughput 454-pyrosequencing for determining influential chemical elements in structuring of the leaf-associated fungal communities, we mined a recent dataset¹⁶ that explored the effects of urbanization on the diversity and composition of the fungal communities associated with a native tree Quercus macrocarpa. From a total list of more than 700 non-singleton fungal OTUs, we selected fifty with highest overall frequency to provide an observationrich dataset for elemental effect assessment; these OTUs accounted for 84.5% of all sequences. Even so, many of these OTUs had a number of zero frequencies (Fig. 1), highlighting one of the difficulties in the use of environmental sequencing data. We omitted one OTU (OTU630 with a likely affinity to Trimmatostroma cordae [Mycosphaerellaceae]) that was strongly affected by the original land use design (urbanization; Wilcoxon rank sum test with a Bonferroni adjustment) and therefore unlikely to be representative for the present analyses of elemental drivers. This OTU was replaced with one with the next highest frequency. The frequencies of these 50 OTUs were investigated in the context of concentrations for 29 elements after the omission of five (Ag, Au, C, δ¹³C, δ¹⁵N) in the final analyses because of their strong association with the land use or the difficulty of finding a biological relevance. Of the remaining elements three (Fe, Cr and Ni) had pairwise correlations exceeding 0.98 between the three pairings; others showed no similar high correlations. To allow comparable evaluation across the broad array of elements, all concentrations were standardized to have a mean equal to zero and a standard deviation equal to one.Open in a separate windowFigure 1Rank-ordered distribution of observed frequencies for those OTU s whose frequency had an extreme slope when associated with the concentrations of one or more chemical elements in the mixed effects model. The asterisk denotes one extreme frequency for OTU 313 with a value 0.8636. Numbers in parentheses indicate the number of observations with a frequency equal to zero. The OTU s were assigned to approximate taxa using BLAST:²⁰ 425: Alternaria alternata (Pleosporaceae); 46: Phoma glomerata (Pleosporaceae); 686: Aureobasidium pullulans (Dothioraceae); 520: Davidiella tassiana (Davidiellaceae); 567: Cladosporioum tenuissimum (Davidiellaceae); 313 Oidium heveae (mitosporic Erysiphaceae); 586: Erysiphe hypogena (Erysiphaceae); 671: Mycosphaerella microsora (Mycosphaerellaceae); 555: Pestalotiopsis sp. (Amphisphaeriaceae); 607: Pleiochaeta setosa (incertae sedis).To rank elements according to their magnitude of association with the abundance of each OTU, a total of 1,450 models (50 OTUs times 29 elements) relating element concentration to OTU abundance were fit to the data. For each model, OTU frequency was the dependent variable, element concentration and time (a factor with three levels) were fixed effects, and—to account for the spatial arrangement of the experimental units—random effects associated with tree nested within site were included in the error structure. Time by element interactions were also investigated and tested using a likelihood ratio test. These mixed effect models were fit using R and the package lme4 (www.rproject.org).Statistical “tests of significance” that produce p-values can be sensitive to assumptions or outliers. Because of this and the fact that our analyses evaluated a total of 1,450 models, p-values themselves were not considered a reliable measure of importance when associating elements with OTU frequency. Instead, we emphasized metrics that highlight extraordinary findings rather than rely on tests of statistical significance. This approach facilitates finding few elements that have the strongest effect on OTU frequency. Note that the use of standardized element concentrations (above) provided slope coefficients that are comparable across all models. “Extreme slopes”, i.e., models where the OTU response to element concentration was strongest, were identified as those with estimated slope coefficients in the lower or upper 2.5 percentile, i.e., those farther than 1.8 standard deviations from the mean across all estimated slopes (Fig. 2). Using this approach, we identified a total of 69 models with extreme slopes (Open in a separate windowFigure 2Distribution of estimated slopes (i.e., the slope for element concentration) for a model relating OTU frequency to element concentration, time and a concentration by time interaction, including a tree-nested-within-site random effect. The mean across all 1,450 OTU s is approximately zero; the two vertical lines identify upper and lower 2.5 percentiles, beyond which the slopes were considered extreme (large black symbols). The horizontal line identifies the cut off maximum leverage (0.24), above which the slopes were considered to have observations with high leverage. Models with observations with a high leverage were tested for extreme slopes by refitting without those observations. Models are ranked from bottom to top in order of increasing leverage and the element for which the high-leverage observations and extreme slopes were recorded are identified on the right y-axis.

Table 1

Slopes identified as extreme in our analyses

Ecosystem services–A tool for sustainable management of human–environment systems. Case study Finnish Forest Lapland

The concept of ecosystem services (ESs) is a relatively new scientific methodology, offering a possible approach to the prevention of ecological problems caused by human action and to the resolution of conflicts arising from land-use questions. Since ESs were launched as a major conceptual tool in the Millennium Ecosystem Assessment (MA, 2005), interest in them has been increasing. Despite the scientific as well as economic and political enthusiasm for the ES approach, only few case studies have as yet been published. We studied the interface between ESs and landscape planning in Forest Lapland, in northern Finland. In the article, we present a methodology and various databases which can be used in applied research on ESs. We classify the ESs offered by various biotopes of the study area, and examine the effects of different land-use forms on the provision of ESs. On the basis of our results, we suggest possible uses of the European CORINE land cover database in case studies.     

Accumulation of DNA damage and reduced levels of nicotine adenine dinucleotide in the brains of Atm-deficient mice.

Ataxia-telangiectasia (A-T) is a human genetic disorder caused by mutational inactivation of the ATM gene. A-T patients display a pleiotropic phenotype, in which a major neurological feature is progressive ataxia due to degeneration of cerebellar Purkinje and granule neurons. Disruption of the mouse Atm locus creates a murine model of A-T that exhibits most of the clinical and cellular features of the human disease, but the neurological phenotype is barely expressed. We present evidence for the accumulation of DNA strand breaks in the brains of Atm(-/-), supporting the notion that ATM plays a major role in maintaining genomic stability. We also show a perturbation of the steady state levels of pyridine nucleotides. There is a significant decrease in both the reduced and the oxidized forms of NAD and in the total levels of NADP(T) and NADP(+) in the brains of Atm(-/-) mice. The changes in NAD(T), NADH, NAD(+), NADP(T), and NADP(+) were progressive and observed primarily in the cerebellum of 4-month-old Atm(-/-) mice. Higher rates of mitochondrial respiration were also recorded in 4-month-old Atm(-/-) cerebella. Taken together, our findings support the hypothesis that absence of functional ATM results in continuous stress, which may be an important cause of the degeneration of cerebellar neurons in A-T.     

The Arabidopsis thaliana cysteine-rich receptor-like kinases CRK6 and CRK7 protect against apoplastic oxidative stress

Receptor-like kinases are important regulators of many different processes in plants. Despite their large number only a few have been functionally characterized. One of the largest subgroups of receptor-like kinases in Arabidopsis is the cysteine-rich receptor like kinases (CRKs). High sequence similarity among the CRKs has been suggested as major cause for functional redundancy. The genomic localization of CRK genes in back-to-back repeats has made their characterization through mutant analysis unpractical. Expression profiling has linked the CRKs with reactive oxygen species, important signaling molecules in plants. Here we have investigated the role of two CRKs, CRK6 and CRK7, and analyzed their role in extracellular ROS signaling. CRK6 and CRK7 are active protein kinases with differential preference for divalent cations. Our results suggest that CRK7 is involved in mediating the responses to extracellular but not chloroplastic ROS production.     

A surplus no more? Variation in krill availability impacts reproductive rates of Antarctic baleen whales

The krill surplus hypothesis of unlimited prey resources available for Antarctic predators due to commercial whaling in the 20th century has remained largely untested since the 1970s. Rapid warming of the Western Antarctic Peninsula (WAP) over the past 50 years has resulted in decreased seasonal ice cover and a reduction of krill. The latter is being exacerbated by a commercial krill fishery in the region. Despite this, humpback whale populations have increased but may be at a threshold for growth based on these human-induced changes. Understanding how climate-mediated variation in prey availability influences humpback whale population dynamics is critical for focused management and conservation actions. Using an 8-year dataset (2013–2020), we show that inter-annual humpback whale pregnancy rates, as determined from skin-blubber biopsy samples (n = 616), are positively correlated with krill availability and fluctuations in ice cover in the previous year. Pregnancy rates showed significant inter-annual variability, between 29% and 86%. Our results indicate that krill availability is in fact limiting and affecting reproductive rates, in contrast to the krill surplus hypothesis. This suggests that this population of humpback whales may be at a threshold for population growth due to prey limitations. As a result, continued warming and increased fishing along the WAP, which continue to reduce krill stocks, will likely impact this humpback whale population and other krill predators in the region. Humpback whales are sentinel species of ecosystem health, and changes in pregnancy rates can provide quantifiable signals of the impact of environmental change at the population level. Our findings must be considered paramount in developing new and more restrictive conservation and management plans for the Antarctic marine ecosystem and minimizing the negative impacts of human activities in the region.     

N-glycolyl GM1 ganglioside as a receptor for simian virus 40

Carbohydrate microarrays have emerged as powerful tools in analyses of microbe-host interactions. Using a microarray with 190 sequence-defined oligosaccharides in the form of natural glycolipids and neoglycolipids representative of diverse mammalian glycans, we examined interactions of simian virus 40 (SV40) with potential carbohydrate receptors. While the results confirmed the high specificity of SV40 for the ganglioside GM1, they also revealed that N-glycolyl GM1 ganglioside [GM1(Gc)], which is characteristic of simian species and many other nonhuman mammals, is a better ligand than the N-acetyl analog [GM1(Ac)] found in mammals, including humans. After supplementing glycolipid-deficient GM95 cells with GM1(Ac) and GM1(Gc) gangliosides and the corresponding neoglycolipids with phosphatidylethanolamine lipid groups, it was found that GM1(Gc) analogs conferred better virus binding and infectivity. Moreover, we visualized the interaction of NeuGc with VP1 protein of SV40 by molecular modeling and identified a conformation for GM1(Gc) ganglioside in complex with the virus VP1 pentamer that is compatible with its presentation as a membrane receptor. Our results open the way not only to detailed studies of SV40 infection in relation to receptor expression in host cells but also to the monitoring of changes that may occur with time in receptor usage by the virus.     

Isolation,patch size and matrix effects on bird assemblages in forest reserves

Habitat fragmentation is one of the most studied topics in ecology but our knowledge is still limited particularly concerning matrix effects on species distribution in a human-dominated landscape. We tested the ability of random sampling hypothesis, colonization–extinction dynamics and matrix-related concepts to explain the variation in species richness, total bird density and community composition in old-forest bird assemblages in two contrasting landscapes. We collected data on breeding bird abundances from 66 old-growth forest reserves in NE Finland and six larger areas in adjacent Russian Karelia using the line transect method. In Finland, protected old-forest patches are embedded in a matrix dominated by young regeneration stands. In Russia, study areas were situated in continuous, old forest dominated landscapes. Bird assemblages in old-forest patches embedded in human-modified matrix in Finland were not random samples from Russian bird assemblages. In the Finnish assemblages, species richness was lower and total bird density higher. Species richness declined with increasing distance (isolation) to Russia. Bird assemblages in large forest reserves in Finland close to Russia were structurally more similar to assemblages in the continuous reference landscape than those in small and more distant reserves. The results support the idea that several mechanisms related to colonisation–extinction dynamics and to matrix resource availability influence species distribution in fragmented landscapes but in species-specific ways. We conclude that even though small and isolated protected areas may foster high relative bird species density their ecological integrity is compromised, and therefore, improving matrix quality around reserves may lead to considerable conservation benefits.     

The longitudinal development of running economy in males and females aged between 13 and 27 years: The Amsterdam Growth and Health Study

The purpose of this study was to describe the longitudinal development of running economy [defined as the oxygen uptake (V˙O₂) at a submaximal running speed] in males and females from teenage to young adult age using data from the Amsterdam Growth and Health Study. Submaximal V˙O₂ (in ml · kg⁻¹· min⁻¹) was measured in 84 males and 98 females while they ran on a treadmill at a constant speed of 8 km · h⁻¹ for 6 min at three different treadmill slopes (0%, 2.5% and 5%). This test was carried out six times, on the same subjects at the ages of 13, 14, 15, 16, 21, and 27 years. The longitudinal development of running economy in males and females was analysed using a two-way analysis of variance for repeated measurements. At all three slopes, a significant decrease in V˙O₂ with increasing age was found for both males and females, implying a significant increase in running economy for both sexes. Males showed significantly higher V˙O₂ values than females at all ages measured and for all three slopes, suggesting that females have a significantly higher running economy than males. In order to make a better comparison of the V˙O₂ of individuals of different sizes, allometric models were used; power function ratios were constructed in which body mass was expressed to an exponential power. Following this analysis the difference in submaximal V˙O₂ and running economy between males and females appeared even larger. Accepted: 20 February 1997