Spiders (Araneae) in the biological and integrated pest management of apple in the Czech Republic

Abstract:  Two management systems, biological and integrated, were compared to control the major pest, codling moth ( Cydia pomonella ) in apple orchards. The aim of the study was to assess the effect of these two systems on arboreal spiders. The biological system was based on the use of biological preparations, whereas in the integrated system selective pesticides were employed. The control plot had no pesticide treatment. The abundance of spiders was similar in all study plots during 3 years of study. Diversity, however, was higher in the biological plot than in the control and the integrated plots, suggesting that the response of spiders to management was guild-specific. Four spider families dominanted in all plots: Araneidae (orb weavers), Theridiidae (space-web weavers), and Philodromidae and Thomisidae (ambushers). While Araneidae and Thomisidae were similarly abundant on all plots, the density of Theridiidae and Philodromidae differed. In the integrated plot there were significantly more theridiid spiders, whereas in the control plot philodromid spiders were significantly more abundant. On the biological plot, the two families were similarly abundant. These differences were attributed to different age of trees in the control and treated plots, different prey spectrum, different susceptibility of the two families to applied chemicals and intraguild predation of theridiids by philodromids.     

Cryptic genetic variation in an inbreeding and cosmopolitan pest,Xylosandrus crassiusculus,revealed using ddRADseq

Each year new exotic species are transported across the world through global commerce, causing considerable economic and ecological damage. An important component of managing invasion pathways is to identify source populations. Some of the most widespread exotic species are haplodiploid ambrosia beetles. The ability to mate with siblings (inbreed) and their transportable food source (symbiotic fungus) have enabled them to colonize most of the world and become pests of plant nurseries, lumber, and forests. One of the fastest spreading ambrosia beetles is Xylosandrus crassiusculus. In order to discover the source populations of this globally invasive species, track its movement around the world, and test biogeographical scenarios, we combined restriction site‐associated DNA sequencing (RADseq) with comprehensive sampling across the species native and introduced range. From 1,365 genotyped SNP loci across 198 individuals, we determined that in its native range, X. crassiusculus is comprised of a population in Southeast Asia that includes mainland China, Thailand, and Taiwan, and a second island population in Japan. North America and Central America were colonized from the island populations, while Africa and Oceania were colonized from the mainland Asia, and Hawaii was colonized by both populations. Populations of X. crassiusculus in North America were genetically diverse and highly structured, suggesting (1) numerous, repeated introductions; (2) introduction of a large founding population; or (3) both scenarios with higher than expected outcrossing. X. crassiusculus, other wood‐boring insects, and indeed many other pests with unusual genetic structure continue to spread around the world. We show that contemporary genetic methods offer a powerful tool for understanding and preventing pathways of future biosecurity threats.     

SNP‐based genetic characterization of the Tulane National Primate Research Center's conventional and specific pathogen‐free rhesus macaque (Macaca mulatta) populations

Background

The rhesus macaque is an important biomedical model organism, and the Tulane National Primate Research Center (TNPRC) has one of the largest rhesus macaque breeding colonies in the United States.

Methods

SNP profiles from 3266 rhesus macaques were used to examine the TNPRC colony genetic composition over time and across conventional or SPF animals of Chinese and Indian ancestry.

Results

Chinese origin animals were the least genetically diverse and the most inbred; however, since their derivation from their conventional forebearers, neither the Chinese nor the Indian SPF animals exhibit any significant loss of genetic diversity or differentiation.

Conclusions

The TNPRC colony managers have successfully minimized loss in genetic variation across generations. Although founder effects and bottlenecks among the Indian animals have been successfully curtailed, the Chinese subpopulation still show some influences from these events.     

Morphological,genetic and symptomatic identification of an invasive jujube pest in Korea,Dasineura jujubifolia Jiao & Bu (Diptera: Cecidomyiidae)

Correct identification of invasive species is an important step for further management of this species. New type of leafcurling gall midge feeding on jujube, Ziziphus jujube Miller (Rhamnaceae), was discovered in Korea in 2011. At that time, this species was not identified as specific level, so it was recorded as Dasineura sp. Since then, the jujube gall midge has become a serious alien insect pest in Korea. In this study, this jujube gall midge was surveyed by collecting damaged leaf samples from different major jujube-producing regions in Korea. Morphological, genetic and symptomatic identification of this invasive jujube pest identified as Dasineura jujubifolia Jiao & Bu. Here we provided its taxonomic status, geographical distribution, morphological characteristics of all stages and infestation symptoms. In addition to classical morphological characters, mitochondrial COI barcoding sequences were generated for several specimens. The possible pathway of invasion and subsequent socioeconomic consequences were discussed.     

