SNP markers for the genetic characterization of Mexican shrimp broodstocks

A high-density linkage map and sex-determination loci in Pacific white shrimp (Litopenaeus vannamei)

BACKGROUND

Expansion of genomic resources for the Pacific white shrimp (Litopenaeus vannamei), such as the construction of dense genetic linkage maps, is crucial for the application of genomic tools in order to improve economically relevant traits. Sexual dimorphism exists in Pacific white shrimp, and the mapping of the sex-determination region in this species may help in future reproductive applications. We have constructed male, female, and sex-averaged high-density genetic maps using a 50 K single-nucleotide polymorphism (SNP) array, followed by a genome-wide association study (GWAS) to identify genomic regions associated with sex in white shrimp.

RESULTS

The genetic map yielded 15,256 SNPs assigned to 44 linkage groups (LG). The lengths of the male, female, and sex-averaged maps were 5,741.36, 5,461.20 and 5,525.26 cM, respectively. LG18 was found to be the largest for both sexes, whereas LG44 was the shortest for males and LG31 for females. A sex-determining region was found in LG31 with 21 statistically significant SNPs. The most important SNP was previously identified as a sex-linked marker and was able to identify 99% of the males and 88% of the females. Although other significant markers had a lower ability to determine sex, putative genes were intercepted or close to them. The oplophorus-luciferin 2-monooxygenase, serine/arginine repetitive matrix protein and spermine oxidase genes were identified as candidates with possible participation in important processes of sexual differentiation in shrimp.

CONCLUSIONS

Our results provide novel genomic resources for shrimp, including a high-density linkage map and new insights into the sex-determining region in L. vannamei, which may be usefulfor future genetics and reproduction applications.

Genetic diversity of the common carp black strain population based on mtDNA (D-loop and cytb)

The common strain black carp (Cyprinus carpio var. baisenensis) is a culturally important carp strain that is raised and cultured in Guangxi Province, China. Its color reflects the interactions between the Burau people and their surrounding environment. The population of the common carp black strain was isolated and cultured in a rice-fish integration system. To explore the genetic diversity and protection of germplasm resources, we analyzed mitochondrial DNA (mtDNA) sequences, specifically the displacement loop (D-loop) and cytochrome b (Cytb), using single-nucleotide polymorphisms (SNP). We compared these sequences with those from four other local common carp populations. The study included a total of 136 adult common carps from five strain populations: the common black carp strain (HJ = 31), Jian (F = 30), Heilongjiang (H = 10), Songpu (S = 31), and Saijiang (SJ = 34). The results of the Cytb and D-loop analyses showed that the Heilongjiang carp (H) and Saijiang (SJ) populations had the highest levels of haplotype diversity (0.867 ± 0.034785) and nucleotide diversity (π = 0.0063 ± 0.000137 and 0.0093 ± 0.000411), respectively. On the other hand, the Common carp black strain population (HJ) exhibited the lowest haplotype diversity in both Cytb and D-loop, with haplotype 2 being the most commonly observed among the populations. Private haplotypes dominated the five common carp populations, which were significantly different at P<0.001. Furthermore, analyzing the coefficient of genetic differentiation (Fst), the highest genetic difference was observed between Saijiang (SJ) and Heilongjiang (H) (Fst = 0.963), whereas the lowest was observed between Songpu (S) and the Common carp black strain population (HJ) (Fst = 0.019) for the Cytb gene sequences. For the D-loop, the Common carp black strain population (HJ) and Songpu (S) (Fst = 0.7) had the highest values, and Heilongjiang (H) and Common black carp strain (HJ) had an Fst of 0.125. Additionally, the AMOVA analysis revealed a higher level of variance for the Cytb and D-loop genes, indicating lower genetic diversity within the local carp community. On the other hand, the phylogenetic tree analysis showed that the five carp populations were closely related and formed a distinct cluster. The distinct cluster of populations suggests a common ancestor or recent gene flow, possibly due to geographic proximity or migration history, and unique genetic characteristics, possibly due to adaptations or selective pressures. The results of this study provide valuable insights into the genetic diversity of the common strain black carp, which can have implications for conservation, breeding programs, evolutionary studies, and fisheries management.

Improved genome assembly of the whiteleg shrimp Penaeus (Litopenaeus) vannamei using long- and short-read sequences from public databases

