Detecting an elusive invasive species: a diagnostic PCR to detect Burmese python in Florida waters and an assessment of persistence of environmental DNA

Recent studies have demonstrated that detection of environmental DNA (eDNA) from aquatic vertebrates in water bodies is possible. The Burmese python, Python bivittatus, is a semi‐aquatic, invasive species in Florida where its elusive nature and cryptic coloration make its detection difficult. Our goal was to develop a diagnostic PCR to detect P. bivittatus from water‐borne eDNA, which could assist managers in monitoring this invasive species. First, we used captive P. bivittatus to determine whether reptilian DNA could be isolated and amplified from water samples. We also evaluated the efficacy of two DNA isolation methods and two DNA extraction kits commonly used in eDNA preparation. A fragment of the mitochondrial cytochrome b gene from P. bivittatus was detected in all water samples isolated with the sodium acetate precipitate and the QIAamp DNA Micro Kit. Next, we designed P. bivittatus‐specific primers and assessed the degradation rate of eDNA in water. Our primers did not amplify DNA from closely related species, and we found that P. bivittatus DNA was consistently detectable up to 96 h. Finally, we sampled water from six field sites in south Florida. Samples from five sites, where P. bivittatus has been observed, tested positive for eDNA. The final site was negative and had no prior documented evidence of P. bivittatus. This study shows P. bivittatus eDNA can be isolated from water samples; thus, this method is a new and promising technique for the management of invasive reptiles.     

Limited genetic structure and evidence for dispersal among populations of the endangered Florida grasshopper sparrow, Ammodramus savannarum floridanus

The Florida grasshopper sparrow, Ammodramus savannarum floridanus, is a non-migratory, endangered subspecies endemic to the prairie region of south-central Florida. It has experienced significant population declines and is currently restricted to five locations. We found substantial levels of variation in microsatellites and mtDNA control region sequences, estimates of inbreeding genetic effective population sizes that were much larger than the estimated census size, and no evidence of inbreeding within five sampled populations (n = 105). We also found a lack of genetic structure among populations (F _ST = 0.0123 for microsatellites and θ = 0.008 for mtDNA), and evidence for dispersal between populations, with 7.6% of all individuals identified as immigrants to their population of capture. We suggest that the subspecies be managed as a single management unit on a regional scale rather than as multiple management units on a local subpopulation scale. There is still a limited opportunity to preserve much of the present genetic variation in this subspecies, if immediate measures are taken to reverse the current population decline before this variation is reduced by genetic drift.     

The Genetic Integrity of the Ex Situ Population of the European Wildcat (Felis silvestris silvestris) Is Seriously Threatened by Introgression from Domestic Cats (Felis silvestris catus)

Studies on the genetic diversity and relatedness of zoo populations are crucial for implementing successful breeding programmes. The European wildcat, Felis s. silvestris, is subject to intensive conservation measures, including captive breeding and reintroduction. We here present the first systematic genetic analysis of the captive population of Felis s. silvestris in comparison with a natural wild population. We used microsatellites and mtDNA sequencing to assess genetic diversity, structure and integrity of the ex situ population. Our results show that the ex situ population of the European wildcat is highly structured and that it has a higher genetic diversity than the studied wild population. Some genetic clusters matched the breeding lines of certain zoos or groups of zoos that often exchanged individuals. Two mitochondrial haplotype groups were detected in the in situ populations, one of which was closely related to the most common haplotype found in domestic cats, suggesting past introgression in the wild. Although native haplotypes were also found in the captive population, the majority (68%) of captive individuals shared a common mtDNA haplotype with the domestic cat (Felis s. catus). Only six captive individuals (7.7%) were assigned as wildcats in the STRUCTURE analysis (at K = 2), two of which had domestic cat mtDNA haplotypes and only two captive individuals were assigned as purebred wildcats by NewHybrids. These results suggest that the high genetic diversity of the captive population has been caused by admixture with domestic cats. Therefore, the captive population cannot be recommended for further breeding and reintroduction.     

Microsatellite evidence for high clonality and limited genetic diversity in Ziziphus celata (Rhamnaceae), an endangered, self-incompatible shrub endemic to the Lake Wales Ridge, Florida, USA

