Characterization of hybridization within a secondary contact region of the inshore fish, Bostrychus sinensis, in the East China Sea

Phylogeography and population structure of Lagocephalus spadiceus (Richardson, 1845) (Tetraodontiformes, Tetraodontidae) in the South China Sea

The climate fluctuations during the Late Pleistocene significantly influenced the phylogeographic structure and historical dynamics of marine fishes in the marginal seas of the western Pacific Ocean. The puffer fish, Lagocephalus spadiceus, holds substantial nutritional and economic value in the South China Sea. To investigate the demographic history and population structure of the L. spadiceus, the mitochondrial DNA COI and Cyt b gene datasets from 300 individuals across eight populations in the South China Sea were sequenced. Our findings revealed high haplotype diversity (0.874 ± 0.013) and low nucleotide diversity (0.00075 ± 0.00058). The phylogenetic tree and haplotype networks revealed no significant genetic differentiation along the northern coast of South China Sea. Neutrality tests, mismatch distribution analyses, and Bayesian skyline plots suggested that L. spadiceus underwent population expansion during the Late Pleistocene. Both ocean currents and climate change significantly influenced the geographical distribution and genetic population structure of L. spadiceus.

Mitochondrial genomes reveal mid-Pleistocene population divergence, and post-glacial expansion, in Australasian snapper (Chrysophrys auratus)

Glacial cycles play important roles in determining the phylogeographic structure of terrestrial species, however, relatively little is known about their impacts on the distribution of marine biota. This study utilised modern (n = 350) and ancient (n = 26) mitochondrial genomes from Australasian snapper (Chrysophrys auratus) sampled in New Zealand to assess their demographic and phylogeographic history. We also tested for changes in genetic diversity using the up to 750-year-old mitochondrial genomes from pre-European archaeological sites to assess the potential impacts of human exploitation. Nucleotide diversity and haplotype diversity was high (π = 0.005, h = 0.972). There was no significant change in nucleotide diversity over the last 750 years (p = 0.343), with no detectable loss of diversity as a result of indigenous and industrial-scale fishing activity. While there was no evidence for contemporary population structure (AMOVA, p = 0.764), phylogeographic analyses identified two distinct mitochondrial clades that diverged approximately 650,000 years ago during the mid-Pleistocene, suggesting the species experienced barriers to gene flow when sea levels dropped over 120 m during previous glacial maxima. An exponential population increase was also observed around 8000 years ago consistent with a post-glacial expansion, which was likely facilitated by increased ocean temperatures and rising sea levels. This study demonstrates that glacial cycles likely played an important role in the demographic history of C. auratus and adds to our growing understanding of how dynamic climatic changes have influenced the evolution of coastal marine species.

Genetic pattern and demographic history of cutlassfish (Trichiurus nanhaiensis) in South China Sea by the influence of Pleistocene climatic oscillations

Trichiurus nanhaiensis is one of the most important commercial fish species in the South China Sea. This study aimed to investigate the level of genetic variation and population genetic structure of T. nanhaiensis in the South China Sea for the first time, using 281 individuals collected from seven locations along the coast of mainland China, Taiwan, and Hainan Island. A high level of haplotype diversity and low nucleotide diversity were detected in the mitochondrial DNA cyt b gene and nuDNA RYR 3 gene. The overall expected heterozygosity (He = 0.693) among the seven populations ranged from 0.681 to 0.706 in microsatellite DNA data, which revealed high levels of genetic diversity. Significant genetic differentiation was found in Taidong populations in Taiwan, revealing the prevention of gene flow caused by the Kuroshio Current. Two major lineages based on the cyt b gene suggested that the Taiwan Strait acted as a geographic barrier for T. nanhaiensis during the glacier periods in the late Pleistocene. The Bayesian skyline plot also revealed that population demographic expansion of T. nanhaiensis was estimated to have occurred in 0.1 Mya. Our results indicated that all populations of T. nanhaiensis had experienced a recent genetic bottleneck following recent expansion based on ABC analysis.

Demographic history and population genetic structure of Anisakis pegreffii in the cutlassfish Trichiurus japonicus along the coast of mainland China and Taiwan

Studying the genetic diversity of nematode parasite populations is crucial to gaining insight into parasite infection dynamics and informing parasite phylogeography. Anisakiasis is a zoonotic disease caused by the consumption of infectious third-stage larvae (L3) of Anisakis spp. carried by marine fish. In the present study, a total of 206 mitochondrial DNA sequences (cytochrome c oxidase 2, cox2) were used to study the genetic diversity, genetic structure, and historical demography of twelve A. pegreffii populations from Trichiurus japonicas along the coast of mainland China and Taiwan. Two distinct evolutionary lineages of A. pegreffii and no significant genealogical structures corresponding to sampling localities suggested that isolation in the marginal seas shaped their patterns of phylogeographic distribution along the coast of mainland China and Taiwan during glaciation with lower sea levels. Furthermore, pairwise F_ST values and AMOVA did not indicate any significant genetic differentiation among groups with no relation to the geographic area, which might be attributed to fewer barriers to gene flow as well as large population sizes. The results of the neutrality test, mismatch distribution, and Bayesian skyline plot analyses showed that entire population underwent population expansion during the late Pleistocene. Analysis of the demographic history revealed that A. pegreffii underwent historical lineage diversification and admixture due to secondary contact based on ABC analysis. The present research represents the first definitive population structure and demographic history across sampling locations of A. pegreffii along the coast of mainland China and Taiwan.

