RESUMEN

BACKGROUND: Microbiota alterations are linked with gastric cancer (GC). However, the relationship between the oral microbiota (especially oral fungi) and GC is not known. In this study, we aimed to apply 2b-RAD sequencing for Microbiome (2b-RAD-M) to characterize the oral microbiota in patients with GC. METHODS: We performed 2b-RAD-M analysis on the saliva and tongue coating of GC patients and healthy controls. We carried out diversity, relative abundance, and composition analyses of saliva and tongue coating bacteria and fungi in the two groups. In addition, indicator analysis, the Gini index, and the mean decrease accuracy were used to identify oral fungal indicators of GC. RESULTS: In this study, fungal imbalance in the saliva and tongue coating was observed in the GC group. At the species level, enriched Malassezia globosa (M. globosa) and decreased Saccharomyces cerevisiae (S. cerevisiae) were observed in saliva and tongue coating samples of the GC group. Random forest analysis indicated that M. globosa in saliva and tongue coating samples could serve as biomarkers to diagnose GC. The Gini index and mean decreases in accuracy for M. globosa in saliva and tongue coating samples were the largest. In addition, M. globosa in saliva and tongue coating samples classified GC from the control with areas under the receiver operating curve (AUCs) of 0.976 and 0.846, respectively. Further ecological analysis revealed correlations between oral bacteria and fungi. CONCLUSION: For the first time, our data suggested that changes in oral fungi between GC patients and controls may help deepen our understanding of the complex spectrum of the different microbiotas involved in GC development. Although the cohort size was small, this study is the first to use 2b-RAD-M to reveal that oral M. globosa can be a fungal biomarker for detecting GC.

Asunto(s)

RESUMEN

Coilia nasus is a threatened migratory species in the Yangtze River Basin. To reveal the genetic diversity of natural and farmed populations of C. nasus and the status of germplasm resources in the Yangtze River, the genetic diversity and structure of two wild populations (Yezhi Lake: YZ; Poyang Lake: PY) and two farmed populations (Zhenjiang: ZJ; Wuhan: WH) of C. nasus were analyzed using 44,718 SNPs obtained via 2b-RAD sequencing. The results indicate that both the wild and farmed populations had low genetic diversity, and germplasm resources have undergone varying degrees of degradation. Population genetic structure analyses indicated that the four populations may have come from two ancestral groups. Different amounts of gene flow were identified among WH, ZJ, and PY populations, but gene flow among YZ and other populations was low. It is speculated that the river-lake isolation of Yezhi Lake is the main cause of this phenomenon. In conclusion, this study revealed that genetic diversity reduction and germplasm resource degradation had occurred in both wild and farmed C. nasus, suggesting that conservation of its resources is of great urgency. This study provides a theoretical basis for the conservation and rational exploitation of germplasm resources for C. nasus.

RESUMEN

Animal pollinators directly affect plant gene flow by transferring pollen grains between individuals. Pollinators with restricted mobility are predicted to limit gene flow within and among populations, whereas pollinators that fly longer distances are likely to promote genetic cohesion. These predictions, however, remain poorly tested. We examined population genetic structure and fine-scale spatial genetic structure (FSGS) in six perennial understory angiosperms in Andean cloud forests of northwestern Ecuador. Species belong to three families (Gesneriaceae, Melastomataceae, and Rubiaceae), and within each family we paired one insect-pollinated with one hummingbird-pollinated species, predicting that insect-pollinated species have greater population differentiation (as quantified with the FST statistic) and stronger FSGS (as quantified with the SP statistic) than hummingbird-pollinated species. We confirmed putative pollinators through a literature review and fieldwork, and inferred population genetic parameters with a genome-wide genotyping approach. In two of the three species pairs, insect-pollinated species had much greater (>2-fold) population-level genetic differentiation and correspondingly steeper declines in fine-scale genetic relatedness. In the Gesneriaceae pair, however, FST and SP values were similar between species and to those of the other hummingbird-pollinated plants. In this pair, the insect pollinators are euglossine bees (as opposed to small bees and flies in the other pairs), which are thought to forage over large areas, and therefore may provide similar levels of gene flow as hummingbirds. Overall, our results shed light on how different animal pollination modes influence the spatial scale of plant gene flow, suggesting that small insects strongly decrease genetic cohesion.

