PacBio Single-molecule Sequencing Research Articles

Chromosome level genome assembly of giant freshwater prawn (Macrobrachium rosenbergii)

The giant freshwater prawn (Macrobrachium rosenbergii) has many advantages in aquaculture, such as fast growth rate, short breeding cycle and good nutrition, which makes it a freshwater shrimp with high economic value. Herein, high-quality chromosome-level genome of both female and male prawns were obtained by combining Illumina paired-end sequencing, PacBio single molecule sequencing technique and High-through chromosome conformation capture (Hi-C) technologies. In ZZ male prawn, a final contig assembly of 3118.58 Mb with a N50 length of 956,237 bp was obtained. In WW female prawn, a final contig assembly of 3333.31 Mb with a N50 length of 1,143,555 bp was obtained. The assembled genome sequences from prawns were anchored to 59 chromosomes. Moreover, the sex chromosomes including W chromosome and Z chromosome were generated in prawn with the length of 36.23 Mb and 27.33 Mb, respectively. The sequence similarity of Z chromosome and W chromosome reached to 74.90%. The high-quality genome resource will be useful for further molecular breeding and functional genomic research of giant freshwater prawns.

Endophytic bacteria with allelopathic potential regulate gene expression and metabolite production in host Casuarina equisetifolia.

Casuarina equisetifolia is a common protective forest in coastal areas. However, artificial C. equisetifolia forests cannot self-renew, mainly due to the accumulation of allelochemicals. Endophytic bacteria may alleviate the root growth inhibition caused by allelochemicals in C. equisetifolia seedlings. B. amyloliquefaciens and B. aryabhattai were endophytic bacteria with strong allelopathy in C. equisetifolia root. The allelopathy mechanism of these two endophytes and their interaction with C. equisetifolia remains to be studied. Whole-genome sequencing of B. amyloliquefaciens and B. aryabhattai isolated from the roots of allelochemical-accumulating C. equisetifolia was performed using Illumina Hiseq and PacBio single-molecule sequencing platforms. Sterile seedlings of C. equisetifolia were treated with either individual or mixed bacterial cultures through root drenching. Transcriptional and metabolomics analyses were conducted after 3 days of infection. Whole-genome sequencing of Bacillus aryabhattai and Bacillus amyloliquefaciens showed that the two strains contained various horizontal gene transfer elements such as insertion sequence, prophage and transposon. In addition, these two strains also contain numerous genes related to the synthesis and catabolism of allelochemicals. After these two strains of bacteria were individually or mixed infected with C. equisetifolia, metabolomics and transcriptomic analysis of C. equisetifolia showed the 11 important secondary metabolite biosynthesis among them alkaloids biosynthesis, phenylpropanoid and terpenes biosynthesis and related genes were putatively regulated. Correlation analysis revealed that 48 differentially expressed genes had strong positive correlations with 42 differential metabolites, and 48 differentially expressed genes had strong negative correlations with 36 differential metabolites. For example, CMBL gene showed positive correlations with the allelochemical (-)-Catechin gallate, while Bp10 gene showed negative correlations with (-)-Catechin gallate. The intergenerational accumulation of allelochemicals may induce horizontal gene transfer in endogenic bacteria of Casuarina equisetifolia root. Endophytic Bacillus plays an allelopathic role by assisting the host in regulating gene expression and the production and/or variety of allelochemicals. This comprehensive study sheds light on the intricate genetic and metabolic interactions between Bacillus endophytes and C. equisetifolia. These findings provide insights into endophyte-mediated allelopathy and its potential uses in plant biology and forest sustainability.

Single-Molecule Sequencing of the C9orf72 Repeat Expansion in Patient iPSCs.

