Genetic Characterization of the Fish Piaractus brachypomus by Microsatellites Derived from Transcriptome Sequencing

Toxic effects of aniline in liver, gills and kidney of freshwater fish Channa punctatus after acute exposure

Aniline (C6H5NH2) is one of the hazardous aromatic amine where an amino group -NH2) is connected to phenyl ring (C6H5). Based on the evaluation of the 96-hour LC50 of aniline, two sublethal concentrations (4.19 mg/l and 8.39 mg/l) were selected for acute exposure tests in freshwater fish Channa punctatus. The liver, gills and kidney of fish being the principal sites of xenobiotic material accumulation, respiration, biotransformation, and excretion are the focus of the present study. Throughout the exposure time, the comet assay revealed increased tail length and tail DNA percentage indicating maximum damage to liver, gills and kidney of treated group after 96 h. After acute exposure, there was a significant (p ≤ 0.05) increase in the enzymatic activity of glutathione-S-transferase (GST) and acetylcholinesterase (AChE), whereas decline in superoxide dismutase (SOD) and catalase (CAT) activity was observed. Meanwhile, levels of malondialdehyde (MDA) increased over the exposure period for both concentrations. After 96 h of exposure, degree of tissue change (DTC) was evaluated in liver, gill and kidney of aniline exposed fish. Additionally, light microscopy revealed multiple abnormalities in liver, gills and kidney of all the treated groups. Significant changes were observed in the levels of biochemical markers viz., glucose, triglyceride, cholesterol, aspartate transaminase, alanine transaminase and urea following a 96-hour exposure to aniline. Studies using ATR-FTIR and transmission electron microscopy (TEM) revealed changes in biomolecules and structural abnormalities in several tissues of the aniline-exposed groups in comparison to the control group respectively.

Microsatellite markers reveal genetic diversity and population structure of Portunus trituberculatus in the Bohai Sea, China

The swimming crab, Portunus trituberculatus, is one of the main aquaculture species in Chinese coastal regions due to its palatability and high economic value. To obtain a better understanding of the genetic diversity of P. trituberculatus in the Bohai Sea, the present study used 40 SSR loci to investigate the genetic diversity and population structure of 420 P. trituberculatus individuals collected from seven populations in the Bohai Sea. Genetic parameters revealed a low level of genetic diversity in the cultured population (SI = 1.374, He = 0.687, and PIC = 0.643) in comparison with wild populations (SI ≥ 1.399, He ≥ 0.692, and PIC ≥ 0.651). The genetic differentiation index (Fst) and gene flow (Nm) ranged from 0.001 to 0.060 (mean: 0.022) and 3.917 to 249.750 (mean: 31.289) respectively, showing a low differentiation among the seven populations of P. trituberculatus. Population structure analysis, phylogenetic tree, and principal component analysis (PCA) demonstrated that the seven groups of P. trituberculatus were divided into four subpopulations (K = 4), but the correlation between genetic structure and geographical distribution was not obvious. These results are expected to provide useful information for the fishery management of wild swimming crabs.

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.

Development of a multi-species SNP array for serrasalmid fish Colossoma macropomum and Piaractus mesopotamicus

Scarce genomic resources have limited the development of breeding programs for serrasalmid fish Colossoma macropomum (tambaqui) and Piaractus mesopotamicus (pacu), the key native freshwater fish species produced in South America. The main objectives of this study were to design a dense SNP array for this fish group and to validate its performance on farmed populations from several locations in South America. Using multiple approaches based on different populations of tambaqui and pacu, a final list of 29,575 and 29,612 putative SNPs was selected, respectively, to print an Axiom AFFYMETRIX (THERMOFISHER) SerraSNP array. After validation, 74.17% (n = 21,963) and 71.25% (n = 21,072) of SNPs were classified as polymorphic variants in pacu and tambaqui, respectively. Most of the SNPs segregated within each population ranging from 14,199 to 19,856 in pacu; and from 15,075 to 20,380 in tambaqui. Our results indicate high levels of genetic diversity and clustered samples according to their hatchery origin. The developed SerraSNP array represents a valuable genomic tool approaching in-depth genetic studies for these species.

Haplotypes traceability and genetic variability of the breeding population of pacu (Piaractus mesopotamicus) revealed by mitochondrial DNA