Asunto(s)

RESUMEN

The presence of Listeria monocytogenes (Lm) in the food processing environment (facilities and products) is a challenging problem in food safety management. Lm is one of the main causes of mortality in foodborne infections, and the trend is continuously increasing. In this study, a collection of 323 Lm strain isolates recovered from food matrices and food industry environments (surfaces and equipment) over four years from 80 food processing facilities was screened using a restriction site-associated tag sequencing (2b-RAD) typing approach developed for Lm. Thirty-six different restriction site-associated DNA (RAD) types (RTs) were identified, most of which correspond to lineage II. RT1, the most represented genotype in our collection and already reported as one of the most prevalent genotypes in the food environment, was significantly associated with meat processing facilities. The sequencing of the genomes of strains belonging to the same RT and isolated in the same facility in different years revealed several clusters of persistence. The definition of the persistent strains (PSs) allowed the identification of the potential source of contamination in the incoming raw meat that is introduced in the facility to be processed. The slaughterhouses, which, according to the European Union (EU) regulation, are not inspected for the presence of Lm could be hotspots for the persistence of Lm PSs.

Asunto(s)

RESUMEN

Certain species within the genus Pseudo-nitzschia are able to produce the neurotoxin domoic acid (DA), which can cause illness in humans, mass-mortality of marine animals, and closure of commercial and recreational shellfisheries during toxic events. Understanding and forecasting blooms of these harmful species is a primary management goal. However, accurately predicting the onset and severity of bloom events remains difficult, in part because the underlying drivers of bloom formation have not been fully resolved. Furthermore, Pseudo-nitzschia species often co-occur, and recent work suggests that the genetic composition of a Pseudo-nitzschia bloom may be a better predictor of toxicity than prevailing environmental conditions. We developed a novel next-generation sequencing assay using restriction site-associated DNA (2b-RAD) genotyping and applied it to mock Pseudo-nitzschia communities generated by mixing cultures of different species in known abundances. On average, 94% of the variance in observed species abundance was explained by the expected abundance. In addition, the false positive rate was low (0.45% on average) and unrelated to read depth, and false negatives were never observed. Application of this method to environmental DNA samples collected during natural Pseudo-nitzschia spp. bloom events in Southern California revealed that increases in DA were associated with increases in the relative abundance of P. australis. Although the absolute correlation across time-points was weak, an independent species fingerprinting assay (Automated Ribosomal Intergenic Spacer Analysis) supported this and identified other potentially toxic species. Finally, we assessed population-level genomic variation by mining SNPs from the environmental 2bRAD dataset. Consistent shifts in allele frequencies in P. pungens and P. subpacifica were detected between high and low DA years, suggesting that different intraspecific variants may be associated with prevailing environmental conditions or the presence of DA. Taken together, this method presents a potentially cost-effective and high-throughput approach for studies aiming to evaluate both population and species dynamics in mixed samples.

Asunto(s)

RESUMEN

The aim of this study was to assess the current status of the germplasm resources of golden-backed carp (Cyprinus carpio var. Jinbei) cultured in paddy fields in Guizhou Province, China. Five populations of golden-backed carp in Liping County, Jinping County, Huangping County, Congjiang County and Duyun City in Guizhou Province were subjected to high-throughput sequencing by 2b-RAD technology, and their genetic diversity and genetic differentiation were analysed. Based on sequencing, 44,896 SNP loci were obtained, and all five population genetic diversity indicators showed low diversity. In the NJ tree, the Congjiang and Liping populations were mixed together, and the other three groups formed a cluster. A cross-validation error box plot and pong cluster plot were constructed to show the K value results. When K = 1, the cross-validation error rate was the lowest. Principal component analysis showed that the Duyun population formed a group separate from the group comprising the other four populations. The genetic differentiation index and genetic distances between the Duyun population and the remaining four populations were greater than 0.05, indicating population differentiation. The genetic diversity of the five populations of golden-backed carp in Guizhou Province was low, the genetic differentiation of the Duyun population was the most significant, and the Duyun population was separate from the other four groups.

RESUMEN

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria-Canary Islands, Faro-Portugal, and Ebro Delta-Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, "regulation of transcription" and "signalling" were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

RESUMEN

