N50 Length Research Articles

Comparative full-length transcriptome analysis by Oxford Nanopore Technologies reveals genes involved in anthocyanin accumulation in storage roots of sweet potatoes (Ipomoea batatasL.).

BackgroundStorage roots of sweet potatoes (Ipomoea batatas L.) with different colors vary in anthocyanin content, indicating different economically agronomic trait. As the newest DNA/RNA sequencing technology, Oxford Nanopore Technologies (ONT) have been applied in rapid transcriptome sequencing for investigation of genes related to nutrient metabolism. At present, few reports concern full-length transcriptome analysis based on ONT for study on the molecular mechanism of anthocyanin accumulation leading to color change of tuberous roots of sweet potato cultivars.ResultsThe storage roots of purple-fleshed sweet potato (PFSP) and white-fleshed sweet potato (WFSP) at different developmental stages were subjected to anthocyanin content comparison by UV-visible spectroscopy as well as transcriptome analysis at ONT MinION platform. UV-visible spectrophotometric measurements demonstrated the anthocyanin content of PFSP was much higher than that of WFSP. ONT RNA-Seq results showed each sample generated average 2.75 GB clean data with Full-Length Percentage (FL%) over 70% and the length of N50 ranged from 1,192 to 1,395 bp, indicating reliable data for transcriptome analysis. Subsequent analysis illustrated intron retention was the most prominent splicing event present in the resulting transcripts. As compared PFSP with WFSP at the relative developmental stages with the highest (PH vs. WH) and the lowest (PL vs. WL) anthocyanin content, 282 and 216 genes were up-regulated and two and 11 genes were down-regulated respectively. The differential expression genes involved in flavonoid biosynthesis pathway include CCoAOMT, PpLDOX, DFR, Cytochrome P450, CHI, and CHS. The genes encoding oxygenase superfamily were significantly up-regulated when compared PFSP with WFSP at the relative developmental stages.ConclusionsComparative full-length transcriptome analysis based on ONT serves as an effective approach to detect the differences in anthocyanin accumulation in the storage roots of different sweet potato cultivars at transcript level, with noting that some key genes can now be closely related to flavonoids biosynthesis. This study helps to improve understanding of molecular mechanism for anthocyanin accumulation in sweet potatoes and also provides a theoretical basis for high-quality sweet potato breeding.

Salinity Fluctuation on the Genetic Regulatory Mechanisms of the Crustacean, Charybdis japonica

The intertidal salinity fluctuations provide pose a significant physical challenge for the crustaceans that live here. In the present study, Charybdis japonica was selected as the research species to explore the genetic regulatory mechanisms of intertidal crustaceans responding to salinity fluctuations. We cultured C. japonica s with three salinity gradients (the salinity of 15‰, 25‰, and 35‰) and 68.90 Gb clean transcriptome reads were obtained by RNA-seq. All clean reads were then de novo assembled to 41,058 unigene with a mean length of 1,179.39 bp and an N50 length of 2,033 bp. Furthermore, a total of 18,100 (24.74%) unigenes were successfully matched with the sequences from at least one database. We further analyzed the transcriptome structure of C. japonica and a total of 26,853 CDSs and 12,190 SSRs were predicted. The gene expression levels of C. japonica at 25‰ were used as control, and 204 and 5,392 DEGs were found at the salinity of 35‰ and 15‰, respectively. Based on the annotated information of DEGs, we speculated that C. japonicas were subjected to greater stress under the salinity of 35‰ and isotonic intracellular regulation was activated. Furthermore, neuronal activity, oxidation response, intracellular regulatory activity and osmotic regulation-related genes were speculated to be critical genes for C. japonica to cope with salinity fluctuations. The present results will provide fundamental information for revealing the salinity fluctuation-related genetic regulatory mechanisms of C. japonica and other intertidal crustaceans, and then help us predict how future salinity changes will affect the survival of crustaceans.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus.

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.

