Metabolomics and Molecular Networking-Guided Screening of Bacillus-Derived Bioactive Compounds Against a Highly Lethal Vibrio Species.

Microorganisms are important sources of bioactive natural products. However, the complexity of microbial metabolites and the low abundance of active compounds render the isolation and purification process laborious and inefficient. During our search for active substances capable of inhibiting the newly discovered highly lethal Vibrio strain vp-HL, we found that the fermentation broth of multiple Bacillus strains exhibited antibacterial activity. However, the substances responsible for the activity remained unclear. Metabolomics, molecular networking (MN), and the Structural similarity Network Annotation Platform for Mass Spectrometry (SNAP-MS) were employed in conjunction with bioactivity screening to predict the antibacterial compounds from Bacillus strains. The analysis of fractions, and their isolation, NMR-based annotation, and bioactivity evaluation of an amicoumacin compound partially confirmed the prediction from these statistical analyses. This work presents the potential of marine Bacillus in producing active substances against Vibrio species. Additionally, it highlighted the significance and feasibility of metabolomics and MN in the dereplication of compounds and the determination of isolation targets.

Galbibacter pacificus sp. nov., isolated from surface seawater of the western Pacific Ocean and transfer of Joostella marina to the genus Galbibacter as Galbibacter orientalis nom. nov. and emended description of the genus Galbibacter.

Two Gram-stain-negative, non-motile, non-spore-forming, strictly aerobic and rod-shaped bacterial strains, CMA-7T and CAA-3, were isolated from surface seawater samples collected from the western Pacific Ocean. Phylogeny of 16S rRNA gene sequences indicated they were related to the genera Galbibacter and Joostella and shared 95.1, 90.9 and 90.8% sequence similarity with G. mesophilus Mok-17T, J. marina DSM 19592T and G. marinus ck-I2-15T, respectively. Phylogenomic analysis showed that the two strains, together with the members of the genera Galbibacter and Joostella, formed a monophyletic clade that could also be considered a monophyletic taxon. This distinctiveness was supported by amino acid identity and percentage of conserved proteins indices, phenotypic and chemotaxonomic characteristics and comparative genomics analysis. Digital DNA‒DNA hybridization values and average nucleotide identities between the two strains and their closest relatives were 18.0-20.8 % and 77.7-79.3 %, respectively. The principal fatty acids were iso-C15 : 0, iso-C17 : 0 3-OH, iso-C15 : 1 G, Summed Feature 3 (C16 : 1 ω7c/C16 : 1 ω6c or C16 : 1 ω6c/C16 : 1 ω7c), Summed Feature 9 (iso-C17 : 1 ω9c or C16 : 0 10-methyl), and C15 : 0 3-OH. The predominant respiratory quinone was MK-6. The polar lipids were phosphatidylethanolamine, aminolipid, aminophospholipid, phospholipid, phosphoglycolipid, glycolipid and unknown polar lipid. The genomic DNA G+C content of strains CMA-7T and CAA-3 was both 38.4 mol%. Genomic analysis indicated they have the potential to degrade cellulose and chitin. Based on the polyphasic evidence presented in this study, the two strains represent a novel species within the genus Galbibacter, for which the name Galbibacter pacificus sp. nov. is proposed. The type strain is CMA-7T (=MCCC M28999T = KCTC 92588T). Moreover, the transfer of Joostella marina to the genus Galbibacter as Galbibacter orientalis nom. nov. (type strain En5T = KCTC 12518T = DSM 19592T=CGMCC 1.6973T) is also proposed.

Alcanivorax xiamenensis sp. nov., a novel marine hydrocarbonoclastic bacterium isolated from surface seawater.

