The count per cage of each variety was averaged to estimate survival and reproduction

Plots of UC 92 were vacuumed in 2019 and monitored with water traps in 2021 to check the difference in L. hesperus densities. In 2019, plots in the sprayed block had on average 0.5 L. hesperus while plots in the unsprayed block had on average 1.5 L. hesperus. In 2021, water traps collected every 3 days for 2 weeks after the insecticide treatment indicated that in UC 92 plots there were an average of 0.125 L. hesperus in the sprayed section and an average of 0.375 L. hesperus in the unsprayed section. These data indicate that there was a difference in L. hesperus densities between the two treatment blocks. High tolerance of damage by L. hesperus was defined by dividing the yield of untreated plots by the yield of insecticide treated plots . If the quotient was high, the variety was highly defended against L. hesperus. The most desirable varieties for breeding are those with high tolerance to L. hesperus and high yield. It should be noted that the term “tolerance” is used here as a place holder for the true mechanism of resilience. It is possible that some varieties may be resistant rather than tolerant. A resistant variety can avoid damage while a tolerant variety can recover from damage. Sixteen commercial cultivars from the United States were included in this study . On average, square pot these elite lines yielded 65% as much in conditions unprotected by insecticide as compared to conditions with protection of insecticide. Pat and Kingstone had the highest L. hesperus tolerance while UC Lee , UC 92 , and Fordhook had the lowest.

UC Haskell performed the best when both yield and L. hesperus tolerance are accounted for.Among the 56 heirloom and international cultivars, the average yield in unprotected conditions compared to protected conditions was 73%. Noir De Kisenyi – G26196 , a baby, black-and-white-seeded, bush type from Rwanda, and Hopi 50 , a baby, white-seeded, vine type from the United States, were the least susceptible to L. hesperus but did not compete well on overall yield. The American heirloom variety, Jackson Wonder, had the best combination of high L. hesperus tolerance and an average yield of 2338.4 g/plot when treated with insecticide and 1945.4 g/plot when unprotected by insecticide. Many of the international and heirloom varieties included in this study have already been incorporated into the UC Davis Lima Bean Breeding Program. These results will offer additional guidance in parent selection and the choice of varieties used in future research on insect defense mechanisms.When comparing total number of L. hesperus per plot and yield per plot, yield did not appear to be affected by the density of L. hesperus, as shown by the nearly horizontal trend line . In this group of varieties, UC 92 appeared to be an outlier with low yields but also a low number of L. hesperus. UC 92 is the only variety in this study from the Andean gene pool of Lima bean. While the yields of Andean and Mesoamerican varieties are comparable, Andean lines tend to have much larger seeds. As a result, the loss of a single seed due to L. hesperus herbivory would constitute a loss of a higher percentage of yield than when a single small seed is lost .It should also be noted that data from this experimental design may over represent resistance or tolerance traits. On a typical California farm, a single variety will be planted across several if not tens of acres. While L. hesperus adults are highly mobile, they have less choice in an actual farm setting than in this experimental field.

Nymphs were not considered in this study because they were crushed by the strong pressure of the vacuumand difficult to identify with the available expertise. Nymphs are much less mobile than adults and therefore would have less ability to choose between varieties. It is interesting to note the strong difference in the average number of L. hesperus per plot between the two highest yielding varieties, UC Haskell and UC Beija Flor. The large difference in average L. hesperus might be attributable to sampling bias caused by the denser canopy of the indeterminate UC Haskell compared to the determinate UC Beija Flor. It could also be that these two varieties have different mechanisms of resistance or tolerance to L. hesperus. UC Haskell and UC Beija Flor share as a parent UC Cariblanco N, but their resistance probably comes from their other respective parents. UC Beija Flor is an F10 progeny from a cross of UC Cariblanco N and CIAT accession G25165 . This accession was included in the diversity panel yield comparison under insecticide treated and untreated conditions described above. It has an average sprayed yield of 1650.8 g per plot and an average unsprayed yield of 1106.4 g per plot. This puts its average yield with insecticide protection above average but its estimated tolerance of L. hesperus below average. UC Haskell is an F10 progeny from a cross of UC Cariblanco N and an accession introduced from CIAT accession “P&T 4255” . This accession has been less well studied at it is unknown what degree or mechanism of resilience it may offer to insect herbivory. These results open several avenues of research to pursue in the future. From further analysis of these and other data it may be possible to distinguish resistance and tolerance traits.

Another area of future exploration should be the differences between the defense phenotypes of Lima beans in the Andean and Mesoamerican gene pools. Unlike the Mesoamerican gene pool, the range of Andean Lima beans falls completely outside the native range of the genus L. hesperus and therefore may have fewer defense adaptations to it and closely related species . In figure 12, the only Andean variety, UC 92, supports very low L. hesperus densities. However, as shown in Figure 11, it is very sensitive to L. hesperus with nearly a 50% drop in yield in the unsprayed plots compared to the sprayed plots. This combination of data may indicate that the variety has some resistance traits but very little tolerance to L. hesperus herbivory. Future studies should include additional Andean lines for a more robust comparison of defense phenotypes between the two gene pools.There was a very significant effect of variety on levels of cyanide in flowers, but the effect was not significant in pods . Large-seeded Lima beans, UC 92 and UC Lee had significantly less cyanide in their flowers than UC Haskell. Henderson Bush and UC Beija Flor had intermediate levels of cyanide in flowers that were not significantly different from UC Haskell, UC 92, or UC Lee . Cyanide was not detectable in the mature seeds. The presence of L. hesperus had no effect on the level of cyanide in flowers but it did affect the level of cyanide in pods . Pods collected from plants with L. hesperus had higher levels of cyanide . This may be due to increased enzymatic activity increasing the cyanogenic capacity rather than potential as this was found in Lima bean leaves . There was no significant difference in cyanide levels when considering the interactions between treatment and variety. The duration of the L. hesperus presence, as measured by “Time” in the study, did not affect the level of cyanide inflowers or pods. Nor was there a significant difference in cyanide levels when considering the interaction of treatment and time, square plastic plant pots or time and variety.From the greenhouse experiment of cyanogenic capacity with and without the presence of L. hesperus, insects from each of the cages were counted four weeks after flowering and three weeks after the L. hesperus had been introduced. For thorough counts, all plants were completely deconstructed with each leaf carefully checked for nymphs and the surface of the potting soil carefully searched. It is interesting to note that the two varieties known to be resistant or tolerant of L. hesperus, UC Haskell and UC Beija Flor, have very different numbers of surviving adults and new nymphs. It is possible that cyanogenesis may be responsible for this difference. In quantified studies of greenhouse grown samples without L. hesperus, UC Beija Flor was estimated to have a cyanogenic capacity of 153.4 nM/30 minutes of volatile HCN released from floral bud samples and 224.1 nM/30 minutes of volatile HCN released from immature pod samples .

In the same study, UC Haskell only had a cyanogenic capacity of 126.86 nM/30 minutes of volatile HCN released from floral bud samples and 63.4 nM nM/30 minutes of volatile HCN released from immature pod samples. Given that survival of adult and nymph L. hesperus is more strongly correlated with cyanide in immature pods than in flowers this difference in cyanogenic capacity could be contributing to the higher number of insects surviving on UC Haskell as compared to UC Beija Flor. It should be noted however that this was a no-choice experiment and in the field collected vacuum samples, in which there were many varieties to choose from, more L. hesperus were captured on UC Beija Flor than UC Haskell. This could be due to sampling bias due to the denser canopy of the indeterminate UC Haskell or it could indicate that a mechanism other than cyanogenesis is contributing to L. hesperus preference and survival.This study has several limitations. One is that no wild or weedy types were included for comparison with the domesticated forms. While growing these accessions can be challenging in confined greenhouse space, including them in a similar analysis could shed light on how the trait of cyanogenesis has been affected by domestication. Another limitation of the study was the use of discrete semiquantitative measurements using Fiegl Anger paper. While this method of measuring cyanogenic capacity is well established, safe, and high throughput, it does not provide a true quantitative measurement, or the continuous sampling provided by the enclosed detection systems used for studies of volatilized cyanide .This study found that the survival and reproduction of L. hesperus was negatively correlated with the cyanogenic capacity of their host plant. This indicates that selecting plants with higher cyanogenic capacity in their flowers and young pods may be an effective way to control L. hesperus. The evidence does not support the hypothesis that cyanogenic capacity was induced by the presence of L. hesperus. For consumer safety, future research should determine if there is a relationship between the cyanogenic glucoside content of flowers, immature pods, and the mature seeds they grow into as a correlation has been found between the cyanogenic content of mature seeds and the cotyledons from the seedlings those seeds grow into .Farmland covers more than 35% of Earth’s ice-free terrestrial area, and agriculture is expanding and intensifying in many regions to meet the growing demands of human populations . This trend threatens biodiversity and the ecosystem services on which agriculture depends, including crop pollination . Indeed, recent reviews have highlighted how multiple anthropogenic pressures lead to a decline in wild pollinators such as bees, flies, beetles, and butterflies . However, practices to enhance wild pollinators in agroecosystems are still in development , and considerable uncertainty remains regarding their effects on crop yield and farmers’ profits. Here we review recent research on the topic, including the impacts of certain practices on wild pollinators, crop pollination, yield, and profits . We focus on practices that enhance the carrying capacity of habitats for wild insect assemblages that may then provide crop pollination services; practices to conserve or manage a particular pollinator species are outside our scope although they have received attention elsewhere . We offer general science-based advice to land managers and policy makers and highlight knowledge gaps. Throughout, we emphasize the need to consider population-level processes, rather than just short-term behavioral responses of pollinators to floral resources.Plant–pollinator interactions are typically very general, with many pollinators being rewarded with pollen, nectar, or other resources from several plant species , and with most angiosperms being pollinated by multiple insect species . Humans benefit from this generalized nature of pollination systems, as exotic crops brought far from their ancestral ranges can find effective pollinators within native insect assemblages .

The biased expression of the maternal nuclear copy would resolve any potential conflicts

Our results were also consistent with those previously reported in goldfish, which suggests the most common mechanism for duplicate gene retention in these allopolyploid cyprinines since their 4R event is due to dosage constraints. However, mechanisms for duplicate gene retention are not strictly inferable because a multilevel set of phenomena that range across WGD.Allopolyploids face the unique challenge of integrating two subgenomes, which evolved independently in the diploid progenitors since their most recent common ancestor, that now reside in a single nucleus. One way to resolve potential genetic or epigenetic conflicts is “subgenome dominance”, which results in one subgenome being dominant over the ‘submissive’ subgenome. The dominant subgenome not only has higher gene expression but also retains a greater number of ohnologs compared to the submissive subgenome. To better understand the dynamics among subgenomes of our three sequenced allotetraploids, we compared their gene loss , gene expression level, the density of TEs near ohnologs, constraint on conserved noncoding sequences , DNA methylation patterns and 3D genome structure. To examine gene fractionation differences among subgenomes, blueberry grow bag size gene retention patterns between the two subgenomes of the three allopolyploid species were examined relative to the diploid references from zebrafish, O. macrolepis and Sc. acanthopterus.

These results revealed that in all cases, the maternal subM showed slightly higher gene retention rates relative to the paternal subP . Compared to the reference zebrafish, subM showed 2.815% higher gene retention in L. capito, 0.427% higher gene retention in P. rabaudi, and 0.819% higher gene retention in S. sinensis relative to subP. However, these patterns are not supportive of strong subgenome dominance pattterns as has been reported in some plant allopolyploids . Ohnolog retention bias of certain sets of genes, including BUSCO genes towards one subgenome has been recently reported for the Prussian carp, goldfish and common carp. Similarly, we found that the number of BUSCO singleton genes in maternal subM was significantly higher than those in subP for all three allotetraploids . For example, subM of S. sinensis has retained 609 complete and single copy BUSCO genes, compared to only 448 in subP . Next, we performed GO analysis of the genes that returned to single copy in subP and subM. Functional enrichment analysis revealed that similar GO term classes were identified for all species, including mitochondrial related processes, nc/rRNA processes and DNA repair . These GO terms were also identified as returning to singleton state post-WGD from a previous analysis of the Prussian carp, goldfish, and common carp genomes.A previous study investigating subgenome dominance in octoploid strawberry revealed that the dominant subgenome retained a significantly greater number of tandem duplicated genes. Here, we uncovered a similar pattern for retained tandem gene duplications being biased towards the maternal subM in all three allotetraploid cyprinids .

