Development and evaluation of 200 novel SNP assays for population genetic studies of westslope cutthroat trout and genetic identification of related taxa

DNA sequence data were collected and screened for single nucleotide polymorphisms (SNPs) in westslope cutthroat trout (Oncorhynchus clarki lewisi) and also for substitutions that could be used to genetically discriminate rainbow trout (O. mykiss) and cutthroat trout, as well as several cutthroat trout subspecies. In total, 260 expressed sequence tag‐derived loci were sequenced and allelic discrimination genotyping assays developed from 217 of the variable sites. Another 50 putative SNPs in westslope cutthroat trout were identified by restriction‐site‐associated DNA sequencing, and seven of these were developed into assays. Twelve O. mykiss SNP assays that were variable within westslope cutthroat trout and 12 previously published SNP assays were also included in downstream testing. A total of 241 assays were tested on six westslope cutthroat trout populations (N = 32 per population), as well as collections of four other cutthroat trout subspecies and a population of rainbow trout. All assays were evaluated for reliability and deviation from Hardy–Weinberg and linkage equilibria. Poorly performing and duplicate assays were removed from the data set, and the remaining 200 assays were used in tests of population differentiation. The remaining markers easily distinguished the various subspecies tested, as evidenced by mean G_ST of 0.74. A smaller subset of the markers (N = 86; average G_ST = 0.40) was useful for distinguishing the six populations of westslope cutthroat trout. This study increases by an order of magnitude the number of genetic markers available for the study of westslope cutthroat trout and closely related taxa and includes many markers in genes (developed from ESTs).     

Fine-scale genetic population structure in a mobile marine mammal: inshore bottlenose dolphins in Moreton Bay, Australia

Highly mobile marine species in areas with no obvious geographic barriers are expected to show low levels of genetic differentiation. However, small‐scale variation in habitat may lead to resource polymorphisms and drive local differentiation by adaptive divergence. Using nuclear microsatellite genotyping at 20 loci, and mitochondrial control region sequencing, we investigated fine‐scale population structuring of inshore bottlenose dolphins (Tursiops aduncus) inhabiting a range of habitats in and around Moreton Bay, Australia. Bayesian structure analysis identified two genetic clusters within Moreton Bay, with evidence of admixture between them (F_ST = 0.05, P = 0.001). There was only weak isolation by distance but one cluster of dolphins was more likely to be found in shallow southern areas and the other in the deeper waters of the central northern bay. In further analysis removing admixed individuals, southern dolphins appeared genetically restricted with lower levels of variation (AR = 3.252, π = 0.003) and high mean relatedness (r = 0.239) between individuals. In contrast, northern dolphins were more diverse (AR = 4.850, π = 0.009) and were mixing with a group of dolphins outside the bay (microsatellite‐based STRUCTURE analysis), which appears to have historically been distinct from the bay dolphins (mtDNA Φ_ST = 0.272, P < 0.001). This study demonstrates the ability of genetic techniques to expose fine‐scale patterns of population structure and explore their origins and mechanisms. A complex variety of inter‐related factors including local habitat variation, differential resource use, social behaviour and learning, and anthropogenic disturbances are likely to have played a role in driving fine‐scale population structure among bottlenose dolphins in Moreton Bay.     

Sequence polymorphisms in porcine homologs of murine coat colour‐related genes

Herein, we report the variability among 57 porcine homologs of murine coat colour‐related genes. We identified single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within 44 expressed gene sequences by aligning eight pig complementary DNA (cDNA) samples. The sequence alignment revealed a total of 485 SNPs and 15 InDels. The polymorphisms were then validated by performing matrix‐assisted laser desorption/ionization time‐of‐flight mass spectrometry (MALDI‐TOF MS) with reference DNA samples obtained from 384 porcine individuals. Of the 384 individuals, three parents of the experimental F₂ family were included to detect polymorphisms between them for linkage mapping. We also genotyped previously reported polymorphisms of 12 genes, and one SNP each in three genes that were detected by performing a BLAST search of the Trace database. A total of 211 SNPs and three InDels were successfully genotyped from our porcine DNA panel. We detected SNPs in 33 of the 44 genes among the parents of an experimental F₂ family and then constructed a linkage map of the 33 genes for this family. The linkage assignment of each gene to the porcine chromosomes was consistent with the location of the BAC clone in the porcine genome and the corresponding gene sequence. We confirmed complete substitutions of EDNRB and MLPH in the Jinhua and Clawn miniature breeds, respectively. Furthermore, we identified polymorphic alleles exclusive to each pig group: 13 for Jinhua, two for Duroc, three for Meishan, four for the Japanese wild boar, one for the Clawn miniature pig and four for the Potbelly pig.     

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.     

