Genetic identification of two sibling species of Lutzomyia longipalpis (Diptera: Psychodidae) that produce distinct male sex pheromones in Sobral,Ceará State,Brazil

Lutzomyia longipalpis, the main sandfly vector for New World visceral leishmaniasis is a complex of an as yet undefined number of sibling species. At present, there is no consensus on the status (single species vs. species complex) of Brazilian populations. We applied five microsatellite loci to test the hypothesis that L. longipalpis occurs as two sympatric cryptic species in Sobral, Ceará State, Brazil as predicted by male sex pheromone chemotypes described previously for field specimens from this site [S-9-methyl-germacrene-B (9MGB) and a cembrene compound]. Abdominal spot morphology corresponds with pheromone type at this locality (9MGB in '1 spot' males and cembrene in '2 spot' males). Genotype data from 190 wild-caught L. longipalpis specimens collected in October 1999 and April 2001 were used to estimate genetic differentiation between the two sex pheromone populations and sampling dates. No significant (P > 0.05) genetic differences were found between the 1999 and 2001 9MGB samples (theta = 0.018; RST = -0.005), and genetic differentiation was low between the cembrene collections (theta = 0.037, P < 0.05; RST = -0.043, P > 0.05). By contrast, highly divergent allelic frequencies (largely at two microsatellite loci) corresponded to significant (P > 0.05) genetic differentiation (theta = 0.221; RST = 0.215) for all comparisons between samples with different pheromones. When pheromone samples were pooled across sample date, genetic differentiation was high (theta = 0.229; P < 0.001; Nem = 0.84). The allele frequency distribution at each of the five microsatellite loci was similar for males and females from the two collection years. Two of these loci showed highly divergent allele frequencies in the two sex pheromone populations. This was reflected in the highly significant genetic differentiation obtained from the male genotypes, between populations producing different pheromones (theta = 0.229-0.268; P < 0.0001 for the 2001 and theta = 0.254-0.558; P < 0.0001 for the 1999 collections, respectively). Similar results were obtained when the females, assigned to a pheromone type, were included in the analysis. Both a Bayesian analysis of the data set and a population assignment test provided strong evidence for two distinct populations corresponding to pheromone type. Given its genotype, the probability of assigning a 9MGB male to the original 9MGB population was 100% once the two years' collections were pooled. For cembrene-producing '2 spot' males this probability although still high, was lower than for 9MGB males, at 86%. This microsatellite data together with previously reported reproductive isolation between the two Sobral populations confirm that premating barriers are important in speciation of L. longipalpis.     

Local scale population genetic structure of Entomophthora muscae epidemics

Isolates of the entomopathogenic fungus Entomophthora muscae were obtained from houseflies sampled at five stables on Zealand during an epidemic in fall 2005. DNA fingerprints were generated from single conidia from 40 E. muscae isolates using the PCR-based method of Inter-Simple Sequence Repeats (ISSR). This resulted in fingerprint patterns consisting of about 50 fragments, of which, 14 were polymorphic. From the polymorphic loci we estimated the reproductive mode, genetic differentiation and gene diversity of E. muscae populations using the statistics of index of association, Weir & Cockerhams theta and Nei's analysis of gene diversity in subdivided populations. Our results revealed no significant differences in allele frequencies among the five populations. The index of association test rejected the null hypothesis of random mating, but the test for paired locus compatibility showed weak or no linkage between loci indicating recombination or homoplasy in the dataset. From this study we cannot exclude the possibility that the population genetic structure underlying E. muscae epidemics could be panmictic consisting of several lineages with a high level of reciprocal migration.     

Influence of northern limit range on genetic diversity and structure in a widespread North American tree,sugar maple (Acer saccharum Marshall)

