Genome-wide mining,characterization, and development of microsatellite markers in Marsupenaeus japonicus by genome survey sequencing

The kuruma prawn, Marsupenaeus japonicus, is one of the most cultivated and consumed species of shrimp. However, very few molecular genetic/genomic resources are publically available for it. Thus, the characterization and distribution of simple sequence repeats(SSRs) remains ambiguous and the use of SSR markers in genomic studies and marker-assisted selection is limited. The goal of this study is to characterize and develop genome-wide SSR markers in M. japonicus by genome survey sequencing for application in comparative genomics and breeding. A total of 326 945 perfect SSRs were identified, among which dinucleotide repeats were the most frequent class(44.08%), followed by mononucleotides(29.67%), trinucleotides(18.96%), tetranucleotides(5.66%), hexanucleotides(1.07%), and pentanucleotides(0.56%). In total, 151 541 SSR loci primers were successfully designed. A subset of 30 SSR primer pairs were synthesized and tested in 42 individuals from a wild population, of which 27 loci(90.0%) were successfully amplified with specific products and 24(80.0%) were polymorphic. For the amplified polymorphic loci, the alleles ranged from 5 to 17(with an average of 9.63), and the average PIC value was 0.796. A total of 58 256 SSR-containing sequences had significant Gene Ontology annotation; these are good functional molecular marker candidates for association studies and comparative genomic analysis. The newly identified SSRs significantly contribute to the M. japonicus genomic resources and will facilitate a number of genetic and genomic studies, including high density linkage mapping, genome-wide association analysis, marker-aided selection, comparative genomics analysis, population genetics, and evolution.     

Quantitative trait loci detection of <Emphasis Type="Italic">Edwardsiella tarda</Emphasis> resistance in Japanese flounder <Emphasis Type="Italic">Paralichthys olivaceus</Emphasis> using bulked segregant analysis

In recent years, Edwardsiella tarda has become one of the most deadly pathogens of Japanese flounder (Paralichthys olivaceus), causing serious annual losses in commercial production. In contrast to the rapid advances in the aquaculture of P. olivaceus, the study of E. tarda resistance-related markers has lagged behind, hindering the development of a disease-resistant strain. Thus, a marker-trait association analysis was initiated, combining bulked segregant analysis (BSA) and quantitative trait loci (QTL) mapping. Based on 180 microsatellite loci across all chromosomes, 106 individuals from the F1333 (♀: F0768 ×♂: F0915) (Nomenclature rule: F+year+family number) were used to detect simple sequence repeats (SSRs) and QTLs associated with E. tarda resistance. After a genomic scan, three markers (Scaffold 404-21589, Scaffold 404-21594 and Scaffold 270-13812) from the same linkage group (LG)-1 exhibited a significant difference between DNA, pooled/bulked from the resistant and susceptible groups (P <0.001). Therefore, 106 individuals were genotyped using all the SSR markers in LG1 by single marker analysis. Two different analytical models were then employed to detect SSR markers with different levels of significance in LG1, where 17 and 18 SSR markers were identified, respectively. Each model found three resistance-related QTLs by composite interval mapping (CIM). These six QTLs, designated qE1–6, explained 16.0%–89.5% of the phenotypic variance. Two of the QTLs, qE-2 and qE-4, were located at the 66.7 cM region, which was considered a major candidate region for E. tarda resistance. This study will provide valuable data for further investigations of E. tarda resistance genes and facilitate the selective breeding of disease-resistant Japanese flounder in the future.     

Characterization and identification of enzyme-producing microflora isolated from the gut of sea cucumber <Emphasis Type="Italic">Apostichopus japonicus</Emphasis>

Gut microorganisms play an important role in the digestion of their host animals. The purpose of this research was to isolate and assess the enzyme-producing microbes from the Apostichopus japonicus gut. Thirty-nine strains that can produce at least one of the three digestive enzymes (protease, amylase, and cellulase) were qualitatively screened based on their extracellular enzyme-producing abilities. The enzyme-producing strains clustered into eight groups at the genetic similarity level of 100% by analyzing the restriction patterns of 16S rDNA amplified with Mbo I. Phylogenetic analysis revealed that 37 strains belonged to the genus Bacillus and two were members of the genus Virgibacillus. Enzyme-producing capability results indicate that the main enzyme-producing microflora in the A. japonicus gut was Bacillus, which can produce protease, amylase, and cellulase. Virgibacillus, however, can only produce protease. The high enzyme-producing capability of the isolates suggests that the gut microbiota play an important role in the sea cucumber digestive process.     

