Using a network modularity analysis to inform management of a rare endemic plant in the northern Great Plains, USA

Microbial communities in tree root-compartment niches under Cd and Zn pollution: Structure, assembly process and co-occurrence relationship

Woody plants have showed great potential in remediating severely contaminated soils by heavy metals (HMs) due to their cost-efficient and ecologically friendly trait. It is believed the root-associated microbiota plays a vital role in phytoremediation for HMs. However, the ecological process controlling the assembly and composition of tree root-associated microbial communities under HMs stress remains poorly understood. Herein, we profiled the bulk soil, rhizosphere and endosphere microbial communities of trees growing in heavily Cd and Zn polluted soil. The microbiota was gradually filtered from bulk soil to the tree roots and was selectively enriched in roots with specific taxa, such as Proteobacteria and Ascomycota. The microbial community assembly along the soil-root continuum was mainly controlled by deterministic processes from bulk soil to the endosphere, with the normalized stochasticity ratio (NST) indices of 67.16-31.05 % and 30.37-15.02 % for bacteria and fungi, respectively. Plant selection pressure sequentially increased from bulk soil to rhizosphere to endosphere, with the reduced bacterial alpha diversity accompanying the consequently reduced complexity of the co-occurrence network. Together, the findings provide new evidence for horizontal transmission of endophytic microbiome from soil to the host, which can shed light on the future screening and application of microbial-assisted phytoremediation.

Edge disturbance shapes liana diversity and abundance but not liana-tree interaction network patterns in moist semi-deciduous forests, Ghana

Edge disturbance can drive liana community changes and alter liana-tree interaction networks, with ramifications for forest functioning. Understanding edge effects on liana community structure and liana-tree interactions is therefore essential for forest management and conservation. We evaluated the response patterns of liana community structure and liana-tree interaction structure to forest edge in two moist semi-deciduous forests in Ghana (Asenanyo and Suhuma Forest Reserves: AFR and SFR, respectively). Liana community structure and liana-tree interactions were assessed in 24 50 × 50 m randomly located plots in three forest sites (edge, interior and deep-interior) established at 0-50 m, 200 m and 400 m from edge. Edge effects positively and negatively influenced liana diversity in forest edges of AFR and SFR, respectively. There was a positive influence of edge disturbance on liana abundance in both forests. We observed anti-nested structure in all the liana-tree networks in AFR, while no nestedness was observed in the networks in SFR. The networks in both forests were less connected, and thus more modular and specialised than their null models. Many liana and tree species were specialised, with specialisation tending to be symmetrical. The plant species played different roles in relation to modularity. Most of the species acted as peripherals (specialists), with only a few species having structural importance to the networks. The latter species group consisted of connectors (generalists) and hubs (highly connected generalists). Some of the species showed consistency in their roles across the sites, while the roles of other species changed. Generally, liana species co-occurred randomly on tree species in all the forest sites, except edge site in AFR where lianas showed positive co-occurrence. Our findings deepen our understanding of the response of liana communities and liana-tree interactions to forest edge disturbance, which are useful for managing forest edge.

Patch utilization and flower visitations by wild bees in a honey bee-dominated, grassland landscape

Understanding habitat needs and patch utilization of wild and managed bees has been identified as a national research priority in the United States. We used occupancy models to investigate patterns of bee use across 1030 transects spanning a gradient of floral resource abundance and richness and distance from apiaries in the Prairie Pothole Region (PPR) of the United States. Estimates of transect use by honey bees were nearly 1.0 during our 3.5-month sampling period, suggesting honey bees were nearly ubiquitous across transects. Wild bees more frequently used transects with higher flower richness and more abundant flowers; however, the effect size of the native flower abundance covariate (β ^ native  = 3.90 ± 0.65 [1SE]) was four times greater than the non-native flower covariate (β ^ n o n - n a t i v e  = 0.99 ± 0.17). We found some evidence that wild bee use was lower at transects near commercial apiaries, but the effect size was imprecise (β ^ distance  = 1.4 ± 0.81). Honey bees were more frequently detected during sampling events with more non-native flowers and higher species richness but showed an uncertain relationship with native flower abundance. Of the 4039 honey bee and flower interactions, 85% occurred on non-native flowers, while only 43% of the 738 wild bee observations occurred on non-native flowers. Our study suggests wild bees and honey bees routinely use the same resource patches in the PPR but often visit different flowering plants. The greatest potential for resource overlap between honey bees and wild bees appears to be for non-native flowers in the PPR. Our results are valuable to natural resource managers tasked with supporting habitat for managed and wild pollinators in agroecosystems.

