Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD-seq for the small yellow croaker (Larimichthys polyactis)

New perspectives on the genetic structure of dotted gizzard shad (Konosirus punctatus) based on RAD-seq

To maintain, develop and rationally utilize marine organisms, understanding their genetic structure and habitat adaptation pattern is necessary. Konosirus punctatus, which is a commercial fish species inhabiting the Indo-west Pacific Ocean, has shown an obvious annual global capture and aquaculture production decline due to climate changes and human activities. In the present study, restriction-site associated DNA sequencing (RAD-seq) was used to describe its genome-wide single nucleotide polymorphisms panel (SNPs). Among 146 individuals collected at nine locations scattered in China, Korea and Japan, a set of 632,090 SNPs were identified. Population genetic analysis showed that K. punctatus individuals were divided into two significant genetic clusters. Meanwhile, potential genetic differentiation between northern and southern population of K. punctatus was found. Treemix results indicated that gene flow existed among sampling locations of K. punctatus, especially from southern Japan to others. Moreover, candidate genes associated with habitat adaptations of K. punctatus were identified, which are involved in diverse physiological processes of K. punctatus including growth and development (e.g., KIDINS220, PAN3), substance metabolism (e.g., PGM5) and immune response (e.g., VAV3, CCT7, HSPA12B). Our findings may aid in understanding the possible mechanisms for the population genetic structure and local adaptation of K. punctatus, which is beneficial to establish the management and conservation units of K. punctatus, guiding the rational use of resources, with reference significance for a profound understanding of the adaptative mechanisms of other marine organisms to the environment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00216-2.

Mixture effects of imidacloprid and difenconazole on enzymatic activity and gene expression in small yellow croakers (Larimichthys polyactis)

Agrochemicals usually exist as mixtures in aqueous ecosystems and have harmful impacts on the natural environment. Nonetheless, the combined effects and underlying mechanisms of agrochemicals on aqueous organisms remain poorly understood. In the present study, the interactive effects of imidacloprid (IMI) and difenconazole (DIF) on the embryos of small yellow croakers (Larimichthys polyactis) were assessed using various toxicological assays, including acute toxicity, enzymatic activity, and gene expression changes. The results showed that DIF (72-h LC₅₀ value of 0.20 mg L^-1) had higher toxicity than IMI (72-h LC₅₀ value of 12.5 mgL^-1). Simultaneously, combinations of IMI and DIF exerted synergistic acute effects on the embryos of L. polyactis. In addition, the SOD, CAT, GST, and CarE activities were noticeably altered in most single and mixed exposures, relative to the untreated control. The expression of four genes (cyp19a1b, ngln2, klf2a, and socs3a) related to the immune system, endocrine system, and neurodevelopment was also surprisingly altered when the embryos of L. polyactis were subjected to individual and combined exposures relative to the untreated control. Changes in enzymatic activity and gene expression might provide early warning indices for the identification of agrochemical co-exposure. The results of this study provide valuable insights into the comprehensive toxicity of agrochemical mixtures to L. polyactis. Further studies on the long-term effects of agrochemical mixtures on marine fish should be conducted to formulate definitive conclusions concerning hazards.

Joint toxic effects of phoxim and lambda-cyhalothrin on the small yellow croaker (Larimichthys polyactis)

Although pesticides commonly exist as combinations in real-life situations of the aquatic ecosystem, the impact of the toxicity of their mixtures has remained largely unclear. In this study, we investigated the combined effects of two neurotoxic pesticides, including one organophosphate insecticide phoxim (PHO) and one pyrethroid insecticide lambda-cyhalothrin (LCY), on the embryos of the small yellow croaker (Larimichthys polyactis), and their potential pathways. LCY exhibited higher toxicity relative to PHO, with a 72-h LC₅₀ value of 0.0074 mg a.i. L^-1, while the corresponding value for PHO was 0.12 mg a.i. L^-1. The mixture of PHO and LCY exerted a synergistic effect on the embryos of L. polyactis. The activities of antioxidant enzyme CAT and apoptotic enzyme caspase 3 were substantially changed in most single and combined exposure groups relative to the baseline value. Under both single and combined exposures, more significant changes were found in the mRNA expression of five genes, including the immunosuppression gene ngln2, the apoptosis gene P53, the endocrine system gene cyp19a1b, as well as neurodevelopment genes of ap and acp2, relative to the baseline value. Furthermore, the non-target metabolomic analysis demonstrated that hundreds of differential metabolites, including two bile acids (taurodeoxycholic acid and tauroursodeoxycholic acid), were significantly increased in the exposure groups. The bile acids were closely associated with the gut microbiota, and 16S rRNA sequencing results demonstrated dysfunction of the gut microbiota after exposure, especially in the combined exposure group. Our findings indicated that there might be a potential risk connected to the co-occurrence of these two pesticides in aquatic vertebrates. Consequently, future ecological risk assessments should incorporate synergistic mixtures because the current risk assessments do not consider them.