Reproductive biology and annual reproductive cycles of two sympatric lineages of Bostrychus sinensis with a natural habitat on southeastern coast of China

The two phylogeographic lineages [South China Sea (SCS) lineage and East China Sea (ESC) lineage] of the four-eyed sleeper (Bostrychus sinensis) occur sympatrically along the southeastern coast of China, where there is a small percentage of hybrids in a natural habitat. To assess the mechanism responsible for the incomplete reproductive isolation between the SCS and ECS lineages of four-eyed sleeper, there was sampling of individuals from the Yueqing Bay from November 2018 to November 2019 to determine whether there are differences in characteristics of the reproductive cycles and reproductive biology. The two lineages varied in reproductive seasonality, body size and shape, and egg number. The SCS lineage spawned annually between May and August, with the absolute fecundity (AF) of the SCS lineage being 9960 to 39,517 eggs per female, and there was positive allometric growth of this species. The ECS lineage spawned annually between June and July, with females producing 4064 to 10,370 eggs per fish and individuals having an isometric growth pattern. Preliminary results indicated that there was partial overlap in the spawning season between the two lineages in the sympatric region, which may not be a primary factor in the incomplete reproductive isolation. Furthermore, the lineage differences in fecundity and body shape and size may lead to variation in mate choice and reproductive strategy. The results from this study provide insight into the biology, reproductive strategies, and speciation of the four-eyed sleepers.

Balancing selection maintains ancient polymorphisms at conserved enhancers for the olfactory receptor genes of a Chinese marine fish

The study of balancing selection, as a selective force maintaining adaptive genetic variation in gene pools longer than expected by drift, is currently experiencing renewed interest due to the increased availability of new data, methods of analysis, and case studies. In this investigation, evidence of balancing selection operating on conserved enhancers of the olfactory receptor (OR) genes is presented for the Chinese sleeper (Bostrychus sinensis), a coastal marine fish that is emerging as a model species for evolutionary studies in the Northwest Pacific marginal seas. Coupled with tests for Gene Ontology enrichment and transcription factor binding, population genomic data allow for the identification of an OR cluster in the sleeper with a downstream flanking region containing three enhancers that are conserved with human and other fish species. Phylogenetic and population genetic analyses indicate that the enhancers are under balancing selection as evidenced by their translineage polymorphisms, excess common alleles, and increased within-group diversities. Age comparisons between the translineage polymorphisms and most recent common ancestors of neutral genealogies substantiate that the former are old, and thus, due to ancient balancing selection. The survival and reproduction of vertebrates depend on their sense of smell, and thereby, on their ORs. In addition to locus duplication and allelic variation of structural genes, this study highlights a third mechanism by which receptor diversity can be achieved for detecting and responding to the huge variety of environmental odorants (i.e., by balancing selection acting on OR gene expression through their enhancer variability).

Nuclear microsatellites reveal population genetic structuring and fine-scale pattern of hybridization in the Japanese mantis shrimp Oratosquilla oratoria

The interplay between historical and contemporary processes can produce complex patterns of genetic differentiation in the marine realm. Recent mitochondrial and nuclear sequence analyses revealed cryptic speciation in the Japanese mantis shrimp Oratosquilla oratoria. Herein, we applied nuclear microsatellite markers to examine patterns and causes of genetic differentiation in this morphotaxon. Population structure analyses revealed two genetically divergent and geographically structured clades in O. oratoria, one dominating the temperate zone of the Northwestern (NW) Pacific and the other occurring in the subtropical and tropical waters where are influenced by the Kuroshio Current. Two sympatric zones, one around the Changjiang Estuary in China coast and the other in the northern Japan Sea, were demonstrated to be hybrid zones where introgressive hybridization occurred asymmetrically. The interaction between historical climate shifts and contemporary factors (e.g., freshwater discharge, temperature gradient and isolation by distance) may contribute to the present-day genetic architecture in the Japanese mantis shrimp. Range shift induced by climate changes and oceanographic factors may promote hybridization and gene flow between the O. oratoria complex. Our results provide insights into the interacting mechanisms that give rise to diversification and speciation of coastal species in the NW Pacific.

Fine-Scale Population Genetic Structure and Parapatric Cryptic Species of Kuruma Shrimp (Marsupenaeus japonicus), Along the Northwestern Pacific Coast of China