First, significantly more tandem duplicates are encoded on the maternal subM compared to the paternal subP . Second, a greater number of tandem gene arrays were observed in the maternal subM compared to paternal subP . Lastly, the maternal subM genomes contained a greater number of larger tandem arrays than the paternal subP . An analysis of protein family domains revealed an enrichment of functions associated with the immune system for retained tandem duplicates in these subgenomes .We also tested the hypothesis that DNA methylation patterns in genes and TEs in the extant relatives of diploid progenitor species,and thus subgenomes within allotetraploids, may explain observed subgenome expression bias patterns. Whole-genome bisulfite sequencing of the muscle tissue from two diploid ancestors and three allotetraploids was performed . Levels of CH methylation were very low in genes of all five species , which was also observed in the common carp genome and is typical of somatic tissues in humans. Therefore, we focused on CG methylation for all subsequent analyses. A similar pattern of CG methylation was observed within the gene body and 2 kb flanking regions in case of the diploid Sc. acanthopterus and the three allotetraploid species . However, for O. macrolepis, there was much lower CG methylation levels ~1 kb upstream up to the transcriptional start site and higher levels throughout the gene body and ~1 kb downstream. No difference in CG methylation was observed among the diploid and subgenomes of tetraploid species . However, this analysis of the entire set of ohnologs may obscure more subtle differences.

To examine this, we next analyzed CG methylation for genes with biased expression in muscle tissue towards either the paternal subP or maternal subM. Interestingly, CG methylation levels of expression biased genes towards the subgenome A were lower from ~1.5 kb upstream to TSS compared with subM levels . Similarly, the upstream region of subM bias genes for all species showed lower CG methylation levels in this same region than those of the corresponding regions of duplicated genes in subP . This suggests that CG methylation levels in upstream regions of genes may have a role in observed expression bias towards a particular subgenome. Further, to determine if there are any significant differences in TE CG methylation between subP and subM of tetraploid species, we investigated CG methylation of TEs that are in 1 kb vicinity of 7040 positionally conserved syntenic ohnologs and at the whole genome level. We found some degree of variation in mCG levels between subgenome TEs that were found in 1 kb vicinity of 7040 duplicate orthologs . However, elevated levels of TE methylation in subP were observed in L. capito which was opposite to what was observed in S. sinensis and P. rabaudi where subM showed higher methylation levels. This phenomenon was also observed for TE methylation at the whole genome level . This opposite trend of TE methylation in L. capito in comparison to S. sinensis and P. rabaudi can be attributed to the difference in the TE density of respective genomes .Cypriniformes represent the largest clade of freshwater fish with ~600 described species in the family Cyprinidae, which has experienced multiple rounds of independent WGD. The phylogenetic relationships, evolutionary history, and the genetic basis of previously reported subgenome dominance of these polyploids has remained poorly understood. In this study, high-quality genomes of twenty-one cyprinid fishes, including subgenome-resolved allotetraploid genomes from three tribes, were de novo assembled and analyzed to investigate subgenome evolution at the genetic and epigenetic levels. Our results are supportive of previous reports for subgenome dominance at both the gene retention and transcriptome level. In addition, we observed that the dominant subgenome retained a greater number of tandem duplicates with a functional bias towards immune related processes. Our phylogenetic analyses revealed that S. sinensis, L. capito, and P. rabaudi are allopolyploids and that observed dominance is consistently towards the subgenome contributed by the maternal parent. Also, the most recent polyploid event in P. rabaudi is likely shared with common carp and goldfish. Functional enrichment analyses revealed similar GO term classes, blueberry box including mitochondrial related processes, for the genes that returned to single copy in all examined allopolyploids. The observed consistent bias towards the maternal subgenome donor, alongside the bias towards mitochondrial functions, suggests that observed subgenome dominance patterns in these allopolyploid fish may be due to maternal dominance. The maternal contributed nuclearencoded genes that interact with mitochondrial encoded genes may be favored to maintain proper cytonuclear interactions. The mitochondrial proteome contains products from over a thousand genes, while the mitochondrial genome encodes approximately only 13 proteins.

The vast majority of genes are now nuclear genome encoded following the horizontal gene transfer from the organellar genome to the nuclear genome over the past hundred million years. However, these nuclear genes might encode dosage-sensitive proteins that function in either organellar signaling networks or macromolecular complexes that must maintain proper stoichiometric balance with interacting partner that are encoded in the organellar genome. Furthermore, the sequence of the proteins encoded by both organellar and nuclear-encoded mitochondrial genes may have diverged among the diploid progenitors. Thus, there’s a possibility for incompatibilities to arise from “mismatches” between the genes contributed by the paternal subgenome and the organellar genomes contributed by the maternal parents in allopolyploids. The model that we are proposing here is that observed dominance patterns in these allopolyploids is to preserve proper cytonuclear interactions, and ultimately, core cellular functions. Nonetheless, we cannot exclude the possibility that some ofthe observed subgenome expression differences, particularly at the individual gene level, is due to differences in DNA methylation and transposable element density differences as hypothesized in previous studies. We observed that methylation levels at CG sites in upstream regions of genes, ~1.5 kb upstream to the transcriptional start site, may have a role in observed expression bias towards the maternal subM genome. Epigenetic factors, including changes in methylation at certain CG sites, have been previously shown to alter gene expression and involved in maternal imprinting including of nuclear encoded mitochondrial and DNA repair genes. We also observed that the dominantly expressed ohnolog, from either subgenome, in some cases, had significantly lower TE densities. This suggests that both maternal dominance and TE differences are likely contributing to observed independently repeated subgenome dominance patterns in allopolyploid cyprinid fishes. To the best of our knowledge, this is the first study to show the potential role of maternal dominance in contributing to subgenome dominance in any allopolyploid animal. Future studies of other allopolyploids are needed to determine if these observed patterns are shared by other polyploid animals or are potentially unique to cyprinids. Furthermore, our multi-species comparisons suggest that genetic divergence of the diploid progenitors, for the allopolyploids and divergence times examined in this study, did not contribute to subgenome expression dominance. However, it is important to note the possibility that the divergence of the diploid species in each allopolypoid wasn’t sufficiently different to observe additive subgenome expression dominance effects. Lastly, we also examined genome organization using Hi-C data and selective constraints on noncoding regulatory sequences, which revealed no significant differences among subgenomes. These new reference genomes and various datasets should serve as a powerful platform for the community to further investigate genome evolution of cyprinids, and as a valuable resource for a wide range of studies including modeling human disease.Genomic DNA degradation and contamination was monitored on agarose gels. DNA purity was checked using the NanoPhotometer spectrophotometer . DNA concentration was measured using Qubit DNA Assay Kit in Qubit 2.0 Flurometer . Microgram genomic DNA spiked with lambda DNA were fragmented by sonication to 200-300 bp with Covaris S220, followed by end repair and adenylation. Cytosine-methylated barcodes were ligated to sonicated DNA as per manufacturer’s instructions. Then these DNA fragments were treated twice with bisulfite using EZ DNA Methylation-GoldTM Kit , before the resulting single-strand DNA fragments were PCR amplificated using KAPA HiFi HotStart Uracil + ReadyMix . Library concentration was quantified by Qubit 2.0 Flurometer and sequenced by Novaseq platform w producing 24.39-55.95 Gb raw bases with a bisulfite conversion rate of 99.57–99.75%. MethylC-seq data for each sample were aligned to their respective genomes and methylation called using the methylpy pipeline v.1.4.6. This pipeline uses Cutadapt v.4.1 for adapter trimming, Bowtie 2 v.2.4.4 for alignment, and Picard v.2.26.10 to mark duplicate reads. Spiked-in unmethylated lambda phage DNA was used as a control to calculate non-conversion rates from bisulfite treatment . Gene and TE metaplots were made as previously done using custom scripts and pybedtools v.0.9.0. Gene/TE bodies were divided into 20 bins, and the weighted methylation level calculated across all genes/TEs. For gene bodies, only exonic cytosines were included. This process was repeated for both 2 kb upstream and downstream regions, and the data plotted in R with ggplot2. To examine the effects of neighboring TEs on genicmethylation, we used bedtools v.2.30.0 to identifying genes with an intersecting TE within 1 kb.Domestication has altered the interactions among crop plants, herbivorous insects, and higher trophic levels of agroecosystems . While selection has resulted in enhanced agronomic traits like yield, defensive traits such as toxic compounds have been reduced or removed . A consistent pattern of reduced defense has not been established across species, but crop plants tend to be more vulnerable to herbivory than their wild relatives . Domesticated beans may be more attractive to insect herbivores than their wild relatives, but they are also more effective at recruiting parasitoids .

Dry samples were then weighed and used to calculate an average dry matter percentage

Genomic selection , first introduced by Meuwissen et al. , has been adopted by alfalfa breeding programs in recent years for improvement of a variety of complex traits . By enabling breeders to perform marker-only selection, a cycle of selection could be completed in less than a year, addressing one of the biggest impediments to faster genetic gain – the need for multi-year evaluations . The majority ofGS applications in alfalfa to date are based on phenotypic and genotypic data at the individual plant level . However, alfalfa is marketed as synthetic cultivars with breeding based on the evaluation of families, not individual plants. This is particularly relevant when attempting to evaluate yield. An alternative approach to GS is pooling individuals and genotyping family bulks rather than individuals . Phenotypic data can then be collected at the family level from densely planted plots, a better representation of commercial yield, and the relationship matrix is built based on allele frequency marker data, rather than individual genotyping calls. High quality phenotypic data is essential for maximizing the predictive ability of a GS model . Field trials are inherently variable, particularly in the case of perennial crops with multiple harvests where spatial and temporal differences can have a large impact the quality of phenotypic data. Randomization, blocking, pe grow bag and complex trial designs are commonly used to account for field variation; however, they are often inadequate to account for all the spatial trends in large experiments .

Residual maximum likelihood has become commonplace in plant breeding to estimate variance components and calculate genetic parameters using linear mixed models . Best linear unbiased predictions of random effects have proven effective in predicting breeding values to guide breeding decisions . Mixed-model analysis has the advantage of handling unbalanced data and can incorporate additional information using covariance matrices, providing more accurate predictions . In field trials, experimental units in close proximity to one-another are expected to be more highly correlated than those far apart, and successive harvests likely will be more highly correlated than harvests separated by a greater time interval. Therefore, a better way of managing spatial and temporal variation is to include covariance matrices in the analysis to account for spatial and temporal trends in the data .The training populations used in this study were part of the UC Davis elite non-dormant alfalfa breeding program. Plants were selected from an evaluation trial in Davis, CA and crossed in the greenhouse in 2018 resulting in two populations – UCAL1960 and UCAL1970. In 2019, seed from both populations was germinated in the greenhouse and transplanted to the field at the UC Davis Plant Sciences Farm. Approximately 200 plants from each population were grown in adjacent areas that were each placed under mesh cages to enable pollination with leaf cutter bees . The seed from each plant was harvested individually and the highest seed yielding families of each population were included in the trial.

The entries consisted of 105 half-sib families from UCAL1960, 88 half-sib families from UCAL1970, balanced bulks from each population, and three cultivar checks: Highline, UC Impalo and CUF. A field trial was established at two locations, one in Yolo County and the other in Solano County , on the UC Davis Plant Sciences Farm near Davis, CA in May 2020. Although the two sites are close together, the Solano location is in a floodplain of the nearby Putah Creek. The locations are in a Mediterranean environment characterized by hot, dry summers, cool winters, and moderate annual rainfall, which falls predominantly in the cooler months from November to March . The soil type is a Yolo silt clay loam . Soil tests were conducted prior to planting to adjust P, K, and pH according to soil test recommendations. Seeds were sown in trays in the greenhouse in March 2020 and seedlings transplanted in early May. The trial at each location was a randomized complete block design consisting of two replicates, each with 7 rows and 29 ranges. Plots consisted of 24 plants laid out in a regular 4 × 6 grid with 20 cm spaces between plants. There was a 30 cm space between rows and a 110 cm space between ranges to allow room for mechanical harvesting. This trial was managed as a high-yielding alfalfa stand; soil tests were conducted each year, with amendments made accordingly. The trial was initially irrigated using sprinklers, but after plants were well established, flood irrigation was used to match crop evapotranspiration. Weeds were managed by a combination of manual removal and herbicides, and insect pests were monitored and controlled with insecticide as needed.Dry matter yield was measured by mechanical harvest when the plants had reached 10% bloom.

