Pacific Biosciences Single-molecule Real-time Sequencing Research Articles

Genome of the Southern Giant Petrel Assembled Using Third-Generation DNA Sequencing and Linked Reads Reveals Evolutionary Traits of Southern Avian.

Simple SummaryThe southern giant petrel Macronectes giganteus is one of key species on the Antarctic and sub-Antarctic regions. We reported the whole genome of M. giganteus for understanding the evolutionary mechanisms on Antarctic environments and for studying a more effective genetic monitoring of threatened species. The M. giganteus genome was 1.248 Gb in size with a scaffold N50 length of 27.4 Mb and a longest scaffold length of 120.4 Mb, and its genome contains 14,993 predicted protein-coding genes. The estimated historical effective population size of southern giant petrel underwent a severe reduction during a period coinciding with the early Pleistocene. The M. giganteus showed genomic expansion related to maintenance of energy homeostasis, being essential for survival and effective functioning in cold environments. Moreover, we employed a classification of microsatellite markers for studying the genetic diversity within and among populations. Genomic research of this first Antarctic bird helps to address the environmental adaptation and evolution of avian species.The southern giant petrel Macronectes giganteus, a large seabird of the southern oceans, is one of only two members of the genus Macronectes and is the largest species in the order Procellariiformes. Although these two families account for the vast majority of the avian fauna inhabiting the Antarctic and sub-Antarctic regions, studies on the status of some populations and the associated genetic data are currently extremely limited. In this study, we assembled the genome of M. giganteus by integrating Pacific Biosciences single-molecule real-time sequencing and the Chromium system developed by 10x Genomics. The final M. giganteus genome assembly was 1.248 Gb in size with a scaffold N50 length of 27.4 Mb and a longest scaffold length of 120.4 Mb. The M. giganteus genome contains 14,993 predicted protein-coding genes and has 11.06% repeat sequences. Estimated historical effective population size analysis indicated that the southern giant petrel underwent a severe reduction in effective population size during a period coinciding with the early Pleistocene. The availability of this newly sequenced genome will facilitate more effective genetic monitoring of threatened species. Furthermore, the genome will provide a valuable resource for gene functional studies and further comparative genomic studies on the life history and ecological traits of specific avian species.

The genome of Magnolia biondii Pamp. provides insights into the evolution of Magnoliales and biosynthesis of terpenoids

Magnolia biondii Pamp. (Magnoliaceae, magnoliids) is a phylogenetically, economically, and medicinally important ornamental tree species widely grown and cultivated in the north-temperate regions of China. Determining the genome sequence of M. biondii would help resolve the phylogenetic uncertainty of magnoliids and improve the understanding of individual trait evolution within the Magnolia genus. We assembled a chromosome-level reference genome of M. biondii using ~67, ~175, and ~154 Gb of raw DNA sequences generated via Pacific Biosciences single-molecule real-time sequencing, 10X Genomics Chromium, and Hi-C scaffolding strategies, respectively. The final genome assembly was ~2.22 Gb, with a contig N50 value of 269.11 kb and a BUSCO complete gene percentage of 91.90%. Approximately 89.17% of the genome was organized into 19 chromosomes, resulting in a scaffold N50 of 92.86 Mb. The genome contained 47,547 protein-coding genes, accounting for 23.47% of the genome length, whereas 66.48% of the genome length consisted of repetitive elements. We confirmed a WGD event that occurred very close to the time of the split between the Magnoliales and Laurales. Functional enrichment of the Magnolia-specific and expanded gene families highlighted genes involved in the biosynthesis of secondary metabolites, plant–pathogen interactions, and responses to stimuli, which may improve the ecological fitness and biological adaptability of the lineage. Phylogenomic analyses revealed a sister relationship of magnoliids and Chloranthaceae, which are sister to a clade comprising monocots and eudicots. The genome sequence of M. biondii could lead to trait improvement, germplasm conservation, and evolutionary studies on the rapid radiation of early angiosperms.

Pacific Biosciences assembly with Hi-C mapping generates an improved, chromosome-level goose genome.