A hexanucleotide GGGGCC repeat expansion in the C9orf72 gene is the most frequent genetic cause of amyotrophic lateral sclerosis (ALS) and frontal temporal dementia (FTD). C9orf72 repeat expansions are currently identified with long-range PCR or Southern blot for clinical and research purposes, but these methods lack accuracy and sensitivity. The GC-rich and repetitive content of the region cannot be amplified by PCR, which leads traditional sequencing approaches to fail. We turned instead to PacBio single-molecule sequencing to detect and size the C9orf72 repeat expansion without amplification. We isolated high molecular weight genomic DNA from patient-derived iPSCs of varying repeat lengths and then excised the region containing the C9orf72 repeat expansion from naked DNA with a CRISPR/Cas9 system. We added adapters to the cut ends, capturing the target region for sequencing on PacBio's Sequel, Sequel II, or Sequel IIe. This approach enriches the C9orf72 repeat region without amplification and allows the repeat expansion to be consistently and accurately sized, even for repeats in the thousands. Key features • This protocol is adapted from PacBio's previous "no-amp targeted sequencing utilizing the CRISPR-Cas9 system." • Optimized for sizing C9orf72 repeat expansions in patient-derived iPSCs and applicable to DNA from any cell type, blood, or tissue. • Requires high molecular weight naked DNA. • Compatible with Sequel I and II but not Revio.

The chromosome level genome assembly of the Asian green mussel, Perna viridis

The Asian green mussel, Perna viridis is an important aquaculture species in the family Mytilidae contributing substantially to molluscan aquaculture. We generated a high-quality chromosome level assembly of this species by combining PacBio single molecule sequencing technique (SMRT), Illumina paired-end sequencing, high-throughput chromosome conformation capture technique (Hi-C) and Bionano mapping. The final assembly resulted in a genome of 723.49 Mb in size with a scaffold N50 of 49.74 Mb with 99% anchored into 15 chromosomes. A total of 49654 protein-coding genes were predicted from the genome. The presence of 634 genes associated with the cancer pathway and 408 genes associated with viral carcinogenesis indicates the potential of this species to be used as a model for cancer studies. The chromosome-level assembly of this species is also a valuable resource for further genomic selection and selective breeding for improving economically important aquaculture traits and augmenting aquaculture productivity.

Chromosome-level genome assembly and annotation of the yellow grouper, Epinephelus awoara

Epinephelus awoara, as known as yellow grouper, is a significant economic marine fish that has been bred artificially in China. However, the genetic structure and evolutionary history of yellow grouper remains largely unknown. Here, this work presents the high-quality chromosome-level genome assembly of yellow grouper using PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The 984.48 Mb chromosome-level genome of yellow grouper was assembled, with a contig N50 length of 39.77 Mb and scaffold N50 length of 41.39 Mb. Approximately 99.76% of assembled sequences were anchored into 24 pseudo-chromosomes with the assistance of Hi-C reads. Furthermore, approximately 41.17% of the genome was composed of repetitive elements. In total, 24,541 protein-coding genes were predicted, of which 22,509 (91.72%) genes were functionally annotated. The highly accurate, chromosome-level reference genome assembly and annotation are crucial to the understanding of population genetic structure, adaptive evolution and speciation of the yellow grouper.

Structure of plant-associated microeukaryotes in roots and leaves of aquatic and terrestrial plants revealed by blocking peptide-nucleic acid (PNA) amplification.

Studies of plant-microbe interactions, including mutualistic, antagonistic, parasitic, or commensal microbes, have greatly benefited our understanding of ecosystem functioning. New molecular identification tools have increasingly revealed the association patterns between microorganisms and plants. Here, we integrated long-read PacBio single-molecule sequencing technology with a blocking protein-nucleic acid (PNA) approach to minimise plant amplicons in a survey of plant-eukaryotic microbe relationships in roots and leaves of different aquatic and terrestrial plants to determine patterns of organ, host, and habitat preferences. The PNA approach reduced the samples' relative amounts of plant reads and did not distort the fungal and other microeukaryotic composition. Our analyses revealed that the eukaryotic microbiomes associated with leaves and roots of aquatic plants exhibit a much larger proportion of non-fungal microorganisms than terrestrial plants, and leaf and root microbiomes are similar. Terrestrial plants had much stronger differentiation of leaf and root microbiomes and stronger partner specificity than aquatic plants.

A Global Survey of the Full-Length Transcriptome of Apis mellifera by Single-Molecule Long-Read Sequencing.

