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.     

不同地区甜菜夜蛾种群的遗传多样性分析

甜菜夜蛾Spodoptera exigua (Hübner)是一种重要的农业害虫,曾给我国农牧业生产造成过重大的经济损失。为了阐明甜菜夜蛾的种群动态规律,改善和提高其预测、防治水平,我们应用扩增片段长度多态性(AFLP)技术对我国7个甜菜夜蛾种群的42头雄蛾的遗传多样性进行了研究分析。结果表明,不同种群个体间遗传相似性指数分布在0.143～0.824之间,同一种群内不同个体间的遗传相似性指数分布在0.250～0.786之间,同一种群内个体间的遗传相似性指数有许多小于不同种群个体间的。江西种群的多态性条带比例最低,为80.7%; 山东种群的多态性条带比例最高,达88.6%。总体而言,北方种群内的遗传多样性高于南方种群。聚类分析结果表明,同一种群的个体并不总能聚在一起,即种群间不存在明显的遗传分化。这些结果为明确我国甜菜夜蛾的迁飞规律及各发生危害区的虫源关系提供了进一步的实验依据。     

Response of multiple generations of beet armyworm,Spodoptera exigua (Hübner), feeding on transgenic Bt cotton

Development, reproduction and food utilization of three successive generations of beet armyworm, Spodoptera exigua (Hübner), fed on transgenic and non‐transgenic Bt cotton were examined. Significantly longer larval life‐span and lower pupal weight were observed in three successive generations of S. exigua fed on transgenic Bt cotton compared with non‐transgenic Bt cotton. Significantly higher survival rate and adult fecundity of S. exigua were found in three successive generations of S. exigua fed on transgenic Bt cotton compared with non‐transgenic Bt cotton. The survival rate and adult fecundity of S. exigua were occurred significant increase in the third generation compared with the first generation after feeding on transgenic Bt cotton. Significantly lower consumption, frass and relative growth rate (RGR) were observed in three successive generations of S. exigua fed on transgenic Bt cotton compared with non‐transgenic Bt cotton. Cotton variety significantly affected all indices of larval consumption and utilization in three successive generations of S. exigua, except for efficiency of conversion of ingested food. However, beet armyworm generation only significantly affected RGR of S. exigua. The results of this study indicated food quality on the diet‐utilization efficiency of S. exigua was different along with beet armyworm generation. Measuring multigenerational development and food utilization of S. exigua at individual and population level in response to Bacillus thuringiensis (Bt) can provide a more meaningful evaluation of long‐term population dynamics than experiments on a single generation. It is imperative to develop an appropriate multigenerational pest management tactic to monitor the field population dynamics of non‐target pests (e.g., beet armyworm) in agricultural Bt cotton ecosystem.     

Dynamic micro-geographic and temporal genetic diversity in vertebrates: the case of lake-spawning populations of brown trout (Salmo trutta)

Conservation of species should be based on knowledge of effective population sizes and understanding of how breeding tactics and selection of recruitment habitats lead to genetic structuring. In the stream‐spawning and genetically diverse brown trout, spawning and rearing areas may be restricted source habitats. Spatio–temporal genetic variability patterns were studied in brown trout occupying three lakes characterized by restricted stream habitat but high recruitment levels. This suggested non‐typical lake‐spawning, potentially representing additional spatio–temporal genetic variation in continuous habitats. Three years of sampling documented presence of young‐of‐the‐year cohorts in littoral lake areas with groundwater inflow, confirming lake‐spawning trout in all three lakes. Nine microsatellite markers assayed across 901 young‐of‐the‐year individuals indicated overall substantial genetic differentiation in space and time. Nested gene diversity analyses revealed highly significant (≤P = 0.002) differentiation on all hierarchical levels, represented by regional lakes (F_LT = 0.281), stream vs. lake habitat within regional lakes (F_HL = 0.045), sample site within habitats (F_SH = 0.010), and cohorts within sample sites (F_CS = 0.016). Genetic structuring was, however, different among lakes. It was more pronounced in a natural lake, which exhibited temporally stable structuring both between two lake‐spawning populations and between lake‐ and stream spawners. Hence, it is demonstrated that lake‐spawning brown trout form genetically distinct populations and may significantly contribute to genetic diversity. In another lake, differentiation was substantial between stream‐ and lake‐spawning populations but not within habitat. In the third lake, there was less apparent spatial or temporal genetic structuring. Calculation of effective population sizes suggested small spawning populations in general, both within streams and lakes, and indicates that the presence of lake‐spawning populations tended to reduce genetic drift in the total (meta‐) population of the lake.     

