High-quality Assembly Research Articles

The high-quality genome of Grona styracifolia uncovers the genomic mechanism of high levels of schaftoside, a promising drug candidate for treatment of COVID-19.

Recent study has evidenced that traditional Chinese medicinal (TCM) plant-derived schaftoside shows promise as a potential drug candidate for COVID-19 treatment. However, the biosynthetic pathway of schaftoside in TCM plants remains unknown. In this study, the genome of the TCM herb Grona styracifolia (Osbeck) H.Ohashi & K.Ohashi (GSO), which is rich in schaftoside, was sequenced, and a high-quality assembly of GSO genome was obtained. Our findings revealed that GSO did not undergo recent whole genome duplication (WGD) but shared an ancestral papilionoid polyploidy event, leading to the gene expansion of chalcone synthase (CHS) and isoflavone 2'-hydroxylase (HIDH). Furthermore, GSO-specific tandem gene duplication resulted in the gene expansion of C-glucosyltransferase (CGT). Integrative analysis of the metabolome and transcriptome identified 13 CGTs and eight HIDHs involved in the biosynthetic pathway of schaftoside. Functional studies indicated that CGTs and HIDHs identified here are bona fide responsible for the biosynthesis of schaftoside in GSO, as confirmed through hairy root transgenic system and in vitro enzyme activity assay. Taken together, the ancestral papilionoid polyploidy event expanding CHSs and HIDHs, along with the GSO-specific tandem duplication of CGT, contributes, partially if not completely, to the robust biosynthesis of schaftoside in GSO. These findings provide insights into the genomic mechanisms underlying the abundant biosynthesis of schaftoside in GSO, highlighting the potential of GSO as a source of bioactive compounds for pharmaceutical development.

Discerning the dissemination mechanisms of antibiotic resistance genes through whole genome sequencing of extended-spectrum beta-lactamase (ESBL)-producing E. coli isolated from veterinary clinics and farms in South Korea

Extended-spectrum beta-lactamase (ESBL)-producing bacteria are resistant to most beta-lactams, including third-generation cephalosporins, limiting the treatment methods against the infections they cause. In this study, we performed whole genome sequencing of ESBL-producing E. coli to determine the mechanisms underlying the dissemination of antibiotic resistance genes. We analyzed 141 ESBL-producing isolates which had been collected from 16 veterinary clinics and 16 farms in South Korea. Long- and short-read sequencing platforms were used to obtain high-quality assemblies. The results showed that blaCTX-M is the dominant ESBL gene type found in South Korea. The spread of blaCTX-M appears to have been facilitated by both clonal spread between different host species and conjugation. Most blaCTX-M genes were found associated with diverse mobile genetic elements that may contribute to the chromosomal integration of the genes. Diverse incompatibility groups of blaCTX-M-harboring plasmids were also observed, which allows their spread among a variety of bacteria. Comprehensive whole genome sequence analysis was useful for the identification of the most prevalent types of ESBL genes and their dissemination mechanisms. The results of this study suggest that the propagation of ESBL genes can occur through clonal spread and plasmid-mediated dissemination, and that suitable action plans should be developed to prevent further propagation of these genes.

Technology-enabled great leap in deciphering plant genomes.

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

High quality assemblies of four indigenous chicken genomes and related functional data resources

Many lines of evidence indicate that red jungle fowl (RJF) is the primary ancestor of domestic chickens. Although multiple versions of RJF (galgal2-galgal5 and GRCg6a) and commercial chickens (GRCg7b/w and Huxu) genomes have been assembled since 2004, no high-quality indigenous chicken genomes have been assembled, hampering the understanding of chicken domestication and evolution. To fill the gap, we sequenced the genomes of four indigenous chickens with distinct morphological traits in southwest China, using a combination of short, long and Hi-C reads. We assembled each genome (~1.0 Gb) into 42 chromosomes with chromosome N50 90.5–90.9 Mb, amongst the highest quality of chicken genome assemblies. To provide resources for gene annotation and functional analysis, we also sequenced transcriptomes of 10 tissues for each of the four chickens. Moreover, we corrected many mis-assemblies and assembled missing micro-chromosomes 29 and 34–39 for GRCg6a. Our assemblies, sequencing data and the correction of GRCg6a can be valuable resources for studying chicken domestication and evolution.