Genetic structure of the small yellow croaker (Larimichthys polyactis) across the Yellow Sea and the East China Sea by microsatellite DNA variation: implications for the division of management units

The small yellow croaker, Larimichthys polyactis, is a commercial fish of the order Perciformes that mainly inhabit estuaries and coastal waters.In recent years, the resources and catch of L. polyactis have undergone huge fluctuations. To detect genetic variations caused by the fluctuation of resources, genetic diversity of L. polyactis in the coastal waters of China were analyzed in this study using microsatellite DNA marker. The results revealed high genetic diversity of this species. The STRUCTURE, DAPC and F ST results all indicated that there was no genetic structure consistent with the distribution pattern. Overall, our main findings are in agreement with previous studies, indicating that L. polyactis showed high genetic diversity and low genetic differentiation. Our results for high genetic connectivity among L. polyactis localities provide insights into the development of management strategies, that is, to manage this species as a single management unit.

The potential function of KIF17 in large yellow croaker (Larimichthys crocea) spermatid remodeling: molecular characterization and expression pattern during spermiogenesis

KIF17, which belongs to the kinesin-2 protein family, plays an indispensable role in mammalian spermiogenesis. However, the role of KIF17 in fish spermatid remodeling during spermiogenesis remains poorly understood. Therefore, we aimed to study the role of KIF17 in spermatid remodeling during Larimichthys crocea (L. crocea) spermiogenesis. The kif17 cDNA sequence, 3247 bp in length, was cloned from L. crocea testis, which consisted of a 347-bp 5'-untranslated region (UTR), 413-bp 3' -UTR, and 2487-bp open reading frame. Bioinformatic analyses revealed that KIF17 obtained from L. crocea (Lc-KIF17) exhibited a high sequence identity compared with those from other teleosts and possessed the structural features of other kinesin-2 proteins. Based on structural similarity, we speculate that the role of Lc-KIF17 may be similar to that of KIF17 in other animals. Lc-kif17 mRNA was diffusely expressed in L. crocea tissues and was highly expressed in the testis, especially at stage IV testicular development. Immunofluorescence analysis revealed that Lc-KIF17 signals colocalized with β-tubulin signals and migrated from the perinuclear cytoplasm to the side of the nucleus where the tail forms during spermiogenesis. These findings revealed that KIF17 may be involved in L. crocea spermiogenesis. In particular, KIF17 may participate in spermatid remodeling by interacting with perinuclear microtubules during L. crocea spermiogenesis. Collectively, this study contributes to an improved understanding of the mechanism underlying L. crocea spermiogenesis and provides a basis for further research on L. crocea reproduction and development.

Characterization and analysis of transcriptome complexity using SMRT-Seq combined with RNA-Seq for a better understanding of Acanthogobius ommaturus in response to temperature stress