ABSTRACTBackgroundThe domestic goose is an economically important and scientifically valuable waterfowl; however, a lack of high-quality genomic data has hindered research concerning its genome, genetics, and breeding. As domestic geese breeds derive from both the swan goose (Anser cygnoides) and the graylag goose (Anser anser), we selected a female Tianfu goose for genome sequencing. We generated a chromosome-level goose genome assembly by adopting a hybrid de novo assembly approach that combined Pacific Biosciences single-molecule real-time sequencing, high-throughput chromatin conformation capture mapping, and Illumina short-read sequencing.FindingsWe generated a 1.11-Gb goose genome with contig and scaffold N50 values of 1.85 and 33.12 Mb, respectively. The assembly contains 39 pseudo-chromosomes (2n = 78) accounting for ∼88.36% of the goose genome. Compared with previous goose assemblies, our assembly has more continuity, completeness, and accuracy; the annotation of core eukaryotic genes and universal single-copy orthologs has also been improved. We have identified 17,568 protein-coding genes and a repeat content of 8.67% (96.57 Mb) in this genome assembly. We also explored the spatial organization of chromatin and gene expression in the goose liver tissues, in terms of inter-pseudo-chromosomal interaction patterns, compartments, topologically associating domains, and promoter-enhancer interactions.ConclusionsWe present the first chromosome-level assembly of the goose genome. This will be a valuable resource for future genetic and genomic studies on geese.

The chromosome-level draft genome of Dalbergia odorifera.

Background Dalbergia odorifera T. Chen (Fabaceae) is an International Union for Conservation of Nature red-listed tree. This tree is of high medicinal and commercial value owing to its officinal, insect-proof, durable heartwood. However, there is a lack of genome reference, which has hindered development of studies on the heartwood formation.FindingsWe presented the first chromosome-scale genome assembly of D. odorifera obtained on the basis of Illumina paired-end sequencing, Pacific Biosciences single-molecule real-time sequencing, 10x Genomics linked reads, and Hi-C technology. We assembled 97.68% of the 653.45 Mb D. odorifera genome with scaffold N50 and contig sizes of 56.16 and 5.92 Mb, respectively. Ten super-scaffolds corresponding to the 10 chromosomes were assembled, with the longest scaffold reaching 79.61 Mb. Repetitive elements account for 54.17% of the genome, and 30,310 protein-coding genes were predicted from the genome, of which ∼92.6% were functionally annotated. The phylogenetic tree showed that D. odorifera diverged from the ancestor of Arabidopsis thaliana and Populus trichocarpa and then separated from Glycine max and Cajanus cajan.ConclusionsWe sequence and reveal the first chromosome-level de novo genome of D. odorifera. These studies provide valuable genomic resources for the research of heartwood formation in D. odorifera and other timber trees. The high-quality assembled genome can also be used as reference for comparative genomics analysis and future population genetic studies of D. odorifera.

Analysis of the Antifungal Properties of Bacillus velezensis B-4 Through a Bioassay and Complete-Genome Sequencing.

The strain B-4, isolated from a field in Changsha (China), presents strong antifungal activities, as identified by the Kirby–Bauer test, especially for pathogens that harm crops. Here, we obtained the complete genome sequence of the strain B-4 by Pacific Biosciences single-molecule real-time sequencing, making it well analyzed for understanding mechanisms and creating biological agents. Its 3,919-kb circular chromosome genome has 3,725 protein-coding genes [coding sequences (CDSs)] and 46.7% guanine–cytosine content. A comparative genome analysis of B-4 with other published strains (including Bacillus velezensis, Bacillus amyloliquefaciens, and Bacillus subtilis) revealed that the strain B-4 is a B. velezensis strain. These different strains have 2,889 CDSs in common, whereas 179 CDSs were found to be unique in the strain B-4, which is a far greater number than that in other strains. Regarding the antifungal activities of B-4, we were specifically concerned with the genes involved in the biosynthesis of secondary metabolites. In total, more than 19.56% of the genome was annotated to 12 gene clusters relating to synthesis of antimicrobial metabolites, which contained various enzyme-encoding operons for non-ribosomal peptide synthetases, polyketide synthases, and lantipeptide synthesis proteins. They were all considered to be related to the production of bacteriostatic substances or stimulation of induced systemic resistance by bacterial metabolites. These situations also present an advantage over those of other strains for biocontrol potential. We provide evidence that the biological control effect of the strain B-4, as demonstrated in antibacterial activity experiments and predicted from the complete genome sequence analysis, provides the basis for research promoting agricultural research on sustainable development, especially the contribution of biotechnology to agriculture.

Chromosome-Level Clam Genome Helps Elucidate the Molecular Basis of Adaptation to a Buried Lifestyle.