Knowledge on correlations between environmental factors and genome divergence between populations of marine species is crucial for sustainable management of fisheries and wild populations. The edible cockle (Cerastoderma edule) is a marine bivalve distributed along the Northeast Atlantic coast of Europe and is an important resource from both commercial and ecological perspectives. We performed a population genomics screening using 2b-RAD genotyping on 9309 SNPs localized in the cockle's genome on a sample of 536 specimens pertaining to 14 beds in the Northeast Atlantic Ocean to analyse the genetic structure with regard to environmental variables. Larval dispersal modelling considering species behaviour and interannual/interseasonal variation in ocean conditions was carried out as an essential background to which compare genetic information. Cockle populations in the Northeast Atlantic displayed low but significant geographical differentiation between populations (F ST = 0.0240; p < 0.001), albeit not across generations. We identified 742 and 36 outlier SNPs related to divergent and balancing selection in all the geographical scenarios inspected, and sea temperature and salinity were the main environmental correlates suggested. Highly significant linkage disequilibrium was detected at specific genomic regions against the very low values observed across the whole genome. Two main genetic groups were identified, northwards and southwards of French Brittany. Larval dispersal modelling suggested a barrier for larval dispersal linked to the Ushant front that could explain these two genetic clusters. Further genetic subdivision was observed using outlier loci and considering larval advection. The northern group was divided into the Irish/Celtic Seas and the English Channel/North Sea, while the southern group was divided into three subgroups. This information represents the baseline for the management of cockles, designing conservation strategies, founding broodstock for depleted beds and producing suitable seed for aquaculture production.

RESUMEN

BACKGROUND: Currently, Chinese laboratory macaques are widely used in biomedical research. Correspondingly, clarity regarding the genetic diversity of Chinese laboratory macaques is important for both vendors and users. METHODS: Genome sequences of 55 laboratory macaques (40 cynomolgus macaques and 15 rhesus macaques) bred in South China were analyzed using 2b-RAD simplified genome sequencing. RESULTS: A total of 115,681 single-nucleotide polymorphisms (SNPs) were found that were distributed in 21 chromosomes and an unplaced scaffold. Genetic diversity indices varied across populations and exhibited low values. The results of principal coordinate analysis (PCA) were consistent with those of the arithmetic mean (UPGMA) clustered tree and supported the structure analysis, demonstrating that the genetic differentiation in rhesus macaques was higher than that in cynomolgus macaques. Introgressive hybridization with the Chinese rhesus macaque was supported in more than 80% (32/40) of cynomolgus macaques. CONCLUSIONS: Chinese laboratory macaques had relatively low genetic diversity at the genomic level, and genetic differentiation in Chinese rhesus macaques was higher than in cynomolgus macaques. The genome of cynomolgus macaques has been shaped by introgression after hybridization with the Chinese rhesus macaques.

Asunto(s)

RESUMEN

Sea urchin (Strongylocentrotus intermedius) is an economically important mariculture species in Asia, and its gonads are the only edible part. The efficiency of genetic breeding in sea urchins is hampered due to the inability to distinguish gender by appearance. In this study, we first identified a sex-associated single nucleotide polymorphism (SNP) by combining type IIB endonuclease restriction site-associated DNA sequencing (2b-RAD-seq) and genome survey. Importantly, this SNP is located within spata4, a gene specifically expressed in male. Knocking down of spata4 by RNA interference (RNAi) in male individuals led to the downregulation of other conserved testis differentiation-related genes and germ cell marker genes. We also revealed that sex ratio in this validated culture population of S. intermedius is not 1:1. Moreover, after a 58-day feeding experiment with estradiol, the expression levels of several conserved genes that are related to testis differentiation, ovary differentiation, and estrogen metabolism were dynamically changed. Taken together, our results will contribute toward improving breeding efficiency, developing sex-controlled breeding, and providing a solid base for understanding sex determination mechanisms in sea urchins.

Asunto(s)

RESUMEN

