Contrasting patterns of population structure and demographic history in cryptic species of Bostrychia intricata (Rhodomelaceae,Rhodophyta) from New Zealand

Spatial patterns of genetic diversity provide insight into the demography and history of species. Morphologically similar but genetically distinct “cryptic” species are increasingly being recognized in marine organisms through molecular analyses. Such species are, on closer inspection, often discovered to display contrasting life histories or occasionally minor morphological differences; molecular tools can thus be useful indicators of diversity. Bostrychia intricata, a marine red alga, is widely distributed throughout the Southern Hemisphere and comprises many cryptic species. We used mitochondrial cytochrome c oxidase I gene sequences to assess the genetic variation, population genetic structure, and demographic history of B. intricata in New Zealand. Our results supported the existence of three cryptic species of B. intricata (N2, N4, and N5) in New Zealand. Cryptic species N4, which was found throughout New Zealand, showed a higher genetic diversity and wider distribution than the other two species, which were only found in the North Island and northern South Island. Our analyses showed low to moderate genetic differentiation among eastern North Island populations for cryptic species N2, but high differentiation among North and South Island populations for N4, suggesting different population structure between these cryptic species. Data also indicated that N2 has recently undergone population expansion, probably since the Last Glacial Maximum (LGM), while the higher genetic diversity in N4 populations suggests persistence in situ through the LGM. The contrasting population structures and inferred demographic histories of these species highlight that life history can vary greatly even among morphologically indistinguishable taxa.     

Population structure,population history and DNA barcoding of fruit fly Bactrocera latifrons (Hendel) (Diptera: Tephritidae)

The fruit fly Bactrocera latifrons (Hendel) is an important pest of commercially significant plants such as chili, tomato and eggplant. The species is native to South and Southeast Asia, but has now invaded Japan, Hawaii and Africa. In this study, mitochondrial DNA sequences were used to infer genetic structure and demographic history of B. latifrons. The efficiency of DNA barcodes for identification of B. latifrons was also tested. Ninety‐three specimens infesting four host‐plant species were obtained from 11 sampling locations in Thailand. The mitochondrial haplotype network revealed no major divergent lineage, which was consistent with a phylogenetic analysis that found strong support for the monophyly of B. latifrons. Population pairwise F_ST revealed that most (65%) comparisons were not significantly different, suggesting a high rate of gene flow. Analysis of molecular variance (amova ) found no significant genetic differentiation among populations from different host‐plant species. Sharing of several haplotypes among flies from different host‐plants indicates that the flies were moved freely across the plant species. Demographic history analysis revealed that the population has undergone recent expansion dating back to the end of the last glaciation. Thus, the results indicate that both ongoing and historical factors have played important roles in determining the genetic structure and diversity of B. latifrons. DNA barcoding analysis revealed that B. latifrons specimens were clearly differentiated from other species with 100% correct identification. Therefore, cytochrome oxidase I (COI) barcoding sequences could be effectively used to identify this important pest species, which could encourage monitoring and control efforts for this species.     

Genetic diversity and population structure of the commercially harvested sea urchin Paracentrotus lividus (Echinodermata, Echinoidea)

The population structure of the edible Atlanto-Mediterranean sea urchin Paracentrotus lividus is described by analysing sequence variation in a fragment of the mitochondrial gene cytochrome c oxidase subunit I in 127 individuals from 12 localities across south-west Europe. The study revealed high levels of genetic diversity but low levels of genetic structure, suggesting a large degree of gene flow between populations and panmixis within each, the Mediterranean and Atlantic basins. However, we found significant genetic differentiation between the two basins probably due to restricted gene flow across the geographical boundary imposed by the area of the Strait of Gibraltar. Populations of P. lividus appeared to have experienced a recent demographic expansion in the late Pleistocene. We provide new evidence on the population structure of this commercial species, predicting a healthy stock of this sea urchin on the Mediterranean and Atlantic coasts.     

Spatial expansion and population structure of the neotropical malaria vector, Anopheles darlingi (Diptera: Culicidae)

