High-density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis

Genetic variation mining of the Chinese mitten crab (Eriocheir sinensis) based on transcriptome data from public databases

At present, public databases house an extensive repository of transcriptome data, with the volume continuing to grow at an accelerated pace. Utilizing these data effectively is a shared interest within the scientific community. In this study, we introduced a novel strategy that harnesses SNPs and InDels identified from transcriptome data, combined with sample metadata from databases, to effectively screen for molecular markers correlated with traits. We utilized 228 transcriptome datasets of Eriocheir sinensis from the NCBI database and employed the Genome Analysis Toolkit software to identify 96 388 SNPs and 20 645 InDels. Employing the genome-wide association study analysis, in conjunction with the gender information from databases, we identified 3456 sex-biased SNPs and 639 sex-biased InDels. The KOG and KEGG annotations of the sex-biased SNPs and InDels revealed that these genes were primarily involved in the metabolic processes of E. sinensis. Combined with SnpEff annotation and PCR experimental validation, a highly sex-biased SNP located in the Kelch domain containing 4 (Klhdc4) gene, CHR67-6415071, was found to alter the splicing sites of Klhdc4, generating two splice variants, Klhdc4_a and Klhdc4_b. Additionally, Klhdc4 exhibited robust expression across the ovaries, testes, and accessory glands. The sex-biased SNPs and InDels identified in this study are conducive to the development of unisexual cultivation methods for E. sinensis, and the alternative splicing event caused by the sex-biased SNP in Klhdc4 may serve as a potential mechanism for sex regulation in E. sinensis. The analysis strategy employed in this study represents a new direction for the rational exploitation and utilization of transcriptome data in public databases.

Multi-omics analysis of oxidative stress and apoptosis in hepatopancreas cells induced by Polyascus gregaria parasitizing the Eriocheir sinensis

Polyascus gregaria, a parasitic barnacle, poses a significant threat to Eriocheir sinensis farms by inhibiting crab growth. However, the molecular and pathological mechanisms behind P. gregaria infection in the hepatopancreas of E. sinensis remain unclear. In this study, we investigated the impact and underlying mechanisms of P. gregaria infection on E. sinensis through analyzing the infected hepatopancreatic tissues by tandem mass tag technology and RNA-Seq high-throughput sequencing. Among the identified 10,693 differentially expressed genes, 294 genes were significantly altered following P. gregaria infection, including 92 upregulated and 202 downregulated genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses further revealed the involvement of these genes in oxidative decomposition, lipid metabolism, inflammation, and hepatopancreas metabolism. Meanwhile, the identified 253 differentially expressed proteins, including 143 upregulated and 110 downregulated proteins, are mainly related to cellular and metabolic processes, catalytic activity, and cell components. The pathway analysis indicated their enrichment in glycolysis/gluconeogenesis, oxidative phosphorylation, endoplasmic reticulum protein processing, and actin cytoskeleton regulation. The involvement of these differentially expressed genes and proteins in the peroxisome proliferator-activated receptors pathway during host immune responses against P. gregaria infection has been highlighted. Furthermore, pathological examinations and biochemical indicators jointly demonstrated the hepatopancreatic damage and increased oxidative stress and apoptosis in the infected E. sinensis. Collectively, our study provides crucial insights into the mechanisms underlying the E. sinensis-P. gregaria interactions, and may contribute to the development of novel strategies for parasite control and reducing economic losses in aquaculture.

Evolution of sex determination in crustaceans

Sex determination (SD) involves mechanisms that determine whether an individual will develop into a male, female, or in rare cases, hermaphrodite. Crustaceans harbor extremely diverse SD systems, including hermaphroditism, environmental sex determination (ESD), genetic sex determination (GSD), and cytoplasmic sex determination (e.g., Wolbachia controlled SD systems). Such diversity lays the groundwork for researching the evolution of SD in crustaceans, i.e., transitions among different SD systems. However, most previous research has focused on understanding the mechanism of SD within a single lineage or species, overlooking the transition across different SD systems. To help bridge this gap, we summarize the understanding of SD in various clades of crustaceans, and discuss how different SD systems might evolve from one another. Furthermore, we review the genetic basis for transitions between different SD systems (i.e., Dmrt genes) and propose the microcrustacean Daphnia (clade Branchiopoda) as a model to study the transition from ESD to GSD.

Whole-genome resequencing reveals genetic differences and the genetic basis of parapodium number in Russian and Chinese Apostichopus japonicus

BACKGROUND

Apostichopus japonicus is an economically important species in the global aquaculture industry. Russian A. japonicus, mainly harvested in the Vladivostok region, exhibits significant phenotypic differentiation, including in many economically important traits, compared with Chinese A. japonicus owing to differences in their habitat. However, both the genetic basis for the phenotypic divergence and the population genetic structure of Russian and Chinese A. japonicus are unknown.

RESULT

In this study, 210 individuals from seven Russian and Chinese A. japonicus populations were sampled for whole-genome resequencing. The genetic structure analysis differentiated the Russian and Chinese A. japonicus into two groups. Population genetic analyses indicated that the Russian population showed a high degree of allelic linkage and had undergone stronger positive selection compared with the Chinese populations. Gene ontology terms enriched among candidate genes with group selection analysis were mainly involved in immunity, such as inflammatory response, antimicrobial peptides, humoral immunity, and apoptosis. Genome-wide association analysis yielded eight single-nucleotide polymorphism loci significantly associated with parapodium number, and these loci are located in regions with a high degree of genomic differentiation between the Chinese and Russia populations. These SNPs were associated with five genes. Gene expression validation revealed that three of these genes were significantly differentially expressed in individuals differing in parapodium number. AJAP08772 and AJAP08773 may directly affect parapodium production by promoting endothelial cell proliferation and metabolism, whereas AJAP07248 indirectly affects parapodium production by participating in immune responses.

CONCLUSIONS

This study, we performed population genetic structure and GWAS analysis on Chinese and Russian A. japonicus, and found three candidate genes related to the number of parapodium. The results provide an in-depth understanding of the differences in the genetic structure of A. japonicus populations in China and Russia, and provide important information for subsequent genetic analysis and breeding of this species.

The Sex-Specific Splicing of Doublesex in Brine Shrimp Artemia franciscana

The understanding of sex determination and differentiation in animals has recently made remarkable strides through the use of advanced research tools. At the gene level, the Mab-3-related transcription factor (Dmrt) gene family, which encodes for the typical DNA-binding doublesex/Mab-3 (DM) domain in their protein, is known for its contribution to sex determination and differentiation in insects. In this study, DNA-binding DM domain screening has identified eight transcripts from Artemia franciscana transcriptomic that encode proteins containing one conserved DNA-binding DM domain. The genome mapping confirmed that these eight transcripts are transcribed from six different loci on the A. franciscana genome assembly. One of those loci, the Af.dsx-4 locus, is closely related to Doublesex, a gene belonging to the Dmrt gene family. This locus could be transcribed into three alternative transcripts, namely Af.dsx⁴, Af.dsx^F and Af.dsx^M. While Af.dsx⁴ and Af.dsx^F could putatively be translated to form an identical Af.dsx^F protein of 186 aa long, Af.dsx^M translates for an Af.dsx^M protein of 289 aa long but shares a DNA-binding DM domain. Interestingly, Af.dsx^F and Af.dsx^M are confirmed as sex-specific transcripts, Af.dsx^F is only present in females, and Af.dsx^M is only present in male individuals. The results suggest that the sex-specific splicing mechanism of the doublesex described in insects is also present in A. franciscana. Af.dxs-4 locus can be used in further studies to clarify the sex determination pathways in A. fracnciscana.

