Menu
July 7, 2019

Genome sequence of the lager-brewing yeast Saccharomyces sp. strain M14, used in the high-gravity brewing industry in China.

Lager-brewing yeasts are mainly used for the production of lager beers. Illumina and PacBio-based sequence analyses revealed an approximate genome size of 22.8 Mb, with a GC content of 38.98%, for the Chinese lager-brewing yeast Saccharomyces sp. strain M14. Based on ab initio prediction, 9,970 coding genes were annotated. Copyright © 2017 Liu et al.


July 7, 2019

Draft genome sequence of the fruiting myxobacterium Nannocystis exedens DSM 71.

In response to starvation, members of the order Myxococcales form morphologically very different fruiting bodies. To determine whether fruiting myxobacteria share a common genetic program that leads to fruiting body formation, we sequenced and assembled the genome of Nannocystis exedens DSM 71 as two contigs with a total GC content of 72%. Copyright © 2017 Treuner-Lange et al.


July 7, 2019

Bioinformatics analysis and characterization of highly efficient polyvinyl alcohol (PVA)-degrading enzymes from the novel PVA degrader Stenotrophomonas rhizophila QL-P4.

Polyvinyl alcohol (PVA) is used widely in industry, and associated environmental pollution is a serious problem. Herein, we report a novel, efficient PVA degrader, Stenotrophomonas rhizophila QL-P4, isolated from fallen leaves from virgin forest in the Qinling Mountains. The complete genome was obtained using single-molecule real-time (SMRT) technology and corrected using Illumina sequencing. Bioinformatics analysis revealed eight PVA/OVA (vinyl alcohol oligomer)-degrading genes. Of these, seven genes were predicted to be involved in the classical intracellular PVA/OVA degradation pathway, and one (BAY15_3292) was identified as a novel PVA oxidase. Five PVA/OVA-degrading enzymes were purified and characterised. Among which, BAY15_1712, a PVA dehydrogenase (PVADH), displayed high catalytic efficiency towards PVA and OVA substrate. All reported PVADHs only have PVA-degrading ability. Most importantly, we discovered a novel PVA oxidase (BAY15_3292) that exhibited highest PVA-degrading efficiency than the reported PVADHs. Further investigation indicated that BAY15_3292 plays a crucial role in PVA degradation in S. rhizophila QL-P4. Knocking out BAY15_3292 resulted in a significant decline in PVA-degrading activity in S. rhizophila QL-P4. Interestingly, we found that BAY15_3292 possesses exocrine activity, which distinguishes it from classical PVADHs. Transparent circle experiments further proved that BAY15_3292 greatly affects extracellular PVA degradation in S. rhizophila QL-P4. The exocrine characteristics of BAY15_3292 facilitate its potential application to PVA bioremediation. In addition, we report three new efficient secondary alcohol dehydrogenases (SADHs) with OVA-degrading ability in S. rhizophila QL-P4, compared with only one OVA-degrading SADH as reported previously.Importance With the widespread application of PVA in industry, PVA-related environmental pollution is an increasingly serious issue. Because PVA is difficult to degrade, it accumulates in aquatic environments and causes chronic toxicity to aquatic organisms. Biodegradation of PVA, as an economical and environment-friendly method, has attracted much interest. To date, effective and applicable PVA-degrading bacteria/enzymes have not been reported. Herein, we report a new efficient PVA degrader (S. rhizophila QL-P4) that has five PVA/OVA-degrading enzymes with high catalytic efficiency, among which BAY15_1712 is the only reported PVADH with both PVA- and OVA-degrading abilities. Importantly, we discovered a novel PVA oxidase (BAY15_3292) that is not only more efficient than other reported PVA-degrading PVADHs, but also has exocrine activity. Overall, our findings provide new insight into PVA-degrading pathways in microorganisms, and suggest S. rhizophila QL-P4 and its enzymes have potential for application to PVA bioremediation to reduce or eliminate PVA-related environmental pollution. Copyright © 2017 American Society for Microbiology.


July 7, 2019

Heat resistance mediated by pLM58 plasmid-borne ClpL in Listeria monocytogenes.