Comparative transcriptome sequencing analysis of the narrow-clawed crayfish Pontastacus leptodactylus (Eschscholtz, 1823) and discovery of candidate sex-related genes

The comparative analysis of the gonads transcriptomes of the narrow-clawed crayfish Pontastacus leptodactylus were conducted to screen for sex-differentiation related genes based on Illumina RNA-seq technology. The transcriptome sequencing of cDNA libraries generated 543 million clean reads, and the N50 and mean length of the identified 59,348 unigenes were 1205 and 746 bp, respectively. A total of 4875 differentially expressed genes (DEGs) were obtained between the two sexes, of which 1627 predicted to be up-regulated in ovaries and 3248 in testes. A total of 86 candidate genes known to be associated with sexual differentiation and development were searched in the transcriptome of P. leptodactylus. Our results showed that the highly expressed male-related genes (dmrt1, amh, gsdf, spata 1, spata 13, spata 2, srm and smox) were significantly up -regulated to a substantial degree, suggesting the potential roles of these genes in gonads of P. leptodactylus. In addition, 11,989 potential simple sequence repeats (SSRs) and 190,655 putative single nucleotide polymorphisms (SNPs) were detected. The putative DEGs between the gonads provide an important source of information for further study of the sex-differentiation related genes in P. leptodactylus as well as development of molecular markers.

The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis.

BackgroundManchurian walnut (Juglans mandshurica Maxim.) is a tree with multiple industrial uses and medicinal properties in the Juglandaceae family (walnuts and hickories). J. mandshurica produces juglone, which is a toxic allelopathic agent and has potential utilization value. Furthermore, the seed of J. mandshurica is rich in various unsaturated fatty acids and has high nutritive value.FindingsHere, we present a high-quality chromosome-scale reference genome assembly and annotation for J. mandshurica (n = 16) with a contig N50 of 21.4 Mb by combining PacBio high-fidelity reads with high-throughput chromosome conformation capture data. The assembled genome has an estimated sequence size of 548.7 Mb and consists of 657 contigs, 623 scaffolds, and 40,453 protein-coding genes. In total, 60.99% of the assembled genome consists of repetitive sequences. Sixteen super-scaffolds corresponding to the 16 chromosomes were assembled, with a scaffold N50 length of 33.7 Mb and a BUSCO complete gene percentage of 98.3%. J. mandshurica displays a close sequence relationship with Juglans cathayensis, with a divergence time of 13.8 million years ago. Combining the high-quality genome, transcriptome, and metabolomics data, we constructed a gene-to-metabolite network and identified 566 core and conserved differentially expressed genes, which may be involved in juglone biosynthesis. Five CYP450 genes were found that may contribute to juglone accumulation. NAC, bZip, NF-YA, and NF-YC are positively correlated with the juglone content. Some candidate regulators (e.g., FUS3, ABI3, LEC2, and WRI1 transcription factors) involved in the regulation of lipid biosynthesis were also identified.ConclusionsOur genomic data provide new insights into the evolution of the walnut genome and create a new platform for accelerating molecular breeding and improving the comprehensive utilization of these economically important tree species.

High-quality genome assembly and pan-genome studies facilitate genetic discovery in mung bean and its improvement