Due to climate change, the ranges of many North American tree species are expected to shift northward. Sugar maple (Acer saccharum Marshall) reaches its northern continuous distributional limit in northeastern North America at the transition between boreal mixed‐wood and temperate deciduous forests. We hypothesized that marginal fragmented northern populations from the boreal mixed wood would have a distinct pattern of genetic structure and diversity. We analyzed variation at 18 microsatellite loci from 23 populations distributed along three latitudinal transects (west, central, and east) that encompass the continuous–discontinuous species range. Each transect was divided into two zones, continuous (temperate deciduous) and discontinuous (boreal mixed wood), based on sugar maple stand abundance. Respective positive and negative relationships were found between the distance of each population to the northern limit (D_north), and allelic richness (A_R) and population differentiation (F_ST). These relations were tested for each transect separately; the pattern (discontinuous–continuous) remained significant only for the western transect. structure analysis revealed the presence of four clusters. The most northern populations of each transect were assigned to a distinct group. Asymmetrical gene flow occurred from the southern into the four northernmost populations. Southern populations in Québec may have originated from two different postglacial migration routes. No evidence was found to validate the hypothesis that northern populations were remnants of a larger population that had migrated further north of the species range after the retreat of the ice sheet. The northernmost sugar maple populations possibly originated from long‐distance dispersal.     

A high-resolution map of Arabidopsis recombinant inbred lines by whole-genome exon array hybridization

Recombinant populations were the basis for Mendel's first genetic experiments and continue to be key to the study of genes, heredity, and genetic variation today. Genotyping several hundred thousand loci in a single assay by hybridizing genomic DNA to oligonucleotide arrays provides a powerful technique to improve precision linkage mapping. The genotypes of two accessions of Arabidopsis were compared by using a 400,000 feature exon-specific oligonucleotide array. Around 16,000 single feature polymorphisms (SFPs) were detected in ~8,000 of the ~26,000 genes represented on the array. Allelic variation at these loci was measured in a recombinant inbred line population, which defined the location of 815 recombination breakpoints. The genetic linkage map had a total length of 422.5 cM, with 676 informative SFP markers representing intervals of ~0.6 cM. One hundred fifteen single gene intervals were identified. Recombination rate, SFP distribution, and segregation in this population are not uniform. Many genomic regions show a clustering of recombination events including significant hot spots. The precise haplotype structure of the recombinant population was defined with unprecedented accuracy and resolution. The resulting linkage map allows further refinement of the hundreds of quantitative trait loci identified in this well-studied population. Highly variable recombination rates along each chromosome and extensive segregation distortion were observed in the population.     

Genetic structure of populations of Rhizoctonia solani AG-3 on potato in eastern North Carolina

A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) method was developed to identify and differentiate genotypes of Rhizoctonia solani anastomosis group 3 subgroup PT (AG-3 PT), a fungal pathogen of potato. Polymorphic co-dominant single-locus PCR-RFLP markers were identified after sequencing of clones from a genomic library and digestion with restriction enzymes. Multilocus genotypes were determined by a combination of PCR product and digestion with a specific restriction enzyme for each of seven loci. A sample of 104 isolates from one commercial field in each of five counties in eastern North Carolina was analyzed, and evidence for high levels of gene flow between populations was revealed. When data were clone-corrected and samples pooled into one single North Carolina population, random associations of alleles were found for all loci or pairs of loci, indicating random mating. However, when all genotypes were analyzed, the observed genotypic diversity deviated from panmixia and alleles within and between loci were not randomly associated. These findings support a model of population structure for R. solani AG-3 PT on potato that includes both recombination and clonality.     

Multiscale spatial genetic structure within and between populations of wild cherry trees in nuclear genotypes and chloroplast haplotypes

Spatial genetic structure (SGS) of plants mainly depends on the effective population size and gene dispersal. Maternally inherited loci are expected to have higher genetic differentiation between populations and more intensive SGS within populations than biparentally inherited loci because of smaller effective population sizes and fewer opportunities of gene dispersal in the maternally inherited loci. We investigated biparentally inherited nuclear genotypes and maternally inherited chloroplast haplotypes of microsatellites in 17 tree populations of three wild cherry species under different conditions of tree distribution and seed dispersal. As expected, interpopulation genetic differentiation was 6–9 times higher in chloroplast haplotypes than in nuclear genotypes. This difference indicated that pollen flow 4–7 times exceeded seed flow between populations. However, no difference between nuclear and chloroplast loci was detected in within‐population SGS intensity due to their substantial variation among the populations. The SGS intensity tended to increase as trees became more aggregated, suggesting that tree aggregation biased pollen and seed dispersal distances toward shorter. The loss of effective seed dispersers, Asian black bears, did not affect the SGS intensity probably because of mitigation of the bear loss by other vertebrate dispersers and too few tree generations after the bear loss to alter SGS. The findings suggest that SGS is more variable in smaller spatial scales due to various ecological factors in local populations.     