Mitochondrial DNA phylogeography of Spodoptera exigua across a broad geographic area in China

The beet armyworm, Spodoptera exigua (Lepidoptera: Noctuidae), is an important agriculture pest in China that causes serious economic losses in some of the main crop‐producing areas. To monitor and manage this pest effectively, it is necessary to investigate its phylogeographic patterns in China. In this study, we used a partial sequence of mitochondrial DNA Cytb gene consisting of 724 bps to investigate the genetic diversity of the beet armyworm. A total of 765 individuals from 47 populations across the main distribution range of the species were collected, and 112 haplotypes were identified. Moderate‐to‐high levels of genetic diversity (Hd = 0.672 ± 0.017, Pi = 0.00268 ± 0.00021) for the total populations were obtained. Phylogenetic and median‐joining network analyses indicated there was no distinct geographic distribution pattern among haplotypes. Overall, the study also revealed significant differentiation among some populations (P < 0.05). The F_ST values of Shenyang population (SY2012–SY 2014), as well as Baoding (BD), Taian (TA), Lucheng (LC), Zhengzhou (ZZ) and Wuhan (WH), were significantly different from those of the populations in most other locations. Hierarchical AMOVA showed there was no significantly genetic structure between populations located in seven geographic regions and four main bioclimatic zones. Finally, unimodal mismatch distribution combined with negative Tajima's D (D = ?2.696, P < 0.001) and Fu's F_S (F_S = ?207.228, P > 0.05) indicated recent population expansion of S. exigua at large spatial scales in China.     

Isotopic discrimination and kinetic parameters of RubisCO from the marine bloom‐forming diatom,Skeletonema costatum

The cosmopolitan, bloom‐forming diatom, Skeletonema costatum, is a prominent primary producer in coastal oceans, fixing CO₂ with ribulose 1,5‐bisphosphate carboxylase/oxygenase (RubisCO) that is phylogenetically distinct from terrestrial plant RubisCO. RubisCOs are subdivided into groups based on sequence similarity of their large subunits (IA–ID, II, and III). ID is present in several major oceanic primary producers, including diatoms such as S. costatum, coccolithophores, and some dinoflagellates, and differs substantially in amino acid sequence from the well‐studied IB enzymes present in most cyanobacteria and in green algae and plants. Despite this sequence divergence, and differences in isotopic discrimination apparent in other RubisCO enzymes, stable carbon isotope compositions of diatoms and other marine phytoplankton are generally interpreted assuming enzymatic isotopic discrimination similar to spinach RubisCO (IB). To interpret phytoplankton δ¹³C values, S. costatum RubisCO was characterized via sequence analysis, and measurement of its K_CO2 and V_max, and degree of isotopic discrimination. The sequence of this enzyme placed it among other diatom ID RubisCOs. Michaelis‐Menten parameters were similar to other ID enzymes (K_CO2 = 48.9 ± 2.8 μm ; V_max = 165.1 ± 6.3 nmol min^?1 mg^?1). However, isotopic discrimination (ε = [¹²k/¹³k ? 1] × 1000) was low (18.5‰; 17.0–19.9, 95% CI) when compared to IA and IB RubisCOs (22–29‰), though not as low as ID from coccolithophore, Emiliania huxleyi (11.1‰). Variability in ε‐values among RubisCOs from primary producers is likely reflected in δ¹³C values of oceanic biomass. Currently, δ¹³C variability is ascribed to physical or chemical factors (e.g. illumination, nutrient availability) and physiological responses to these factors (e.g. carbon‐concentrating mechanisms). Estimating the importance of these factors from δ¹³C measurements requires an accurate ε‐value, and a mass‐balance model using the ε‐value for S. costatum RubisCO is presented. Clearly, appropriate ε‐values must be included in interpreting δ¹³C values of environmental samples.     

Epigenetic variation reflects dynamic habitat conditions in a rare floodplain herb