Acanthogobius ommaturus is a eurythermic fish, which is widely distributed in coastal, estuarine and bay waters of China, Japan and Korea. Due to the lack of whole genomic information, full-length transcriptome of A. ommaturus was firstly generated by single molecule real-time sequencing (SMRT-seq) in this study. A total of 49,833 full-length non-redundant transcripts (FLNRTs), 2255 alternative splices, 46,856 simple sequence repeats, 5094 long non-coding RNAs and 2708 transcription factors were obtained. Additionally, FLNRTs were used as reference sequences for the following transcriptome analysis of the temperature stress (7 °C, 14 °C, 21 °C (control), 28 °C and 35 °C). GO and KEGG enrichment analysis using GSEA were performed on all genes in 10 response modules which were screened out by WGCNA. Enrichment analysis showed that protein degradation, immune response and energy metabolism play an active role in the temperature stress of A. ommaturus. The differentially expressed hub genes (DEHGs) in response modules were closely related to adhesion, vascular remodeling and disease. The results of this study provided the first systematical full-length transcriptome profile of A. ommaturus and characterized its temperature stress responses, which will serve as the foundation for further exploring the molecular mechanism of the temperature stress in fish.

Exploring the genomic resources and analysing the genetic diversity and population structure of Chinese indigenous rabbit breeds by RAD-seq

Background

Chinese indigenous rabbits have distinct characteristics, such as roughage resistance, stress resistance and environmental adaptability, which are of great significance to the sustainable development of the rabbit industry in China. Therefore, it is necessary to study the genetic diversity and population structure of this species and develop genomic resources.

Results

In this study, we used restriction site-associated DNA sequencing (RAD-seq) to obtain 1,006,496 SNP markers from six Chinese indigenous rabbit breeds and two imported rabbit breeds. Jiuyishan and Fujian Yellow rabbits showed the highest nucleotide diversity (π) and decay of linkage disequilibrium (LD), as well as higher observed heterozygosity (Ho) and expected heterozygosity (He), indicating higher genetic diversity than other rabbits. The inbreeding coefficient (F_IS) of New Zealand rabbits and Belgian rabbits was higher than that of other rabbits. The neighbour-joining (NJ) tree, principal component analysis (PCA), and population structure analysis of autosomes and Y chromosomes showed that Belgian, New Zealand, Wanzai, Sichuan White, and Minxinan Black rabbits clustered separately, and Fujian Yellow, Yunnan Colourful, and Jiuyishan rabbits clustered together. Wanzai rabbits were clearly separated from other populations (K = 3), which was consistent with the population differentiation index (F_ST) analysis. The selection signature analysis was performed in two populations with contrasting coat colours. With Sichuan White and New Zealand rabbits as the reference populations and Minxinan Black and Wanzai rabbits as the target populations, 408, 454, 418, and 518 genes with a selection signature, respectively, were obtained. Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the genes with a selection signature. The results showed that the genes with a selection signature were enriched in the melanogenesis pathway in all four sets of selection signature analyses.

Conclusions

Our study provides the first insights into the genetics and genomics of Chinese indigenous rabbit breeds and serves as a valuable resource for the further effective utilization of the species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07833-6.

A Transient Hermaphroditic Stage in Early Male Gonadal Development in Little Yellow Croaker, Larimichthys polyactis

Animal taxa show remarkable variability in sexual reproduction, where separate sexes, or gonochorism, is thought to have evolved from hermaphroditism for most cases. Hermaphroditism accounts for 5% in animals, and sequential hermaphroditism has been found in teleost. In this study, we characterized a novel form of the transient hermaphroditic stage in little yellow croaker (Larimichthys polyactis) during early gonadal development. The ovary and testis were indistinguishable from 7 to 40 days post-hatching (dph). Morphological and histological examinations revealed an intersex stage of male gonads between 43 and 80 dph, which consist of germ cells, somatic cells, efferent duct, and early primary oocytes (EPOs). These EPOs in testis degenerate completely by 90 dph through apoptosis yet can be rescued by exogenous 17-β-estradiol. Male germ cells enter the mitotic flourishing stage before meiosis is initiated at 180 dph, and they undergo normal spermatogenesis to produce functional sperms. This transient hermaphroditic stage is male-specific, and the ovary development appears to be normal in females. This developmental pattern is not found in the sister species Larimichthys crocea or any other closely related species. Further examinations of serum hormone levels indicate that the absence of 11-ketotestosterone and elevated levels of 17-β-estradiol delineate the male intersex gonad stage, providing mechanistic insights on this unique phenomenon. Our research is the first report on male-specific transient hermaphroditism and will advance the current understanding of fish reproductive biology. This unique gonadal development pattern can serve as a useful model for studying the evolutionary relationship between hermaphroditism and gonochorism, as well as teleost sex determination and differentiation strategies.