A novel Alcanivorax-related strain, designated 6-D-6T, was isolated from the surface seawater collected around Xiamen Island. The novel strain is Gram-stain-negative, rod-shaped and motile, and grows at 10-45 °C, pH 6.0-9.0 and in the presence of 0.5-15.0 % (w/v) NaCl. Phylogenetic analysis based on the 16S rRNA gene sequences indicated that it belongs to the genus Alcanivorax, with the highest sequence similarity to Alcanivorax dieselolei B5T (99.9 %), followed by Alcanivorax xenomutans JC109T (99.5 %), Alcanivorax balearicus MACL04T (99.3 %) and other 13 species of the genus Alcanivorax (93.8 %-95.6 %). The digital DNA-DNA hybridization and average nucleotide identity values between strain 6-D-6T and three close type strains were 40.1-42.9/90.6-91.4 %, and others were below 22.9/85.1 %, respectively. The novel strain contained major cellular fatty acids of C16 : 0 (31.0 %), C19 : 0 ω8c cyclo (23.5 %), C17 : 0 cyclo (9.7 %), C12 : 0 3OH (8.6 %), summed feature 8 (7.6 %) and C12 : 0 (5.4 %). The genomic G+C content of strain 6-D-6T was 61.38 %. Phosphatidylethanolamine, phosphatidylglycerol, diphosphatidylglycerol, two unidentified phospholipids and one amino-group-containing phospholipid were detected. On the basis of phenotypic and genotypic characteristics, strain 6-D-6T represents a novel species within the genus Alcanivorax, for which the name Alcanivorax xiamenensis sp. nov. is proposed. The type strain is 6-D-6T (=MCCC 1A01359T=KCTC 92480T).

Apatinib remodels the immunosuppressive tumor ecosystem of gastric cancer enhancing anti-PD-1 immunotherapy.

Apatinib has been shown to clinically enhance anti-PD-1 immunotherapy for advanced gastric cancer (GC). However, the complexity of GC immunosuppression remains a challenge for precision immunotherapy. Here, we profile the transcriptomes of 34,182 single cells from GC patient-derived xenografts of humanized mouse models treated with vehicle, nivolumab, or nivolumab plus apatinib. Notably, excessive expression of CXCL5 in the CellCycle malignant epithelium, induced by anti-PD-1 immunotherapy and blocked by combined apatinib treatment, is found to be a key driver of tumor-associated neutrophil (TAN) recruitment in the tumor microenvironment through the CXCL5/CXCR2 axis. We further show that the protumor TAN signature is associated with anti-PD-1 immunotherapy-related progressive disease and poor cancer prognosis. Molecular and functional analyses in cell-derived xenograft models confirm the positive in vivo therapeutic effect of targeting the CXCL5/CXCR2 axis during anti-PD-1 immunotherapy. Altogether, our study elucidates the GC immunosuppressive landscape in anti-PD-1 immunotherapy and highlights potential targets for overcoming checkpoint immunotherapy resistance.

Biodegradation of polystyrene (PS) by marine bacteria in mangrove ecosystem.

Plastics pollution poses a new threat to marine ecosystems. Mangrove locating at estuary worldwide is probably the most heavily polluted area trapping various plastics transported from terrestrial and nearby marine aquaculture. Expanded polystyrene (EPS) is one of most common plastic debris therein and even in the plastic garbage. Here we showed the bacterial diversity of the polystyrene-degrading microbial community from EPS waste sites from a subtropical mangrove area. After enrichment with EPS, the degradation consortia were obtained. They shared a similar community structure dominated by bacteria of Sphingomonadaceae, Rhodanobacteraceae, Rhizobiaceae, Dermacoccaceae, Rhodocyclaceae, Hyphomicrobiaceae, and Methyloligellaceae. Diverse bacteria standing for the first member of the genera of Novosphingobium, Gordonia, Stappia, Mesobacillus, Alcanivorax, Flexivirga, Cytobacillus, Thioclava, and Thalassospira showed PS degradation capability as a pure culture. Further, PS biodegradation of Gordonia sp. and Novosphingobium sp. was quantified by weight loss, in addition to obvious morphological and structural changes of the PS films observed by SEM, ATR-FTIR, and contact angle analysis. The formation of new oxygen-containing functional groups implied the degradation pathway of oxidation. Although the degradation rates ranged from 2.7% to 7.7% after one month in lab and possibly lower in situ, their role in EPS removal is unneglectable.

Characterization of a bioactive meroterpenoid isolated from the marine-derived fungus Talaromyces sp.