The main objective of this study was to estimate the genetic diversity levels and haplotype traceability in pacu Piaractus mesopotamicus from the breeding program located in Brazil by analyses of the mitochondrial DNA control region (mtDNA). Moreover, broodstocks from eight commercial fish farms were used for comparative evaluation, four from Brazil (Br1-Br4) and four from Argentina (Ar1-Ar4). The descriptive results revealed 47 polymorphic sites and 51 mutations, which evidenced 34 haplotypes. Ten haplotypes were shared among fish farms and 24 were exclusive. The nucleotide diversity (π) ranged from 0.00031 to 0.01462 and haplotype diversity (Hd) from 0.125 to 0.868. The analysis of molecular variance (AMOVA) indicated high structure present in the analyzed stocks (FST = 0.13356 and ФST = 0.52707). The genetic diversity was high in most of the commercial broodstocks, especially those from Brazil. We observed seven haplotypes in the genetic breeding population, of which four were exclusive and three shared among the commercial fish farms. The genetic diversity was moderate (π = 0.00265 and Hd = 0.424) and considered appropriated for this breeding population of pacu. Our results provide support for the genetic diversity maintenance and mtDNA traceability of pacu commercial broodstocks.

Applications of Next-Generation Sequencing Technologies and Computational Tools in Molecular Evolution and Aquatic Animals Conservation Studies: A Short Review

Aquatic ecosystems that form major biodiversity hotspots are critically threatened due to environmental and anthropogenic stressors. We believe that, in this genomic era, computational methods can be applied to promote aquatic biodiversity conservation by addressing questions related to the evolutionary history of aquatic organisms at the molecular level. However, huge amounts of genomics data generated can only be discerned through the use of bioinformatics. Here, we examine the applications of next-generation sequencing technologies and bioinformatics tools to study the molecular evolution of aquatic animals and discuss the current challenges and future perspectives of using bioinformatics toward aquatic animal conservation efforts.

Assessing Genetic Diversity for a Pre-Breeding Program in Piaractus mesopotamicus by SNPs and SSRs

The pacu (Piaractus mesopotamicus) is a Neotropical fish with remarkable productive performance for aquaculture. Knowledge of genetic resources in Neotropical fish is essential for their applications in breeding programs. The aim of this study was to characterize the genetic diversity of seven farmed populations of pacu which will constitute the basis for a broodstock foundation for coming breeding programs in Brazil. Analysis of one wild population (Paraná River) was used as a reference to compare genetic parameters in the farmed populations. The analyses were performed using 32 single-nucleotide polymorphisms (SNP) and 8 simple sequence repeat (SSR) markers. No significant differences in genetic diversity between populations estimated through the number of alleles and allelic richness, observed heterozygosity, expected heterozygosity, and minimum allele frequency were detected (p > 0.05). Low genetic diversity was observed in all farmed stocks and the wild population. Moreover, we detected low genetic structure when comparing farmed and wild populations for SNPs (F_ST = 0.07; K = 3) and SSRs (F_ST = 0.08; K = 2). Analysis of molecular variance (AMOVA) demonstrated that genetic variation was mostly within populations. Kinship analysis showed that most fish farms included related individuals at a proportion of at least 25%. Our results suggest that the basal broodstock for pacu breeding programs should be founded with individuals from different fish farms for higher genetic diversity and to avoid inbreeding risks.

Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing

Background

The C. macropomum is a characiform fish from the Amazon basin that has been hybridized with other pacu species to produce commercial hybrids, such as the tambacu. However, little is known of the functional genomics of the parental species or these hybrid forms. The transcriptome of C. macropomum and tambacu were sequenced using 454 Roche platform (pyrosequencing) techniques to characterize the domains of Gene Ontology (GO) and to evaluate the levels of gene expression in the two organisms.

Results

The 8,188,945 reads were assembled into 400,845 contigs. A total of 58,322 contigs were annotated with a predominance of biological processes for both organisms, as determined by Gene Ontology (GO). Similar numbers of metabolic pathways were identified in both the C. macropomum and the tambacu, with the metabolism category presenting the largest number of transcripts. The BUSCO analysis indicated that our assembly was more than 40% complete. We identified 21,986 genes for the two fishes. The P and Log2FC values indicated significant differences in the levels of gene expression, with a total of 600 up-regulated genes.

Conclusion

In spite of the lack of a reference genome, the functional annotation was successful, and confirmed a considerable difference in the specificity and levels of gene expression between the two organisms. This report provides a comprehensive baseline for the genetic management of these commercially important fishes, in particular for the identification of specific genes that may represent markers involved in the immunity, growth, and fertility of these organisms, with potential practical applications in aquaculture management.