Extensive population structuring is known to occur in Anopheles darlingi , the primary malaria vector of the Neotropics. We analysed the phylogeographic structure of the species using the mitochondrial cytochrome oxidase I marker. Diversity is divided into six main population groups in South America: Colombia, central Amazonia, southern Brazil, south-eastern Brazil, and two groups in north-east Brazil. The ancestral distribution of the taxon is hypothesized to be central Amazonia, and there is evidence of expansion from this region during the late Pleistocene. The expansion was not a homogeneous front, however, with at least four subgroups being formed due to geographic barriers. As the species spread, populations became isolated from each other by the Amazon River and the coastal mountain ranges of south-eastern Brazil and the Andes. Analyses incorporating distances around these barriers suggest that the entire South American range of An. darlingi is at mutation–dispersal–drift equilibrium. Because the species is distributed throughout such a broad area, the limited dispersal across some landscape types promotes differentiation between otherwise proximate populations. Moreover, samples from the An. darlingi holotype location in Rio de Janeiro State are substantially derived from all other populations, implying that there may be additional genetic differences of epidemiological relevance. The results obtained contribute to our understanding of gene flow in this species and allow the formulation of human mosquito health protocols in light of the potential population differences in vector capacity or tolerance to control strategies.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 854–866.     

Ecological speciation by temporal isolation in a population of the stonefly Leuctra hippopus (Plecoptera,Leuctridae)

Stream dwelling invertebrates are ideal candidates for the study of ecological speciation as they are often adapted to particular environmental conditions within a stream and inhabit only certain reaches of a drainage basin, separated by unsuitable habitat. We studied an atypical population of the stonefly Leuctra hippopus at a site in central Norway, the Isterfoss rapids, in relation to three nearby and two remote conspecific populations. Adults of this population emerge about a month earlier than those of nearby populations, live on large boulders emerging from the rapids, and are short‐lived. This population also has distinct morphological features and was studied earlier during the period 1975–1990. We reassessed morphological distinctness with new measurements and added several analyses of genetic distinctness based on mitochondrial and nuclear sequence markers, as well as AFLP fingerprinting and SNPs mined from RAD sequences. The Isterfoss population is shown to be most closely related to its geographical neighbors, yet clearly morphologically and genetically distinct and homogeneous. We conclude that this population is in the process of sympatric speciation, with temporal isolation being the most important direct barrier to gene flow. The shift in reproductive season results from the particular temperature and water level regime in the Isterfoss rapids. The distinct adult body shape and loss of flight are hypothesized to be an adaptation to the unusual habitat. Ecological diversification on small spatial and temporal scales is one of the likely causes of the high diversity of aquatic insects.     

Species status and population structure of the Australian Eucalyptus pest Paropsis atomaria Olivier (Coleoptera: Chrysomelidae)

1 Paropsis atomaria Olivier represents an emergent pest of Eucalyptus plantations in Queensland and New South Wales, Australia. Most prior studies on the biology and control of P. atomaria have centred on populations from Canberra in the Australian Capital Territory, but the biological relationship between beetles from Canberra and those from up to 1500 km further north are unknown. 2 DNA markers were used to determine whether P. atomaria from Canberra are the same biological species as those from Eucalyptus forestry plantations in northern New South Wales and Queensland, where the beetle has become an important pest. Using the mitochondrial gene, cytochrome c oxidase I (COI), individuals collected from across the distribution of P. atomaria were investigated for haplotype diversity and levels of mitochondrial divergence. 3 Within P. atomaria, genetic distance averaged 0.5% across 23 unique haplotypes for 93 individuals, with an average of 14% difference between P. atomaria and the outgroup species, Paropsis obsoleta. Significant genetic structure was observed relative to geographical distribution, but not with respect to host plant species of origin. Greatest divergence was between the southern‐most sample site (Canberra) and northern sites in New South Wales and Queensland, indicating reduced gene flow between these regions. 4 Individuals from across eastern Australia belong to the same genetic species with population substructuring evident. Consequently, there is no evidence to suggest cryptic species complexes exist within the currently defined taxon. Continued implementation of control strategies for P. atomaria across its distribution is appropriate.     