Mung bean is an economically important legume crop species that is used as a food, consumed as a vegetable, and used as an ingredient and even as a medicine. To explore the genomic diversity of mung bean, we assembled a high-quality reference genome (Vrad_JL7) that was ∼479.35 Mb in size, with a contig N50 length of 10.34 Mb. A total of 40,125 protein-coding genes were annotated, representing ∼96.9% of the genetic region. We also sequenced 217 accessions, mainly landraces and cultivars from China, and identified 2,229,343 high-quality single-nucleotide polymorphisms (SNPs). Population structure revealed that the Chinese accessions diverged into two groups and were distinct from non-Chinese lines. Genetic diversity analysis based on genomic data from 750 accessions in 23 countries supported the hypothesis that mung bean was first domesticated in south Asia and introduced to east Asia probably through the Silk Road. We constructed the first pan-genome of mung bean germplasm and assembled 287.73 Mb of non-reference sequences. Among the genes, 83.1% were core genes and 16.9% were variable. Presence/absence variation (PAV) events of nine genes involved in the regulation of the photoperiodic flowering pathway were identified as being under selection during the adaptation process to promote early flowering in the spring. Genome-wide association studies (GWASs) revealed 2,912 SNPs and 259 gene PAV events associated with 33 agronomic traits, including a SNP in the coding region of the SWEET10 homolog (jg24043) involved in crude starch content and a PAV event in a large fragment containing 11 genes for color-related traits. This high-quality reference genome and pan-genome will provide insights into mung bean breeding.

Phasing analysis of lung cancer genomes using a long read sequencer

Chromosomal backgrounds of cancerous mutations still remain elusive. Here, we conduct the phasing analysis of non-small cell lung cancer specimens of 20 Japanese patients. By the combinatory use of short and long read sequencing data, we obtain long phased blocks of 834 kb in N50 length with >99% concordance rate. By analyzing the obtained phasing information, we reveal that several cancer genomes harbor regions in which mutations are unevenly distributed to either of two haplotypes. Large-scale chromosomal rearrangement events, which resemble chromothripsis events but have smaller scales, occur on only one chromosome, and these events account for the observed biased distributions. Interestingly, the events are characteristic of EGFR mutation-positive lung adenocarcinomas. Further integration of long read epigenomic and transcriptomic data reveal that haploid chromosomes are not always at equivalent transcriptomic/epigenomic conditions. Distinct chromosomal backgrounds are responsible for later cancerous aberrations in a haplotype-specific manner.

Global leaf and root transcriptome in response to cadmium reveals tolerance mechanisms in Arundo donax L

The expected increase of sustainable energy demand has shifted the attention towards bioenergy crops. Due to their know tolerance against abiotic stress and relatively low nutritional requirements, they have been proposed as election crops to be cultivated in marginal lands without disturbing the part of lands employed for agricultural purposes. Arundo donax L. is a promising bioenergy crop whose behaviour under water and salt stress has been recently studied at transcriptomic levels. As the anthropogenic activities produced in the last years a worrying increase of cadmium contamination worldwide, the aim of our work was to decipher the global transcriptomic response of A. donax leaf and root in the perspective of its cultivation in contaminated soil. In our study, RNA-seq libraries yielded a total of 416 million clean reads and 10.4 Gb per sample. De novo assembly of clean reads resulted in 378,521 transcripts and 126,668 unigenes with N50 length of 1812 bp and 1555 bp, respectively. Differential gene expression analysis revealed 5,303 deregulated transcripts (3,206 up- and 2,097 down regulated) specifically observed in the Cd-treated roots compared to Cd-treated leaves. Among them, we identified genes related to “Protein biosynthesis”, “Phytohormone action”, “Nutrient uptake”, “Cell wall organisation”, “Polyamine metabolism”, “Reactive oxygen species metabolism” and “Ion membrane transport”. Globally, our results indicate that ethylene biosynthesis and the downstream signal cascade are strongly induced by cadmium stress. In accordance to ethylene role in the interaction with the ROS generation and scavenging machinery, the transcription of several genes (NADPH oxidase 1, superoxide dismutase, ascorbate peroxidase, different glutathione S-transferases and catalase) devoted to cope the oxidative stress is strongly activated. Several small signal peptides belonging to ROTUNDIFOLIA, CLAVATA3, and C-TERMINALLY ENCODED PEPTIDE 1 (CEP) are also among the up-regulated genes in Cd-treated roots functioning as messenger molecules from root to shoot in order to communicate the stressful status to the upper part of the plants. Finally, the main finding of our work is that genes involved in cell wall remodelling and lignification are decisively up-regulated in giant reed roots. This probably represents a mechanism to avoid cadmium uptake which strongly supports the possibility to cultivate giant cane in contaminated soils in the perspective to reserve agricultural soil for food and feed crops.