Sex-specific markers play an important role in revealing sex-determination mechanism. Sea urchin (Mesocentrotus nudus) is an economically important mariculture species in several Asian countries and its gonads are the sole edible parts for people. However, growth rate and immunocompetence differ by sex in this species, sex-specific markers have not been identified, and the sex-determination mechanism of sea urchin remains undetermined. In this study, type IIB endonuclease restriction-site associated DNA sequencing (2b-RAD-seq) and a genome survey of M. nudus were performed, and three female-specific markers and three female heterogametic single nucleotide polymorphism (SNP) loci were identified. We validated these sex-specific markers via PCR amplification in a large number of individuals, including wild and artificially bred populations. Several open reading frames (ORFs) were predicted, although there are no potential genes known for sex determination and sex differentiation within the scaffold in which the sex-specific markers are located. Importantly, the female-specific sequences and female heterozygous SNP loci indicate that a female heterogametic and male homogametic ZW/ZZ sex-determination system should exist in M. nudus. The results provide a solid basis for revealing the sex-determination mechanism of this species, and open up new possibilities for developing sex-control breeding in sea urchin.

RESUMEN

Because of their extraordinary flower and leaf morphology, passion flowers (Passifloraceae) have fascinated naturalists since their discovery. Within the large, diverse (600 species) genus Passiflora is an especially enigmatic and species-rich (120 spp.) subclade, Section Decaloba, which occurs in the Neotropics and has its center of diversity in Andean montane forests. A recent phylogenetic study of Passifloraceae showed that Section Decaloba was monophyletic, but was unable to resolve relationships within the clade, thus preventing inferences of evolutionary history and biogeography. The goal of this study was to elucidate the phylogeny and biogeography of Section Decaloba. We sampled 206 accessions representing 91 of the ~ 120 known species in section Decaloba and four outgroups, with samples derived predominantly from herbarium specimens. We generated DNA sequences using a high-throughput DNA sequencing technique called 2b-RAD, reconstructed the phylogeny, and conducted ancestral area reconstructions to infer the biogeographic history of the group. We recovered predominantly well-supported trees in which species were grouped into two main clades: 1) the Central American clade, within which the majority of nodes well supported and species were monophyletic and 2) the South American clade, a large clade that showed overall lower resolution and included several polyphyletic species and species complexes that need additional research. RASP analysis showed that section Decaloba originated in Central America around 10.4 Ma, and then dispersed to South America, the Greater Antilles, and the Bahamas. The South American clade diversified in the Northern Andes and then dispersed to the rest of South America, and Lesser Antilles. Results suggest that both long-distance dispersal and colonization of newly available habitats (i.e., in the Andes) likely promoted diversification of this clade. This study also illustrates how using herbarium specimens and a RAD-seq approach can produce phylogenies for broadly distributed, highly diverse, and poorly accessible groups of plants where field collections would be unfeasible.

Asunto(s)

RESUMEN

Despite that the ballast water management (BWM) convention has come into force to prevent the spread of harmful aquatic organisms, to date, very few bacteria can be identified through microbial culture method. In this study, we explored a reduced-representation sequencing of 2b-RAD approach to investigate the bacterial diversity in ballast water and sediments (BWS). Our results indicated a large amount of bacteria species (1496) detected in BWS up to now, including 13 pathogens that are seriously concerning in marine environment and aquaculture like the most harmful Vibrio harveyi and Aurantimonas coralicida. We showed that the ballast water had relative lower species, which was dominated by Proteobacteria. In contrast, the sediments had richer species, which was dominated by Bacteroidetes. Although BWS differed significantly in species composition, sediments shared most of the concerned pathogens with ballast water, highlighting the importance of sediment management. In conclusion, 2b-RAD sequencing shows promise in future BWM.

Asunto(s)

RESUMEN

