De novo assembly and characterization of foot transcriptome and microsatellite marker development for Paphia textile

Population Genetics of Haliotis discus hannai in China Inferred Through EST-SSR Markers

BACKGROUND/OBJECTIVES

The Pacific abalone Haliotis discus hannai originated in cold waters and is an economically important aquaculture shellfish in China. Our goal was to clarify the current status of the genetic structure of Pacific abalone in China.

METHODS

In this study, eighteen polymorphic EST-SSR loci were successfully developed based on the hemolymph transcriptome data of Pacific abalone, and thirteen highly polymorphic EST-SSR loci were selected for the genetic variation analysis of the six populations collected.

RESULTS

The results showed that the average number of observed alleles was 8.0769 (RC)-11.3848 (DQ) in each population. The number of observed alleles in the DQ, NH, and TJ populations was significantly higher than that in the RC population. The cultivated population outside the Changshan Islands has experienced a 22.79% reduction in allele diversity compared to the wild population of DQ. The pairwise Fst values and analysis of molecular variance (AMOVA) revealed significant population differentiation among all populations except DQ and NH populations, with RC and ZZ cultured populations exhibiting the largest population differentiation (Fst = 0.1334). The phylogenetic tree and structural analysis divided the six populations into two groups (group 1: NH, DQ, and ZZ; group 2: DL, TJ, and RC), and there was no relationship between geographical distance and genetic distance.

CONCLUSIONS

These results may reflect the large-scale culture from different populations in China and the exchange of juveniles between hatcheries. Different breeding conditions have led to a higher degree of genetic differentiation between the RC and ZZ populations. This study enables a better understanding of the genetic diversity and structure of current Pacific abalone populations.

Reproductive studies on the carpet clam Paphia textile (Paratapes textilis) (Gmelin 1791) (Family: Veneridae): a guide of aquaculture management along the Egyptian coasts of the Red Sea and Suez Canal

BACKGROUND

Most aquatic biota's reproductive biology and life cycle are essential to the sustainable management and development of coastal ecosystems and aquaculture. The bivalve Paphia textile (Gmelin 1791), also known as Paratapes textilis, has an economic value in Indo-Pacific waters, including the Red Sea and the Suez Canal lakes, the Egyptian coasts. However, P. textile suffers from extensive fishing and exploitation.

AIM

The present work aims to study the Paphia textile's reproductive cycle on the Egyptian coasts of the Red Sea for the first time. It helps to manage and develop the coastal ecosystems and aquaculture.

METHODOLOGY

Samples were collected monthly from two saline lakes in the Suez Gulf from December 2019 to November 2020. As part of the comprehensive research study, sex ratio, condition index, sexuality, histological analysis of gonads, shell size, and gonad index were used to investigate the reproductive cycle.

RESULTS

The results reveal a male-biased sex ratio, possibly due to anthropogenic stressors. The Paphia textile is dioecious. No hermaphrodite cases were observed in the studied specimens. The condition index in winter and spring indicates periods dominated by mature individuals. Five reproductive maturity stages were assigned for both P. textile males and females. Due to the simultaneous development of several developmental stages monthly throughout the sampling year, warm water may be responsible for non-sequential gametogenic cycles. As measured environmental parameters correlate with maturity stages, temperature, salinity, and chlorophyll a play important role in gonad growth. The size at first sexual maturity at which 50% of the Paphia textile population reached maturity ranged from 28.60 to 31.50 mm for females, and between 31.70 and 34.10 mm for males. As the gonad index increases during the ripe stages, this index decreases during the resting, spawning, and spent phases.

CONCLUSIONS

The findings suggest the most suitable temperature for aquaculture spawning is between 20 °C and 30 °C in subtropical waters. Fishing should generally be prohibited at sizes less than 28.60 mm for better management and sustainability of this valuable aquatic resource on the Egyptian coasts of the Red Sea.

Development of Microsatellite Markers Based on Transcriptome Sequencing and Evaluation of Genetic Diversity in Swimming Crab (Portunus trituberculatus)

