Total Assembly Length Research Articles

Assessing the de novo assemblers: a metaviromic study of apple and first report of citrus concave gum-associated virus, apple rubbery wood virus 1 and 2 infecting apple in India

BackgroundThe choice of de novo assembler for high-throughput sequencing (HTS) data remains a pivotal factor in the HTS-based discovery of viral pathogens. This study assessed de novo assemblers, namely Trinity, SPAdes, and MEGAHIT for HTS datasets generated on the Illumina platform from 23 apple samples, representing 15 exotic and indigenous apple varieties and a rootstock. The assemblers were compared based on assembly quality metrics, including the largest contig, total assembly length, genome coverage, and N50.ResultsMEGAHIT was most efficient assembler according to the metrics evaluated in this study. By using multiple assemblers, near-complete genome sequences of citrus concave gum-associated virus (CCGaV), apple rubbery wood virus 1 (ARWV-1), ARWV-2, apple necrotic mosaic virus (ApNMV), apple mosaic virus, apple stem pitting virus, apple stem grooving virus, apple chlorotic leaf spot virus, apple hammerhead viroid and apple scar skin viroid were reconstructed. These viruses were further confirmed through Sanger sequencing in different apple cultivars. Among them, CCGaV, ARWV-1 and ARWV-2 were recorded from apples in India for the first time. The analysis of virus richness revealed that ApNMV was dominant, followed by ARWV-1 and CCGaV. Moreover, MEGAHIT identified novel single-nucleotide variants.ConclusionsOur analyses highlight the crucial role of assembly methods in reconstructing near-complete apple virus genomes from the Illumina reads. This study emphasizes the significance of employing multiple assemblers for de novo virus genome assembly in vegetatively propagated perennial fruit crops.

Chromosome-level genome provides insights into evolution and diving adaptability in the vulnerable common pochard (Aythya ferina)

The common pochard (Aythya ferina) is a freshwater diving duck found in the Palearctic region that has been classified as vulnerable by the IUCN due to continuous and rapid population declines across their distribution. To gain a better understanding of its genetic mechanism of adaptive evolution, we successfully sequenced and assembled the first high-quality chromosome-level genome of A. ferina using Illumina, Nanopore and Hi-C sequencing technologies. A total assembly length of 1,130.78 Mbp was obtained, with over 98.81% (1,117.37Mbp) of sequence anchored to 35 pseudo-chromosomes. We predicted 17,232 protein-coding genes, 95.9% of which were functionally annotated. We identified 339 expanded and 937 contracted gene families in the genome of A. ferina, and detected 95 genes that have been positively selected. The significantly enriched Gene Ontology and enriched pathways were related to energy metabolism, immune, nervous, and sensory systems, suggests that these factors likely played an important role in its evolution. Importantly, we recovered signatures of positive selection on genes related to vasoconstriction that may be associated with thermoregulatory adaptations of A. ferina for underwater diving. Overall, the high-quality genome assembly and annotation in this study provides valuable genomic resources for ecological and evolutionary studies, as well as toward the conservation of A. ferina.

Chromosome-scale Genome Assembly of the Rough Periwinkle Littorina saxatilis.

The intertidal gastropod Littorina saxatilis is a model system to study speciation and local adaptation. The repeated occurrence of distinct ecotypes showing different levels of genetic divergence makes L. saxatilis particularly suited to study different stages of the speciation continuum in the same lineage. A major finding is the presence of several large chromosomal inversions associated with the divergence of ecotypes and, specifically, the species offers a system to study the role of inversions in this divergence. The genome of L. saxatilis is 1.35 Gb and composed of 17 chromosomes. The first reference genome of the species was assembled using Illumina data, was highly fragmented (N50 of 44 kb), and was quite incomplete, with a BUSCO completeness of 80.1% on the Metazoan dataset. A linkage map of one full-sibling family enabled the placement of 587 Mbp of the genome into 17 linkage groups corresponding to the haploid number of chromosomes, but the fragmented nature of this reference genome limited the understanding of the interplay between divergent selection and gene flow during ecotype formation. Here, we present a newly generated reference genome that is highly contiguous, with a N50 of 67 Mb and 90.4% of the total assembly length placed in 17 super-scaffolds. It is also highly complete with a BUSCO completeness of 94.1% of the Metazoa dataset. This new reference will allow for investigations into the genomic regions implicated in ecotype formation as well as better characterization of the inversions and their role in speciation.