A self-propelled forage plot harvester with a flail chopper and weigh bin was used to cut plots uniformly at 7.5 cm and measure the fresh weight. There was a total of 12 harvests, three in the establishment year, seven in the first full year of production , and two in the second production year . Each harvest, 20 hand grab samples were collected throughout the day of harvest and weighed before being dried at 60°C in a forced-air drying oven for 7 days. Plot yield was recorded on a dry matter basis after adjusting the fresh weights with the average dry matter percentage.The first fully expanded trifoliolate leaf from the growing tip of the plant was collected from every plant within each family from a single replication of the trial and pooled for DNA extraction. Subsequently, tissue samples were lyophilized and ground. DNA was extracted from the ground tissue using DNeasy 96 plant kits quantified using a BioTek Synergy H1 Microplate Reader with Take3 Micro-volume plate. Genotyping-by-sequencing libraries were constructed using the protocol of Elshire et al. modified as described previously in Li et al. . Briefly, 100 ng of DNA from each sample was digested with PstI enzyme. Two different types of adapters, a barcode adapter and common adapter, were ligated with the DNA fragment using T4 DNA ligase. Equal volumes of the ligation product from each sample were pooled together and purified using AMPure XP beads . The ligation product was then amplified with the Phusion High-Fidelity PCR Master mix , using 50 ng of DNA as the template and 25 nmole of each PCR primer in a 50 ul reaction system. PCR was performed under the following thermal cycling conditions: 72°C for 5 min, denaturation at 98°C for 30 sec, followed by 12 cycles with 10 sec of denaturation at 98°C, 30 sec of primer annealing at 65°C and 30 sec extension at 72°C. The resulting library was again purified using the AMPure XP beads . The DNA fragment sizes were checked using a Bioanalyzer then all families were pooled and sequenced with Illumina Nextseq. The raw reads were compressed using clumpify software from the BBTools package , then trimmed using Cutadapt . The Tassel 5.0 GBS v2 pipeline was used for SNP discovery . Default parameters were used for all steps unless otherwise stated. First, tags and taxa are identified from the raw FastQ files and stored in a local database using GBSSeqToTagDBPlugin. Next, distinct tags from the database were retrieved and reformatted to be aligned with the reference genome using the TagExportToFastqPlugin. Bowtie2 v2.5.1was used to align the filtered reads with the Medicago sativa refence genome developed by Shen et al. using the –very-sensitive-local input. The resulting SAM file created from this step was then used to update the local database with the SAMToGBSdbPlugin. DiscoverySNPCallerPluginV2 then identified SNPs from the aligned tags and updated the database. Finally, ProductionSNPCallerPluginV2 was used to generate a VCF formatted genotype file. SNPs were further filtered with VCFTools using the following criteria: biallelic variants only, growing bags minor allele frequency > 0.05, minimum sequencing depth for each data point > 4, minimum mean read depth over all samples > 30, and maximum missing data of 10%. After filtering, 39,229 SNPs remained for downstream analysis and genomic prediction.

Alfalfa has suffered from a lack of yield improvement over the last 30 years of breeding and research. Traditional recurrent phenotypic selection, although useful for the improvement of simply inherited traits, has been unsuccessful in translating yield gains to commercial producers. In this study, genomic selection was incorporated into the UC Davis alfalfa breeding program across two elite non-dormant populations in an effort to select families for improved yield potential. Whereas previous predictive models in alfalfa have been developed based on genotyping single plants and phenotypic information from space plants or short family rows, this study used bulk genotyping of families, as demonstrated by Andrade et al. , and uses phenotypic data on densely transplanted mini-sward plots. The goal was to better capture yield that could be expected in a commercial planting, while reducing genotyping costs, and aligning the genotyping with current breeding methodology in alfalfa, which is also performed on a family level. To further improve the quality of phenotypic data, spatial and temporal variation were modeled using variance/covariance structures. The predictive ability of GS models is affected by a range of factors, incluing the number, density and distribution of genetic markers, population size and structure, the degree of LD, and the quality of phenotypic data. The LD decay in this study was moderate , but our marker density was sufficient to provide coverage of the genome. In a study optimizing whole-genome sequencing in autotetraploid blueberry , the predictive ability of genomic selection increased rapidly with the inclusion of more markers; however, a plateau was reached with 10,000 markers or more. GBS therefore should provide sufficient markers for accurate GS in alfalfa, as we noted previously . Population structure can induce bias in genomic predictions if not incorporated into the prediction model. Organized breeding programs can develop significant population structure, as observed in this study, where the two populations included in the training population were related, yet clearly distinct. Here we accounted for population structure by including the first two principal components from a PCA in the model. Population size is perhaps the most important factor influencing the predictive ability of GS. A significant increase in predictive ability for fruit yield in blueberry occurred by increasing population size from 120 to ~1500, with an ideal size of 1000 or more . We evaluated relatively few families in this study which may have impacted the predictive ability of our model, even though it represented a maximum number of families we could reasonably harvest and manage. High throughput phenotyping could be a solution to greatly increase the number of families evaluated for DMY in alfalfa without a significant increase in cost to the breeder . Remotely assessing biomass from drone-based cameras and harvesting a subset of plots would allow a breeder to evaluate significantly more families with a similar cost to current phenotyping methods. High quality phenotypic data is essential for producing accurate genomic predictions. The correlation between harvests was significant and positive in this study, and spatial variation is inherent in field trials. Accounting for spatial variation and temporal correlation between repeated measures improves the quality of phenotypic data, thus allowing for more accurate predictions which may translate to genetic gain. Moving forward, we intend to generate a range of populations based on both phenotypic and genomic selection from this yield evaluation trial. These will then be evaluated alongside other elite material for DMY to assess if there has been any improvement. Due to the lengthy breeding process of perennial forages, it will take several years to determine whether these methods have been successful. We will continue to develop the predictive model as more material is evaluated and adjust accordingly. To improve the predictive ability of the model moving forward a combination of evaluating a greater number of families and improving the quality of phenotypic data through better modeling will be imposed. Increasing the size of the training population could be facilitated without a significant increase in costs by using modern high-throughput phenotyping tools, such as drone based remote sensing. With decreasing costs of genotyping, improved computational software and the availability of genomic resources , genomic selection is becoming increasingly available to more resource limited breeding programs like alfalfa. There is still much research required to assess whether actual yield gain can be achieved; however, these studies provide a baseline for future studies to investigate potential yield improvement.

This is an important factor for how much area could even be potentially converted to pollinator habitat

The main correction for this type of mapping project would showing suburban and residential areas do have some ecological habitat value for bees. However, showing how much habitat value a residential parcel has is not simple, but subsequent student designers tried with other techniques shown in section 5.9, below.The goal for the studio students was to create an extensive poster emphasizing designing for bees in a city of their choice. The final project had them concentrate at the scale of up to a few kilometers or miles maximum. One of the best projects include Figure 5, which shows a regional bee habitat design for the SLO area. This student was careful to use both GIS land cover data for more precise habitat mapping as well as mapping information . In their map titled, “The Network Map”, they were able to precisely prescribe where to place new habitat for “pollinator usability” which in ecological terms relates to ensuring foraging habitats are close enough together for continuity between pollinator patches. This student also does a wonderful job demonstrating three simple design techniques in their category “Pollinator Pockets” which could be implemented at various landscape types. So, for example, with major streets, slim narrow planters could be installed; with residential areas, seed packet programs could be implemented; and finally with parks, plastic plant pot suitable foraging flowers could be planted in groups. All three of these design techniques add together for the desired result of increased pollinator habitat circuity and connectivity.

This design solution is akin to adding sunflowers throughout a town to bolster bee network connectivity, as seen in The Great Sunflower Project in Sonoma, California, which was very well received. Overall, Roa’s project is an outstanding example of assessing a city’s pollinator habitat and offering design solutions. This is the sort of spatial resolution landscape designers should be striving for to create bee habitat network plans. This project is a shining example of how bee or pollinator habitat analysis and design could be accomplished. Another final project worth discussing is a pollinator habitat project for densely urban, Glendora, California. This student opted to trace every portion of vegetation from an aerial imagery map using Illustrator . This process, at this scale was painstaking and took several days, but the results are fascinating. They were able to capture the distribution of vegetation through this highly urbanized area. We realized that even the proportion of green spaces were visible between neighborhoods at this scale . Moreover, this project helps to show the proportion of vegetated space versus hardscape in very urbanized areas. This project could have been stronger by identifying how far their focal bee could travel and/or possibly occupy habitat patches within this map. Additionally, this project lacks a vision that the previous one, Figure 5, showed for implementation in typical design situations. However, the attention to detail in mapping green space was outstanding and should not be overshadowed by these constructive criticisms.

Finally, this sort of mapping technique, if used in GIS could be very helpful for improving the urban pollinator network and quantifying existing conditions as well as potential habitat modifications. By applying landscape categorizations to real landscapes, it was possible to graphically demonstrate how bees perceive the greater landscape. Only then, once habitat deficiencies are identified can designers go about fixing habitat short-falls. How different bee genera experience the same landscape will also be important to understand and see spatially, as dynamic landscapes create geographic isolation or continuity depending on bee genera. Future studies should work to untangle how urban habitat performs; provides habitat for various focal bees.This research focuses on developing the field of landscape architecture to better aid in bee conservation design at multiple scales. My goal is to build resilient pollination landscapes for the future. As a profession we must invest in helping stabilize landscapes and their functionality or the effects climate change could be much worse. Climate change is causing ecological stress for bees, both in terms of environmental, including phenological, mismatches , but also in terms of increased physiological issues, including bee body temperature. Bees are keystone species and conserving them, saves so much more than just bees. Pollination functionality underlines the very basis for nearly all terrestrial ecosystems on Earth. Every day, landscape architects and designers create plans to change Earth’s lands and with every plant they choose, they either provide foraging habitat for a bee… or not. Planting designs need to serve functional pollination purpose, to feed the local fauna. Landscape designers, should put an emphasis on bee conservation now to hopefully preserve and prioritize pollination ecosystem services.

Strategic solutions to habitat deficiencies, education and also raising public awareness are all important for resilient landscapes. Now is the time to act. Its time to make great strides in educating not only the new generations of designers, but also, with the greater public. Through implementing designs that celebrate local pollinator diversity, biology, and ecology, we can help to support nature’s and our own, human, future as well.For the last decade or so, horticultural and landscape researchers have been developing sustainable California plant varieties and palettes. Most research on sustainable themes has been completed on low water use landscape plants, but more recently, wildlife friendly too. Field trials for future suitable plants have been completed to help identify, promote and produce climate suitable plantings . For over a decade the UC Davis Arboretum and Public Garden has placed a strong emphasis on using drought tolerant plants, regardless of their geographic origin, the published list being called “The Arboretum All-Stars.” Notably, many of these plants, of non-California origin exhibit traits which this research project has found are also attractive to and functional for foraging native bees. Importantly, this plant palette has been carefully crafted to first prioritize low water use, which is an essential quality facing California’s Central Valley future facing climate change . Consider that in similar urban butterfly studies, it is estimated that around 40% of native host plants which once existed in the California Central Valley habitat areas no longer exist where they once grew and exotic species now provide essential host plant habitat . Thus, it is inappropriate to attempt to eradicate all non-native plants, as many are currently native bee foraging resources. The non-invasive exotic plants that native bees currently use may provide important habitat resources in the future if native plants cannot survive extreme climate change, while the exotic species may be more resilient.Preserving and enhancing pollination ecosystem services are extremely important in landscapes facing great uncertainties with climate change. Ecosystem services are a tool which can be used to help buffer the negative effects of drought, extreme storms, increased or decreased precipitation. Pollination ecosystem services are infrastructure and essential not only for conservation, nature and biodiversity, but also for much of the food humans depend on . Between one-third and two-thirds of the food we eat is the result of insect pollination. Aside from that, the next top priority in California should be water conservation. While pollinator plant lists exist, it is important to explore their success in a real-world setting and examine which ecological aspects could still be improved upon. Ecologically, landscape design can be used to provide and promote bee habitat and connectivity of habitat.Tree death is a natural part of forest dynamics , nursery pots but increasing rates of mortality can result when climatic conditions exceed a species’ physiological threshold . Although directional climate change has historically resulted in shifts in the distributions of species and ecosystems , comparatively rapid shifts in tree distributions attributed to anthropogenic climate change have been documented on all six plant-covered continents . Recent research has focused predominantly on causal mechanisms of tree death, feedbacks to the climate system, and predictive modeling . Ecologists generally agree that trees and forests in temperate regions will shift to higher latitudes and upward in elevation due to warming trends . However, understanding how forests will behave at the ‘‘trailing ends’’ is limited . Stand development patterns following forest mortality events are of considerable interest because they indicate future structure and composition of affected forests, and the ability of these forests to maintain biodiversity and other ecosystem services .