Variation of DNA methylation is thought to play an important role for rapid adjustments of plant populations to dynamic environmental conditions, thus compensating for the relatively slow response time of genetic adaptations. However, genetic and epigenetic variation of wild plant populations has not yet been directly compared in fast changing environments. Here, we surveyed populations of Viola elatior from two adjacent habitat types along a successional gradient characterized by strong differences in light availability. Using amplified fragment length polymorphisms (AFLP) and methylation‐sensitive amplification polymorphisms (MSAP) analyses, we found relatively low levels of genetic (H′_gen = 0.19) and epigenetic (H′_epi = 0.23) diversity and high genetic (?_ST = 0.72) and epigenetic (?_ST = 0.51) population differentiation. Diversity and differentiation were significantly correlated, suggesting that epigenetic variation partly depends on the same driving forces as genetic variation. Correlation‐based genome scans detected comparable levels of genetic (17.0%) and epigenetic (14.2%) outlier markers associated with site specific light availability. However, as revealed by separate differentiation‐based genome scans for AFLP, only few genetic markers seemed to be actually under positive selection (0–4.5%). Moreover, principal coordinates analyses and Mantel tests showed that overall epigenetic variation was more closely related to habitat conditions, indicating that environmentally induced methylation changes may lead to convergence of populations experiencing similar habitat conditions and thus may play a major role for the transient and/or heritable adjustment to changing environments. Additionally, using a new MSAP‐scoring approach, we found that mainly the unmethylated (?_ST = 0.60) and CG‐methylated states (?_ST = 0.46) of epiloci contributed to population differentiation and putative habitat‐related adaptation, whereas CHG‐hemimethylated states (?_ST = 0.21) only played a marginal role.     

Immunochip analysis identifies association of the RAD50/IL13 region with human longevity

Human longevity is characterized by a remarkable lack of confirmed genetic associations. Here, we report on the identification of a novel locus for longevity in the RAD50/IL13 region on chromosome 5q31.1 using a combined European sample of 3208 long‐lived individuals (LLI) and 8919 younger controls. First, we performed a large‐scale association study on 1458 German LLI (mean age 99.0 years) and 6368 controls (mean age 57.2 years) by targeting known immune‐associated loci covered by the Immunochip. The analysis of 142 136 autosomal single nucleotide polymorphisms (SNPs) revealed an Immunochip‐wide significant signal (P_Immunochip = 7.01 × 10^–9) for the SNP rs2075650 in the TOMM40/APOE region, which has been previously described in the context of human longevity. To identify novel susceptibility loci, we selected 15 markers with P_Immunochip < 5 × 10^–4 for replication in two samples from France (1257 LLI, mean age 102.4 years; 1811 controls, mean age 49.1 years) and Denmark (493 LLI, mean age 96.2 years; 740 controls, mean age 63.1 years). The association at SNP rs2706372 replicated in the French study collection and showed a similar trend in the Danish participants and was also significant in a meta‐analysis of the combined French and Danish data after adjusting for multiple testing. In a meta‐analysis of all three samples, rs2706372 reached a P‐value of P_{Immunochip+Repl} = 5.42 × 10^?7 (OR = 1.20; 95% CI = 1.12–1.28). SNP rs2706372 is located in the extended RAD50/IL13 region. RAD50 seems a plausible longevity candidate due to its involvement in DNA repair and inflammation. Further studies are needed to identify the functional variant(s) that predispose(s) to a long and healthy life.     

SNP discovery in nonmodel organisms: strand bias and base‐substitution errors reduce conversion rates

Single nucleotide polymorphisms (SNPs) have become the marker of choice for genetic studies in organisms of conservation, commercial or biological interest. Most SNP discovery projects in nonmodel organisms apply a strategy for identifying putative SNPs based on filtering rules that account for random sequencing errors. Here, we analyse data used to develop 4723 novel SNPs for the commercially important deep‐sea fish, orange roughy (Hoplostethus atlanticus), to assess the impact of not accounting for systematic sequencing errors when filtering identified polymorphisms when discovering SNPs. We used SAMtools to identify polymorphisms in a velvet assembly of genomic DNA sequence data from seven individuals. The resulting set of polymorphisms were filtered to minimize ‘bycatch’—polymorphisms caused by sequencing or assembly error. An Illumina Infinium SNP chip was used to genotype a final set of 7714 polymorphisms across 1734 individuals. Five predictors were examined for their effect on the probability of obtaining an assayable SNP: depth of coverage, number of reads that support a variant, polymorphism type (e.g. A/C), strand‐bias and Illumina SNP probe design score. Our results indicate that filtering out systematic sequencing errors could substantially improve the efficiency of SNP discovery. We show that BLASTX can be used as an efficient tool to identify single‐copy genomic regions in the absence of a reference genome. The results have implications for research aiming to identify assayable SNPs and build SNP genotyping assays for nonmodel organisms.     