A new meroterpenoid, taladrimanin A (1), was isolated from a marine-derived fungus Talaromyces sp. HM6-1-1, together with eleven biogenetically related compounds (2-12). A plausible biosynthetic pathway for the meroterpenoids (1-4) was proposed. The planar structure of 1 was assigned by HRESIMS and NMR. Its relative configuration was established by quantum chemical NMR calculation of two possible isomers and analyzed by DP4 + method. Finally, X-ray diffraction unambiguously confirmed the relative configuration and revealed the absolute configuration of compound 1. 2-12 were assigned by comparing their NMR data with those reported in the literature. 1 was the first drimane-type meroterpenoid with a C10 polyketide unit bearing an 8R-configuration. In the bioactive assay, 1 exhibited antitumor activity against gastric cancer cells MGC803 and MKN28; it also inhibited the colony formation and induced apoptosis in MGC803 cells both in a concentration-dependent manner. Additionally, 1 displayed selective antibacterial activity against Staphylococcus aureus 6538P, and low activities towards strains of Vibrio parahaemolyticus and Escherichia coli in this study. KEY POINTS: • Twelve compounds were obtained from Talaromyces sp., including four meroterpenoids, one of which was new. • The new compound taladrimanin A (1) inhibits the growth of gastric cancer cells MGC803 and MKN28 as well as the pathogenic bacteria Staphylococcus aureus 6538P. • The biosynthetic pathway of the meroterpenoids was proposed.

Bacillus pumilus Group Comparative Genomics: Toward Pangenome Features, Diversity, and Marine Environmental Adaptation.

BACKGROUND: Members of the Bacillus pumilus group (abbreviated as the Bp group) are quite diverse and ubiquitous in marine environments, but little is known about correlation with their terrestrial counterparts. In this study, 16 marine strains that we had isolated before were sequenced and comparative genome analyses were performed with a total of 52 Bp group strains. The analyses included 20 marine isolates (which included the 16 new strains) and 32 terrestrial isolates, and their evolutionary relationships, differentiation, and environmental adaptation. RESULTS: Phylogenomic analysis revealed that the marine Bp group strains were grouped into three species: B. pumilus, B. altitudinis and B. safensis. All the three share a common ancestor. However, members of B. altitudinis were observed to cluster independently, separating from the other two, thus diverging from the others. Consistent with the universal nature of genes involved in the functioning of the translational machinery, the genes related to translation were enriched in the core genome. Functional genomic analyses revealed that the marine-derived and the terrestrial strains showed differences in certain hypothetical proteins, transcriptional regulators, K+ transporter (TrK) and ABC transporters. However, species differences showed the precedence of environmental adaptation discrepancies. In each species, land specific genes were found with possible functions that likely facilitate survival in diverse terrestrial niches, while marine bacteria were enriched with genes of unknown functions and those related to transcription, phage defense, DNA recombination and repair. CONCLUSION: Our results indicated that the Bp isolates show distinct genomic features even as they share a common core. The marine and land isolates did not evolve independently; the transition between marine and non-marine habitats might have occurred multiple times. The lineage exhibited a priority effect over the niche in driving their dispersal. Certain intra-species niche specific genes could be related to a strains adaptation to its respective marine or terrestrial environment(s). In summary, this report describes the systematic evolution of 52 Bp group strains and will facilitate future studies toward understanding their ecological role and adaptation to marine and/or terrestrial environments.

Characterization of Sulfurimonas hydrogeniphila sp. nov., a Novel Bacterium Predominant in Deep-Sea Hydrothermal Vents and Comparative Genomic Analyses of the Genus Sulfurimonas.