Change of digestive physiology in sea cucumber <Emphasis Type="Italic">Apostichopus japonicus</Emphasis> (Selenka) induced by corn kernels meal and soybean meal in diets

The present study was conducted to determine the change of digestive physiology in sea cucumber Apostichopus japonicus(Selenka) induced by corn kernels meal and soybean meal in diets. Four experimental diets were tested, in which Sargassum thunbergii was proportionally replaced by the mixture of corn kernels meal and soybean meal. The growth performance, body composition and intestinal digestive enzyme activities in A. japonicus fed these 4 diets were examined. Results showed that the sea cucumber exhibited the maximum growth rate when 20% of S. thunbergii in the diet was replaced by corn kernels meal and soybean meal, while 40% of S. thunbergii in the diet can be replaced by the mixture of corn kernels meal and soybean meal without adversely affecting growth performance of A. japonicus. The activities of intestinal trypsin and amylase in A. japonicus can be significantly altered by corn kernels meal and soybean meal in diets. Trypsin activity in the intestine of A. japonicus significantly increased in the treatment groups compared to the control, suggesting that the supplement of corn kernels meal and soybean meal in the diets might increase the intestinal trypsin activity of A. japonicus. However, amylase activity in the intestine of A. japonicus remarkably decreased with the increasing replacement level of S. thunbergii by the mixture of corn kernels meal and soybean meal, suggesting that supplement of corn kernels meal and soybean meal in the diets might decrease the intestinal amylase activity of A. japonicus.     

Characterization of genome-wide microsatellites of <Emphasis Type="Italic">Saccharina japonica</Emphasis> based on a preliminary assembly of Illumina sequencing reads

Microsatellites or simple sequence repeats (SSR) function widely and locate dependently in genome. However, their characteristics are often ignored due to the lack of genomic sequences of most species. Kelp (Saccharina japonica), a brown macroalga, is extensively cultured in China. In this study, the genome of S. japonica was surveyed using an Illumina sequencing platform, and its microsatellites were characterized. The preliminarily assembled genome was 469.4 Mb in size, with a scaffold N₅₀ of 20529 bp. Among the 128370 identified microsatellites, 90671, 25726 and 11973 were found in intergenic regions, introns and exons, averaging 339.3, 178.8 and 205.4 microsatellites per Mb, respectively. These microsatellites distributed unevenly in S. japonica genome. Mononucleotide motifs were the most abundant in the genome, while trinucleotide ones were the most prevalent in exons. The microsatellite abundance decreased significantly with the increase of motif repeat numbers, and the microsatellites with a small number of repeats accounted for a higher proportion of the exons than those of the intergenic regions and introns. C/G-rich motifs were more common in exons than in intergenic regions and introns. These characteristics of microsatellites in S. japonica genome may associate with their functions, and ultimately their adaptation and evolution. Among the 120140 pairs of designed microsatellite primers, approximately 75% were predicted to be able to amplify S. japonica DNA. These microsatellite markers will be extremely useful for the genetic breeding and population evolution studies of kelp.     

Profile of candidate microsatellite markers in <Emphasis Type="Italic">Sebastiscus marmoratus</Emphasis> using 454 pyrosequencing

Sebastiscus marmoratus is an important sedentary ovoviparous fish distributed in near-shore coastal waters from the coast of China to Japan. Candidate S. marmoratus microsatellite markers were developed in the present study using 454 pyrosequencing, and the marker profile was analyzed. A total of 2 000 000 raw sequence reads were assembled to reduce redundancy. Among them, 1 043 dinucleotide, 925 trinucleotide, 692 tetranucleotide, and 315 pentanucleotide repeats were detected. AC repeats were the most frequent motifs among the dinucleotide repeats, and AAT was the most abundant among the trinucleotide repeats. AAAT, ATAG, and ATCC were the three most common tetranucleotide motifs, and AAGAT and AATAT were the most dominant pentanucleotide motifs. The greatest numbers of loci and potentially amplifiable loci were found in dinucleotide repeats, whereas trinucleotide repeats had the fewest. In summary, a wide range of candidate microsatellite markers were identified in the present study using a rapid and efficient 454 pyrosequencing approach.     

Succession and seasonal variation in epilithic biofilms on artificial reefs in culture waters of the sea cucumber <Emphasis Type="Italic">Apostichopus japonicus</Emphasis>