A new single nucleotide polymorphism in the rabbit (Oryctolagus cuniculus) myostatin (MSTN) gene is associated with carcass composition traits

This study aimed at the identification of genetic variations in the myostatin (MSTN) gene and testing their effects on carcass quality traits. We comparatively sequenced Giant Grey (GG) and New Zealand White (NZW) rabbits that were founders of a cross‐bred population. Alignment of our sequence data with the GenBank sequence of the rabbit MSTN gene (Ensembl Gene ID ENSOCUG00000012663) identified three single nucleotide polymorphisms (SNPs). The two novel SNPs (c.?125T>C, c.373+234G>A) and one known SNP (c.747+34C>T) were subsequently analysed for linkage with carcass composition traits in 363 F₂ animals of the cross GG × NZW. Significant linkage was found between c.373+234G>A and nine carcass composition traits (P < 0.05). No significant effects were found for c.?125T>C and c.747+34C>T. Because the linked SNP is located in intron 1 and no genetic variation was found in the coding region, further investigations are necessary to understand the functional effect of the c.373+234G>A variant on the variability of the traits.     

Isolation and development of microsatellite markers for the Japanese dormouse, Glirulus japonicus

Eight microsatellite markers were developed for the Japanese dormouse (Glirulus japonicus), a natural monument and near‐threatened species in Japan. The markers amplify in individuals from all of the mitochondrial lineages detected in a previous study. Numerous polymorphisms were detected in specimens from a local population in central Honshu (11–21 alleles per locus; n = 31) and from the entire distribution range of the species (19–41 alleles per locus; n = 152). These microsatellites will be useful in conservation genetic studies of G. japonicus.     

Age determination and growth estimates of the white‐spotted conger eel,Conger myriaster (Brevoort, 1856) in marine waters of South Korea

The age and growth of conger eel, Conger myriaster, were investigated by measuring transversely sectioned sagittal otoliths samples from 635 individuals. Sample ages ranged from 1 to 13 years in the female data. Parameters of the von Bertalanffy growth function were estimated using nonlinear regression from back‐calculation, mean length of samples at age relationships, and otolith weight‐at‐age relationships. Best‐fitting value of the three methods was the otolith weight‐at‐age relationship (r² = .87). Parameters of otolith weight‐at‐age were estimated as L_∞ = 143.76 cm, K = 0.081, and t₀ = ?1.285. Maximum oocyte diameter (MOD) ranged from 50 to 430 μm. Reproductive traits of ovaries showed a positive relationship between GSI and MOD (r² = .8515). It is suggested that oogenesis begins to develop from 4 years of age and at lengths of about 45 cm TL. In conclusion, these data provide reliable fundamental data for the fish stock management of Conger myriaster in South Korea.     

Age,growth and mortality of Lutjanus alexandrei in estuarine and coastal waters of the tropical south‐western Atlantic

Otolith‐based methods were used to determine life history traits of the endemic Brazilian snapper (Lutjanus alexandrei) in estuarine and coastal environments in the south‐western Atlantic. Fishes were caught as juveniles inside mangrove‐bordered estuaries by traditional corral fisheries whereas adults were captured at sea using motorboats with trap and gill nets. Fish were sampled during landings and 331 otolith pairs were extracted from L. alexandrei. Inshore mangroves comprised individuals of 0–4 years (mean: 2 years), while individuals in deeper reef environments were older (range: 3–22; mean: 8 years), indicating an ontogenetic shift at approximately age 3 or 4. Edge analysis was used to validate the annual deposition in the otoliths, suggesting that opaque growth rings were formed between April and September. Age‐at‐length data were used to predict L. alexandrei growth rates using the von Bertalanffy growth model from where the parameters were calculated: L_∞ = 31 cm, k = 0.24, t₀ = ?1.26, r² = 0.97. Mortality rates were estimated for coastal habitats, with Z = 0.22 and S = 0.78 year^?1, based on ages 7–17. Additionally, evidence of ontogenetic migration is provided by age and size structure.     