Pleistocene isolation caused by sea-level fluctuations shaped genetic characterization of Pampus minor over a large-scale geographical distribution

The southern lesser pomfret (Pampusminor) is an economically important fish, and its numbers are declining because of overfishing and environmental pollution. In addition, owing to the similarities of its external morphological characteristics to other species in the genus Pampus, it is often mistaken for grey pomfret (P.cinereus) or silver pomfret (P.argenteus) juveniles. In this study, the genetic diversity and structure of 264 P.minor individuals from 11 populations in China and Malaysia coastal waters were evaluated for the first time, to the best of our knowledge, using mitochondrial cytochrome b fragments. The results showed that P.minor had moderate haplotype diversity and low nucleotide diversity. Furthermore, two divergent lineages were detected within the populations, but the phylogenetic structure corresponded imperfectly with geographical location; thus, the populations may have diverged in different glacial refugia during the Pleistocene low sea levels. Analysis of molecular variation (AMOVA) showed that genetic variation originated primarily from individuals within the population. Pairwise F_ST results showed significant differentiation between the Chinese and Malaysian populations. Except for the Xiamen population, which was classified as a marginal population, the genetic differentiation among the other Chinese populations was not significant. During the Late Pleistocene, P.minor experienced a population expansion event starting from the South China Sea refugium that expanded outward, and derivative populations quickly occupied and adapted to the new habitat. The results of this study will provide genetic information for the scientific conservation and management of P.minor resources.

Development of Genomic Resources and Identification of Genetic Diversity and Genetic Structure of the Domestic Bactrian Camel in China by RAD Sequencing

The domestic Bactrian camel is indispensable to agricultural production in the desertification area of China owning to its endurance to hunger and thirst, cold resistance, drought resistance, and good long-distance transportation. Therefore, it is necessary to investigate the genetic diversity, genetic structure, and genes with important roles in the evolution of this species. In this study, 1,568,087 SNPs were identified in 47 domestic Bactrian camels inhabiting four regions of China, namely Inner Mongolia, Gansu, Qinghai, and Xinjiang, by restriction site associated DNA sequencing (RAD-seq). The SNP data were used for nucleotide diversity analysis (π) and linkage disequilibrium (LD) attenuation analysis to elucidate the genetic diversity of the domestic Bactrian camel in the four regions studied. Results showed that Xinjiang camels had the highest nucleotide diversity and the fastest decay rate of the LD coefficient; therefore, Xinjiang camels had the highest genetic diversity. Structure analysis, principal component analysis (PCA), and phylogenetic tree construction by the neighbor-joining (NJ) method showed that Qinghai camels clustered separately, at a larger phylogenetic distance from camels in the other regions. Through analyses of selection signals, it was found that the number of selected genes shared by Inner Mongolia camels, Qinghai camels, Xinjiang camels, and Gansu camels was 7, 24, 25, and 113, respectively. The shared selected genes of the domestic Bactrian camel in the four regions were further analyzed, and three shared genes (GRIA3, XIAP, and THOC2) of the domestic Bactrian camel in China were identified. Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the shared selected genes of the domestic Bactrian camel in all four regions studied. Across all regions, genes involved in the cellular process were the most abundant subcategory under biological process. Cell and cell part represented the main proportion of genes under cellular component. Binding represented the main molecular function. In addition, the shared selected genes of the domestic Bactrian camel in the four regions of China were significantly enriched in the long-term depression pathway. The research should enable further study of the genetic resources of the domestic Bactrian camel, as well as the conservation of these resources.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Expression and potential functions of KIF3A/3B to promote nuclear reshaping and tail formation during Larimichthys polyactis spermiogenesis