RAD genotyping reveals fine‐scale genetic structuring and provides powerful population assignment in a widely distributed marine species,the American lobster (Homarus americanus)

Deciphering genetic structure and inferring connectivity in marine species have been challenging due to weak genetic differentiation and limited resolution offered by traditional genotypic methods. The main goal of this study was to assess how a population genomics framework could help delineate the genetic structure of the American lobster (Homarus americanus) throughout much of the species’ range and increase the assignment success of individuals to their location of origin. We genotyped 10 156 filtered SNPs using RAD sequencing to delineate genetic structure and perform population assignment for 586 American lobsters collected in 17 locations distributed across a large portion of the species’ natural distribution range. Our results revealed the existence of a hierarchical genetic structure, first separating lobsters from the northern and southern part of the range (F_CT = 0.0011; P‐value = 0.0002) and then revealing a total of 11 genetically distinguishable populations (mean F_ST = 0.00185; CI: 0.0007–0.0021, P‐value < 0.0002), providing strong evidence for weak, albeit fine‐scale population structuring within each region. A resampling procedure showed that assignment success was highest with a subset of 3000 SNPs having the highest F_ST. Applying Anderson's (Molecular Ecology Resources, 2010, 10, 701) method to avoid ‘high‐grading bias’, 94.2% and 80.8% of individuals were correctly assigned to their region and location of origin, respectively. Lastly, we showed that assignment success was positively associated with sample size. These results demonstrate that using a large number of SNPs improves fine‐scale population structure delineation and population assignment success in a context of weak genetic structure. We discuss the implications of these findings for the conservation and management of highly connected marine species, particularly regarding the geographic scale of demographic independence.     

Population structure and gene flow in the global pest,Helicoverpa armigera

Helicoverpa armigera is a major agricultural pest that is distributed across Europe, Asia, Africa and Australasia. This species is hypothesized to have spread to the Americas 1.5 million years ago, founding a population that is at present, a distinct species, Helicoverpa zea. In 2013, H. armigera was confirmed to have re‐entered South America via Brazil and subsequently spread. The source of the recent incursion is unknown and population structure in H. armigera is poorly resolved, but a basic understanding would highlight potential biosecurity failures and determine the recent evolutionary history of region‐specific lineages. Here, we integrate several end points derived from high‐throughput sequencing to assess gene flow in H. armigera and H. zea from populations across six continents. We first assemble mitochondrial genomes to demonstrate the phylogenetic relationship of H. armigera with other Heliothine species and the lack of distinction between populations. We subsequently use de novo genotyping‐by‐sequencing and whole‐genome sequences aligned to bacterial artificial chromosomes, to assess levels of admixture. Primarily, we find that Brazilian H. armigera are derived from diverse source populations, with strong signals of gene flow from European populations, as well as prevalent signals of Asian and African ancestry. We also demonstrate a potential field‐caught hybrid between H. armigera and H. zea, and are able to provide genomic support for the presence of the H. armigera conferta subspecies in Australasia. While structure among the bulk of populations remains unresolved, we present distinctions that are pertinent to future investigations as well as to the biosecurity threat posed by H. armigera.     

Population genomics of wild and laboratory zebrafish (Danio rerio)

Understanding a wider range of genotype–phenotype associations can be achieved through ecological and evolutionary studies of traditional laboratory models. Here, we conducted the first large‐scale geographic analysis of genetic variation within and among wild zebrafish (Danio rerio) populations occurring in Nepal, India, and Bangladesh, and we genetically compared wild populations to several commonly used lab strains. We examined genetic variation at 1832 polymorphic EST‐based single nucleotide polymorphisms (SNPs) and the cytb mitochondrial gene in 13 wild populations and three lab strains. Natural populations were subdivided into three major mitochondrial DNA clades with an average among‐clade sequence divergence of 5.8%. SNPs revealed five major evolutionarily and genetically distinct groups with an overall F_ST of 0.170 (95% CI 0.105–0.254). These genetic groups corresponded to discrete geographic regions and appear to reflect isolation in refugia during past climate cycles. We detected 71 significantly divergent outlier loci (3.4%) and nine loci (0.5%) with significantly low F_ST values. Valleys of reduced heterozygosity, consistent with selective sweeps, surrounded six of the 71 outliers (8.5%). The lab strains formed two additional groups that were genetically distinct from all wild populations. An additional subset of outlier loci was consistent with domestication selection within lab strains. Substantial genetic variation that exists in zebrafish as a whole is missing from lab strains that we analysed. A combination of laboratory and field studies that incorporates genetic variation from divergent wild populations along with the wealth of molecular information available for this model organism provides an opportunity to advance our understanding of genetic influences on phenotypic variation for a vertebrate species.     