GWAS Reveal Novel Sex-Related Markers and Candidate Genes in Sea Urchin Mesocentrotus nudus

Mesocentrotus nudus is an economically important mariculture species. Identification of sex-related markers and candidate genes has potential implications for sex-control breeding of Mesocentrotus nudus. The aim of this study was to identify the molecular markers and genes associated with the sex of M. nudus. Initial GWAS analysis based on 80 individuals genotyped by using GBS identified 22 sex-related SNPs located on 9 GBS tags. Further targeted sequencing in another population of 124 individuals confirmed that 7 SNPs located on 7 GBS tags displayed complete associations with sex, suggesting a ZW/ZZ sex-determination system in M. nudus. Besides, genome and transcriptome annotations presented that the heat shock protein 75 kDa, mitochondrial (trap1), and protein furry homolog-like (fryl) may be important candidate genes involved in sex determination and differentiation in M. nudus. However, further studies are needed to elucidate the functions of these genes. Overall, the current results not only provide molecular markers that may contribute to the sex-control breeding for M. nudus, but also provide new insights to understand the regulatory mechanism of sea urchin sex.

The Chinese mitten crab genome provides insights into adaptive plasticity and developmental regulation

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general.

QTL mapping and marker identification for sex determination in the ridgetail white prawn, Exopalaemon carinicauda

Sex determination is an important and intriguing research topic in the field of evolutionary and developmental biology. Quantitative trait locus (QTL) mapping for sex is helpful in clarifying the sex determination system of species. In this study, a second high-resolution genetic linkage map was constructed for the ridgetail white prawn, Exopalaemon carinicauda, which included 9280 markers, covering 99.98% of the complete genome. Based on the linkage map, a highly significant sex-related QTL was first mapped to a single linkage group (LG3, LOD > 55.6). Fifty-two markers in the QTL region were significantly associated with sex (p ≤ 10^-40), of which heterogametic genotypes in females supported the ZW sex determination mechanism. Six markers were verified to be significantly associated with sex in the wild population. Some sex-related genes were identified, including phospholipase D, protein kinase shaggy, and longitudinals lacking protein. These results inform our understanding of the mechanisms of sex determination in E. carinicauda.

High accuracy of pooled DNA genotyping by 2b-RAD sequencing in the Pacific white shrimp, Litopenaeus vannamei

Using pooled DNA genotyping to estimate the proportional contributions from multiple families in a pooled sample is of particular interest for selective breeding in aquaculture. We compared different pooled libraries with separate 2b-RAD sequencing of Litopenaeus vannamei individuals to assess the effect of different population structures (different numbers of individuals and families) on pooled DNA sequencing, the accuracy of parent sequencing of the DNA pools and the effect of SNP numbers on pooled DNA sequencing. We demonstrated that small pooled DNA genotyping of up to 53 individuals by 2b-RAD sequencing could provide a highly accurate assessment of population allele frequencies. The accuracy increased as the number of individuals and families increased. The allele frequencies of the parents from each pool were highly correlated with those of the pools or the corresponding individuals in the pool. We chose 500-28,000 SNPs to test the effect of SNP number on the accuracy of pooled sequencing, and no linear relationship was found between them. When the SNP number was fixed, increasing the number of individuals in the mixed pool resulted in higher accuracy of each pooled genotyping. Our data confirmed that pooled DNA genotyping by 2b-RAD sequencing could achieve higher accuracy than that of individual-based genotyping. The results will provide important information for shrimp breeding programs.

Sex-Specific Transcriptome Sequencing of Zoea I Larvae and Identification of Sex-Linked Genes Using Bulked Segregant Analysis in Pacific White Shrimp Litopenaeus vannamei

In order to learn the sex determination and differentiation process in shrimp, we performed sex-specific transcriptome sequencing of the Zoea I larvae in Litopenaeus vannamei (L. vannamei) using a DNA/RNA co-extraction method. The sex-specific genes and the sex-linked SNPs were identified. De novo assembly of all the clean reads generated 41,270 unigenes with an average length of 1026 bp and an N50 of 1922 bp. A total of 60 differentially expressed genes (DEGs) between females and males at Zoea I stage were obtained, in which 41 unigenes had annotations. DEGs that might be related to sex development in L. vannamei were analyzed. Besides, paralogous genes of some known sex-related genes in model organisms such as sex-lethal (Sxl) showed no sex-biased difference, which indicated they may not play roles in sex development at Zoea I stage. Candidate sex-linked SNPs were validated in six populations of L. vannamei, and SNPs in Unigene0020898 and Unigene0020336 were fully linked to genders, suggesting the two genes located in the sex-determining region of L. vannamei. Our report provided the sex-specific gene expression in early developmental stage of L. vannamei, and found some sex-biased genes that may participate in the sex development process. We also detected two sex-linked genes which may be located in the sex determination region of L. vannamei through bulked segregant analysis, and the SNPs on these genes also provided a method for genetic sex identification at transcriptional level. Our findings will facilitate further researches on molecular mechanism of sex determination and differentiation in shrimps.

Accurate Sex Identification of Ancient Elephant and Other Animal Remains Using Low-Coverage DNA Shotgun Sequencing Data

Sex identification of ancient animal biological remains can benefit our understanding of historical population structure, demography and social behavior. Traditional methods for sex identification (e.g., osteological and morphometric comparisons) may be ineffective when animal remains are not well preserved, when sex distinguishing characteristics have not yet developed, or where organisms do not exhibit sex-associated phenotypic dimorphisms. Here we adapt a method developed for human sex determination so that it can be used to identify the sex of ancient and modern animal taxa. The method identifies sex by calculating the ratio of DNA reads aligning to the X chromosome to DNA reads aligning to autosomes (termed the Rx ratio). We tested the accuracy of this method using low coverage genomes from 15 modern elephants (Loxodonta africana) for which sex was known. We then applied this method to ancient elephant ivory samples for which sex was unknown, and describe how this method can be further adapted to the genomes of other taxa. This method may be especially useful when only low-coverage genomic data are obtainable. Furthermore, because this method relies on only the X and not the Y chromosome, it can be used to determine the sex of organisms for which a reference genome was obtained from a female or for which only the X chromosome is reported. Such taxa include the domestic cat, sheep, goat, and horse; and non-domesticated animals such as the Sumatran orangutan, western lowland gorilla and meerkat.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Validation of growth-related quantitative trait loci markers in different Exopalaemon carinicauda families for marker-assisted selection

We detected growth-related QTL and associated markers from the backcross population of Exopalaemon carinicauda in the previous study. Based on our previous study, the 47 SNP markers associated with candidate growth trait QTL were selected to analyze the association between these markers and body weight (BW), body length and abdominal segment length traits in four different populations including wild population, a full-sib family, a half-sib family and a backcross population for evaluating their potential application of marker-assisted selection in E. carinicauda. The general linear model (GLM) and mixed linear model were applied and the associations between SNP loci and three growth-related traits verified. The results showed that the Marker79268 and Marker100644 were significantly associated with the BW trait in more than three populations by the GLM method. The Marker100644 was significantly associated with BW in the full-sib family, half-sib family and backcross populations by the GLM and mixed linear model methods. Our findings will provide useful SNP markers to go forward to improve growth performance through more refined marker-assisted selection in E. carinicauda.