Listeria monocytogenes is one of the most heat-resistant non-spore-forming food-borne pathogens and poses a notable risk to food safety, particularly when mild heat treatments are used in food processing and preparation. While general heat stress properties and response mechanisms of L. monocytogenes have been described, accessory mechanisms providing particular L. monocytogenes strains with the advantage of enhanced heat resistance are unknown. Here, we report plasmid-mediated heat resistance of L. monocytogenes for the first time. This resistance is mediated by the ATP-dependent protease ClpL. We tested the survival of two wild-type L. monocytogenes strains-both of serotype 1/2c, sequence type ST9, and high sequence identity-at high temperatures and compared their genome composition in order to identify genetic mechanisms involved in their heat survival phenotype. L. monocytogenes AT3E was more heat resistant (0.0 CFU/ml log10 reduction) than strain AL4E (1.4 CFU/ml log10 reduction) after heating at 55°C for 40 min. A prominent difference in the genome compositions of the two strains was a 58-kb plasmid (pLM58) harbored by the heat-resistant AT3E strain, suggesting plasmid-mediated heat resistance. Indeed, plasmid curing resulted in significantly decreased heat resistance (1.1 CFU/ml log10 reduction) at 55°C. pLM58 harbored a 2,115-bp open reading frame annotated as an ATP-dependent protease (ClpL)-encoding clpL gene. Introducing the clpL gene into a natively heat-sensitive L. monocytogenes strain (1.2 CFU/ml log10 reduction) significantly increased the heat resistance of the recipient strain (0.4 CFU/ml log10 reduction) at 55°C. Plasmid-borne ClpL is thus a potential predictor of elevated heat resistance in L. monocytogenes. IMPORTANCEListeria monocytogenes is a dangerous food pathogen causing the severe illness listeriosis that has a high mortality rate in immunocompromised individuals. Although destroyed by pasteurization, L. monocytogenes is among the most heat-resistant non-spore-forming bacteria. This poses a risk to food safety, as listeriosis is commonly associated with ready-to-eat foods that are consumed without thorough heating. However, L. monocytogenes strains differ in their ability to survive high temperatures, and comprehensive understanding of the genetic mechanisms underlying these differences is still limited. Whole-genome-sequence analysis and phenotypic characterization allowed us to identify a novel plasmid, designated pLM58, and a plasmid-borne ATP-dependent protease (ClpL), which mediated heat resistance in L. monocytogenes. As the first report on plasmid-mediated heat resistance in L. monocytogenes, our study sheds light on the accessory genetic mechanisms rendering certain L. monocytogenes strains particularly capable of surviving high temperatures-with plasmid-borne ClpL being a potential predictor of elevated heat resistance.


July 7, 2019

FKBP12-dependent inhibition of calcineurin mediates immunosuppressive antifungal drug action in Malassezia.

The genus Malassezia includes yeasts that are commonly found on the skin or hair of animals and humans as commensals and are associated with a number of skin disorders. We have previously developed an Agrobacterium tumefaciens transformation system effective for both targeted gene deletion and insertional mutagenesis in Malassezia furfur and M. sympodialis In the present study, these molecular resources were applied to characterize the immunophilin FKBP12 as the target of tacrolimus (FK506), ascomycin, and pimecrolimus, which are calcineurin inhibitors that are used as alternatives to corticosteroids in the treatment of inflammatory skin disorders such as those associated with Malassezia species. While M. furfur and M. sympodialis showed in vitro sensitivity to these agents, fkb1? mutants displayed full resistance to all three of them, confirming that FKBP12 is the target of these calcineurin inhibitors and is essential for their activity. We found that calcineurin inhibitors act additively with fluconazole through an FKBP12-dependent mechanism. Spontaneous M. sympodialis isolates resistant to calcineurin inhibitors had mutations in the gene encoding FKBP12 in regions predicted to affect the interactions between FKBP12 and FK506 based on structural modeling. Due to the presence of homopolymer nucleotide repeats in the gene encoding FKBP12, an msh2? hypermutator of M. sympodialis was engineered and exhibited an increase of more than 20-fold in the rate of emergence of resistance to FK506 compared to that of the wild-type strain, with the majority of the mutations found in these repeats.IMPORTANCEMalassezia species are the most abundant fungal components of the mammalian and human skin microbiome. Although they belong to the natural skin commensal flora of humans, they are also associated with a variety of clinical skin disorders. The standard treatment for Malassezia-associated inflammatory skin infections is topical corticosteroids, although their use has adverse side effects and is not recommended for long treatment periods. Calcineurin inhibitors have been proposed as a suitable alternative to treat patients affected by skin lesions caused by Malassezia Although calcineurin inhibitors are well-known as immunosuppressive drugs, they are also characterized by potent antimicrobial activity. In the present study, we investigated the mechanism of action of FK506 (tacrolimus), ascomycin (FK520), and pimecrolimus in M. furfur and M. sympodialis and found that the conserved immunophilin FKBP12 is the target of these drugs with which it forms a complex that directly binds calcineurin and inhibits its signaling activity. We found that FKBP12 is also required for the additive activity of calcineurin inhibitors with fluconazole. Furthermore, the increasing natural occurrence in fungal pathogen populations of mutator strains poses a high risk for the rapid emergence of drug resistance and adaptation to host defense. This led us to generate an engineered hypermutator msh2? mutant strain of M. sympodialis and genetically evaluate mutational events resulting in a substantially increased rate of resistance to FK506 compared to that of the wild type. Our study paves the way for the novel clinical use of calcineurin inhibitors with lower immunosuppressive activity that could be used clinically to treat a broad range of fungal infections, including skin disorders caused by Malassezia. Copyright © 2017 Ianiri et al.