Periphytic biofilms in aquaculture waters are thought to improve water quality, provide an additional food source, and improve the survival and growth of some reared animals. In the Asia- Pacific region, particularly in China, artificial reefs are commonly used in the commercial farming of sea cucumbers. However, few studies have examined the epilithic biofilms on the artificial reefs. To gain a better understanding of the succession of epilithic biofilms and their ecological processes in sea cucumber culture waters, two experiments were conducted in culture waters of the sea cucumber Apostichopus japonicus in Rongcheng, China, using artificial test panels. On the test panels of succession experiment, more than 67 species were identified in the biofilms. On the test panels of seasonal variation experiment, more than 46 species were recorded in the biofilms. In both experiments, communities of epilithic biofilms were dominated by diatoms, green algae and the annelid Spirorbis sp. In the initial colonization, the dominant diatoms were Cocconeis sp., Amphora spp. and Nitzschia closterium in June, which were succeeded by species of Navicula, Cocconeis and Nitzschia (July to September), and then by Licmophora abbreviata, Nitzschia closterium and Synedra spp. in the following months. A diatom bloom in the autumn and filamentous green algae burst in the summer were also observed. Ecological indices well annotated the succession and seasonal changes in epilithic communities. Multidimensional scaling (MDS) analysis found significant differences in diatom community composition among months and seasons. Fast growth of biofilms was observed in the summer and autumn, whereas the biomass of summer biofilms was largely made up of filamentous green algae. Present results show that the components of epilithic biofilms are mostly optimal foods of A. japonicus, suggesting that biofilms on artificial reefs may contribute important nutritional sources for sea cucumbers during their growth seasons. Future works should include quantitative determination of the contribution of epilithic biofilms to the diet of A. japonicus, potential roles of epilithic biofilms in regulating the water quality of sea cucumber ponds, and the regulation of epilithic biofilms in sea cucumber culture ponds.     

Characterization of bacterial communities associating with larval development of Yesso Scallop (<Emphasis Type="Italic">Patinopecten yessoensisis</Emphasis> Jay, 1857) by high-throughput sequencing

Bacterial community presumably plays an essential role in inhibiting pathogen colonization and maintaining the health of scallop larvae, but limiting data are available for Yesso scallop (Patinopecten yessoensisis Jay, 1857) larval development stages. The aim of this study was to characterize and compare the bacterial communities associating with Yesso scallop larval development at fertilized egg S1, trochophora S2, D-shaped larvae S3, umbo larvae S4, and juvenile scallop S5 stages by Illumina high-throughput sequencing. Genomic DNA was extracted from the larvae and their associating bactera, and a gene segment covering V3-V4 region of 16S rRNA gene was amplified and sequenced using an Illumina Miseq sequencer. Overall, 106760 qualified sequences with an average length of 449 bp were obtained. Sequences were compared with those retrieved from 16S rRNA gene databases, and 4 phyla, 7 classes, 15 orders, 21 families, 31 genera were identified. Proteobacteria was predominant phylum, accounting for more than 99%, at all 5 larval development stages. At genus level, Pseudomonas was dominant at stages S1 (80.60%), S2 (87.77%) and S5 (68.71%), followed by Photobacterium (17.06%) and Aeromonas (1.64%) at stage S1, Serratia (6.94%), Stenotrophomonas (3.08%) and Acinetobacter (1.2%) at stage S2, Shewanella (25.95%) and Pseudoalteromonas (4.57%) at stage S5. Moreover, genus Pseudoalteromonas became dominant at stages S3 (44.85%) and S4 (56.02%), followed by Photobacterium (29.82%), Pseudomonas (11.86%), Aliivibrio (8.60%) and Shewanella (3.39%) at stage S3, Pseudomonas (18.16%), Aliivibrio (14.29%), Shewanella (4.11%), Psychromonas (4.04%) and Psychrobacter (1.81%) at stage S4. From the results, we concluded that the bacterial community changed significantly at different development stages of Yesso Scallop larvae.     

Development and validation of 89 novel expressed sequence tag-derived microsatellite markers in blood clam, <Emphasis Type="Italic">Tegillarca granosa</Emphasis>

Blood clam, Tegillarca granosa, is an important shellfish in Chinese mariculture industry. Investigative research in this species, such as genetic linkage mapping, requires a large panel of molecular markers. In present study, a total of 89 polymorphic microsatellite markers were developed in T. granosa using the sequence database of Life Sciences Technology 454 next generation sequencing technology. All 89 loci were characterized in 20 individual clams from a natural population inhabiting Yueqing Gulf, Zhejiang Province, China. The number of alleles per polymorphic locus varied between 2 and 15, while the observed heterozygosity, expected heterozygosity and polymorphic information content varied between 0.000 and 1.000, 0.102 and 0.921, and 0.048 and 0.886, respectively. Of the 89 loci identified, 32 loci deviated significantly from Hardy-Weinberg equilibrium following Bonferroni correction. Thirty nine markers, which were shown to be polymorphic in a full-sibling family, were tested in Mendelian segregations. As expected, 32 loci were co-dominantly segregated in a Mendelian fashion. These novel developed microsatellite markers represent useful research tools for investigation of population genetic structure and genetic diversity in this species.     