SummaryBivalve mollusks are economically important invertebrates that exhibit marked diversity in benthic lifestyle and provide valuable resources for understanding the molecular basis of adaptation to benthic life. In this report, we present a high-quality, chromosome-anchored reference genome of the Venus clam, Cyclina sinensis. The chromosome-level genome was assembled by Pacific Bioscience single-molecule real-time sequencing, Illumina paired-end sequencing, 10× Genomics, and high-throughput chromosome conformation capture technologies. The final genome assembly of C. sinensis is 903.2 Mb in size, with a contig N50 size of 2.6 Mb and a scaffold N50 size of 46.5 Mb. Enrichment analyses of significantly expanded and positively selected genes suggested evolutionary adaptation of this clam to buried life. In addition, a change in shell color represents another mechanism of adaptation to burial in sediment. The high-quality genome generated in this work provides a valuable resource for investigating the molecular mechanisms of adaptation to buried lifestyle.

Single-molecule long-read sequencing of the full-length transcriptome of Rhododendron lapponicum L.

Rhododendron lapponicum L. is a familiar ornamental plant worldwide with important ornamental and economic value. However, a full-length R. lapponicum transcriptome is still lacking. In the present study, we used the Pacific Biosciences single-molecule real-time sequencing technology to generate the R. lapponicum transcriptome. A total of 346,270 full-length non-chimeric reads were generated, from which we obtained 75,002 high-quality full-length transcripts. We identified 55,255 complete open reading frames, 7,140 alternative splicing events and 2,011 long non-coding RNAs. In gene annotation analyses, 71,155, 33,653, 30,359 and 31,749 transcripts were assigned to the Nr, GO, COG and KEGG databases, respectively. Additionally, 3,150 transcription factors were detected. KEGG pathway analysis showed that 96 transcripts were identified coding for the enzymes associated with anthocyanin synthesis. Furthermore, we identified 64,327 simple sequence repeats from 45,319 sequences, and 150 pairs of primers were randomly selected to develop SSR markers. This study provides a large number of full-length transcripts, which will facilitate the further study of the genetics of R. lapponicum.

Profiling of the viable bacterial and fungal microbiota in fermented feeds using single-molecule real-time sequencing.

Fermented concentrated feed has been widely recognized as an ideal feed in the animal industry. In this study, we used a powerful method, coupling propidium monoazide (PMA) pretreatment with single-molecule real-time (SMRT) sequencing technology to compare the bacterial and fungal composition of feeds before and after fermentation with four added lactic acid bacteria (LAB) inoculants (one Lactobacillus casei strain and three L. plantarum strains). Five feed samples consisting of corn, soybean meal, and wheat bran were fermented with LAB additives for 3 d. Following anaerobic fermentation, the pH rapidly decreased, and the mean numbers of LAB increased from 106 to 109 colony-forming units (cfu)/g fresh matter. SMRT sequencing results showed that the abundance and diversity of bacteria and fungi in the feed were significantly higher before fermentation than after fermentation. Fifteen bacterial species and eight fungal genera were significantly altered following fermentation, and L. plantarum was the dominant species (relative abundance 88.94%) in the post-fermentation group. PMA treatment revealed that the bacteria Bacillus cereus, B. circulans, Alkaliphilus oremlandii, Cronobacter sakazakii, Paenibacillus barcinonensis, and P. amylolyticus (relative abundance >1%) were viable in the raw feed. After fermentation, their relative abundances decreased sharply to <0.2%; however, viable L. plantarum was still the dominant species post fermentation. We inferred that our LAB additives grew rapidly and inhibited harmful microorganisms and further improved feed quality. In addition, coupling PMA treatment with the Pacific Biosciences SMRT sequencing technology was a powerful tool for providing accurate live microbiota profiling data in this study.

Chromosome-level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818).

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high-quality chromosome-level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single-molecule real-time sequencing, Illumina paired-end sequencing, 10X Genomics linked-reads and Hi-C reads. The genome size was 1,220.85Mb, containing scaffold N50 of 65.93Mb and contig N50 of 976.94Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single-Copy Orthologs were determined in this genome assembly. And Hi-C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein-coding genes were predicted in the S.constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S.constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S.constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.

A high-quality genome assembly for the endangered golden snub-nosed monkey (Rhinopithecus roxellana).