Identifying the true oysters (Bivalvia: Ostreidae) with mitochondrial phylogeny and distance-based DNA barcoding

Oysters (family Ostreidae), with high levels of phenotypic plasticity and wide geographic distribution, are a challenging group for taxonomists and phylogenetics. As a useful tool for molecular species identification, DNA barcoding might offer significant potential for oyster identification and taxonomy. This study used two mitochondrial fragments, cytochrome c oxidase I (COI) and the large ribosomal subunit (16S rDNA), to assess whether oyster species could be identified by phylogeny and distance-based DNA barcoding techniques. Relationships among species were estimated by the phylogenetic analyses of both genes, and then pairwise inter- and intraspecific genetic divergences were assessed. Species forming well-differentiated clades in the molecular phylogenies were identical for both genes even when the closely related species were included. Intraspecific variability of 16S rDNA overlapped with interspecific divergence. However, average intra- and interspecific genetic divergences for COI were 0-1.4% (maximum 2.2%) and 2.6-32.2% (minimum 2.2%), respectively, indicating the existence of a barcoding gap. These results confirm the efficacy of species identification in oysters via DNA barcodes and phylogenetic analysis.     

Genetic diversity of Paramecium species on the basis of multiple loci analysis and ITS secondary structure models

Among ciliates, Paramecium has become a privileged model for the study of “species problem” particularly in the case of the “Paramecium aurelia complex” that has been intensely investigated. Despite extensive studies, the taxonomy of Paramecium is still challenging. The major problem is an uneven sampling of Paramecium with relatively few representatives of each species. To investigate species from the less discovered region (Pakistan), 10 isolates of Paramecium species including a standing-alone FT8 strain previously isolated by some of us were subjected to molecular characterization. Fragments of 18S recombinant DNA (rDNA), ITS1-5.8S-ITS2-5′LSU rDNA, cytochrome c oxidase subunit II, and hsp70 genes were used as molecular markers for phylogenetic analysis of particular isolates. The nucleotide sequences of polymerase chain reaction products of all markers were compared with the available sequences of relevant markers of other Paramecium species from GenBank. Phylogenetic trees based on all molecular markers showed that all the nine strains had a very close relationship with Paramecium primaurelia except for the FT8 strain. FT8 consistently showed its unique position in comparison to all other species in the phylogenetic trees. Available sequences of internal transcribed spacer 1 (ITS1) and ITS2 and some other ciliate sequences from GenBank were used for the construction of secondary models. Two highly conserved helices supported by compensatory base changes among all ciliates of ITS2 secondary structures were found similar to other eukaryotes. Therefore, the most conserved 120 to 180 base pairs regions were identified for their comparative studies. We found that out of the three helices in ITS1 structure, helix B was more conserved in Paramecium species. Despite various substitutions in the primary sequence, it was observed that secondary structures of ITS1 and ITS2 could be helpful in interpreting the phylogenetic relationships both at species as well as at generic level.     

Life cycles,phenology and genetic structure of endangered Megacrania tsudai Shiraki (Phasmatodea: Phasmatidae): Male individuals from a geographic parthenogenesis species

Megacrania tsudai, a peripherally distributed member of Megacrania, requires conservation in Taiwan; it has limited distribution in Taiwan and its eastern offshore islands. It feeds on screw pines (Pandanus odoratissimus) in nature and has demonstrated a specific defensive mechanism involving actinidine secretion from the prothoracic gland. However, details of its distribution area, life cycle and developmental phenology remain largely unknown. In this study, a field survey and review of published works revealed M. tsudai distribution in coastal zones and along river shores near estuaries. At room temperature, the egg period was 128 days. The development of the first to sixth instars required 17, 26, 27, 26, 34 and 43 days, respectively, on average; and a generation cycle required approximately 204 days. The phenology of the mesonotal granules was recorded. Moreover, genetic analysis of the mitochondrial cytochrome oxidase I (COI), 16S ribosomal DNA (16S rDNA) and the nuclear ribosomal spacer indicated the occurrence of genetic drift. Therefore, the rearing procedures proposed in this study for the primary and last instars of M. tsudai can facilitate its conservation. Megacrania tsudai was previously recorded as parthenogenetic; however, two male individuals were fostered unexpectedly. The male body length was 91 mm, which is shorter than the female length (120 mm). During mating, the male climbs onto the female's back and protrudes its genitalia downward. Geographical parthenogenesis is likely the reproductive strategy among peripheral M. tsudai; however, the rarely found M. tsudai male could be an intermediate link of reproductive strategy in the transition from tychoparthenogenesis to parthenogenesis.     