Coflowering invasive plants and a congener have neutral effects on fitness components of a rare endemic plant

Network analyses rarely include fitness components, such as germination, to tie invasive plants to population-level effects on the natives. We address this limitation in a previously studied network of flower visitors around a suite of native and invasive plants that includes an endemic plant at Badlands National Park, South Dakota, USA. Eriogonum visheri coflowers with two abundant invasive plants, Salsola tragus and Melilotus officinalis, as well as a common congener, E. pauciflorum. Network analyses had suggested strong linkages between E. visheri and S. tragus and E. pauciflorum, with a weaker link to M. officinalis. We measured visitation, pollen deposited on stigmas, achene weight and germination over three field seasons (two for germination) in four populations (two in the final season) of E. visheri and applied in situ pollen treatments to E. visheri, adding pollen from other flowers on the same plant; flowers on other E. visheri plants; S. tragus, M. officinalis, or E. pauciflorum; open pollination; or excluding pollinators. Insect visitation to E. visheri was not affected by floral abundance of any of the focal species. Most visitors were halictid bees; one of these (Lasioglossum packeri) was the only identified species to visit E. visheri all three years. Ninety-seven percent of pollen on collected E. visheri stigmas was conspecific, but 22% of flowers had >1 grain of E. pauciflorum pollen on stigmas and 7% had >1 grain of S. tragus pollen; <1% of flowers had M. officinalis pollen on stigmas. None of the pollen treatments produced significant differences in weight or germination of E. visheri achenes. We conclude that, in contrast to the results of the network analysis, neither of the invasive species poses a threat, via heterospecific pollen deposition, to pollination of the endemic E. visheri, and that its congener provides alternative pollen resources to its pollinators.

A Comparison of Honey Bee-Collected Pollen From Working Agricultural Lands Using Light Microscopy and ITS Metabarcoding

Taxonomic identification of pollen has historically been accomplished via light microscopy but requires specialized knowledge and reference collections, particularly when identification to lower taxonomic levels is necessary. Recently, next-generation sequencing technology has been used as a cost-effective alternative for identifying bee-collected pollen; however, this novel approach has not been tested on a spatially or temporally robust number of pollen samples. Here, we compare pollen identification results derived from light microscopy and DNA sequencing techniques with samples collected from honey bee colonies embedded within a gradient of intensive agricultural landscapes in the Northern Great Plains throughout the 2010-2011 growing seasons. We demonstrate that at all taxonomic levels, DNA sequencing was able to discern a greater number of taxa, and was particularly useful for the identification of infrequently detected species. Importantly, substantial phenological overlap did occur for commonly detected taxa using either technique, suggesting that DNA sequencing is an appropriate, and enhancing, substitutive technique for accurately capturing the breadth of bee-collected species of pollen present across agricultural landscapes. We also show that honey bees located in high and low intensity agricultural settings forage on dissimilar plants, though with overlap of the most abundantly collected pollen taxa. We highlight practical applications of utilizing sequencing technology, including addressing ecological issues surrounding land use, climate change, importance of taxa relative to abundance, and evaluating the impact of conservation program habitat enhancement efforts.

Exotic Plant Infestation Is Associated with Decreased Modularity and Increased Numbers of Connectors in Mixed-Grass Prairie Pollination Networks