July 7, 2019

Genome-wide discovery of genes required for capsule production by uropathogenic Escherichia coli.

Uropathogenic Escherichia coli (UPEC) is a major cause of urinary tract and bloodstream infections and possesses an array of virulence factors for colonization, survival, and persistence. One such factor is the polysaccharide K capsule. Among the different K capsule types, the K1 serotype is strongly associated with UPEC infection. In this study, we completely sequenced the K1 UPEC urosepsis strain PA45B and employed a novel combination of a lytic K1 capsule-specific phage, saturated Tn5 transposon mutagenesis, and high-throughput transposon-directed insertion site sequencing (TraDIS) to identify the complement of genes required for capsule production. Our analysis identified known genes involved in capsule biosynthesis, as well as two additional regulatory genes (mprA and lrhA) that we characterized at the molecular level. Mutation of mprA resulted in protection against K1 phage-mediated killing, a phenotype restored by complementation. We also identified a significantly increased unidirectional Tn5 insertion frequency upstream of the lrhA gene and showed that strong expression of LrhA induced by a constitutive Pcl promoter led to loss of capsule production. Further analysis revealed loss of MprA or overexpression of LrhA affected the transcription of capsule biosynthesis genes in PA45B and increased sensitivity to killing in whole blood. Similar phenotypes were also observed in UPEC strains UTI89 (K1) and CFT073 (K2), demonstrating that the effects were neither strain nor capsule type specific. Overall, this study defined the genome of a UPEC urosepsis isolate and identified and characterized two new regulatory factors that affect UPEC capsule production.IMPORTANCE Urinary tract infections (UTIs) are among the most common bacterial infections in humans and are primarily caused by uropathogenic Escherichia coli (UPEC). Many UPEC strains express a polysaccharide K capsule that provides protection against host innate immune factors and contributes to survival and persistence during infection. The K1 serotype is one example of a polysaccharide capsule type and is strongly associated with UPEC strains that cause UTIs, bloodstream infections, and meningitis. The number of UTIs caused by antibiotic-resistant UPEC is steadily increasing, highlighting the need to better understand factors (e.g., the capsule) that contribute to UPEC pathogenesis. This study describes the original and novel application of lytic capsule-specific phage killing, saturated Tn5 transposon mutagenesis, and high-throughput transposon-directed insertion site sequencing to define the entire complement of genes required for capsule production in UPEC. Our comprehensive approach uncovered new genes involved in the regulation of this key virulence determinant. Copyright © 2017 Goh et al.


July 7, 2019

Insights from the genome sequence of Mycobacterium lepraemurium: Massive gene decay and reductive evolution.