Constructing telomere-to-telomere diploid genome by polishing haploid nanopore-based assembly.

Draft genomes generated from Oxford Nanopore Technologies (ONT) long reads are known to have a higher error rate. Although existing genome polishers can enhance their quality, the error rate (including mismatches, indels and switching errors between paternal and maternal haplotypes) can be significant. Here, we develop two polishers, hypo-short and hypo-hybrid to address this issue. Hypo-short utilizes Illumina short reads to polish an ONT-based draft assembly, resulting in a high-quality assembly with low error rates and switching errors. Expanding on this, hypo-hybrid incorporates ONT long reads to further refine the assembly into a diploid representation. Leveraging on hypo-hybrid, we have created a diploid genome assembly pipeline called hypo-assembler. Hypo-assembler automates the generation of highly accurate, contiguous and nearly complete diploid assemblies using ONT long reads, Illumina short reads and optionally Hi-C reads. Notably, our solution even allows for the production of telomere-to-telomere diploid genomes with additional manual steps. As a proof of concept, we successfully assembled a fully phased telomere-to-telomere diploid genome of HG00733, achieving a quality value exceeding 50.

Chromosome-level genome assembly of Platycarya strobilacea

Platycarya strobilacea belongs to the walnut family (Juglandaceae), is commonly known as species endemic to East Asia, and is an ecologically important, wind pollinated, woody deciduous tree. To facilitate this ancient tree for the ecological value and conservation of this ancient tree, we report a new high-quality genome assembly of P. strobilacea. The genome size was 677.30 Mb, with a scaffold N50 size of 45,791,698 bp, and 98.43% of the assembly was anchored to 15 chromosomes. We annotated 32,246 protein-coding genes in the genome, of which 96.30% were functionally annotated in six databases. This new high-quality assembly of P. strobilacea provide valuable resource for the phylogenetic and evolutionary analysis of the walnut family and angiosperm.

A de novo chromosome-scale assembly of the Lablab purpureus genome.

Lablab (Lablab purpureus (L.) Sweet), an underutilized tropical legume crop, plays a crucial role in global food and nutritional security. To enhance our understanding of its genetic makeup towards developing elite cultivars, we sequenced and assembled a draft genome of L. purpureus accession PK2022T020 using a single tube long fragment read (stLFR) technique. The preliminary assembly encompassed 367 Mb with a scaffold N50 of 4.3 Mb. To improve the contiguity of our draft genome, we employed a chromatin contact mapping (Hi-C) approach to obtain a pseudochromosome-level assembly containing 366 Mb with an N50 length of 31.1 Mb. A total of 327.4 Mb had successfully been anchored into 11 pseudomolecules, corresponding to the haploid chromosome number in lablab. Our gene prediction recovered 98.4% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Comparative analyses utilizing sequence information from single-copy orthologous genes demonstrated that L. purpureus diverged from the last common ancestor of the Phaseolus/Vigna species approximately 27.7 million years ago. A gene family expansion analysis revealed a significant expansion of genes involved in responses to biotic and abiotic stresses. Our high-quality chromosome-scale reference assembly provides an invaluable genomic resource for lablab genetic improvement and future comparative genomics studies among legume species.

Genome Assembly and Annotation of the Dark-Branded Bushbrown Butterfly Mycalesis mineus (Nymphalidae: Satyrinae).

We report a high-quality genome draft assembly of the dark-branded bushbrown, Mycalesis mineus, a member of the Satyrinae subfamily of nymphalid butterflies. This species is emerging as a promising model organism for investigating the evolution and development of phenotypic plasticity. Using 45.99 Gb of long-read data (N50 = 11.11 kb), we assembled a genome size of 497.4 Mb for M. mineus. The assembly is highly contiguous and nearly complete (96.8% of Benchmarking Universal Single-Copy Orthologs lepidopteran genes were complete and single copy). The genome comprises 38.71% of repetitive elements and includes 20,967 predicted protein-coding genes. The assembled genome was super-scaffolded into 28 pseudo-chromosomes using a closely related species, Bicyclus anynana, with a chromosomal-level genome as a template. This valuable genomic tool will advance both ongoing and future research focused on this model organism.