Although widespread mortality events can have negative impacts to ecosystem services , there may be benefits that are also important for adaptation in the human dimension . A global overview of climate induced forest mortality provides a detailed assessment of events driven by climatic water/heat stress since 1970; few of these documented dieback events provide opportunity to examine vegetation changes that occur over a longer time frame. Yellow-cedar , a species distributed from the northern Klamath Mountains of California to Prince William Sound in Alaska, has been dying in southeast Alaska since the late 1800s with intensifying rates observed in the 1970s and 1980s . Recent research reveals a complex ‘‘tree injury pathway’’ where climate change plays a key role in a web of interactions leading to widespread yellow-cedar mortality, referred to as yellow-cedar decline . Prominent factors in this injury pathway include cold tolerance of roots, timing of dehardening, and regional trends of reduced snowpack at low elevations . Early springtime thaws trigger dehardening and reduce snow cover that insulates soil and shallow fine roots from periodic extreme cold events; this can lead to injury of yellow-cedar roots to initiate tree mortality, which is predominantly limited to lower elevations . Despite the extent of research on the mechanisms of decline, overstory and understory dynamics in declining stands are not well understood . The direct loss of yellow-cedar has important ecological, economic, and cultural implications; however, other changes are also relevant in these forests that emerge in response to decline. Researchers are just beginning to understand the influence of dead cedars on watershed nutrient export . Economically and culturally, yellow-cedar trees are important because they provide valuable products for Alaska Native communities and the forest industry . These coastal forests also provide forage for the Sitka black-tailed deer , an important game animal throughout the region. Since the 1980s, much forest-related research in southeast Alaska has addressed the implications of various active forest management regimes on habitat of this commonly hunted species and biodiversity ; aspects of this research centered on old growth habitat and the effects of land use practices, such as clearcutting or partial cutting on forage . To date, researchers have not addressed the effects of yellow-cedar decline on the availability of key forage species. Death of yellow-cedar and the shifts in plant community dynamics in forests affected by decline can have cascading effects on the human-natural system by affecting the ecosystem services these forests provide . We studied the process of forest development using a chronosequence to compare forests unaffected by widespread mortality with those affected at different time points over approximately one century. Considering size classes from seedlings to large trees across the chronosequence, our analysis of the conifer species populations at various life history stages, including death, documented changes occurring in forests affected by decline, and extended a view of forest composition and structure into the future. We hypothesized that: western hemlock and other conifers increase in importance as the contribution of yellow-cedar to the conifer community structure is reduced over time, seedling and sapling regeneration increases as yellow-cedars die and the canopy opens, community composition of understory plants changes over time such that shrubs increase in abundance, and the volume of key forage species for the Sitka black-tailed deer increases in forests affected by decline. Our study illustrates the long-term consequences for many plant species when a single tree species suffers from climate-induced mortality.Modern climate in the southeast region of Alaska is mild and hypermaritime with year round precipitation, absence of prolonged dry periods, and comprised of comparatively mild season conditions than continental climates at similar latitudes . Mean annual rainfall measured in Sitka and Gustavus, the two closest towns to the remote, outer coast study area, measure 2200 and 1700 mm, respectively. The high rainfall that occurs throughout the Alexander Archipelago, combined with its unique island geography, geologic history, and absence of fires maintain some of the most expansive old-growth forests found in North America.

This study shows WHR models are feasible and can be constructed from existing literature

These results suggest that the existing literature identifies few actual foraging associations for the vast majority of bee genera in California, This dearth of information shows the importance of this study and the need for further research in this area. The model independence tests for each respective bee genus observed in the Arboretum reveals that 15 of the 28 models are highly significant at levels below 0.01 with degrees of freedom ranging from 1-237. The remaining models were found to be lacking independence mostly due to a low frequency of use from our observations for those bee genera.According to these results, it is possible to build foraging models for bees based on floral foraging preferences. However, while the aggregate model evaluation scored quite well overall with Table 1 plants being utilized for forage by bees, analysis of the bee-to-plant associations demonstrated that the individual models were not as predictive at the finer scale of bees to specific plant genera, mostly due to the high omission rates. Thus, bees were found to be using many more plants than the predictive plants alone in Table 1. Future research should focus on better understanding bee-to-plant associations. Furthermore, plastic flower buckets wholesale investigating the spatial consequences of bee-to-plant associations is essential in studying habitat connectivity and fragmentation for bee genera.

Studying bees at the genus level was effective for observation of their foraging preferences and trends, which require specific conservation strategies. The variability of results between bee genera emphasizes that, for more effective conservation, beegenera should be studied individually from each other, using an autecological approach, and not aggregated to achieve more effective conservation. Overall, there was much variability in the degree to which bee genera utilized Table 1 plants or not . In contrast to Frankie’s estimate of 17 common bee genera in California , we found 27 bee genera observed to be foraging in the Arboretum, five of which were not predicted to have been there . The 22 predicted bee genera included: Agapostemon, Andrena, Anthidium, Anthophora, Apis, Ashmeadiella, Bombus, Ceratina, Coelioxys, Diadasia, Eucera, Habropoda, Halictus, Hoplitus, Hylaeus, Lasioglossum, Megachile, Melissodes, Osmia, Peponapis, Svastra, and Xylocopa. At 27 observed bee genera, our findings demonstrated only one more bee genus than Robbin Thorp’s confirmed personal collection of 26 bee genera found in Davis, California. It is a good sign that our findings closely confirm the local expert’s specimen collection. Map accuracy was very helpful in determining the presence or absence of bees. Interestingly, as a whole, bees utilized 84 of the 134 plants they were predicted to use in the Arboretum per Table 1. Since bees utilized 297 forage plants, 213 plants were novel plants not included in Table 1.

Error analysis demonstrated variability in predictive success between bee genera. The genera Anthophora, Diadasia, and Habropoda demonstrated the highest true positive fractions; theywere the most predictable in terms of their feeding preferences, ranging from 0.3-0.5. On average, for all bees, the sensitivity score was low, at 0.14, meaning only 14% of the literature’s plants were predicted ‘correctly’ with the corresponding bee association. Conversely, the false negative fraction of 0.86 means that 86% of bee foraging was not predicted via the existing literature models. Therefore, it is valuable to compare separately multiple bee genera and also their foraging preference associations within the same study. Studies which aggregate all bees risk erroneous conclusions regarding the relationships between bee and plant genera. Similarly, studies of only one bee genera are unlikely to be comprehensive enough to encompass conservation of a diversity of bees. We have demonstrated that different angles of model analysis can yield varied ranges of success. Our findings indicate the existing plant lists currently lack effectiveness to comprehensively predict bee genera foraging in California. These results emphasize how more needs to be learned about bee genera foraging patterns. We suggest more research should be done to collect and publish bee association data. Moreover, these findings also show how strongly plant selection can influence garden habitats for bee genera. Table 3 shows how the Arboretum gardens vary greatly in their ability to attract bees and also sustain them over time.

Gardens which facilitate the most foraging over time represent models for how urban planting schemes could attempt to accommodate and conserve California’s bees. The Mary Wattis Brown California native garden and Ruth Storer Arboretum All-Star garden perform the best at providing bee habitat. Therefore, future garden design to support bees should be based off of the planting characteristics of both garden types,maximizing plants best at supporting bee foraging. Furthermore, almost all bees were found to forage on exotic plant species which were not emphasized on some of the recommended suitable plant lists . The observed opportunistic foraging nature of bees shows that reconciliation ecology seems in part suitable for maximizing bee habitat conservation and design. This study presents compelling evidence that that current suggested plant lists have significant inadequacies as habitat for bees. Furthermore, this research has identified many plants within the Arboretum’s novel plant communities that are highly functional as habitat for California’s native bees . To further improve bee habitat relationship model information, habitat elements required to support bee nesting need to be included. In future research, it would be valuable to add further layers of information, such as nesting resources, to increase the helpfulness of the model toward building or predicting comprehensive suitable habitat, not just foraging. For example, the primary means of reproduction by most native bees is usually one of three strategies : ground nesting types, cavity nesting in trees or plant stems, and plant stems, both woody and herbaceous. There may be also be garden physical form factors regarding pollinator attractiveness and this should be addressed in future research. For example, the age of the plants, degree of under or overstory, sun aspect ratio, and more. Detailed bee landscape design is explored in Chapter 2 in depth.Since the existing bee-plant lists seemed to be under performing, we conducted further investigations into these list inadequacies. The existing literature plant list included 61%native California plants and 38% non-native, while our foraging study showed nearly an inverse in the actual feeding trends, with 43% native California plants and 57% non-native. In addition, among plants which were not listed in Table 1, “unexpected novel” foraging plants, only 36% were California natives, while 64% were non-natives. Scientists may have overvalued the forage attractiveness of native plants for bees and undervalued the positive function of nonnative plants. For example, an urban ecology study in the Central Valley of California by Shapiro found 40% of lepidopteran faunal host plants are non-native plants. While California native plants contribute many evolutionary and ecological aspects to a site, black flower buckets our study showed that bees in horticultural environments exhibit trends of feeding opportunistically and without majority preference for the native status of a plant. As climate change increasingly puts stresses on native plants forcing range contractions and expansions, exotic plant species may ultimately need to be utilized to maintain and conserve native bees.

The UC Davis Arboretum and Public Garden hosts a variety of specialty gardens, boasting plants curated from geographically distant locations. However, the Arboretum plantings have been consciously selected to suit the local Mediterranean climate. The Arboretum tends not to grow plants which require intensive care and high water demands, rather they take low maintenance approach, with policies about not fertilizing and/or spraying for pests, for example, and this is beneficial to bees. Furthermore, the Arboretum has a plant promotion program called the “AllStars” and one hundred Arboretum All-Star plants were selected for climate suitability, some native, some not, with an emphasis on low-water use plants and also attractiveness from the human perspective . Many of these plants were also found to be utilized by bees. This study allowed a unique look at how California native bees perceive unique foraging opportunities, some of which were previously unknown, but seemingly beneficial. Among nonnative plants there are consistent trends of plant origin . Based on observed foraging results, from a bee’s perspective, the Arboretum would improve bee foraging by adding more drought tolerant plants from Africa, Australia, Europe, South Africa, South America, the Canary Islands, and New Zealand.It is essential to better understand bee foraging from a bees’ perspective in order to design suitable habitat. This research has demonstrated the limitations of currently understood bee-toplant associations. The basic connection of a bee to its foraging plants is essential to understand which types of habitats are appropriate or not. WHR accuracy can be improved with extensive empirical testing. This study elucidates the need for further quantitative studies on bee WHR models. Through accommodating native bees in novel ecosystems, such as in gardens and in hedgerows near agricultural crops, valuable pollination services can be conserved for the benefit of people, agriculture, and natural communities.Reconciliation ecology provides one of the best frameworks for resilient bee habitat design. The data from this study have shown there are benefits to using both California native plants as well as non-natives to achieve peak foraging habitat for California native bees, as well as naturalized European honey bees. As seen in the Arboretum’s themed gardens, bees are utilizing the novel plantings which also were selected for drought tolerance and aesthetic beauty. By combining the best foraging plants, landscape designers can significantly improve foraging habitats for native bees. Best management practices and design responses for improving bee habitat are likely to be very specialized to increase foraging optimization at designed garden sites. For example, in California’s Central Valley drought tolerance is a very important plant attribute since water will likely become less available in the future due to climactic water deficit projections . Therefore, a design framework should be responsive to a site’s prioritized ecological needs.Although the data collected for this study is extensive in time and space, it represents just one year of bee-plant visitation in the Arboretum. Thus, interannual variability of weather and climate could produce other visitation patterns. However, we feel we captured much intraannual variability with our sampling approach. This study was conducted in the Central Valley ecoregion of California and this could produce biased results towards the bees associated with this ecoregion. As noted in the Introduction WHR models tend to be regional in nature and therefore results may not transfer to all other areas of California. All methods for bee fieldwork have gone under scrutiny in recent years, which are not unique to this study. It also should be pointed out that this study did not explicitly evaluate nesting habitat at the study site.Development of WHR models for native bees is an important step in conservation planning for these species. They can be significantly improved with additional field studies regarding bee foraging plants. In particular, it is important to study bee-to-plant associations, which are currently not well understood or, in the case of novel associations, underrepresented. More quantitative foraging studies should be undertaken to make accurate foraging models. Results from this study show that more optimal pollinator plantings could be designed and constructed in a more strategic and scientific manner. Using data from this study, plant palettes for ideal garden design for native bee conservation could be created from a bee’s point of view. Many bees in this research project were found to be opportunistic foragers—neither exclusively utilizing native or non-native plants. Thus, bees may have more opportunities to thrive if they are given a broader spectrum of plants, which can be done through strategic conservation actions in both space and time. This study elucidates the feeding preferences of native bees, which can be used to better manage and conserve them in the California landscape. How many other foraging plants would native bees use or prefer? This can only be answered with further empirical studies and assessment of possible plant candidates for native bees.Bees are a diverse suite of insects which provide the greatest percentage of plant pollination in the world . The ecosystem services that bees provide are essential to people and the world’s ecosystems. Bees also pollinate the majority of food plants that humans consume . It is estimated that up to two-thirds of food crops for humans require bee pollination of which native bees have been estimated to contribute, around 50 percent of bee foraging activity .

One major finding from floral studies is that microbes are dispersal limited at regional scales