Mycobacterium lepraemurium is the causative agent of murine leprosy, a chronic, granulomatous disease similar to human leprosy. Due to the similar clinical manifestations of human and murine leprosy and the difficulty of growing both bacilli axenically, Mycobacterium leprae and M. lepraemurium were once thought to be closely related, although it was later suggested that M. lepraemurium might be related to Mycobacterium avium In this study, the complete genome of M. lepraemurium was sequenced using a combination of PacBio and Illumina sequencing. Phylogenomic analyses confirmed that M. lepraemurium is a distinct species within the M. avium complex (MAC). The M. lepraemurium genome is 4.05 Mb in length, which is considerably smaller than other MAC genomes, and it comprises 2,682 functional genes and 1,139 pseudogenes, which indicates that M. lepraemurium has undergone genome reduction. An error-prone repair homologue of the DNA polymerase III a-subunit was found to be nonfunctional in M. lepraemurium, which might contribute to pseudogene formation due to the accumulation of mutations in nonessential genes. M. lepraemurium has retained the functionality of several genes thought to influence virulence among members of the MAC.IMPORTANCEMycobacterium lepraemurium seems to be evolving toward a minimal set of genes required for an obligatory intracellular lifestyle within its host, a niche seldom adopted by most mycobacteria, as they are free-living. M. lepraemurium could be used as a model to elucidate functions of genes shared with other members of the MAC. Its reduced gene set can be exploited for studying the essentiality of genes in related pathogenic species, which might lead to discovery of common virulence factors or clarify host-pathogen interactions. M. lepraemurium can be cultivated in vitro only under specific conditions and even then with difficulty. Elucidating the metabolic (in)capabilities of M. lepraemurium will help develop suitable axenic media and facilitate genetic studies. Copyright © 2017 Benjak et al.


July 7, 2019

Complete genome sequence of Citrobacter werkmanii strain BF-6 isolated from industrial putrefaction.

In our previous study, Citrobacter werkmanii BF-6 was isolated from an industrial spoilage sample and demonstrated an excellent ability to form biofilms, which could be affected by various environmental factors. However, the genome sequence of this organism has not been reported so far.We report the complete genome sequence of C. werkmanii BF-6 together with the description of the genome features and its annotation. The size of the complete chromosome is 4,929,789 bp with an average coverage of 137×. The chromosome exhibits an average G + C content of 52.0%, and encodes 4570 protein coding genes, 84 tRNA genes, 25 rRNA operons, 3 microsatellite sequences and 34 minisatellite sequences. A previously unknown circular plasmid designated as pCW001 was also found with a length of 212,549 bp and a G + C content of 48.2%. 73.5%, 75.6% and 92.6% of the protein coding genes could be assigned to GO Ontology, KEGG Pathway, and COG (Clusters of Orthologous Groups) categories respectively. C. werkmanii BF-6 and C. werkmanii NRBC 105721 exhibited the closest evolutionary relationships based on 16S ribosomal RNA and core-pan genome assay. Furthermore, C. werkmanii BF-6 exhibits typical bacterial biofilm formation and development. In the RT-PCR experiments, we found that a great number of biofilm related genes, such as bsmA, bssR, bssS, hmsP, tabA, csgA, csgB, csgC, csgD, csgE, and csgG, were involved in C. werkmanii BF-6 biofilm formation.This is the first complete genome of C. werkmanii. Our work highlights the potential genetic mechanisms involved in biofilm formation and paves a way for further application of C. werkmanii in biofilms research.


July 7, 2019

Complete genome sequence of Paenibacillus yonginensis DCY84(T), a novel plant symbiont that promotes growth via induced systemic resistance.

This article reports the full genome sequence of Paenibacillus yonginensis DCY84(T) (KCTC33428, JCM19885), which is a Gram-positive rod-shaped bacterium isolated from humus soil of Yongin Forest in Gyeonggi Province, South Korea. The genome sequence of strain DCY84(T) provides greater understanding of the Paenibacillus species for practical use. This bacterium displays plant growth promotion via induced systemic resistance of abiotic stresses.


July 7, 2019

Trypanosoma cruzi specific mRNA amplification by in vitro transcription improves parasite transcriptomics in host-parasite RNA mixtures.