Replicate divergence between and within sounds in a marine fish: the copper rockfish (Sebastes caurinus)

Understanding the factors that influence larval dispersal and connectivity among marine populations is critical to the conservation and sustainable management of marine resources. We assessed genetic subdivision among ten populations of copper rockfish (Sebastes caurinus) representing paired samples of outer coast and the heads of inlets in five replicate sounds on the west coast of Vancouver Island, British Columbia, using 17 microsatellite DNA loci. Overall, subdivision (F_ST) was low (F_ST = 0.031, P < 0.001), but consistently higher between paired coast and head of inlet sites (mean F_ST = 0.047, P < 0.001) compared to among the five coast sites (mean F_ST = ?0.001, P > 0.5) or among the five head of inlet sites (mean F_ST = 0.026, P < 0.001). Heterozygosity, allelic richness and estimates of effective population size were also lower in head of inlet sites than in coast sites. Bayesian analysis identified two genetic groups across all samples, a single genetic group among only coast samples, two genetic groups among head of inlet samples and two genetic groups within each sound analysed separately. Head of inlet copper rockfish tended to be shorter with lower condition factors and grew more slowly than coast sites fish. Reduced physical connectivity and selection against immigrants in contrasting outer coast–head of inlet environments likely contribute to the evolution of population structure of copper rockfish. Based on genetic connectivity, coast sites appear to be better served by existing marine protected areas than are head of inlet sites.     

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep‐sea shark (Centroscymnus coelolepis)

Knowledge of the mechanisms limiting connectivity and gene flow in deep‐sea ecosystems is scarce, especially for deep‐sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep‐sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497‐bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ_ST = ?0.0091, F_ST = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ_ST = 0.0744, P < 0.01; F_ST = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ_ST = 0.3808, F_ST = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.     

Genetically distinct populations of northern shrimp,Pandalus borealis,in the North Atlantic: adaptation to different temperatures as an isolation factor

The large‐scale population genetic structure of northern shrimp, Pandalus borealis, was investigated over the species’ range in the North Atlantic, identifying multiple genetically distinct groups. Genetic divergence among sample localities varied among 10 microsatellite loci (range: F_ST = ?0.0002 to 0.0475) with a highly significant average (F_ST = 0.0149; P < 0.0001). In contrast, little or no genetic differences were observed among temporal replicates from the same localities (F_ST = 0.0004; P = 0.33). Spatial genetic patterns were compared to geographic distances, patterns of larval drift obtained through oceanographic modelling, and temperature differences, within a multiple linear regression framework. The best‐fit model included all three factors and explained approximately 29% of all spatial genetic divergence. However, geographic distance and larval drift alone had only minor effects (2.5–4.7%) on large‐scale genetic differentiation patterns, whereas bottom temperature differences explained most (26%). Larval drift was found to promote genetic homogeneity in parts of the study area with strong currents, but appeared ineffective across large temperature gradients. These findings highlight the breakdown of gene flow in a species with a long pelagic larval phase (up to 3 months) and indicate a role for local adaptation to temperature conditions in promoting evolutionary diversification and speciation in the marine environment.     

Rise and fall of a wolf population: genetic diversity and structure during recovery,rapid expansion and drastic decline