While there are no studies looking at the connection between floral topography and seed microbial transmission, experiments with flowers and leaves have demonstrated that bacterial dispersal is influenced by plant surface topography and surface water distribution. Conducting similar micro-scale inoculation experiments like these in flower to-seed systems will illuminate how microbes actually move. Seeing that microbes can be florally transmitted to seed, we need to consider studies on the dispersal of floral microbes to understand seed microbial communities at the meso-and macroscales. For example, Belisle et al., 2012 found that yeast frequency in nectar communities of Mimulus aurantiacus was correlated with flower proximity, and they inferred that dispersal limitation was controlled by pollinator behavior. In a study on the floral microbiome across wildflower species of California, Vannette et al. observed that fungi were more dispersal limited between individual flowers and plant species than bacteria. Another major finding has been that pollinators can vector microbes between flowers and influence microbial community patterns. For example, Vannette and Fukami explored the variable effects of dispersal limitation on beta diversity in the nectar microbiome. Using a pollinator exclusion experiment, procona valencia buckets they found that increased dispersal by pollinators raised beta diversity and hypothesized that this increase was due to the stochasticity of dispersal timing which strengthens priority effects .

A pollinator exclusion experiment in B. napus demonstrated that pollinators can also vector bacteria to seeds through flowers, impacting the local and regional diversity . These experiments indicate that dispersal may have unique effects on diversity in flower microbiome meta communities via arrival history. However, all of these studies were performed only on the macro-scale, and they did not characterize dispersal traits. Furthermore, the association between dispersal patterns in floral microbial communities and those in seed communities has yet to be studied. Future experiments should explore if dispersal traits and arrival history consistently enhance beta diversity in flower and seed microbial communities among spatial scales .Ecological drift is defined as random fluctuations in species abundances over time, and can be driven by random birth, death, and migration events . Drift is particularly important when local communities are small and filtering is weak . This is key to note for seed microbial communities because they typically have low population sizes and low species richness . Random migration events may be particularly important for seed microbes, such as those vectored by rain or wind . However, while there is a lot of interest in drift and stochasticity in seed microbe research , drift as a process is difficult to study because it is hard to manipulate. Meta community ecologists have also generally found it difficult to get direct evidence of drift, with limited examples from experiments testing the coexistence of ecologically equivalent taxa . One alternative approach to direct observation used in plant microbiota studies is to fit community data to neutral models, where community members are assumed to be ecologically equivalent, and non-significant variation in community composition across samples is explained by neutral processes.

Rezki et al. took this approach when studying the seed microbiota of R. sativus by fitting fungal and bacterial community data to a Sloan neutral model . This model accounts for neutral birth, death, and immigration rates, and estimates immigration rates into communities based on species frequencies across samples . Immigration rates within the confidence interval of the predicted values imply that drift is structuring the community . Based on the model, they found that bacterial community assembly was driven primarily by drift, while fungal communities were driven more by dispersal . This study indicated that drift is important for some seed microbes, and more model-fitting studies or coexistence experiments are needed.At its core, meta community ecology emphasizes not only how the processes described above play out individually, but also how they interact with each other to produce emergent community patterns across scales. In plant microbiome research, the interaction between abiotic and host filters, also known as genotype-by-environment interactions, has been of growing interest because it provides a more holistic explanation for microbiome variation . Such an explanation can be applied to seed microbial communities, which may vary with seed nutrient profiles, osmotic stress, and water availability. However, as previously mentioned, GxE studies on plant microbiota face a scale problem where genotype and environment become synonymous at the micro-scale. Taking a plant trait-based approach to these studies may make the role of these effects more clear, and can connect micro-and macro-scales via host local adaptation. While not emphasized as much as GxE interactions, the interaction between dispersal and filtering is also important during seed microbiome assembly.

At the micro-scale, variation in the plant surface landscape can create differences in dispersal limitation between taxa. Doan et al. demonstrated this interaction on synthetic leaf surfaces, finding that surface water acted as a conduit for bacterial dispersal. This effect may also be present in floral stigmas, which are highly heterogeneous landscapes . Indeed, in their work on transmission of the pathogen A. citrulli from watermelon flowers and fruit to seeds, Dutta et al. found that inoculum from the flower dispersed more frequently and ended up in deeper seed tissues than inoculum from the fruit. While these examples suggest that heterogeneity in the plant landscape impacts dispersal limitation to seeds, more studies are needed. Dispersal also intersects with species interactions, most clearly through historical contingency or priority effects . In this phenomenon, the arrival order of community members dictates assembly outcomes, typically with an advantage to taxa that arrive first . Priority effects can occur either through niche preemption, where the first colonizers fill all available niches, or by niche modification, where the first colonizers alter the environment and its resulting niches . These effects are often cited as important in seed communities because they have few members . However, priority effect experiments in plant microbiota have typically been done in leaf and wood communities . As such, there is a need to understand the role of priority effects in seed communities. An exciting new approach for studying the multiple, interactive processes of dispersal, filtering, drift, and species interaction is with Joint Species Distribution Models , which extend single-species distributions to community-level dynamics . Leibold et al. used these models in tandem with variation partitioning to explain the internal structure of simulated meta communities. They found that this approach was a promising way to connect meta community pattern data to multiple assembly processes . In the seed microbiology literature, Fort et al. used JSDMs to infer how maternal filtering and abiotic filtering contributed to seed mycobiome variation in Q. petraea seeds . They found that fungal guild influenced which taxa varied with abiotic filters, with elevation selecting saprotrophs and seed specialists, and all taxon co-occurrences were positive associations . While JSDMs were not used in a meta community context for this study, and they are limited in their omission of abundance data, these models provide an integrative approach for looking at seed microbiome assembly.Multiple tools exist for exploring and exposing the effects and interactions of filtering, species interactions, dispersal, and drift on microbial community assembly of individual seeds at multiple spatial scales. However, future work can do a better job of integrating and connecting meta community ecology models to traditional seed microbial ecology studies at micro-, procona buckets meso-and macro-scales. One technical challenge of taking this approach pertains to the interrogation of microbial communities in individual seeds. Culture-based studies of individual seeds report low isolation frequencies, with most seeds containing zero or one microbial taxon . Additionally, most sequence-based studies to date pool seeds by fields or other groupings . As exceptions, Bergmann and Busby and Fort et al. sequenced fungi from individual tree seeds and found that sequencing depth was fairly high. However, the tree species in these studies produce large seeds; sequence-based detection of microbiota might be more difficult in small-seeded species .

Additionally, it is often difficult or impossible to treat individual seeds as independent since experimental treatments or predictors are often applied at the fruit or plant level. To resolve these issues, future work could focus on species where seeds are fertilized independently , or one seed per fruit/plant could be sampled for large-seeded species. Alternatively, seeds could be pooled at the fruit or plant levels, since these are the levels where treatments are often applied and they sufficiently capture the variation in seed microbiota while still allowing for a meta community approach at the meso-and macro-scales. The appropriate level of pooling should be selected based on the transmission pathway of interest . Finally, seeds of large-fruited species could be pooled by parts of the fruit/pod for spatially explicit sampling at the meso-and micro-scales. These scale-explicit pooling approaches, along with the use of additional methods at the meso-and micro-scales , will allow for characterization of microbiota at or near the individual seed level while mitigating issues of low DNA amounts and cross-contamination. At the micro-scale, there are many opportunities to take a traits-based approach to host filtering of seed microbiota. Experiments can go beyond studying if plant traits have an effect to testing what these effects are . These experiments could also take a microbial trait-based approach to host filtering and identify the genes required for successful transmission, which are still largely unknown . This could provide valuable insights into the genes required for transmission across the different pathways . Furthermore, meta genomic analyses across plants, populations, species, etc., could determine if these transmission-associated genes are common across meta communities. Such information could show if there is functional conservation across microbial communities, even if they are taxonomically variable. Finally, micro-scale experiments can also test how microbial community assembly is impacted by the interplay between deterministic and stochastic processes. In addition to these tests of microbial and plant trait impacts, experiments testing the role of dispersal in seed microbial community assembly among spatial scales should be conducted. At the micro-scale, experiments using synthetic microbial communities on stigmas with varying chemistry and topography can demonstrate how dispersal and selection occur between flowers and seeds, and what the role is of plant genetics and microbial adaptations. At the macro-scale, pollinator exclusion experiments similar to those in Prado et al. could be conducted across sites in natural landscapes. By using sites at varying distances and connectivity levels from each other, and analyzing both within- and among-site seed microbial community variation, one may obtain new information about how pollinators and patch connectivity impact multi-scale dispersal ability. These proposed studies would elucidate how dispersal contributes to meta community assembly among spatial scales. Along with these single-process studies, we envision studying the interactions between processes through both observational and experimental studies. As JSDMs continue to be refined to model nested and continuous meta communities, they will provide a way to analyze seed microbiome patterns and their associated assembly processes that is more sophisticated than previous modeling approaches. Additionally, priority effect experiments conducted at multiple points in the seed life cycle may reveal how historical contingencies impact seed microbiome assembly throughout the seed life cycle. Such experiments would also test the Primary Symbiont Hypothesis , which argues that seed communities are dominated by a single microbe with significant functional consequences for the plant. Finally, questions will need to be asked about seed microbiome assembly that go beyond just testing for spatial mechanisms. Primary among these questions is: what fitness benefit does transmission into seeds provide to microbes and their host plants? Such a question gets at the eco-evolutionary dynamics in these microbial meta communities, which can have long-term consequences for both microbes and plants. Because microbial communities behave and evolve at shorter time-scales than macro-organisms , it is feasible to design simple experiments testing how microbes evolve in response to plant defenses, nutrient availability, and micromorphology. Such eco-evolutionary studies may have applications in understanding microbial community shifts with crop domestication . Additionally, both microbes and seeds have dormant stages, which can impact meta community dynamics through tradeoffs with dispersal and delayed responses to environmental conditions . The role of dormancy in seed and plant microbial meta community assembly has yet to be explored, so there is much room to study how dormancy impacts these systems over longer temporal scales. Finally, a hot topic in plant microbiome research is how to modify plant microbial communities for climate resilience and other beneficial traits . However, the impacts of climate change-associated disturbances on plant microbiomes have been limited to pattern-based studies in leaves and roots .

Tomato provides a good example of the relationship between pollination and yield

Comparing the pollen-provision microbial communities of closely related specialist and generalist bees will help address how diet breadth affects microbe exposure and establishment. To help understand how diet breadth affects bee nest microbial communities, we characterized the pollen provision microbial communities of a pair of closely related bees collected across four sites: Osmia lignaria, a polylectic bee that prefers to forage on orchard trees in the family Roseaceae; and O. ribifloris, a specialist on blueberries and its relatives with Berberis serving as an alternate host. We address the hypothesis that, due to the importance of floral transmission for microbes found in the pollen provisions of megachilid bees, diet breadth affects pollen provision microbial composition. We predicted that the specialist bee would harbor fewer species of microbes and show less variation across sites compared to the generalist bee.We collected pollen provisions from the nest cells of two bee species native to the area of Logan, Utah, USA: Osmia lignaria propinqua and Osmia ribifloris biedermannii in Spring 2015. We chose four locations near Logan that varied in floral composition and interaction with other bee species: Site 1 was a cultivated floral garden located in a suburban area with a combination of native and introduced floral species, plastic planters wholesale including Rubus idaeus L. and Calluna vulgaris .

Site 2 was a bee-tight screenhouse that contained older, in-ground plantings of R. idaeus and Berberis fremontii , along with potted C. vulgaris. Sites 3 and 4 were natural areas of Logan Canyon that contained limestone cliffs and a mixed forest, and were approximately 1836 m apart, which is likely outside of the flight range of Osmia species. We note that the locations differed in the composition of their pollinator communities: Site 1 had a diverse pollinator community present, including numerous honey bee and bumble bee colonies. Site 2 was a research screenhouse that had no honey bees nor bumble bees present for at least five preceding years. Sites 3 and 4 had diverse pollinator communities without honey bees present, and we did not collect any O. ribifloris from Site 4. See supplemental Table S1 for the number of each bee species collected at each site. There were also differences in the sources of sampled Osmia nests between locations: We sampled natural Osmia populations from Sites 1, 3, and 4, while we sampled nests made by commercially acquired bees for Site 2. In this location, O. lignaria were sourced from Watts Solitary Bees and O. ribifloris were purchased from NativeBees.com . Within each location, we sampled the pollen provisions from cavities within wooden nesting blocks that were lined with paper straws and placed in each location for one week. As both species of bee will nest in these cavities at the same time, we were able to concurrently sample newly made nests. We used 1–6 nesting straws per bee species, then removed the straws from the nest blocks, numbered them sequentially, and X-rayed the straws to confirm that eggs were present and had not yet hatched.