Historical demography and genetic differentiation of the giant freshwater prawn Macrobrachium rosenbergii in Bangladesh based on mitochondrial and ddRAD sequence variation

Macrobrachium rosenbergii, the giant freshwater prawn, is an important source of high quality protein and occurs naturally in rivers as well as commercial farms in South and South‐East Asia, including Bangladesh. This study investigated the genetic variation and population structure of M. rosenbergii sampled from four rivers in Bangladesh (sample size ranged from 19 to 20), assessing sequence variation, both in the mitochondrial cytochrome oxidase subunit 1 (CO1) gene and in 106 single nucleotide polymorphisms (SNPs) sampled randomly from the genome with double digest RAD sequencing (ddRADseq). The mitochondrial variation presented a shallow genealogy with high haplotype diversity (h = 0.95), reflecting an expansion in population size for the last ~82 kyr. Based on the CO1 variation the current effective population size (N_e) was 9.7 × 10⁶ (CI: 1.33 × 10⁶ – 35.84 × 10⁶) individuals. A significant population differentiation was observed with the mitochondrial CO1 sequence variation and based on the ddRADseq variation, which could be traced to the divergence of the population in the Naf River in the South‐East border with Myanmar from the other populations. A differentiation in mtDNA haplotype frequencies was also observed between the Biskhali River and the Karnaphuli Rivers in eastern Bangladesh. This study demonstrated the use of high‐throughput genotyping based on the ddRADseq method to reveal population structure at a small geographical scale for an important freshwater prawn. The information from this study can be utilized for management and conservation of this species in Bangladesh.     

Dissecting the association between a gall midge,Asteromyia carbonifera,and its symbiotic fungus,Botryosphaeria dothidea

The Ambrosia gall midge [Asteromyia carbonifera (Osten Sacken) (Diptera: Cecidomyiidae: Alycaulini)] consists, in part, of a complex of genetically differentiated populations that have diverged in gall morphology on the host plant Solidago altissima L. (Asteraceae). This divergence appears to be an incipient adaptive radiation that may be driven by parasitoid pressure. Understanding the mechanisms driving this genetic and phenotypic diversification requires a close examination of the relationship between the midge and its fungal associate Botryosphaeria dothidea (Moug.) Ces. & De Not. (Ascomycota: Dothideomycetes), whose mycelia actually form the protective gall structure. We used manipulative experiments to test the degree of interdependency of the fungus and the midge, and we employed field and laboratory studies to gain insight into the source of fungal conidia, which our data and observations indicate are collected by females and stored in specialized pockets (mycangia) on the ovipositor. Manipulative experiments demonstrate that fungal proliferation on the host plant is dependent on the midge larvae and larvae exhibit significant growth on the fungus alone. Field observations and experiments were unable to identify the source of mycangial conidia; however, analyses of conidia shape suggest a biotrophic source. We conclude that this association is an obligatory mutualism with respect to successful gall formation. These findings corroborate recent findings that the primary food source of the midge is the gall fungus.     

Secondary contact and changes in coastal habitat availability influence the nonequilibrium population structure of a salmonid (Oncorhynchus keta)

Numerous empirical studies have reported lack of migration–drift equilibrium in wild populations. Determining the causes of nonequilibrium population structure is challenging because different evolutionary processes acting at a variety of spatiotemporal scales can produce similar patterns. Studies of contemporary populations in northern latitudes suggest that nonequilibrium population structure is probably caused by recent colonization of the region after the last Pleistocene ice age ended ~13 000 years ago. The chum salmon's (Oncorhynchus keta) range was fragmented by dramatic environmental changes during the Pleistocene. We investigated the population structure of chum salmon on the North Alaska Peninsula (NAP) and, using both empirical data and simulations, evaluated the effects of colonization timing and founder population heterogeneity on patterns of genetic differentiation. We screened 161 single nucleotide polymorphisms and found evidence of nonequilibrium population structure when the slope of the isolation‐by‐distance relationship was examined at incremental spatial scales. In addition, simulations suggested that this pattern closely matched models of recent colonization of the NAP by secondary contact. Our results agree with geological and archaeological data indicating that the NAP was a dynamic landscape that may have been more recently colonized than during the last deglaciation because of dramatic changes in coastal hydrology over the last several thousand years.     