The grey wolves (Canis lupus) of Finland have had a varied history, with a period of rapid population expansion after the mid‐1990s followed by a decline with a current census size of about 140 wolves. Here, we investigate the impact of unstable population size and connectivity on genetic diversity and structure in a long‐term genetic study of 298 Finnish wolves born in 1995–2009 and genotyped for 17 microsatellite loci. During the initial recovery and prior to population expansion, genetic diversity was high (1995–1997: LD‐N_e = 67.2; H_o = 0.749; H_e = 0.709) despite a small census size and low number of breeders (N_c < 100; N_b < 10) likely reflecting the status of the Russian source population. Surprisingly, observed heterozygosity decreased significantly during the study period (t = ?2.643, P = 0.021) despite population expansion, likely a result of an increase in inbreeding (F_IS = 0.108 in 2007–2009) owing to a low degree of connectivity with adjacent Russian wolf population (m = 0.016–0.090; F_ST = 0.086, P < 0.001) and population crash after 2006. However, population growth had a temporary positive impact on N_e and number of family lines. This study shows that even strong population growth alone might not be adequate to retain genetic diversity, especially when accompanied with low amount of subsequent gene flow and population decline.     

Geographic variation in sex‐chromosome differentiation in the common frog (Rana temporaria)

In sharp contrast with birds and mammals, sex‐determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex‐determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500‐km transect across Sweden. Genetic differentiation between sexes is strongest (F_ST = 0.152) in a northern‐boreal population, where male‐specific alleles and heterozygote excesses (F_IS = ?0.418 in males, +0.025 in females) testify to a male‐heterogametic system and lack of X‐Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (F_ST = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male‐specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X‐Y recombination, co‐option of an alternative sex‐chromosome pair, or a mixed sex‐determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the ‘sexual races’ described in common frogs in the 1930s.     

Selection of nuclear polyhedrosis viruses as biological control agents ofSpodoptera exigua [Lep.: Noctuidae]

The virulence of 5 nuclear polyhedrosis viruses infectious for larvae of beet armyworm,Spodoptera exigua, was studied and their potential as biological control agents of this accidentally introduced pest in Dutch greenhouse crops is discussed. Three of the virus isolates were collected from deceased beet armyworm larvae found in Dutch greenhouses. Based on restriction endonuclease patterns of their DNA they appeared to be closely related toMamestra brassicae nuclear polyhedrosis virus (MbMNPV) and therefore were named MbMNPV-NL80, MbMNPV-NL82 and MbMNPV-NL83. These isolates were not related toAutographa californica MNPV (AcMNPV) or toSpodoptera exigua MNPV (SeMNPV), both originating from the USA. Comparison of the oiological activity of these 5 isolates showed that the SeMNPV was more virulent against beet armyworm than the other isolates. There was no significant difference in virulence between MbMNPV-NL80, NL82, NL83 and AcMNPV forS. exigua. The LD-50 values of the 5 isolates for 2nd instar larvae were 3, 26, 14, 17 and 18 polyhedra, respectively. Despite compensating qualities of the other MNPVs, such as a broader host range and potential production in alternate hosts or cell-lines, SeMNPV is considered to be the most suitable candidate as biological control agent of beet armyworm.      

Microsatellite‐based analysis of genetic structure and gene flow of Mythimna separata (Walker) (Lepidoptera: Noctuidae) in China

The oriental armyworm, Mythimna separata, is a serious agricultural pest in China. Seasonal and roundtrip migration has recently led to sudden, localized outbreaks and crop losses. To evaluate genetic differentiation between populations in eastern and western China and elucidate gene flow, the genetic structure of 20 natural populations from nine provinces was examined using seven microsatellite markers. The results indicated high genetic diversity. However, little to moderate (0 < F_ST < 0.15) genetic differentiation was detected, and there was no correlation between genetic distance and geographical distance. Bayesian clustering analysis identified three groups whereas discriminant analysis of principal components identified ten clusters that were considered as two clear‐cut clusters and one admixed group. Gene flow occurred frequently in most population pairs, and an asymmetrical migration rate was detected in several pairwise population comparisons. The bottleneck test showed that few populations had experienced recent bottlenecks. Correspondingly, large‐scale and long‐distance migration of M. separata has caused low genetic differentiation and frequent gene exchange. Our findings are important for studying genetic evolution and help to improve predictions of M. separata outbreaks in China.     

Evidence for regional diversity and host adaptation in Iranian populations of the Russian wheat aphid

Effective pest management is greatly facilitated by knowledge of the genetic structure and host adaptation of the pest species in question. The Russian wheat aphid (RWA), Diuraphis noxia (Mordvilko) (Homoptera: Aphididae: Macrosiphini), is an important economic pest in many cereal‐growing areas of the world, and in this study we investigated these aspects of its populations, using microsatellite markers and host plant response assays. Diuraphis noxia was sampled from 38 locations in Iran and genotyped at four polymorphic microsatellite loci that had been isolated from various Sitobion species. We identified 50 multilocus genotypes in 376 individuals. The overall observed heterozygosity was 0.134. F‐statistics showed a regional partitioning in D. noxia populations with overall F_ST = 0.231. In addition, there was a significant correlation between genetic and geographic distances. In order to test for the ecological consequences of genetic variability in D. noxia, biotypic variation amongst the isolates collected from wheat (Triticum aestivum L.) and barley (Hordeum vulgare L.) was evaluated on a number of resistant and susceptible wheat varieties. The plant variables we measured were damage rating (based on leaf chlorosis, leaf rolling, wilting, and death of the host plant), host plant dry weight, and root length. Damage rating was the best criterion for detecting biotypic variation in D. noxia. Discriminant analysis correctly classified the isolates in respective groups in 80–91.8% of the cases. The barley isolate showed no differences in performance on resistant and susceptible wheat, indicating a lack of gene‐by‐gene relationship with wheat plants. In contrast, wheat isolates differentially damaged the resistant and susceptible plants and showed moderate to severe virulence.     

Genomic analysis of the differential response to experimental infection with porcine circovirus 2b

Porcine circovirus type 2 (PCV2) is the etiological agent of a group of associated diseases (PCVAD) that affect production efficiency and can lead to mortality. Using different crossbred lines of pigs, we analyzed host genetic variation of viral load, immune response and weight change following experimental infection with a PCV2b strain (n = 386). Pigs expressed variation in the magnitude and initiation of viremia and immune response recorded weekly until 28 days post‐infection. A higher viral load was correlated with weight gain (r = ?0.26, P < 0.0001) and presence of PCV2‐specific antibodies (IgM, r = 0.26–0.34, P < 0.0001; IgG, r = 0.17–0.20, P < 0.01). In genome‐wide association analyses of the responses at different time points, the proportions of phenotypic variation explained by combined effects of 56 433 SNPs were 34.8–59.4% for viremia, 10.1–59.5% for antibody response and 5.6–14.9% for weight change. Relationships between genomic prediction of overall viral load and weight gain during the first weeks of challenge were negative (?0.21 and ?0.24 respectively, P < 0.0001). Individuals that carried more favorable alleles across three SNPs on SSC9 (0.60 Mb) and SSC12 (6.8 and 18.2 Mb) partially explained this relationship, having lower viral load (P < 0.0001); lower viremia at day 14 (P < 0.0001), day 21 (P < 0.01) and day 28 (P < 0.05) and greater overall average daily gain during infection (ADG_i; P < 0.01), ADG_i at week 3 (P < 0.001) and week 4 (P < 0.01). These additive genetic relationships could lead to molecular solutions to improve animal health and reduce production costs.     

Identification and characterization of nuclear microsatellite loci in the aquatic moss Platyhypnidium riparioides (Brachytheciaceae)

Eight microsatellite loci from the aquatic moss Platyhypnidium riparioides were identified using the method of microsatellite‐enriched libraries. Polymorphism was assessed in a sample of four populations of 20 individuals each from four streams of the Meuse hydrographic basin in southern Belgium. The markers amplified three to seven alleles per locus. Comparison of observed and expected heterozygosities as well as F‐statistics (F_ST = 0.62) reveals a significant genetic differentiation among populations. These markers will be useful for further investigation of population genetic structure and diversity at different nested spatial scales.     

Multiple origins and incursions of the Atlantic barnacle Chthamalus proteus in the Pacific

Chthamalus proteus, a barnacle native to the Caribbean and western Atlantic, was introduced to the Pacific within the last few decades. Using direct sequencing of mitochondrial DNA (COI), we characterized genetic variation in native and introduced populations and searched for genetic matches between regions to determine if there were multiple geographical sources and introduction points for this barnacle. In the native range, we found great genetic differences among populations (max. F_ST = 0.613) encompassing four lineages: one endemic to Panama, one endemic to Brazil, and two occurring Caribbean-wide. All four lineages were represented in the Pacific, but not equally; the Brazilian lineage was most prevalent and the Panamanian least common. Twenty-one individuals spread among nearly every island from where the barnacle is known in the Pacific, exactly matched six haplotypes scattered among Curaçao, the Netherlands Antilles; St John, US Virgin Islands; Puerto Rico; and Brazil, confirming a multigeographical origin for the Pacific populations. Significant genetic differences were also found in introduced populations from the Hawaiian Islands (F_CT = 0.043, P < 0.001), indicating introduction events have occurred at more than one locality. However, the sequence, timing and number of arrival events remains unknown. Possible reasons for limited transport of this barnacle through the Panama Canal are discussed. This and a preponderance of Brazilian-type individuals in the Pacific suggest an unexpected route of entry from around Cape Horn, South America. Unification in the Pacific of historically divergent lineages of this barnacle raises the possibility for selection of ‘hybrids’ with novel ecological adaptations in its new environment.     

Nitrogen fertilization rate affects feeding, larval performance, and oviposition preference of the beet armyworm, Spodoptera exigua, on cotton

Nitrogen (N) is one of the most critical chemical elements for plant and animal growth, exerting a variety of bottom‐up effects. Development and oviposition of the beet armyworm, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae), were studied in relation to varying N fertilization levels (42, 112, 196, and 280 p.p.m.) in cotton [Gossypium hirsutum L. (Malvaceae)]. Low N fertilization of cotton plants led to reduced plant biomass and a lower percentage of N in leaf blades and in leaf petioles. Development of S. exigua larvae fed with plants with reduced N applications (42 and 112 p.p.m.) was prolonged relative to treatments receiving higher N fertilization. Almost all larvae reared on artificial diets underwent only five instars before pupation. However, most larvae reared on cotton plants, irrespective of N levels, experienced a supernumerary sixth larval instar. Furthermore, significantly more larvae reared on lower N cotton plants underwent supernumerary development compared to larvae reared on higher N cotton plants. Life‐time feeding damage per larva ranged from 55 to 65 cm², depending on the nutritional quality of the food plant, although the differences were not statistically significant. Larvae distinguished between cotton plants with various nutritional qualities and fed preferentially on higher N plants. Female moth oviposition choice was also affected by host plant nutritional quality: cotton plants with higher N levels were preferentially chosen by S. exigua females for oviposition. The mechanisms of these effects are unclear, but they can have important implications for population dynamics and pest status of beet armyworms in the field.     