Once the straws were collected, we carefully slit open the straws, removed the pollen provisions from each initially formed cells , then sterilely removed the eggs from nest cells before DNA extraction. We carefully excluded cell wall partitions from the collection of the pollen provision. Each cell was treated as a separate sample.We extracted DNA from the pollen provisions based on a modified protocol from Engel et al. 2013, Rothman et al., and Pennington et al. 2018. We used a Qiagen DNeasy Blood and Tissue kit for the DNA extractions, with slight modifications. We sterilely transferred entire pollen provisions to a 96-well tissue lysis plate , then added 100 µL glass beads, one 3.2 mm steel-chrome bead, and 180 µL buffer ATL to eachwell. We bead-beat the mixture at 30 Hz for three minutes, turned the plates over, then bead-beat for three more minutes. We incubated the mixture at 50 ◦C overnight, then followed the rest of the manufacturer’s protocol to finish the DNA extraction. We used the extracted DNA to prepare 16S rRNA gene libraries as in McFrederick and Rehan 2016, Rothman et al. 2019, and Pennington et al. 2017. Briefly, we used a dual-indexing approach to build an amplicon construct consisting of the universal primers 799F-mod3 and 1115R, a unique 8-mer barcode, and the Illumina adapter sequence as in Hanshew et al. 2013. We built the libraries in two rounds of PCR amplification: First, we used 4 µL of template DNA, 0.5 µL of 10 µM forward and reverse barcoded primers, 10 µL water, and 10 µL of 2× Pfusion DNA polymerase with an annealing temperature of 52 ◦C for 30 cycles. Next, we cleaned the PCR product with a MoBio UltraClean PCR cleanup kit . We then used the cleaned amplicons as template for another PCR reaction, with the following conditions: 1 µL template DNA, 0.5 µL of 10 µM primers PCR2F and PCR2R, 13 µL water, and 10 µL 2× Pfusion DNA polymerase with an annealing temperature of 58 ◦C for 15 cycles.

We cleaned and normalized the libraries with a SequalPrep Normalization kit , pooled 5 µL of each library, then cleaned and concentrated the libraries with a MoBio UltraClean PCR cleanup kit . Lastly, we used an Illumina MiSeq to sequence the libraries at 2 × 300 cycles in the UC Riverside Genomics Core Facility.Other studies have found that bee species identity, foraging patterns, and geography can affect the microbes found in pollen provisions. McFrederick and Rehan showed that the pollen provision microbiome of Ceratina calcarata co-varies with pollen usage across habitats, especially when considering fungi. Voulgari-Kokota et al. showed that bee species and foraging patterns drive the bacterial communities found in pollen provisions. Our findings add to the consensus of these other studies that floral transmission and pollen usage influence the composition of the pollen provision microbiome. Our study also extends these previous studies by adding an understanding of how diet breadth does and does not affect the pollen provision microbial community in these two species of Osmia. Whether the different bacterial variants found in the host species and sites across studies is due solely from different transmission networks, filtering of microbes via differential floral and pollen provision chemistry, or a combination of both requires further study. As in other bee species’ pollen provisions, Osmia spp. pollen provisions contain a wide diversity of microbes. In agreement with these previous studies, the bacteria identified here can also be found in the environment, in both flowers and soil. For example, four of the 18 ASVs that were found at all sites in our study belonged to the genus Acinetobacter. This genus is ubiquitous in the environment and has been found in association with plants and animals, in floral nectar, in sewage, and in water and soil. The small fragment of the 16S rRNA gene that we use for Illumina sequencing rarely allows species-level resolution of these bacteria, so we are unable to determine whether these taxa were sourced from the mud that the bees used to partition their brood cells or from the nectar they mixed into their pollen provisions. Acinetobacter, however, has been reported in association with pollen provisions of several different bee species, and it is therefore not surprising that it occurs at all study sites. One conspicuous group of bacteria that are found in many pollen provision microbiomes but were uncommon here are bacteria from the Apilactobacillus micheneri clade. These bacteria have been found in the pollen provisions of other megachilid species in North America and Europe, but it is becoming clear that they are unevenly distributed across species. For example, in Germany, Apilactobacillus spp. are abundant in the pollen provisions of Megachile spp., plastic plant pot at variable abundances in O. caerulescens provisions, at low abundance in Heriades truncorum provisions, and absent in O. bicornis and O. leaiana provisions. In Texas , we detected these same lactobacilli at high relative abundances in pollen provisions of Osmia chalybea, Osmia subfasciata and Megachile policaris. Here, we report that lactobacilli are present only at very low relative abundances in the pollen provisions of O. lignaria and O. ribifloris, and that the A. micheneri clade lactobacilli are absent. As these bacteria have been isolated from flowers in both Texas and California, we hypothesize that either foraging preferences or pollen provision chemistry drives the presence or absence of these lactobacilli, and understanding the apparently cosmopolitan phenomenon of uneven distribution of lactobacilli across wild and solitary bee species should be a priority for pollen provision microbial community studies.

Many of the bacteria that we identified in O. lignaria pollen provisions have also been found in association with O. lignaria adults. Cohen et al. found that adults that had been foraging in the environment had a different or more variable microbiome compared to bees that emerged under sterile conditions in the lab, again supporting the importance of environmental transmission for the wild and solitary bee microbiome. Many of the bacteria found in adult O. lignaria adults are found in the environment and in the pollen provisions that we studied here. For example, Massilia is a root-colonizing soil bacterium that has been found in O. bicornis nests. This bacterium may be found in adult O. lignaria and their pollen provisions due to the adult’s habit of collecting mud to build partitions between brood cells in their nests. Pantoea is a plant-associated microbe that is abundant in the environment but has also been reported in association with other solitary bees and honey bees. Surprisingly, Cohen et al. found lactobacilli associated with adult O. lignaria, and additionally found that the abundance of flowers at a site positively correlated with the relative abundance of lactobacilli associated with adult O. lignaria. Future studies examining the adult and pollen provision microbiome of O. lignaria will help unravel how the environment shapes the microbiome of these separate but connected niches. The sole bacterial family that we found to be differentially abundant by host species was Micrococcacea. Micrococcacea is a diverse family that occurs in the environment, includes a commensal but opportunistic human pathogen, and has been classified as a pathogen honey bees. As Micrococcacea represented around 4% of the reads in seemingly healthy O. ribofloris nests, it is unlikely that this bacterium is pathogenic to Osmia. The differential abundance between the two Osmia species may mean that it is somehow important for bee health, but may also be due to the differential abundance between sampling site too. Along this thought, we found that many bacterial families were differentially abundant between geographical location. This provides further support to the hypothesis of environmental and floral transmission for solitary bee microbes as has been shown in previous studies.Insects pollinate many plant species, including several major crops. Bees are the single most important insect pollinator group and can be a limiting factor for the success of plant reproduction. Consequently, there is strong inter- and intra-specific competition among plants for the attention of pollinators. With respect to insect-pollinated crops, pollinator visitation is required to obtain maximal seed and fruit production. Consequently, pollination facilitates higher yields even when a crop plant is self-compatible. Bumblebees are important pollinators of tomato and other Solanum species that utilize an unusual pollination system called ‘buzz-pollination’. Buzz-pollinated flowers provide excess pollen as a reward to foraging bumblebees that feed it to their developing larvae. Although domesticated tomato is to a large extent ‘self-fertilizing’, buzz-pollination by bumblebees or by manual application of mechanical vibration ‘wands’ is required for maximal seed production, which in turn promotes increased fruit yield . Cucumber mosaic virus , one of the major viral pathogens of tomato, is a positive sense RNA virus that encodes five proteins including the 2b protein, which is a viral suppressor of RNA silencing. Bees do not transmit CMV but the virus is vectored by several aphid species. Virus infection causes dramatic changes in plant host metabolism . CMV-induced metabolic changes include qualitative and quantitative alterations in the emission of volatile compounds and in certain host species this makes infected hosts more attractive to aphid vectors.It is not known if the virus-induced alterations in host volatile emission that influence aphid behavior can also affect plant-pollinator interactions. Most bee-plant interaction studies have focussed on the effects of visual cues. Therefore, the influences of floral and non-floral volatiles on bee-mediated pollination are less well understood.

Health-protective effects have also been described for other phytochemicals

Oregon weather data were used to illustrate the impact of the different levels of parasitism on D. suzukii populations.In the U.S.A. and Italy, the current suite of parasitoids attack only late larval and early pupal life stages of D. suzukii, and field rates of parasitism are estimated to be around 2 % . In the biological control model runs, we made the simplifying assumptions that parasitism remains constant and that parasitism is attributable only to limited biological control agents based on field observations . For Oregon, we also assumed abundance of alternate hosts in surrounding vegetation from March 1, and that oviposition is possible as soon as the reproductive thresholds are met . Collections from South Korea show parasitism rates as high as 17 % . For these reasons, 2 and 15 % parasitism on the late larvalearly pupal life stages of D. suzukii were incorporated into the model on a season-long basis and compared to populations exhibiting no parasitism. These model runs would illustrate current and potential future population impact because of classical biological control using parasitoids.This research suggests that DD are very useful for estimating physiological time, and their use extends beyond phenology models. In addition, we describe the impacts of IPM strategies on early-season D. suzukii populations presumed to consist mainly of adults. The data presented are an alternative to the model developed by Wiman et al. to estimate populations of D. suzukii within environmentally known thermal thresholds. Previously, 10 plastic plant pots insects were forced to progress through the age-structured matrix with mandatory daily advancement without taking into consideration physiological age.

An advantage of the current model structure is that it allows flies to proceed through the physiological age-structured matrix contingent on accumulation of physiological time as approximated by DD. This refinement to the model allows survival and reproduction of concurrent cohorts comprising a population to be more precisely estimated. The two functions that describe mortality and fecundity within thermal thresholds in physiological time provided good fit using standard population survival and fecundity fitting techniques. We realize that a partial or complete shift of phenotypes inresponse to temperature extremes is not taken into account. However, future modeling efforts will be aimed at capturing phenotypic plasticity and its implications for survival and reproduction of populations experiencing extreme environmental conditions. We provided parameters of D. suzukii survival outside of thermal development thresholds, simulating a sudden transition to extreme cold. The fitted exponential function describes the survival rates of un-acclimated D. suzukii populations under prolonged environmental conditions. When exposed to suboptimal temperatures, un-acclimated D. suzukii populations clearly display high mortality rates. Theoretically, the horizontal asymptote of the fitted exponential function will not result in the extinction of winter or summer populations under these extreme conditions, and in the field, this phenomenon is observed as populations build during periods of favorable temperatures. Investigation of this observation, however, falls outside the scope of the current study where we focused on the growing season.

We therefore realize that the ability to estimate population response to a wide range of changing environmental conditions will further improve population estimation of adaptive invasive species such as D. suzukii. One step was made to more clearly understand how improving environmental conditions impact winter-surviving D. suzukii reproductive potential. Dissections of latedormant female D. suzukii from Oregon and Washington, U.S.A., and Trentino, Italy, displayed increasing levels of reproductive maturity as DD accumulated. Our laboratory data show that egg laying is initiated at 210 DD, whereas field increases in reproductive potential range from 50 to 800 DD. There may be various explanations for these differences, including year, microclimate, genetic variability, and trapping techniques. Although the timing of reproductive maturity with DD in these studies differed, the data nevertheless show a clear relationship between reproductive potential and physiological time. These data illustrate the importance of more suitable environmental conditions as a factor contributing to increased reproductive potential of D. suzukii. We realize that temperature is not the only factor contributing toward such increased reproductive potential, as the roles of humidity, alternate food sources, and host media during the late dormant period of D. suzukii may also need consideration . Our data suggest that in some regions, female sexual maturity may occur very early in the growing season before hosts become widely available. The lack of suitable fruit hosts in certain areas early in the growing season likely makes D. suzukii more dependent on alternate nutrient sources. These sources include pollen and nectar, which may be utilized by D. suzukii to increase reproductive potential and survival levels.

Such resource availability may be a factor resulting in variability of our data. For these reasons, it is important in the future to describe the contribution of these different factors to D. suzukii population dynamics. Clear parameters should be developed to more accurately model early-season D. suzukii population increase. We used a model to illustrate the importance of key periods when pest population structures can be exploited to the advantage of IPM . To the best of our current understanding of the structure of D. suzukii populations, adult females comprise the majority of overwintering individuals . Thus, early spring is a key period when only sexually maturing adults are present. Elimination of these initial adults before the population becomes established and spread among different ages would be ideal. In Parlier, California, U.S.A., a second key period exists during the summer when suboptimal hot temperatures prevail and D. suzukii populations decline. Insecticide A had a longer-lasting residual against adult stages and a shorter residual against immature stages as compared with Insecticide B, which explains the strongest effect on the early-season population. Insecticide B model runs simulated compounds designed to cause mortality at all life stages of D. suzukii . These compounds are of increasing importance as the population age structure becomes more diverse. In the early season, as in our simulation, adults dominate. During latter portions of the season, larger portions of D. suzukii populations are expected to be in immature life stages. During such latter portions of the growing season, there should be increased focus on pesticides that target all life stages of D. suzukii. The use of adulticides to prevent fruit rejections due to infestation before harvest when all ages of flies are present, however, will remain a key component of IPM programs. Additional considerations include the development of D. suzukii insecticide protocols to minimize the development of insecticide resistance. Genetic control using RNAi biopesticide technology with seven-day residual periods against all life stages resulted in minimal population reduction compared with untreated populations. Again, targeting of larvae makes more of a difference on structurally diverse populations. The utility of this technology to target immature life stages of D. suzukii may make it more effective at curatively managing pest populations compared to the adulticides during periods when immature life stages dominate the population structure. An additional advantage of such genetic pest management tools is that they can be designed to be species-specific and target D. suzukii only. Disadvantages currently include regulatory obstacles, wary public perception, and potential incompatibility with organic production practices. Both levels of biological control inputted into model runs resulted in D. suzukii suppression. Higher levels of parasitism will probably result in concomitantly lower levels of pest pressure during the harvest period. Currently, only low levels of biological control are found in most production regions . Model runs indicate that levels of parasitism close to 15 % will result in significant reductions of pest populations during the early portion of the season due to the loss of a portion of the initial crop-infesting population. Biological control is only effective, however, plastic pots large if it can suppress pest populations as the crop ripens. Model outputs estimate lower levels of suppression during the earlier portion of the season at the 2 % level of biological control. Clearly parasitism, particularly at the higher rate, helped to destabilize the D. suzukii population and is an effect that may be enhanced by additional compatible control measures.