Microsatellite loci isolated from the tropical tree Hymenaea courbaril L. (Fabaceae)

We developed 11 microsatellite markers for Hymenaea courbaril for the purpose of studying spatial genetic structure and gene flow. The microsatellite loci were screened in 44 trees from two populations. All loci were polymorphic, exhibiting between two and 16 alleles, and levels of expected heterozygosity from 0.174 to 0.909. Departures from Hardy-Weinberg equilibrium were detected for all loci in one population. The estimated null allele frequency is low or moderate. No locus combinations exhibited linkage disequilibrium.     

Spatial genetic structure among and within populations of Primula sieboldii growing beside separate streams

We investigated the hierarchical genetic structure of SSR (simple sequence repeats) and cpDNA (chloroplast DNA) polymorphisms among and within populations of Primula sieboldii, a heterostylous clonal herb. Seven out of eight populations at the study site, located in a mountainous region of Nagano Prefecture, had each developed alongside a different stream, and the other occurred on a flat area 70 m from the nearest stream. The magnitude of genetic differentiation among streamside populations in maternally inherited cpDNA (Phi = 0.341) was much higher than that in biparentally inherited SSRs (Phi = 0.011). This result suggests that seed dispersal among streams was restricted, and pollen was the primary agent of gene flow among streamside populations. In contrast, genetic differentiation among subpopulations within streams were low at both markers (Phi = 0.053 for cpDNA, Phi = 0.025 for SSR). This low differentiation among subpopulations in cpDNA compared with that among streamside populations suggest that seed dispersal occur along the stream probably during flooding. This hypothesis was supported by the fact that in cpDNA haplotypes, no clear genetic structure was detected within the streamside population, while a significant genetic structure was found within 20 m in the nonstreamside population. Furthermore, within the streamside populations, two pairs of ramets with identical multilocus genotypes for eight SSR loci were distantly (> 50 m) distributed along the same streamside, suggesting dispersal of clonal propagule. Our study showed that the heterogeneity of the landscape can influence gene flow and hence spatial genetic structure.     

Nuclear gene sequences and DNA variation of Cryptomeria japonica samples from the postglacial period

Genomic DNA was extracted from heartwood blocks of six Cryptomeria japonica individuals that had been buried (in an area now covered by rice fields) for about 3600 years. Attempts were made to determine the sequences of five nuclear genes following polymerase chain reaction amplification, using previously obtained C. japonica expressed sequence tag (EST) information. We detected 15 nucleotide substitutions and four insertion/deletions (indels) in a partial GapC gene sequence among 13 individuals of the buried and an extant population, which allowed us to estimate the extent of DNA variation within the buried populations, and the level of genetic differentiation between the buried population and the extant population growing in a neighbouring area. For the entire haplotypes of the GapC region, pi and theta nucleotide diversity estimates were 0.0063 and 0.0010, respectively, when both populations were included, while corresponding figures for the buried population alone were 0.0009 and 0.0017. Estimates of DNA divergence statistics (dXY = 0.0062, dA = 0.0005, FST = 0.0832 and KST = 0.0935) suggest that differentiation between the two populations was not great. However, permutation tests gave FST and KST values rejecting the null hypothesis (that populations were not differentiated) at the 5% and 1% probability levels, respectively. The significant genetic differentiation between the two populations was mainly caused by differences in haplotype diversity. The significant level of haplotype diversity in the extant population compared to the buried population might be the result of gene flow from neighbouring artificial forests. Alternatively, it is possible that we failed to detect all the DNA variation in the buried population because of clonal growth in the buried population.     

Genetic linkage map of Coffea canephora: effect of segregation distortion and analysis of recombination rate in male and female meioses.