Draft Genome Sequence of Bacillus thuringiensis INTA 103-23 Reveals Its Insecticidal Properties: Insights from the Genomic Sequence

The genome of Bacillus thuringiensis strain INTA 103-23 was sequenced, revealing a high-quality draft assembly comprising 243 contigs with a total size of 6.30 Mb and a completeness of 99%. Phylogenetic analysis classified INTA 103-23 within the Bacillus cereus sensu stricto cluster. Genome annotation identified 6993 genes, including 2476 hypothetical proteins. Screening for pesticidal proteins unveiled 10 coding sequences with significant similarity to known pesticidal proteins, showcasing a potential efficacy against various insect orders. AntiSMASH analysis predicted 13 biosynthetic gene clusters (BGCs), including clusters with 100% similarity to petrobactin and anabaenopeptin NZ857/nostamide A. Notably, fengycin exhibited a 40% similarity within the identified clusters. Further exploration involved a comparative genomic analysis with ten phylogenetically closest genomes. The ANI values, calculated using fastANI, confirmed the closest relationships with strains classified under Bacillus cereus sensu stricto. This comprehensive genomic analysis of B. thuringiensis INTA 103-23 provides valuable insights into its genetic makeup, potential pesticidal activity, and biosynthetic capabilities. The identified BGCs and pesticidal proteins contribute to our understanding of the strain’s biocontrol potential against diverse agricultural pests.

A combination of conserved and diverged responses underlies Theobroma cacao’s defense response to Phytophthora palmivora

BackgroundPlants have complex and dynamic immune systems that have evolved to resist pathogens. Humans have worked to enhance these defenses in crops through breeding. However, many crops harbor only a fraction of the genetic diversity present in wild relatives. Increased utilization of diverse germplasm to search for desirable traits, such as disease resistance, is therefore a valuable step towards breeding crops that are adapted to both current and emerging threats. Here, we examine diversity of defense responses across four populations of the long-generation tree crop Theobroma cacao L., as well as four non-cacao Theobroma species, with the goal of identifying genetic elements essential for protection against the oomycete pathogen Phytophthora palmivora.ResultsWe began by creating a new, highly contiguous genome assembly for the P. palmivora-resistant genotype SCA 6 (Additional file 1: Tables S1-S5), deposited in GenBank under accessions CP139290-CP139299. We then used this high-quality assembly to combine RNA and whole-genome sequencing data to discover several genes and pathways associated with resistance. Many of these are unique, i.e., differentially regulated in only one of the four populations (diverged 40 k–900 k generations). Among the pathways shared across all populations is phenylpropanoid biosynthesis, a metabolic pathway with well-documented roles in plant defense. One gene in this pathway, caffeoyl shikimate esterase (CSE), was upregulated across all four populations following pathogen treatment, indicating its broad importance for cacao’s defense response. Further experimental evidence suggests this gene hydrolyzes caffeoyl shikimate to create caffeic acid, an antimicrobial compound and known inhibitor of Phytophthora spp.ConclusionsOur results indicate most expression variation associated with resistance is unique to populations. Moreover, our findings demonstrate the value of using a broad sample of evolutionarily diverged populations for revealing the genetic bases of cacao resistance to P. palmivora. This approach has promise for further revealing and harnessing valuable genetic resources in this and other long-generation plants.

Nonlinear Tolerancing: Variation Simulation and Assembly Analysis with Regard to Contact Interaction of Parts

The variation analysis is a key tool for ensuring the high quality assembly in the process of developing the technology for manufacturing of aircraft parts. One of the main factors in variations is the deviations in the positioning procedure. This paper is devoted to the development of an approach that allows taking into account the variations during positioning and merging it with the special algorithm of contact problem solving. The impact of varied boundary conditions is incorporated into an additional vector of forces that can be interpreted as reactions to the shift of supports. The obtained results are illustrated with a case of wing-to-fuselage assembly.

The genome of the ant Tetramorium bicarinatum reveals a tandem organization of venom peptides genes allowing the prediction of their regulatory and evolutionary profiles