P. trituberculatus is an economically important mariculture species in China. Evaluating its genetic diversity and population structure can contribute to the exploration of germplasm resources and promote sustainable aquaculture production. In this study, a total of 246,243 SSRs were generated by transcriptome sequencing of P. trituberculatus. Among the examined 254,746 unigenes, 66,331 had more than one SSR. Among the different SSR motif types, dinucleotide repeats (110,758, 44.98%) were the most abundant. In 173 different base repeats, A/T (96.86%), AC/GT (51.46%), and ACC/GGT (26.20%) were dominant in mono-, di-, and trinucleotide, respectively. GO annotations showed 87,079 unigenes in 57 GO terms. Cellular process, cell, and binding were the most abundant terms in biological process, cellular component, and molecular function categories separately. A total of 34,406 annotated unigenes were classified into 26 functional categories according to the functional annotation analysis of KOG, of which “general function prediction only” was the biggest category (6,028 unigenes, 17.52%). KEGG pathway annotations revealed the clustering of 34,715 unigenes into 32 different pathways. Nineteen SSRs were identified as polymorphic and, thus, used to assess the genetic diversity and structure of 240 P. trituberculatus individuals from four populations in the Bohai Sea. Genetic parameter analysis showed a similar level of genetic diversity within wild populations, and the cultured population indicated a reduction in genetic diversity compared with wild populations. The pairwise F_ST values were between 0.001 and 0.04 with an average of 0.0205 (p < 0.05), suggesting a low but significant level of genetic differentiation among the four populations. Structure analysis demonstrated that the four populations were classified into two groups including the cultured group and other populations. The phylogenetic tree and PCA revealed that a vast number of samples were clustered together and that cultivated individuals were distributed more centrally than wild individuals. The findings contribute to the further assessment of germplasm resources and assist to provide valuable SSRs for marker-assisted breeding of P. trituberculatus in the future.

Analysis of Genetic Diversity of Banana Weevils (Cosmopolites sordidus) (Coleoptera: Curculionidae) Using Transcriptome-Derived Simple Sequence Repeat Markers

The banana weevil, Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) is an economically important insect pest of bananas. It causes up to 100% yield losses and substantial lifespan reduction in bananas. Advances in genomics, proteomics, and sequencing technologies have provided powerful pathways to genotyping disastrous pests such as C. sordidus. However, such technologies are often not available to the majority of rural subtropical African banana growers and pest control managers. This study was therefore motivated by the need to create cheap and easily accessible C. sordidus genotyping methods that could be deployed by banana pest control managers to the benefit of C. sordidus control programs in the tropics where such advanced technologies are not readily accessible. We used an in-house C. sordidus transcriptome from the an-ongoing study from which we mined an array of simple sequence repeat (SSR) markers. Of these, six highly polymorphic transcriptome-derived SSR markers were used to successfully genotype within and among banana weevil population genetic diversity of 12 C. sordidus populations collected from four banana-growing agro-ecological zones (AEZs) in Uganda. The developed transcriptome-derived SSR markers can be used by researchers in population genetics for characterization of the C. sordidus and identification of new genes that are linked to traits of particular interest. The significant genetic diversity revealed in C. sordidus provides pertinent information for integrated pest management strategies.

De novo assembly, transcriptome characterization and marker discovery in Indian major carp, Labeo rohita through pyrosequencing

Labeo rohita, one of the Indian major carps, is the most popular culture species in Indian subcontinent due to its consumer preference and delicacy. A selective breeding program for harvest body weight has resulted in an average genetic gain of 17% per generation. Transcriptome resource for this species is scanty. Here, we have characterized the liver and muscle transcriptomes of rohu using Roche 454 GS-FLX next generation sequencing platform. In total, 1.2 million reads were generated, de novo assembly and clustering resulted in 4171 transcripts. Out of these, 4171 had significant blast hit against NCBI nr database, and 2130 transcripts were successfully annotated. In total, 289 SSRs were identified with an identification rate of 5.8%, and dinucleotide repeat motifs were observed to be the most abundant SSRs. Further, 2231 putative SNPs were identified with high confidence. Validation of eight putative SNPs using Sanger sequencing resulted in 100% true SNPs. Significant allelic imbalance of M1, M4 and M5 loci between growth selected and control individual were observed. Furthermore, 13 transcription factors were identified in the present study belonging to six different transcription factor families. The present study demonstrated the utility of RNAseq to develop genomics resources in non-model fish species, and the marker resources developed would support the genetic improvement program of this species.

Molecular characterization of AwSox2 from bivalve Anodonta woodiana: Elucidating its player in the immune response