Bacteria of the genus Sulfurimonas within the class Campylobacteria are predominant in global deep-sea hydrothermal environments and widespread in global oceans. However, only few bacteria of this group have been isolated, and their adaptations for these extreme environments remain poorly understood. Here, we report a novel mesophilic, hydrogen- and sulfur-oxidizing bacterium, strain NW10T, isolated from a deep-sea sulfide chimney of Northwest Indian Ocean.16S rRNA gene sequence analysis showed that strain NW10T was most closely related to the vent species Sulfurimonas paralvinellae GO25T with 95.8% similarity, but ANI and DDH values between two strains were only 19.20 and 24.70%, respectively, indicating that strain NW10 represents a novel species. Phenotypic characterization showed strain NW10T is an obligate chemolithoautotroph utilizing thiosulfate, sulfide, elemental sulfur, or molecular hydrogen as energy sources, and molecular oxygen, nitrate, or elemental sulfur as electron acceptors. Moreover, hydrogen supported a better growth than reduced sulfur compounds. During thiosulfate oxidation, the strain can produce extracellular sulfur of elemental α-S8 with an unknown mechanism. Polyphasic taxonomy results support that strain NW10T represents a novel species of the genus Sulfurimonas, and named as Sulfurimonas hydrogeniphila sp. nov. Genome analyses revealed its diverse energy metabolisms driving carbon fixation via rTCA cycling, including pathways of sulfur/hydrogen oxidation, coupled oxygen/sulfur respiration and denitrification. Comparative analysis of the 11 available genomes from Sulfurimonas species revealed that vent bacteria, compared to marine non-vent strains, possess unique genes encoding Type V Sqr, Group II, and Coo hydrogenase, and are selectively enriched in genes related to signal transduction and inorganic ion transporters. These phenotypic and genotypic features of vent Sulfurimonas may explain their thriving in hydrothermal environments and help to understand the ecological role of Sulfurimonas bacteria in hydrothermal ecosystems.

Genome sequencing of deep-sea hydrothermal vent snails reveals adaptions to extreme environments.

BACKGROUND: The scaly-foot snail (Chrysomallon squamiferum) is highly adapted to deep-sea hydrothermal vents and has drawn much interest since its discovery. However, the limited information on its genome has impeded further related research and understanding of its adaptation to deep-sea hydrothermal vents. FINDINGS: Here, we report the whole-genome sequencing and assembly of the scaly-foot snail and another snail (Gigantopelta aegis), which inhabits similar environments. Using Oxford Nanopore Technology, 10X Genomics, and Hi-C technologies, we obtained a chromosome-level genome of C. squamiferum with an N50 size of 20.71 Mb. By constructing a phylogenetic tree, we found that these 2 deep-sea snails evolved independently of other snails. Their divergence from each other occurred ∼66.3 million years ago. Comparative genomic analysis showed that different snails have diverse genome sizes and repeat contents. Deep-sea snails have more DNA transposons and long terminal repeats but fewer long interspersed nuclear elements than other snails. Gene family analysis revealed that deep-sea snails experienced stronger selective pressures than freshwater snails, and gene families related to the nervous system, immune system, metabolism, DNA stability, antioxidation, and biomineralization were significantly expanded in scaly-foot snails. We also found 251 H-2 Class II histocompatibility antigen, A-U α chain-like (H2-Aal) genes, which exist uniquely in the Gigantopelta aegis genome. This finding is important for investigating the evolution of major histocompatibility complex (MHC) genes. CONCLUSION: Our study provides new insights into deep-sea snail genomes and valuable resources for further studies.

Identification of a Tumor Microenvironment-relevant Gene set-based Prognostic Signature and Related Therapy Targets in Gastric Cancer.