Two complementary segregating plant populations of Coffea canephora were produced from the same clone. One population (DH) comprised 92 doubled haploids derived from female gametes, while the other population (TC) was a test cross consisting of 44 individuals derived from male gametes. Based on the DH population, a genetic linkage map comprising 160 loci was constructed. Eleven linkage groups that putatively correspond to the 11 gametic chromosomes of C. canephora were identified. The mapped loci included more than 40 specific sequence-tagged site markers, either single-copy RFLP probes or microsatellites, that could serve as standard landmarks in coffee-genome analyses. Furthermore, comparisons for segregation distortion and recombination frequency between the two populations were performed. Although segregation distortions were observed in both populations, the frequency of loci exhibiting a very pronounced degree of distortion was especially high in the DH population. This observation is consistent with the hypothesis of strong zygotic selection among the DH population. The recombination frequencies in both populations were found to be almost indistinguishable. These results offer evidence in favour of the lack of significant sex differences in recombination in C. canephora.     

Mating patterns and genetic diversity in the wild daffodil Narcissus longispathus (Amaryllidaceae)

Despite the importance of Narcissus to ornamental horticulture, there have been no population genetic studies of wild species, many of which have narrow distributions. Here, we measure selfing rates and levels of genetic diversity at allozyme loci in six populations of Narcissus longispathus, a self-compatible daffodil endemic to a few mountain ranges in southeastern Spain. The populations were distributed among four distinct river valleys encompassing two main watersheds in the Sierra de Cazorla mountains. Selfing rates averaged 0.37 (range 0.23-0.46), resulting in significant inbreeding coefficients for the progeny (f = 0.324). In contrast, estimates of inbreeding in parental genotypes were not significantly different from zero (f = 0.001), indicating that few selfed offspring survive to maturity because of inbreeding depression. Species-wide estimates of genetic diversity for the six populations were P(s) = 0.38, H(es) = 0.119 and A(s) = 1.27 with significant genetic differentiation among populations theta = 0.15. The observed patterns of genetic differentiation among populations are likely influenced by the mating system, and a combination of local topography, watershed affinities and gene flow.     

Population Genetics of Euphydryas Butterflies. I. Genetic Variation and the Neutrality Hypothesis

Twenty-one populations of the checkerspot butterfly, Euphydryas editha, and ten populations of Euphydryas chalcedona were sampled for genetic variation at eight polymorphic enzyme loci. Both species possessed loci that were highly variable from population to population and loci that were virtually identical across all populations sampled. Our data indicate that the neutrality hypothesis is untenable for the loci studied, and therefore selection is indicated as the major factor responsible for producing these patterns. Thorough ecological work allowed gene flow to be ruled out (in almost all instances) as a factor maintaining similar gene frequencies across populations. The Lewontin-Krakauer test indicated magnitudes of heterogeneity among standardized variances of gene frequencies inconsistent with the neutrality hypothesis. The question of whether or not to correct this statistic for sample size is discussed. Observed equitability of gene frequencies of multiple allelic loci was found to be greater than that predicted under the neutrality hypothesis. Genetic differentiation persisting through two generations was found between the one pair of populations known to exchange significant numbers of individuals per generation. Two matrices of genetic distance between populations, based on the eight loci sampled, were found to be significantly correlated with a matrix of environmental distance, based on measures of fourteen environmental parameters. Correlations between gene frequencies and environmental parameters, results of multiple regression analysis, and results of principle component analysis showed strong patterns of association and of "explained" variation. The correlation analyses suggest which factors might be further investigated as proximate selective agents.     

Beyond the point of no return? A comparison of genetic diversity in captive and wild populations of two nearly extinct species of Goodeid fish reveals that one is inbred in the wild

The relative importance of genetic and non-genetic factors in extinction liability has been extensively debated. Here, we examine the levels of genetic variability at 13 (seven informative) loci in wild and captive populations of two endangered species of Mexican Goodeid fish, Ameca splendens and Zoogoneticus tequila. Allelic diversity was higher in the wild populations, and F(IS) lower. Values of theta (=4Nemu) were estimated using a coalescent approach. These implied that the effective population size of all captive populations of A. splendens were smaller than that of the wild population; qualitatively similar results were obtained using an analytical method based on within-population gene identity disequilibrium. However, the wild population of Z. tequila did not show a significantly greater estimate of theta. We used the Beaumont approach to infer population declines, and found that both species showed clear evidence of a decline in effective population size, although this was stronger and probably occurred over a longer period of time in Z. tequila than in A. splendens. The decline in Z. tequila probably occurred before captive populations were established. We discuss implications for the conservation of critically endangered populations.     