Quantitative Trait Loci Mapping and Marker Identification for Low Salinity Tolerance Trait in the Swimming Crab (Portunus trituberculatus)

Low salinity is one of the most important abiotic factors that directly affect the abundance of the swimming crab, Portunus trituberculatus. Quantitative trait loci (QTL) mapping could be helpful in identifying the markers and genes involved in low salinity tolerance. In this study, two QTLs of low salt tolerance were mapped on linkage group 17 (LG17, 2.6-5.2 cM) based on a high-density linkage map. Ninety-five markers related to low salinity tolerance were identified via association analysis, and seventy-nine low salt-related candidate genes (including ammonium transport, aldehyde dehydrogenase, and glucosyltransferase) were screened from draft genome of the species via these markers. This represents the first report of QTL mapping for low salinity tolerance in the swimming crab, which may be useful to elucidate salinity adaptation mechanisms.

Sex identification from distinctive gene expression patterns in Antarctic krill (Euphausia superba)

Antarctic krill (Euphausia superba) is a highly abundant keystone species of the Southern Ocean ecosystem, directly connecting primary producers to high-trophic level predators. Sex ratios of krill vary remarkably between swarms and this phenomenon is poorly understood, as identification of krill sex relies on external morphological differences that appear late during development. Sex determination mechanisms in krill are unknown, but could include genetic, environmental or parasitic mechanisms. Similarly, virtually nothing is known about molecular sex differentiation. The krill genome has to date not been sequenced, and due to its enormous size and large amount of repetitive elements, it is currently not feasible to develop sex-specific DNA markers. To produce a reliable molecular marker for sex in krill and to investigate molecular sex differentiation we therefore focused on identifying sex-specific transcriptomic differences. Through transcriptomic analysis, we found large gene expression differences between testes and ovaries and identified three genes exclusively expressed in female whole krill from early juvenile stages onwards. The sex-specific expression of these three genes persisted through sexual regression, although our regressed samples originated from a krill aquarium and may differ from wild-regressed krill. Two slightly male-biased genes did not display sufficient expression differences to clearly differentiate sexes. Based on the expression of the three female-specific genes we developed a molecular test that for the first time allows the unambiguous sex determination of krill samples lacking external sex-specific features from juvenile stages onwards, including the sexually regressed krill we examined.

Molecular Characterization and Functional Study of Insulin-Like Androgenic Gland Hormone Gene in the Red Swamp Crayfish, Procambarus clarkii

The androgenic gland (AG) is a male-specific endocrine organ that controls the primary and secondary sexual characteristics in male crustaceans. More evidence indicates that the insulin-like androgenic gland hormone gene (IAG) is the key male sexual differentiation factor, particularly the application of RNA interference (RNAi) technology on IAG. In this study, the full-length cDNA of IAG (termed PcIAG) was isolated from the red swamp crayfish, Procambarusclarkii. Tissue distribution analysis showed that in addition to its expression in the AG of male P. clarkii, PcIAG was widely expressed in female tissues and other male tissues. The PcIAG protein was detected in the reproductive and nervous systems of adult male P. clarkii. Additionally, RNAi results showed that the PcIAG expression could be silenced efficiently, and the male sperm maturation and release possibly present a transient adverse interference at lower doses (0.1 μg/g and 1 μg/g) of PcIAG–dsRNA (PcIAG double-stranded RNA). Dramatically, the expression level of PcIAG increased sharply shortly after the injection of higher doses (5 μg/g and 10 μg/g) of PcIAG–dsRNA, which might accelerate the maturation and release of sperm. Moreover, the expression of PcSxl (P. clarkii Sex-lethal) was detected by Quantitative Real-Time PCR (qPCR) after the injection of PcIAG–dsRNA to explore whether the PcIAG gene regulates the PcSxl gene, and we found that the PcIAG did not directly regulate PcSxl in P. clarkii. The study could help accelerate the progress of PcIAG functional research and provide a useful reference for the single-sex selective breeding of P. clarkii.

Genome survey and high-resolution backcross genetic linkage map construction of the ridgetail white prawn Exopalaemon carinicauda applications to QTL mapping of growth traits

BACKGROUND

High-resolution genetic linkage map is critical for QTL mapping, genome sequence assembly and marker-assisted selection in aquaculture species. The ridgetail white prawn Exopalaemon carinicauda is one of the most economic shrimp species naturally distributed in the coasts of eastern China and western Korea. However, quite limited genomics and genetics information have been exploited for genetic improvement of economic traits in this species.

RESULTS

In the present study, we conducted genome survey and constructed high-resolution genetic linkage maps of the ridgetail white prawn with reciprocal-cross mapping family genotyped using next-generation sequencing approaches. The estimated genome size was 9.33 Gb with a heterozygosity of 0.26% and a repeat sequence ratio of 76.62%. 65,772 protein-coding genes were identified by genome annotation. A total of 10,384 SNPs were used to high-throughput genotyping and assigned to 45 linkage groups (LGs) from reciprocal backcross families of E. carinicauda, and the average marker distances were 0.73 cM and 0.55 cM, respectively. Based on the high-resolution linkage map, twenty-three QTLs related to five growth traits were detected. All QTLs could explain 8.8-15.7% of the total growth-traits variation.

CONCLUSIONS

The genome size of E. carinicauda was estimated more accurately by genome survey analysis, which revealed basic genomic architecture. The first high-resolution backcross genetic linkage map and QTLs related to growth traits will provide important information for QTL fine mapping, genome assembly and genetic improvement of E. carinicauda and other palaemon shrimps.

Sex-Biased Gene Expression and Dosage Compensation on the Artemia franciscana Z-Chromosome

Males and females of Artemia franciscana, a crustacean commonly used in the aquarium trade, are highly dimorphic. Sex is determined by a pair of ZW chromosomes, but the nature and extent of differentiation of these chromosomes is unknown. Here, we characterize the Z chromosome by detecting genomic regions that show lower genomic coverage in female than in male samples, and regions that harbor an excess of female-specific SNPs. We detect many Z-specific genes, which no longer have homologs on the W, but also Z-linked genes that appear to have diverged very recently from their existing W-linked homolog. We assess patterns of male and female expression in two tissues with extensive morphological dimorphism, gonads, and heads. In agreement with their morphology, sex-biased expression is common in both tissues. Interestingly, the Z chromosome is not enriched for sex-biased genes, and seems to in fact have a mechanism of dosage compensation that leads to equal expression in males and in females. Both of these patterns are contrary to most ZW systems studied so far, making A. franciscana an excellent model for investigating the interplay between the evolution of sexual dimorphism and dosage compensation, as well as Z chromosome evolution in general.

High-Density Genetic Linkage Maps Provide Novel Insights Into ZW/ZZ Sex Determination System and Growth Performance in Mud Crab (Scylla paramamosain)