A high-quality de novo genome assembly based on nanopore sequencing of a wild-caught coconut rhinoceros beetle (Oryctes rhinoceros)

BackgroundAn optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle’s mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation.ResultsHigh-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles.ConclusionsThe genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

Transcriptome Analysis of Pennisetum glaucum (L.) R. Br. Provides Insight Into Heat Stress Responses.

Pennisetum glaucum (L.) R. Br., being widely grown in dry and hot weather, frequently encounters heat stress at various stages of growth. The crop, due to its inherent capacity, efficiently overcomes such stress during vegetative stages. However, the same is not always the case with the terminal (flowering through grain filling) stages of growth, where recovery from stress is more challenging. However, certain pearl millet genotypes such as 841-B are known to overcome heat stress even at the terminal growth stages. Therefore, we performed RNA sequencing of two contrasting genotypes of pearl millet (841-B and PPMI-69) subjected to heat stress (42°C for 6 h) at flowering stages. Over 274 million high quality reads with an average length of 150 nt were generated, which were assembled into 47,310 unigenes having an average length of 1,254 nucleotides, N50 length of 1853 nucleotides, and GC content of 53.11%. Blastx resulted in the annotation of 35,628 unigenes, and functional classification showed 15,950 unigenes designated to 51 Gene Ontology terms. A total of 13,786 unigenes were allocated to 23 Clusters of Orthologous Groups, and 4,255 unigenes were distributed to 132 functional Kyoto Encyclopedia of Genes and Genomes database pathways. A total of 12,976 simple sequence repeats and 305,759 SNPs were identified in the transcriptome data. Out of 2,301 differentially expressed genes, 10 potential candidate genes were selected based on log2 fold change and adjusted p value parameters for their differential gene expression by qRT-PCR. We were able to identify differentially expressed genes unique to either of the two genotypes, and also, some DEGs common to both the genotypes were enriched. The differential expression patterns suggested that 841-B 6 h has better ability to maintain homeostasis during heat stress as compared to PPMI-69 6 h. The sequencing data generated in this study, like the SSRs and SNPs, shall serve as an important resource for the development of genetic markers, and the differentially expressed heat responsive genes shall be used for the development of transgenic crops.

Transcriptome analysis of Mythimna separata: De novo assembly and detection of genes related to beta-cypermethrin resistance.

The oriental armyworm Mythimna separata (Walker) is a devastating pest of cereal crops mainly in Asia and Oceania and recently become resistant to beta-cypermethrin (beta-CP). However, molecular biological studies of its response to beta-CP are scarce, and related genomic information is not available. In this study, we sequenced and de novo assembled the transcriptome of beta-CP susceptible M. separata (MsS-whole, abbr. MsS-W). A total of 30,486 unigenes were obtained, with an N50 length of 506 bp. A number of 21,051 unigenes were matched to public databases, of which 10,107 were classified into 59 gene ontology subcategories, 5792 were assigned into 25 clusters of orthologous groups of proteins subgroups and 12,123 were assigned to 257 Kyoto Encyclopedia of Genes and Genomes pathways. A total of 729 simple sequence repeatswere detected. In addition, a total of 323 cytochrome P450-associated sequences from nine lepidopterous species, of which 130 were from M. separata, were analyzed using the maximum likelihood method and Bayesian inference. Among the 130 cytochrome P450-associated sequences from M. separata, 60 were dropped into CYP3 clan, which is associated with metabolizing xenobiotics and plant natural compounds. Furthermore, the beta-CP susceptible (MsS-2) and resistant (MsR-2) M. separata population transcriptomes were sequenced. Certain critical genes involved in beta-CP detoxification were detected and verified by quantitative real-time polymerase chain reaction. Collectively, our results provided a basis for further studies of the molecular mechanism of insecticide resistance in M. separata.