Trypanosomatids are a group of protozoan parasites that includes the etiologic agents of important human illnesses as Chagas disease, sleeping sickness and leishmaniasis. These parasites have a significant distinction from other eukaryotes concerning mRNA structure, since all mature mRNAs have an identical species-specific sequence of 39 nucleotides at the 5′ extremity, named spliced leader (SL). Considering this peculiar aspect of trypanosomatid mRNA, the aim of the present work was to develop a Trypanosoma cruzi specific in vitro transcription (IVT) linear mRNA amplification method in order to improve parasite transcriptomics analyses.We designed an oligonucleotide complementary to the last 21 bases of T. cruzi SL sequence, bearing an upstream T7 promoter (T7SL primer), which was used to direct the synthesis of second-strand cDNA. Original mRNA was then amplified by IVT using T7 RNA polymerase. T7SL-amplified RNA from two distinct T. cruzi stages (epimastigotes and trypomastigotes) were deep sequenced in SOLiD platform. Usual poly(A) + RNA and and T7-oligo(dT) amplified RNA (Eberwine method) were also sequenced. RNA-Seq reads were aligned to our new and improved T. cruzi Dm28c genome assembly (PacBio technology) and resulting transcriptome pattern from these three RNA preparation methods were compared, mainly concerning the conservation of mRNA transcritional levels and DEGs detection between epimastigotes and trypomastigotes.T7SL IVT method detected more potential differentially expressed genes in comparison to either poly(A) + RNA or T7dT IVT, and was also able to produce reliable quantifications of the parasite transcriptome down to 3 ng of total RNA. Furthermore, amplification of parasite mRNA in HeLa/epimastigote RNA mixtures showed that T7SL IVT generates transcriptome quantification with similar detection of differentially expressed genes when parasite RNA mass was only 0.1% of the total mixture (R = 0.78 when compared to poly(A) + RNA).The T7SL IVT amplification method presented here allows the detection of more potential parasite differentially expressed genes (in comparison to poly(A) + RNA) in host-parasite mixtures or samples with low amount of RNA. This method is especially useful for trypanosomatid transcriptomics because it produces less bias than PCR-based mRNA amplification. Additionally, by simply changing the complementary region of the T7SL primer, the present method can be applied to any trypanosomatid species.


July 7, 2019

Whole-genome assembly of Babesia ovata and comparative genomics between closely related pathogens.

Babesia ovata, belonging to the phylum Apicomplexa, is an infectious parasite of bovids. It is not associated with the manifestation of severe symptoms, in contrast to other types of bovine babesiosis caused by B. bovis and B. bigemina; however, upon co-infection with Theileria orientalis, it occasionally induces exacerbated symptoms. Asymptomatic chronic infection in bovines is usually observed only for B. ovata. Comparative genomic analysis could potentially reveal factors involved in these distinguishing characteristics; however, the genomic and molecular basis of these phenotypes remains elusive, especially in B. ovata. From a technical perspective, the current development of a very long read sequencer, MinION, will facilitate the obtainment of highly integrated genome sequences. Therefore, we applied next-generation sequencing to acquire a high-quality genome of the parasite, which provides fundamental information for understanding apicomplexans.The genome was assembled into 14,453,397 bp in size with 5031 protein-coding sequences (91 contigs and N50 = 2,090,503 bp). Gene family analysis revealed that ves1 alpha and beta, which belong to multigene families in B. bovis, were absent from B. ovata, the same as in B. bigemina. Instead, ves1a and ves1b, which were originally specified in B. bigemina, were present. The B. ovata and B. bigemina ves1a configure one cluster together even though they divided into two sub-clusters according to the spp. In contrast, the ves1b cluster was more dispersed and the overlap among B. ovata and B. bigemina was limited. The observed redundancy and rapid evolution in sequence might reflect the adaptive history of these parasites. Moreover, same candidate genes which potentially involved in the distinct phenotypes were specified by functional analysis. An anamorsin homolog is one of them. The human anamorsin is involved in hematopoiesis and the homolog was present in B. ovata but absent in B. bigemina which causes severe anemia.Taking these findings together, the differences demonstrated by comparative genomics potentially explain the evolutionary history of these parasites and the differences in their phenotypes. Besides, the draft genome provides fundamental information for further characterization and understanding of these parasites.


July 7, 2019

Public health surveillance in the UK revolutionises our understanding of the invasive Salmonella Typhimurium epidemic in Africa.