Mud crab, Scylla paramamosain is one of the most important crustacean species in global aquaculture. To determine the genetic basis of sex and growth-related traits in S. paramamosain, a high-density genetic linkage map with 16,701 single nucleotide polymorphisms (SNPs) was constructed using SLAF-seq and a full-sib family. The consensus map has 49 linkage groups, spanning 5,996.66 cM with an average marker-interval of 0.81 cM. A total of 516 SNP markers, including 8 female-specific SNPs segregated in two quantitative trait loci (QTLs) for phenotypic sex were located on LG32. The presence of female-specific SNP markers only on female linkage map, their segregation patterns and lower female: male recombination rate strongly suggest the conformation of a ZW/ZZ sex determination system in S. paramamosain. The QTLs of most (90%) growth-related traits were found within a small interval (25.18–33.74 cM) on LG46, highlighting the potential involvement of LG46 in growth. Four markers on LG46 were significantly associated with 10–16 growth-related traits. BW was only associated with marker 3846. Based on the annotation of transcriptome data, 11 and 2 candidate genes were identified within the QTL regions of sex and growth-related traits, respectively. The newly constructed high-density genetic linkage map with sex-specific SNPs, and the identified QTLs of sex- and growth-related traits serve as a valuable genetic resource and solid foundation for marker-assisted selection and genetic improvement of crustaceans.

Construction of a Genomic Bacterial Artificial Chromosome (BAC) Library for the Prawn Macrobrachium rosenbergii and Initial Analysis of ZW Chromosome-Derived BAC Inserts

Knowledge on sex determination has proven valuable for commercial production of the prawn Macrobrachium rosenbergii due to sex dimorphism of the male and female individuals. Previous studies indicated that prawn sex is determined by a ZW-ZZ chromosomal system, but no genomic information is available for the sex chromosome. Herein, we constructed a genomic bacterial artificial chromosome (BAC) library and identified the ZW-derived BAC clones for initial analysis of the sex chromosomal DNA sequence. The arrayed BAC library contains 200,448 clones with average insert size of 115.4 kb, corresponding to ∼ 4× coverage of the estimated 5.38 Gb genome. Based on a short female-specific marker, a Z- and a W-fragment were retrieved with the genomic walking method. Screening the BAC library using a ZW-specific marker as probe resulted in 12 positive clones. From these, a Z-derived (P331M17) and a W-derived (P122G2) BAC clones were randomly selected and sequenced by PacBio method. We report the construction of a large insert, deep-coverage, and high-quality BAC library for M. rosenbergii that provides a useful resource for positional cloning of target genes, genomic organization, and comparative genomics analysis. Our study not only confirmed the ZW/ZZ system but also discovered sex-linked genes on ZW chromosomes for the first time, contributing to a comprehensive understanding of the genomic structure of sex chromosomes in M. rosenbergii.

Construction of a High-Density Linkage Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Channel Catfish (Ictalurus punctatus)

A high-density genetic linkage map is of particular importance in the fine mapping for important economic traits and whole genome assembly in aquaculture species. The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world. Outside of the United States, China has become the major producer and consumer of channel catfish after experiencing rapid development in the past three decades. In this study, based on restriction site associated DNA sequencing (RAD-seq), a high-density genetic linkage map of channel catfish was constructed by using single nucleotide polymorphisms (SNPs) in a F₁ family composed of 156 offspring and their two parental individuals. A total of 4,768 SNPs were assigned to 29 linkage groups (LGs), and the length of the linkage map reached 2,480.25 centiMorgans (cM) with an average distance of 0.55 cM between loci. Based on this genetic linkage map, 223 genomic scaffolds were anchored to the 29 LGs of channel catfish, and a total length of 704.66 Mb was assembled. Quantitative trait locus (QTL) mapping and genome-wide association analysis identified 10 QTLs of sex-related and six QTLs of growth-related traits at LG17 and LG28, respectively. Candidate genes associated with sex dimorphism, including spata2, spata5, sf3, zbtb38, and fox, were identified within QTL intervals on the LG17. A sex-linked marker with simple sequence repeats (SSR) in zbtb38 gene of the LG17 was validated for practical verification of sex in the channel catfish. Thus, the LG17 was considered as a sex-related LG. Potential growth-related genes were also identified, including important regulators such as megf9, npffr1, and gas1. In a word, we constructed the high-density genetic linkage map and developed the sex-linked marker in channel catfish, which are important genetic resources for future marker-assisted selection (MAS) of this economically important teleost.

Identification of male-specific SNP markers and development of PCR-based genetic sex identification technique in crucifix crab (Charybdis feriatus) with implication of an XX/XY sex determination system

In this study, we first identified male-specific SNP markers using restriction site-associated DNA sequencing, and further developed a PCR-based sex identification technique for Charybdis feriatus. A total of 296.96 million clean reads were obtained, with 114.95 and 182.01 million from females and males. After assembly and alignment, 10 SNP markers were identified being heterozygous in males but homozygous in females. Five markers were further confirmed to be male-specific in a large number of individuals. Moreover, two male-specific sense primers and a common antisense primer were designed, using which, a PCR-based genetic sex identification method was successfully developed and used to identify the sex of 103 individuals, with a result of 49 females and 54 males. The presence of male-specific SNP markers suggests an XX/XY sex determination system for C. feriatus. These findings should be helpful for better understanding sex determination mechanism, and drafting artificial breeding program in crustaceans.

Genome survey, high-resolution genetic linkage map construction, growth-related quantitative trait locus (QTL) identification and gene location in Scylla paramamosain

Scylla paramamosain is one of the most economically important crabs in China. In this study, the first genome survey sequencing of this crab was performed, and the results revealed that the estimated genome size was 1.21 Gb with high heterozygosity (1.3%). Then, RAD technology was used to construct a high-resolution linkage map for this species. A total of 24,444 single nucleotide polymorphism (SNP) makers were grouped into 47 linkage groups. The total length of the linkage groups was 3087.53 cM with a markers interval of 0.92 cM. With the aid of transcriptome and genome scaffold data, 4,271 markers were linked to genes, including several important growth-related genes such as transforming growth factor-beta regulator I, immune related-gene C-type lectin and ecdysone pathway gene broad-complex-like protein. Further, 442 markers, representing 279 QTLs, associated with 24 traits were identified, and of these markers, 78 were linked to genes. Some interesting genes, such as dedicator of cytokinesis protein 3, tenascin-X and DNA helicase MCM8, were believed to have important relationship with specific traits and merit further exploration. The results of this study will accelerate the genetic improvement and genome sequencing analysis of the mud crab.

Female-specific SNP markers provide insights into a WZ/ZZ sex determination system for mud crabs Scylla paramamosain, S. tranquebarica and S. serrata with a rapid method for genetic sex identification

BACKGROUND

Mud crabs, Scylla spp., are commercially important large-size marine crustaceans in the Indo-West Pacific region. As females have the higher growth rate and economic value, the production of all female stocks is extremely essential in aquaculture. However, the sex determination mechanism is still unclear. Development of sex-specific genetic markers based on next-generation sequencing proved to be an effective tool for discovering sex determination system in various animals.

RESULTS

Restriction-site associated DNA sequencing (RAD-seq) was employed to isolate sex-specific SNP markers for S. paramamosain. A total of 335.6 million raw reads were obtained from 20 individuals, of which 204.7 million were from 10 females and 130.9 million from 10 males. After sequence assembly and female-male comparison, 20 SNP markers were identified to be sex-specific. Furthermore, ten SNPs in a short sequence (285 bp) were confirmed heterozygous in females and homozygous in males in a large population by PCR amplification and sequencing. Subsequently, a female-specific primer was successfully designed according to the female-specific nucleotide which could amplify an expected band from females but not from males. Thus, a rapid and effective method for molecular sexing in S. paramamosain was developed, meanwhile, this method could successfully identify the sex of S. tranquebarica and S. serrata. Finally, nine and four female-specific SNP markers were detected in S. tranquebarica and S. serrata, respectively.