Sox2 is an embryonal stem cell Ag essential for early embryonic development, tissue homeostasis and immune regulation. In the current study, one complete Sox2 cDNA sequence was cloned from freshwater bivalve Anodonta woodiana and named AwSox2. Histological changes of testis derived from Bisphenol A (BPA) treatment were analyzed by hematoxylin and eosin staining. Expressions of AwSox2 derived from BPA, LPS and polyinosinic:polycytidylic (Poly I:C) challenge were measured by quantitative real-time PCR. The full-length cDNA of AwSox2 contained an open reading frame of 927 nucleotides bearing the typical structural features of Sox2 family. Obvious degeneration, irregular arrangement of spermatids, and clotted dead and intertwined spermatids were observed in BPA-treated groups. Administration of BPA could result in a dose-dependent up-regulation of AwSox2 expression in the male gonadal tissue of A. woodiana. In addition, expression of AwSox2 was significantly induced by LPS and Poly I:C treatment in the hepatopancreas, gill and hemocytes, compared with that of control group. These results indicated that up-regulations of AwSOx2 are closely related to apoptosis of spermatogonial stem cells derived from BPA treatment as well as enhancement of immune defense against LPS and Poly I:C challenge in A. woodiana.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

Transcriptome analysis of juvenile genetically improved farmed tilapia (Oreochromis niloticus) livers by dietary resveratrol supplementation

Here we used RNA-Seq to explore the transcriptomic response and specific involvement of hepatic mRNA of juvenile Oreochromis niloticus (GIFT) as a result of dietary resveratrol supplementation (0.05 g/kg RES). More than 24,513,018 clean reads were reference genome guided assembly into 23,417 unigenes. 12,596 unigenes (29.64%) were annotated to GO database. There were 5, 179 and 1526 genes significantly differentially expressed genes at 15, 30 and 45 d respectively, and 8 KEGG pathways were enriched associated with this immune response. Hyperemia and compressed hepatic sinusoid, fibrosis of liver cell and abnormal hepatic epidermal cell revealed by H&E and SEM analysis respectively. Genes related with cytokine production (il12rb2, scfr), immune system (ig8l, hlfl, cd226, prf1l), autophagy regulation (atg4b), foxo signaling (ccnb2), steroid hormone biosynthesis (cyp3a40), fatty acid metabolism (scd1), metabolism (cacna1b) have been significantly decreased, while genes associated with such pathways above (leap-2, prdx4, mb, homer1, mif, sat1, cytbc1_8) and the pathway of protein processing in endoplasmic reticulum (cne1, tram1) have been significantly increased. These findings suggested RES activated some immune and biological process-related genes to enhance GIFT's innate immunity. It also suggested high concentration addition or long-time administration may bring negative effect in tilapia liver.

Mining of simple sequence repeats (SSRs) loci and development of novel transferability-across EST-SSR markers from de novo transcriptome assembly of Angelica dahurica

Angelica dahurica is a widely grown plant species with multiple uses, especially in the medical field. However, the frequent introduction of A. dahurica to new areas has made it difficult to distinguish between varieties. Simple sequence repeats (SSRs) detected based on transcriptome analyses are very useful for constructing genetic maps and analyzing genetic diversity. They are also relevant for the molecular marker-assisted breeding of A. dahurica. We identified 33,724 genic SSR loci based on transcriptome sequencing data. A total of 114 primer pairs were designed for the SSR loci and were tested for their specificity and diversity. Ten SSR loci in untranslated regions were ultimately selected. Subsequently, 56 A. dahurica ecotypes collected from different regions were analyzed. The SSR loci comprised 2–8 alleles, with a mean of 5.2 alleles per locus. The polymorphic information content value and Shannon’s information index were 0.6274–0.2702 (average of 0.4091) and 1.3040–0.5618 (average of 0.8475), respectively. Thus, the 10 novel SSRs identified in this study were almost in accordance with Harvey-Weinberg equilibrium and will be useful for analyzing A. dahurica genetic relationships. The results of this study confirm the potential value of transcriptome databases for the development of new SSR markers.

De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea

Crassadoma gigantea is an important commercial marine bivalve species in Baja California and Mexico. In this study, we have applied RNA-Seq technology to profile the transcriptome of the C. gigantea for the first time. A total of 80,832,518 raw reads were produced from a Illumina HiSeq4000 platform, and 77,306,198 (95.64%) clean reads were generated after trimming the adaptor sequences. The transcriptome assembled into 158,855 transcripts with an N50 size of 1995 bp and an average size of 1008 bp. A number of DNA repair related genes, such as MSH3, EGF, TGF, IGF, FGF, encoding different groups of growth factors were found in the transcriptome data of C. gigantean. In addition, immune related genes Toll-like receptor (TLR) including TLR1, TLR2, TLR3, TLR4, TLR5, TLR6, TLR7, TLR8, and TLR9 was also observed in C. gigantean. A set of 12 polymorphic microsatellite loci were firstly developed and characterized in C. gigantea. The results show that the number of alleles and expected heterozygosity ranged from 3 to 9 and from 0.254 to 0.820, respectively. The average polymorphic information content was 0.790. These microsatellite loci will facilitate future studies of population structure and conservation genetics in this species.