Methylocystis suflitae sp. nov., a novel type II methanotrophic bacterium isolated from landfill cover soil.

A bacterial strain, designated NLS-7T, was isolated through enrichment of landfill cover soil in methane-oxidizing conditions. Strain NLS-7T is a Gram-stain negative, non-motile rod, approximately 0.8 µm wide by 1.3 µm long. Phylogenetic analysis based on 16S rRNA gene sequencing places it within the genus Methylocystis, with its closest relatives being M. hirsuta, M. silviterrae and M. rosea, with 99.9, 99.7 and 99.6 % sequence similarity respectively. However, average nucleotide identity and average amino acid identity values below the 95 % threshold compared to all the close relatives and digital DNA-DNA hybridization values between 20.9 and 54.1 % demonstrate that strain NLS-7T represents a novel species. Genome sequencing generated 4.31 million reads and genome assembly resulted in the generation of 244 contigs with a total assembly length of 3 820 957 bp (N50, 37 735 bp; L50, 34). Genome completeness is 99.5 % with 3.98 % contamination. It is capable of growth on methane and methanol. It grows optimally at 30 °C between pH 6.5 and 7.0. Strain NLS-7T is capable of atmospheric dinitrogen fixation and can use ammonium (as NH4Cl), l-aspartate, l-arginine, yeast extract, nitrate, l-leucine, l-proline, l-methionine, l-lysine and l-alanine as nitrogen sources. The major fatty acids are C18:1 ω8c and C18:1 ω7c. Based upon this polyphasic taxonomic study, strain NLS-7T represents a novel species of the genus Methylocystis, for which the name Methylocystis suflitae sp. nov. is proposed. The type strain is NLS-7T (=ATCC TSD-256T=DSM 112294T). The 16S rRNA gene and genome sequences of strain NLS-7T have been deposited in GenBank under accession numbers ON715489 and GCA_024448135.1, respectively.

High-speed whole-genome sequencing of a Whippet: Rapid chromosome-level assembly and annotation of an extremely fast dog's genome.

The time required for genome sequencing and de novo assembly depends on the interaction between laboratory work, sequencing capacity, and the bioinformatics workflow, often constrained by external sequencing services. Bringing together academic biodiversity institutes and a medical diagnostics company with extensive sequencing capabilities, we aimed at generating a high-quality mammalian de novo genome in minimal time. We present the first chromosome-level genome assembly of the Whippet, using PacBio long-read high-fidelity sequencing and reference-guided scaffolding. The final assembly has a contig N50 of 55 Mbp and a scaffold N50 of 65.7 Mbp. The total assembly length is 2.47 Gbp, of which 2.43 Gpb were scaffolded into 39 chromosome-length scaffolds. Annotation using mammalian genomes and transcriptome data yielded 28,383 transcripts, 90.9% complete BUSCO genes, and identified 36.5% repeat content. Sequencing, assembling, and scaffolding the chromosome-level genome of the Whippet took less than a week, adding another high-quality reference genome to the available sequences of domestic dog breeds.

Comparative Transcriptome Analysis of the Skin and the Peritoneal Wall Layer of Triplophysa stenura Distributed in High Elevations.

A total of 81,868 All-Unigenes were sequenced and assembled by the transcriptome in the dorsal skin, the lateral skin, and the peritoneal wall layer of Triplophysa stenura with a total assembly length of 123,827,585 bp, and 68,750 unigenes were annotated to seven functional databases. A total of 588 DEGs were screened between the dorsal and lateral skin, 17,097 DEGs were screened between the dorsal skin and the peritoneal wall layer, and 16,598 DEGs were screened between the lateral skin and the peritoneal wall layer. Most of DEGs in three tissues were annotated to GO terms related to cellular structures, binding, cellular processes, and catalytic activity. They were also annotated to KEGG pathways such as the MAPK signaling pathway, PI3K-Akt signaling pathway, Wnt signaling pathway, melanogenesis, tyrosine metabolism, and cell cycle. A total of twenty-three DEGs were found to be enriched in the melanin synthesis pathway by a local Blast comparison, of which nine DEGs were significantly upregulated in the peritoneal wall layer and six DEGs were significantly upregulated in the dorsal and lateral skin. The results suggest that these genes may be associated with the molecular mechanism of melanin synthesis in T. stenura, and the differential regulation of genes may be related to the differences in UVR intensity and tissue sites of melanin synthesis. Further investigation is needed on how these genes specifically regulate melanin synthesis.