The kuruma shrimp (Marsupenaeus japonicus) includes two cryptic species, which are distributed mostly allopatrically but co-occur in the northern South China Sea (from Huilai to Beihai). To obtain a better understanding of the fine-scale genetic structure and parapatric diversification of these two varieties in the northwestern Pacific region, we used a genotyping-by-sequencing (GBS) and comparative transcriptomics approach to establish their phylogenetic relationships. Using the GBS technique, we genotyped 28891 SNPs in 160 individuals in the Northwest Pacific. The results supported two highly diverged evolutionary lineages of kuruma shrimp (var. I and II). The ND and XM populations showed complex genetic patterns, which might be affected by the complex environment of the Taiwan Strait. In addition, the migration rates and inbreeding coefficients of XM and BH were much lower than those of the other populations, which might be related to the land-sea changes and complex ocean currents in the Taiwan Strait and Qiongzhou Strait. Based on the synonymous substitution rates (ds) of 2,491 candidate orthologs, we estimated that the divergence time between the two varieties was 0.26~0.69 Mya. Choice and no-choice interbreeding experiments provided support for the biological species concept, by showing the existence of reproductive isolation or incompatibility. In view of these differences between the two Marsupenaeus species, we believe that it is essential and urgent to establish a genetic database for each and reevaluate their ecological suitable conditions in order to improve species-specific culturing techniques. Moreover, this research can serve as a case study for future research on speciation and hybridization.

Evolutionary adaptation analysis of immune defense and hypoxia tolerance in two closely related Marsupenaeus species based on comparative transcriptomics

Kuruma shrimp, a major farmed shrimp species in the world, includes two cryptic or sibling species, Form I (Marsupenaeus japonicus) and Form II (Marsupenaeus pulchricaudatus). Due to the lack of genomic resources, little is known about the molecular mechanisms associated with immune defense and hypoxia tolerance. Here, we sequenced the transcriptomes of two closely related Marsupenaeus species and compared genomic divergence. This study obtained 77049 and 84561 unigenes with N50 values of 1281bp and 1244bp for M. japonicus and M. pulchricaudatus, respectively, and 5036 pairs of putative orthologs were identified between two Marsupenaeus species. Estimation of Ka/Ks ratios indicated that 165 orthologous genes may be under positive selection (Ka/Ks > 0.5), including 49 pairs with a Ka/Ks ratio >1. According to the peak of synonymous rates, the divergence time between M. japonicus and M. pulchricaudatus was about 0.26-0.69 Mya. These positively selected orthologous genes related to the immune process mainly comprised single VWC domain protein, legumain, ras-related C3 botulinum, caspase, C-type lectin and were enriched in functions related to immune (Toll-like receptor and PI3K-Akt signaling) and hypoxia signaling (HIF-1 signaling and VEGF signaling). In this study, dozens of caspase-like unigenes were screened from two Marsupenaeus transcriptomes. Among these, the PjCaspase orthologous gene was subjected to positive selection (Ka/Ks = 1.22), which had different secondary and three-dimensional structure prediction. Based on the single copy caspase gene, eight populations of Marsupenaeus species were divided into two phylogeographic lineages from the East and South China. We characterized the transcriptomes of the two Marsupenaeus species and obtained several key orthologs associated with immune defense and hypoxia tolerance, which provides new insights into the immunity and genetic divergence of the two varieties. Moreover, this study will facilitate further comparative genomic studies of the two varieties.

Genetic, phenotypic and ecological differentiation suggests incipient speciation in two Charadrius plovers along the Chinese coast

Background

Speciation with gene flow is an alternative to the nascence of new taxa in strict allopatric separation. Indeed, many taxa have parapatric distributions at present. It is often unclear if these are secondary contacts, e.g. caused by past glaciation cycles or the manifestation of speciation with gene flow, which hampers our understanding of how different forces drive diversification. Here we studied genetic, phenotypic and ecological aspects of divergence in a pair of incipient shorebird species, the Kentish (Charadrius alexandrinus) and the White-faced Plovers (C. dealbatus), shorebirds with parapatric breeding ranges along the Chinese coast. We assessed divergence based on molecular markers with different modes of inheritance and quantified phenotypic and ecological divergence in aspects of morphometric, dietary and climatic niches.

Results

Our integrative analyses revealed small to moderate levels of genetic and phenotypic distinctiveness with symmetric gene flow across the contact area at the Chinese coast. The two species diverged approximately half a million years ago in dynamic isolation with secondary contact occurring due to cycling sea level changes between the Eastern and Southern China Sea in the mid-late Pleistocene. We found evidence of character displacement and ecological niche differentiation between the two species, invoking the role of selection in facilitating divergence despite gene flow.

Conclusion

These findings imply that ecology can indeed counter gene flow through divergent selection and thus contributes to incipient speciation in these plovers. Furthermore, our study highlights the importance of using integrative datasets to reveal the evolutionary history and assist the inference of mechanisms of speciation.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1449-5) contains supplementary material, which is available to authorized users.