Overall, these data suggest that biological control, as it currently stands, will not be effective as a standalone management tactic but will result in additive contributions to IPM programs targeting D. suzukii. Increased benefit will undoubtedly be gained from additional classic biological control introductions. We believe that the refined model presented here can be used as a comparative tool for practitioners and scientists, and such models will allow for the integration and optimization of multiple IPM technologies. This approach also illustrates that IPM practitioners should take advantage of environmental conditions that create vulnerability of the pest to management activities. During more suitable summer conditions, alternative factors such as pesticide use and biological control are considered as key management techniques. We realize that manipulations of D. suzukii populations are not the only factors that can be used in IPM strategies. Behavioral techniques including push–pull strategies, cultural methods, and insect barriers could also contribute to sustainable management of D. suzukii. In practice, IPM strategies often focus on single technologies including biological control, host plant resistance, chemical pesticides, or biopesticides. Typically, very little attention is given to the interaction or compatibility of the different technologies used and their timing . Future studies using stage-structured models such as the one presented here should be conducted to investigate these interactions. We realize that population modeling is only one approach for understanding how to manage this damaging pest. However, it can provide powerful insights into the relative performance of different tactics and combinations, and we expect that population modeling will allow rapid assessment of different integrated control programs and their expected performance under different environmental conditions.Cardiometabolic diseases encompass a cluster of cardiovascular, metabolic, prothrombotic, and inflammatory abnormalities that are recognized as disease states by the American Society of Endocrinology, the National Cholesterol Education Program, and the WHO . Food intake plays a key role in reducing the risk of cardiometabolic diseases, with data suggesting that >30% of all deaths could be prevented through dietary changes, particularly by increased consumption of plant based foods . Plant foods are rich sources of fiber and essential micronutrients, such as vitamins and minerals. They are also sources of a large group of bio-active compounds, which might not be essential throughout life or cause clinically manifested deficiencies, but when consumed with the diet, these phytochemicals may promote health and wellbeing in adulthood and the elderly population by reducing the risk of age-related chronic diseases . The major categories of dietary phytochemicals include polyphenols, such as flavonoids or phenolic acids, carotenoids, or plant sterols.Growing evidence from mechanistic studies, clinical trials, and prospective cohort studies suggest that these bio-active compounds may help in promoting health when consumed as part of the habitual diet. Polyphenols are secondary metabolites of plants and are found in fruits, vegetables, and their products . Carotenoids, including b-carotene, lycopene, lutein, and zeaxanthin, are phytochemicals found in many fruits and vegetables and account for the brilliant colors of these foods . Phytosterols are cholesterol-like molecules found in all plant foods, with the highest concentrations occurring in vegetable oils, but they can also be found in nuts, breads, or whole vegetables . Furthermore, caffeine ranks as one of the most commonly consumed dietary micronutrients; it is found in coffee beans, cacao beans, kola nuts, guarana berries, and tea leaves including yerba mate . Accumulating evidence from cohort studies suggests that an increased intake of polyphenols, which are the most abundant category of phytochemicals present in our foods, may reduce the risk of cardiovascular diseases . This evidence is supported by animal and clinical studies that have reported beneficial effects of the intake of polyphenol-rich foods or purified compounds on intermediate risk factors for CVD, including LDL cholesterol, blood pressure, and endothelial function . The most convincing clinical evidence for the cardioprotective benefits of the consumption of dietary polyphenols relates to their observed beneficial effect on endothelial function . For example, the consumption of suitable doses of plant sterols has repeatedly been shown in randomized controlled trials to reduce LDL cholesterol concentrations and thus reduce risk of subsequent CVD . Another phytochemical-rich source is coffee, one of the most widely consumed beverages worldwide. Coffee consumption may reduce the risk of type 2 diabetes mellitus and hypertension, as well as other conditions associated with cardiovascular risk , and epidemiological studies suggest that regular coffee drinkers have reduced mortality, predominantly as a result of their reduced risk of developing CVD .

Grapevine is the crop associated with the highest number of intracellular pathogens

Viruses belonging to the recently established genus Coguvirus also infect plants and, although this genus has not been included yet in the family Phenuiviridae, many structural and molecular features and phylogenetic relationships with other phenuiviruses support its classification in this family. With respect to nsRNA viruses, there is only one preliminary report of tomato spotted wilt virus in this host ; however, later studies failed to transmit TSWV to grapevine. During the screening of two selections of Garan dmak and Muscat rose grapevines via HTS, partial sequences with homology to apple rubbery wood viruses 1 and 2 were identified. ARWV1 and ARWV2 are two novel nsRNA viruses with tri-segmented genomes recently discovered in apple trees displaying rubbery wood symptoms in Canada. Based on sequence identity and phylogenetic analysis, ARWV1 and ARWV2 were proposed to be representative members of a new genus, tentatively named Rubodvirus, within the family Phenuiviridae. Here, we report the natural occurrence and discovery of two novel nsRNA viruses, blueberry container size named grapevine Muscat rose virus and grapevine Garan dmak virus , in grapevine plants. Even more, this is the first evidence for the natural occurrence of phenui-like viruses in grapevine.

Later, reverse transcription PCR -based assays for the specific detection of both viruses were designed to investigate the prevalence in different grapevine populations. Phylogenetics and transmission of GMRV and GGDV were also investigated.In 2015, a new selection of grapevine, cultivar Garan dmak, was received as dormant cuttings from Armenia under a USDA Animal and Plant Health Inspection Service controlled import permit for inclusion in the Foundation Plant Services collection. Cuttings were propagated under mist and later transferred to single pots; all this within an insect-proof greenhouse enclosure. Six months after bud break, leaf tissue from the Garan dmak plant was collected for HTS as part of the plant certification program at FPS. Additionally, a Muscat rose grapevine at the USDA National Clonal Germplasm Repository but originated from Argentina was sampled and analyzed by HTS during a separate study about the genetic diversity of grapevine leafroll-associated virus 3. At the time of sampling, both grapevine cultivars did not display any known symptom associated with virus infection. Total nucleic acid extracts were prepared using a MagMax Plant RNA Isolation kit as per manufacturer’s protocol, but adjusting the amount of plant material based on the following molecular analysis. For RT-PCR, 0.2 g of tissue was homogenized in 2 mL of guanidine isothiocyanate lysis buffer PVP-40 using a Homex grinder. In the case of HTS, 0.7 g of tissue was homogenized in 7 mL of guanidine isothiocyanate lysis buffer. Subsequently, the quality of TNA extracts was verified using an 18S ribosomal RNA assay.

HTS analysis was performed as described by Al Rwahnih et al.. Briefly, aliquots of TNA from source samples were subjected to rRNA depletion and complementary DNA library construction. Later, cDNA libraries were sequenced using the Illumina NextSeq 500 platform using asingle-end 75-bp regime. Illumina reads were demultiplexed and adapter trimmed prior to analysis using Illumina bcl2fastq v2.20.0.422. Trimmed reads were subsequently de novo assembled into contigs using SPAdes v3.13. Generated contigs were compared against the viral database of the National Center for Biotechnology Information using tBLASTx. Novel virus contigs were initially identified by their conserved protein domains. Open reading frames were annotated using HMMER v3.1 to look for Pfam protein domains associated with viruses infecting land plants. In the case of Phenuiviridae-like contigs, the large RNAs were annotated with the bunyavirus RNA-dependent RNA polymerase domain. The medium RNAs were annotated with the viral movement protein domain. The small RNAs were annotated with the Tenuivirus/Phlebovirus nucleocapsid protein domain. The association between the small and the medium RNAs was investigated by BLASTn sequence similarity of their 50 ends. These contig sequences were subsequently confirmed by BLASTx hits to the NCBI nucleotide database which produced top hits to different accessions of ARWV2 and ARWV1. To complete the 50 end of each RNA segment present in GMRV and GGDV, the SMARTer RACE 5 0 /3 0 Kit was employed following the instructions provided by the manufacturer. In the case of the 30 ends, the methodology described by Navarro et al. was used, which involved the addition of a poly tail to the RNA template.In the last few years, as a consequence of the increasing application of HTS, many novel nsRNA viruses have been identified, most of which are from invertebrates. Thus, the classification of these viruses has been recently reassessed, with plant-infecting viruses now classified in the order Bunyavirales, Serpentovirales and Mononegavirales .

Recently, the tentative genus Rubodvirus has been officially proposed to classify two novel nsRNA viruses from apple trees, with a suggested species demarcation criteria of <95% aa identity for the RdRp. In the present study, HTS allowed the identification of two novel plant-infecting viruses with segmented nsRNA genome, GMRV, and GGDV, which are also the first nsRNA viruses identified and transmitted in grapevine. Although these viruses infect the same host species, the aa sequence identity between the putative proteins encoded by their genomic RNAs is always below 75% , indicating that GMRV and GGDV are two different viruses. Their genomes encode proteins showing the highest sequence identity with ARWV1 and ARWV2. GGDV and GMRV also share other traits with these viruses, including the number of genomic components, limited to three nsRNAs, identical terminal nucleotides in the genomic RNAs, and the lack of any ORF coding for glycoproteins. Moreover, close phylogenetic relationships between these four viruses are supported by the ML phylogenetic trees reported here, in which, independently of the considered protein , they always clustered in the same clade, which is significantly separated from all the other nsRNA viruses included in the analyses. Altogether these data support the classification of GMRV and GGDV as two novel species in the tentative genus Rubodvirus. When the other bunyavirales are considered, the ML phylograms inferred from the RdRps or the NPs show a close phylogenetic relationship between rubodviruses and coguviruses, which together with the arthropod-infecting LLV, form a superclade nested at a basal node, in closer proximity to arthropod-infecting viruses than other plant-infecting viruses. These data are consistent with the hypothesis, previously advanced for the coguviruses, that all the members of this superclade evolved from a common ancestor virus infecting arthropods. In this evolutionary scenario, the acquisition of the MP gene appears to be the key step in the adaptation of the ancestor virus to plants, an event that likely happened through the typical modular genome evolution process proposed for most eukaryotic viruses. However, the clustering of rubodviruses and coguviruses in two distant clades, observed in the ML tree inferred with MPs , supports the independent acquisition of the MP gene by the ancestor of the viruses included in these two taxa. These data are consistent with the hypothesis that the adaptation of invertebrate-infecting nsRNA ancestor viruses to plants happened several times through independent events during the evolutionary history of nsRNA viruses infecting plants. It is worthy of note that rubodviruses, laulaviruses, and coguviruses, although phylogenetically related, have divergent genome structures and gene expression strategies. In fact, the members of the first two genera have a genome composed of three monocistronic nsRNAs encoding different proteins;RdRp, NP, and putative MP in the case of rubodviruses, and RdRp, NP and a protein of unknown function in the case of laulaviruses. Instead, coguviruses have a bipartite genome consisting of one nsRNA encoding the RdRp , and one ambisense RNA , in which the ORFs encoding the NP and the MP are separated by a long intergenic region .