Hepatic transcriptome analysis of juvenile GIFT tilapia (Oreochromis niloticus), fed diets supplemented with different concentrations of resveratrol

The GIFT (Genetically Improved Farmed Tilapia) tilapia, Oreochromis niloticus, is cultured widely for the production of freshwater fish in China. Streptococcosis, which is related to pathogenic infections, occurs frequently in juvenile and adult female GIFT individuals. Resveratrol (RES) has been used in feed to control these infections in freshwater tilapia. To address the effects of RES on tilapia, we used high-throughput RNA sequencing technology (RNA-Seq, HiSeq. 2500) to explore the global transcriptomic response and specific involvement of hepatic mRNA of juvenile O. niloticus fed with diets containing different concentrations of (0, 0.025, 0.05, and 0.1g/kg) RES. A total of > 24,513,018 clean reads were generated and then assembled into 23,244 unigenes. The unigenes were annotated by comparing them against non-redundant protein sequence (Nr), non-redundant nucleotide (Nt), Swiss-Prot, Pfam, Gene Ontology database (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and 12,578 unigenes were annotated to the GO database. A total of 1444 (0.025g/kg RES), 1526 (0.05g/kg RES), and 3135 (0.1g/kg RES) genes were detected as significant differentially expressed genes (DEGs), when compared with the controls. A total of 6 (0.025 vs 0.05g/kg RES), 19 (0.025 vs 0.1g/kg RES), and 124 (0.05 vs 0.1g/kg RES) genes were detected as significant DEGs. Six genes, including dnah7x1, sox4, fam46a, hsp90a, ddit4, and nmrk2, were associated with an immune response. These findings provide information on the innate immunity of GIFT and might contribute to the development of strategies for the effective management of diseases and long-term sustainability of O. niloticus culture.

Development of Novel Polymorphic EST-SSR Markers in Bailinggu (Pleurotus tuoliensis) for Crossbreeding

Identification of monokaryons and their mating types and discrimination of hybrid offspring are key steps for the crossbreeding of Pleurotus tuoliensis (Bailinggu). However, conventional crossbreeding methods are troublesome and time consuming. Using RNA-seq technology, we developed new expressed sequence tag-simple sequence repeat (EST-SSR) markers for Bailinggu to easily and rapidly identify monokaryons and their mating types, genetic diversity and hybrid offspring. We identified 1110 potential EST-based SSR loci from a newly-sequenced Bailinggu transcriptome and then randomly selected 100 EST-SSRs for further validation. Results showed that 39, 43 and 34 novel EST-SSR markers successfully identified monokaryons from their parent dikaryons, differentiated two different mating types and discriminated F₁ and F₂ hybrid offspring, respectively. Furthermore, a total of 86 alleles were detected in 37 monokaryons using 18 highly informative EST-SSRs. The observed number of alleles per locus ranged from three to seven. Cluster analysis revealed that these monokaryons have a relatively high level of genetic diversity. Transfer rates of the EST-SSRs in the monokaryons of closely-related species Pleurotuseryngii var. ferulae and Pleurotus ostreatus were 72% and 64%, respectively. Therefore, our study provides new SSR markers and an efficient method to enhance the crossbreeding of Bailinggu and closely-related species.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

De novo assembly, characterization and annotation for the transcriptome of Sarcocheilichthys sinensis

The Chinese lake gudgeon Sarcocheilichthys sinensis is a small cyprinid fish with great aquaculture potential both for its edible and ornamental values. Nevertheless, available genomic and transcriptomic information for this fish is extremely deficient. In this study, a normalized cDNA library was constructed using 13 mixed tissues of an adult male S. sinensis, and was sequenced by the Illumina HiSeq2500 platform. De novo assembly was performed using 38,911,511 obtained clean reads, and a total of 147,282 unigenes with an average length of 900 bp were finally achieved. 96.2% of these unigenes were annotated in 9 public databases, and 16 segments of growth-related genes were identified for future studies. In addition, 28,493 unigenes were assigned to 61 subcategories of Gene Ontology (GO), and 10,483 unigenes were assigned to 25 categories of Cluster of Orthologous Group (COG). Moreover, 14,943 unigenes were classified into 225 pathways of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 30,666 microsatellites were detected from 17,627 unigenes with an average distribution density of 1:2405 bp. This transcriptome data set will be valuable for researches on discovery, expression and evolution on genes of interest. Meanwhile, the identified microsatellites would be useful tools for genetic and genomic studies in S. sinensis.