Contrasting mtDNA diversity and population structure in a direct-developing marine gastropod and its trematode parasites

The comparative genetic structure of hosts and their parasites has important implications for their coevolution, but has been investigated in relatively few systems. In this study, we analysed the genetic structure and diversity of the New Zealand intertidal snail Zeacumantus subcarinatus ( n  =   330) and two of its trematode parasites, Maritrema novaezealandensis ( n  =   269) and Philophthalmus sp. ( n  =   246), using cytochrome c oxidase subunit I gene ( COI ) sequences. Snails and trematodes were examined from 11 collection sites representing three regions on the South Island of New Zealand. Zeacumantus subcarinatus displayed low genetic diversity per geographic locality, strong genetic structure following an isolation by distance pattern, and low migration rates at the scale of the study. In contrast, M. novaezealandensis possessed high genetic diversity, genetic homogeneity among collection sites and high migration rates. Genetic diversity and migration rates were typically lower for Philophthalmus sp. compared to M. novaezealandensis and it displayed weak to moderate genetic structure. The observed patterns likely result from the limited dispersal ability of the direct developing snail and the utilization of bird definitive hosts by the trematodes. In addition, snails may occasionally experience long-distance dispersal. Discrepancies between trematode species may result from differences in their effective population sizes and/or life history traits.     

Distribution and ecology of mosquitoes in a focus of dirofilariasis in northwestern Iran, with the first finding of filarial larvae in naturally infected local mosquitoes

A study of the distribution and ecology of mosquitoes (Diptera: Culicidae) in Ardebil Province of northwestern Iran, with special reference to the known dirofilariasis focus in Meshkinshahr County, was carried out in July‐August 2005 and April‐October 2006. In total, 17 533 adult mosquitoes and 3090 third‐ and fourth‐instar larvae were collected on 14 occasions and identified using morphological characters and DNA sequence data. Twenty species belonging to seven genera were found: Anopheles claviger (Meigen), An. hyrcanus (Pallas), An. maculipennis Meigen, An. pseudopictus Grassi*, An. sacharovi Favre, An. superpictus Grassi, Aedes vexans (Meigen)*, Coquillettidia richiardii (Ficalbi)*, Culex hortensis Ficalbi, Cx. modestus Ficalbi, Cx. pipiens Linnaeus, Cx. theileri Theobald, Cx. torrentium Martini*, Cx. tritaeniorhynchus Giles, Culiseta longiareolata (Macquart), Cs. annulata (Schrank)*, Cs. subochrea (Edwards), Ochlerotatus caspius (Pallas) s.l.* (= Aedes caspius sensu auctorum), Oc. geniculatus (Olivier)* (= Aedes geniculatus sensu auctorum) and Uranotaenia unguiculata Edwards (asterisks indicate new occurrence records for the province). The most prevalent species in adult catches were An. maculipennis (52%), Cx. theileri (45%) and Cx. hortensis (1%); the most prevalent species caught as larvae were Cx. theileri (27%), Cx. hortensis (21%) and An. maculipennis (19%). Anopheles maculipennis, Cx. pipiens and Cx. theileri were most widely distributed in the province. The occurrence of Cx. torrentium in Iran is verified based on differential characters of fourth‐instar larvae. Anopheles maculipennis and An. sacharovi of the Maculipennis Group were identified from their diagnostic ITS2 sequences. For the first time, cytochrome c oxidase I (COI) sequences were obtained from Iranian specimens of An. hyrcanus, An. pseudopictus, Cx. theileri and Oc. caspiuss.l. Culex theileri and An. maculipennis were found naturally infected with third‐stage (infective) larvae of Dirofilaria immitis (Leidy) and Setaria labiatopapillosa (Alessandrini) (Spirurida: Onchocercidae), respectively, for the first time in Iran.     