Rationale: The prognosis of gastric cancer (GC) patients is poor, and there is limited therapeutic efficacy due to genetic heterogeneity and difficulty in early-stage screening. Here, we developed and validated an individualized gene set-based prognostic signature for gastric cancer (GPSGC) and further explored survival-related regulatory mechanisms as well as therapeutic targets in GC. Methods: By implementing machine learning, a prognostic model was established based on gastric cancer gene expression datasets from 1699 patients from five independent cohorts with reported full clinical annotations. Analysis of the tumor microenvironment, including stromal and immune subcomponents, cell types, panimmune gene sets, and immunomodulatory genes, was carried out in 834 GC patients from three independent cohorts to explore regulatory survival mechanisms and therapeutic targets related to the GPSGC. To prove the stability and reliability of the GPSGC model and therapeutic targets, multiplex fluorescent immunohistochemistry was conducted with tissue microarrays representing 186 GC patients. Based on multivariate Cox analysis, a nomogram that integrated the GPSGC and other clinical risk factors was constructed with two training cohorts and was verified by two validation cohorts. Results: Through machine learning, we obtained an optimal risk assessment model, the GPSGC, which showed higher accuracy in predicting survival than individual prognostic factors. The impact of the GPSGC score on poor survival of GC patients was probably correlated with the remodeling of stromal components in the tumor microenvironment. Specifically, TGFß and angiogenesis-related gene sets were significantly associated with the GPSGC risk score and poor outcome. Immunomodulatory gene analysis combined with experimental verification further revealed that TGFß1 and VEGFB may be developed as potential therapeutic targets of GC patients with poor prognosis according to the GPSGC. Furthermore, we developed a nomogram based on the GPSGC and other clinical variables to predict the 3-year and 5-year overall survival for GC patients, which showed improved prognostic accuracy than clinical characteristics only. Conclusion: As a tumor microenvironment-relevant gene set-based prognostic signature, the GPSGC model provides an effective approach to evaluate GC patient survival outcomes and may prolong overall survival by enabling the selection of individualized targeted therapy.

Scallop genome reveals molecular adaptations to semi-sessile life and neurotoxins.

Bivalve molluscs are descendants of an early-Cambrian lineage superbly adapted to benthic filter feeding. Adaptations in form and behavior are well recognized, but the underlying molecular mechanisms are largely unknown. Here, we investigate the genome, various transcriptomes, and proteomes of the scallop Chlamys farreri, a semi-sessile bivalve with well-developed adductor muscle, sophisticated eyes, and remarkable neurotoxin resistance. The scallop's large striated muscle is energy-dynamic but not fully differentiated from smooth muscle. Its eyes are supported by highly diverse, intronless opsins expanded by retroposition for broadened spectral sensitivity. Rapid byssal secretion is enabled by a specialized foot and multiple proteins including expanded tyrosinases. The scallop uses hepatopancreas to accumulate neurotoxins and kidney to transform to high-toxicity forms through expanded sulfotransferases, probably as deterrence against predation, while it achieves neurotoxin resistance through point mutations in sodium channels. These findings suggest that expansion and mutation of those genes may have profound effects on scallop's phenotype and adaptation.

Scallop genome provides insights into evolution of bilaterian karyotype and development.

Reconstructing the genomes of bilaterian ancestors is central to our understanding of animal evolution, where knowledge from ancient and/or slow-evolving bilaterian lineages is critical. Here we report a high-quality, chromosome-anchored reference genome for the scallop Patinopecten yessoensis, a bivalve mollusc that has a slow-evolving genome with many ancestral features. Chromosome-based macrosynteny analysis reveals a striking correspondence between the 19 scallop chromosomes and the 17 presumed ancestral bilaterian linkage groups at a level of conservation previously unseen, suggesting that the scallop may have a karyotype close to that of the bilaterian ancestor. Scallop Hox gene expression follows a new mode of subcluster temporal co-linearity that is possibly ancestral and may provide great potential in supporting diverse bilaterian body plans. Transcriptome analysis of scallop mantle eyes finds unexpected diversity in phototransduction cascades and a potentially ancient Pax2/5/8-dependent pathway for noncephalic eyes. The outstanding preservation of ancestral karyotype and developmental control makes the scallop genome a valuable resource for understanding early bilaterian evolution and biology.

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.

An integrated genetic and cytogenetic map for Zhikong scallop, Chlamys farreri, based on microsatellite markers.