The Pacific whiteleg shrimp Penaeus (Litopenaeus) vannamei is a highly relevant species for the world's aquaculture development, for which an incomplete genome is available in public databases. In this work, PacBio long-reads from 14 publicly available genomic libraries (131.2 Gb) were mined to improve the reference genome assembly. The libraries were assembled, polished using Illumina short-reads, and scaffolded with P. vannamei, Feneropenaeus chinensis, and Penaeus monodon genomes. The reference-guided assembly, organized into 44 pseudo-chromosomes and 15,682 scaffolds, showed an improvement from previous reference genomes with a genome size of 2.055 Gb, N50 of 40.14 Mb, L50 of 21, and the longest scaffold of 65.79 Mb. Most orthologous genes (92.6%) of the Arthropoda_odb10 database were detected as "complete," and BRAKER predicted 21,816 gene models; from these, we detected 1,814 single-copy orthologues conserved across the genomic references for Marsupenaeus japonicus, F. chinensis, and P. monodon. Transcriptomic-assembly data aligned in more than 99% to the new reference-guided assembly. The collinearity analysis of the assembled pseudo-chromosomes against the P. vannamei and P. monodon reference genomes showed high conservation in different sets of pseudo-chromosomes. In addition, more than 21,000 publicly available genetic marker sequences were mapped to single-site positions. This new assembly represents a step forward to previously reported P. vannamei assemblies. It will be helpful as a reference genome for future studies on the evolutionary history of the species, the genetic architecture of physiological and sex-determination traits, and the analysis of the changes in genetic diversity and composition of cultivated stocks.

The anaplerotic pyruvate carboxylase from white shrimp Litopenaeus vannamei: Gene structure, molecular characterization, protein modelling and expression during hypoxia

Hypoxic zones are spreading worldwide in marine environments affecting many organisms. Shrimp and other marine crustaceans can withstand environmental hypoxia using several strategies, including the regulation of energy producing metabolic pathways. Pyruvate carboxylase (PC) catalyzes the first reaction of gluconeogenesis to produce oxaloacetate from pyruvate. In mammals, PC also participates in lipogenesis, insulin secretion and other processes, but this enzyme has been scarcely studied in marine invertebrates. In this work, we characterized the gene encoding PC in the white shrimp Litopenaeus vannamei, modelled the protein structure and evaluated its gene expression in hepatopancreas during hypoxia, as well as glucose and lactate concentrations. The PC gene codes for a mitochondrial protein and has 21 coding exons and 4 non-coding exons that generate three transcript variants with differences only in the 5'-UTR. Total PC expression is more abundant in hepatopancreas compared to gills or muscle, indicating tissue-specific expression. Under hypoxic conditions of 1.53 mg/L dissolved oxygen, PC expression is maintained in hepatopancreas, indicating its key role even in energy-limited conditions. Finally, both glucose and lactate concentrations were maintained under hypoxia for 24-48 h in hepatopancreas.

Characterization of Penaeus vannamei mitogenome focusing on genetic diversity

The diversity of the Penaeus vannamei mitochondrial genome has still been poorly characterized, there are no validated mitochondrial markers available for populational studies, and the heteroplasmy has not yet been investigated in this species. In this study, metagenomic reads extracted from the muscle of a single individual were used to assemble the mitochondrial genome (mtDNA). These data associated with mitochondrial genomes previously described allowed to evaluate the inter-individual variability and heteroplasmy. Comparison among 45 mtDNA control regions led to the detection of conserved and variable segments and the characterization of two hypervariable regions. The analysis of diversity revealed mostly low frequency polymorphisms, and heteroplasmy was found in practically all mitochondrial genes, with a high occurrence of indels. These results indicate that the design of mitochondrial markers for P. vannamei must be done with caution. The mapping of conserved and variable regions and the characterization of heteroplasmy presented here will contribute to increasing the efficiency of mitochondrial markers for population or individual studies.

Single nucleotide polymorphism associated with disease resistance in Penaeus vannamei

Despite the considerable number of genetic markers published for Penaeus vannamei, the classification of these markers and their standardization in specific databases is still insufficient. As a consequence, access to these markers is difficult, hampering their application in genetic association studies. In this study, all previously described single nucleotide polymorphisms (SNPs) related to resistance for P. vannamei were revised, and 512 SNPs were identified and classified in detail. We observed that most of the SNPs occurred in the proteins including Toll like receptors 1 and 3, hemocyanin large and small subunits, and anti-lipopolysaccharide factors 1 and 2, allowing to propose to use them as targets in association studies involving resistance in P. vannamei. Additionally, the potential effects of the most frequent non-synonymous coding SNPs in the secondary structure of the main target proteins were evaluated using an in silico approach. These data can serve as the starting point for the development of new genetic and computational tools as well as for the design of new association studies that involve resistance in P. vannamei.

Patterns of Geographical and Potential Adaptive Divergence in the Genome of the Common Carp (Cyprinus carpio)

The common carp, Cyprinus carpio, is a cyprinid fish species cultured in Europe and Asia. It accounts for >70% of freshwater aquaculture production worldwide. We conducted a population genomics analysis on C. carpio using high-throughput SNP genotyping of 2,198 individuals from 14 populations worldwide to determine the genetic architecture of common carp populations and the genetic bases for environmental adaptation. Structure analyses including phylogeny and principal component analysis were also conducted, showing distinct geographical patterns in European and Asian populations. The linkage disequilibrium block average lengths of the 14 populations ranged from 3.94 kb to 36.67 kb. Genes within selective sweep regions were identified by genome scanning among the different populations, including gdf6a, bmpr1b, and opsin5. Gene Ontology and KEGG enrichment analyses revealed potential trait-related loci and genes associated with body shape, scaling patterns, and skin color. This population genomics analysis may provide valuable clues for future genome-assisted breeding of C. carpio.