Genetic Diversity in Introduced Golden Mussel Populations Corresponds to Vector Activity

We explored possible links between vector activity and genetic diversity in introduced populations of Limnoperna fortunei by characterizing the genetic structure in native and introduced ranges in Asia and South America. We surveyed 24 populations: ten in Asia and 14 in South America using the mitochondrial cytochrome c oxidase subunit I (COI) gene, as well as eight polymorphic microsatellite markers. We performed population genetics and phylogenetic analyses to investigate population genetic structure across native and introduced regions. Introduced populations in Asia exhibit higher genetic diversity (H _E = 0.667–0.746) than those in South America (H _E = 0.519–0.575), suggesting higher introduction effort for the former populations. We observed pronounced geographical structuring in introduced regions, as indicated by both mitochondrial and nuclear markers based on multiple genetic analyses including pairwise Ф_ST, F _ST, Bayesian clustering method, and three-dimensional factorial correspondence analyses. Pairwise F _ST values within both Asia (F _ST = 0.017–0.126, P = 0.000–0.009) and South America (F _ST = 0.004–0.107, P = 0.000–0.721) were lower than those between continents (F _ST = 0.180–0.319, P = 0.000). Fine-scale genetic structuring was also apparent among introduced populations in both Asia and South America, suggesting either multiple introductions of distinct propagules or strong post-introduction selection and demographic stochasticity. Higher genetic diversity in Asia as compared to South America is likely due to more frequent propagule transfers associated with higher shipping activities between source and donor regions within Asia. This study suggests that the intensity of human-mediated introduction vectors influences patterns of genetic diversity in non-indigenous species.     

Plasmodium vivax Diversity and Population Structure across Four Continents

Plasmodium vivax is the geographically most widespread human malaria parasite. To analyze patterns of microsatellite diversity and population structure across countries of different transmission intensity, genotyping data from 11 microsatellite markers was either generated or compiled from 841 isolates from four continents collected in 1999–2008. Diversity was highest in South-East Asia (mean allelic richness 10.0–12.8), intermediate in the South Pacific (8.1–9.9) Madagascar and Sudan (7.9–8.4), and lowest in South America and Central Asia (5.5–7.2). A reduced panel of only 3 markers was sufficient to identify approx. 90% of all haplotypes in South Pacific, African and SE-Asian populations, but only 60–80% in Latin American populations, suggesting that typing of 2–6 markers, depending on the level of endemicity, is sufficient for epidemiological studies. Clustering analysis showed distinct clusters in Peru and Brazil, but little sub-structuring was observed within Africa, SE-Asia or the South Pacific. Isolates from Uzbekistan were exceptional, as a near-clonal parasite population was observed that was clearly separated from all other populations (F _ST>0.2). Outside Central Asia F _ST values were highest (0.11–0.16) between South American and all other populations, and lowest (0.04–0.07) between populations from South-East Asia and the South Pacific. These comparisons between P. vivax populations from four continents indicated that not only transmission intensity, but also geographical isolation affect diversity and population structure. However, the high effective population size results in slow changes of these parameters. This persistency must be taken into account when assessing the impact of control programs on the genetic structure of parasite populations.     

Genetic Population Structure of Macridiscus multifarius (Mollusca: Bivalvia) on the Basis of Mitochondrial Markers: Strong Population Structure in a Species with a Short Planktonic Larval Stage

The clam Macridiscus multifarius with a planktonic larval stage of about 10 days is an ecologically and economically important species in the coastal regions of China. In this study, 3 mt-DNA markers (COI, 12S rRNA, and ND1) were used to investigate the population structure and demography of wild M. multifarius populations in 3 coastal localities of the East China Sea (ZS and ZP populations) and Beibu Gulf in the South China Sea (BH population). Sequences of 685 bp in COI, 350 bp in 12S rRNA, and 496 bp in ND1 were determined. High level and significant F _ST values were obtained among the different localities on the basis of either COI (F _ST = 0.100–0.444, p < 0.05) or 12S rRNA (F _ST = 0.199–0.742, p < 0.05) gene, indicating a high degree of genetic differentiation among the populations. F _ST values were significant but weak for the ND1 gene because it is highly conservative. The median-joining network suggested an obvious genetic differentiation between ZS and BH populations, and the finding is consistent with the results of our demographic analyses using the unweighted pair group method with arithmetic mean. Our study unraveled the extant population genetic structure of M. multifarius and explained the strong population structure of a species with a short planktonic larval stage species; this information could be useful for fishery management measures, including artificial breeding and conservation.     