CONCLUSIONS

Sex-specific SNP markers were firstly identified in crab species and showed female heterogamety and male homogamety, which provided strong genetic evidence for a WZ/ZZ sex determination system in mud crabs S. paramamosain, S. tranquebarica and S. serrata. These findings will lay a solid foundation for the study of sex determination mechanism, sex chromosome evolution, and the development of mono-sex population in crustaceans.

High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing

BACKGROUND

Flavonoids are important components that confer upon tea plants a unique flavour and health functions. However, the traditional breeding method for selecting a cultivar with a high or unique flavonoid content is time consuming and labour intensive. High-density genetic map construction associated with quantitative trait locus (QTL) mapping provides an effective way to facilitate trait improvement in plant breeding. In this study, an F1 population (LJ43×BHZ) was genotyped using 2b-restriction site-associated DNA (2b-RAD) sequencing to obtain massive single nucleotide polymorphism (SNP) markers to construct a high-density genetic map for a tea plant. Furthermore, QTLs related to flavonoids were identified using our new genetic map.

RESULTS

A total of 13,446 polymorphic SNP markers were developed using 2b-RAD sequencing, and 4,463 of these markers were available for constructing the genetic linkage map. A 1,678.52-cM high-density map at an average interval of 0.40 cM with 4,217 markers, including 427 frameset simple sequence repeats (SSRs) and 3,800 novel SNPs, mapped into 15 linkage groups was successfully constructed. After QTL analysis, a total of 27 QTLs related to flavonoids or caffeine content (CAF) were mapped to 8 different linkage groups, LG01, LG03, LG06, LG08, LG10, LG11, LG12, and LG13, with an LOD from 3.14 to 39.54, constituting 7.5% to 42.8% of the phenotypic variation.

CONCLUSIONS

To our knowledge, the highest density genetic map ever reported was constructed since the largest mapping population of tea plants was adopted in present study. Moreover, novel QTLs related to flavonoids and CAF were identified based on the new high-density genetic map. In addition, two markers were located in candidate genes that may be involved in flavonoid metabolism. The present study provides valuable information for gene discovery, marker-assisted selection breeding and map-based cloning for functional genes that are related to flavonoid content in tea plants.

A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

Divergent Evolutionary Trajectories of Two Young, Homomorphic, and Closely Related Sex Chromosome Systems

There exists extraordinary variation among species in the degree and nature of sex chromosome divergence. However, much of our knowledge about sex chromosomes is based on comparisons between deeply diverged species with different ancestral sex chromosomes, making it difficult to establish how fast and why sex chromosomes acquire variable levels of divergence. To address this problem, we studied sex chromosome evolution in two species of African clawed frog (Xenopus), both of whom acquired novel systems for sex determination from a recent common ancestor, and both of whom have female (ZW/ZZ) heterogamy. Derived sex chromosomes of one species, X. laevis, have a small region of suppressed recombination that surrounds the sex determining locus, and have remained this way for millions of years. In the other species, X. borealis, a younger sex chromosome system exists on a different pair of chromosomes, but the region of suppressed recombination surrounding an unidentified sex determining gene is vast, spanning almost half of the sex chromosomes. Differences between these sex chromosome systems are also apparent in the extent of nucleotide divergence between the sex chromosomes carried by females. Our analyses also indicate that in autosomes of both of these species, recombination during oogenesis occurs more frequently and in different genomic locations than during spermatogenesis. These results demonstrate that new sex chromosomes can assume radically different evolutionary trajectories, with far-reaching genomic consequences. They also suggest that in some instances the origin of new triggers for sex determination may be coupled with rapid evolution sex chromosomes, including recombination suppression of large genomic regions.

QTL Mapping and Marker Identification for Sex-Determining: Indicating XY Sex Determination System in the Swimming Crab (Portunus trituberculatus)

Sex determination is an important area of research, which has always had an intriguing aspect in evolutionary and developmental biology. Quantitative trait locus (QTL) mapping for sex will be helpful in clarifying the sex determination system. In this study, the sex QTL mapping of the swimming crab (Portunus trituberculatus) was performed based on a high-density linkage map, and a highly significant QTL specifically mapped on a single linkage group (LG) was firstly identified (LG24, LOD > 14). Twenty markers in the QTL region showed significant associations with sex by association analysis, of which heterogametic genotypes in males supported the XY sex determination mechanism. Two sex-specific markers at the family level were identified via segregation distortion analysis, which were known to be the most closely linked to the sex of P. trituberculatus. Based on sex marker sequences (Marker3840, Marker20320, and Marker10494), three potential sex-related genes were identified, and the quantitative real-time PCR results suggested that these genes were important in spermatogenesis or sex characteristics in males. Our results will contribute to the fine-mapping of sex-determining genes and clarify the sex determination mechanism of P. trituberculatus.

A high-resolution genetic linkage map and QTL fine mapping for growth-related traits and sex in the Yangtze River common carp (Cyprinus carpio haematopterus)

Background

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection for economic traits in aquaculture species. The Yangtze River common carp (Cyprinus carpio haematopterus) is one of the most important aquacultured strains in China. However, quite limited genetics and genomics resources have been developed for genetic improvement of economic traits in such strain.

Results

A high-resolution genetic linkage map was constructed by using 7820 2b-RAD (2b-restriction site-associated DNA) and 295 microsatellite markers in a F2 family of the Yangtze River common carp (C. c. haematopterus). The length of the map was 4586.56 cM with an average marker interval of 0.57 cM. Comparative genome mapping revealed that a high proportion (70%) of markers with disagreed chromosome location was observed between C. c. haematopterus and another common carp strain (subspecies) C. c. carpio. A clear 2:1 relationship was observed between C. c. haematopterus linkage groups (LGs) and zebrafish (Danio rerio) chromosomes. Based on the genetic map, 21 QTLs for growth-related traits were detected on 12 LGs, and contributed values of phenotypic variance explained (PVE) ranging from 16.3 to 38.6%, with LOD scores ranging from 4.02 to 11.13. A genome-wide significant QTL (LOD = 10.83) and three chromosome-wide significant QTLs (mean LOD = 4.84) for sex were mapped on LG50 and LG24, respectively. A 1.4 cM confidence interval of QTL for all growth-related traits showed conserved synteny with a 2.06 M segment on chromosome 14 of D. rerio. Five potential candidate genes were identified by blast search in this genomic region, including a well-studied multi-functional growth related gene, Apelin.

Conclusions

We mapped a set of suggestive and significant QTLs for growth-related traits and sex based on a high-density genetic linkage map using SNP and microsatellite markers for Yangtze River common carp. Several candidate growth genes were also identified from the QTL regions by comparative mapping. This genetic map would provide a basis for genome assembly and comparative genomics studies, and those QTL-derived candidate genes and genetic markers are useful genomic resources for marker-assisted selection (MAS) of growth-related traits in the Yangtze River common carp.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4613-1) contains supplementary material, which is available to authorized users.

Comparative cytogenetic analysis of marine Palaemon species reveals a X1X1X2X2/X1X2Y sex chromosome system in Palaemon elegans

Background

The maintenance of species and the promotion of speciation are closely related to chromosomal rearrangements throughout evolution. Decapoda represents the most species-rich order among crustaceans and, despite its ecological and economic importance, little is known about decapod karyology. We aim at cytogenetically characterizing two sympatric prawn species.

Results