Absence of population genetic differentiation in the New Zealand greenshell mussel Perna canaliculus (Gmelin 1791) as assessed by allozyme variation

Genetic variation in the endemic New Zealand greenshell mussel, Perna canaliculus (Gmelin 1791), was examined using starch-gel electrophoresis at seven protein-coding loci (Idh; Acon-1; Acon-2; Gpd; Pgi; Pgm; Pgd) in 35 populations (N=1038 mussels). For all loci and all populations, Fisher's exact tests indicated highly significant departures from Hardy-Weinberg equilibrium (HWE), but this overall result was caused by significant heterozygote deficiencies at only two loci (Pgm and Pgd), and in only three northern populations (Kuaotunu, Te Haumi and Days Bay). Allelic and genotypic differentiation between population pairs at individual loci and across all loci were nonsignificant, and genotypic disequilibrium at each locus pair was also nonsignificant for all populations. Genetic variation in all populations was high (mean heterozygosity, 0.210+/-0.113), while Nei's D among populations was very low (0.002+/-0.002). Low population subdivision (θ=-0.001+/-0.002) and high levels of gene flow (Nm(p)=10.18; Nm(θ)=infinity) also indicated that the single panmictic unit model best explains population genetic homogeneity in P. canaliculus over a north-south distance >2000 km. Lack of genetic subdivision in this species is discussed in light of two previous allozyme studies, with differing results: one suggested that a north-south division exists between greenshell mussel stocks, and the other suggested that population structure in this species can be explained through isolation by distance model modified by local hydrology.     

Genetic analysis reveals population structuring and a bottleneck in the black-faced lion tamarin (Leontopithecus caissara)

The ability of a population to evolve in a changing environment may be compromised by human-imposed barriers to gene flow. We investigated the population structure and the possible occurrence of a genetic bottleneck in two isolated populations of the black-faced lion tamarin (Leontopithecus caissara), a species with very reduced numbers (less than 400) in a very restricted range in the Atlantic Forest of southeast Brazil. We determined the genotypes of 52 individuals across 9 microsatellite loci. We found genetic divergence between the populations, each exhibiting low genetic diversity. Analysis revealed broad- and fine-scale population structuring. Both populations have evidently experienced population reduction and a genetic bottleneck without presenting any apparent detrimental effect. Anyway, measures should be taken to effectively protect the forests where L. caissara occurs in order to allow its populations to increase and counteract the eventual effects of genetic impoverishment.     

Genetic diversity and population structure of two important medicinal plant species <Emphasis Type="Italic">Schisandra chinensis</Emphasis> and <Emphasis Type="Italic">Schisandra sphenanthera</Emphasis> revealed by nuclear microsatellites

Seven polymorphic and transferable nuclear microsatellites were used to investigate the population structure of genetic diversity of Schisandra chinensis and Schisandra sphenanthera for facilitating their conservation and sustainable utilization. High levels of gene diversity were revealed in these two medicinal species, the majority of genetic diversity was harbored within populations, and population structure was might due to restricted gene flow among populations. Isolation by distance was close to significance in S. chinensis but not in S. sphenanthera. In S. chinensis, null alleles were identified as a cause for excess of homozygotes at loci G24 and WGA60, but inbreeding might also be partly responsible for the positive F _IS values in this species. In contrast, null allele frequencies were high at all the seven loci in S. sphenanthera and resulted in overestimation of fixation index. The strategy for ex situ conservation of these two medicinal species is discussed based on the genetic results.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

尼罗罗非鱼ghrelin基因的多态性及其与生长性状相关SNP位点的筛选