The majority of pollinating insects are generalists whose lifetimes overlap flowering periods of many potentially suitable plant species. Such generality is instrumental in allowing exotic plant species to invade pollination networks. The particulars of how existing networks change in response to an invasive plant over the course of its phenology are not well characterized, but may shed light on the probability of long-term effects on plant-pollinator interactions and the stability of network structure. Here we describe changes in network topology and modular structure of infested and non-infested networks during the flowering season of the generalist non-native flowering plant, Cirsium arvense in mixed-grass prairie at Badlands National Park, South Dakota, USA. Objectives were to compare network-level effects of infestation as they propagate over the season in infested and non-infested (with respect to C. arvense) networks. We characterized plant-pollinator networks on 5 non-infested and 7 infested 1-ha plots during 4 sample periods that collectively covered the length of C. arvense flowering period. Two other abundantly-flowering invasive plants were present during this time: Melilotus officinalis had highly variable floral abundance in both C. arvense-infested and non-infested plots and Convolvulus arvensis, which occurred almost exclusively in infested plots and peaked early in the season. Modularity, including roles of individual species, and network topology were assessed for each sample period as well as in pooled infested and non-infested networks. Differences in modularity and network metrics between infested and non-infested networks were limited to the third and fourth sample periods, during flower senescence of C. arvense and the other invasive species; generality of pollinators rose concurrently, suggesting rewiring of the network and a lag effect of earlier floral abundance. Modularity was lower and number of connectors higher in infested networks, whether they were assessed in individual sample periods or pooled into infested and non-infested networks over the entire blooming period of C. arvense. Connectors typically did not reside within the same modules as C. arvense, suggesting that effects of the other invasive plants may also influence the modularity results, and that effects of infestation extend to co-flowering native plants. We conclude that the presence of abundantly flowering invasive species is associated with greater network stability due to decreased modularity, but whether this is advantageous for the associated native plant-pollinator communities depends on the nature of perturbations they experience.

Taxonomic Characterization of Honey Bee (Apis mellifera) Pollen Foraging Based on Non-Overlapping Paired-End Sequencing of Nuclear Ribosomal Loci

Identifying plant taxa that honey bees (Apis mellifera) forage upon is of great apicultural interest, but traditional methods are labor intensive and may lack resolution. Here we evaluate a high-throughput genetic barcoding approach to characterize trap-collected pollen from multiple North Dakota apiaries across multiple years. We used the Illumina MiSeq platform to generate sequence scaffolds from non-overlapping 300-bp paired-end sequencing reads of the ribosomal internal transcribed spacers (ITS). Full-length sequence scaffolds represented ~530 bp of ITS sequence after adapter trimming, drawn from the 5’ of ITS1 and the 3’ of ITS2, while skipping the uninformative 5.8S region. Operational taxonomic units (OTUs) were picked from scaffolds clustered at 97% identity, searched by BLAST against the nt database, and given taxonomic assignments using the paired-read lowest common ancestor approach. Taxonomic assignments and quantitative patterns were consistent with known plant distributions, phenology, and observational reports of pollen foraging, but revealed an unexpected contribution from non-crop graminoids and wetland plants. The mean number of plant species assignments per sample was 23.0 (+/- 5.5) and the mean species diversity (effective number of equally abundant species) was 3.3 (+/- 1.2). Bray-Curtis similarities showed good agreement among samples from the same apiary and sampling date. Rarefaction plots indicated that fewer than 50,000 reads are typically needed to characterize pollen samples of this complexity. Our results show that a pre-compiled, curated reference database is not essential for genus-level assignments, but species-level assignments are hindered by database gaps, reference length variation, and probable errors in the taxonomic assignment, requiring post-hoc evaluation. Although the effective per-sample yield achieved using custom MiSeq amplicon primers was less than the machine maximum, primarily due to lower “read2” quality, further protocol optimization and/or a modest reduction in multiplex scale should offset this difficulty. As small quantities of pollen are sufficient for amplification, our approach might be extendable to other questions or species for which large pollen samples are not available.