Analysis of mitotic metaphases and meiotic diakinesis of the common prawn Palaemon serratus and the rockpool prawn P. elegans, revealed considerable differences between their karyotypes including chromosome numbers and sex determination systems. The cytogenetic data for P. serratus showed a diploid number of 56 and the putative absence of heteromorphic sex chromosomes. However, the diploid chromosome number in P. elegans was 90 for females and 89 for males. The karyotype of the females consisted of the three largest acrocentric pairs and 42 submetacentric and metacentric pairs, while the karyotype of the males comprised a clearly identifiable large metacentric chromosome and two acrocentric pairs as well as the smaller 42 pairs. These results highlight the presence of the X₁X₁X₂X₂/X₁X₂Y multiple sex chromosome system in P. elegans, which constitute the only sexual system for Decapoda reported cytogenetically using modern techniques. The origin of this sex chromosome system is discussed. We hypothesize that the chromosome evolution within the genus could involve several fusion events giving rise to a reduction on the chromosome number in P. serratus. In both species, the major ribosomal genes were located in two chromosome pairs and hybridization signals of the telomeric sequences (TTAGGG)_n were visualized at the telomeres of all chromosomes. C-banding revealed that, when present, constitutive heterochromatin had a predominantly telomeric distribution and no centromeric constitutive heterochromatin was observed.

Conclusions

Although more comparative cytogenetic analyses are needed to clarify our hypotheses, the findings of this work indicate that the prawns of the genus Palaemon represent a promising model among Decapoda representatives to investigate the karyotype evolution and the patterns of sex chromosome differentiation.

Linkage mapping aided by de novo genome and transcriptome assembly in Portunus trituberculatus: applications in growth-related QTL and gene identification

A high-resolution genetic linkage map is an essential tool for decoding genetics and genomics in non-model organisms. In this study, a linkage map was constructed for the swimming crab (Portunus trituberculatus) with 10,963 markers; as far as we know, this number of markers has never been achieved in any other crustacean. The linkage map covered 98.85% of the whole genome with a mean marker interval of 0.51 cM. The de novo assembly based on genome and transcriptome sequencing data enabled 2,378 explicit annotated markers to be anchored to the map. Quantitative trait locus (QTL) mapping revealed 10 growth-related QTLs with a phenotypic variance explained (PVE) range of 12.0-35.9. Eight genes identified from the growth-related QTL regions, in particular, RE1-silencing transcription factor and RNA-directed DNA polymerase genes with nonsynonymous substitutions, were considered important growth-related candidate genes. We have demonstrated that linkage mapping aided by de novo assembly of genome and transcriptome sequencing could serve as an important platform for QTL mapping and the identification of trait-related genes.

A High-Density Genetic Linkage Map and QTL Fine Mapping for Body Weight in Crucian Carp (Carassius auratus) Using 2b-RAD Sequencing

A high-resolution genetic linkage map is essential for a wide range of genetics and genomics studies such as comparative genomics analysis and QTL fine mapping. Crucian carp (Carassius auratus) is widely distributed in Eurasia, and is an important aquaculture fish worldwide. In this study, a high-density genetic linkage map was constructed for crucian carp using 2b-RAD technology. The consensus map contains 8487 SNP markers, assigning to 50 linkage groups (LGs) and spanning 3762.88 cM, with an average marker interval of 0.44 cM and genome coverage of 98.8%. The female map had 4410 SNPs, and spanned 3500.42 cM (0.79 cM/marker), while the male map had 4625 SNPs and spanned 3346.33 cM (0.72 cM/marker). The average recombination ratio of female to male was 2.13:1, and significant male-biased recombination suppressions were observed in LG47 and LG49. Comparative genomics analysis revealed a clear 2:1 syntenic relationship between crucian carp LGs and chromosomes of zebrafish and grass carp, and a 1:1 correspondence, but extensive chromosomal rearrangement, between crucian carp and common carp, providing evidence that crucian carp has experienced a fourth round of whole genome duplication (4R-WGD). Eight chromosome-wide QTL for body weight at 2 months after hatch were detected on five LGs, explaining 10.1-13.2% of the phenotypic variations. Potential candidate growth-related genes, such as an EGF-like domain and TGF-β, were identified within the QTL intervals. This high-density genetic map and QTL analysis supplies a basis for genome evolutionary studies in cyprinid fishes, genome assembly, and QTL fine mapping for complex traits in crucian carp.

Two spliced isoforms of the sex-determination gene fruitless in the Chinese mitten crab Eriocheir sinensis

The fruitless (fru) gene plays an important role in the sex determination pathway and courtship behavior of Drosophila melanogaster. In the present study, two fru isoforms (Esfru1 and Esfru2) were identified from the Chinese mitten crab Eriocheir sinensis. Sequence analysis showed that Esfru1 and Esfru2 were encoded by the same genomic locus and generated by alternative splicing of pre-mRNA. Esfru1 had all introns completely spliced out, while Esfru2 had a longer exon1 with an additional 78bp sequence. They both contained a conserved BTB domain which was also found in D. melanogaster Fru isoforms. Analysis on temporal expression profiles of Esfru1-2 (a common region of Esfru1 and Esfru2) and Esfru2 (a specific region of Esfru2) showed that they expressed in similar patterns during embryonic development but in different patterns during larval development. The expression of Esfru1-2 decreased gradually from zoea III stage, however, Esfru2 increased from zoea IV stage and reached the peak at megalopa stage. Expression distribution in tissues and in situ hybridization analysis revealed that they showed sexually dimorphic expression in gonads, hepatopancreas and brains. Esfru1-2 showed significantly higher expression in female gonads and hepatopancreas than in males, but highly expressed in male brains than that in females. Interestingly, Esfru2 was displayed in a male-specific manner in each tissue, especially in gonads and hepatopancreas. Our results indicate that Esfru1 might be involved in both sexual brain neuronal structure development and related to female-specific character development. Esfru2 may participate in male-specific character development. This is the first report that characterizes two spliced variants of fru in crustaceans and provides basic information for further functional studies of the crab sex-determination mechanism.

A second generation SNP and SSR integrated linkage map and QTL mapping for the Chinese mitten crab Eriocheir sinensis

The Chinese mitten crab Eriocheir sinensis is the most economically important cultivated crab species in China, and its genome has a high number of chromosomes (2n = 146). To obtain sufficient markers for construction of a dense genetic map for this species, we employed the recently developed specific-locus amplified fragment sequencing (SLAF-seq) method for large-scale SNPs screening and genotyping in a F1 full-sib family of 149 individuals. SLAF-seq generated 127,677 polymorphic SNP markers, of which 20,803 valid markers were assigned into five segregation types and were used together with previous SSR markers for linkage map construction. The final integrated genetic map included 17,680 SNP and 629 SSR markers on the 73 linkage groups (LG), and spanned 14,894.9 cM with an average marker interval of 0.81 cM. QTL mapping localized three significant growth-related QTL to a 1.2 cM region in LG53 as well as 146 sex-linked markers in LG48. Genome-wide QTL-association analysis further identified four growth-related QTL genes named LNX2, PAK2, FMRFamide and octopamine receptors. These genes are involved in a variety of different signaling pathways including cell proliferation and growth. The map and SNP markers described here will be a valuable resource for the E. sinensis genome project and selective breeding programs.

Serial sequencing of isolength RAD tags for cost-efficient genome-wide profiling of genetic and epigenetic variations