Red sage (Salvia miltiorrhiza) is a widely used medicinal plant for treatment of cardiovascular and cerebrovascular diseases. Because of excessive excavation by huge market demand and habitat loss by human activities, the wild population resources of S. miltiorrhiza have reduced drastically in recent years. Meanwhile, population status of two closely related species S. bowleyana and S. paramiltiorrhiza were in a trend of decreasing due to their potential replacement of S. miltiorrhiza. Particularly, S. paramiltiorrhiza was threatened and endemic to a small region in eastern China. However, to date there has been no conservation genetic research reported for wild S. miltiorrhiza population and its endangered relatives. Assess the wild germplasm diversity for S. miltiorrhiza and its related species would provide fundamental genetic background for cultivation and molecular breeding of this medicinally important species. In the present study, we investigated the genetic diversity, population structure, and intra/inter-specific differentiation of S. miltiorrhiza and above two relatives using 2b-RAD genome-wide genotyping method. By investigating 81 individuals of S. miltiorrhiza, 55 individuals of S. bowleyana and 15 individuals of S. paramiltiorrhiza from 23 locations in China, we obtained 23,928 SNPs in total. A comparatively high genetic diversity was observed in S. miltiorrhiza (π = 0.0788, H e = 0.0783 ± 0.0007). The observed and expected heterozygosity in populations of these three species ranged from 0.0297 to 0.1481 and 0.0251 to 0.831, respectively. Two major lineage groups were detected in the examined S. miltiorrhiza populations. The results indicated that Dabie Mountain as a genetic diversity center of S. miltiorrhiza and possible complex inter-specific genetic exchange/hybridization occurred between S. miltiorrhiza and the two relatives. We suggest that strategic conservation and germplasm preservation should be considered not only for wild populations of S. miltiorrhiza, but also for its related S. bowleyana and S. paramiltiorrhiza.

RESUMEN

Little is known about the mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different habitat temperature. Here, the genome-wide genetic information of three O. oratoria populations were obtained by IIB restriction site-associated DNA (2b-RAD) sequencing and 2403 single-nucleotide polymorphisms (SNPs) were identified. Based on the 2403 SNPs, we found a remarkable genetic differentiation between the Yellow Sea and the East China Sea groups of O. oratoria. Furthermore, 63 SNPs are thought to be associated with different sea temperatures. Based on the 63 SNPs, it is hypothesised that the long-term temperature differences may contribute to the variation of genes associated with multiple biological functions, such as material metabolism, cytoskeleton, cellular processes, inflammatory response and hormonal regulation. This study provides new information for elucidating the molecular mechanisms underlying the relationship between genetic variation and the adaptation of Oratosquilla oratoria populations to different temperature.

Asunto(s)

RESUMEN

BACKGROUND: The irruption of Next-generation sequencing (NGS) and restriction site-associated DNA sequencing (RAD-seq) in the last decade has led to the identification of thousands of molecular markers and their genotyping for refined genomic screening. This approach has been especially useful for non-model organisms with limited genomic resources. Many building-loci pipelines have been developed to obtain robust single nucleotide polymorphism (SNPs) genotyping datasets using a de novo RAD-seq approach, i.e. without reference genomes. Here, the performances of two building-loci pipelines, STACKS 2 and Meyer's 2b-RAD v2.1 pipeline, were compared using a diverse set of aquatic species representing different genomic and/or population structure scenarios. Two bivalve species (Manila clam and common edible cockle) and three fish species (brown trout, silver catfish and small-spotted catshark) were studied. Four SNP panels were evaluated in each species to test both different building-loci pipelines and criteria for SNP selection. Furthermore, for Manila clam and brown trout, a reference genome approach was used as control. RESULTS: Despite different outcomes were observed between pipelines and species with the diverse SNP calling and filtering steps tested, no remarkable differences were found on genetic diversity and differentiation within species with the SNP panels obtained with a de novo approach. The main differences were found in brown trout between the de novo and reference genome approaches. Genotyped vs missing data mismatches were the main genotyping difference detected between the two building-loci pipelines or between the de novo and reference genome comparisons. CONCLUSIONS: Tested building-loci pipelines for selection of SNP panels seem to have low influence on population genetics inference across the diverse case-study scenarios here studied. However, preliminary trials with different bioinformatic pipelines are suggested to evaluate their influence on population parameters according with the specific goals of each study.

Asunto(s)

RESUMEN

Finding optimal sample sizes is critical for the accurate estimation of genetic diversity of large invasive populations. Based on previous studies, we hypothesized that a minimal sample size of 3-8 individuals is sufficient to dissect the population architecture of the harlequin lady beetle, Harmonia axyridis, a biological control agent and an invasive alien species. Here, equipped with a type IIB endonuclease restriction site-associated (2b-RAD) DNA sequencing approach, we identified 13,766 and 13,929 single nucleotide polymorphisms (SNPs), respectively, among native and invasive H. axyridis populations. With this information we simulated populations using a randomly selected 3000 SNPs and a subset of individuals. From this simulation we finally determined that six individuals is the minimum sample size required for the accurate estimation of intra- and inter-population genetic diversity within and across H. axyridis populations. Our findings provide an empirical advantage for population genomic studies of H. axyridis in particular and suggest useful tactics for similar studies on multicellular organisms in general.