Chromosome-level genome assembly of Chouioia cunea Yang, the parasitic wasp of the fall webworm

Chouioia cunea Yang 1989 is a parasitic wasp of many lepidopteran insects during their pupal stage, and has been successfully used to control pests such as the fall webworm Hyphantria cunea. Here we reported the chromosome-level genome of C. cunea by using short (MGI-SEQ), long (Oxford Nanopore), chromatin-linked (Hi-C) sequencing reads and transcriptomic data, representing the first chromosome-level genome of parasitic wasps of the family Eulophidae. The total assembly length is 171.99 Mb, containing 6 pesudo-chromosomes with a GC content of 36.89% and the scaffold/contig N50 length of 31.70/26.52 Mb. The BUSCO completeness of the assembly was estimated to be 98.7%. A total of 12,258 protein-coding genes (PCGs), 10,547 3′-UTRs, and 10,671 5′-UTRs were annotated. This high-quality genome is an important step toward a better understanding of the genomes of the Eulophidae (Chalcidoidea), and will serve as a valuable resource for analyses of phylogenetic relationships and the evolution of Hymenoptera.

Benchmarking genome assembly methods on metagenomic sequencing data.

Metagenome assembly is an efficient approach to reconstruct microbial genomes from metagenomic sequencing data. Although short-read sequencing has been widely used for metagenome assembly, linked- and long-read sequencing have shown their advancements in assembly by providing long-range DNA connectedness. Many metagenome assembly tools were developed to simplify the assembly graphs and resolve the repeats in microbial genomes. However, there remains no comprehensive evaluation of metagenomic sequencing technologies, and there is a lack of practical guidance on selecting the appropriate metagenome assembly tools. This paper presents a comprehensive benchmark of 19 commonly used assembly tools applied to metagenomic sequencing datasets obtained from simulation, mock communities or human gut microbiomes. These datasets were generated using mainstream sequencing platforms, such as Illumina and BGISEQ short-read sequencing, 10x Genomics linked-read sequencing, and PacBio and Oxford Nanopore long-read sequencing. The assembly tools were extensively evaluated against many criteria, which revealed that long-read assemblers generated high contig contiguity but failed to reveal some medium- and high-quality metagenome-assembled genomes (MAGs). Linked-read assemblers obtained the highest number of overall near-complete MAGs from the human gut microbiomes. Hybrid assemblers using both short- and long-read sequencing were promising methods to improve both total assembly length and the number of near-complete MAGs. This paper also discussed the running time and peak memory consumption of these assembly tools and provided practical guidance on selecting them.

Draft Genome Sequence of Plesiomonas shigelloides Strain zfcc0051 (Phylum Proteobacteria).

ABSTRACTHere, we report a draft genome sequence of Plesiomonas shigelloides strain zfcc0051, an isolate derived from zebrafish (Danio rerio) feces. The genome consists of 115 contigs (>500 bp) and has a total assembly length of 4,041,537 bases.

Improvements in the sequencing and assembly of plant genomes.

Advances in DNA sequencing have made it easier to sequence and assemble plant genomes. Here, we extend an earlier study, and compare recent methods for long read sequencing and assembly. Updated Oxford Nanopore Technology software improved assemblies. Using more accurate sequences produced by repeated sequencing of the same molecule (Pacific Biosciences HiFi) resulted in less fragmented assembly of sequencing reads. Using data for increased genome coverage resulted in longer contigs, but reduced total assembly length and improved genome completeness. The original model species, Macadamia jansenii, was also compared with three other Macadamia species, as well as avocado (Persea americana) and jojoba (Simmondsia chinensis). In these angiosperms, increasing sequence data volumes caused a linear increase in contig size, decreased assembly length and further improved already high completeness. Differences in genome size and sequence complexity influenced the success of assembly. Advances in long read sequencing technology continue to improve plant genome sequencing and assembly. However, results were improved by greater genome coverage, with the amount needed to achieve a particular level of assembly being species dependent.

Chromosome-Level Reference Genome Assembly for the American Pika (Ochotona princeps).

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.

Genome sequence of the coffee root-knot nematode Meloidogyne exigua.