Isolength restriction site-associated DNA (isoRAD) sequencing is a very simple but powerful approach that was originally developed for genome-wide genotyping at minimal labor and cost, and it has recently extended its applicability to allow quantification of DNA methylation levels. The isoRAD method is distinct from other genotyping-by-sequencing (GBS) methods because of its use of special restriction enzymes to produce isolength tags (32-36 bp), and sequencing of these uniform tags can bring many benefits. However, the relatively short tags produced by the original protocol are mostly suited to single-end (SE) sequencing (36-50 bp), and therefore they cannot efficiently match the gradually increased sequencing capacity of next-generation sequencing (NGS) platforms. To address this issue, we describe an advanced protocol that allows the preparation of five concatenated isoRAD tags for Illumina paired-end (PE) sequencing (100-150 bp). The configuration of the five concatenated tags is highly flexible, and can be defined by users to work with a desired combination of samples and/or restriction enzymes to suit specific research purposes. In comparison with the original protocol, the advanced protocol has an additional digestion and ligation step, and library preparation can be completed in ∼8 h.

MethylRAD: a simple and scalable method for genome-wide DNA methylation profiling using methylation-dependent restriction enzymes

Characterization of dynamic DNA methylomes in diverse phylogenetic groups has attracted growing interest for a better understanding of the evolution of DNA methylation as well as its function and biological significance in eukaryotes. Sequencing-based methods are promising in fulfilling this task. However, none of the currently available methods offers the 'perfect solution', and they have limitations that prevent their application in the less studied phylogenetic groups. The recently discovered Mrr-like enzymes are appealing for new method development, owing to their ability to collect 32-bp methylated DNA fragments from the whole genome for high-throughput sequencing. Here, we have developed a simple and scalable DNA methylation profiling method (called MethylRAD) using Mrr-like enzymes. MethylRAD allows for de novo (reference-free) methylation analysis, extremely low DNA input (e.g. 1 ng) and adjustment of tag density, all of which are still unattainable for most widely used methylation profiling methods such as RRBS and MeDIP. We performed extensive analyses to validate the power and accuracy of our method in both model (plant Arabidopsis thaliana) and non-model (scallop Patinopecten yessoensis) species. We further demonstrated its great utility in identification of a gene (LPCAT1) that is potentially crucial for carotenoid accumulation in scallop adductor muscle. MethylRAD has several advantages over existing tools and fills a void in the current epigenomic toolkit by providing a universal tool that can be used for diverse research applications, e.g. from model to non-model species, from ordinary to precious samples and from small to large genomes, but at an affordable cost.

Applying the Power of Transcriptomics: Understanding Male Sexual Development in Decapod Crustacea

The decapod Crustacea are the most species-rich order of the Crustacea and include some of the most charismatic and highly valued commercial species. Thus the decapods draw a significant research interest in relation to aquaculture, as well as gaining a broader understanding of these species' biology. However, the diverse physiology of the group considered with the lack of a model species have presented an obstacle for comparative analyses. In reflection of this, the recent integration of comparative transcriptomics has rapidly advanced our understanding of key regulatory pathways and developmental phenomena, an example being our understanding of sexual development. We discuss our work in the Eastern spiny lobster, Sagmariasus verreauxi, in the context of what is currently known about male sexual development in the decapods, highlighting the importance of transcriptomic techniques in achieving our recent advancements. We describe the progression made in our understanding of male sexual differentiation and maturation, as mediated by the insulin-like androgenic gland hormone (IAG), integrating the role of regulatory binding proteins (IGFBPs), a tyrosine kinase insulin receptor (TKIR), as well as the upstream effect of neuroendocrine hormones (GIH and MIH). We then consider the less well understood mechanism of male sex determination, with an emphasis on what we believe to be the key regulatory factors, the Dsx- and mab-3-related transcription factors (Dmrts). Finally, we discuss the function of the antennal gland (AnG) in sexual development, relating to the emergence of male-biased upregulation in the AnG in later sexual maturation and the sexually dimorphic expression of two key genes Sv-TKIR and Sv-Dmrt1 We then present the AnG as a case study to illustrate how comparative transcriptomic techniques can be applied to guide preliminary analyses, like the hypothesis that the AnG may function in pheromone biosynthesis. In summary, we describe the power of transcriptomics in facilitating the progress made in our understanding of male sexual development, as illustrated by the commercial decapod species, S. verreauxi Considering future directions, we suggest that the integration of multiple omics-based techniques offers the most powerful tool to ensure we continue to piece together the biology of the important group of decapod Crustacea.

Similarities between decapod and insect neuropeptidomes

Background. Neuropeptides are important regulators of physiological processes and behavior. Although they tend to be generally well conserved, recent results using trancriptome sequencing on decapod crustaceans give the impression of significant differences between species, raising the question whether such differences are real or artefacts. Methods. The BLAST+ program was used to find short reads coding neuropeptides and neurohormons in publicly available short read archives. Such reads were then used to find similar reads in the same archives, and the DNA assembly program Trinity was employed to construct contigs encoding the neuropeptide precursors as completely as possible. Results. The seven decapod species analyzed in this fashion, the crabs Eriocheir sinensis, Carcinus maenas and Scylla paramamosain, the shrimp Litopenaeus vannamei, the lobster Homarus americanus, the fresh water prawn Macrobrachium rosenbergii and the crayfish Procambarus clarkii had remarkably similar neuropeptidomes. Although some neuropeptide precursors could not be assembled, in many cases individual reads pertaining to the missing precursors show unambiguously that these neuropeptides are present in these species. In other cases, the tissues that express those neuropeptides were not used in the construction of the cDNA libraries. One novel neuropeptide was identified: elongated PDH (pigment dispersing hormone), a variation on PDH that has a two-amino-acid insertion in its core sequence. Hyrg is another peptide that is ubiquitously present in decapods and is likely a novel neuropeptide precursor. Discussion. Many insect species have lost one or more neuropeptide genes, but apart from elongated PDH and hyrg all other decapod neuropeptides are present in at least some insect species, and allatotropin is the only insect neuropeptide missing from decapods. This strong similarity between insect and decapod neuropeptidomes makes it possible to predict the receptors for decapod neuropeptides that have been deorphanized in insects. This includes the androgenic insulin-like peptide that seems to be homologous to drosophila insulin-like peptide 8.

Evaluation of the 2b-RAD method for genomic selection in scallop breeding

The recently developed 2b-restriction site-associated DNA (2b-RAD) sequencing method provides a cost-effective and flexible genotyping platform for aquaculture species lacking sufficient genomic resources. Here, we evaluated the performance of this method in the genomic selection (GS) of Yesso scallop (Patinopecten yessoensis) through simulation and real data analyses using six statistical models. Our simulation analysis revealed that the prediction accuracies obtained using the 2b-RAD markers were slightly lower than those obtained using all polymorphic loci in the genome. Furthermore, a small subset of markers obtained from a reduced tag representation (RTR) library presented comparable performance to that obtained using all markers, making RTR be an attractive approach for GS purpose. Six GS models exhibited variable performance in prediction accuracy depending on the scenarios (e.g., heritability, sample size, population structure), but Bayes-alphabet and BLUP-based models generally outperformed other models. Finally, we performed the evaluation using an empirical dataset composed of 349 Yesso scallops that were derived from five families. The prediction accuracy for this empirical dataset could reach 0.4 based on optimal GS models. In summary, the genotyping flexibility and cost-effectiveness make 2b-RAD be an ideal genotyping platform for genomic selection in aquaculture breeding programs.

Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing

BACKGROUND

Undaria pinnatifida is an important economic brown alga in East Asian countries. However, its genetic and genomic information is very scarce, which hinders further research in this species. A high-density genetic map is a basic tool for fundamental and applied research such as discovery of functional genes and mapping of quantitative trait loci (QTL). In this study the recently developed specific length amplified fragment sequencing (SLAF-seq) technology was employed to construct a high-density genetic linkage map and locate a sex determining locus for U. pinnatifida.

RESULTS

A total of 28.06 Gb data including 140.31 M pair-end reads was obtained. After linkage analysis 4626 SLAF markers were mapped onto the genetic map. After adding the sex linked simple sequence repeat (SSR) marker [GenBank:AY738602.1], the final genetic map was 1816.28 cM long, consisting of 30 linkage groups with an average distance of 0.39 cM between adjacent markers. The length of LGs ranged from 20.12 to 106.95 cM. A major sex associated QTL was mapped to LG22 within a window starting at 29.01 cM and ending at 68.81 cM with a total of 68 SLAF markers. The SSR marker and five SLAF markers (Marker6556, 19020, 43089, 60771 and 26359) were identified as tightly sex-linked markers, as indicated by the absence of recombination between them and the sex phenotype. These markers were located at the position of 59.50 cM, which was supposed to be the sex determining region.

CONCLUSIONS

A high-density genetic linkage map was constructed using SLAF-seq technique and F1 gametophyte population for the first time in the economically important brown alga U. pinnatifida. For the first time, a major sex associated QTL suggesting a sex determining region was mapped to a single LG. This map will facilitate the further fundamental and applied research such as QTL mapping and map-based gene clone in U. pinnatifida and provide a reference for studies in other kelp species.

Construction of a High-Density Genetic Map and Quantitative Trait Locus Mapping in the Sea Cucumber Apostichopus japonicus

Genetic linkage maps are critical and indispensable tools in a wide range of genetic and genomic research. With the advancement of genotyping-by-sequencing (GBS) methods, the construction of a high-density and high-resolution linkage maps has become achievable in marine organisms lacking sufficient genomic resources, such as echinoderms. In this study, high-density, high-resolution genetic map was constructed for a sea cucumber species, Apostichopus japonicus, utilizing the 2b-restriction site-associated DNA (2b-RAD) method. A total of 7839 markers were anchored to the linkage map with the map coverage of 99.57%, to our knowledge, this is the highest marker density among echinoderm species. QTL mapping and association analysis consistently captured one growth-related QTL located in a 5 cM region of linkage group (LG) 5. An annotated candidate gene, retinoblastoma-binding protein 5 (RbBP5), which has been reported to be an important regulator of cell proliferation, was recognized in the QTL region. This linkage map represents a powerful tool for research involving both fine-scale QTL mapping and marker assisted selection (MAS), and will facilitate chromosome assignment and improve the whole-genome assembly of sea cucumber in the future.

Transcriptomic variation of hepatopancreas reveals the energy metabolism and biological processes associated with molting in Chinese mitten crab, Eriocheir sinensis

Molting is a critical developmental process for crustaceans, yet the underlying molecular mechanism is unknown. In this study, we used RNA-Seq to investigate transcriptomic profiles of the hepatopancreas and identified differentially expressed genes at four molting stages of Chinese mitten crab (Eriocheir sinensis). A total of 97,398 transcripts were assembled, with 31,900 transcripts annotated. Transcriptomic comparison revealed 1,189 genes differentially expressed amongst different molting stages. We observed a pattern associated with energy metabolism and physiological responses during a molting cycle. In specific, differentially expressed genes enriched in postmolt were linked to energy consumption whereas genes enriched in intermolt were related to carbohydrates, lipids metabolic and biosynthetic processes. In premolt, a preparation stage for upcoming molting and energy consumption, highly expressed genes were enriched in response to steroid hormone stimulus and immune system development. The expression profiles of twelve functional genes detected via RNA-Seq were corroborated through real-time RT-PCR assay. Together, our results, including assembled transcriptomes, annotated functional elements and enriched differentially expressed genes amongst different molting stages, provide novel insights into the functions of the hepatopancreas in energy metabolism and biological processes pertaining to molting in crustaceans.

Comparative Transcriptome Analysis Reveals Sex-Biased Gene Expression in Juvenile Chinese Mitten Crab Eriocheir sinensis

Sex-biased genes are considered to account for most of phenotypic differences between males and females. In order to explore the sex-biased gene expression in crab, we performed the whole-body transcriptome analysis in male and female juveniles of the Chinese mitten crab Eriocheir sinensis using next-generation sequencing technology. Of the 23,349 annotated unigenes, 148 were identified as sex-related genes. A total of 29 candidate genes involved in primary sex determination pathways were detected, indicating the sex determination cascade of the mitten crab might be more complex than previously supposed. Differential expression analysis showed 448 differentially expressed genes (DEGs) between the two transcriptomes. Most of DEGs were involved in processes such as metabolism and immunity, and not associated with obvious sexual function. The pathway predominantly enriched for DEGs were related to lysosome, which might reflect the differences in metabolism between males and females. Of the immune DGEs, 18 up-regulated genes in females were humoral immune factors, and eight up-regulated genes in males were pattern recognition receptors, suggesting sex differences of immune defense might exist in the mitten crab. In addition, two reproduction-related genes, vitellogenin and insulin-like androgenic gland factor, were identified to express in both sexes but with significantly higher level in males. Our research provides the first whole-body RNA sequencing of sex-specific transcriptomes for juvenile E. sinensis and will facilitate further studies on molecular mechanisms of crab sexual dimorphism.

Molecular characterization and expression profile of three Fem-1 genes in Eriocheir sinensis provide a new insight into crab sex-determining mechanism

The FEM-1 protein of Caenorhabditis elegans plays a crucial role in the nematode sex-determination pathway. Here, we reported the characterization of three members of Fem-1 gene family in Eriocheir sinensis (designated EsFem-1a, EsFem-1b, and EsFem-1c), which were homologs of the nematode FEM-1 protein. The amino acid sequences of EsFem-1a, EsFem-1b, and EsFem-1c contained eight, nine, and eight ankyrin repeats, respectively. None of the ankyrin repeats had its own specific signature, and the evolution of ankyrin repeat was not completely independent. The predicted three-dimensional structure of EsFem-1 proteins exhibited highly similar superhelical conformation, especially the N-terminal six contiguous ankyrin repeats, which provided a binding surface for the protein-protein interaction. Phylogenetic tree based on the amino acid sequences revealed that EsFem-1a, EsFem-1b, and EsFem-1c were divided into three obvious separated clades. EsFem-1 genes were highly expressed in fertilized egg, 2-4 cell and blastula stage comparing with larval stage (P<0.01), which suggested they might be maternal genes. They also showed a certain degree of sexually dimorphic expression in some tissues. Notably, the highest expression of EsFem-1a was in the hepatopancreas, with EsFem-1b in testes and EsFem-1c in muscle (P<0.05), which indicated their potential role in a broad array of tissues. In addition, the genes initially involved in sex differentiation were not limited to those specifically expressed in the developing gonad. Taken together, these results suggested that EsFem-1 might function in crab early sex determination and late gonad development. The identification of Fem-1 gene family in E. sinensis provides a new insight into crab sex-determination mechanism.