The effects of subinhibitory concentrations of costus oil on virulence factor production in Staphylococcus aureus

Aim: To determine the antimicrobial activity of costus (Saussurea lappa) oil against Staphylococcus aureus, and to evaluate the influence of subinhibitory concentrations of costus oil on virulence‐related exoprotein production in staph. aureus. Methods and Results: Minimal inhibitory concentrations (MICs) were determined using a broth microdilution method, and the MICs of costus oil against 32 Staph. aureus strains ranged from 0.15 to 0.6 μl ml^?1. The MIC₅₀ and MIC₉₀ were 0.3 and 0.6 μl ml^?1, respectively. Western blot, haemolytic, tumour necrosis factor (TNF) release and real‐time RT‐PCR assays were performed to evaluate the effects of subinhibitory concentrations of costus oil on virulence‐associated exoprotein production in Staph. aureus. The data presented here show that costus oil dose dependently decreased the production of α‐toxin, toxic shock syndrome toxin 1 (TSST‐1) and enterotoxins A and B in both methicillin‐sensitive Staph. aureus (MSSA) and methicillin‐resistant Staph. aureus (MRSA). Conclusion: Costus oil has potent antimicrobial activity against Staph. aureus, and the production of α‐toxin, TSST‐1 and enterotoxins A and B in Staph. aureus was decreased by costus oil. Significance and Impact of the Study: The data suggest that costus oil may deserve further investigation for its potential therapeutic value in treating Staph. aureus infections. Furthermore, costus oil could be rationally applied in food products as a novel food preservative both to inhibit the growth of Staph. aureus and to repress the production of exotoxins, particularly staphylococcal enterotoxins.     

Isolation and characterization of microsatellite loci for the Potentilla core group (Rosaceae) using 454 sequencing

Microsatellites are valuable markers for the analysis of genetic diversity, linkage mapping or genotyping. The limited availability of microsatellites for the genus Potentilla (Rosaceae) stipulated the isolation of markers from a representative (Potentilla pusilla Host) of the Potentilla core group that constitutes the most species‐rich evolutionary lineage within the genus. Thousand four hundred and seventy‐six simple sequence repeat (SSR) containing candidate sequences were isolated from a single‐type line using 454 sequencing. Seventy‐four functional microsatellite markers were developed from 200 sequences selected for suitable priming sites flanking microsatellite repeats referring to a 37% primer‐to‐marker conversion ratio. Seventy‐two markers were polymorphic. These numbers confirm the increased efficiency of pyrosequencing over traditional isolation techniques in the development of microsatellites. Amplification primer sequences and the sequences of corresponding target fragments are provided for all functional markers, and molecular polymorphisms estimated for four accessions of P. pusilla and among seven core group species represented by 14 individuals are reported. Cross‐species transferability ranged between 86.4% and 97.3% among the studied taxa, and 57, 11 and six of the selected primer pairs amplified fragments of expected size and number in seven, six and five of the species, respectively. Reproducibility of the molecular phenotypes was 97.0%, which was inferred using a replicate sample of P. pusilla.     

Dynamics of cell cycle phase perturbations by trabectedin (ET-743) in nucleotide excision repair (NER)-deficient and NER-proficient cells, unravelled by a novel mathematical simulation approach

Abstract. Objectives: Trabectedin (ET‐743, Yondelis^®) is a natural marine product, with antitumour activity, currently in phase II/III clinical trials. Previous studies have shown that cells hypersensitive to ultraviolet (UV)‐rays because of nucleotide excision repair (NER) deficiency, were resistant to trabectedin. The purpose of this study was to investigate whether this resistance was associated with different drug‐induced cell cycle perturbations. Materials and Methods: An isogenic NER‐proficient cellular system (CHO‐AA8) and a NER‐deficient one (CHO‐UV‐96), lacking functional ERCC‐1, were studied. Flow cytometric assays showed progressive accumulation of cells in G₂ + M phase in NER‐proficient but not in NER‐deficient cells. Applying a computer simulation method, we realized that the dynamics of the cell cycle perturbations in all phases were complex. Results: Cells exposed to trabectedin during G₁ and G₂ + M first experienced a G₁ block, while those exposed in S phase were delayed in S and G₂ + M phases but eventually divided. In the presence of functional NER, exit from the G₁ block was faster; then, cells progressed slowly through S phase and were subsequently blocked in G₂ + M phase. This G₂ + M processing of trabectedin‐induced damage in NER‐proficient cells was unable to restore cell cycling, suggesting a difficulty in repairing the damage. Conclusions: This might be due either to important damage left unrepaired by previous G₁ repair, or that NER activity itself caused DNA damage, or both. We speculate that in UV‐96 cells repair mechanisms other than NER are activated both in G₁ and G₂ + M phases.     