The reliability of genome analysis and proficiency of genetic manipulation requires knowledge of the correspondence between the genetic and cytogenetic maps. In the present study, we integrated cytogenetic and microsatellite-based linkage maps for Zhikong scallop, Chlamys farreri. Thirty-eight marker-anchored BAC clones standing for the 19 linkage groups were used to be FISH probes. Of 38 BAC clones, 30 were successfully located on single chromosome by FISH and used to integrate the genetic and cytogenetic map. Among the 19 linkage groups, 12 linkage groups were physically anchored by 2 markers, 6 linkage groups were anchored by 1 marker, and one linkage group was not anchored any makers by FISH. In addition, using two-color FISH, six linkage groups were distinguished by different chromosomal location; linkage groups LG6 and LG16 were placed on chromosome 10, LG8 and LG18 on chromosome 14. As a result, 18 of 19 linkage groups were localized to 17 pairs of chromosomes of C. farreri. We first integrated genetic and cytogenetic map for C. farreri. These 30 chromosome specific BAC clones in the cytogenetic map could be used to identify chromosomes of C. farreri. The integrated map will greatly facilitate molecular genetic studies that will be helpful for breeding applications in C. farreri and the upcoming genome projects of this species.

Genome-wide analysis of DNA methylation in five tissues of Zhikong scallop, Chlamys farreri.

DNA methylation plays a vital role in tissue development and differentiation in eukaryotes. Epigenetic studies have been seldom conducted in the extremely diverse and evolutionarily highly successful bilaterian lineage Mollusca. In the present study, we conducted the genome-wide profiling of DNA methylation for five tissues of a bivalve mollusc, Chlamys farreri using the methylation-sensitive amplification polymorphism (MSAP) technique. The methylation levels were quite similar among tissues, ranging from 20.9% to 21.7%. CG methylation was the dominant type (14.9%-16.5%) in the C. farreri genome, but CHG methylation also accounted for a substantial fraction of total methylation (5.1%-6.3%). Relatively high methylation diversity was observed within tissues. Methylation differentiation between tissues was evaluated and 460 tissue-specific epiloci were identified. Kidney differs from the other tissues in DNA methylation profiles. Our study presents the first look at the tissue-specific DNA methylation patterns in a bivalve mollusc and represents an initial step towards understanding of epigenetic regulatory mechanism underlying tissue development and differentiation in bivalves.

Identification of two secreted ferritin subunits involved in immune defense of Yesso scallop Patinopecten yessoensis.

As an important iron storage protein, ferritin plays a crucial role in the iron-withholding defense system. In this study, two secreted ferritin subunits (PyFerS1 and PyFerS2) were identified from the Yesso scallop, Patinopecten yessoensis. The complete DNA sequences of the two ferritins are 7101 and 5359 bp, consisting of seven and five exons, respectively. The full-length cDNAs of PyFerS1 and PyFerS2 are 960 and 956 bp in length, encoding 228 and 220 amino acids, respectively. They have typical ferritin structures, with four long α-helices, one short α-helix and an L-loop. Signal peptides were found at the N-terminus of both ferritins, and phylogenetic analysis showed that they both clustered with secreted mollusc ferritins. PyFerS1 possesses all seven conserved residues of the ferroxidase center, whereas PyFerS2 only has two. Real-time PCR analysis indicated high expression level of PyFerS2 in the D-shaped larvae, and PyFerS1 in both D-shaped larvae and fertilized eggs. In adult scallops, PyFerS1 was only detected in the hepatopancreas, whereas PyFerS2 was detected in both hepatopancreas and mantle. After the scallops were challenged by iron ion or bacteria Vibrio anguillarum, the expression of both PyFerS1 and PyFerS2 was significantly elevated, suggesting they may play a role in scallop innate immune defense. For the first time, secreted ferritins were cloned and comprehensively characterized in bivalve molluscs. It will assist in better understanding of the role of secreted ferritins in bivalve innate immunity.

Large-scale development of gene-associated single-nucleotide polymorphism markers for molluscan population genomic, comparative genomic, and genome-wide association studies.

Mollusca is the second most diverse group of animals in the world. Despite their perceived importance, omics-level studies have seldom been applied to this group of animals largely due to a paucity of genomic resources. Here, we report the first large-scale gene-associated marker development and evaluation for a bivalve mollusc, Chlamys farreri. More than 21,000 putative single-nucleotide polymorphisms (SNPs) were identified from the C. farreri transcriptome. Primers and probes were designed and synthesized for 4500 SNPs, and 1492 polymorphic markers were successfully developed using a high-resolution melting genotyping platform. These markers are particularly suitable for population genomic analysis due to high polymorphism within and across populations, a low frequency of null alleles, and conformation to neutral expectations. Unexpectedly, high cross-species transferability was observed, suggesting that the transferable SNPs may largely represent ancestral genetic variations that have been preserved differentially among subfamilies of Pectinidae. Gene annotations were available for 73% of the markers, and 65% could be anchored to the recently released Pacific oyster genome. Large-scale association analysis revealed key candidate genes responsible for scallop growth regulation, and provided markers for further genetic improvement of C. farreri in breeding programmes.