As important pollinators, honey bees play a crucial role in both maintaining the ecological balance and providing products for humans. Although several versions of the western honey bee genome have already been published, its transcriptome information still needs to be refined. In this study, PacBio single-molecule sequencing technology was used to sequence the full-length transcriptome of mixed samples from many developmental time points and tissues of A. mellifera queens, workers and drones. A total of 116,535 transcripts corresponding to 30,045 genes were obtained. Of these, 92,477 transcripts were annotated. Compared to the annotated genes and transcripts on the reference genome, 18,915 gene loci and 96,176 transcripts were newly identified. From these transcripts, 136,554 alternative splicing (AS) events, 23,376 alternative polyadenylation (APA) sites and 21,813 lncRNAs were detected. In addition, based on the full-length transcripts, we identified many differentially expressed transcripts (DETs) between queen, worker and drone. Our results provide a complete set of reference transcripts for A. mellifera that dramatically expand our understanding of the complexity and diversity of the honey bee transcriptome.

Chromosome-level genome assembly of largemouth bass (Micropterus salmoides) using PacBio and Hi-C technologies

The largemouth bass (Micropterus salmoides) has become a cosmopolitan species due to its widespread introduction as game or domesticated fish. Here a high-quality chromosome-level reference genome of M. salmoides was produced by combining Illumina paired-end sequencing, PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. Ultimately, the genome was assembled into 844.88 Mb with a contig N50 of 15.68 Mb and scaffold N50 length of 35.77 Mb. About 99.9% assembly genome sequences (844.00 Mb) could be anchored to 23 chromosomes, and 98.03% assembly genome sequences could be ordered and directed. The genome contained 38.19% repeat sequences and 2693 noncoding RNAs. A total of 26,370 protein-coding genes from 3415 gene families were predicted, of which 97.69% were functionally annotated. The high-quality genome assembly will be a fundamental resource to study and understand how M. salmoides adapt to novel and changing environments around the world, and also be expected to contribute to the genetic breeding and other research.

Chromosome-Level Assembly of the Southern Rock Bream (Oplegnathus fasciatus) Genome Using PacBio and Hi-C Technologies.

The Rock Bream (Oplegnathus fasciatus) is an economically important rocky reef fish of the Northwest Pacific Ocean. In recent years, it has been cultivated as an important edible fish in coastal areas of China. Despite its economic importance, genome-wide adaptions of domesticated O. fasciatus are largely unknown. Here we report a chromosome-level reference genome of female O. fasciatus (from the southern population in the subtropical region) using the PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The genome was assembled into 120 contigs with a total length of 732.95 Mb and a contig N50 length of 27.33 Mb. After chromosome-level scaffolding, 24 chromosomes with a total length of 723.22 Mb were constructed. Moreover, a total of 27,015 protein-coding genes and 5,880 ncRNAs were annotated in the reference genome. This reference genome of O. fasciatus will provide an important resource not only for basic ecological and population genetic studies but also for dissect artificial selection mechanisms in marine aquaculture.

A chromosome-level genome assembly of the yellowfin seabream (Acanthopagrus latus; Hottuyn, 1782) provides insights into its osmoregulation and sex reversal

The yellowfin seabream Acanthopagrus latus is the economically most important Sparidae fish in the northern South China Sea. As euryhaline fish, they are perfect model for investigating osmoregulatory mechanisms in teleosts. Moreover, the reproductive biology of hermaphrodites has long been intriguing; however, little information is known about the molecular pathways underlying their sex change. Here, we report a chromosome level reference genome of A. latus generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The draft genome of yellowfin seabream was 806 Mb, with 732 Mb scaffolds anchored on 24 chromosomes. The contig N50 and scaffold N50 were 2.6 Mb and 30.17 Mb, respectively. The assembly is of high integrity and includes 92.23% universal single-copy orthologues based on benchmarking universal single-copy orthologs (BUSCO) analysis. A total of 19,631 protein-coding genes were functionally annotated in the reference genome. Moreover, ARRDC3 and GSTA gene families which related to osmoregulation underwent an extensive expansion in two euryhaline sparids fish genomes compared to other teleost genomes. Moreover, integrating sex-specific transcriptome analyses, several genes related to the transforming growth factor beta (TGF-β) signalling pathway involved in sex differentiation and development. This genomic resource will not only be valuable for studying the osmoregulatory mechanisms in estuarine fish and sex determination in hermaphrodite vertebrate species, but also provide useful genomic tools for facilitating breeding of the yellowfin seabream.

Repeat-Driven Generation of Antigenic Diversity in a Major Human Pathogen,Trypanosoma cruzi.