Root-knot nematodes (Meloidogyne spp.) cause serious damages on most crops. Here, we report a high-quality genome sequence of Meloidogyne exigua (population Mex1, Costa Rica), a major pathogen of coffee. Its mitogenome (20,974 bp) was first assembled and annotated. The nuclear genome was then constructed consisting of 206 contigs, with an N50 length of 1.89 Mb and a total assembly length of 42.1 Mb.

Chromosome-level genome assembly of a benthic associated Syngnathiformes species: the common dragonet, Callionymus lyra.

The common dragonet, Callionymus lyra, is one of three Callionymus species inhabiting the North Sea. All three species show strong sexual dimorphism. The males show strong morphological differentiation, e.g., species-specific colouration and size relations, while the females of different species have few distinguishing characters. Callionymus belongs to the 'benthic associated clade' of the order Syngnathiformes. The 'benthic associated clade' so far is not represented by genome data and serves as an important outgroup to understand the morphological transformation in 'long-snouted' syngnatiformes such as seahorses and pipefishes. Here, we present the chromosome-level genome assembly of C.lyra. We applied Oxford Nanopore Technologies' long-read sequencing, short-read DNBseq, and proximity-ligation-based scaffolding to generate a high-quality genome assembly. The resulting assembly has a contig N50 of 2.2 Mbp and a scaffold N50 of 26.7 Mbp. The total assembly length is 568.7 Mbp, of which over 538 Mbp were scaffolded into 19 chromosome-length scaffolds. The identification of 94.5% complete BUSCO genes indicates high assembly completeness. Additionally, we sequenced and assembled a multi-tissue transcriptome with a total length of 255.5 Mbp that was used to aid the annotation of the genome assembly. The annotation resulted in 19,849 annotated transcripts and identified a repeat content of 27.7%. The chromosome-level assembly of C.lyra provides a high-quality reference genome for future population genomic, phylogenomic, and phylogeographic analyses.

Regional sequence expansion or collapse in heterozygous genome assemblies.

High levels of heterozygosity present a unique genome assembly challenge and can adversely impact downstream analyses, yet is common in sequencing datasets obtained from non-model organisms. Here we show that by re-assembling a heterozygous dataset with variant parameters and different assembly algorithms, we are able to generate assemblies whose protein annotations are statistically enriched for specific gene ontology categories. While total assembly length was not significantly affected by assembly methodologies tested, the assemblies generated varied widely in fragmentation level and we show local assembly collapse or expansion underlying the enrichment or depletion of specific protein functional groups. We show that these statistically significant deviations in gene ontology groups can occur in seemingly high-quality assemblies, and result from difficult-to-detect local sequence expansion or contractions. Given the unpredictable interplay between assembly algorithm, parameter, and biological sequence data heterozygosity, we highlight the need for better measures of assembly quality than N50 value, including methods for assessing local expansion and collapse.

Genome sequence of the root-knot nematode Meloidogyne luci.

Root-knot nematodes from the genus Meloidogyne are polyphagous plant endoparasites and agricultural pests of global importance. Here, we report the high-quality genome sequence of Meloidogyne luci population SI-Smartno V13. The resulting genome assembly of M. luci SI-Smartno V13 consists of 327 contigs, with an N50 contig length of 1,711,905 bp and a total assembly length of 209.16 Mb.

First Insights into the Guar (Cyamopsis tetragonoloba (L.) Taub.) Genome of the ‘Vavilovskij 130’ Accession, Using Second and Third-Generation Sequencing Technologies

Guar (Cyamopsis tetragonoloba (L.) Taub.) is becoming a popular industrial crop in response to industry demand for the guar gum extracted from seeds’ endosperm. Breeding efforts of new guar varieties would greatly benefit from genomic resources developed for marker assisted selection (MAS) purposes. We have undertaken the first steps to establish a whole-genome assembly of the guar ‘Vaviloskij 130’ accession, bred at VIR. Using a combination of second (Illumina short reads) and third generation (Oxford Nanopore long reads) sequencing methods, a dataset of approx. 5X of genome coverage was obtained. We tested assemblers for short reads, namely SOAPdenovo, AbySS and SGA, based on different algorithms. For short reads (Illumina MiSeq and HiSeq data), the better result in terms of total number of scaffolds and total assembly length were obtained with SGA (String Graph Assembler). For Oxford Nanopore dataset, we used the combination of minimap + miniASM assembly, then corrected the assembly with raw Illumina and Nanopore data. The current preliminary de novo assembly of the guar genome covers 1.2 Gb, corresponding to 50% of the genome. The data confirm the phylogeny position of C. tetragoloba as being highly related to the genus Vigna, Abrus, Glycine and Lupinus genomes. This preliminary reference genome paves the way to further detailed diversity and genetic analyses into that important agro-industrial crop.