RESUMEN

Possessing powerful adaptive capacity and a pleasant taste, spotted sea bass (Lateolabrax maculatus) has a broad natural distribution and is one of the most popular mariculture fish in China. However, the genetic improvement program for this fish is still in its infancy. Growth is the most economically important trait and is controlled by quantitative trait loci (QTL); thus, the identification of QTLs and genetic markers for growth-related traits is an essential step for the establishment of marker-assisted selection (MAS) breeding programs. In this study, we report the first high-density linkage map of spotted sea bass constructed by sequencing 333 F1 generation individuals in a full-sib family using 2b-RAD technology. A total of 6883 SNP markers were anchored onto 24 linkage groups, spanning 2189.96 cM with an average marker interval of 0.33 cM. Twenty-four growth-related QTLs, including 13 QTLs for body weight and 11 QTLs for body length, were successfully detected, with phenotypic variance explained (PVE) ranging from 5.1 to 8.6%. Thirty potential candidate growth-related genes surrounding the associated SNPs were involved in cell adhesion, cell proliferation, cytoskeleton reorganization, calcium channels, and neuromodulation. Notably, the fgfr4 gene was detected in the most significant QTL; this gene plays a pivotal role in myogenesis and bone growth. The results of this study may facilitate marker-assisted selection for breeding populations and establish the foundation for further genomic and genetic studies investigating spotted sea bass.

Asunto(s)

RESUMEN

Global change forecasts in ecosystems require knowledge of within-species diversity, particularly of dominant species within communities. We assessed site-level diversity and capacity for adaptation in Bouteloua gracilis, the dominant species in the Central US shortgrass steppe biome. We quantified genetic diversity from 17 sites across regional scales, north to south from New Mexico to South Dakota, and local scales in northern Colorado. We also quantified phenotype and plasticity within and among sites and determined the extent to which phenotypic diversity in B. gracilis was correlated with climate. Genome sequencing indicated pronounced population structure at the regional scale, and local differences indicated that gene flow and/or dispersal may also be limited. Within a common environment, we found evidence of genetic divergence in biomass-related phenotypes, plasticity, and phenotypic variance, indicating functional divergence and different adaptive potential. Phenotypes were differentiated according to climate, chiefly median Palmer Hydrological Drought Index and other aridity metrics. Our results indicate conclusive differences in genetic variation, phenotype, and plasticity in this species and suggest a mechanism explaining variation in shortgrass steppe community responses to global change. This analysis of B. gracilis intraspecific diversity across spatial scales will improve conservation and management of the shortgrass steppe ecosystem in the future.

Asunto(s)

RESUMEN

Analysis of population genetics provides insights into the evolutionary processes, among which the sample size choice is per se a crucial issue in the analysis. Genome-wide high-throughput techniques based on RADseq have been increasingly used in studies on the population genomics of invasive species. However, there is little information available regarding optimal sample sizes for analyzing population genomics of invasive species. In this study, we first use type IIB endonucleases restriction site-associated DNA (2b-RAD) to mine thousands of single nucleotide polymorphisms (SNPs) for native and introduced populations in Q1 clade (SPB and 17JN) and Q2 clade (ISQ and UAS0601) of the whitefly, Bemisia tabaci (Gennadius) MED (also known as B. tabaci biotype Q). Then, we used resampling techniques to create simulated populations with a random subset of individuals and 3,000 SNPs to determine how many individuals should be sampled for accurate estimates of intra- and interpopulation genetic diversity. We calculated the intrapopulation genetic diversity parameters (unbiased expected heterozygosity, observed heterozygosity, and the number of effect alleles) and pairwise genetic differentiation F ST; finally, an ad hoc statistic, ΔK, was used to determine the optimal value. Our results showed that a sample size greater than four individuals (n ≥ 4) has little impact on estimates of genetic diversity within whitefly populations; moreover, precise estimate of F ST can be easily achieved at a very small simple size (n = 3 or 4). Our results will provide in-depth understanding of the optimization of sampling schemes in population genomics of invasive species.