Genetic data are often crucial for designing management strategies for rare and endangered species. Ziziphus celata is an endangered sandhill shrub endemic to the Lake Wales Ridge of central Florida. This self-incompatible clonal species is known from only 14 wild populations, most of which are small (under 100 plants). Focusing on the five populations discovered in 2007, we evaluate the level of genetic diversity and identify clonal lineages within the wild populations of the species with a set of microsatellite loci. To account for somatic mutations and genotyping errors, we identified clonal lineages using a threshold cutoff for pair-wise genetic distances among samples. The microsatellites had up to 18 alleles/locus, and, consistent with outcrossing, samples were highly heterozygous (average population level H _o  = 0.69). Most populations of Z. celata consist of a single clone, and the most diverse population has only 10 clones. Overall Z. celata comprises 41 multi-locus genotypes, and 30 clonal lineages. With nearly 1,000 recorded plants (595 genotyped) and only 30 clonal lineages, Ziziphus celata is highly clonal: clonal richness, R = 0.049. The pair-wise distance method facilitates identification of clonal lineages, avoiding overestimation of clonal diversity. In most cases, the samples that grouped into a lineage were one to four plants differing from a surrounding genotype by a single microsatellite repeat insertion/deletion mutation, consistent with these having arisen via somatic mutations. Our data will enable managers to incorporate extant diversity from wild populations into ex situ collections. Additionally, our research demonstrates the utility of microsatellites for conservation of imperiled species, identifying genotypes of high priority for preservation.     

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.     

High connectivity on a global scale in the pelagic wahoo,Acanthocybium solandri (tuna family Scombridae)

The population genetic structure and phylogeography of wahoo, Acanthocybium solandri, were investigated on a global scale with intron six of lactate dehydrogenase‐A (ldhA6, 8 locations, N = 213) and mtDNA cytochrome b (Cytb, 10 locations, N = 322). Results show extensive sharing of haplotypes across the wahoo's entire global range, and analyses were unable to detect significant structure (nuclear F_ST = 0.0125, P = 0.106; mtDNA Φ_ST < 0.0001, P = 0.634). Power analyses indicated 95% confidence in detecting nuclear F_ST ≥ 0.0389 and mtDNA Φ_ST ≥ 0.0148. These findings appear unique, as most other tunas, billfishes, and oceanic sharks exhibit significant population structure on the scale of East–West Atlantic, Atlantic vs. Indian‐Pacific, or East–West Pacific. Overall nuclear heterozygosity (H = 0.714) and mtDNA haplotype diversity (h = 0.918) are both high in wahoo, while overall mtDNA nucleotide diversity (π= 0.006) and nuclear nucleotide diversity (π=0.004) are uniformly low, indicating a recent increase in population size. Coalescence analyses yield an estimate of effective female population size (N_eF) at ~816 000, and a population bottleneck ~690 000 years ago. However, conclusions about population history from our Cytb data set are not concordant with a control region survey, a finding that will require further investigation. This is the first example of a vertebrate with a single globally distributed population, a finding we attribute to extensive dispersal at all life stages. The indications of a worldwide stock for wahoo reinforce the mandate for international cooperation on fisheries issues.     

Phylogeography of the Temminck’s Stint (Calidris temminckii): historical vicariance but little present genetic structure in a regionally endangered Palearctic wader

Aim We study the population differentiation and phylogeography of the Temminck’s Stint (Calidris temminckii). Specifically, we seek signs of past and present population size changes and dispersal events and evaluate management and conservation unit status of the populations. We also study the possibility of introgression as the origin of two mitochondrial DNA (mtDNA) lineages found and estimate the divergence time of the lineages. Location Northern Eurasia. Methods We analysed 583 bp of mtDNA control region domains I and II and 11 microsatellite loci from 13 localities throughout the breeding range. In addition, we used mitochondrial cytochrome c oxidase subunit I (COI), a barcoding gene, to search for signs of introgression. Results More population differentiation was found from microsatellites than from mtDNA, although differentiation was weak in both markers. Signs of past population growth were observed, in addition to more recent decline in some areas. Both control region and COI sequences revealed two maternal lineages coexisting in Fennoscandia and in north‐west Siberia. No signs of introgression were detected. Lineage divergence time was estimated to have occurred during the glacial periods of Pleistocene. Main conclusions Slight differences in mtDNA and microsatellite differentiation and diversity may reflect different features – such as the mutation rate and effective population size – of the markers used, or female‐biased dispersal pattern and high male site‐fidelity of the species. The coexistence of the two mitochondrial lineages is most likely a consequence of post‐glacial mixing of two refugial Pleistocene populations. Based on genetic information alone, global conservation concerns are not imminent. However, fast decline of a marginal Bothnian Bay population and the smallness and remoteness of a Central Yakutian population warrant conservation actions.     