BackgroundThe golden snub-nosed monkey (Rhinopithecus roxellana) is an endangered colobine species endemic to China, which has several distinct traits including a unique social structure. Although a genome assembly for R. roxellana is available, it is incomplete and fragmented because it was constructed using short-read sequencing technology. Thus, important information such as genome structural variation and repeat sequences may be absent.FindingsTo obtain a high-quality chromosomal assembly for R. roxellana qinlingensis, we used 5 methods: Pacific Bioscience single-molecule real-time sequencing, Illumina paired-end sequencing, BioNano optical maps, 10X Genomics link-reads, and high-throughput chromosome conformation capture. The assembled genome was ∼3.04 Gb, with a contig N50 of 5.72 Mb and a scaffold N50 of 144.56 Mb. This represented a 100-fold improvement over the previously published genome. In the new genome, 22,497 protein-coding genes were predicted, of which 22,053 were functionally annotated. Gene family analysis showed that 993 and 2,745 gene families were expanded and contracted, respectively. The reconstructed phylogeny recovered a close relationship between R. rollexana and Macaca mulatta, and these 2 species diverged ∼13.4 million years ago.ConclusionWe constructed a high-quality genome assembly of the Qinling golden snub-nosed monkey; it had superior continuity and accuracy, which might be useful for future genetic studies in this species and as a new standard reference genome for colobine primates. In addition, the updated genome assembly might improve our understanding of this species and could assist conservation efforts.

The genome assembly and annotation of yellowhorn (Xanthoceras sorbifolium Bunge).

BackgroundYellowhorn (Xanthoceras sorbifolium Bunge), a deciduous shrub or small tree native to north China, is of great economic value. Seeds of yellowhorn are rich in oil containing unsaturated long-chain fatty acids that have been used for producing edible oil and nervonic acid capsules. However, the lack of a high-quality genome sequence hampers the understanding of its evolution and gene functions.FindingsIn this study, a whole genome of yellowhorn was sequenced and assembled by integration of Illumina sequencing, Pacific Biosciences single-molecule real-time sequencing, 10X Genomics linked reads, Bionano optical maps, and Hi-C. The yellowhorn genome assembly was 439.97 Mb, which comprised 15 pseudo-chromosomes covering 95.42% (419.84 Mb) of the assembled genome. The repetitive fractions accounted for 56.39% of the yellowhorn genome. The genome contained 21,059 protein-coding genes. Of them, 18,503 (87.86%) genes were found to be functionally annotated with ≥1 "annotation" term by searching against other databases. Transcriptomic analysis showed that 341, 135, 125, 113, and 100 genes were specifically expressed in hermaphrodite flower, staminate flower, young fruit, leaf, and shoot, respectively. Phylogenetic analysis suggested that yellowhorn and Dimocarpus longan diverged from their most recent common ancestor ∼46 million years ago.ConclusionsThe availability and subsequent annotation of the yellowhorn genome, as well as the identification of tissue-specific functional genes, provides a valuable reference for plant comparative genomics, evolutionary studies, and molecular design breeding.

Draft genome of Glyptosternon maculatum, an endemic fish from Tibet Plateau.

sBackgroundMechanisms for high-altitude adaption have attracted widespread interest among evolutionary biologists. Several genome-wide studies have been carried out for endemic vertebrates in Tibet, including mammals, birds, and amphibians. However, little information is available about the adaptive evolution of highland fishes. Glyptosternon maculatum (Regan 1905), also known as Regan or barkley and endemic to the Tibetan Plateau, belongs to the Sisoridae family, order Siluriformes (catfishes). This species lives at an elevation ranging from roughly 2,800 m to 4,200 m. Hence, a high-quality reference genome of G. maculatum provides an opportunity to investigate high-altitude adaption mechanisms of fishes.FindingsTo obtain a high-quality reference genome sequence of G. maculatum, we combined Pacific Bioscience single-molecule real-time sequencing, Illumina paired-end sequencing, 10X Genomics linked-reads, and BioNano optical map techniques. In total, 603.99 Gb sequencing data were generated. The assembled genome was about 662.34 Mb with scaffold and contig N50 sizes of 20.90 Mb and 993.67 kb, respectively, which captured 83% complete and 3.9% partial vertebrate Benchmarking Universal Single-Copy Orthologs. Repetitive elements account for 35.88% of the genome, and 22,066 protein-coding genes were predicted from the genome, of which 91.7% have been functionally annotated.ConclusionsWe present the first comprehensive de novo genome of G. maculatum. This genetic resource is fundamental for investigating the origin of G. maculatum and will improve our understanding of high-altitude adaption of fishes. The assembled genome can also be used as reference for future population genetic studies of G. maculatum.

A single-molecule sequencing assay for the comprehensive profiling of T4 DNA ligase fidelity and bias during DNA end-joining.