KIF3A and KIF3B are homologous motor subunits of the Kinesin II protein family. KIF3A, KIF3B, and KAP3 form a heterotrimeric complex and play a significant role in spermatogenesis. Here, we first cloned full-length kif3a/3b cDNAs from Larimichthys polyactis. Lp-kif3a/3b are highly related to their homologs in other animals. The proteins are composed of three domains, an N-terminal head domain, a central stalk domain, and a C-terminus tail domain. Lp-kif3a/3b mRNAs were found to be ubiquitously expressed in the examined tissues, with high expression in the testis. Fluorescence in situ hybridization (FISH) was used to analyze the expression of Lp-kif3a/3b mRNAs during spermiogenesis. The results showed that Lp-kif3a/3b mRNAs had similar expression pattern and were continuously expressed during spermiogenesis. From middle spermatid to mature sperm, Lp-kif3a/3b mRNAs gradually localized to the side of the spermatid where the midpiece and tail form. In addition, we used immunofluorescence (IF) to observe that Lp-KIF3A protein co-localizes with tubulin during spermiogenesis. In early spermatid, Lp-KIF3A protein and microtubule signals were randomly distributed in the cytoplasm. In middle spermatid, however, the protein was detected primarily around the nucleus. In late spermatid, the protein migrated primarily to one side of the nucleus where the tail forms. In mature sperm, Lp-KIF3A and microtubules accumulated in the midpiece. Moreover, Lp-KIF3A co-localized with the mitochondria. In mature sperm, Lp-KIF3A and mitochondria were present in the midpiece. Therefore, Lp-KIF3A/KIF3B may be involved in spermiogenesis in L. polyactis, particularly during nuclear reshaping and tail formation.

Transcriptomic Responses of Two Ecologically Divergent Populations of Japanese Mantis Shrimp (Oratosquilla oratoria) under Thermal Stress

Simple Summary

Rising ocean temperature would change the seawater chemistry and affect the external and internal physiology of crustaceans due to their lack of certain efficient temperature regulators. In addition, the infraspecific populations of crustaceans might also have different response strategies to the rising of temperature. Therefore, we identified the transcriptomic variations to the same thermal stress between ecologically divergent populations of Oratosquilla oratoria. The aim of this study was to investigate the population-specific function genes and relevant pathways in response to thermal stress in O. oratoria. The results showed that gene-expressed variation was in a population-specific pattern, which indicated that the local environment could lead to the evolvement of changes in gene regulation, ultimately leading to adaptive divergences. Additionally, we found several genes with large pleiotropic effects in the Zhoushan population, which might indicate that the regulation mechanisms of the Zhoushan population were more efficient than those of the Qingdao population under same thermal stress. The results provided some novel insights into the local adaptive differences of the infraspecific populations of O. oratoria and other crustaceans.

Abstract

Crustaceans are generally considered more sensitive to ocean warming due to their lack of certain efficient regulators. However, the alterations in the physiology and behavior of crustaceans in response to thermal stress differ vastly even among the infraspecific populations of heterogeneous landscapes. Consequently, understanding the impact of temperature fluctuation on crustacean infraspecific populations might be essential for maintaining a sustainable persistence of populations at existing locations. In the present study, we chose the Japanese mantis shrimp (Oratosquilla oratoria) as the representative crustacean population, and conducted transcriptome analyses in two divergent O. oratoria populations (the Zhoushan and Qingdao populations) under same thermal stress (20–28 °C) to identify the population-specific expression response to thermal stress. The results showed significant differences in gene expressions, GO terms and metabolic pathways between the two populations. We hypothesized that intraspecific mutations in the same or different genes might lead to thermal adaptive divergences. Temperature increases from 20–28 °C produced significant enrichment in GO terms and altered the metabolic pathways in the Zhoushan population despite the lack of differentially expressed unigenes. Therefore, several functional genes with large pleiotropic effects may underlie the response to thermal stress in the Zhoushan population. Furthermore, the most significantly enriched biological processes of the Qingdao population were associated with the state or activity of cells and its significant enriched pathways with genetic information processing as well as immune and environmental information processing. In contrast, the differentially regulated unigenes of the Zhoushan population were primarily involved in the regulatory cellular and transcription processes and the most significant pathways found were metabolic and digestive. Consequently, the regulatory mechanisms of the Zhoushan population are probably more efficient than those of the Qingdao population under the same thermal stress.

Analysis of evolution and genetic diversity of sweetpotato and its related different polyploidy wild species I. trifida using RAD-seq

BACKGROUND