Trypanosoma cruzi, a zoonotic kinetoplastid protozoan parasite, is the causative agent of American trypanosomiasis (Chagas disease). Having a very plastic, repetitive and complex genome, the parasite displays a highly diverse repertoire of surface molecules, with pivotal roles in cell invasion, immune evasion and pathogenesis. Before 2016, the complexity of the genomic regions containing these genes impaired the assembly of a genome at chromosomal level, making it impossible to study the structure and function of the several thousand repetitive genes encoding the surface molecules of the parasite. We here describe the genome assembly of the Sylvio X10/1 genome sequence, which since 2016 has been used as a reference genome sequence for T. cruzi clade I (TcI), produced using high coverage PacBio single-molecule sequencing. It was used to analyze deep Illumina sequence data from 34 T. cruzi TcI isolates and clones from different geographic locations, sample sources and clinical outcomes. Resolution of the surface molecule gene distribution showed the unusual duality in the organization of the parasite genome, a synteny of the core genomic region with related protozoa flanked by unique and highly plastic multigene family clusters encoding surface antigens. The presence of abundant interspersed retrotransposons in these multigene family clusters suggests that these elements are involved in a recombination mechanism for the generation of antigenic variation and evasion of the host immune response on these TcI strains. The comparative genomic analysis of the cohort of TcI strains revealed multiple cases of such recombination events involving surface molecule genes and has provided new insights into T. cruzi population structure.

A high-quality de novo genome assembly of one swamp eel (Monopterus albus) strain with PacBio and Hi-C sequencing data.

The swamp eel (Monopterus albus) is one economically important fish in China and South-Eastern Asia and a good model species to study sex inversion. There are different genetic lineages and multiple local strains of swamp eel in China, and one local strain of M. albus with deep yellow and big spots has been selected for consecutive selective breeding due to superiority in growth rate and fecundity. A high-quality reference genome of the swamp eel would be a very useful resource for future selective breeding program. In the present study, we applied PacBio single-molecule sequencing technique (SMRT) and the high-throughput chromosome conformation capture (Hi-C) technologies to assemble the M. albus genome. A 799 Mb genome was obtained with the contig N50 length of 2.4 Mb and scaffold N50 length of 67.24 Mb, indicating 110-fold and ∼31.87-fold improvement compared to the earlier released assembly (∼22.24 Kb and 2.11 Mb, respectively). Aided with Hi‐C data, a total of 750 contigs were reliably assembled into 12 chromosomes. Using 22,373 protein-coding genes annotated here, the phylogenetic relationships of the swamp eel with other teleosts showed that swamp eel separated from the common ancestor of Zig-zag eel ∼49.9 million years ago, and 769 gene families were found expanded, which are mainly enriched in the immune system, sensory system, and transport and catabolism. This highly accurate, chromosome-level reference genome of M. albus obtained in this work will be used for the development of genome-scale selective breeding.

Community Structure of Arbuscular Mycorrhizal Fungi in Soils of Switchgrass Harvested for Bioenergy.

Learning more about the biodiversity and composition of arbuscular mycorrhizal fungi (AMF) under alternative agricultural management scenarios may be important to the sustainable intensification of switchgrass grown as a bioenergy crop. Using PacBio single-molecule sequencing and taxonomic resolution to the level of amplicon sequence variant (ASV), we assessed the effects of nitrogen amendment on AMF associating with switchgrass and explored relationships between AMF and switchgrass yield across three sites of various productivities in Wisconsin. Nitrogen amendment had little effect on AMF diversity metrics or community composition. While AMF ASV diversity was not correlated with switchgrass yield, AMF family richness and switchgrass yield had a strong, positive relationship at one of our three sites. Each of our sites was dominated by unique ASVs of the species Paraglomus brasilianum, indicating regional segregation of AMF at the intraspecific level. Our molecular biodiversity survey identified putative core members of the switchgrass microbiome, as well as novel clades of AMF, especially in the order Paraglomerales and the genus Nanoglomus Furthermore, our phylogenies unite the cosmopolitan, soil-inhabiting clade deemed GS24 with Pervetustaceae, an enigmatic family prevalent in stressful environments. Future studies should isolate and characterize the novel genetic diversity found in switchgrass agroecosystems and explore the potential yield benefits of AMF richness.IMPORTANCE We assessed the different species of beneficial fungi living in agricultural fields of switchgrass, a large grass grown for biofuels, using high-resolution DNA sequencing. Contrary to our expectations, the fungi were not greatly affected by fertilization. However, we found a positive relationship between plant productivity and the number of families of beneficial fungi at one site. Furthermore, we sequenced many species that could not be identified with existing reference databases. One group of fungi was highlighted in an earlier study for being widely distributed but of unknown taxonomy. We discovered that this group belonged to a family called Pervetustaceae, which may benefit switchgrass in stressful environments. To produce higher-yielding switchgrass in a more sustainable manner, it could help to study these undescribed fungi and the ways in which they may contribute to greater switchgrass yield in the absence of fertilization.