Annual cycle of the legume pod borer Maruca vitrata Fabricius (Lepidoptera: Crambidae) in southwestern Burkina Faso

Maruca vitrata is an economically significant insect pest of cowpea in sub-Saharan Africa. Understanding the seasonal population patterns of M. vitrata is essential for the establishment of effective pest management strategies. M. vitrata larval populations on cultivated cowpea and adult flying activities were monitored in addition to scouting for host plants and parasitoids during 2 consecutive years in 2010 and 2011 in southwestern Burkina Faso. Our data suggest that M. vitrata populations overlapped on cultivated cowpea and alternate host plants during the rainy season. During the cowpea off-season, M. vitrata maintained a permanent population on the wild host plants Mucuna poggei and Daniella oliveri. The parasitoid fauna include three species, Phanerotoma leucobasis Kri., Braunsia kriegeri End. and Bracon sp. Implications of these finding for pest management strategies are discussed.     

Worldwide Distribution of the MYH9 Kidney Disease Susceptibility Alleles and Haplotypes: Evidence of Historical Selection in Africa

MYH9 was recently identified as renal susceptibility gene (OR 3–8, p<10⁻⁸) for major forms of kidney disease disproportionately affecting individuals of African descent. The risk haplotype (E-1) occurs at much higher frequencies in African Americans (≥60%) than in European Americans (<4%), revealing a genetic basis for a major health disparity. The population distributions of MYH9 risk alleles and the E-1 risk haplotype and the demographic and selective forces acting on the MYH9 region are not well explored. We reconstructed MYH9 haplotypes from 4 tagging single nucleotide polymorphisms (SNPs) spanning introns 12–23 using available data from HapMap Phase II, and by genotyping 938 DNAs from the Human Genome Diversity Panel (HGDP). The E-1 risk haplotype followed a cline, being most frequent within sub-Saharan African populations (range 50–80%), less frequent in populations from the Middle East (9–27%) and Europe (0–9%), and rare or absent in Asia, the Americas, and Oceania. The fixation indexes (F_ST) for pairwise comparisons between the risk haplotypes for continental populations were calculated for MYH9 haplotypes; F_ST ranged from 0.27–0.40 for Africa compared to other continental populations, possibly due to selection. Uniquely in Africa, the Yoruba population showed high frequency extended haplotype length around the core risk allele (C) compared to the alternative allele (T) at the same locus (rs4821481, iHs = 2.67), as well as high population differentiation (F_{ST(CEU vs. YRI)} = 0.51) in HapMap Phase II data, also observable only in the Yoruba population from HGDP (F_ST = 0.49), pointing to an instance of recent selection in the genomic region. The population-specific divergence in MYH9 risk allele frequencies among the world''s populations may prove important in risk assessment and public health policies to mitigate the burden of kidney disease in vulnerable populations.     

Plio-Pleistocene history and phylogeography of Acacia senegal in dry woodlands and savannahs of sub-Saharan tropical Africa: evidence of early colonisation and recent range expansion

Drylands are extensive across sub-Saharan Africa, socio-economically and ecologically important yet highly sensitive to environmental changes. Evolutionary history, as revealed by contemporary intraspecific genetic variation, can provide valuable insight into how species have responded to past environmental and population changes and guide strategies to promote resilience to future changes. The gum arabic tree (Acacia senegal) is an arid-adapted, morphologically diverse species native to the sub-Saharan drylands. We used variation in nuclear sequences (internal transcribed spacer (ITS)) and two types of chloroplast DNA (cpDNA) markers (PCR-RFLP, cpSSR) to study the phylogeography of the species with 293 individuals from 66 populations sampled across its natural range. cpDNA data showed high regional and rangewide haplotypic diversity (h_T(cpSSR)=0.903–0.948) and population differentiation (G_ST(RFLP)=0.700–0.782) with a phylogeographic pattern that indicated extensive historical gene flow via seed dispersal. Haplotypes were not restricted to any of the four varieties, but showed significant geographic structure (G_ST(cpSSR)=0.392; R_ST=0.673; R_ST>R_ST (permuted)), with the major division separating East and Southern Africa populations from those in West and Central Africa. Phylogenetic analysis of ITS data indicated a more recent origin for the clade including West and Central African haplotypes, suggesting range expansion in this region, possibly during the Holocene humid period. In conjunction with paleobotanical evidence, our data suggest dispersal to West Africa, and across to the Arabian Peninsula and Indian subcontinent, from source populations located in the East African region during climate oscillations of the Plio-Pleistocene.     

Geographic distribution of phylogenetically-distinct legume pod borer, Maruca vitrata (Lepidoptera: Pyraloidea: Crambidae)