Sweetpotato (Ipomoea batatas (L.) Lam.) is one of the most important crops from the family of Convolvulaceae. It is widely reported that cultivated sweetpotato was originated from Ipomoea trifida. However, diploid, tetraploid and hexaploid I. trifida were found in nature. The relationship, between them, and among them and sweetpotato, is remaining unclear.

RESULTS

In the present study, we detected the genome diversity and relationship of sweetpotato and different polyploidy types I. trifida using Restriction-site Associated DNA Sequencing (RAD-seq). A total of 38,605 RAD-tags containing 832,204 SNPs had been identified. These tags were annotated using five public databases, about 11,519 tags were aligned to functional genes in various pathways. Based on SNP genotype, phylogenetic relation analysis results confirmed that cultivated sweetpotato has a closer relationship with I. trifida 6× than with I. trifida 4X and I. trifida 2×. Besides, 5042 SSRs were detected in I. trifida 6×, and 3202 pairs of high-quality SSR primers were developed. A total of 68 primers were randomly selected and synthesized, of which 61 were successfully amplified.

CONCLUSION

These results provided new evidence that cultivated sweetpotato originated from I. trifida 6×, and that I. trifida 6× evolved from I. trifida 4X and I. trifida 2×. Therefore, using I. trifida 6× as the model plant of sweetpotato research should be more practical than using I. trifida 2× in the future. Meanwhile, sequence information and markers from the present study will be helpful for sweetpotato and I. trifida studies in the future.

Clipperton Atoll as a model to study small marine populations: Endemism and the genomic consequences of small population size

Estimating population sizes and genetic diversity are key factors to understand and predict population dynamics. Marine species have been a difficult challenge in that respect, due to the difficulty in assessing population sizes and the open nature of such populations. Small, isolated islands with endemic species offer an opportunity to groundtruth population size estimates with empirical data and investigate the genetic consequences of such small populations. Here we focus on two endemic species of reef fish, the Clipperton damselfish, Stegastes baldwini, and the Clipperton angelfish, Holacanthus limbaughi, on Clipperton Atoll, tropical eastern Pacific. Visual surveys, performed over almost two decades and four expeditions, and genetic surveys based on genomic RAD sequences, allowed us to estimate kinship and genetic diversity, as well as to compare population size estimates based on visual surveys with effective population sizes based on genetics. We found that genetic and visual estimates of population numbers were remarkably similar. S. baldwini and H. limbaughi had population sizes of approximately 800,000 and 60,000, respectively. Relatively small population sizes resulted in low genetic diversity and the presence of apparent kinship. This study emphasizes the importance of small isolated islands as models to study population dynamics of marine organisms.

Inferring spatial patterns and drivers of population divergence of Neolitsea sericea (Lauraceae), based on molecular phylogeography and landscape genomics

The relative roles of geography, climate and ecology in driving population divergence and (incipient) speciation has so far been largely neglected in studies addressing the evolution of East Asia's island flora. Here, we employed chloroplast and ribosomal DNA sequences and restriction site-associated DNA sequencing (RADseq) loci to investigate the phylogeography and drivers of population divergence of Neolitsea sericea. These data sets support the subdivision of N. sericea populations into the Southern and Northern lineages across the 'Tokara gap'. Two distinct sublineages were further identified for the Northern lineage of N. sericea from the RADseq data. RADseq was also used along with approximate Bayesian computation to show that the current distribution and differentiation of N. sericea populations resulted from a combination of relatively ancient migration and successive vicariant events that likely occurred during the mid to late Pleistocene. Landscape genomic analyses showed that, apart from geographic barriers, barrier, potentially local adaptation to different climatic conditions appears to be one of the major drivers for lineage diversification of N. sericea.

Assessing the genetic diversity of the critically endangered Chinese sturgeon Acipenser sinensis using mitochondrial markers and genome-wide single-nucleotide polymorphisms from RAD-seq