DNA ligases are key enzymes in molecular and synthetic biology that catalyze the joining of breaks in duplex DNA and the end-joining of DNA fragments. Ligation fidelity (discrimination against the ligation of substrates containing mismatched base pairs) and bias (preferential ligation of particular sequences over others) have been well-studied in the context of nick ligation. However, almost no data exist for fidelity and bias in end-joining ligation contexts. In this study, we applied Pacific Biosciences Single-Molecule Real-Time sequencing technology to directly sequence the products of a highly multiplexed ligation reaction. This method has been used to profile the ligation of all three-base 5′-overhangs by T4 DNA ligase under typical ligation conditions in a single experiment. We report the relative frequency of all ligation products with or without mismatches, the position-dependent frequency of each mismatch, and the surprising observation that 5′-TNA overhangs ligate extremely inefficiently compared to all other Watson–Crick pairings. The method can easily be extended to profile other ligases, end-types (e.g. blunt ends and overhangs of different lengths), and the effect of adjacent sequence on the ligation results. Further, the method has the potential to provide new insights into the thermodynamics of annealing and the kinetics of end-joining reactions.

Genome sequencing and analysis of the first complete genome of Lactobacillus kunkeei strain MP2, an Apis mellifera gut isolate.

Background. The honey bee (Apis mellifera) is the most important pollinator in agriculture worldwide. However, the number of honey bees has fallen significantly since 2006, becoming a huge ecological problem nowadays. The principal cause is CCD, or Colony Collapse Disorder, characterized by the seemingly spontaneous abandonment of hives by their workers. One of the characteristics of CCD in honey bees is the alteration of the bacterial communities in their gastrointestinal tract, mainly due to the decrease of Firmicutes populations, such as the Lactobacilli. At this time, the causes of these alterations remain unknown. We recently isolated a strain of Lactobacillus kunkeei (L. kunkeei strain MP2) from the gut of Chilean honey bees. L. kunkeei, is one of the most commonly isolated bacterium from the honey bee gut and is highly versatile in different ecological niches. In this study, we aimed to elucidate in detail, the L. kunkeei genetic background and perform a comparative genome analysis with other Lactobacillus species.Methods. L. kunkeei MP2 was originally isolated from the guts of Chilean A. mellifera individuals. Genome sequencing was done using Pacific Biosciences single-molecule real-time sequencing technology. De novo assembly was performed using Celera assembler. The genome was annotated using Prokka, and functional information was added using the EggNOG 3.1 database. In addition, genomic islands were predicted using IslandViewer, and pro-phage sequences using PHAST. Comparisons between L. kunkeei MP2 with other L. kunkeei, and Lactobacillus strains were done using Roary.Results. The complete genome of L. kunkeei MP2 comprises one circular chromosome of 1,614,522 nt. with a GC content of 36,9%. Pangenome analysis with 16 L. kunkeei strains, identified 113 unique genes, most of them related to phage insertions. A large and unique region of L. kunkeei MP2 genome contains several genes that encode for phage structural protein and replication components. Comparative analysis of MP2 with other Lactobacillus species, identified several unique genes of L. kunkeei MP2 related with metabolism, biofilm generation, survival under stress conditions, and mobile genetic elements (MGEs).Discussion. The presence of multiple mobile genetic elements, including phage sequences, suggest a high degree of genetic variability in L. kunkeei. Its versatility and ability to survive in different ecological niches (bee guts, flowers, fruits among others) could be given by its genetic capacity to change and adapt to different environments. L. kunkeei could be a new source of Lactobacillus with beneficial properties. Indeed, L. kunkeei MP2 could play an important role in honey bee nutrition through the synthesis of components as isoprenoids.

Denoising DNA deep sequencing data-high-throughput sequencing errors and their correction.

Characterizing the errors generated by common high-throughput sequencing platforms and telling true genetic variation from technical artefacts are two interdependent steps, essential to many analyses such as single nucleotide variant calling, haplotype inference, sequence assembly and evolutionary studies. Both random and systematic errors can show a specific occurrence profile for each of the six prominent sequencing platforms surveyed here: 454 pyrosequencing, Complete Genomics DNA nanoball sequencing, Illumina sequencing by synthesis, Ion Torrent semiconductor sequencing, Pacific Biosciences single-molecule real-time sequencing and Oxford Nanopore sequencing. There is a large variety of programs available for error removal in sequencing read data, which differ in the error models and statistical techniques they use, the features of the data they analyse, the parameters they determine from them and the data structures and algorithms they use. We highlight the assumptions they make and for which data types these hold, providing guidance which tools to consider for benchmarking with regard to the data properties. While no benchmarking results are included here, such specific benchmarks would greatly inform tool choices and future software development. The development of stand-alone error correctors, as well as single nucleotide variant and haplotype callers, could also benefit from using more of the knowledge about error profiles and from (re)combining ideas from the existing approaches presented here.