为了研究尼罗罗非鱼(Oreochromis niloticus)生长激素促分泌素基因(ghrelin)的多态性及其与生长的相关性, 研究以两个尼罗罗非鱼群体(快长群体和基础群体)的DNA样本各40份为模板, 通过PCR扩增和测序获得ghrelin基因序列。通过Dnasp v5和MEGA 5.0分析序列多态性、筛选有效SNP 位点; 采用Snapshot法对两个群体子代ghrelin基因中SNP位点进行基因分型, 然后分析SNP位点基因型与生长性状的相关性。结果表明, 快长群体ghrelin基因中的单核苷酸变异位点数(S)比基础群体要少, 而核苷酸多态性(Pi)和平均核苷酸差异数(K)要略高于基础群体。共筛得3个有效SNP 位点(S1、S2和S3), 均分布于第1个内含子中。遗传结构分析表明, 3个SNP 位点在两个群体的子代中均为低度多态性位点(PIC0.25), 但处于Hardy-Weinberg平衡(P0.05);快长群体子代中3个SNP 位点的观测杂合度、期望杂合度和多态信息含量等遗传多样性参数均小于基础群体子代的相应值, 3个SNP 位点的遗传多样性参数、基因型和基因频率在同一群体中高度一致, SNP 位点之间完全连锁。两个群体子代中3个SNP 位点处的优势基因型相同, 但快长群体子代中优势基因型频率要明显大于基础群体子代中相应基因型频率。对两个群体子代的生长性状与SNP基因型进行关联性分析的结果表明,尼罗罗非鱼个体的多项生长指标(体重、体长、体高、头长和尾柄高等)在不同基因型中存在显著差异(S1:GG AG, S2:TT AT, S3:AA AT)(P0.05)。D1双倍型(S1:GG, S2:TT, S3:AA)所对应的尼罗罗非鱼个体的多项生长指标(体重、体长、体高、头长和尾柄高等)显著高于D2双倍型(S1:AG, S2:AT, S3:AT)。以上结果表明, 尼罗罗非鱼ghrelin基因3个SNP 位点完全连锁, D1双倍型与快长性状密切相关, 可作为尼罗罗非鱼分子标记辅助育种的候选标记。     

Variability of the Ramularia collo‐cygni Population in Central Europe

The genetic structure of the fungal barley pathogen Ramularia collo‐cygni (Rcc) population in Central Europe involving the isolates from the Czech Republic, the Slovak Republic, Germany and Swiss was determined using amplified fragment length polymorphism (AFLP) analysis. One hundred and eighty‐four markers were chosen to determine genetic and genotypic diversity and to test the hypothesis of random mating and population differentiation of Rcc isolates. Among the 337 isolates collected, the overall gene diversity was moderate ( = 0.216). The level of multilocus genotypic diversity was higher within populations than among them. All individuals had unique multilocus genotypes. Genetic differentiation was significant among populations in localities, but at a moderate level (θ = 0.12; P < 0.001), suggesting that gene flow is occurring among populations. The isolates from all twelve clusters produced by Structure were found in all local populations, although at different frequencies. Therefore, the inferred clusters did not represent geographical populations. Although the null hypothesis of random mating in Rcc populations was rejected, the high level of genotypic diversity suggests that the Rcc population structure appears to be generated by a mixed reproductive system including both asexual and sexual reproduction, along with a rather high migration rate.     

Inter‐island divergence within Drosophila mauritiana,a species of the D. simulans complex: Past history and/or speciation in progress?

Speciation with gene flow may be more common than generally thought, which makes detailed understanding of the extent and pattern of genetic divergence between geographically isolated populations useful. Species of the Drosophila simulans complex provide a good model for speciation and evolutionary studies, and hence understanding their population genetic structure will increase our understanding of the context in which speciation has occurred. Here, we describe genetic diversity and genetic differentiation of two distant populations of D. mauritiana (Mauritius and Rodrigues Islands) at mitochondrial and nuclear loci. We surveyed the two populations for their mitochondrial haplotypes, eight nuclear genes and 18 microsatellite loci. A new mitochondrial type is fixed in the Rodrigues population of D. mauritiana. The two populations are highly differentiated, their divergence appears relatively ancient (100,000 years) compared to the origin of the species, around 0.25MYA, and they exhibit very limited gene flow. However, they have similar levels of divergence from their sibling, D. simulans. Both nuclear genes and microsatellites revealed contrasting demographic histories between the two populations, expansion for the Mauritius population and stable population size for the Rodrigues Island population. The discovery of pronounced geographic structure within D. mauritiana combined to genetic structuring and low gene flow between the two island populations illuminates the evolutionary history of the species and clearly merits further attention in the broad context of speciation.