A 66‐bp deletion in growth hormone releasing hormone gene 5′‐flanking region with largemouth bass recessive embryonic lethal

Growth hormone releasing hormone (GHRH) regulates the secretion of growth hormone (GH) in the pituitary gland. A 66‐bp deletion (c.‐923_‐858del) was detected in the 5′‐flanking sequence of the largemouth bass (Micropterus salmoides) GHRH gene. In two cultured random populations of adult individuals (A: n = 170 and B: n = 150), the genotype ratios of +/+:+/? were 2.5:1 and 2.8:1 respectively. Only one ?/? fish was detected. A Largemouth bass family was constructed with two heterozygous individuals (+/?) as parents. The genotype ratio of +/+:+/?:?/? in the filial generation embryos was 1:1.6:0.1 at the neurula and 1:2:0 at hatched larvae stages. This indicated that the 66‐bp deletion was a recessive lethal site and that homozygous individuals (?/?) died off in embryonic development. The growth traits (body weight, body length and body depth) were measured, and the GHRH mRNA expression levels in brain tissue were detected using real‐time PCR. The effects of genotype (+/?) on growth traits and GHRH mRNA expression were not significant. Although the cause of death was not clear, the results hint that the 66‐bp deletion site in GHRH 5′‐flanking sequence significantly affects the livability in largemouth bass embryonic development.     

Development of EST‐derived microsatellite markers for Arabidopsis lyrata subspecies petraea (L.)

A set of expressed sequence tag–simple sequence repeat (EST‐SSR) loci has been developed for Arabidopsis lyrata ssp. petraea. From 768 root cDNA clones, 126 microsatellites, including di‐, tri‐, tetra‐ and pentanucleotide repeat motifs were identified and primers were designed to 24 EST‐SSRs. Eleven loci were subsequently screened on 150 individuals sampled from five natural populations, which revealed three to nine alleles per locus (mean 5.36) and expected heterozygosity (H_E) estimates ranging from 0.046 to 0.698. Significant deviations from random mating were observed at 10 EST‐SSR loci, likely due to inbreeding (global F_IS = 0.151) and population structure (global F_ST = 0.246).     

Genetic variation,phylogeographic structure of Spodoptera exigua in the welsh onion‐producing areas of North China

The beet armyworm, Spodoptera exigua (Lepidiptera: Noctuidae), is an important agricultural pest in China, causing serious economic losses in the main welsh onion‐producing areas of North China. In order to effectively monitor and manage this pest, it is necessary to investigate its genetic variation and intraspecific phylogeographic structure in North China. In this study, we used a DNA fragment of cytochrome oxidase subunit I (COI) consisting of 518 bps to investigate its genetic diversity. A total of 500 individuals of this species were collected from 20 geographic locations, and 8 haplotypes were identified. Low levels of genetic diversity (h = 0.256 ± 0.025, π = 0.00058 ± 0.00008) for the total populations in these study regions were obtained. Phylogenetic and median‐joining network analyses indicated that there was no distinct geographic distribution pattern among haplotypes. The results also indicated that the mean pairwise sequence difference between the populations was 0.1%, ranged from zero to 0.2%. The genetic differentiation (F_st) between the populations varied from ?0.198 to 0.656, and there was a significant differentiation between the population of Kunshan (KUNS), Shenyang (SY), Lucheng (LUC), Binzhou (BZ) and the remaining populations, respectively. The analysis of molecular variance (amova ) showed that the percentage of variation within a population (69.05%) was greater than that between the populations (30.95%), and the level of population differentiation (F_ST = 0.310, P < 0.001) was highly significant after 1000 random permutations. Finally, unimodal mismatch distribution combined with negative, Tajima's D (D = ?1.956, P < 0.05) and Fu's Fs (Fs = ?5.875, P > 0.05) indicated that S. exigua experienced recent population expansion in North China.