Conservation genetics of the alligator snapping turtle: cytonuclear evidence of range-wide bottleneck effects and unusually pronounced geographic structure

A previous mtDNA study indicated that female-mediated gene flow was extremely rare among alligator snapping turtle populations in different drainages of the Gulf of Mexico. In this study, we used variation at seven microsatellite DNA loci to assess the possibility of male-mediated gene flow, we augmented the mtDNA survey with additional sampling of the large Mississippi River System, and we evaluated the hypothesis that the consistently low within-population mtDNA diversity reflects past population bottlenecks. The results show that dispersal between drainages of the Gulf of Mexico is rare (F _STmsat  = 0.43, Φ_STmtDNA = 0.98). Past range-wide bottlenecks are indicated by several genetic signals, including low diversity for microsatellites (1.1–3.9 alleles/locus; H _e = 0.06–0.53) and mtDNA (h = 0.00 for most drainages; π = 0.000–0.001). Microsatellite data reinforce the conclusion from mtDNA that the Suwannee River population might eventually be recognized as a distinct taxonomic unit. It was the only population showing fixation or near fixation for otherwise rare microsatellite alleles. Six evolutionarily significant units are recommended on the basis of reciprocal mtDNA monophyly and high levels of microsatellite DNA divergence.     

Genetic diversity of mitochondrial D-LOOP sequences in the spotted scat (Scatophagus argus) from different geographical populations along the northern coast of the South China Sea

The genetic diversity of 289 spotted scat (Scatophagus argus) from seven populations along the northern coast of the South China Sea was studied by analyzing the full-length sequences of the mitochondrial control region (D-LOOP). The S. argus D-LOOP sequence was 1,004–1,010 bp long and contained 156 variant sites. The seven studied S. argus populations had a high degree of genetic diversity (haplotype diversity [Hd] = 0.99135; nucleotide diversity (π) = 0.01313). There was no obvious genetic differentiation among the seven geographical populations and gene exchange was frequent (Fst = −0.01867–0.01117, p > .05). Four distinct mitochondrial lineages were identified in the phylogenetic tree and the haplotype network. The between-lineage Fst was 0.71690–0.84940 (p < .001), but these lineages showed no obvious phylogeographic pattern. Based on D-LOOP mutation rates, we estimated that the four lineages diverged approximately 513,800–93,600 years ago, during the Eocene ice age, at which time falling sea levels may have led to population segregation. We estimated that S. argus population expansion occurred approximately 2.29–0.68 million years ago, during the late Pleistocene. During this period, sea levels rose again, allowing previously separated lineages to come into sympatry, which eventually gave rise to a highly genetically diverse population without pyhlogeographic structure. Here, we characterized the genetic structure and differentiation of seven S. argus populations from the northern coast of the South China Sea. Our results suggested that the seven S. argus populations from the northern coast of the South China Sea have a relatively low level of genetic variation and can be considered a single unit for the purposes of fishery development, utilization, and management.     

Current status of the Sardinian partridge (<Emphasis Type="Italic">Alectoris barbara</Emphasis>) assessed by molecular markers

Four Alectoris species inhabit the Mediterranean area, where they represent important gamebirds subject to human manipulations. The Sardinian partridge is peculiar in Europe, in that it belongs to the African species Alectoris barbara. Nevertheless, no comprehensive study has as yet investigated its genetic status as regards both the extant levels of genetic diversity and the possible contamination due to introgressive hybridization with other Mediterranean species. For the purposes of this study, we analyzed 65 samples of Sardinian partridges, 40 of which came from the wild population and 25 from captive stocks. No one of them showed a mtDNA polymerase chain reaction restriction fragment length polymorphism haplotype assigned to another species than A. barbara, thus, ruling out a possible introgression in the maternal line. In addition, we compared these samples with 94 partridges from other circum-Mediterranean populations using a set of eight chicken (Gallus gallus) microsatellites. A low level of genetic variation was observed in the Sardinian population (H _E = 0.310; k _AR = 2.69), comparable only to that observed in the Sicilian rock partridge (A. graeca). The comparison with the Tunisian population showed that its present genetic composition is consistent with a historical introduction from North Africa, showing possible effects of a post-introductional genetic drift. Bayesian tests assigned all but one individuals with >90% probability to A. barbara, thus, providing evidence that no or only a few exotic Alectoris genes have introgressed into Sardinian partridges.     