BackgroundVenoms have evolved independently over a hundred times in the animal kingdom to deter predators and/or subdue prey. Venoms are cocktails of various secreted toxins, whose origin and diversification provide an appealing system for evolutionary researchers. Previous studies of the ant venom of Tetramorium bicarinatum revealed several Myrmicitoxin (MYRTX) peptides that gathered into seven precursor families suggesting different evolutionary origins. Analysis of the T. bicarinatum genome enabling further genomic approaches was necessary to understand the processes underlying the evolution of these myrmicitoxins.ResultsHere, we sequenced the genome of Tetramorium bicarinatum and reported the organisation of 44 venom peptide genes (vpg). Of the eleven chromosomes that make up the genome of T. bicarinatum, four carry the vpg which are organized in tandem repeats. This organisation together with the ML evolutionary analysis of vpg sequences, is consistent with evolution by local duplication of ancestral genes for each precursor family. The structure of the vpg into two or three exons is conserved after duplication events while the promoter regions are the least conserved parts of the vpg even for genes with highly identical sequences. This suggests that enhancer sequences were not involved in duplication events, but were recruited from surrounding regions. Expression level analysis revealed that most vpg are highly expressed in venom glands, although one gene or group of genes is much more highly expressed in each family. Finally, the examination of the genomic data revealed that several genes encoding transcription factors (TFs) are highly expressed in the venom glands. The search for binding sites (BS) of these TFs in the vpg promoters revealed hot spots of GATA sites in several vpg families.ConclusionIn this pioneering investigation on ant venom genes, we provide a high-quality assembly genome and the annotation of venom peptide genes that we think can fosters further genomic research to understand the evolutionary history of ant venom biochemistry.

Haplotype-resolved assembly of a tetraploid potato genome using long reads and low-depth offspring data

Potato is one of the world’s major staple crops, and like many important crop plants, it has a polyploid genome. Polyploid haplotype assembly poses a major computational challenge. We introduce a novel strategy for the assembly of polyploid genomes and present an assembly of the autotetraploid potato cultivar Altus. Our method uses low-depth sequencing data from an offspring population to achieve chromosomal clustering and haplotype phasing on the assembly graph. Our approach generates high-quality assemblies of individual chromosomes with haplotype-specific sequence resolution of whole chromosome arms and can be applied in common breeding scenarios where collections of offspring are available.

Chromosome-Level Assemblies for the Pine Pitch Canker Pathogen Fusarium circinatum.

The pine pitch canker pathogen, Fusarium circinatum, is globally regarded as one of the most important threats to commercial pine-based forestry. Although genome sequences of this fungus are available, these remain highly fragmented or structurally ill-defined. Our overall goal was to provide high-quality assemblies for two notable strains of F. circinatum, and to characterize these in terms of coding content, repetitiveness and the position of telomeres and centromeres. For this purpose, we used Oxford Nanopore Technologies MinION long-read sequences, as well as Illumina short sequence reads. By leveraging the genomic synteny inherent to F. circinatum and its close relatives, these sequence reads were assembled to chromosome level, where contiguous sequences mostly spanned from telomere to telomere. Comparative analyses unveiled remarkable variability in the twelfth and smallest chromosome, which is known to be dispensable. It presented a striking length polymorphism, with one strain lacking substantial portions from the chromosome's distal and proximal regions. These regions, characterized by a lower gene density, G+C content and an increased prevalence of repetitive elements, contrast starkly with the syntenic segments of the chromosome, as well as with the core chromosomes. We propose that these unusual regions might have arisen or expanded due to the presence of transposable elements. A comparison of the overall chromosome structure revealed that centromeric elements often underpin intrachromosomal differences between F. circinatum strains, especially at chromosomal breakpoints. This suggests a potential role for centromeres in shaping the chromosomal architecture of F. circinatum and its relatives. The publicly available genome data generated here, together with the detailed metadata provided, represent essential resources for future studies of this important plant pathogen.

Chromosome-Level Genome Assembly and Annotation of a Periodical Cicada Species: Magicicada septendecula.