As a living fossil, the endangered Chinese sturgeon (Acipenser sinensis) has been considered a national treasure in China. Here, the famous Gezhouba Dam and Three Gorges Dam on the Yangtze River were built in 1988 and 2006, for economic purposes. The natural population of Chinese sturgeon has declined since then, as these dams block its migratory route to the original spawning grounds in the middle reaches of the Yangtze River. In 2013 and 2014, there was an absence of spawning where it typically happened near the Gezhouba Dam. Nevertheless, from April to June in 2015, over 1,000 larvae with different body lengths (10-35 cm) were detected along the Shanghai Yangtze Estuary; but only little is currently known about the population genetic structure of the Chinese sturgeon. Herein, we inferred population genetic parameters from 462 available Chinese sturgeon specimens based on a 421-bp fragment of the mitochondrial DNA (mtDNA) D-loop region and 1,481,620 SNPs (single-nucleotide polymorphisms) generated by restriction site-associated DNA sequencing (RAD-seq). For the D-loop dataset, 15 haplotypes were determined. Randomly picked 23 individuals, representing the 15 D-loop haplotype groups, were subsequently used for further RAD-seq validation. The average nucleotide diversity calculated from the mtDNA and RAD datasets was 0.0086 and 0.000478, respectively. The overall effective female population size was calculated to be 1,255 to 2,607, and the long-term effective population size was estimated to range from 11,950 to 119,500. We observed that the genetic variability and the effective female population size of the current population in the Yangtze River are severely low, which are similar to the data reported over 10 years ago. The deduced relatively small effective population of female fish, limiting the genetic connectivity among Chinese sturgeon, should be considered a serious threat to this endangered species.

Genome-wide single nucleotide polymorphisms (SNPs) for a model invasive ascidian Botryllus schlosseri

Invasive species cause huge damages to ecology, environment and economy globally. The comprehensive understanding of invasion mechanisms, particularly genetic bases of micro-evolutionary processes responsible for invasion success, is essential for reducing potential damages caused by invasive species. The golden star tunicate, Botryllus schlosseri, has become a model species in invasion biology, mainly owing to its high invasiveness nature and small well-sequenced genome. However, the genome-wide genetic markers have not been well developed in this highly invasive species, thus limiting the comprehensive understanding of genetic mechanisms of invasion success. Using restriction site-associated DNA (RAD) tag sequencing, here we developed a high-quality resource of 14,119 out of 158,821 SNPs for B. schlosseri. These SNPs were relatively evenly distributed at each chromosome. SNP annotations showed that the majority of SNPs (63.20%) were located at intergenic regions, and 21.51% and 14.58% were located at introns and exons, respectively. In addition, the potential use of the developed SNPs for population genomics studies was primarily assessed, such as the estimate of observed heterozygosity (H _O ), expected heterozygosity (H _E ), nucleotide diversity (π), Wright's inbreeding coefficient (F _IS ) and effective population size (Ne). Our developed SNP resource would provide future studies the genome-wide genetic markers for genetic and genomic investigations, such as genetic bases of micro-evolutionary processes responsible for invasion success.

An optimized approach for local de novo assembly of overlapping paired-end RAD reads from multiple individuals

Restriction site-associated DNA (RAD) sequencing is revolutionizing studies in ecological, evolutionary and conservation genomics. However, the assembly of paired-end RAD reads with random-sheared ends is still challenging, especially for non-model species with high genetic variance. Here, we present an efficient optimized approach with a pipeline software, RADassembler, which makes full use of paired-end RAD reads with random-sheared ends from multiple individuals to assemble RAD contigs. RADassembler integrates the algorithms for choosing the optimal number of mismatches within and across individuals at the clustering stage, and then uses a two-step assembly approach at the assembly stage. RADassembler also uses data reduction and parallelization strategies to promote efficiency. Compared to other tools, both the assembly results based on simulation and real RAD datasets demonstrated that RADassembler could always assemble the appropriate number of contigs with high qualities, and more read pairs were properly mapped to the assembled contigs. This approach provides an optimal tool for dealing with the complexity in the assembly of paired-end RAD reads with random-sheared ends for non-model species in ecological, evolutionary and conservation studies. RADassembler is available at https://github.com/lyl8086/RADscripts.

Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing

The KwaZulu‐Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu‐Natal Province, South Africa. In this study, we developed a single‐nucleotide polymorphism (SNP) dataset from double‐digest restriction site‐associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focusing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high‐quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well‐supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north–south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.