Fully Phased Sequence of a Diploid Human Genome Determined de Novo from the DNA of a Single Individual.

In recent years, improved sequencing technology and computational tools have made de novo genome assembly more accessible. Many approaches, however, generate either an unphased or only partially resolved representation of a diploid genome, in which polymorphisms are detected but not assigned to one or the other of the homologous chromosomes. Yet chromosomal phase information is invaluable for the understanding of phenotypic trait inheritance in the cases of compound heterozygosity, allele-specific expression or cis-acting variants. Here we use a combination of tools and sequencing technologies to generate a de novo diploid assembly of the human primary cell line WI-38. First, data from PacBio single molecule sequencing and Bionano Genomics optical mapping were combined to generate an unphased assembly. Next, 10x Genomics linked reads were combined with the hybrid assembly to generate a partially phased assembly. Lastly, we developed and optimized methods to use short-read (Illumina) sequencing of flow cytometry-sorted metaphase chromosomes to provide phase information. The final genome assembly was almost fully (94%) phased with the addition of approximately 2.5-fold coverage of Illumina data from the sequenced metaphase chromosomes. The diploid nature of the final de novo genome assembly improved the resolution of structural variants between the WI-38 genome and the human reference genome. The phased WI-38 sequence data are available for browsing and download at wi38.research.calicolabs.com. Our work shows that assembling a completely phased diploid genome de novo from the DNA of a single individual is now readily achievable.

A high-quality reference genome of wild Cannabis sativa

Cannabis sativa is a well-known plant species that has great economic and ecological significance. An incomplete genome of cloned C. sativa was obtained by using SOAPdenovo software in 2011. To further explore the utilization of this plant resource, we generated an updated draft genome sequence for wild-type varieties of C. sativa in China using PacBio single-molecule sequencing and Hi-C technology. Our assembled genome is approximately 808 Mb, with scaffold and contig N50 sizes of 83.00 Mb and 513.57 kb, respectively. Repetitive elements account for 74.75% of the genome. A total of 38,828 protein-coding genes were annotated, 98.20% of which were functionally annotated. We provide the first comprehensive de novo genome of wild-type varieties of C. sativa distributed in Tibet, China. Due to long-term growth in the wild environment, these varieties exhibit higher heterozygosity and contain more genetic information. This genetic resource is of great value for future investigations of cannabinoid metabolic pathways and will aid in promoting the commercial production of C. sativa and the effective utilization of cannabinoids. The assembled genome is also a valuable resource for intensively and effectively investigating the C. sativa genome further in the future.

Single molecule sequencing reveals response of manganese-oxidizing microbiome to different biofilter media in drinking water systems

Rapid sand biofiltration (RSBF) is widely used for the removal of contaminants from drinking water treatment systems. Biofilm microbiomes in the biofilter media play essential roles in biotransformation of contaminants, but is not comprehensively understood. This study reports on Mn(II) oxidation and the core microbiomes in magnetite sand RSBF (MagS–RSBF) and manganese sand RSBF (MnS–RSBF). MnS–RSBF showed a relatively higher Mn(II) removal rate (40–91.2%) than MagS–RSBF during the start-up. MagS–RSBF and MnS–RSBF had similar Mn(II) removal rates (94.13% and 99.16%) over stable operation for 80 days. Mn(II) removal rates at different depths in the MnS–RSBF reactor significantly changed with operation time, and the filter in the upper layer of MnS–RSBF made the largest contribution to Mn(II) oxidation once operation had stabilized. PacBio single molecule sequencing of full-length 16S rRNA gene indicated that biofilter medium had a significant impact on the core microbiomes of the biofilms from the two biofilters. The magnetite sand biofilter facilitated the enrichment of Mn(II)-oxidizing biofilms. The dominant populations consisted of Pedomicrobium, Pseudomonas, and Hyphomicrobium in the RSBF, which have been affiliated with putative manganese-oxidizing bacteria (MnOB). The relative abundance of Pedomicrobium manganicum increased with operation time in both RSBF reactors. In addition, Nordella oligomobilis and Derxia gummosa were statistically correlated with Mn(II) oxidation. Species-species co-occurrence networks indicated that the microbiome of MnS–RSBF had more complex correlations than that of MagS–RSBF, implying that biofilter medium substantially shaped the microbial community in the RSBF. Hyphomicrobium and nitrite-oxidizing Nitrospira moscoviensis were positively correlated. The core microbiomes’ composition of both RSBF reactors converged over operation time. A hybrid biofilter medium with magnetite sand and manganese sand may therefore be best in rapid sand filtration for Mn(II) oxidation.