Genome‐wide SNP analysis unveils genetic structure and phylogeographic history of snow sheep (Ovis nivicola) populations inhabiting the Verkhoyansk Mountains and Momsky Ridge (northeastern Siberia)

Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise F_ST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.     

Revisiting comparisons of genetic diversity in stable and declining species: assessing genome‐wide polymorphism in North American bumble bees using RAD sequencing

Genetic variation is of key importance for a species’ evolutionary potential, and its estimation is a major component of conservation studies. New DNA sequencing technologies have enabled the analysis of large portions of the genome in nonmodel species, promising highly accurate estimates of such population genetic parameters. Restriction site‐associated DNA sequencing (RADseq) is used to analyse thousands of variants in the bumble bee species Bombus impatiens, which is common, and Bombus pensylvanicus, which is in decline. Previous microsatellite‐based analyses have shown that gene diversity is lower in the declining B. pensylvanicus than in B. impatiens. RADseq nucleotide diversities appear much more similar in the two species. Both species exhibit allele frequencies consistent with historical population expansions. Differences in diversity observed at microsatellites thus do not appear to have arisen from long‐term differences in population size and are either recent in origin or may result from mutational processes. Additional research is needed to explain these discrepancies and to investigate the best ways to integrate next‐generation sequencing data and more traditional molecular markers in studies of genetic diversity.     

Assessing polar bear (Ursus maritimus) population structure in the Hudson Bay region using SNPs

Defining subpopulations using genetics has traditionally used data from microsatellite markers to investigate population structure; however, single‐nucleotide polymorphisms (SNPs) have emerged as a tool for detection of fine‐scale structure. In Hudson Bay, Canada, three polar bear (Ursus maritimus) subpopulations (Foxe Basin (FB), Southern Hudson Bay (SH), and Western Hudson Bay (WH)) have been delineated based on mark–recapture studies, radiotelemetry and satellite telemetry, return of marked animals in the subsistence harvest, and population genetics using microsatellites. We used SNPs to detect fine‐scale population structure in polar bears from the Hudson Bay region and compared our results to the current designations using 414 individuals genotyped at 2,603 SNPs. Analyses based on discriminant analysis of principal components (DAPC) and STRUCTURE support the presence of four genetic clusters: (i) Western—including individuals sampled in WH, SH (excluding Akimiski Island in James Bay), and southern FB (south of Southampton Island); (ii) Northern—individuals sampled in northern FB (Baffin Island) and Davis Strait (DS) (Labrador coast); (iii) Southeast—individuals from SH (Akimiski Island in James Bay); and (iv) Northeast—individuals from DS (Baffin Island). Population structure differed from microsatellite studies and current management designations demonstrating the value of using SNPs for fine‐scale population delineation in polar bears.     

Genome-wide population structure and evolutionary history of the Frizarta dairy sheep

In the present study, we used genomic data, generated with a medium density single nucleotide polymorphisms (SNP) array, to acquire more information on the population structure and evolutionary history of the synthetic Frizarta dairy sheep. First, two typical measures of linkage disequilibrium (LD) were estimated at various physical distances that were then used to make inferences on the effective population size at key past time points. Population structure was also assessed by both multidimensional scaling analysis and k-means clustering on the distance matrix obtained from the animals’ genomic relationships. The Wright’s fixation F_ST index was also employed to assess herds’ genetic homogeneity and to indirectly estimate past migration rates. The Wright’s fixation F_IS index and genomic inbreeding coefficients based on the genomic relationship matrix as well as on runs of homozygosity were also estimated. The Frizarta breed displays relatively low LD levels with r² and |Dʹ| equal to 0.18 and 0.50, respectively, at an average inter-marker distance of 31 kb. Linkage disequilibrium decayed rapidly by distance and persisted over just a few thousand base pairs. Rate of LD decay (β) varied widely among the 26 autosomes with larger values estimated for shorter chromosomes (e.g. β=0.057, for OAR6) and smaller values for longer ones (e.g. β=0.022, for OAR2). The inferred effective population size at the beginning of the breed’s formation was as high as 549, was then reduced to 463 in 1981 (end of the breed’s formation) and further declined to 187, one generation ago. Multidimensional scaling analysis and k-means clustering suggested a genetically homogenous population, F_ST estimates indicated relatively low genetic differentiation between herds, whereas a heat map of the animals’ genomic kinship relationships revealed a stratified population, at a herd level. Estimates of genomic inbreeding coefficients suggested that most recent parental relatedness may have been a major determinant of the current effective population size. A denser than the 50k SNP panel may be more beneficial when performing genome wide association studies in the breed.     