Draft genome sequence of a clinical Acinetobacter haemolyticus isolate from South Africa.

Here we describe the draft genome sequence of a clinical Acinetobacter haemolyticus isolate (A109R1B4) from a rectal swab of a hospitalised patient in South Africa. Genomic DNA from the isolate was sequenced using an Illumina MiSeq platform. Generated reads were de novo assembled using Qiagen CLC Genomics Workbench. The assembled contigs were annotated and antimicrobial resistance genes and sequence type were identified. The genome comprised 281 contigs with a total assembly length of 3 371 389bp, a G+C content of 39.9% and an N50 value of 58 196bp, and reference coverage was 93.08 while the coverage was 88.08. A total of 3387 genes, 3292 coding sequences (CDS), 3175 coding genes and 95 RNA genes were detected. The antimicrobial resistance genes blaOXA-264 and aac(6')-Ig were detected. The genome sequence reported will serve as a reference point for molecular epidemiological studies of antibiotic-resistant A. haemolyticus in Africa.

High-Quality Draft Genome Sequence of Fusarium oxysporum f. sp. cubense Strain 160527, a Causal Agent of Panama Disease.

Fusarium oxysporum f. sp. cubense is the causal agent of banana Fusarium wilt, also known as Panama disease. Here, we present a high-quality genome sequence of F. oxysporum f. sp. cubense strain 160527. The genome assembly is composed of 12 contigs with a total assembly length of 51,139,495 bp (N 50 contig length, 4,884,632 bp).

Hypoxic and Cold Adaptation Insights from the Himalayan Marmot Genome

SummaryThe Himalayan marmot (Marmota himalayana) is a hibernating mammal that inhabits the high-elevation regions of the Himalayan mountains. Here we present a draft genome of the Himalayan marmot, with a total assembly length of 2.47 Gb. Phylogenetic analyses showed that the Himalayan marmot diverged from the Mongolian marmot approximately 1.98 million years ago. Transcriptional changes during hibernation included genes responsible for fatty acid metabolism in liver and genes involved in complement and coagulation cascades and stem cell pluripotency pathways in brain. Two selective sweep genes, Slc25a14 and ψAamp, showed apparent genotyping differences between low- and high-altitude populations. As a processed pseudogene, ψAamp may be biologically active to influence the stability of Aamp through competitive microRNA binding. These findings shed light on the molecular and genetic basis underlying adaptation to extreme environments in the Himalayan marmot.

Genome sequencing and conservation genomics in the Scandinavian wolverine population.

Genetic approaches have proved valuable to the study and conservation of endangered populations, especially for monitoring programs, and there is potential for further developments in this direction by extending analyses to the genomic level. We assembled the genome of the wolverine (Gulo gulo), a mustelid that in Scandinavia has recently recovered from a significant population decline, and obtained a 2.42Gb draft sequence representing >85% of the genome and including >21,000protein-coding genes. We then performed whole-genome resequencing of 10Scandinavian wolverines for population genomic and demographic analyses. Genetic diversity was among the lowest detected in a red-listed population (mean genome-wide nucleotide diversity of 0.05%). Results of the demographic analyses indicated a long-term decline of the effective population size (Ne ) from 10,000well before the last glaciation to <500 after this period. Current Ne appeared even lower. The genome-wide FIS level was 0.089 (possibly signaling inbreeding), but this effect was not observed when analyzing a set of highly variable SNP markers, illustrating that such markers can give a biased picture of the overall character of genetic diversity. We found significant population structure, which has implications for population connectivity and conservation. We used an integrated microfluidic circuit chip technology to develop an SNP-array consisting of 96 highly informative markers that, together with a multiplex pre-amplification step, was successfully applied to low-quality DNA from scat samples. Our findings will inform management, conservation, and genetic monitoring of wolverines and serve as a genomic roadmap that can be applied to other endangered species. The approach used here can be generally utilized in other systems, but we acknowledge the trade-off between investing in genomic resources and direct conservation actions.