The sequencing and de novo assembly of the Larimichthys crocea genome using PacBio and Hi-C technologies

Larimichthys crocea is an endemic marine fish in East Asia that belongs to Sciaenidae in Perciformes. L. crocea has now been recognized as an “iconic” marine fish species in China because not only is it a popular food fish in China, it is a representative victim of overfishing and still provides high value fish products supported by the modern large-scale mariculture industry. Here, we report a chromosome-level reference genome of L. crocea generated by employing the PacBio single molecule sequencing technique (SMRT) and high-throughput chromosome conformation capture (Hi-C) technologies. The genome sequences were assembled into 1,591 contigs with a total length of 723.86 Mb and a contig N50 length of 2.83 Mb. After chromosome-level scaffolding, 24 scaffolds were constructed with a total length of 668.67 Mb (92.48% of the total length). Genome annotation identified 23,657 protein-coding genes and 7262 ncRNAs. This highly accurate, chromosome-level reference genome of L. crocea provides an essential genome resource to support the development of genome-scale selective breeding and restocking strategies of L. crocea.

The sequence and de novo assembly of Takifugu bimaculatus genome using PacBio and Hi-C technologies

Takifugu bimaculatus is a native teleost species of the southeast coast of China where it has been cultivated as an important edible fish in the last decade. Genetic breeding programs, which have been recently initiated for improving the aquaculture performance of T. bimaculatus, urgently require a high-quality reference genome to facilitate genome selection and related genetic studies. To address this need, we produced a chromosome-level reference genome of T. bimaculatus using the PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. The genome was assembled into 2,193 contigs with a total length of 404.21 Mb and a contig N50 length of 1.31 Mb. After chromosome-level scaffolding, 22 chromosomes with a total length of 371.68 Mb were constructed. Moreover, a total of 21,117 protein-coding genes and 3,471 ncRNAs were annotated in the reference genome. The highly accurate, chromosome-level reference genome of T. bimaculatus provides an essential genome resource for not only the genome-scale selective breeding of T. bimaculatus but also the exploration of the evolutionary basis of the speciation and local adaptation of the Takifugu genus.

Genome Assembly and Annotation of the Trichoplusia ni Tni-FNL Insect Cell Line Enabled by Long-Read Technologies.

Background: Trichoplusia ni derived cell lines are commonly used to enable recombinant protein expression via baculovirus infection to generate materials approved for clinical use and in clinical trials. In order to develop systems biology and genome engineering tools to improve protein expression in this host, we performed de novo genome assembly of the Trichoplusia ni-derived cell line Tni-FNL. Methods: By integration of PacBio single-molecule sequencing, Bionano optical mapping, and 10X Genomics linked-reads data, we have produced a draft genome assembly of Tni-FNL. Results: Our assembly contains 280 scaffolds, with a N50 scaffold size of 2.3 Mb and a total length of 359 Mb. Annotation of the Tni-FNL genome resulted in 14,101 predicted genes and 93.2% of the predicted proteome contained recognizable protein domains. Ortholog searches within the superorder Holometabola provided further evidence of high accuracy and completeness of the Tni-FNL genome assembly. Conclusions: This first draft Tni-FNL genome assembly was enabled by complementary long-read technologies and represents a high-quality, well-annotated genome that provides novel insight into the complexity of this insect cell line and can serve as a reference for future large-scale genome engineering work in this and other similar recombinant protein production hosts.

Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.