text:
Research | Open | Published: 05 September 2017
Enrichment of beneficial bacteria in the skin microbiota of bats persisting with white-nose syndrome
Virginie Lemieux-Labonté ORCID: orcid.org/0000-0001-9964-44951,4,
Anouk Simard2,4,
Craig K. R. Willis3 &
François-Joseph Lapointe1,4
Infectious diseases of wildlife are increasing worldwide with implications for conservation and human public health. The microbiota (i.e. microbial community living on or in a host) could influence wildlife disease resistance or tolerance. White-nose syndrome (WNS), caused by the fungus Pseudogymnoascus destructans (Pd), has killed millions of hibernating North American bats since 2007. We characterized the skin microbiota of naïve, pre-WNS little brown bats (Myotis lucifugus) from three WNS-negative hibernation sites and persisting, previously exposed bats from three WNS-positive sites to test the hypothesis that the skin microbiota of bats shifts following WNS invasion.
Using high-throughput 16S rRNA gene sequencing on 66 bats and 11 environmental samples, we found that hibernation site strongly influenced the composition and diversity of the skin microbiota. Bats from WNS-positive and WNS-negative sites differed in alpha and beta diversity, as well as in microbiota composition. Alpha diversity was reduced in persisting, WNS-positive bats, and the microbiota profile was enriched with particular taxa such Janthinobacterium, Micrococcaceae, Pseudomonas, Ralstonia, and Rhodococcus. Some of these taxa are recognized for their antifungal activity, and specific strains of Rhodococcus and Pseudomonas are known to inhibit Pd growth. Composition of the microbial community in the hibernaculum environment and the community on bat skin was superficially similar but differed in relative abundance of some bacterial taxa.
Our results are consistent with the hypothesis that Pd invasion leads to a shift in the skin microbiota of surviving bats and suggest the possibility that the microbiota plays a protective role for bats facing WNS. The detection of what appears to be enrichment of beneficial bacteria in the skin microbiota of persisting bats is a promising discovery for species re-establishment. Our findings highlight not only the potential value of management actions that might encourage transmission, growth, and establishment of beneficial bacteria on bats, and within hibernacula, but also the potential risks of such management actions.
Infectious diseases of wildlife are on the rise worldwide with dramatic consequences for wildlife conservation and human public health [1,2,3]. In North America, insectivorous bats provide important ecosystem services by limiting insect pests and potentially saving billions of dollars annually for agriculture [4, 5]. However, a number of ecologically important species are threatened by white-nose syndrome (WNS). This skin disease, caused by the fungus Pseudogymnoascus destructans (Pd) [6, 7], has killed millions of North American bats since 2006 [8].
White-nose syndrome involves invasion of exposed skin by Pd, and the disease is defined by cup-shaped erosions and ulcerations on the tissue of the flight membranes (wings and tail), ears, and muzzle [9]. Infection of the flight membranes is thought to be the most pathologically significant aspect of the infection because this tissue is involved in fluid balance, thermoregulation, and gas exchange [10]. Pd invades hair follicles and sebaceous and apocrine glands [9]. This likely disrupts secretions that contribute to skin integrity [11, 12] with consequences for defense against pathogens and important skin commensal microorganisms [13]. Hibernating bats survive the winter on just a few grams of stored fat by using prolonged energy-saving bouts of torpor characterized by dramatically reduced body temperatures and metabolism [14,15,16]. Pd is adapted for growth at the low temperature characteristic of bat skin during torpor [6], and infection causes hibernating bats to warm up too frequently during winter and deplete their fat reserves [17, 18]. The immune system is downregulated during hibernation [19,20,21] which, in turn, facilitates infection.
Seven species of bats have suffered impacts from WNS in North America [22] but not all bat species are equally affected [23, 24]. It has been suggested that environmental conditions inside hibernacula, physiology, and behavior could all play a role in the variable tolerance of, or resistance to, infection with Pd among species [22, 24,25,26]. In Canada, the northern long-eared bat (Myotis septentrionalis), the little brown bat (Myotis lucifugus), and the tricolored bat (Perimyotis subflavus) are listed as federally endangered [22] due to mortality rates of 75–90% during the several-year invasion stage of the disease [27].
Despite extremely high mortality during the epidemic stage of WNS, some hibernating colonies of at least one highly vulnerable species (e.g., M. lucifugus) seem to have persisted following disease invasion [24, 28,29,30] with colony counts stabilizing at about 5 to 30% of their initial size [24, 28]. Recently, it was observed that intensity of infection with Pd, based on swabs of bat forearms and quantitative PCR, was significantly lower for persisting colonies in which Pd had become established, compared to colonies in the midst of the epidemic phase and massive declines [31]. One mechanism that could explain this pattern is a fundamental shift in the microbial community living on bat skin due to selection for Pd antagonists. Strong selection for microbial taxa that inhibit Pd could provide resistance to the fungus and increase bat survival.
Animal skin is an ecosystem inhabited by highly variable and complex communities of microorganisms [13]. This community, called microbiota, can be divided into a resident flora, defined as a relatively stable assemblage in size and composition, and a transient flora, acquired from the local environment and that only temporarily colonizes the skin [32]. A healthy skin microbiota can directly contribute to host fitness by occupying pathogen adhesion sites and producing pathogen inhibitors [13, 33]. Competitive interactions between beneficial and pathogenic skin microbes are hypothesized to play a role in disease dynamics for wild animals [34]. For example, the bacterium Janthinobacterium lividum, which lives on salamander skin, appears to confer resistance to the devastating fungal pathogen Batrachochytrium dendrobatidis [35] and could explain why some salamander populations decline while others do not. Recently, a strain of the bacterium Pseudomonas fluorescens isolated from the skin of a bat species thought to be resistant to WNS (Eptesicus fuscus) was shown to inhibit Pd growth in vitro [36] as well as in vivo for M. lucifugus [37]. It has been hypothesized that WNS could cause a shift in microbiota communities of the skin [38], and this could be one mechanism underlying resistance in persisting bats. However, it could also have negative consequences for bat populations if a shift in the microbiota makes it easier for opportunistic pathogens other than Pd to invade the skin. A detailed characterization of the skin microbiota for WNS-positive and WNS-negative bats is, therefore, needed to fully understand potential implications of skin microbial communities in the context of WNS.
Due to its direct exposure to the local environment, the skin microbiota is more dynamic and should be more strongly influenced by the environment, than the gut microbiota [39]. Environment and host species are strong predictors of variation in the skin microbiota among bats [40,41,42]. However, one study [43] found a strong influence of site and habitat type on the skin microbiota of 13 bat species in the western USA but was not able to detect an influence of host species or sex. For bats and amphibians, the local environment appears to act as a reservoir for skin microbiota, while conditions on the skin may lead to selection favoring or enriching particular taxa [40, 44, 45]. Consequently, host and local environmental factors appear to interact closely to shape the skin microbiota. This suggests that the skin microbiota could exhibit dramatic temporal variation for species characterized by seasonal shifts in physiology and habitat selection. Bats exhibit enormous changes in metabolism and habitat selection between winter and summer. Therefore, to fully characterize the skin microbiota and its relevance to WNS, bats must be sampled at the appropriate time during hibernation.
Several studies have reported on the skin microbiota of North American bats, but, to date, these have involved individuals not yet affected by WNS [42, 43] or have been based on relatively small sample sizes [38, 40]. Our objective was to understand the potential interaction between Pd and the skin microbiota of bats by comparing individuals from WNS-positive and WNS-negative regions. We used high-throughput 16S amplicon sequencing to characterize the composition and diversity of the skin microbiota of M. lucifugus sampled from WNS-positive (Québec) and WNS-negative (Manitoba) hibernacula in the northern part of this species' range, in Canada. We tested two predictions of the hypothesis that WNS is causing selection favoring Pd antagonists on the skin microbiota of bats in the affected region. First, we predicted that bats persisting in WNS-affected sites would exhibit reduced diversity of their microbiota consistent with strong selection for a subset of pre-WNS microbial species [38]. Second, we predicted that the microbiota of persisting bats from WNS-affected sites would show a proportional increase in antifungal/anti-Pd bacterial species such as those identified in previous studies [36, 37, 46]. We also tested the third hypothesis that variation in the skin microbiota of hibernating bats relates to environmental variation in the microbial community of a given cave. We predicted that the diversity and composition of the microbial community living on individual bats would be similar to that found on substrates in the local environment of their hibernaculum and would differ from that on bats and the local environment in other hibernacula.
Sampling and ethics
During winter 2015–2016, we sampled the skin microbiota of 33 M. lucifugus from three WNS-negative hibernacula in central Manitoba (Canada) about 50 km north of the town of Grand Rapids (53° 30′ N; 99° 24′ W) and another 33 individuals from three sites known to be WNS-positive since 2010 in Québec (Canada) within 60 km north of Gatineau city (45° 28′ N; 75° 42′ W). The temperature within sites ranged from −3 to 7 °C at sampling time. Site and population information are specified in Table 1.
Table 1 M. lucifugus hibernaculum sites information in Manitoba and Québec provinces
Bats in a given hibernaculum were always sampled from within the same area (i.e., room, gallery, corridor). We selected bats at random from among those we could reach from the ground. Little brown bats are highly gregarious during hibernation, and most individuals spend at least part of their time huddling or clustering during hibernation. We defined bats in direct contact with each other as being members of the same cluster. Sixty-four of the 66 bats we sampled were clustering with other bats, and cluster sizes ranged in size from 2 to 11 individuals. We sampled two bats at Emerald that were roosting solitarily (Additional file 1). We swabbed 11 individual bats per site. Samples were collected by swabbing in linear strokes the back and forearm of each bat for 20 s with a sterile Whatman Omniswab (Fisher Scientific) soaked in sterile 0.15 M NaCl [41]. Swab tips were ejected into MoBio Powersoil DNA isolation Kit tubes (MoBio Laboratories), which were transferred to −20 °C within 24 h of sampling until DNA extraction. Local environment samples were also collected by swabbing cave walls adjacent to clusters of sampled bats for 20 s in linear strokes (approx. 5 cm). As a negative control, a humidified sterile swab was exposed to open air for 20 s, prior to ejecting its tip into a MoBio tube.
Bats are vulnerable to disturbance during hibernation, and we were careful to minimize the impact of our visits. Only two people entered hibernacula for sampling, and bats were not handled during swabbing so we did not determine their sex. A previous study established that sex was not a significant predictor of the external microbiota of bats [43]. Therefore, differences among hibernacula are likely to reflect the influence of the local habitat (e.g., differences in temperature, humidity, and environmental bacteria), rather than difference in sex ratio among sites. All methods were approved by the Animal Welfare and Ethics Committee at Université de Montréal (Protocol Number #16-015) and the University of Winnipeg Animal Care Committee (Protocol Number AEO5639).
DNA extraction, amplification, and sequencing
Bacterial genomic DNA was extracted from each swab using the MoBio Powersoil DNA isolation Kit according to the manufacturer's protocol. Extractions were randomized for site and region to avoid detecting false patterns [47]. Extraction, amplification blanks, and the HM-782D Human Microbiome Project mock community (BEI Resources) were also included to detect possible contamination and assess sequencing accuracy [47, 48]. Amplification and sequencing were then performed as previously described [49]. Libraries were prepared using a two-step PCR. The first PCR amplified the hypervariable region V4 of the 16S small subunit ribosomal gene with forward primer U515_f: ACACGACGCTCTTCCGATCTYRYRGTGCCA GCMGCCGCGGTAA and reverse primer E786_R: CGGCATTCCTGCTGAACCGCTCTTCC GATCTGGACTACHVGGGTWTCTAAT [50]. Two microliters of extracted DNA (equivalent DNA amount by sample) was added to the PCR reaction containing 14.25 μl of sterile water, 5 μl HF buffer, 0.5 μl DNTPs, 0.25 μl Phusion High-Fidelity DNA Polymerase (New England Biolabs Inc.), and 1.5 μl of forward and reverse primers. Amplifications were performed with a Mastercycler Nexus GSX1 (Eppendorf) under the following conditions: initial denaturation at 98 °C for 30 s; 30 cycles alternating 98 °C for 25 s, 40 s at 54 °C, 35 s at 72 °C, and final elongation step for 1 min at 72 °C. Each sample was amplified in quadruplicate and pooled to limit possible PCR artifacts. All PCR products were then purified by PCR purification Agencourt AMPure XP (Beckman Coulter). The second PCR step consisted of adding primers containing a barcode (index) and Illumina adapter sequences to each DNA amplicon. To do so, 4 μl of the first step amplification product was added to a PCR reaction containing 10.25 μl of sterile water, 5 μl HF buffer, 0.5 μl DNTPs, 0.25 μl Phusion High-Fidelity DNA Polymerase, and 2.5 μl of forward primer PE-III-PCR-F:AATGATACGGCGACCACCGAGATCTACACTCTTT CCCTACACGACGCTCTTCCGATCT and reverse primer PE-III-PCR-001-096:CAAGCAGA AGACGGCATACGAGATNNNNNNNNNCGGTCTCGGCATTCCTGCTGAACCGCTCTTCCGATCT (N indicating the unique barcode) [51]. Indexing was performed under the following thermal conditions: initial denaturation at 98 °C for 30 s, 7 cycles alternating 98 °C for 30 s, 30 s at 83 °C, and finally 30 s at 72 °C. This second amplification was performed in triplicate. Samples were pooled and purified with the PCR purification Agencourt AMPure XP (Beckman Coulter). Qubit 2.0 Fluorometer (Invitrogen) was used to measure the DNA concentration of each sample. Indexed samples were then pooled to obtain a final concentration range between 10 and 20 ng/μl. DNA was next diluted and denatured according to the manufacturer's protocol for paired-end sequencing using MiSeq Reagent Kit v2 (500 cycles) 2 × 250 bp on MiSeq (Illumina).
We amplified 4,072,792 sequences classified into 13,224 operational taxonomic units (OTUs) from the 66 swabs of bat skin and 11 environmental samples (one or two per site). A total of 3,729,096 sequences classified in 11,812 OTUs were amplified from bat samples with a mean of 56,501 sequences per sample (range 9920–100,812). A total of 343,696 sequences classified in 9302 OTUs were obtained from the 11 environmental samples, with a mean of 31,245 sequences per sample (range 10,325–73,877). We were able to match all expected sequences in the mock positive control, except for Helicobacter pylori, which genus was nonetheless the most abundant in the compositional data (see Additional file 2). The genera or families of the 20 expected mock taxa were also the most abundant in the mock profile. We detected 36 false positives, with very low abundances (< 0.3%) (see Additional file 2). After filtering out OTUs with abundance values smaller than 3, sampling controls, extraction controls, and library negative controls were dominated by Halomonas (5–75%, mean of 56%) and Shewanella genera (1–26%, mean 18%) (see Additional file 3).
Preclustering, quality filtering, primer removal, merging of raw sequences, and postclustering dereplication were performed with the SmileTrain scripts [52] for 16S data processing using USEARCH v. 7.0.1090 [53]. Distribution-based clustering using the dbOTUcaller algorithm was performed to cluster sequences into OTUs by considering the distribution of DNA sequences across samples and distances between sequences [51]. The corresponding OTU table, providing abundances of bacterial taxa in the different samples was assigned with QIIME version 1.9.0 [54] using GreenGenes database release 13_5 [55]. For alpha diversity and compositional analysis of bat samples, mitochondrial and chloroplastic DNA sequences, as well as OTUs with abundance values smaller than 3, were filtered out, leaving 3,716,672 sequences classified into 9897 OTUs. In addition, the genera Halomonas and Shewanella, present in negative controls, were filtered out from all bat samples for compositional analysis, resulting in 3,145,399 sequences classified into 9575 OTUs.
The diversity of the skin microbial community (alpha diversity) of each sample was computed using the Shannon index [56]. The Shannon index, which includes both OTU richness and evenness, was selected due to its reduced sensitivity to sample depth differences [49, 57] (Additional file 4). R version 3.1.3 [58] was used for all statistical analyses. Log-transformed alpha diversity values were compared between WNS-positive and WNS-negative regions, using a linear mixed-effect model (lme() function), and significance was tested with anova.lme() of the nlme package [59]. Hibernaculum and clusters were included as a random effect. Variation in diversity among sites within the WNS-positive and WNS-negative regions was tested using a one-way ANOVA (function aov()) and post hoc Tukey test (function TukeyHSD()) of the package stats [58].
The change in diversity among skin microbial community (beta diversity) was calculated among skin microbiota samples and environmental samples. Two distinct phylogenetic distances, unweighted UniFrac (qualitative) and weighted UniFrac (quantitative) [60, 61], were computed on rarefied data, as such measures could be sensitive to differences in sequencing depth [62, 63]. UniFrac distances were computed from bat samples rarefied at 9886 sequences/sample and from environmental + bat samples rarefied at 9898 sequences/sample after retrieving OTUs in low abundance (<3 sequences). Computations were performed with the phyloseq package [64]. All beta diversity results were visualized with principal coordinates analysis (PCoA) [65] using the ordinate() function. The UniFrac distance matrix was checked with is.euclid() function of the ade4 package [66] prior to the ordination to ensure that all distances were Euclidian and properly represented by PCoA [67]. When required, square-root transformations were applied to obtain distance matrices satisfying the Euclidian condition. All phylogeny-based UniFrac distances were calculated using a phylogenetic tree constructed with FastTree 2.1.8 [68].
To assess the influence of explanatory variables on the microbiota composition, we used distance-based redundancy analysis (db-RDA), a method intended to conduct a redundancy analysis (RDA) on distance matrices [69]. It is computed by first decomposing UniFrac distances (weighted or unweighted) into principal coordinates and then applying RDA to the corresponding principal coordinates using the capscale() function of the R package vegan [70]. Four distinct models were constructed to test the relative importance of (1) WNS status (i.e., WNS-negative vs. WNS-positive), (2) sampling sites (i.e., the six different hibernacula), (3) types of samples (i.e., bat samples vs. local environment samples), and (4) clusters (i.e., bat clusters within each hibernaculum). To better understand the relationships among explanatory models in the variation of the microbial assemblages, partial db-RDA was also computed [71]. This form of RDA allows for exploration of the contribution of an explanatory variable model while controlling for other explanatory models. Adjusted R-squared (R 2) values [72] were calculated to compare the explanatory power of such models containing different numbers of variables. Significance of db-RDA and partial db-RDA was tested via 9999 permutations with the anova.cca() function of the R package vegan .
The microbiota composition was explored down to genus level to assess differences among hibernacula and between WNS-positive and WNS-negative sites. To emphasize these differences, Indicator Value tests (IndVal) [73] were performed on relative abundance data, using the 26 taxa with a relative abundance larger than 1% for the analysis. The IndVal indicator value is based on the comparison of occurrences and abundances of taxa across predefined groups of bats (e.g., grouped by sites or WNS status). The analysis for any given taxon is not influenced by other taxa present in the dataset. It provides an index ranging between 0 and 1, the maximum value indicating a taxon exclusively present in one group. IndVal is calculated as the product of A (specificity, i.e., the probability that a site belongs to the group given the fact that a given species is found in that site) and B (fidelity, i.e., the probability of finding a given taxon at a site when the site belongs to that group) [73]. The multipatt() function of the R package indicspecies [74] was used to compute indicator values, and significance was assessed with 9999 permutations of object between groups. The p.adjust() function of the R package stats was used to correct p values for multiple comparisons [75]. A corrected p value threshold of 0.05 was considered significant in all tests, and only significant taxa with a specificity of A ≥ 0.4 were retained as indicators.
Alpha diversity in WNS-positive and WNS-negative regions
After controlling for sites and bat clusters using a linear mixed-effects model, we found significant differences in Shannon diversity between our pooled set of WNS-positive hibernacula in Québec and WNS-negative hibernacula in Manitoba (F 1,4 = 16.27, p ≤ 0.05) (Fig. 1a). WNS-positive sites had significantly lower Shannon diversity than WNS-negative sites (Fig. 1a). We also found significant variation in alpha diversity between some of the hibernacula in the WNS-negative region (ANOVA: F 2,30 = 22.84, p ≤ 0.001), whereas all the three WNS-positive sites in Québec were statistically indistinguishable from each other and had relatively low alpha diversity (ANOVA: F 2,30 = 0.96, p = 0.395) (Fig. 1b). Within the WNS-negative region in Manitoba, Abyss cave harbored particularly high alpha diversity and was significantly different from Dale's (Tukey's; p ≤ 0.001) and Microwave (Tukey's: p ≤ 0.001).
Alpha diversity of M. lucifugus skin microbiota in WNS-positive and WNS-negative sites in Canada. Distribution of alpha diversity within groups as estimated by the Shannon index for a hibernacula pooled by WNS status (positive vs. negative) and b all six hibernacula sampled in the study. Error bars represent standard deviations. Significant differences in alpha diversity among groups are indicated by different letters according to linear mixed model effect, ANOVA, and Tukey's test (p ≤ 0.05)
Beta diversity analysis of microbial community assemblage
We first used beta diversity analysis to explore compositional differences among skin microbiota samples alone, that is, after removing all environmental samples from the analysis. The PCoA, based on unweighted UniFrac, revealed a clear separation between WNS-positive sites in Québec and WNS-negative sites in Manitoba, and also grouped samples from within the same hibernaculum (Fig. 2a). This pattern was not observed with weighted UniFrac (Fig. 2b), which implies that accounting for differential abundances (weighted UniFrac), and not just the presence/absence of bacterial OTUs between samples (unweighted UniFrac), affected our results. However, the first principal axes, accounting for 20.1% of the variation in the data, support the separation of microbiota samples according to WNS status (Fig. 2a).
Principal coordinate analysis of M. lucifugus skin microbiota in WNS-positive and WNS-negative sites. a Principal coordinate analysis of unweighted UniFrac distances. b Principal coordinate analysis of weighted UniFrac distances. Each point represents a sample from an individual bat hibernating in one of the six different hibernacula that differed in WNS status
In order to better relate these patterns to different variables, we used a distance-based redundancy analysis (db-RDA) to compute from UniFrac distances among skin microbiota samples using three distinct explanatory models: (1) WNS status, (2) sampling sites (hibernacula), and (3) clusters (bat clusters within each hibernaculum). Unweighted UniFrac distances revealed that each of these models explained a significant fraction of microbiota community variation (Table 2). The WNS status model explained 8%, sites explained 22%, and bat cluster explained 28% of microbiota community variation among samples. Weighted UniFrac distances accounting for abundance of taxa revealed similar patterns with WNS status explaining 14%, sites explaining 26%, and cluster explaining 30% of the variation in the microbiota samples (Table 2).
Table 2 db-RDA of unweighted and weighted UniFrac distances of M. lucifugus skin microbiota samples
In light of these results, we conducted partial RDA to better distinguish the relative influence of our three explanatory models (Table 2). This analysis revealed that WNS status explained no variation after controlling for sites and/or clustering. Similarly, the sites model explained none of the variation after controlling for bat cluster. The cluster model, on the other hand, explained a significant fraction of variation in microbiota community after controlling for sites, both with weighted and unweighted UniFrac distances. These results suggest that bat clustering within each hibernaculum exerts strong influence on the composition of the skin microbiota. Controlling for WNS status, however, greatly reduced the variation explained by sites and cluster models alone, regardless of UniFrac distances (unweighted or weighted). Taken together, the results of simple and partial db-RDA analyses suggest that sites, combined with a local effect of bat clustering within sites, contribute to shaping the skin microbiota, whereas WNS status have much less influence.
We next used unweighted and weighted UniFrac to explore the relationship between skin microbiota samples and environmental samples collected at each site. The first PCoA, based on unweighted UniFrac distances, grouped skin microbiota samples and environmental samples by sites, with some overlap between them (Fig. 3a). The second PCoA based on weighted UniFrac distances revealed a different pattern, however. In that case, the PCoA plot clearly distinguished between environmental samples and skin microbiota samples collected at all sites (Fig. 3b), except for three bat samples from Abyss and one from Dale's. The results of PCoA suggest that the presence/absence of OTUs in bat skin samples is influenced by the local environment within each hibernaculum. Yet, the same analyses also reveal differences in abundance patterns of some OTUs in bat skin samples compared to local environmental samples.
Principal coordinate analysis comparing local environment sites samples and M. lucifugus skin microbiota in Canada. a Principal coordinate analysis of unweighted UniFrac distances. b Principal coordinate analysis of weighted UniFrac distances. Each point represents a single sample
We then used RDA to explore the influence of site (hibernaculum) and sample type (bat vs. environment samples) on variation in microbial community assemblage. Both of these models were significant (Table 3), but the sites model accounted for the most variation in the data, explaining 18% of the microbial community variation in both UniFrac distances employed (Table 3). The sample type model only explained 5% of microbial variation for unweighted UniFrac distances and 8% for weighted distances. The higher explanatory power of the weighted UniFrac model indicates that differences observed between environmental samples and bat samples partly depend on the relative abundance of each taxon within the corresponding microbial communities. This is consistent with patterns revealed by the PCoA plots (Fig. 3).
Table 3 db-RDA of unweighted and weighted UniFrac distances among local environment and bat skin microbiota samples
We used a partial db-RDA to better understand the relationship between the two explanatory models, and their influence on the composition of microbial communities. In both cases, when variation of one model was controlled for using partial db-RDA, the ability of the models to explain variation in microbial community composition was slightly increased by 1% (Table 3). These results suggest that both sample type and sites models had a non-redundant influence on microbial community variation and that the local environment is an important factor explaining skin microbiota patterns of hibernating bats in our study sites.
Taxonomic indicators of WNS status
We found that the most abundant bacterial taxa were shared among all hibernacula, but we also identified indicator taxa present more often and more abundant at particular sites. At the phylum level, the dominant taxa accounting together for 86 to 98% of overall profiles in a given cave were Actinobacteria (23 to 53%), Proteobacteria (24 to 51%), and Bacteroidetes (6 to 38%) (Additionalfile 5). At the class level, the six principal taxa accounting for 80 to 98% of the total abundance were Actinobacteria, Gammaproteobacteria, Flavobacteriia, Sphingobacteriia, Alphaproteobacteria, and Betaproteobacteria (Additional file 6).
Generalist genera such as Arthrobacter, Chryseobacterium, Flavobacterium, Intrasporangiaceae, Pedobacter, Mycoplana, Pseudonocardiaceae, Ralstonia, Rhodococcus, Sinobacteraceae, and Sphingobacterium were identified from all sites (Additional file 7). Significant representatives were found among the 26 more abundant taxa representing more than 1% of the total composition profile (Fig. 4, Additional file 8). Among the more abundant taxa identified at three WNS-positive sites in Québec, only Pseudomonas and Acinetobacter were indicators of one site (Fig. 4, Additional file 8). On the other hand, the more abundant taxa at three WNS-negative sites in Manitoba, Knoellia, Brucellaceae:Other, Microbacterium, and Pseudomonadaceae were all indicators of the Microwave site (Fig. 4, Additional file 8). The largest indicator value was obtained for Nitrosovibrio at the Abyss site. Cytophagaceae and Flavobacteriaceae were also associated with Abyss. Representative taxa were identified from all hibernacula, except for Emerald, Laflèche and Dale's.
M. lucifugus skin microbiota taxa indicator of the six hibernacula of different WNS status in Canada. The sites from WNS-negative (Manitoba, Canada) and WNS-positive (Québec) regions are presented. The significant indicators were identified by IndVal analysis among the 26 more abundant taxa representing more than 1% of total abundance. Stars indicate hibernacula with significant representative taxa. *IndVal > 0.60, **IndVal L ≥ 0.75, ***IndVal ≥ 0.89
We compared skin microbiota profiles based on WNS status in order to highlight possible differences in microbial composition related to the fungal disease. Here, again, some of the most abundant taxa such as Pedobacter and Intrasporangiaceae were not significant representatives as they were identified in both areas in similar relative abundance (Additional file 7). However, a large number of significant indicators were detected, some with high indicator values. At WNS-negative sites in Manitoba, significant indicators were Knoellia, Brucellaceae:Other, Nitrosovibrio, Flavobacteriaceae, Enterobacteriaceae, Microbacterium, Sphingobacterium, Cytophagaceae, Chryseobacterium, and Xanthomonadaceae (Fig. 5, Additional file 9). On the other hand, significant indicators of WNS-positive sites in Québec were Ralstonia, Janthinobacterium, Rhodococcus, Micrococcaceae, and Pseudomonas (Fig. 5, Additional file 9).
M. lucifugus skin microbiota taxa indicator of WNS-positive (Québec) and WNS-negative regions (Manitoba) in Canada. Significant indicators were found among the 26 more abundant taxa representing more than 1% of total abundance with IndVal analysis. Stars indicate regions with significant representative taxa. **IndVal ≥ 0.75, ***IndVal ≥ 0.89
We compared the skin microbiota of bats from WNS-positive and WNS-negative sites to better understand the role of the microbiota as a factor in the host-pathogen interaction associated with WNS. We found support for the hypotheses that WNS has contributed changes in the skin microbiota for bats that are persisting in affected regions and that the skin microbiota is strongly influenced by the local environment within hibernacula. Although we cannot rule out the role of geographic variation confounded here with WNS status we believe that our results are more consistent with the proposal that WNS has led to a shift in the microbiota of bats inhabiting WNS-positive sites. For one, our analyses based on weighted UniFrac distances show no clustering environmental samples by province when analyzed together with bat samples, supporting that region is not a major driver of microbiota communities. Second, because Pd affects and interacts with the skin so directly and because higher levels of bacteria known to inhibit Pd growth in vitro [36, 46] and in vivo [37] were observed, it seems more likely that WNS status, and not geography, explains the compositional patterns. As predicted by our first hypothesis, the diversity of the skin microbiota was indeed smaller at WNS-positive sites compared to WNS-negative sites, which is consistent with a shift in microbiota caused by Pd. In addition, WNS status was a strong predictor of variation in Shannon diversity values across sites. A previous study on tricolored bats (P. subflavus) affected by WNS also revealed a trend for lower diversity values at WNS-positive sites [38], as shown here for little brown bats persisting after WNS invasion. Phylogenetic beta diversity analysis was also consistent with selection on the microbiota by Pd and WNS. Future studies, assessing diversity of the microbiota on bats from the same sites, before and after Pd invasion, would help resolve the WNS influence in microbiota diversity patterns of persisting bats.
At the compositional level, the skin microbiota of hibernating little brown bats is dominated by the classes Actinobacteria, Gammaproteobacteria, Flavobacteriia, Alphaproteobacteria, and Betaproteobacteria, a pattern consistent with previous investigation of the skin microbiota in several species of bats [42] and particularly M. lucifugus [40]. Our analysis also identified Sphingobacteriia as a predominant class. IndVal analysis revealed that one interesting genus, Rhodococcus, was significantly more abundant in skin microbiota samples collected at WNS-positive sites in Québec. This genus has previously been identified on bats [38, 76] and is known for its antifungal activity [77, 78]. Most interesting, a volatile organic chemical produced by R. rhodochrous strain DAP 96253 has been shown to inhibit Pd growth in vitro [46]. Several other genera, reported as antifungal agents, were also identified as significant indicators of bat samples collected at WNS-positive sites. Namely, Pseudomonas was enriched at all WNS-positive sites, whereas Acinetobacter was enriched at a single WNS-positive site. Both taxa are known for their antifungal activity [79, 80] and have been previously identified on the skin of North American bats [38, 40]. Moreover, one strain of Pseudomonas fluorescens has been shown to inhibit Pd growth in vitro and reduce disease severity and improve survival of bats with WNS in a laboratory challenge experiment [36, 37]. Another lesser-known antifungal bacterial genus, Janthinobacterium, was also identified as a significant representative at WNS-positive sites in Québec. Some species from the same genus isolated from the skin of wild amphibians confer resistance against the fungal pathogen Batrachochytrium dendrobatidis [35, 79]. Enrichment of multiple taxa with potential antifungal and anti-Pd activity in bats persisting following WNS invasion is consistent with our second hypothesis that the skin microbiota of bats provides a mechanism for resistance to, or tolerance of, Pd infection. Further studies should focus on any functional influence of these bacteria on the host-pathogen interaction between bats and Pd.
Consistent with our third hypothesis, we found that the skin microbiota of hibernating little brown bats is related to the microbial community composition of the nearby environmental substrates. That is, bacteria living on bats and bacteria living on adjacent cave walls are very likely exchanged by contact. However, bat skin samples and local environmental samples were by no means identical in their compositional profiles, indicating that microbial communities on the skin of hibernating bats are probably not regulated in the same way as in the environment (Additional file 10). These results are consistent with other studies of bats [40] or frogs [44, 45] showing that skin microbiota assemblages do not exactly mirror the microbial communities in the immediate environment. Although we did not detect any related variation within sites of microbial communities of the skin and that of the substrates, we found that individual bats strongly differ across sites. Moreover, the tendency for hibernating little brown bats to cluster, often in large groups of hundreds to thousands of individuals, is likely to reduce variation in the microbiota among individuals because of transfer within clusters, as shown by our analysis. Homogenization of the skin microbiota by close contact among individuals has also been observed in previous studies of bats and humans [41, 81]. Taken together, these results suggest that bat populations could differ in their susceptibility to WNS depending on the microbial community in their immediate environment, their reliance on clustering behavior, and the potential for clustering to homogenize the bacterial community. In this study, we did not attempt to quantify the potential influence of abiotic variables, such as pH, temperature, and humidity, and considered these factors as possible contributors to site effects. It would be interesting in future studies to analyze these factors separately to understand their relative influence on the microbial community on bats and in the environment. Temporal variation may also have influenced the compositional patterns observed in this study, but our experimental protocol was designed to reduce this effect as much as possible. All sites were sampled within a relatively short period of time (less than 3.5 months), and we avoided the start of hibernation when the experience of bats prior to hibernation might be expected to more strongly influence their skin microbiota. Moreover, the microbiota on bats or in the environment for the single WNS-positive site we sampled in March (Emerald) was not different from the two WNS-positive sites we sampled in November (Lames and Laflèche).
This study highlights the role of skin microbiota for wildlife population health, conservation, and management in the face of emerging infectious diseases. The enrichment of potentially beneficial bacteria in skin microbiota samples collected at WNS-positive hibernacula is an encouraging discovery for the prospect of bat population recovery after WNS becomes endemic in a given region. This finding highlights the potential value of management actions that might encourage transmission, growth, and establishment of beneficial bacterial taxa on bats and within hibernacula [37]. However, our findings also highlight a potential risk of some proposed management actions. Considerable funding and time is currently being devoted to development and testing of potential chemical or biological treatments for WNS that could be applied to hibernating bats or hibernaculum substrates. Our results not only support previous work highlighting the potential of some bacteria as biological control agents for Pd (e.g., [36, 37, 46]) but also highlight a potential risk of biological or chemical treatments. Treatments that disrupt the skin microbiota or attenuate selection for a beneficial skin community could cause more harm than good for recovery of bat populations and the establishment of stable, long-term resistance to WNS in the wild. Thus, an important component of testing any potential treatment for WNS should be to confirm that it is selective for Pd and has minimal negative impacts on the whole non-target microbiota on bats or in hibernacula.
Pd :
Pseudogymnoascus destructans
WNS:
White-nose syndrome
Fisher MC, Henk DA, Briggs CJ, Brownstein JS, Madoff LC, McCraw SL, et al. Emerging fungal threats to animal, plant and ecosystem health. Nature. 2012;484(7393):186–94. https://doi.org/10.1038/nature10947.
Daszak P, Cunningham AA, Hyatt AD. Emerging infectious diseases of wildlife-threats to biodiversity and human health. Science. 2000;287(5452):443–9. https://doi.org/10.1126/science.287.5452.443.
Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, et al. Global trends in emerging infectious diseases. Nature. 2008;451(7181):990–3. https://doi.org/10.1038/nature06536.
Boyles JG, Cryan PM, McCracken GF, Kunz TH. Economic importance of bats in agriculture. Science. 2011;332(6025):41–2. https://doi.org/10.1126/science.1201366.
Maine JJ, Boyles JG. Bats initiate vital agroecological interactions in corn. Proc Natl Acad Sci. 2015;112(40):12438–43. https://doi.org/10.1073/pnas.1505413112.
Gargas A, Trest MT, Christensen M, Volk TJ, Blehert DS. Geomyces destructans sp. nov. associated with bat white-nose syndrome. Mycotaxon. 2009;108(1):147–54. https://doi.org/10.5248/108.147.
Lorch JM, Meteyer CU, Behr MJ, Boyles JG, Cryan PM, Hicks AC, et al. Experimental infection of bats with Geomyces destructans causes white-nose syndrome. Nature. 2011;480:376–8. https://doi.org/10.1038/nature10590.
U.S. Fish & Wildlife Service. North American bat death toll exceeds 5.5 million from white-nose syndrome. In: News release. 2012. http://www.batcon.org/pdfs/USFWS_WNS_Mortality_2012_NR_FINAL.pdf. Accessed 28 Dec 2016.
Meteyer CU, Buckles EL, Blehert DS, Hicks AC, Green DE, Shearn-Bochsler V, et al. Histopathologic criteria to confirm white-nose syndrome in bats. J Vet Diagn Investig. 2009;21(4):411–4. https://doi.org/10.1177/104063870902100401.
Cryan PM, Meteyer CU, Boyles JG, Blehert DS. Wing pathology of white-nose syndrome in bats suggests life-threatening disruption of physiology. BMC Biol. 2010;8:135. https://doi.org/10.1186/1741-7007-8-135.
Cortese TA, Nicoll PA. In vivo observations of skin appendages in the bat wing. J Invest Dermatol. 1970;54:1–10. https://doi.org/10.1111/1523-1747.ep12551469.
Sisk MO. A study of the sudoriparous glands of the little brown bat, Myotis lucifugus lucifugus. J Morphol. 1957;101:425–55. https://doi.org/10.1002/jmor.1051010303.
Grice EA, Segre JA. The skin microbiome. Nat Rev Microbiol. 2011;9(4):244–53. https://doi.org/10.1038/nrmicro2537.
Czenze ZJ, Willis CKR. Warming up and shipping out: cues for arousal and emergence in hibernating little brown bats (Myotis lucifugus). J Comp Physiol B. 2015;185:575–86. https://doi.org/10.1007/s00360-015-0900-1.
Czenze ZJ, Park AD, Willis CKR. Staying cold through dinner: cold-climate bats rewarm with conspecifics but not sunset during hibernation. J Comp Physiol B. 2013;183(6):859–66. https://doi.org/10.1007/s00360-013-0753-4.
Jonasson KA, Willis CKR. Hibernation energetics of little brown bats. J Exp Biol. 2012;215:2141–9. https://doi.org/10.1242/jeb.066514.
Reeder DM, Frank CL, Turner GG, Meteyer CU, Kurta A, Britzke ER, et al. Frequent arousal from hibernation linked to severity of infection and mortality in bats with white-nose syndrome. PLoS One. 2012;7(6):e38920. https://doi.org/10.1371/journal.pone.0038920.
Warnecke L, Turner JM, Bollinger TK, Lorch JM, Misra V, Cryan PM, et al. Inoculation of bats with European Geomyces destructans supports the novel pathogen hypothesis for the origin of white-nose syndrome. Proc Natl Acad Sci. 2012;109(18):6999–7003. https://doi.org/10.1073/pnas.1200374109.
Moore MS, Reichard JD, Murtha TD, Zahedi B, Fallier RM, Kunz TH. Specific alterations in complement protein activity of little brown myotis (Myotis lucifugus) hibernating in white-nose syndrome affected sites. PLoS One. 2011;6(11):e27430. https://doi.org/10.1371/journal.pone.0027430.
Bouma HR, Carey HV, Kroese FG. Hibernation: the immune system at rest? J of Leukoc Biol. 2010;88(4):619–24. https://doi.org/10.1189/jlb.0310174.
Geiser F. Metabolic rate and body temperature reduction during hibernation and daily torpor. Annu Rev Physiol. 2004;66:239–74. https://doi.org/10.1146/annurev.physiol.66.032102.115105.
Frick WF, Puechmaille SJ, Willis CKR. White-nose syndrome in bats. In: Voigt CC, Kingston T, editors. Bats in the Anthropocene: conservation of bats in a changing world. Cham, Heidelberg, New York, Dordrecht, London: Springer International Publishing; 2016. p. 245–62.
Turner GG, Reeder DM, Coleman JTH. A five-year assessment of mortality and geographic spread of white-nose syndrome in North American bats, with a look to the future. Bat Res News. 2011;52:13–27.
Langwig KE, Frick WF, Bried JT, Hicks AC, Kunz TH, Kilpatrick AM. Sociality, density-dependence and microclimates determine the persistence of populations suffering from a novel fungal disease, white-nose syndrome. Ecol Lett. 2012;15(9):1050–7. https://doi.org/10.1111/j.1461-0248.2012.01829.x.
Willis CKR, Wilcox A. Hormones and hibernation: possible links between hormone systems, winter energy balance and white-nose syndrome in bats. Horm Behav. 2014;66(1):66–73. https://doi.org/10.1016/j.yhbeh.2014.04.009.
Puechmaille SJ, Frick WF, Kunz TH, Racey PA, Voigt CC, Wibbelt G, et al. White-nose syndrome: is this emerging disease a threat to European bats? Trends Ecol Evol. 2011;26:570–6. https://doi.org/10.1016/j.tree.2011.06.013.
COSEWIC. COSEWIC Assessment and Status Report on the little brown Myotis lucifugus, Northern Myotis septentrionalis and tri-colored bat Perimyotis subflavus in Canada. Committee on the Status of Endangered Wildlife in Canada. 2013. http://www.registrelep-sararegistry.gc.ca/document/default_e.cfm?documentID=1323. Accessed 4 Jan 2017.
Frick WF, Puechmaille SJ, Hoyt JR, Nickel BA, Langwig KE, Foster JT, et al. Disease alters macroecological patterns of North American bats. Glob Ecol Biogeogr. 2015;24(7):741–9. https://doi.org/10.1111/geb.12290.
Dobony CA, Hicks AC, Langwig KE, von Linden RI, Okoniewski JC, Rainbolt RE. Little brown myotis persist despite exposure to white-nose syndrome. J of Fish Wildl Manag. 2011;2(2):190–5. https://doi.org/10.3996/022011-JFWM-014.
Maslo B, Valent M, Gumbs JF, Frick WF. Conservation implications of ameliorating survival of little brown bats with white-nose syndrome. Ecol Appl. 2015;25(7):1832–40. https://doi.org/10.1890/14-2472.1.
Langwig KE, Hoyt JR, Parise KL, Frick WF, Foster JT, Kilpatrick AM. Resistance in persisting bat populations after white-nose syndrome invasion. Philos Trans R Soc B Biol Sci. 2017;372(1712):20160044. https://doi.org/10.1098/rstb.2016.0044.
Price PB. The bacteriology of normal skin; a new quantitative test applied to a study of the bacterial flora and the disinfectant action of mechanical cleansing. J Infect Dis. 1938;63(3):301–18. http://www.jstor.org/stable/30088420.
Roth RR, James WD. Microbial ecology of the skin. Annu Rev Microbiol. 1988;42:441–64. https://doi.org/10.1146/annurev.mi.42.100188.002301.
Belden LK, Harris RN. Infectious diseases in wildlife: the community ecology context. Front Ecol Environ. 2007;5(10):533–9. https://doi.org/10.1890/060122.
Brucker RM, Harris RN, Schwantes CR, Gallaher TN, Flaherty DC, Lam BA, et al. Amphibian chemical defense: antifungal metabolites of the microsymbiont Janthinobacterium lividum on the salamander Plethodon cinereus. J Chem Ecol. 2008;34(11):1422–9. https://doi.org/10.1007/s10886-008-9555-7.
Hoyt JR, Cheng TL, Langwig KE, Hee MM, Frick WF, Kilpatrick AM. Bacteria isolated from bats inhibit the growth of Pseudogymnoascus destructans, the causative agent of white-nose syndrome. PLoS One. 2015;10(4):e0121329. https://doi.org/10.1371/journal.pone.0121329.
Cheng TL, Mayberry H, McGuire LP, Hoyt JR, Langwig KE, Nguyen H, et al. Efficacy of a probiotic bacterium to treat bats affected by the disease white-nose syndrome. J Appl Ecol. 2016; https://doi.org/10.1111/1365-2664.12757.
Lueschow SR. Effect of Pseudogymnoascus destructans on microbial community composition on bats. Ph.D. diss., Western Illinois University. 2015. http://search.proquest.com/docview/1758623705?accountid=12543. Accessed 30 Dec 2016.
Romano-Bertrand S, Licznar-Fajardo P, Parer S, Jumas-Bilak E. Impact de l'environnement sur les microbiotes: focus sur l'hospitalisation et les microbiotes cutanés et chirurgicaux. Revue Francophone des Laboratoires. 2015;2015(469):75–82. https://doi.org/10.1016/S1773-035X(15)72824-8.
Avena CV, Parfrey LW, Leff JW, Archer HM, Frick WF, Langwig KE, et al. Deconstructing the bat skin microbiome: influences of the host and the environment. Front Microbiol. 2016;7:1753. https://doi.org/10.3389/fmicb.2016.01753.
Lemieux-Labonté V, Tromas N, Shapiro BJ, Lapointe FJ. Environment and host species shape the skin microbiome of captive neotropical bats. PeerJ. 2016;4:e2430. https://doi.org/10.7717/peerj.2430.
Winter AS, Kimble JC, Young JM, Buecher DC, Valdez EW, Hathaway JJM, et al. External bacterial diversity on bats in the southwestern United States: changes in bacterial community structure above and below ground. PeerJ Preprints. 2016;4:e2493v1. https://doi.org/10.7287/peerj.preprints.2493v1.
Kooser A, Kimble JC, Young JM, Buecher DC, Valdez EW, Porras-Alfaro A, et al. External microbiota of western United States bats: does it matter where you are from? bioRxiv. 2015; 017319. doi: https://doi.org/10.1101/017319.
Loudon AH, Woodhams DC, Parfrey LW, Archer H, Knight R, McKenzie V, et al. Microbial community dynamics and effect of environmental microbial reservoirs on red-backed salamanders (Plethodon cinereus). ISME J. 2014;8(4):830–40. https://doi.org/10.1038/ismej.2013.200.
Walke JB, Becker MH, Loftus SC, House LL, Cormier G, Jensen RV, et al. Amphibian skin may select for rare environmental microbes. ISME J. 2014;8(11):2207–17. https://doi.org/10.1038/ismej.2014.77.
Cornelison CT, Keel MK, Gabriel KT, Barlament CK, Tucker TA, Pierce GE, et al. A preliminary report on the contact-independent antagonism of Pseudogymnoascus destructans by Rhodococcus rhodochrous strain DAP96253. BMC Microbiol. 2014;14:246. https://doi.org/10.1186/s12866-014-0246-y.
Glassing A, Dowd SE, Galandiuk S, Davis B, Chiodini RJ. Inherent bacterial DNA contamination of extraction and sequencing reagents may affect interpretation of microbiota in low bacterial biomass samples. Gut Pathog. 2016;8:24. https://doi.org/10.1186/s13099-016-0103-7.
Preheim SP, Perrotta AR, Friedman J, Smilie C, Brito I, Smith MB, et al. Computational methods for high-throughput comparative analyses of natural microbial communities. Methods Enzymol. 2013;531:353–70. https://doi.org/10.1016/B978-0-12-407863-5.00018-6.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011;108:4516–22. https://doi.org/10.1073/pnas.1000080107.
Preheim SP, Perrotta AR, Martin-Platero AM, Gupta A, Alm EJ. Distribution-based clustering: using ecology to refine the operational taxonomic unit. Appl Environ Microbiol. 2013;79(21):6593–603. https://doi.org/10.1128/AEM.00342-13.
Alm Laboratories. SmileTrain. https://github.com/almlab/SmileTrain/wiki/. Accessed 1 July 2016.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1. https://doi.org/10.1093/bioinformatics/btq461.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6. https://doi.org/10.1038/nmeth.f.303.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72(7):5069–72. https://doi.org/10.1128/AEM.03006-05.
Shannon CE. A mathematical theory of communications. Bell Syst Tech J. 1948;27(3):379–423, 623–56. https://doi.org/10.1002/j.1538-7305.1948.tb01338.x.
Haegeman B, Hamelin J, Moriarty J, Neal P, Dushoff J, Weitz JS. Robust estimation of microbial diversity in theory and in practice. ISME J. 2013;7(6):1092–101. https://doi.org/10.1038/ismej.2013.10.
R Developement Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2016.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: linear and nonlinear mixed effects models. R package version 3.1–131. 2017. https://CRAN.R-project.org/package=nlme. Accessed 19 July 2017.
Lozupone CA, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35. https://doi.org/10.1128/AEM.71.12.8228-8235.2005.
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85. https://doi.org/10.1128/AEM.01996-06.
Lozupone CA, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5(2):169–72. https://doi.org/10.1038/ismej.2010.133.
Weiss S, Xu ZZ, Peddada S, Amir A, Bittinger K, Gonzalez A, et al. Normalization and microbial differential abundance strategies depend upon data characteristics. Microbiome. 2017;5:27. https://doi.org/10.1186/s40168-017-0237-y.
McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217. https://doi.org/10.1371/journal.pone.0061217.
Gower JC. Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika. 1966;53(3/4):325–38. https://doi.org/10.2307/2333639.
Dray S, Dufour AB. The ade4 package: implementing the duality diagram for ecologists. J Stat Softw. 2007;22(4):1–20. 10.18637/jss.v022.i04.
Gower JC, Legendre P. Metric and Euclidean properties of dissimilarity coefficients. J Classif. 1986;3(1):5–48. https://doi.org/10.1007/BF01896809.
Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490. https://doi.org/10.1371/journal.pone.0009490.
Legendre P, Anderson MJ. Distance-based redundancy analysis: testing multispecies responses in multifactorial ecological experiments. Ecol Monogr. 1999;69(1):1–24. https://doi.org/10.2307/2657192.
Oksanen J, Blanchet G, Kindt R, Legendre P, Minchin PR, O'Hara RB, et al. vegan: Community Ecology Package. R package version 2.4–1. 2016. http://CRAN.R-project.org/package=vegan. Accessed 29 Oct 2016.
Davies PT, Tso MKS. Procedures for reduced-rank regression. Appl Stat. 1982;31(3):244–55. https://doi.org/10.2307/2347998.
Ezekiel M. Methods of correlation analysis. New York and London: Wiley; 1930.
Dufrêne M, Legendre P. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecol Monogr. 1997;67(3):345–66. https://doi.org/10.2307/2963459.
De Cáceres M, Legendre P. Associations between species and groups of sites: indices and statistical inference. Ecology. 2009;90(12):3566–74. https://doi.org/10.1890/08-1823.1.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Statist. 1979;6(2):65–70. http://www.jstor.org/stable/4615733.
Voigt CC, Caspers B, Speck S. Bats, bacteria, and bat smell: sex-specific diversity of microbes in a sexually selected scent organ. J Mammal. 2005; 86(4):745-749. doi: http://dx.doi.org/10.1644/1545-1542(2005)086[0745:BBABSS]2.0.CO;2.
Chiba H, Agematu H, Kaneto R, Terasawa T, Sakai K, Dobashi K, et al. Rhodopeptins (Mer-N1033), novel cyclic tetrapeptides with antifungal activity from Rhodococcus sp. I. Taxonomy, fermentation, isolation, physico-chemical properties and biological activities. J Antibiot. 1999;52(8):695–9.
Nakayama K, Kawato HC, Inagaki H, Nakajima R, Kitamura A, Someya K, et al. Synthesis and antifungal activity of rhodopeptin analogues. 2. Modification of the west amino acid moiety. Org Lett. 2000;2(7):977–80. https://doi.org/10.1021/ol005630k.
Lauer A, Simon MA, Banning JL, Lam BA, Harris RN. Diversity of cutaneous bacteria with antifungal activity isolated from female four-toed salamanders. ISME J. 2008;2(2):145–57. https://doi.org/10.1038/ismej.2007.110.
Liu CH, Chen X, Liu TT, Lian B, Gu Y, Caer V, et al. Study of the antifungal activity of Acinetobacter baumannii LCH001 in vitro and identification of its antifungal components. Appl Microbiol Biotechnol. 2007;76(2):459–66. https://doi.org/10.1007/s00253-007-1010-0.
Song SJ, Lauber C, Costello EK, Lozupone CA, Humphrey G, Berg-Lyons D, et al. Cohabiting family members share microbiota with one another and with their dogs. elife. 2013;2:e00458. https://doi.org/10.7554/eLife.00458.
We thank Pascal Samson, Jocelyn Caron, Valérie Simard, Ariane Massé, Guillaume Tremblay, Kaleigh Norquay, Ana Breit, Drew Sippell, Manon Gagné, and Armand Yargeau for their assistance with the sampling, Julie Marleau for help with the experiments, Nicolas Tromas for help with the sequence treatment, and Quinn Fletcher for assistance in the data analysis.
The following reagent was obtained through BEI Resources, NIAID, NIH as part of the Human Microbiome Project: Genomic DNA from Microbial Mock Community B (Even, Low Concentration), v5.1L, for 16S rRNA Gene Sequencing, HM-782D.
This research was funded by Discovery Grants from the Natural Sciences and Engineering Research Council (NSERC Canada) OGP0155251 to FJL and RGPIN-2015-04437 to CKRW and excellence awards from the Quebec Centre for Biodiversity Science. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
The datasets generated and analyzed during the current study are available in the Figshare repository. Raw sequence data and metadata are available at DOI: https://figshare.com/s/623a1e47b4bed20459a7 and DOI: https://figshare.com/s/74d9497a792f9c0c76df.
Département de Sciences Biologiques, Université de Montréal, CP 6182, Succursale Centre-ville, Montréal, Québec, H2V 2S9, Canada
Virginie Lemieux-Labonté
& François-Joseph Lapointe
Direction de l'expertise sur la faune terrestre, l'herpétofaune et l'avifaune, Ministère des Forêts, de la Faune et des Parcs, Québec, Canada
Anouk Simard
Department of Biology and Centre for Forest Interdisciplinary Research, University of Winnipeg, Winnipeg, Manitoba, Canada
Craig K. R. Willis
Quebec Centre for Biodiversity Science, CP 6182, Succursale Centre-ville, Montréal, Québec, H2V 2S9, Canada
, Anouk Simard
Search for Virginie Lemieux-Labonté in:
Search for Anouk Simard in:
Search for Craig K. R. Willis in:
Search for François-Joseph Lapointe in:
VLL designed the experiments, performed the sampling, carried out the experiments, analyzed the data, and wrote the manuscript. CKRW and AS helped in the sampling logistic and wrote the manuscript. FJL designed the experiments and wrote the manuscript. All authors read and approved the final manuscript.
Correspondence to Virginie Lemieux-Labonté.
All methods were approved by the Animal Welfare and Ethics Committee at Université de Montréal (Protocol Number #16-015) and the University of Winnipeg Animal Care Committee (Protocol Number AEO5639).
Clusters of bats sampled within each hibernaculum and coded as dummy variables for db-RDA analysis.
Positive control mock community analysis. Sequence set comparisons of the mock community to what is expected. File 1 shows the matching sequences and related taxa identified. File 2 shows the taxa composition of the mock in relative abundance, the matching taxa at the genus or family level, and the false positive taxa.
Main taxa relative abundance (> 0.1%) in negative control samples. File 1 presents DNA extraction control samples, file 2 the library control, file 3 the negative site controls, and file 4 presents all controls together.
Rarefaction curves of alpha diversity calculated on multiple rarefied data table for each of the 66 bat skin microbiota samples and 11 environmental samples. (A) Shannon diversity of bat skin samples. (B) Overall richness (OTUs observed) of bat skin samples. (C) Shannon diversity of environmental samples. (D) Overall richness (OTUs observed) of environmental samples.
Main phyla identified in bat skin microbiota samples. The 8 more abundant phyla across all hibernacula are provided.
Main classes identified in bat skin microbiota samples. The 6 more abundant classes across all hibernacula are provided.
Major bacterial taxa identified in bat skin microbiota samples. The 16 more abundant taxa across all hibernacula are provided. Stars represent significant indicator taxa. *IndVal < 0.50, **IndVal ≥ 0.50, ***IndVal ≥ 0.89.
M. lucifugus skin microbiota taxa indicator test and related association measure (A, B) of six hibernaculum groups with different WNS status in Canada. Indicator value tests were computed with the multipatt() function of the indicspecies package in R. Only taxa with A ≥ 0.4 were retained as indicators. A, the specificity, is the probability that a site belongs to the group given the fact that the species is found and B, the fidelity, is the probability of finding a given taxon when the sites belong to that group. *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.
M. lucifugus skin microbiota taxa indicator and related association measure (A, B) of WNS-positive (Québec) and WNS-negative (Manitoba) sites in Canada. Indicator value tests were computed with the multipatt() function of the indicspecies package in R. Only taxa with A ≥ 0.4 were retained as indicators. A, the specificity, is the probability that a site belongs to the group given the fact that the species is found and B, the fidelity, is the probability of finding a given taxon when the sites belong to that group. *p ≤ 0.05, **p ≤ 0.01, ***p ≤ 0.001.
OTU table resulting from the analysis of 66 bat skin microbiota samples and 11 environmental samples.
https://doi.org/10.1186/s40168-017-0334-y
Skin microbiota
Myotis lucifugus
Conservation and management
16S rRNA


Good gift ideas for a males are demanding to locate. You have take into account a lot of things like his interest and his hobbies. You also need to make sure which he will enjoy and treasure it. Are less expensive to suit your budget too. A gift one particular special to be able to show appreciation for one person. Try to finish your finest in giving info about the subject one. There are some things that you can remember a person personally to capacity to offer best gift for a man close for.
There are different forms of concerts actually also be discovered on such sites. A lot of time there are outdoor concerts as well as runners who are held inside stadiums and circles. Most of the time these outdoor concerts are accessible when the next wind storm is nice out like during the spring and summer months, but there is the odd time that there will be one held during also of last year.
You may have heard actor James Woods in the number for the movies he's been in, most recently starring being a lawyer a TV show Shark. He is the involving person that can a evaluate the charges on his visa card at the end of the month and observed that there were thousands of dollars in charges he or she had not put on the card. On the list of items purchased were established very pricey VIP Concert Tickets for hamilton musical.
It's an indisputable incontrovertible fact that your response and your profits from whatever a person promoting will multiply anyone promote it via Residence website or splash fan page.
Brown the natural dancer since the age of more. But as a Vocalist, Brown didn't realize his talent until much later as he turn ten. Before then Brown aspired to develop into a rapper.
He tries on the mask and it takes over him. He decides to go back to the party, on the other hand without creating a lot of noise for that landlord incredibly. He then goes to auto shop and destroys area that it hurts. He wakes up your next morning without remembering a specific thing. A cop comes by and asks him about the disturbance that night. They know it was him so he throws the mask out belonging to the window. It finds its way right back.
They can contact the producers as a legitimate supply of tickets. Certainly they can get sales and commission, and voila; immediate cash. The more tickets sold, better chances they've got for an encouraging profit.
This innovative software a person to to speak while the computer types all the you say making just about of the actual 10 best Christmas gifts 2010 for anyone that hate typing as well as the related typos from typing are eliminated using your pc.


Please Note: This document on the program Rmatrix is accompanied by some notes on the basic theory of waves in layered media. Familiarity with the equations will assist the user with notation, which is rather similar in the FORTRAN code. Realistically, the theory needs to be understood in order to use the program correctly. The code is written with the same cartesian coordinates (x1,x2,z) or (x,y,z) in mind, with z positive downwards. SI units are used.
These notes are not exhaustive and they will grow and become more useful as users' questions are answered.
Rmatrix is a subroutine which returns the reflection and transmission matrices RU, RD, TU and TD for a stack of arbitrarily anisotropic homogeneous plane layers (see theory notes). These matrices embody everything about waves in the stack and may be used to model all kinds of propagation problems. For example, if many thin layers are used to simulate a gradient one may study back-scattered waves from the gradient or quasi-shear wave forward coupling. One can model plane waves incident from below a crustal/upper-mantle structure (e.g. the SKS and receiver function techniques). Surface waves in anisotropic media can be explored.
and for the given slownesses and frequency it returns the 3x3 reflection/transmission matrices TU, RU, TD and RD for a stack of anisotropic layers between two half spaces. The variables ifree and iter are the least obvious. Normally ifree is zero and then the values of TU and RD are those just above the first interface in the model (i.e. at the bottom of the first layer, so the thickness of that layer is irrelevant). When ifree is nonzero the phase shift across the top layer is included just before rmatrix is exited (useful when free-surface effects are to be included at the top of the first layer). The bottom layer of the model is taken as a lower half space, so its thickness is always irrelevant. When iter>1, the frequency-independent R/T coffecients at each single interface are not recalculated, but are assumed to have been stored for the given (px,py) in a previous call to rmatrix with iter=1.
Once the R/T matrices are found one may use them to solve various initial/boundary-value problems. For example, the propagator matrix for displacement/stress vectors can be constructed from equations (4.2) in the theory notes. Subroutines to add, multiply and invert 3x3 matrices are part of the program and these can be exploited.
Important point: Care should be taken with the ordering of the eigenvalues and eigenvectors. In the isotropic case the six columns are easily ordered into P-down, SV-down, SH-down, P-up, SV-up and SH-up (see subroutine isotroc). In anisotropic media the polarisations can sometimes behave in a complicated way, so little effort has been expended yet to provide a general sorting routine. The ordering currently depends on which (numerical) method is used to find the vertical slownesses (i.e. eigenvalues -- see comments in subroutine evset). If the most general routine is used (aniso6cg) or Georg Ruempker's routines (called from aniso3c) are used, the vertical slownesses are ordered in decreasing size (of real or imaginary part as appropriate), so the columns of the eigenvector matrix are possibly like P-down, S1-down, S2-down, S2-up, S1-up and P-up. If Sean Guest's routines are used (via aniso6c1) the columns are ordered more like that in the isotropic case. The user can order the columns as desired, so long as the first three correspond to downward propagating or decaying waves and the last three upward propagating or decaying waves. Awareness of the eigenvector orderings is important for the correct interpretation of the R/T coefficents, ray expansions, etc., and for anisotropic media the investigator will no doubt be looking closely at these quantities and the polarisations.
The program is divided up into 8 files. The first file (rmatrix1.f) has a main program, which calls rmatrix, two input subroutines and rmatrix itself. The second file (rmatrix2.f) contains the core routines of the reflection matrix calculation and its recursion scheme. The next 4 files (isotroc.f, aniso6cg.f, aniso6c1.f & aniso3c.f) contain subroutines to find the eigenvalues and vectors for a given layer. The user may supply personal versions of these if desired (e.g. efficient ones specific to, say, orthorhombic media). Finally there are two files (cg.f & deigv.f) containing public-domain matrix eigen-problem software, which can also be replaced at will by the user.
The travel time across the uppermost layer is included if ifree=1.
The first of these lines has the layer type index (2...fluid, 1...isotropic, 0...anisotropic), the thickness and the density. The second line has the P and S wavespeeds. For an anisotropic layer the first line the same followed by the 81 elements of aijkl=cijkl/rho (where the index i varies most rapidly, then j and so on). These indices refer to the cartesian coordinates, so 3 is definitely depth, for example. If a layer has a particular orientation of, say, olivine, then the corresponding rotated elastic constants must be provided to Rmatrix. Note the elements of aijkl are the elastic parameters cijkl divided by density. Subdirectory [inputs] at the ftp site (see below) contains some example input files (e.g. tranverse isotropy, silicon, sapphire, the Gutenberg model).
contains the eigenvalues and vectors in the layers for the current value of slowness.
Only one output file is produced (rmatrix.out) and by changing various ``iwrite'' variables throughout the program greater or lesser amounts of information can be dumped (e.g. eigenvalues and vectors in the layers).
There are three types of symmetry which can be tested. This is best discussed once the theory is understood.
The program can be obtained via anonymous ftp and information on how to get it can be obtained by sending email to thomson@geol.queensu.ca.
The paper is available in PostScript (77 kB) and GZIPped PostScript (21 kB).
In: Seismic Waves in Complex 3-D Structures, Report 7, pp. 163-167, Dep. Geophys., Charles Univ., Prague, 1998.


text:
"I Was There!" First-Person Accounts of the World
What is Witnify?
Witnify » Garfield
Creator Jim Davis Demonstrates How He Draws Garfield
By Phoebe Goldenberg • June 11, 2014 Uncategorized
Share on Facebook Tweet Add My Story
Jim Davis, creator of the "Garfield" comic strip, reveals how he brings its starring orange tabby cat to life on the page from conception to illustration. Since its debut in 1978, "Garfield" has become a pop culture sensation and currently holds the Guinness World Record for the world's most widely syndicated comic strip.
Submitted by: Phoebe Goldenberg
More Garfield Stories - Total:
Tagged 1970s, Animation, Art, Arts and Culture, Garfield, Jim Davis, TV.
« Son Discusses the Execution of Julius and Ethel Rosenberg
'Not Labeling Anything as Educational. It's All Entertainment.' »
Garfield is a comic strip created by Jim Davis. Published since June 19, 1978, it chronicles the life of the title character, the cat Garfield (named after the grandfather of Davis); his owner, Jon Arbuckle; and Jon's dog, Odie.
Best of Witnify
1940s 1950s 1960s 1970s 1980s 1990s 2000s 2010s Acting Arts and Culture Assassination Assassination of John F. Kennedy Athlete Baseball black history Business Championship Civil Rights Crime Film Football Gun Violence Health historical footage Hockey Invasion of Normandy Mass Tragedy Military Nature People Politics and History Protest rock music Science and Technology Space Speech Sports Terrorism The March on Washington TV U.S. Veteran Women World War II
Tweets by @Witnify
I agree to Witnify Terms of Service
Login Register Lost your password?
Not Logged in, login here

redpajama_set_name: RedPajamaCommonCrawl


text:
Gambling, Tourism Take Spotlight Ahead Of Planned Reopening Of Atlantic City Casinos
The shutdown of Atlantic City casinos due to the COVID-19 pandemic has already made an impact on the state.
AC's nine casinos are forced to rely on just NJ online gambling to bring in any bit of income.
While online play has grown exponentially during the past three months, the casinos have seen an 86% year-over-year total income drop.
And it's not just the casinos that are feeling the pain from this shutdown. The state of New Jersey as a whole is taking a hit from the steep loss of tax revenue.
But there are plenty of reasons to worry about the future of entertainment and tourism in the state, too.
And the big question boils down to safety.
How can Atlantic City, the hub of tourism in the state, stay safe enough so that people want to travel there in the future?
NJ Assembly panel discusses tourism and COVID-19
To that end, a panel discussion was held in Trenton on June 25 where all of these concerns were addressed.
Casino representatives along with people from the summer rental, live entertainment, and tourism industries met to discuss the impact that the coronavirus pandemic has had so far and how they can keep visitors safe.
Assemblyman Ralph Caputo, chair of the state assembly Tourism, Gaming and the Arts Committee, shared the concerns that were raised during the meeting in a statement:
"Tourism and hospitality are vital industries in New Jersey, bringing in over 100 million visitors every year, producing (billions) of dollars in revenue and providing nearly 10 percent of the jobs throughout our state."
Luckily for Atlantic City, Gov. Phil Murphy announced that casinos can reopen at 25% capacity on July 2. All patrons will have to wear a mask and screenings will be done as well.
This news was followed on Monday by the governor canceling plans to reopen indoor dining on July 2.
Caputo goes on to say that "gaming and the arts have made significant contributions to our state's culture and economy." The state has a wide range of beaches, casinos, cuisine, and more, but "all of these have been greatly impacted by COVID-19."
'Successful steps' to keep casinos safe
New Jersey has been one of the more cautious states in regards to reopening places due to the pandemic.
State officials are not messing around about maintaining the downward trend of cases, even going as far as suing a city to prevent them from having indoor dining before it was approved.
Murphy shared some of those concerns during a radio station appearance.
"We're spending a lot of time with the owners and operators as well as the workers and their representatives on casinos," said Murphy to WOND's AC Mike Lopez.
"The bad news is casinos sort of have the attributes that are hardest to deal with this virus. It's indoors, no ventilation. You're sedentary, you're in close proximity. Having said that, they're big footprints and we think there are successful steps that can be taken to address the challenges."
State officials probably also saw what happened when Las Vegas opened up at the start of June without a mask order in place.
There was a sharp uptick in new COVID-19 cases. A required mask order went into effect last week.
AC casinos put safety plans in place
Some of Atlantic City's casinos have already come forward with plans to ensure safety and cleanliness for visitors.
For example, Hard Rock Casino plans to use thermal cameras to check guests' temperature. The casino also will space out slot machines and install plexiglass barriers for table games. All AC casinos plan to implement similar procedures.
Once Hard Rock reopens on July 2, it will also have a new protocol called the Safe + Sound Clean Team. The team consists of 100 employees whose sole job is to deep clean all gaming surfaces, chips, hotel rooms, restaurants, pools, and public areas.
Casinos will also be providing hand sanitizer stations on the casino floor along with enforcing the state's mandates regarding masks and social distancing.
But these plans are all dependent on what guidelines the state and Division of Gaming Enforcement require of casinos, too. Those have yet to be made public.
Cashless gambling in Atlantic City?
Another, more innovative tactic that Atlantic City's casinos could possibly employ is cashless gaming.
Just approved by the Nevada Gaming Commission this past week, cashless gaming would be not only a great way to help prevent COVID-19 spread but also gives players more options.
The new regulations approved by Nevada would allow players to withdraw or deposit money from either table games and/or slot machines with a digital wallet such as PayPal, Venmo, or Apple Pay.
This new way of "money handling" would drastically cut back on the use of chips by players, casino cashiers, and dealers alike.
Implementing cashless gambling, however, is only one part of any plan to make Atlantic City safe for tourists. Safety plans and protocols as well as following state guidelines are essential to AC's gambling future.
The fact remains, though, that the flow of tax revenue into the state's coffers shrank while casinos were closed. Failing to bring Atlantic City back safely will deter visitors.
And if the pandemic worsens in NJ, casinos will likely shut down once more.


text:
Boomtimes Again: Twentieth-Century Mining in the Mojave Desert
Spanish conquistadors sought gold in the desert of what is now the Southwest US, and the prospectors who hurried to California in the post-1849 gold rush also mined this area. However, the Mojave desert truly boomed as a destination for miners when the price of silver escalated quickly in the 1890s. Boomtowns such as Tonopah, Nevada, in the eastern Mojave, became the nation's prime source of this precious metal. Some of these boomtowns, such as Ballarat, California, became ghost towns when mines closed after the panic of 1907, but others thrived through the 1930s. In addition to silver and copper, vast deposits of minerals useful to the nation's burgeoning industries were discovered and mined throughout the Mojave Desert of California and Nevada in the late nineteenth and early twentieth century. Large reserves of potash, soda ash, sodium sulfate, salt, and borax found in the earth around the dry Searles Lake gave rise to the twentieth-century boomtown of Trona, California. Rail lines were constructed through the desert to transport goods to markets and industrial centers around the US, as well as to ports where minerals could be shipped around the world.
A second, mid-twentieth-century boom period was fueled by the need for raw materials during the World Wars and the Cold War. Mining in these years was "big business" rather than the work of smaller, independent prospectors. Post-1920 mining in the Mojave Desert typically saw corporations running mines and processing centers, and company-run towns developing around these industrial centers.
Show full overview
Kerry Dunne, Weston Public Schools, Massachusetts
The Emergence of Modern America (1890-1930)
The Great Depression and World War II (1929-1945)
Cite this set
Chicago citation style
Kerry Dunne. Boomtimes Again: Twentieth-Century Mining in the Mojave Desert. 2016. Retrieved from the Digital Public Library of America, http://dp.la/primary-source-sets/boomtimes-again-twentieth-century-mining-in-the-mojave-desert?subject=us-history&timePeriod=the-emergence-of-modern-america-1890-1930. (Accessed January 17, 2020.)
Kerry Dunne, (2016) Boomtimes Again: Twentieth-Century Mining in the Mojave Desert. Retrieved from the Digital Public Library of America, http://dp.la/primary-source-sets/boomtimes-again-twentieth-century-mining-in-the-mojave-desert?subject=us-history&timePeriod=the-emergence-of-modern-america-1890-1930
MLA citation style
Kerry Dunne. Retrieved from the Digital Public Library of America <http://dp.la/primary-source-sets/boomtimes-again-twentieth-century-mining-in-the-mojave-desert?subject=us-history&timePeriod=the-emergence-of-modern-america-1890-1930>.
Source Set
A 1905 photograph of Tonopah, Nevada.
A photograph of Tonopah, Nevada ca. 1915.
A photograph of late-nineteenth-century silver mining in Tonopah, Nevada.
A photograph of people and mining equipment in early-twentieth-century Tonopah, Nevada.
A postcard image of downtown Tonopah, Nevada ca. 1907.
A 1912 map of mining in Tonopah, Nevada.
A page from a 1912 newspaper, the Tonopah Miner.
A 1906 document listing taxpayers in Tonopah, Nevada.
A 1957 letter from a former mining employee to a friend in Trona, California, describing memories of working in Goldfield, Nevada.
A photograph of a twenty-mule team used to transport metals and minerals prior to the development of railroads in the Mojave Desert.
A 1925 photograph of magnesium bagged for transport by train near Death Valley, California.
An excerpt from a 1915 booklet called "Borax the Magic Crystal."
A 1946 photo of a mineral processing plant in Trona, California.
A 1989 photo of dry lakes used for mineral harvesting in Trona, California.
A 1989 photo of the Kerr McGee plant in Trona, California.
Related Primary Source Sets
The Scopes Trial
Mexican Labor and World War II: The Bracero Program
Environmental Preservation in the Progressive Era
When Miners Strike: West Virginia Coal Mining and Labor History
Eugenics Movement in the United States
Reservations, Resistance, and the Indian Reorganization Act, 1900-1940
Full Steam Ahead: The Steam Engine and Transportation in the Nineteenth Century
Frank Lloyd Wright and Modern American Architecture
Theodore Dreiser's Sister Carrie and the Urbanization of Chicago
Fake News in the 1890s: Yellow Journalism
The Columbian Exchange
Truth, Justice, and the Birth of the Superhero Comic Book
Blackface Minstrelsy in Modern America
The Rise of Italian Fascism and Its Influence on Europe
The United Farm Workers and the Delano Grape Strike
The Invention of the Telephone
Japanese American Internment During World War II
The Yellow Fever Epidemic of 1878
Settlement Houses in the Progressive Era
Ida B. Wells and Anti-Lynching Activism
Second Ku Klux Klan and The Birth of a Nation
World War I: America Heads to War
Beginnings of the American Red Cross
There is No Cure for Polio
Feeding the Hungry with Food Stamp Programs
Immigration through Angel Island
African American Soldiers in World War I
The Homestead Strike
The New Woman
The Populist Movement
Attacks on American Soil: Pearl Harbor and September 11
Electrifying America
The New Deal
Women's Suffrage: Campaign for the Nineteenth Amendment
American Imperialism: The Spanish-American War
Cotton Gin and the Expansion of Slavery
World War II: Women on the Home Front
Treaty of Versailles and the End of World War I
The Panama Canal
The Atomic Bomb and the Nuclear Age


text:
Arvo Valton (właśc. Arvo Vallikivi, ur. 14 grudnia 1935 w Märjamaa) – estoński pisarz, poeta, scenarzysta i krytyk literacki.

Arvo Valton był synem przedsiębiorcy w północno-zachodniej Estonii. W latach 1943-1949 uczęszczał do szkoły podstawowej w Märjamaa (prowincja Rapla). Zesłany wraz z rodziną na Syberię, do obwodu nowosybirskiego, skąd do rodzinnej Estonii mógł wrócić dopiero w 1954 r. W latach 1954-1959 studiował na politechnice w Tallinnie, na wydziale chemii i górnictwa. Uczelnię ukończył z tytułem inżyniera górnictwa. Później przez dwa lata (1959-1961) pracował w kombinacie chemicznym w Maardu (północna Estonia), a następnie do 1968 r. w wytwórni urządzeń pomiarowych w Tallinnie.

W latach 60. w Moskwie zaocznie ukończył kurs scenopisania. W 1965 r. wstąpił do Związku Pisarzy Estońskich, a w latach 1989-1992 pełnił w nim funkcję przewodniczącego. Między 1968 a 1975 r. utrzymywał się z pisania. Następnie zatrudniony był w państwowym przedsiębiorstwie filmowym Tallinnfilm. W 1992 r., krótko po odzyskaniu przez Estonię niepodległości, został posłem do parlamentu estońskiego (Riigikogu). W 1995 r. wycofał się z polityki.

Twórczość
Arvo Valton jest autorem licznych opowiadań, nowel i powieści. Zajmował się również krytyką literacką, pisał wiersze, aforyzmy, reportaże z podróży oraz sztuki teatralne. Zadebiutował w 1960 r. na łamach prasy. Jest autorem m.in. scenariusza do jednego z najbardziej popularnych w Estonii filmów, pt. Viimne reliikvia (1969). Napisał również libretto do opery Lend Eino Tamberga, która została wystawiona po raz pierwszy w Tallinnie w 1980 r. Część jego twórczości, która w krytyczny sposób odnosi się do rzeczywistości radzieckiej lub pisana była pod wpływem niechętnie widzianego w ZSRR egzystencjalizmu, nie mogła się ukazać drukiem ze względu na ingerencję cenzury. W latach 80. podnosił kwestie zniszczeń, wyrządzanych przez przemysł ciężki środowisku naturalnemu; angażował się w działania na rzecz niepodległości państwa estońskiego, wspierał również inicjatywy podejmujące kwestię praw ludów i narodów ugrofińskich w ZSRR.

Dwukrotnie został wyróżniony nagrodą im. Friedeberta Tuglasa: w 1973 r. za opowiadanie Ohtlik leiutis oraz w 1979 r. za opowiadanie Mustamäe armastus, opublikowane w tomie opowiadań pod tym samym tytułem.

W przekładzie na język polski ukazał się tom opowiadań Valtona, zatytułowany Spiritual w wyborze (z różnych tomów) i opracowaniu Aarne Puu oraz w przekładzie Aarne Puu i Walerii Pawłowicz, Kraków 1988. Prócz tego tłumaczenie dziesięciu opowiadań Valtona (w tym wyróżniona nagrodą im. Frideberta Tuglasa Miłość w Mustamäe), znalazło się w tomie Kochanka diabła. Opowiadania estońskie, opracowanym i przełożonym przez Aarne Puu, Warszawa 1984. Opowiadanie Beczka ukazało się na łamach czasopisma Literatura radziecka, nr 1/1989 w przekładzie E. Hessen.

Dzieła
 1963 - Veider soov (opowiadania)
 1965 - Rataste vahel (opowiadania)
 1968 - Kaheksa jaapanlannat
 1972 - Sõnumitooja (opowiadania)
 1972 - Õukondlik mäng (opowiadania)
 1974 - Pööriöö külaskäik (opowiadania)
 1975 - Läbi unemaastike (opowiadania)
 1977 - Uksed kriuksuvad öösiti (aforyzmy)
 1978 - Tee lõpmatuse teise otsa (powieść)
 1978 - Mustamäe armastus (opowiadania)
 1978 - Õndsusse kulgev päev (opowiadanie)
 1980 - Märklaud kilbiks (aforyzmy)
 1981 - Ajaprintsess (książka dla dzieci)
 1981 - Võõras linnas (opowiadania)
 1982 - Loomise iga (powieść)
 1983 - Põhjanaela paine (powieść)
 1984 - Arvid Silberi maailmareis (powieść surrealistyczna)
 1983-1985 - Üksildased ajas (opowiadania w dwóch tomach)
 1989 - Masendus ja lootus (powieść)
 1998 - Kogutud teosed (utwory zebrane)

Bibliografia
 Arvo Valton (notka biograficzna) [w:] Kochanka diabła. Opowiadania estońskie, wybór i przekład Aarne Puu, Warszawa 1984.

Linki zewnętrzne
 Krótki biogram Valtona na stronie estlit.ee

Przypisy

Estońscy pisarze
Estońscy poeci
Urodzeni w 1935
Estońskie ofiary prześladowań komunistycznych 1940–1991

redpajama_set_name: RedPajamaWikipedia


text:
require "rails51/version"

module Rails51
end

redpajama_set_name: RedPajamaGithub


The medicamental selector "Alfa-02" is designed for electropunctural diagnostics by the method of dr. R.Voll with the support of the leading experts of National Medical University and Kiev Medical Academy. According to information of WHO this method receives the wide recognition and is successfully used in 34 countries of the world, as well as in Ukraine and Russia (was recommended for using in medical practice by the Ministry of the Health Protection of USSR in decree N211 from 6.06.89). "For the deserts to suffering humanity and the elaboration of the electroacupuncture diagnostics and therapy method..." dr. R. Voll was rewarded by gold medal of Vatican and by medal of Gufaland.
The condition of organs of the digestion, respiration, movement, circulation of the blood, secretion of urine; nervous, limphatic, endocrinic systems; liver and spleen - without using of X-ray, laboratories tests, ultrasound's examinations and other traditional methods - right up to revealing of ethiological factors of the illness even in early stages.
The acute conditions, where the surgical intervention is required, as sharpening of the appendicitis, sharpening of the chronic adnecsitis, sharpening of the chronic colitis, etc; the serious intestinal illnesses: dysentery, hepatitis, salmonelesis, etc.
The presence in the body of the radio-active elements, nitric fertilizers, pesticides, herbicides, salts of heavy metals and other kinds of industrial and ecological toxins.
Hidden and latent sources of infection and their influence on the body's reaction.
The availability or absence of the allergical or bioenergetic influence on the patient of some kinds of vegetable and animal foods-stuffs, food supplements, solvents, jewelry adornments, precious metals, etc.
The bioenergetic compatibility with the human body of dental materials and cosmetics.
Selector in base configuration contain about 2000 bioenergetic analogies of various preparations - preparations of body's organs and herbs, etiological, allopathic, homeopathic, neoplastic, pathomorphologic preparations, preparations of firms "Vecroma" and "Heel", of allergens, microelements, secretions, radio-active elements, dental materials, ecological toxins.
Therapy influence on the acupuncture points of the patient by the bioenergetic analogies of preparations and reproduce their properties and properties of any remedy on some substances - distilled water, wax and milk sugar. After the reproduction substances becomes biologically active and can be taken inside by certain scheme or used as application on certain acupuncture points for therapy (so named treatment without medicaments.
Verifying the efficiency of the therapy.


text:
PPCE Expands Partnership With Medidata
PPCE expands its partnership with Medidata by including Coder to its suite of services for Clinical Data Management
Pune. India -- Jan 8, 2018 - Phoenix Progressive Certifications Enterprise Pvt. Ltd. (PPCE), a provider of Clinical Data Management, Biostatistics and Medical Writing services for global pharmaceutical and medical device companies, announced today that it has furthered its partnership with Medidata (NASDAQ: MDSO), the leading global provider of cloud-based technology and data analytics for clinical research. The Company has accreditations to implement and market Medidata products, including Medidata Coder®, which streamlines and centralizes medical coding.
"In 2016, PPCE made the decision to partner with Medidata to provide our clients with study design and conduct services on Medidata Rave®, the world's leading solution for capturing, managing and reporting patient data. In October 2017, we completed another accreditation program for the Medical Coding module. Our plan now is to integrate Rave with the Medidata Centralized Statistical Analysis model to further improve the efficiency of the clinical development process by giving our clients the advantage of timely information to aid decision making. Not many companies have access to the sort of meaningful information which can facilitate real-time decisions on, altering the design, for example, or even prematurely terminating a trial if outcomes clearly cannot be achieved." said Khushnoor Dastoor, executive director & head QA at PPCE.
"We're thrilled to see the progress PPCE is making with our platform," said Glen de Vries, Medidata president, and co-founder. "Our partnership highlights how life sciences research can be supported around the world."
About PPCE:
Phoenix Progressive Certifications Enterprise Pvt. Ltd. (PPCE) was founded in 1997 with a vision to add value to the industry by preventing problems arising out of poor quality.
PPCE is a functional service provider of Clinical Data Management, Biostatistics and Medical Writing services for global pharmaceutical and medical device companies. Fully validated platforms and strong quality assurance systems enable PPCE to offer world-class tools for high-quality data resulting in high customer satisfaction and long-term relationships with its partners.
Known for their customized approach designed around clients' specific requirements, PPCE's teams have an in-depth knowledge of international regulations and flexible business models to suit different operational requirements. PPCE has experience in Oncology, Urology, Autoimmune Disorders, Respiratory, Cardiovascular and Immunology studies.
About Medidata:
Medidata is reinventing global drug and medical device development by creating the industry's leading cloud-based solutions for clinical research. Through our advanced applications and intelligent data analytics, Medidata helps advance the scientific goals of life sciences customers worldwide, including over 950 global pharmaceutical companies, biotech, diagnostic and device firms, leading academic medical centers, and contract research organizations.
The Medidata Clinical Cloud® brings a new level of quality and efficiency to clinical trials that empower our customers to make more informed decisions earlier and faster. Our unparalleled clinical trial data assets provide deep insights that pave the way for future growth. The Medidata Clinical Cloud is the primary technology solution powering clinical trials for 18 of the world's top 25 global pharmaceutical companies and is used by 18 of the top 25 medical device developers—from study design and planning through execution, management, and reporting.
Phoenix Progressive Certifications Enterprise Pvt. Ltd.
502, The Chambers, 4/12/13 Viman Nagar
Pune, Maharashtra-411014
Website: http://www.ppceworld.com