Integrative taxonomy of the Clavus canalicularis species complex (Drilliidae,Conoidea, Gastropoda) with description of four new species

ABSTRACT

The conoidean family Drilliidae Olsson, 1964 is a species-rich lineage of marine gastropods, showing a high degree of diversification in comparison to other families of Conoidea. Despite intensive molecular phylogenetic studies during the last decade that have led to notable rearrangements of conoidean systematics, the genus- and species-level taxonomy of Drilliidae has not thus far been affected and remains entirely based on shell features. In the current study we revisit species delimitation in a morphological cluster of species from the Indo-Pacific referred to as the Clavus canalicularis complex, using an integrative taxonomy approach. The species in the complex possess robust thick-walled shells typically over 15?mm in height with sculpture of prominent rounded nodules located at the whorl’s shoulder, sometimes sharp and squamiform, or producing long spines. We find that in addition to five known species, the complex comprises four new species. These are described as Clavus brianmayi n. sp. (New Caledonia), Clavus davidgilmouri n. sp. (the Philippines), Clavus andreolbrichi n. sp. (Vanuatu and New Ireland) and Clavus kirkhammetti n. sp. (Madagascar). Clavus exasperatus (Reeve, 1843), which was previously considered widely distributed in Indo-Pacific, is shown to be confined to the western Indian Ocean.     

Bottlenecks, drift and differentiation: the population structure and demographic history of sika deer (Cervus nippon) in the Japanese archipelago

We assessed genetic differentiation and diversity in 14 populations of sika deer (Cervus nippon) from Japan and four populations of sika deer introduced to the UK, using nine microsatellite loci. We observed extreme levels of differentiation and significant differences in diversity between populations. Our results do not support morphological subspecies designations, but are consistent with previous mitochondrial DNA analyses which suggest the existence of two genetically distinct lineages of sika deer in Japan. The source of sika introduced to the UK was identified as Kyushu. The underlying structure of Japanese populations probably derives from drift in separate glacial refugia and male dispersal limited by distance. This structure has been perturbed by bottlenecks and habitat fragmentation, resulting from human activity from the mid-nineteenth century. Most current genetic differentiation and differences in diversity among populations probably result from recent drift. Coalescent model analysis suggests sika on each of the main Japanese islands have experienced different recent population histories. Hokkaido, which has large areas of continuous habitat, has maintained high levels of gene flow. In Honshu the population is highly fragmented and is likely to have been evolving by drift alone. In Kyushu there has been a balance between gene flow and drift but all the populations have experienced high levels of drift. Habitat fragment size was not significantly associated with genetic diversity in populations but there was a significant correlation between habitat fragment size and effective population size.     

DNA barcoding of the genus Lepidion (Gadiformes: Moridae) with recognition of Lepidion eques as a junior synonym of Lepidion lepidion

DNA sequences of cytochrome c oxidase I gene (COI) from Lepidion spp. were employed to test the efficiency of species identification. A sample of 32 individuals from five Lepidion species was sequenced and combined with 26 sequences from other BOLD projects. As a result, 58 Lepidion DNA sequences of the COI gene belonging to eight of the nine recognized Lepidion species were analysed. Sequences were aligned and formed seven clades in a Bayesian phylogenetic tree, where Lepidion lepidion and Lepidion eques grouped jointly. The Kimura 2‐parameter genetic distances, among congeners were, on average, 4.28%, 16 times greater than among conspecifics (0.27%). The main diagnostic meristic data of Lepidion spp. were compiled and a detailed morphological revision of the congeneric species L. eques and L. lepidion was made. The eye diameter was significantly different between L. eques and L. lepidion (P < 0.001). The number of anal fin rays ranged from 45 to 51 in L. lepidion and from 47 to 54 in L. eques, but no significant differences were obtained in the mean values of this variable (P = 0.07). According to the morphological and genetic analyses, the results strongly suggest that the Mediterranean codling L. lepidion and the North Atlantic codling L. eques are conspecific, making L. eques a junior synonym of L. lepidion.     