Maruca vitrata Fabricius is a pantropical lepidopteran pest of legumes. Phylogenetic analysis of a mitochondrial cytochrome c oxidase-I gene (cox1) fragment indicates that three Maruca sp. mitochondrial lineages have unique geographic distributions [lineages 1 and 2: Australia, Taiwan, and West Africa (Niger, Nigeria, and Burkina Faso), and lineage 3: Puerto Rico]. The haplotype (T30, T114) is specific to lineages 1&2 and was assayed by NsiI and SacI polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) within population samples; it was not observed in the Puerto Rican samples, but was nearly fixed among samples from West Africa, Australia and Taiwan (85.5?C100%). Re-sequencing and phylogenetic analyses of PCR-RFLP defined cox1 haplotypes indicate that nucleotide diversity is highest among samples from West Africa. Phylogenetic reconstruction based upon ribosomal DNA (rDNA) internal transcribed spacer-2 (ITS-2) sequences provided additional evidence for three Maruca sp. clades. These data suggest that multiple unique Maruca species or subspecies are present worldwide, which has implications for the management of this pest species-complex.     

Plasmodium falciparum mitochondrial genetic diversity exhibits isolation-by-distance patterns supporting a sub-Saharan African origin

The geographical distribution of single nucleotide polymorphism (SNP) in the mitochondrial genome of the human malaria parasite Plasmodium falciparum was investigated. We identified 88 SNPs in 516 isolates from seven parasite populations in Africa, Southeast Asia and Oceania. Analysis of the SNPs postulated a sub-Saharan African origin and recovered a strong negative correlation between within-population SNP diversity and geographic distance from the putative African origin over Southeast Asia and Oceania. These results are consistent with those previously obtained for nuclear genome-encoded housekeeping genes, indicating that the pattern of inheritance does not substantially affect the geographical distribution of SNPs.     

Features of evolution and expansion of modern humans,inferred from genomewide microsatellite markers

We study data on variation in 52 worldwide populations at 377 autosomal short tandem repeat loci, to infer a demographic history of human populations. Variation at di-, tri-, and tetranucleotide repeat loci is distributed differently, although each class of markers exhibits a decrease of within-population genetic variation in the following order: sub-Saharan Africa, Eurasia, East Asia, Oceania, and America. There is a similar decrease in the frequency of private alleles. With multidimensional scaling, populations belonging to the same major geographic region cluster together, and some regions permit a finer resolution of populations. When a stepwise mutation model is used, a population tree based on TD estimates of divergence time suggests that the branches leading to the present sub-Saharan African populations of hunter-gatherers were the first to diverge from a common ancestral population (approximately 71-142 thousand years ago). The branches corresponding to sub-Saharan African farming populations and those that left Africa diverge next, with subsequent splits of branches for Eurasia, Oceania, East Asia, and America. African hunter-gatherer populations and populations of Oceania and America exhibit no statistically significant signature of growth. The features of population subdivision and growth are discussed in the context of the ancient expansion of modern humans.     

Multi-Gene Analysis Reveals a Lack of Genetic Divergence between Calanus agulhensis and C. sinicus (Copepoda; Calanoida)

The discrimination and taxonomic identification of marine species continues to pose a challenge despite the growing number of diagnostic metrics and approaches. This study examined the genetic relationship between two sibling species of the genus Calanus (Crustacea; Copepoda; Calanidae), C. agulhensis and C. sinicus, using a multi-gene analysis. DNA sequences were determined for portions of the mitochondrial cytochrome c oxidase I (mtCOI); nuclear citrate synthase (CS), and large subunit (28S) rRNA genes for specimens collected from the Sea of Japan and North East (NE) Pacific Ocean for C. sinicus and from the Benguela Current and Agulhas Bank, off South Africa, for C. agulhensis. For mtCOI, C. sinicus and C. agulhensis showed similar levels of haplotype diversity (H_d = 0.695 and 0.660, respectively) and nucleotide diversity (π = 0.003 and 0.002, respectively). Pairwise F_ST distances for mtCOI were significant only between C. agulhensis collected from the Agulhas and two C. sinicus populations: the Sea of Japan (F_ST = 0.152, p<0.01) and NE Pacific (F_ST = 0.228, p<0.005). Between the species, F_ST distances were low for both mtCOI (F_ST = 0.083, p = 0.003) and CS (F_ST = 0.050, p = 0.021). Large subunit (28S) rRNA showed no variation between the species. Our results provide evidence of the lack of genetic distinction of C. sinicus and C. agulhensis, raise questions of whether C. agulhensis warrants status as a distinct species, and indicate the clear need for more intensive and extensive ecological and genetic analysis.     