Evaluation of cytochrome b mtDNA sequences in genetic diversity studies of <Emphasis Type="Italic">Channa marulius</Emphasis> (Channidae: Perciformes)

Channa marulius (Hamilton, 1822) is a commercially important freshwater fish and a potential candidate species for aquaculture. The present study evaluated partial Cytochrome b gene sequence of mtDNA for determining the genetic variation in wild populations of C. marulius. Genomic DNA extracted from C. marulius samples (n = 23) belonging to 3 distant rivers; Mahanadi, Teesta and Yamuna was analyzed. Sequencing of 307 bp Cytochrome b mtDNA fragment revealed the presence of 5 haplotypes with haplotype diversity value of 0.763 and nucleotide diversity value of 0.0128. Single population specific haplotype was observed in Mahanadi and Yamuna samples and 3 haplotypes in Teesta samples. The analysis of data demonstrated the suitability of partial Cytochrome b sequence in determining the genetic diversity in C. marulius population.     

Environmental DNA (eDNA) Sampling Improves Occurrence and Detection Estimates of Invasive Burmese Pythons

Environmental DNA (eDNA) methods are used to detect DNA that is shed into the aquatic environment by cryptic or low density species. Applied in eDNA studies, occupancy models can be used to estimate occurrence and detection probabilities and thereby account for imperfect detection. However, occupancy terminology has been applied inconsistently in eDNA studies, and many have calculated occurrence probabilities while not considering the effects of imperfect detection. Low detection of invasive giant constrictors using visual surveys and traps has hampered the estimation of occupancy and detection estimates needed for population management in southern Florida, USA. Giant constrictor snakes pose a threat to native species and the ecological restoration of the Florida Everglades. To assist with detection, we developed species-specific eDNA assays using quantitative PCR (qPCR) for the Burmese python (Python molurus bivittatus), Northern African python (P. sebae), boa constrictor (Boa constrictor), and the green (Eunectes murinus) and yellow anaconda (E. notaeus). Burmese pythons, Northern African pythons, and boa constrictors are established and reproducing, while the green and yellow anaconda have the potential to become established. We validated the python and boa constrictor assays using laboratory trials and tested all species in 21 field locations distributed in eight southern Florida regions. Burmese python eDNA was detected in 37 of 63 field sampling events; however, the other species were not detected. Although eDNA was heterogeneously distributed in the environment, occupancy models were able to provide the first estimates of detection probabilities, which were greater than 91%. Burmese python eDNA was detected along the leading northern edge of the known population boundary. The development of informative detection tools and eDNA occupancy models can improve conservation efforts in southern Florida and support more extensive studies of invasive constrictors. Generic sampling design and terminology are proposed to standardize and clarify interpretations of eDNA-based occupancy models.     

Extensive sampling of polar bears (Ursus maritimus) in the Northwest Passage (Canadian Arctic Archipelago) reveals population differentiation across multiple spatial and temporal scales

As global warming accelerates the melting of Arctic sea ice, polar bears (Ursus maritimus) must adapt to a rapidly changing landscape. This process will necessarily alter the species distribution together with population dynamics and structure. Detailed knowledge of these changes is crucial to delineating conservation priorities. Here, we sampled 361 polar bears from across the center of the Canadian Arctic Archipelago spanning the Gulf of Boothia (GB) and M'Clintock Channel (MC). We use DNA microsatellites and mitochondrial control region sequences to quantify genetic differentiation, estimate gene flow, and infer population history. Two populations, roughly coincident with GB and MC, are significantly differentiated at both nuclear (F_ST = 0.01) and mitochondrial (Φ_ST = 0.47; F_ST = 0.29) loci, allowing Bayesian clustering analyses to assign individuals to either group. Our data imply that the causes of the mitochondrial and nuclear genetic patterns differ. Analysis of mtDNA reveals the matrilineal structure dates at least to the Holocene, and is common to individuals throughout the species’ range. These mtDNA differences probably reflect both genetic drift and historical colonization dynamics. In contrast, the differentiation inferred from microsatellites is only on the scale of hundreds of years, possibly reflecting contemporary impediments to gene flow. Taken together, our data suggest that gene flow is insufficient to homogenize the GB and MC populations and support the designation of GB and MC as separate polar bear conservation units. Our study also provide a striking example of how nuclear DNA and mtDNA capture different aspects of a species demographic history.     