The first draft genome assembly and data analysis of the Malaysian mahseer (Tor tambroides)

The Malaysian mahseer (Tor tambroides), one of the most valuable freshwater fish in the world, is mainly targeted for human consumption. The mitogenomic data of this species is available to date, but the genomic information is still lacking. For the first time, we sequenced the whole genome of an adult fish on both Illumina and Nanopore platforms. The hybrid genome assembly had resulted in a sum of 1.23 Gb genomic sequence from the 44,726 contigs found with 44 kb N50 length and BUSCO genome completeness of 87.6%. Four types of SSRs had been detected and identified within the genome with a greater AT abundance than that of GC. Predicted protein sequences had been functionally annotated to public databases, namely GO, KEGG and COG. A maximum likelihood phylogenomic tree containing 52 Actinopterygii species and one Sarcopterygii species as outgroup was constructed, providing first insights into the genome-based evolutionary relationship of T. tambroides with other ray-finned fish. These data are crucial in facilitating the study of population genomics, species identification, morphological variations, and evolutionary biology, which are helpful in the conservation of this species.

Chromosome-level genome assembly of Indian mangrove (Ceriops tagal) revealed a genome-wide duplication event predating the divergence of Rhizophoraceae mangrove species.

Mangrove ecosystems are unique, highly diverse, provide benefits to humans, and aid in coastal protection. The Indian mangrove, or spurred mangrove, [Ceriops tagal (Perr.) C. B. Rob.] is a member of the Rhizophoraceae family and is commonly found along the intertidal zones in tropical regions in Southeast Asia, southern Asia, and Africa. Here, we present the first high-quality reference genome assembly of the Ceriops species. A preliminary draft assembly, generated from the 10× Genomics linked-read library, was scaffolded using the proximity ligation chromatin contact mapping technique (Hi-C) to obtain a chromosome-scale assembly of 231,919,005 bases with an N50 length of 11,408,429 bases. The benchmarking universal single-copy orthologs (BUSCO) analysis revealed that C. tagal gene predictions recovered 95.8% of the highly conserved orthologs. Phylogenetic analyses suggested that C. tagal diverged from the last common ancestor of flat-leaf spurred mangrove [C. decandra (Griff.) Ding Hou] and C. zippeliana Blume ∼10.4 million yr ago (MYA), and the last common ancestor of genera Ceriops, Kandelia, and Rhizophora diverged from that of genus Bruguiera ∼49.4 MYA. In addition, our analysis of the transversion rate at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a recent whole-genome duplication in C. tagal. The STRUCTURE and principal component analyses illustrated that C. tagal individuals investigated in this study were the admixture of two subpopulations, the genetic background of which was influenced primarily by location. The availability of genomic and transcriptomic resources and biodiversity data reported in this work will be useful for future studies that may shed light on adaptive evolutions of mangrove species.

A High-Quality Genome Assembly of Striped Catfish (Pangasianodon hypophthalmus) Based on Highly Accurate Long-Read HiFi Sequencing Data.

The HiFi sequencing technology yields highly accurate long-read data with accuracies greater than 99.9% that can be used to improve results for complex applications such as genome assembly. Our study presents a high-quality chromosome-scale genome assembly of striped catfish (Pangasianodon hypophthalmus), a commercially important species cultured mainly in Vietnam, integrating HiFi reads and Hi-C data. A 788.4 Mb genome containing 381 scaffolds with an N50 length of 21.8 Mb has been obtained from HiFi reads. These scaffolds have been further ordered and clustered into 30 chromosome groups, ranging from 1.4 to 57.6 Mb, based on Hi-C data. The present updated assembly has a contig N50 of 14.7 Mb, representing a 245-fold and 4.2-fold improvement over the previous Illumina and Illumina-Nanopore-Hi-C based version, respectively. In addition, the proportion of repeat elements and BUSCO genes identified in our genome is remarkably higher than in the two previously released striped catfish genomes. These results highlight the power of using HiFi reads to assemble the highly repetitive regions and to improve the quality of genome assembly. The updated, high-quality genome assembled in this work will provide a valuable genomic resource for future population genetics, conservation biology and selective breeding studies of striped catfish.