Global lack of flyway structure in a cosmopolitan bird revealed by a genome wide survey of single nucleotide polymorphisms

Knowledge about population structure and connectivity of waterfowl species, especially mallards (Anas platyrhynchos), is a priority because of recent outbreaks of avian influenza. Ringing studies that trace large‐scale movement patterns have to date been unable to detect clearly delineated mallard populations. We employed 363 single nucleotide polymorphism markers in combination with population genetics and phylogeographical approaches to conduct a population genomic test of panmixia in 801 mallards from 45 locations worldwide. Basic population genetic and phylogenetic methods suggest no or very little population structure on continental scales. Nor could individual‐based structuring algorithms discern geographical structuring. Model‐based coalescent analyses for testing models of population structure pointed to strong genetic connectivity among the world's mallard population. These diverse approaches all support the conclusion that there is a lack of clear population structure, suggesting that the world's mallards, perhaps with minor exceptions, form a single large, mainly interbreeding population.     

Genomic distribution and estimation of nucleotide diversity in natural populations: perspectives from the collared flycatcher (Ficedula albicollis) genome

Properly estimating genetic diversity in populations of nonmodel species requires a basic understanding of how diversity is distributed across the genome and among individuals. To this end, we analysed whole‐genome resequencing data from 20 collared flycatchers (genome size ≈1.1 Gb; 10.13 million single nucleotide polymorphisms detected). Genomewide nucleotide diversity was almost identical among individuals (mean = 0.00394, range = 0.00384–0.00401), but diversity levels varied extensively across the genome (95% confidence interval for 200‐kb windows = 0.0013–0.0053). Diversity was related to selective constraint such that in comparison with intergenic DNA, diversity at fourfold degenerate sites was reduced to 85%, 3′ UTRs to 82%, 5′ UTRs to 70% and nondegenerate sites to 12%. There was a strong positive correlation between diversity and chromosome size, probably driven by a higher density of targets for selection on smaller chromosomes increasing the diversity‐reducing effect of linked selection. Simulations exploring the ability of sequence data from a small number of genetic markers to capture the observed diversity clearly demonstrated that diversity estimation from finite sampling of such data is bound to be associated with large confidence intervals. Nevertheless, we show that precision in diversity estimation in large outbred population benefits from increasing the number of loci rather than the number of individuals. Simulations mimicking RAD sequencing showed that this approach gives accurate estimates of genomewide diversity. Based on the patterns of observed diversity and the performed simulations, we provide broad recommendations for how genetic diversity should be estimated in natural populations.     

Clones or clans: the genetic structure of a deep‐sea sponge,Aphrocallistes vastus,in unique sponge reefs of British Columbia,Canada

Understanding patterns of reproduction, dispersal and recruitment in deep‐sea communities is increasingly important with the need to manage resource extraction and conserve species diversity. Glass sponges are usually found in deep water (>1000 m) worldwide but form kilometre‐long reefs on the continental shelf of British Columbia and Alaska that are under threat from trawling and resource exploration. Due to their deep‐water habitat, larvae have not yet been found and the level of genetic connectivity between reefs and nonreef communities is unknown. The genetic structure of Aphrocallistes vastus, the primary reef‐building species in the Strait of Georgia (SoG) British Columbia, was studied using single nucleotide polymorphisms (SNPs). Pairwise comparisons of multilocus genotypes were used to assess whether sexual reproduction is common. Structure was examined 1) between individuals in reefs, 2) between reefs and 3) between sites in and outside the SoG. Sixty‐seven SNPs were genotyped in 91 samples from areas in and around the SoG, including four sponge reefs and nearby nonreef sites. The results show that sponge reefs are formed through sexual reproduction. Within a reef and across the SoG basin, the genetic distance between individuals does not vary with geographic distance (r = ?0.005 to 0.014), but populations within the SoG basin are genetically distinct from populations in Barkley Sound, on the west coast of Vancouver Island. Population structure was seen across all sample sites (global F_ST = 0.248), especially between SoG and non‐SoG locations (average pairwise F_ST = 0.251). Our results suggest that genetic mixing occurs across sponge reefs via larvae that disperse widely.