We present a high-quality assembly and annotation of the periodical cicada species, Magicicada septendecula (Hemiptera: Auchenorrhyncha: Cicadidae). Periodical cicadas have a significant ecological impact, serving as a food source for many mammals, reptiles, and birds. Magicicada are well known for their massive emergences of 1 to 3 species that appear in different locations in the eastern United States nearly every year. These year classes ("broods") emerge dependably every 13 or 17 yr in a given location. Recently, it has become clear that 4-yr early or late emergences of a sizeable portion of a population are an important part of the history of brood formation; however, the biological mechanisms by which they track the passage of time remain a mystery. Using PacBio HiFi reads in conjunction with Hi-C proximity ligation data, we have assembled and annotated the first whole genome for a periodical cicada, an important resource for future phylogenetic and comparative genomic analysis. This also represents the first quality genome assembly and annotation for the Hemipteran superfamily Cicadoidea. With a scaffold N50 of 518.9 Mb and a complete BUSCO score of 96.7%, we are confident that this assembly will serve as a vital resource toward uncovering the genomic basis of periodical cicadas' long, synchronized life cycles and will provide a robust framework for further investigations into these insects.

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

High-quality metagenome assembly from long accurate reads with metaMDBG

We introduce metaMDBG, a metagenomics assembler for PacBio HiFi reads. MetaMDBG combines a de Bruijn graph assembly in a minimizer space with an iterative assembly over sequences of minimizers to address variations in genome coverage depth and an abundance-based filtering strategy to simplify strain complexity. For complex communities, we obtained up to twice as many high-quality circularized prokaryotic metagenome-assembled genomes as existing methods and had better recovery of viruses and plasmids.

High-quality assembly of the T2T genome for Isodon rubescens f. lushanensis reveals genomic structure variations between 2 typical forms of Isodon rubescens.

Rabdosiae rubescentis herba (Isodon rubescens) is widely used as a folk medicine to treat esophageal cancer and sore throat in China. Its germplasm resources are abundant in China, with I. rubescens (Hemsl.) Hara and I. rubescens f. lushanensis as 2 typical forms. I. rubescens (Hemsl.) Hara is featured by biosynthesis of the diterpenoid oridonin with strong anticancer activity, while I. rubescens f. lushanensis produces another diterpenoid with anticancer activity, lushanrubescensin. However, the biosynthetic pathways of both still need to be fully understood. In particular, little is known about the genetic background of I. rubescens f. lushanensis. We used Pacific Biosciences (PacBio) single-molecule real-time and Nanopore Ultra-long sequencing platforms, respectively, and obtained 139.07 Gb of high-quality data, with a sequencing depth of about 328×. We also obtained a high-quality reference genome for I. rubescens f. lushanensis, with a genome size of 349 Mb and a contig N50 of 28.8 Mb. The heterozygosity of the genome is 1.7% and the repeatability is 83.43%. In total, 34,865 protein-coding genes were predicted. Moreover, we found that most of the variant or unique genes in the diterpenoid synthesis pathways of I. rubescens f. lushanensis and I. rubescens (Hemsl.) Hara were enriched in diterpene synthases. We provide the first genome sequence and gene annotation for the I. rubescens f. lushanensis, which provides molecular evidence for understanding the chemotypic differences of I. rubescens.

DNA Isolation Long-Read Genomic Sequencing in Ctenophores.

Long-read sequencing has proven the necessity for high-quality genomic assemblies of reference species, including enigmatic ctenophores. Obtaining high-molecular-weight genomic DNA is pivotal to this process and has proven highly problematic for many species. Here, we discuss different methodologies for gDNA isolation and present a protocol for isolating gDNA for several members of the phylum Ctenophora. Specifically, we describe a Pacific Biosciences library construction method used in conjunction with gDNA isolation methods that have proven successful in obtaining high-quality genomic assemblies in ctenophores.

A high-quality whole-genome sequence, assembly, and gene annotation of Fusarium oxysporum f. sp. vasinfectum (Fov) race 1 from California.

Fusarium wilt [Fusarium oxysporum f. sp. vasinfectum (FOV)] in cotton is a widespread soilborne pathogen that causes vascular plant disease and is responsible for substantial crop losses worldwide. FOV race 1 (FOV1) is well established across almost all cotton production regions, especially in the USA. Herein, we report a high-quality whole-genome sequence, assembly, and gene annotation of a FOV1 isolate from California.