A chromosome-level genome assembly of the pollinating fig wasp Valisia javana.

Fig wasp has always been thought the species-specific pollinator for their host fig (Moraceae, Ficus) and constitute a model system with its host to study co-evolution and co-speciation. The availability of a high-quality genome will help to further reveal the mechanisms underlying these characteristics. Here, we present a high-quality chromosome-level genome for Valisa javana developed by a combination of PacBio long-read and Illumina short-read. The assembled genome size is 296.34 Mb from 13 contigs with a contig N50 length of 26.76 kb. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid fig wasps living in syconium of its highly specific host. Protein-coding genes associated with chemosensory, detoxification and venom genes were identified. Several differentially expressed genes in transcriptome data of V. javana between odor-stimulated samples and the controls have been identified in some olfactory signal transduction pathways, e.g. olfactory transduction, cAMP, cGMP-PKG, Calcim, Ras and Rap1. This study provides a valuable genomic resource for a fig wasp, and sheds insight into further revealing the mechanisms underlying their adaptive traits to their hosts in different places and co-speciation with their host.

Chromosome-level assembly, annotation and phylome of Pelobates cultripes, the western spadefoot toad.

Genomic resources for amphibians are still hugely under-represented in vertebrate genomic research, despite being a group of major interest for ecology, evolution and conservation. Amphibians constitute a highly threatened group of vertebrates, present a vast diversity in reproductive modes, are extremely diverse in morphology, occupy most ecoregions of the world, and present the widest range in genome sizes of any major group of vertebrates. We combined Illumina, Nanopore and Hi-C sequencing technologies to assemble a chromosome-level genome sequence for an anuran with a moderate genome size (assembly span 3.09 Gb); Pelobates cultripes, the western spadefoot toad. The genome has an N50 length of 330 Mb with 98.6% of the total sequence length assembled into 14 super scaffolds, and 87.7% complete BUSCO genes. We use published transcriptomic data to provide annotations, identifying 32,684 protein-coding genes. We also reconstruct the P. cultripes phylome and identify 2,527 gene expansions. We contribute the first draft of the genome of the western spadefoot toad, P. cultripes. This species represents a relatively basal lineage in the anuran tree with an interesting ecology and a high degree of developmental plasticity, and thus is an important resource for amphibian genomic research.

TrBase: A genome and transcriptome database of Temnopleurus reevesii.