Significant population genetic structure detected in the rock bream <Emphasis Type="Italic">Oplegnathus fasciatus</Emphasis> (Temminck &amp; Schlegel, 1844) inferred from fluorescent-AFLP analysis

Oplegnathus fasciatus (rock bream) is a commercial rocky reef fish species in East Asia that has been considered for aquaculture. We estimated the population genetic diversity and population structure of the species along the coastal waters of China using fluorescent-amplified fragment length polymorphisms technology. Using 53 individuals from three populations and four pairs of selective primers, we amplified 1 264 bands, 98.73% of which were polymorphic. The Zhoushan population showed the highest Nei’s genetic diversity and Shannon genetic diversity. The results of analysis of molecular variance (AMOVA) showed that 59.55% of genetic variation existed among populations and 40.45% occurred within populations, which indicated that a significant population genetic structure existed in the species. The pairwise fixation index F _st ranged from 0.20 to 0.63 and were significant after sequential Bonferroni correction. The topology of an unweighted pair group method with arithmetic mean tree showed two significant genealogical branches corresponding to the sampling locations of North and South China. The AMOVA and STRUCTURE analyses suggested that the O. fasciatus populations examined should comprise two stocks.     

Location of <Emphasis Type="Italic">Vibrio anguillarum</Emphasis> resistance-associated trait loci in half-smooth tongue sole <Emphasis Type="Italic">Cynoglossus semilaevis</Emphasis> at its microsatellite linkage map

A cultured female half-smooth tongue sole (Cynoglossus semilaevis) was crossed with a wild male, yielding the first filial generation of pseudo-testcrossing from which 200 fish were randomly selected to locate the Vibrio anguillarum resistance trait in half-smooth tongue sole at its microsatellite linkage map. In total, 129 microsatellites were arrayed into 18 linkage groups, ≥4 each. The map reconstructed was 852.85 cM in length with an average spacing of 7.68 cM, covering 72.07% of that expected (1 183.35 cM). The V. anguillarum resistance trait was a composite rather than a unit trait, which was tentatively partitioned into Survival time in Hours After V. anguillarum Infection (SHAVI) and Immunity of V. Anguillarum Infection (IVAI). Above a logarithm of the odds (LOD) threshold of 2.5, 18 loci relative to SHAVI and 3 relative to IVAI were identified. The 3 loci relative to IVAI explained 18.78%, 5.87% and 6.50% of the total phenotypic variation in immunity. The microsatellites bounding the 3 quantitative trait loci (QTLs) of IVAI may in future aid to the selection of V. anguillarum-immune half-smooth tongue sole varieties, and facilitate cloning the gene(s) controlling such immunity.     

Analysis of SSRs derived from an EST library of half-smooth tongue sole <Emphasis Type="Italic">Cynoglossus semilaevis</Emphasis>

A complementary DNA (cDNA) library was constructed from half-smooth tongue sole spleen. A long-read expressed sequence tag (EST) database was generated, containing 3100 cDNA clones, of which 220 clones were fully sequenced. A total of 1060 non-redundant simple sequence repeats (SSRs) were obtained from the cDNA library. An average of 5 kb sequence generates 1 SSR in the half-smooth tongue sole spleen cDNA library. The proportion of the SSR unit size was different in the cDNA library. The monomeric repeats (51.4%) are the most abundant class of SSR in the dataset. The dimeric, trimeric, tetrameric and hexameric repeats are represented in decreasing proportions of 27.2%, 16.0%, 2.8% and 1.9%, respectively. The frequency of pentameric repeats was observed the least (only 0.7%). Most of the monomeric and dimeric repeats are distributed in 3′ and 5′ un-translation region. If translation regions are considered merely, trimeric repeats are the highest, accounting for 57% of the total microsatellites.     

Genome size of <Emphasis Type="Italic">Alexandrium catenella</Emphasis> and <Emphasis Type="Italic">Gracilariopsis lemaneiformis</Emphasis> estimated by flow cytometry