F(ST) and Q(ST) under neutrality

A commonly used test for natural selection has been to compare population differentiation for neutral molecular loci estimated by F_ST and for the additive genetic component of quantitative traits estimated by Q_ST. Past analytical and empirical studies have led to the conclusion that when averaged over replicate evolutionary histories, Q_ST = F_ST under neutrality. We used analytical and simulation techniques to study the impact of stochastic fluctuation among replicate outcomes of an evolutionary process, or the evolutionary variance, of Q_ST and F_ST for a neutral quantitative trait determined by n unlinked diallelic loci with additive gene action. We studied analytical models of two scenarios. In one, a pair of demes has recently been formed through subdivision of a panmictic population; in the other, a pair of demes has been evolving in allopatry for a long time. A rigorous analysis of these two models showed that in general, it is not necessarily true that mean Q_ST = F_ST (across evolutionary replicates) for a neutral, additive quantitative trait. In addition, we used finite-island model simulations to show there is a strong positive correlation between Q_ST and the difference Q_ST − F_ST because the evolutionary variance of Q_ST is much larger than that of F_ST. If traits with relatively large Q_ST values are preferentially sampled for study, the difference between Q_ST and F_ST will also be large and positive because of this correlation. Many recent studies have used tests of the null hypothesis Q_ST = F_ST to identify diversifying or uniform selection among subpopulations for quantitative traits. Our findings suggest that the distributions of Q_ST and F_ST under the null hypothesis of neutrality will depend on species-specific biology such as the number of subpopulations and the history of subpopulation divergence. In addition, the manner in which researchers select quantitative traits for study may introduce bias into the tests. As a result, researchers must be cautious before concluding that selection is occurring when Q_ST ≠ F_ST.     

The joint effects of selection and dominance on the QST - FST contrast

Q_ST measures the differentiation of quantitative traits between populations. It is often compared to F_ST, which measures population differentiation at neutral marker loci due to drift, migration, and mutation. When Q_ST is different from F_ST, it is usually taken as evidence that selection has either restrained or accelerated the differentiation of the quantitative trait relative to neutral markers. However, a number of other factors such as inbreeding, dominance, and epistasis may also affect the Q_ST − F_ST contrast. In this study, we examine the effects of dominance, selection, and inbreeding on Q_ST − F_ST. We compare Q_ST with F_ST at selected and neutral loci for populations at equilibrium between selection, drift, mutation, and migration using both analytic and simulation approaches. Interestingly, when divergent selection is acting on a locus, inbreeding and dominance generally inflate Q_ST relative to F_ST when they are both measured at the quantitative locus at equilibrium. As a consequence, dominance is unlikely to hide the signature of divergent selection on the Q_ST − F_ST contrast. However, although in theory dominance and inbreeding affect the expectation for Q_ST − F_ST, of most concern is the very large variance in both Q_ST and F_ST, suggesting that we should be cautious in attributing small differences between Q_ST and F_ST to selection.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

The Distribution and Host Shifts of Cotton-Melon Aphids in Northern China

Aphis gossypii Glover (Hemiptera: Aphididae) is a serious pest of cotton in northern China. A microsatellite analysis was used to characterize the genetic structure of A. gossypii populations from different geographic, host plant, and seasonal populations in 2014. Among 906 individuals, 507 multilocus genotypes were identified, with genotypic richness values of 0.07–1.00 for the populations. We observed moderate levels of genetic differentiation among geographic populations (F_ST = 0.103; 95% confidence interval: 0.065–0.145) and host plant populations (F_ST = 0.237; 95% confidence interval: 0.187–0.296). A Mantel test of isolation by distance revealed no significant correlations between Slatkin’s linearized F_ST and the natural logarithm of geographic distance. A Bayesian analysis of population genetic structures identified three clusters. An analysis of molecular variance revealed significant differences among the three clusters (F = 0.26596, P < 0.0001), among seasons (F = 0.04244, P = 0.00381), and among host populations (F = 0.12975, P = 0.0029). Thus, the A. gossypii populations in northern China exhibit considerable genotypic diversity. Additionally, our findings indicated that the 31 analyzed populations could be classified as one of three host biotypes (i.e., cotton, cucumber, and pomegranate biotypes). There were also clear seasonal effects on population genetic structure diversity among aphids collected from Anyang.     

Phylogenetics,delimitation and historical biogeography of the pantropical tree genus Thespesia (Malvaceae,Gossypieae)