Sea urchins have a long history as model organisms in biology, but their use in genetics is limited because of their long breeding cycle. In sea urchin genetics, genome editing technology was first established in Hemicentrotus pulcherrimus, whose genome has already been published. However, because this species also has a long breeding cycle, new model sea urchins that are more suitable for genetics have been sought. Here, we report a draft genome of another Western Pacific species, Temnopleurus reevesii, which we established as a new model sea urchin recently since this species has a comparable developmental process to other model sea urchins but a short breeding cycle of approximately half a year. The genome of T. reevesii was assembled into 28,742 scaffold sequences with an N50 length of 67.6kb and an estimated genome size of 905.9Mb. In the assembled genome, 27,064 genes were identified, 23,624 of which were expressed in at least one of the seven developmental stages. To provide genetic information, we constructed the genome database TrBase(https://cell-innovation.nig.ac.jp/Tree/). We also constructed the Western Pacific Sea Urchin Genome Database (WestPac-SUGDB) (https://cell-innovation.nig.ac.jp/WPAC/) with the aim of establishing a portal site for genetic information on sea urchins in the West Pacific. This site contains genomic information on two species, T. reevesii and H. pulcherrimus, and is equipped with homology search programs for comparing the two datasets. Therefore, TrBase and WestPac-SUGDB are expected to contribute not only to genetic research using sea urchins but also to comparative genomics and evolutionary research.

Full-Length Transcriptome Analysis Provides New Insights Into the Diversity of Immune-Related Genes in Portunus trituberculatus.

Generally, invertebrates were thought to solely rely on their non-specific innate immune system to fight against invading microorganisms. However, increasing studies have implied that the innate immune response of invertebrates displayed diversity and specificity owing to the hyper-variable immune molecules in organisms. In order to get an insight into the diversity of immune-related genes in Portunus trituberculatus, a full-length transcriptome analysis of several immune-related tissues (hemocytes, hepatopancreas and gills) in P. trituberculatus was performed and the diversity of several immune-related genes was analyzed. The full-length transcriptome analysis of P. trituberculatus was conducted using a combination of SMRT long-read sequencing and Illumina short-read sequencing. A total of 17,433 nonredundant full-length transcripts with average length of 2,271 bp and N50 length of 2,841 bp were obtained, among which 13,978 (80.18%) transcripts were annotated. Moreover, numerous transcript variants of various immune-related genes were identified, including pattern recognition receptors, antimicrobial peptides, heat shock proteins (HSPs), antioxidant enzymes and vital molecules in prophenoloxidase (proPO)-activating system. Based on the full-length transcriptome analysis, open reading frames (ORFs) of four C-type lectins (CTLs) were cloned, and tissue distributions showed that the four CTLs were ubiquitously expressed in all the tested tissues, and mainly expressed in hepatopancreas and gills. The transcription of the four CTLs significantly increased in several immune-related tissues (hemocytes, hepatopancreas and gills) of P. trituberculatus challenged with Vibrio alginolyticus and displayed different profiles. Moreover, the four CTLs displayed distinct bacterial binding and antibacterial activities. The recombinant protein PtCTL-1 (rPtCTL-1) and rPtCTL-3 displayed bacterial binding and antibacterial activities against all tested bacteria. rPtCTL-2 only showed bacterial binding and antibacterial activities against V. alginolyticus. No obvious bacterial binding or antibacterial activities for PtCTL-4 was observed against the tested bacteria. This study enriches the transcriptomic information on P. trituberculatus and provides new insights into the innate immune system of crustaceans. Additionally, our study provided candidates of antibiotic agents for the prevention and treatment of bacteriosis.

Whole-Genome Sequence Resource of Fusarium oxysporum Strain TH15, a Plant Growth Promoting Endophytic Fungus Isolated from Tetrastigma hemsleyanum

Whole-Genome Sequence Resource of <i>Fusarium oxysporum</i> Strain TH15, a Plant Growth Promoting Endophytic Fungus Isolated from <i>Tetrastigma hemsleyanum</i>

Chromosomal-Level Genome Assembly of the Springtail Tomocerus qinae (Collembola: Tomoceridae).

The family Tomoceridae is among the earliest derived collembolan lineages, thus is of key importance in understanding the evolution of Collembola. Here, we assembled a chromosome-level genome of one tomocerid species Tomocerus qinae by combining Nanopore long reads and Hi-C data. The final genome size was 334.44 Mb with the scaffold/contig N50 length of 71.85/13.94 Mb. BUSCO assessment indicated that 96.80% of complete arthropod universal single-copy orthologs (n = 1,013) were present in the assembly. The repeat elements accounted for 26.11% (87.26 Mb) and 494 noncoding RNAs were identified in the genome. A total of 20,451 protein-coding genes were predicted, which captured 96.0% (973) BUSCO genes. Gene family evolution analyses identified 4,825 expanded gene families of T. qinae, among them, 47 experienced significant expansions, and these significantly expanded gene families mainly involved in proliferation and growth. This study provides an important genomic resource for future evolution and comparative genomics analyses of Collembola.