Evaluation of the causality of the low‐density lipoprotein receptor gene (LDLR) for serum lipids in pigs

A significant quantitative trait locus (QTL) for low‐density lipoprotein cholesterol (LDL‐C) and total cholesterol (TC) was identified around the LDLR gene on chromosome 2 (SSC2) in a White Duroc × Erhualian F₂ resource population and Sutai pigs in our previous study. However, in previous reports, the causality of LDLR with serum lipids is controversial in pigs. To systematically assess the causality of LDLR with serum lipids, association analyses were successively performed in three populations: Sutai pigs, a White Duroc × Erhualian F₂ resource population and a Duroc × (Landrace × Large White) population. We first performed a haplotype‐based association study with 60K SNP genotyping data and evidenced the significant association with LDL‐C and TC around the LDLR gene region. We also found that there is more than one QTL for LDL‐C and TC on SSC2. Then, we evaluated the causalities of two missense mutations, c.1812C>T and c.1520A>G, with LDL‐C and TC. We revealed that the c.1812C>T SNP showed the strongest association with LDL‐C (P = 5.40 × 10^?11) and TC (P = 3.64 × 10^?8) and explained all the QTL effect in Sutai pigs. Haplotype analysis found that two missense SNPs locate within a 1.93‐Mb haplotype block. One major haplotype showed the strongest significant association with LDL‐C (P = 4.62 × 10^?18) and TC (P = 1.06 × 10^?9). However, the c.1812C>T SNP was not identified in the White Duroc × Erhualian intercross, and the association of c.1520A>G with both LDL‐C and TC did not achieve significance in this F₂ population, suggesting population heterogeneity. Both missense mutations were identified in the Duroc × (Landrace × Large White) population and showed significant associations with LDL‐C and TC. Our data give evidence that the LDLR gene should be a candidate causative gene for LDL‐C and TC in pigs, but heterogeneity exists in different populations.     

Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano–Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d’Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6–3.0 alleles/locus; H_E = 0.25–0.42) and the fixation index (F_IS = −0.227–0.186) was not significantly different from zero suggesting that Hardy–Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q_ST = 0.055–0.283, Q_STmean = 0.189) strongly exceeds that in eight microsatellite loci (F_ST = −0.011–0.142, F_STmean = 0.047). Global and pairwise F_ST values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.     