Barcoding,molecular taxonomy,and exploration of the diversity of shrews (Soricomorpha: Soricidae) on Mount Nimba (Guinea)

We tested the efficiency of cytochrome oxidase I (COI)‐barcoding as a taxonomic tool to discriminate and identify sympatric shrew species on Mount Nimba (Guinea). We identified 148 specimens at the species level using morphological characters and comparison with type specimens, including several taxa from Mount Nimba. We identified ten morphospecies and tested aspects of genetic diversity and monophyly using genetic data from three mitochondrial (16S, cytochrome b, and COI) and one nuclear marker (the breast cancer gene, BRCA). Nine morphospecies were validated under the phylogenetic and genetic species concepts, including the recently diverged species Crocidura buettikoferi, Crocidura theresae, and Crocidura grandiceps. Under the same concepts, our analyses revealed the presence of two cryptic species amongst animals identified as Crocidura muricauda. We then tested the efficiency of barcoding thanks to commonly used phenetic methods, with the 148 specimens representing 11 potentially valid species based on morphological and molecular data. We show that COI‐barcoding is a powerful tool for shrew identification and can be used for taxonomic surveys. The comparison of genetic divergence values shows the presence of a barcoding gap (i.e. difference between the highest intraspecific and the lowest interspecific genetic divergence values). Given that only a few COI sequences are available for Afrotropical shrews, our work is an important step forward toward their enrichment. We also tested the efficiency of the three other sequenced markers and found that cytochrome b is as efficient as COI for barcoding shrews. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 166 , 672–687.     

Contrasting demographic history and population structure in Capsella rubella and Capsella grandiflora, two closely related species with different mating systems

Both mating system and population history can have large impacts on genetic diversity and population structure. Here, we use multilocus sequence data to investigate how these factors impact two closely related Brassicaceae species: the selfing Capsella rubella and the outcrossing C. grandiflora. To do this, we have sequenced 16 loci in approximately 70 individuals from 7 populations of each species. Patterns of population structure differ strongly between the two species. In C. grandiflora, we observe an isolation-by-distance pattern and identify three clearly delineated genetic groups. In C. rubella, where we estimate the selfing rate to be 0.90-0.94, the pattern is less clear with some sampling populations forming separate genetic clusters while others are highly mixed. The two species also have divergent histories. Our analysis gives support for a bottleneck approximately 73 kya (20-139 kya) in C. rubella, which most likely represents speciation from C. grandiflora. In C. grandiflora, there is moderate support for the standard neutral model in 2 of 3 genetic clusters, while the third cluster and the total data set show evidence of expansion. It is clear that mating system has an impact on these two species, for example affecting the level of genetic variation and the genetic structure. However, our results also clearly show that a combination of past and present processes, some of which are not affected by mating system, is needed to explain the differences between C. rubella and C. grandiflora.     

Low mtDNA variation and shallow population structure of the Chinese pomfret Pampus chinensis along the China coast

In the present study, population genetic structure and genetic diversity of the Chinese pomfret Pampus chinensis, along the China coast were investigated and compared with that from Indonesia using mitochondrial DNA cytochrome b gene sequences. A total of 28 variable sites (including 18 singleton sites and 10 parsimony information sites) were observed and 23 haplotypes were defined in 330 individuals from 11 localities. The haplotype diversity (H_D) of the populations ranged from 0·540 to 0·828, the nucleotide diversity (π) ranged from 0·081 to 0·295%. Pairwise F_ST statistics showed that significant genetic divergence occurred among populations from different geographical regions. The high dispersal capabilities, geographic segregation and ocean currents may be responsible for the present population genetic structure in this species. In addition, a population expansion event during the late Pleistocene period was inferred. The time of population expansion was estimated to occur about 117 000–169 000 years ago.