Influence of drift and admixture on population structure of American black bears (Ursus americanus) in the Central Interior Highlands,USA, 50 years after translocation

Bottlenecks, founder events, and genetic drift often result in decreased genetic diversity and increased population differentiation. These events may follow abundance declines due to natural or anthropogenic perturbations, where translocations may be an effective conservation strategy to increase population size. American black bears (Ursus americanus) were nearly extirpated from the Central Interior Highlands, USA by 1920. In an effort to restore bears, 254 individuals were translocated from Minnesota, USA, and Manitoba, Canada, into the Ouachita and Ozark Mountains from 1958 to 1968. Using 15 microsatellites and mitochondrial haplotypes, we observed contemporary genetic diversity and differentiation between the source and supplemented populations. We inferred four genetic clusters: Source, Ouachitas, Ozarks, and a cluster in Missouri where no individuals were translocated. Coalescent models using approximate Bayesian computation identified an admixture model as having the highest posterior probability (0.942) over models where the translocation was unsuccessful or acted as a founder event. Nuclear genetic diversity was highest in the source (A_R = 9.11) and significantly lower in the translocated populations (A_R = 7.07–7.34; P = 0.004). The Missouri cluster had the lowest genetic diversity (A_R = 5.48) and served as a natural experiment showing the utility of translocations to increase genetic diversity following demographic bottlenecks. Differentiation was greater between the two admixed populations than either compared to the source, suggesting that genetic drift acted strongly over the eight generations since the translocation. The Ouachitas and Missouri were previously hypothesized to be remnant lineages. We observed a pretranslocation remnant signature in Missouri but not in the Ouachitas.     

Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.     

The mitochondrial genome sequence analysis of Ophidascaris baylisi from the Burmese python (Python molurus bivittatus)

Ophidascaris species are parasitic roundworms that inhabit the python gut, resulting in severe granulomatous lesions or even death. However, the classification and nomenclature of these roundworms are still controversial. Our study aims to identify a snake roundworm from the Burmese python (Python molurus bivittatus) and analyze the mitochondrial genome. We identified this roundworm as Ophidascaris baylisi based on the morphology and cytochrome c oxidase subunit I (cox1) sequence. Ophidascaris baylisi complete mitochondrial genome was 14,784 bp in length, consisting of two non-coding regions and 36 mitochondrial genes (12 protein-coding genes, 22 tRNA genes, and two rRNA genes). The protein-coding genes used TTG, ATG, ATT, or TTA as start codons and TAG, TAA, or T as stop codons. All tRNA genes showed a TV-loop structure, except trnS1^AGN and trnS2^UCN revealed a D-loop structure. The mitochondrial large ribosomal subunit 16S (rrnL) and small ribosomal subunit 12S (rrnS) were 956 bp and 700 bp long, respectively. Phylogenetic analysis based on O. baylisi mitochondrial protein-coding genes demonstrated that O. baylisi clustered with the family Ascarididae members and was most closely related to Ophidascaris wangi. These results may enhance the nematode mitochondrial genome database and provide valuable molecular markers for further research on the taxonomy, phylogeny, and genetic relationships of Ophidascaris nematodes.     

Genetic diversity within a threatened,endemic North American species, <Emphasis Type="Italic">Euphorbia telephioides</Emphasis> (Euphorbiaceae)