Thespesia consists of 16 species of trees and shrubs from Southeast Asia–Oceania, Africa and America, the most well known being T. populnea, a small tree of tropical coastal areas around the world. Phylogenetic relationships in the genus and among its allies in tribe Gossypieae were inferred using three plastid and two nuclear regions to ascertain its generic delimitation and explore its biogeographical history. Maximum‐likelihood and Bayesian analyses confirmed that Thespesia is not monophyletic and, based on these results, Azanza is reinstated to accommodate the two species previously placed in Thespesia section Lampas. Dating analyses and ancestral range estimation indicated that Thespesia s.s. most likely originated in Southeast Asia–Oceania c. 30 Mya, but extant species did not begin to differentiate until the late Miocene. Two dispersal events, one into Africa c. 11 Mya and another into America (Antilles) c. 9 Mya, gave rise to the African and the Greater Antillean endemics, respectively. The two most widespread hydrochorous species, T. populnea and T. populneoides, originated in Southeast Asia–Oceania from where they spread to other parts of the world. Our analysis also indicated a much earlier origin than previously reported for the eumalvoid clade and its tribes Gossypieae, Malveae and Hibisceae suggesting that vicariance might have had an important role early in the history of these groups.     

Genetic differentiation and species cohesion in two widespread Central American Begonia species

Begonia is one of the ten largest plant genera, with over 1500 species. This high species richness may in part be explained by weak species cohesion, which has allowed speciation by divergence in allopatry. In this study, we investigate species cohesion in the widespread Central American Begonia heracleifolia and Begonia nelumbiifolia, by genotyping populations at microsatellite loci. We then test for post-zygotic reproductive barriers using experimental crosses, and assess whether sterility barriers are related to intraspecific changes in genome size, indicating major genome restructuring between isolated populations. Strong population substructure was found for B. heracleifolia (F_ST=0.364, F′_ST=0.506) and B. nelumbiifolia (F_ST=0.277, F′_ST=0.439), and Bayesian admixture analysis supports the division of most populations into discrete genetic clusters. Moderate levels of inferred selfing (B. heracleifolia s=0.40, B. nelumbiifolia s=0.62) and dispersal limitation are likely to have contributed to significant genetic differentiation (B. heracleifolia Jost''s D=0.274; B. nelumbiifolia D=0.294). Interpopulation crosses involving a divergent B. heracleifolia population with a genome size ∼10% larger than the species mean had a ∼20% reduction in pollen viability compared with other outcrosses, supporting reproductive isolation being polymorphic within the species. The population genetic data suggest that Begonia populations are only weakly connected by gene flow, allowing reproductive barriers to accumulate between the most isolated populations. This supports allopatric divergence in situ being the precursor of speciation in Begonia, and may also be a common speciation mechanism in other tropical herbaceous plant groups.     

Geographic Influence on Genetic Structure in the Widespread Neotropical Tree Simarouba amara (Simaroubaceae)

This study aimed to assess the population genetic structure of a widespread Neotropical tree species, Simarouba amara, at local, regional and continental spatial scales. We used five microsatellite loci to examine genetic variation in 14 natural populations (N?=?478 individuals) of this vertebrate-dispersed rain forest tree species in Panama, Ecuador, and French Guiana. Estimates of genetic differentiation (F _st and R _st) were significant among all but one population pair and global differentiation was moderate (F _st?=?0.25, R _st?=?0.33) with 94% of genetic variation ascribed to differences among three main geographic regions (Central America, Western Ecuador, Amazon basin). There was no evidence of isolation by distance within regions. Allele-size mutations contributed significantly (R _ST > F _ST) to the divergences between cis- and trans-Andean populations, highlighting the role of the northern Andean cordilleras as an important geographic barrier for this species.     

Human Population Differentiation Is Strongly Correlated with Local Recombination Rate

Allele frequency differences across populations can provide valuable information both for studying population structure and for identifying loci that have been targets of natural selection. Here, we examine the relationship between recombination rate and population differentiation in humans by analyzing two uniformly-ascertained, whole-genome data sets. We find that population differentiation as assessed by inter-continental F _ST shows negative correlation with recombination rate, with F _ST reduced by 10% in the tenth of the genome with the highest recombination rate compared with the tenth of the genome with the lowest recombination rate (P≪10⁻¹²). This pattern cannot be explained by the mutagenic properties of recombination and instead must reflect the impact of selection in the last 100,000 years since human continental populations split. The correlation between recombination rate and F _ST has a qualitatively different relationship for F _ST between African and non-African populations and for F _ST between European and East Asian populations, suggesting varying levels or types of selection in different epochs of human history.