Global mtDNA genetic structure and hypothesized invasion history of a major pest of citrus,Diaphorina citri (Hemiptera: Liviidae)

The Asian citrus psyllid Diaphorina citri Kuwayama is a key pest of citrus as the vector of the bacterium causing the “huanglongbing” disease (HLB). To assess the global mtDNA population genetic structure, and possible dispersal history of the pest, we investigated genetic variation at the COI gene collating newly collected samples with all previously published data. Our dataset consists of 356 colonies from 106 geographic sites worldwide. High haplotype diversity (H‐mean = 0.702 ± 0.017), low nucleotide diversity (π‐mean = 0.003), and significant positive selection (Ka/Ks = 32.92) were observed. Forty‐four haplotypes (Hap) were identified, clustered into two matrilines: Both occur in southeastern and southern Asia, North and South America, and Africa; lineages A and B also occur in eastern and western Asia, respectively. The most abundant haplotypes were Hap4 in lineage A (35.67%), and Hap9 in lineage B (41.29%). The haplotype network identified them as the ancestral haplotypes within their respective lineages. Analysis of molecular variance showed significant genetic structure (F_ST = 0.62, p < .0001) between the lineages, and population genetic analysis suggests geographic structuring. We hypothesize a southern and/or southeastern Asia origin, three dispersal routes, and parallel expansions of two lineages. The hypothesized first route involved the expansion of lineage B from southern Asia into North America via West Asia. The second, the expansion of some lineage A individuals from Southeast Asia into East Asia, and the third involved both lineages from Southeast Asia spreading westward into Africa and subsequently into South America. To test these hypotheses and gain a deeper understanding of the global history of D. citri, more data‐rich approaches will be necessary from the ample toolkit of next‐generation sequencing (NGS). However, this study may serve to guide such sampling and in the development of biological control programs against the global pest D. citri.     

Phylogeography of Chrysanthemum indicum L. (Compositae) in China based on trnL-F sequences

Using chloroplast DNA (cpDNA) trnL-F sequences, we studied the phylogeographic pattern and demographic history of Chrysanthemum indicum with 81 individuals from 27 populations representing the partial geographic range in China. Fifteen haplotypes were identified based on nucleotide variation. Relatively high levels of haplotype diversity (Hd = 0.783) and low levels of nucleotide diversity (Pi = 0.00521) were detected in C. indicum. An analysis of molecular variance (AMOVA) indicated that most of the genetic variation (59.9%, P < 0.001) was partitioned within populations. Genetic differentiation among populations was high (F_ST = 0.35), suggesting that there is little gene flow among populations. Nested clade analysis showed that 15 haplotypes formed twelve 1-step clades, five 2-step clades and two 3-step clades, with four clades showing significant geographical association. We concluded that contiguous range expansion and restricted gene flow/dispersal but with some long distance dispersal seemed to have shaped the present-day distribution of the haplotypes of C. indicum in China.