The southeastern United States and Florida support an unusually large number of endemic plant species, many of which are threatened by anthropogenic habitat disturbance. As conservation measures are undertaken and recovery plans designed, a factor that must be taken into consideration is the genetic composition of the species of concern. Here we describe the levels, and partitioning, of genetic diversity in 17 populations of the rare and threatened Florida endemic, Euphorbia telephioides (telephus spurge). Species-wide genetic diversity was high (P_s = 91%, AP_s = 3.81, A_s = 3.57 and H_es = 0.352) as was mean population genetic diversity (P_p = 81%, AP_p = 2.98, A_p = 2.59 and H_ep = 0.320) which ranks it among the highest 10% of plant species surveyed. Partitioning of genetic variation (G_st = 0.106) was low compared to other herbaceous outcrossing perennials indicating high historical gene flow across its limited geographic range. Among population G_st values within the three Florida counties in which it occurs, Gulf (0.084), Franklin (0.059) and Bay Counties (0.033), were also quite low. Peripheral populations did not generally have reduced genetic variation although there was significant isolation by distance. Rarefaction analysis showed a non-significant relationship between allelic richness and actual population sizes. Our data suggest that E. telephioides populations were probably more continuously distributed in Bay, Gulf and Franklin Counties and that their relative contemporary isolation is a recent phenomenon. These results are important for developing a recovery plan for this species.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Changes in genetic diversity and differentiation in Red‐cockaded woodpeckers (Dryobates borealis) over the past century

Red‐cockaded woodpeckers (RCW; Dryobates borealis) declined after human activities reduced their fire‐maintained pine ecosystem to <3% of its historical range in the southeastern United States and degraded remaining habitat. An estimated 1.6 million RCW cooperative breeding groups declined to about 3,500 groups with no more than 10,000 birds by 1978. Management has increased RCW population abundances since they were at their lowest in the 1990s. However, no range‐wide study has been undertaken since then to investigate the impacts of this massive bottleneck or infer the effects of conservation management and recent demographic recoveries. We used mitochondrial DNA sequences (mtDNA) and nine nuclear microsatellite loci to determine if range‐wide demographic declines resulted in changes to genetic structure and diversity in RCW by comparing samples collected before 1970 (mtDNA data only), between 1992 and 1995 (mtDNA and microsatellites), and between 2010 and 2014 (mtDNA and microsatellites). We show that genetic diversity has been lost as detected by a reduction in the number of mitochondrial haplotypes. This reduction was apparent in comparisons of pre‐1970 mtDNA data with data from the 1992–1995 and 2010–2014 time points, with no change between the latter two time points in mtDNA and microsatellite analyses. The mtDNA data also revealed increases in range‐wide genetic differentiation, with a genetically panmictic population present throughout the southeastern United States in the pre‐1970s data and subsequent development of genetic structure that has remained unchanged since the 1990s. Genetic structure was also uncovered with the microsatellite data, which like the mtDNA data showed little change between the 1992–1995 and 2010–2014 data sets. Temporal haplotype networks revealed a consistent, star‐like phylogeny, suggesting that despite the overall loss of haplotypes, no phylogenetically distinct mtDNA lineages were lost when the population declined. Our results may suggest that management during the last two decades has prevented additional losses of genetic diversity.     

Population genetic structure and intraspecific genetic distance of Periplaneta americana (Blattodea: Blattidae) based on mitochondrial and nuclear DNA markers

The American cockroach (Periplaneta americana) is a globally invasive pest that can cause significant economic loss and threaten human health. Although it is abundant and lives in close proximity to humans, few studies have investigated the genetic diversity of P. americana. Our study analyzed 1,053 P. americana and other Periplaneta species' samples from different locations in China and the United States. A traditional tree‐based method using 17 unique mitochondrial COI haplotypes of P. americana and 20 haplotypes of the other Periplaneta species accurately identified P. americana with a barcoding threshold of 5.1%. To identify the population genetic structure of P. americana, we investigated wingless gene and pooled them with obtained mtDNA data for a combined analysis. Although the genetic diversity of the USA group was relatively higher than the China group, the number of haplotypes and alleles of both groups was small. The analysis of molecular variance (AMOVA), intraspecific phylogeny, and haplotype networks indicated that P. americana had very little global genetic differentiation. The weak geographic genetic structure might reflect the human‐mediated dispersal of P. americana. Despite no apparent phylogeographic assignment of mtDNA and nuclear lineages was observed in both BI trees, the integrated COI sequence data identified four distinct P. americana haplotype groups, showing four ancient maternal lineages of P. americana in China and the United States.