The ST313 sequence type of Salmonella Typhimurium causes invasive non-typhoidal salmonellosis and was thought to be confined to sub-Saharan Africa. Two distinct phylogenetic lineages of African ST313 have been identified.We analysed the whole genome sequences of S. Typhimurium isolates from UK patients that were generated following the introduction of routine whole-genome sequencing (WGS) of Salmonella enterica by Public Health England in 2014.We found that 2.7% (84/3147) of S. Typhimurium from patients in England and Wales were ST313 and were associated with gastrointestinal infection. Phylogenetic analysis revealed novel diversity of ST313 that distinguished UK-linked gastrointestinal isolates from African-associated extra-intestinal isolates. The majority of genome degradation of African ST313 lineage 2 was conserved in the UK-ST313, but the African lineages carried a characteristic prophage and antibiotic resistance gene repertoire. These findings suggest that a strong selection pressure exists for certain horizontally acquired genetic elements in the African setting. One UK-isolated lineage 2 strain that probably originated in Kenya carried a chromosomally located bla CTX-M-15, demonstrating the continual evolution of this sequence type in Africa in response to widespread antibiotic usage.The discovery of ST313 isolates responsible for gastroenteritis in the UK reveals new diversity in this important sequence type. This study highlights the power of routine WGS by public health agencies to make epidemiologically significant deductions that would be missed by conventional microbiological methods. We speculate that the niche specialisation of sub-Saharan African ST313 lineages is driven in part by the acquisition of accessory genome elements.


July 7, 2019

GRIDSS: sensitive and specific genomic rearrangement detection using positional de Bruijn graph assembly.

The identification of genomic rearrangements with high sensitivity and specificity using massively parallel sequencing remains a major challenge, particularly in precision medicine and cancer research. Here, we describe a new method for detecting rearrangements, GRIDSS (Genome Rearrangement IDentification Software Suite). GRIDSS is a multithreaded structural variant (SV) caller that performs efficient genome-wide break-end assembly prior to variant calling using a novel positional de Bruijn graph-based assembler. By combining assembly, split read, and read pair evidence using a probabilistic scoring, GRIDSS achieves high sensitivity and specificity on simulated, cell line, and patient tumor data, recently winning SV subchallenge #5 of the ICGC-TCGA DREAM8.5 Somatic Mutation Calling Challenge. On human cell line data, GRIDSS halves the false discovery rate compared to other recent methods while matching or exceeding their sensitivity. GRIDSS identifies nontemplate sequence insertions, microhomologies, and large imperfect homologies, estimates a quality score for each breakpoint, stratifies calls into high or low confidence, and supports multisample analysis.© 2017 Cameron et al.; Published by Cold Spring Harbor Laboratory Press.


July 7, 2019

The Babesia bovis hap2 gene is not required for blood stage replication, but expressed upon in vitro sexual stage induction.

Babesia bovis, is a tick borne apicomplexan parasite responsible for important cattle losses globally. Babesia parasites have a complex life cycle including asexual replication in the mammalian host and sexual reproduction in the tick vector. Novel control strategies aimed at limiting transmission of the parasite are needed, but transmission blocking vaccine candidates remain undefined. Expression of HAP2 has been recognized as critical for the fertilization of parasites in the Babesia-related Plasmodium, and is a leading candidate for a transmission blocking vaccine against malaria. Hereby we identified the B. bovis hap2 gene and demonstrated that it is widely conserved and differentially transcribed during development within the tick midgut, but not by blood stage parasites. The hap2 gene was disrupted by transfecting B. bovis with a plasmid containing the flanking regions of the hap2 gene and the GPF-BSD gene under the control of the ef-1a-B promoter. Comparison of in vitro growth between a hap2-KO B. bovis clonal line and its parental wild type strain showed that HAP2 is not required for the development of B. bovis in erythrocytes. However, xanthurenic acid-in vitro induction experiments of sexual stages of parasites recovered after tick transmission resulted in surface expression of HAP2 exclusively in sexual stage induced parasites. In addition, hap2-KO parasites were not able to develop such sexual stages as defined both by morphology and by expression of the B. bovis sexual marker genes 6-Cys A and B. Together, the data strongly suggests that tick midgut stage differential expression of hap2 is associated with the development of B. bovis sexual forms. Overall these studies are consistent with a role of HAP2 in tick stages of the parasite and suggest that HAP2 is a potential candidate for a transmission blocking vaccine against bovine babesiosis.


July 7, 2019

The sea cucumber genome provides insights into morphological evolution and visceral regeneration.

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.


Talk with an expert

If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.