Flow cytometry(FCM) technique has been widely applied to estimating the genome size of various higher plants. However, there is few report about its application in algae. In this study, an optimized procedure of FCM was exploited to estimate the genome size of two eukaryotic algae. For analyzing Alexandrium catenella, an important red tide species, the whole cell instead of isolated nucleus was studied, and chicken erythrocytes were used as an internal reference. The genome size of A. catenella was estimated to be 56.48 ± 4.14 Gb(1C), approximately nineteen times larger than that of human genome. For analyzing Gracilariopsis lemaneiformis, an important economical red alga, the purified nucleus was employed, and Arabidopsis thaliana and Chondrus crispus were used as internal references, respectively. The genome size of Gp. lemaneiformis was 97.35 ± 2.58 Mb(1C) and 112.73 ± 14.00 Mb(1C), respectively, depending on the different internal references. The results of this research will promote the related studies on the genomics and evolution of these two species.     

Development of three multiplex PCR primer sets for ark shell (<Emphasis Type="Italic">Scapharca broughtonii</Emphasis>) and their validation in parentage assignment

Scapharca broughtonii is a commercially important and over-exploited species. In order to investigate its genetic diversity and population structure, 43 novel polymorphic microsatellites were isolated and characterized. The number of alleles per locus ranged from 3 to 22 with an average of 6.93, and the observed and expected heterozygosities varied between 0.233 and 1.000, and 0.250 and 0.953, with an average of 0.614 and 0.707, respectively. Three highly informative multiplex PCRs were developed from nine of those microsatellites for S. broughtonii. We evaluated and validated these multiplex PCRs in 8 full-sib families. The average polymorphism information content (PIC) was 0.539. The frequency of null alleles was estimated as 3.13% of all the alleles segregation based on a within-family analysis of Mendelian segregation patterns. Parentage analysis of real offspring demonstrated that 100% of all offspring were unambiguously allocated to a pair of parents based on 3 multiplex sets. Those 43 microsatellite loci with high variability will be helpful for the analysis of population genetics and conservation of wild stock of S. broughtonii. The 3 sets of multiplex PCRs could be an important tool of pedigree reconstruction, population genetic analysis and brood stock management.     

QTL Detection for Albinism-Related Loci in Chinese Tongue Sole (<Emphasis Type="Italic">Cynoglossus semilaevis</Emphasis>)

The Chinese tongue sole (Cynoglossus semilaevis) is widely cultured in the coastal region of East Asia and has excellent economic value. However, the high albino rate of the breeding population has caused a significant loss to the aquaculture industry. To study the molecular mechanism of albinism, the present study used an albino Chinese tongue sole family to construct three simple sequence repeat (SSR) linkage groups, and draft a preliminary linkage map related to albinism. After albinism-related loci mapping, 18 albinism-related loci were detected under two models (containing 2407 genes) compared to the Chinese tongue sole genome. One of these loci, the tyrosinase related protein (tyrp2), which has been reported previously as an important gene regulating both eumelanin and phaeomelanin levels, was indicated to be the possible cause of albinism. Thirty-five Gene Ontology (GO) terms and 14 Kyoto Encyclopedia of Genes and Genome pathways were annotated via bioinformatic analyses. One GO term with protein tyrosine kinase activity, which contained 10 genes, was previously suggested to affect fish albinism. These results establish a foundation for further in-depth study of albinism in Chinese tongue sole.     

Analysis of DNA Methylation of <Emphasis Type="Italic">Gracilariopsis lemaneiformis</Emphasis> Under Temperature Stress Using the Methylation Sensitive Amplification Polymorphism (MSAP) Technique

Gracilariopsis lemaneiformis is an economically important agarophyte, which contains high quality gel and shows a high growth rate. Wild population of G. lemaneiformis displayed resident divergence, though with a low genetic diversity as was revealed by amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) analyses. In addition, different strains of G. lemaneiformis are diverse in morphology. The highly inconsistence between genetic background and physiological characteristics recommends strongly to the regulation at epigenetic level. In this study, the DNA methylation change in G. lemaneiformis among different generation branches and under different temperature stresses was assessed using methylation sensitive amplified polymorphism (MSAP) technique. It was shown that DNA methylation level among different generation branches was diverse. The full and total methylated DNA level was the lowest in the second generation branch and the highest in the third generation. The total methylation level was 61.11%, 60.88% and 64.12% at 15°C, 22°C and 26°C, respectively. Compared with the control group (22°C), the fully methylated and totally methylated ratios were increased in both experiment groups (15°C and 26°C). All of the cytosine methylation/demethylation transform (CMDT) was further analyzed. High temperature treatment could induce more CMDT than low temperature treatment did.