It has been shown that such an IR is AU-rich and self-complementary, thus assuming in both polarity strands a compact conformation containing a long hairpin predicted to serve as a transcription termination signal during the expression of viral genes. Taking this into consideration, growing raspberries in container the question arises of how viruses with such different genomic organizations may have evolved from the same ancestor virus. In this respect, it can be speculated that a recombination event between the viral and the vc strand of two genomic RNAs with long, AU rich, and almost identical 50 UTRs could generate ambisense RNAs containing an IR similar to those of coguviruses. Since nsRNA viruses with genomic RNAs showing structural features compatible with this possibility were not known previously, such a possibility appeared unlikely. However, the very long, AU-rich, and highly conserved 50 UTRs reported here for RNA 2 and RNA 3 of GMRV and GGDV, and also observed in the corresponding RNAs of ARWV1 and ARWV2, are the first clear evidence that nsRNA viruses with the structural features compatible with this evolutionary scenario may exist. Based on these considerations, the possibility that the bipartite genomes of coguviruses originated from a tripartite ancestor with genomic RNAs containing 50 UTRs similar to those observed in the rubodviruses appears feasible. No glycoprotein is encoded by rubodviruses , a feature previously reported also for coguviruses that, according to electron microscopy observations, are flexuous, non-enveloped viruses. In contrast, glycoproteins are expressed by most nsRNA plant viruses transmitted by arthropods . The lack of glycoprotein in the genome of rubodviruses and coguviruses opens the question on the existence of vectors, if any, involved in their transmission. In this respect, it is worthy of note that ophioviruses and varicosaviruses, which are plant-restricted or transmitted by fungi, also do not code for any glycoprotein. Vegetative propagation has been proposed as the prevalent transmission mechanism for ARWV1 and ARWV2. Whether, this is also the case for GMRV and GGDV needs further investigation. Most plant viruses code for viral suppressor proteins counteracting the plant antiviral defense mechanisms based on RNA silencing. Further studies are needed to ascertain whether one or more of the three proteins encoded by GGDV and GMRV and other rubodviruses may interfere with RNA silencing, thus showing multifunctional role, as already reported for VSRs of other viruses. GMRV and GGDV were identified in two different grapevines cultivars, Muscat rose and Garan dmak, that were tested by HTS; moreover, both viruses were found in association with other viruses and viroids. Interestingly, no obvious symptoms were observed in the two grapevines. Although infectivity of both viruses was ascertained by graft-transmission, only a Cabernet franc grapevine infected by GGDV developed symptoms post-grafting; further HTS analysis on this indicator plant reveals the presence of GRSPaV, GYSVd-1, and HSVd. Therefore, it was not possible to ascertain whether GGDV is associated with symptoms and additional studies on its pathogenicity are needed, likewise GMRV. An initial survey using several accessions of grapevine located in three different collections in California resulted in the identification of two and one plants infected by GGDV and GMRV, respectively; like the original sources of the viruses , these grapevines were symptomless. A more extensive survey, including other grapevine-growing regions in California and the USA, is necessary to determine the real distribution of these novel viruses. In that sense, the detection method developed in this study could be useful for virus testing and certification programs. Finally, during the review process of this manuscript, two novel mycoviruses related to coguviruses and rubodviruses were reported, which extends the host range of phenui-like viruses.Molecular networking1 , introduced in 2012, was one of the first data organization approaches to visualize the relationships between tandem mass spectrometry fragmentation spectra. In molecular networking, relationships between similar MS/MS spectra are visualized as edges. As MS/MS spectral similarity implies chemical structural similarity , chemical structural information can thus be represented as a network and chemical relationships can be visualized. This approach forms the basis for the web-based mass spectrometry infrastructure, Global Natural Products Social Molecular Networking which sees ~200,000 new accessions per month. Molecular networking has successfully been used for a range of applications in drug discovery, natural products research, environmental monitoring, medicine, and agriculture. To tap into the chemistry of complex samples through metabolomics, a subset of MS/MS spectra can be annotated by spectral library matching or by using in silico approaches. While molecular networking facilitates the visualization of closely related molecules in molecular families, the inference of chemical relationships at a dataset-wide level and in the context of diverse sample metadata requires complementary representation strategies. To address this need, we developed an approach that uses fragmentation trees and machine learning to calculate all pairwise chemical relationships.

We assessed resource use in each site in both a control and a manipulated state

Tongue length appears to be important for structuring foraging preferences but we are lacking experimental work that evaluates how traits such as tongue length influence pollinators’ response to competition. We systematically manipulated interspecific competition in bumble bee communities through targeted single species removals and examined patterns of species-specific plasticity in resource use in the remaining pollinators. Specifically we examined to what extent tongue length explains species-specific difference in the pollinators’ foraging behavior in response to release from competition. Our sites are self assembled communities with natural diversity varying in plant and bee community composition. Utilizing removals in these natural communities allowed us to examine if the identity of the most abundant bee species, influences bee foraging behavior. We focused on the foraging response of the remaining bees in the community with respect to floral fidelity, or, within plant species movements within a single foraging bout. Floral fidelity is critical for many plants species’ reproductive success because transfer of conspecific pollen must occur in order for fertilization to take place.Study Sites: We worked in 28 subalpine meadow sites in the landscape surrounding the Rocky Mountain Biological Laboratory , plastic potting pots in the Gunnison National Forest, western Colorado, United States.

Eachsite consisted of a 20 × 20-m plot, all with the same dominant plant species to minimize, as much as possible, plant-community-driven differences in foraging behavior. A minimum distance of 1 km separated any two sites. We collected data over three summer growing seasons , in 2010, 2011 and 2013. Manipulations. We assessed each plot in a control state, waited one day, and then assessed in a manipulated state. We kept the interval between control and manipulated states short because of the rapid turnover in flowering phenology in our high-altitude system, allowing us to keep the plant community constant in our control–manipulation comparisons . Manipulations reduce interspecific competition through the temporary, nondestructive removal of the most abundant bumblebee species in each plot. We determined the most abundant bee via inventory of Bombus species richness and abundance on the control day using nondestructive aerial netting, with two field team members netting for a 20-min period, not including handling time . To avoid double-counting, we kept each bee in an individual glass vial, identified to species, and kept in a cool, dark cooler until the inventory time period was over, at which point bees were released. On the manipulation day we removed the most abundant bee species . The removals were accomplished through targeted hand-netting, and we minimized disturbance of other bees and vegetation by carefully placing the insect net over entire inflorescence and allowing bee individuals to fly up into the net . Captured bees were transferred to vials and placed in a cooler during the manipulation and released unharmed afterward.

We used as much time as necessary to remove essentially all individuals of the target species from the sample plot and immediately adjacent area . We left a period of at least 30 minutes between manipulative bee removals and subsequent sampling to minimize the impact of the disturbance on the foraging activities of other bees. We recorded both the abundance of removed individuals, as well as the number of un-captured “escapees” that were observed during bee sampling. Each site was only used once in a manipulated and controlled state per year . Foraging observations: We directly followed the foraging sequences of Bombus individuals in both the control and manipulated states. We recorded the identity of each plant species visited in a foraging sequence. We discontinued an observation when the bee was lost from sight, when it ventured more than 5 m outside of the plot, when it had been observed for 10 full minutes, or when we had tallied 100 individual plants visited. We discarded observations of bees that visited fewer than five plants. The number of individuals observed per/state/site varied due to bee abundance .Tongue Length measurements: The data on proboscis lengths of workers was taken from published measurements of bumblebees collected on the Front Range of the Colorado Rocky Mountains and overlap with the sites and species that we used for this study . Measurements indicate the sum of the individual lengths of the prementum and glossa. The mean tongue length of each bumble bee species was assigned to each individual bee and used in the trait analysis .

Floral fidelity was measured as the binomial counts of individual bee foraging movements that were conspecific vs. heterospecific . We used GLMMs with binomial errors using the logit link in the lme4 package for R to model the floral fidelity response variable . Data from individual bees foraging within a site cannot be considered independent and therefore we used site as a random effect . Relative to a binomial distribution, our data were overdispersed, which we corrected by including an individual-level random effect . We used the R statistical programming language for all models.Our results demonstrate that bees vary in their floral fidelity and that tongue length explains a large part of this variation. Bees with shorter tongues move between plant species more often than bees with longer tongues. We did not find significant variation in the response of bee species to a reduction in interspecific competition, but rather saw a guild-wide reduction in floral fidelity in response to the removal of the dominant bee species . Finally, our results suggest that tongue length of the most abundant bee species, a site-level attribute, explains much of the site-to-site variation in pollinator foraging behavior. In particular, we found that as the tongue length of the most abundant bee increases, the site level foraging fidelity decreases. We found that bumble bee species vary in the degree to which they move between different plant species within a single foraging bout, and tongue length explains much of the variation. Some suggest that long tongued bees should exhibit broader resource usage patterns because their traits permit them access to a wider range of flower types . In contrast, short-tongued bees should act as specialists, with a more restricted range of resource use options, rarely able to access the nectar at the base of the flowers with long corollas . Our results suggest the opposite pattern, that shorter tongue bees are more labile with their foraging patterns and on average move between plant species within a single foraging bout more often than longer tongue bees. We suggest the following interpretation, because long tongues enable bees to access flowers with better rewards and maintain a monopoly on those rewards, they may have less incentive than short-tongued species to move between plant species while foraging. While longer tongue bumble bees are capable of foraging on flowers with short corollas it would provide less energetic gain, potentially making behavioral plasticity less profitable . In contrast, the shorter-tongued bees in our system tend to have smaller bodies and are more likely to depend on resources within a more restricted foraging range . These limitations could favor a labile foraging habit, with shorter tongue bees constantly assessing the resource availability and competitive context in their community. As such, shorter-tongued bees more readily switch between plant species. While we saw little variation in species-specific foraging responses to our manipulations, we found an overall reduction in floral fidelity across sites after the removal of the most abundant bee species. Our results build on the findings of Brosiand Briggs , reaffirming a guild-wide reduction in floral fidelity in response to a reduction in interspecific competition. This study adds an additional two years, 8 sites and 165 bee individuals to our previous study, confirming that the guild-wide results found in Brosi and Briggs are robust. Variation in bumble bee floral fidelity is largely explained by the tongue length of the most abundant bee species in each site. This means that pollinator foraging behavior is context dependent and is determined by the most abundant bee species. In general, raspberry container growing short tongued bees exhibit lower floral fidelity than long tongued bees, and when they are in a site that has a long tongued bee removal, their reduction in floral fidelity is magnified. Bumble bees are large bodied insects that require many floral resources to keep their hive growing throughout the growing season.

As such, we might expect strong competition between these species, and tongue length, arguably the trait most relevant for resource acquisition, could dictate how resources are partitioned within a community, ultimately drive the assembly of bumble bees within communities . Pyke proposed that bumble-bee species with similar tongue lengths could not exist in his altitudinal alpine transects presumably because the bees compete for floral resources. But later studies did not support this pattern , and, as in our study, found that bumble bee species with similar tongue lengths co-ocurred within a community. This has left researchers to wonder if the coexistence of many bee species with substantial overlap in their life history requirements is possible because bumble bee species compete for something other than flower resources allowing so many similar species to co-occur . Our study suggests that in our system, bumble bees do in fact compete for floral resources and that longer tongue bees seem to elicit competition that is experienced across the range of trait values seen in our sites . The willingness of short tongue bees to exhibit behavioral plasticity may allow for such a large number of seemingly similar bee species to coexist in a community. Future work should examine the extent to which this plasticity is adaptive and assess the fitness costs that may result from the willingness to switch floral resources in response to a reduction in interspecific competition. Most pollinators are generalist foragers that can switch between plant species within a single foraging bout . When pollinators move between plant species, they can transfer heterospecific pollen to plant stigmas which in turn can reduce plant reproduction . We found an effect in which competition from a long tongued bee changes the foraging behavior of the rest of the bees in a way that could be detrimental to plant reproduction. From a plant’s perspective, not only do short tongue bees exhibit behavior that likely results in the transfer of heterospecific pollen, but when short tongue bees are in communities in which a longer tongue bee is most abundant, they exhibit even greater floral infidelity, making the likelihood of heterospecific pollen deposition even greater .Pollinator species are on the decline globally . Bumble bees in particular are experiencing population range contractions due to climate change as well reductions in abundance due to disease , agricultural intensification and pesticide use . In a changing world where we are likely to experience an emergence of new interactions via range shifts, introduced species and climate change, exploring how the competitive landscape shapes foraging plasticity will help us generalize to other plant pollinator systems and begin to better predict the functional implications of competitive interactions.Pollinator losses are increasing across the globe which could have potentially strong negative effects on the plants that rely on them for pollination . Network-based simulations suggest that plant communities will be very robust to pollinator extinctions . This robustness is likely driven by two features of network structure. First, pollination networks are dominated by generalist interactions; most plants and pollinators each interact with several species from the other group over the course of their lives . Second, pollination networks have a nested structure in which specialist plants and pollinators tend to interact with a subset of the species in the other group that generalists interact with . Nestedness leads to an asymmetric interaction structure where specialists from one group tend to act with generalists, not specialists, from the other, which could reduce linked extinctions if specialists are vulnerable to stochastic extinctions . One feature of binary-graph plant-pollinator simulation models that may overestimate network robustness is that all interactions in a binary-graph network are positive. A typical binary-graph simulation modeling approach is that if at least one link remains between a plant and a pollinator, the plant will continue to persist . While the assumption of “all interactions positive” is a reasonable starting point in a mutualistic network model, empirical evidence suggests at least two ways in which pollinators can have negative consequences for the reproduction of plants they interact with . First, the benefit that pollinators have on plant reproduction is sensitive to how “faithful” pollinators are to particular plant species in a single foraging bout.