High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc.

Genetic linkage maps are indispensable tools in genetic and genomic studies. Recent development of genotyping-by-sequencing (GBS) methods holds great promise for constructing high-resolution linkage maps in organisms lacking extensive genomic resources. In the present study, linkage mapping was conducted for a bivalve mollusc (Chlamys farreri) using a newly developed GBS method-2b-restriction site-associated DNA (2b-RAD). Genome survey sequencing was performed to generate a preliminary reference genome that was utilized to facilitate linkage and quantitative trait locus (QTL) mapping in C. farreri. A high-resolution linkage map was constructed with a marker density (3806) that has, to our knowledge, never been achieved in any other molluscs. The linkage map covered nearly the whole genome (99.5%) with a resolution of 0.41 cM. QTL mapping and association analysis congruously revealed two growth-related QTLs and one potential sex-determination region. An important candidate QTL gene named PROP1, which functions in the regulation of growth hormone production in vertebrates, was identified from the growth-related QTL region detected on the linkage group LG3. We demonstrate that this linkage map can serve as an important platform for improving genome assembly and unifying multiple genomic resources. Our study, therefore, exemplifies how to build up an integrative genomic framework in a non-model organism.

RADtyping: an integrated package for accurate de novo codominant and dominant RAD genotyping in mapping populations.

Genetic linkage maps are indispensable tools in genetic, genomic and breeding studies. As one of genotyping-by-sequencing methods, RAD-Seq (restriction-site associated DNA sequencing) has gained particular popularity for construction of high-density linkage maps. Current RAD analytical tools are being predominantly used for typing codominant markers. However, no genotyping algorithm has been developed for dominant markers (resulting from recognition site disruption). Given their abundance in eukaryotic genomes, utilization of dominant markers would greatly diminish the extensive sequencing effort required for large-scale marker development. In this study, we established, for the first time, a novel statistical framework for de novo dominant genotyping in mapping populations. An integrated package called RADtyping was developed by incorporating both de novo codominant and dominant genotyping algorithms. We demonstrated the superb performance of RADtyping in achieving remarkably high genotyping accuracy based on simulated and real mapping datasets. The RADtyping package is freely available at http://www2.ouc.edu.cn/mollusk/ detailen.asp?id=727.

Development of a rapid and efficient method for non-lethal DNA sampling and genotyping in scallops.

Non-lethal DNA sampling has long appealed to researchers studying population and conservation genetics, as it does not necessitate removing individuals permanently from their natural environment or destroying valuable samples. However, such an approach has not yet been well established in bivalves. In this study, we demonstrate that the gill represents a good source of tissue for non-lethal sampling in scallops. Removal of a few gill filaments caused no noticeable behavioral abnormalities or increased mortality rates in Zhikong scallop (Chlamys farreri) during a three-month period of observation. To facilitate rapid gill-based DNA extraction, six methods (MA-MF) were designed and evaluated, each requiring less than one hour of processing time. The optimal method was identified as MF, in terms of maintaining DNA integrity and genotyping accuracy. Further optimization of MF method by orthogonal experimental design suggested that the utilization of gills could be limited to 2 mg of sample, which is sufficient for performing up to 20,000 PCR reactions. We also demonstrate the excellent cross-species utility of MF in two additional scallop species, Yesso scallop (Patinopecten yessoensis) and bay scallop (Argopecten irradians). Taken together, our study provides a rapid and efficient approach for applying non-lethal DNA sampling in bivalve species, which would serve as a valuable tool for maintaining bivalve populations and conservation genetics, as well as in breeding studies.