Population genetics analysis of the Nujiang catfish Creteuchiloglanis macropterus through a genome-wide single nucleotide polymorphisms resource generated by RAD-seq

Advances in genome scanning using high-throughput sequencing technologies has led to a revolution in studies of non-model organisms. The glyptosternoid fish Creteuchiloglanis macropterus, is widely distributed in the main stem and tributaries of the Nujiang River basin. Here, we analyzed IIB restriction-site-associated DNA (2b-RAD) sequences and mitochondrial DNA sequences, to assess the genomic signature of adaptation by detecting and estimating the degree of genetic differentiation among ten Creteuchiloglanis macropterus populations from the Nujiang River. The analyses revealed significant population differentiation among the up-tributaries, main stem, mid-tributary and low-tributary. Annotation of contigs containing outlier SNPs revealed that the candidate genes showed significant enrichment in several important biological process terms between up-tributaries and low-tributary, and exhibited prominent enrichment in the term macromolecular metabolic process between all tributaries and the main stem. Population dynamics analyses indicated that the Late Pleistocene glaciations strongly influenced the demographic history of C. macropterus. Our results provide strong evidence for the utility of RAD-seq in population genetics studies, and our generated SNP resource should provide a valuable tool for population genomics studies of C. macropterus in the future.

Genomic evidence for local adaptation in the ovoviviparous marine fish Sebastiscus marmoratus with a background of population homogeneity

Advances in next-generation sequencing techniques have allowed for the generation of genome-wide sequence data, to gain insight into the dynamics influencing genetic structure and the local adaptation of marine fish. Here, using genotyping-by-sequencing (GBS) technique, we identified 31,119 single nucleotide polymorphisms (SNPs) for Sebastiscus marmoratus in 59 individuals from three populations in Chinese coastal waters. Based on all SNPs, there was little evidence of genetic differentiation among populations. However, outlier tests revealed 329 SNPs putatively under divergent selection across populations. Structural and phylogenetic topology analyses based on the outliers showed clear genetic differentiation among populations. Gene Ontology (GO) annotation results revealed that most of these outliers are known or hypothesized to be involved in metabolic process. Together with previous work using mitochondrial cytochrome b sequences, the present results further suggest that the population structure is strongly influenced by locally adaptive pressure. Overall, adaptive evolution in a heterogeneous environment plays an important role in inducing genetic differentiation among local populations. This study increases understanding of the factors (including gene flow and local adaptation) promoting and constraining population genetic differentiation in marine organisms.

Genome-Wide SNP Discovery, Genotyping and Their Preliminary Applications for Population Genetic Inference in Spotted Sea Bass (Lateolabrax maculatus)

Next-generation sequencing and the collection of genome-wide single-nucleotide polymorphisms (SNPs) allow identifying fine-scale population genetic structure and genomic regions under selection. The spotted sea bass (Lateolabrax maculatus) is a non-model species of ecological and commercial importance and widely distributed in northwestern Pacific. A total of 22 648 SNPs was discovered across the genome of L. maculatus by paired-end sequencing of restriction-site associated DNA (RAD-PE) for 30 individuals from two populations. The nucleotide diversity (π) for each population was 0.0028±0.0001 in Dandong and 0.0018±0.0001 in Beihai, respectively. Shallow but significant genetic differentiation was detected between the two populations analyzed by using both the whole data set (F_ST = 0.0550, P < 0.001) and the putatively neutral SNPs (F_ST = 0.0347, P < 0.001). However, the two populations were highly differentiated based on the putatively adaptive SNPs (F_ST = 0.6929, P < 0.001). Moreover, a total of 356 SNPs representing 298 unique loci were detected as outliers putatively under divergent selection by F_ST-based outlier tests as implemented in BAYESCAN and LOSITAN. Functional annotation of the contigs containing putatively adaptive SNPs yielded hits for 22 of 55 (40%) significant BLASTX matches. Candidate genes for local selection constituted a wide array of functions, including binding, catalytic and metabolic activities, etc. The analyses with the SNPs developed in the present study highlighted the importance of genome-wide genetic variation for inference of population structure and local adaptation in L. maculatus.