SLR: a scaffolding algorithm based on long reads and contig classification

Phased chromosome-scale genome assembly of an asexual, allopolyploid root-knot nematode reveals complex subgenomic structure

We present the chromosome-scale genome assembly of the allopolyploid root-knot nematode Meloidogyne javanica. We show that the M. javanica genome is predominantly allotetraploid, comprising two subgenomes, A and B, that most likely originated from hybridisation of two ancestral parental species. The assembly was annotated using full-length non-chimeric transcripts, comparison to reference databases, and ab initio prediction techniques, and the subgenomes were phased using ancestral k-mer spectral analysis. Subgenome B appears to show fission of chromosomal contigs, and while there is substantial synteny between subgenomes, we also identified regions lacking synteny that may have diverged in the ancestral genomes prior to or following hybridisation. This annotated and phased genome assembly forms a significant resource for understanding the origins and genetics of these globally important plant pathogens.

Geography, not lifestyle, explains the population structure of free-living and host-associated deep-sea hydrothermal vent snail symbionts

BACKGROUND

Marine symbioses are predominantly established through horizontal acquisition of microbial symbionts from the environment. However, genetic and functional comparisons of free-living populations of symbionts to their host-associated counterparts are sparse. Here, we assembled the first genomes of the chemoautotrophic gammaproteobacterial symbionts affiliated with the deep-sea snail Alviniconcha hessleri from two separate hydrothermal vent fields of the Mariana Back-Arc Basin. We used phylogenomic and population genomic methods to assess sequence and gene content variation between free-living and host-associated symbionts.

RESULTS

Our phylogenomic analyses show that the free-living and host-associated symbionts of A. hessleri from both vent fields are populations of monophyletic strains from a single species. Furthermore, genetic structure and gene content analyses indicate that these symbiont populations are differentiated by vent field rather than by lifestyle.

CONCLUSION

Together, this work suggests that, despite the potential influence of host-mediated acquisition and release processes on horizontally transmitted symbionts, geographic isolation and/or adaptation to local habitat conditions are important determinants of symbiont population structure and intra-host composition. Video Abstract.

SLHSD: hybrid scaffolding method based on short and long reads

In genome assembly, scaffolding can obtain more complete and continuous scaffolds. Current scaffolding methods usually adopt one type of read to construct a scaffold graph and then orient and order contigs. However, scaffolding with the strengths of two or more types of reads seems to be a better solution to some tricky problems. Combining the advantages of different types of data is significant for scaffolding. Here, a hybrid scaffolding method (SLHSD) is present that simultaneously leverages the precision of short reads and the length advantage of long reads. Building an optimal scaffold graph is an important foundation for getting scaffolds. SLHSD uses a new algorithm that combines long and short read alignment information to determine whether to add an edge and how to calculate the edge weight in a scaffold graph. In addition, SLHSD develops a strategy to ensure that edges with high confidence can be added to the graph with priority. Then, a linear programming model is used to detect and remove remaining false edges in the graph. We compared SLHSD with other scaffolding methods on five datasets. Experimental results show that SLHSD outperforms other methods. The open-source code of SLHSD is available at https://github.com/luojunwei/SLHSD.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

Cave Thiovulum (Candidatus Thiovulum stygium) differs metabolically and genomically from marine species

Thiovulum spp. (Campylobacterota) are large sulfur bacteria that form veil-like structures in aquatic environments. The sulfidic Movile Cave (Romania), sealed from the atmosphere for ~5 million years, has several aqueous chambers, some with low atmospheric O₂ (~7%). The cave's surface-water microbial community is dominated by bacteria we identified as Thiovulum. We show that this strain, and others from subsurface environments, are phylogenetically distinct from marine Thiovulum. We assembled a closed genome of the Movile strain and confirmed its metabolism using RNAseq. We compared the genome of this strain and one we assembled from public data from the sulfidic Frasassi caves to four marine genomes, including Candidatus Thiovulum karukerense and Ca. T. imperiosus, whose genomes we sequenced. Despite great spatial and temporal separation, the genomes of the Movile and Frasassi Thiovulum were highly similar, differing greatly from the very diverse marine strains. We concluded that cave Thiovulum represent a new species, named here Candidatus Thiovulum stygium. Based on their genomes, cave Thiovulum can switch between aerobic and anaerobic sulfide oxidation using O₂ and NO₃^- as electron acceptors, the latter likely via dissimilatory nitrate reduction to ammonia. Thus, Thiovulum is likely important to both S and N cycles in sulfidic caves. Electron microscopy analysis suggests that at least some of the short peritrichous structures typical of Thiovulum are type IV pili, for which genes were found in all strains. These pili may play a role in veil formation, by connecting adjacent cells, and in the motility of these exceptionally fast swimmers.

Genome of Lindera glauca provides insights into the evolution of biosynthesis genes for aromatic compounds

Lindera glauca is a crucial source of diverse industrial oil and medicines. The spicy aroma of tender leaves is caused by the presence of abundant aromatic compounds. Here, we present its chromosome-level genome assembly comprising 12 pseudochromosomes (2,092.2 Mb; scaffold N50: 186.5 Mb), which was predicted to have 65,145 protein-coding genes. Comparative genomic analyses indicated two whole-genome duplication (WGD) events in the Lauraceae family, contributing to the production of numerous terpene synthase (TPS) genes. We identified 138 TPS genes in L. glauca. Comparative transcriptomic analyses revealed high expression of genes Lg03G2346 and Lg08G140 in TPS-a and Lg07G2961 and Lg12G971 in TPS-b subfamilies, which regulated the biosynthesis of the monoterpenoid β-ocimene and sesquiterpenoid D-germacrene in L. glauca. The results suggested a molecular basis for species-specific terpenoid biosynthesis and provided a foundation for molecular breeding to produce desired characteristics and a valuable reference genome.

Flemish soils contain rhizobia partners for Northwestern Europe-adapted soybean cultivars

In Europe, soybean (Glycine max) used for food and feed has to be imported, causing negative socioeconomic and environmental impacts. To increase the local production, breeding generated varieties that grow in colder climates, but the yield using the commercial inoculants is not satisfactory in Belgium because of variable nodulation efficiencies. To look for indigenous nodulating strains possibly adapted to the local environment, we initiated a nodulation trap by growing early-maturing cultivars under natural and greenhouse conditions in 107 garden soils in Flanders. Nodules occurred in 18 and 21 soils in the garden and greenhouse experiments respectively. By combining 16S rRNA PCR on single isolates with HiSeq 16S metabarcoding on nodules, we found a large bacterial richness and diversity from different soils. Furthermore, using Oxford Nanopore Technologies sequencing of DNA from one nodule, we retrieved the entire genome of a Bradyrhizobium species, not previously isolated, but profusely present in that nodule. These data highlight the need of combining diverse identification techniques to capture the true nodule rhizobial community. Eight selected rhizobial isolates were subdivided by whole-genome analysis in three genera containing six genetically distinct species that, except for two, aligned with known type strains and were all able to nodulate soybean in the laboratory.

ONT-Based Draft Genome Assembly and Annotation of Alternaria atra

Species of Alternaria (phylum Ascomycota, family Pleosporaceae) are known as serious plant pathogens, causing major losses on a wide range of crops. Alternaria atra (previously known as Ulocladium atrum) can grow as a saprophyte on many hosts and causes Ulocladium blight on potato. It has been reported that it can also be used as a biocontrol agent against Botrytis cinerea. Here, we present a scaffold-level reference genome assembly for A. atra. The assembly contains 43 scaffolds with a total length of 39.62 Mbp, with scaffold N₅₀ of 3,893,166 bp, L₅₀ of 4, and the longest 10 scaffolds containing 89.9% of the assembled data. RNA-sequencing-guided gene prediction using BRAKER resulted in 12,173 protein-coding genes with their functional annotation. This first high-quality reference genome assembly and annotation for A. atra can be used as a resource for studying evolution in the highly complicated Alternaria genus and might help in understanding the mechanisms defining its role as pathogen or biocontrol agent.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Macrosynteny analysis between Lentinula edodes and Lentinula novae-zelandiae reveals signals of domestication in Lentinula edodes

The basidiomycete fungus Lentinula novae-zelandiae is endemic to New Zealand and is a sister taxon to Lentinula edodes, the second most cultivated mushroom in the world. To explore the biology of this organism, a high-quality chromosome level reference genome of L. novae-zelandiae was produced. Macrosyntenic comparisons between the genome assembly of L. novae-zelandiae, L. edodes and a set of three genome assemblies of diverse species from the Agaricomycota reveal a high degree of macrosyntenic restructuring within L. edodes consistent with signal of domestication. These results show L. edodes has undergone significant genomic change during the course of its evolutionary history, likely a result of its cultivation and domestication over the last 1000 years.

A comprehensive review of scaffolding methods in genome assembly

In the field of genome assembly, scaffolding methods make it possible to obtain a more complete and contiguous reference genome, which is the cornerstone of genomic research. Scaffolding methods typically utilize the alignments between contigs and sequencing data (reads) to determine the orientation and order among contigs and to produce longer scaffolds, which are helpful for genomic downstream analysis. With the rapid development of high-throughput sequencing technologies, diverse types of reads have emerged over the past decade, especially in long-range sequencing, which have greatly enhanced the assembly quality of scaffolding methods. As the number of scaffolding methods increases, biology and bioinformatics researchers need to perform in-depth analyses of state-of-the-art scaffolding methods. In this article, we focus on the difficulties in scaffolding, the differences in characteristics among various kinds of reads, the methods by which current scaffolding methods address these difficulties, and future research opportunities. We hope this work will benefit the design of new scaffolding methods and the selection of appropriate scaffolding methods for specific biological studies.

Correction to: SLR: a scaffolding algorithm based on long reads and contig classification

Following publication of the original article [1], the author reported that there is an error in the original article.

TGS-GapCloser: A fast and accurate gap closer for large genomes with low coverage of error-prone long reads

BACKGROUND

Analyses that use genome assemblies are critically affected by the contiguity, completeness, and accuracy of those assemblies. In recent years single-molecule sequencing techniques generating long-read information have become available and enabled substantial improvement in contig length and genome completeness, especially for large genomes (>100 Mb), although bioinformatic tools for these applications are still limited.

FINDINGS

We developed a software tool to close sequence gaps in genome assemblies, TGS-GapCloser, that uses low-depth (∼10×) long single-molecule reads. The algorithm extracts reads that bridge gap regions between 2 contigs within a scaffold, error corrects only the candidate reads, and assigns the best sequence data to each gap. As a demonstration, we used TGS-GapCloser to improve the scaftig NG50 value of 3 human genome assemblies by 24-fold on average with only ∼10× coverage of Oxford Nanopore or Pacific Biosciences reads, covering with sequence data up to 94.8% gaps with 97.7% positive predictive value. These improved assemblies achieve 99.998% (Q46) single-base accuracy with final inserted sequences having 99.97% (Q35) accuracy, despite the high raw error rate of single-molecule reads, enabling high-quality downstream analyses, including up to a 31-fold increase in the scaftig NGA50 and up to 13.1% more complete BUSCO genes. Additionally, we show that even in ultra-large genome assemblies, such as the ginkgo (∼12 Gb), TGS-GapCloser can cover 71.6% of gaps with sequence data.

CONCLUSIONS

TGS-GapCloser can close gaps in large genome assemblies using raw long reads quickly and cost-effectively. The final assemblies generated by TGS-GapCloser have improved contiguity and completeness while maintaining high accuracy. The software is available at https://github.com/BGI-Qingdao/TGS-GapCloser.

The honeysuckle genome provides insight into the molecular mechanism of carotenoid metabolism underlying dynamic flower coloration

Lonicera japonica is a widespread member of the Caprifoliaceae (honeysuckle) family utilized in traditional medical practices. This twining vine honeysuckle also is a much-sought ornamental, in part due to its dynamic flower coloration, which changes from white to gold during development. The molecular mechanism underlying dynamic flower coloration in L. japonica was elucidated by integrating whole genome sequencing, transcriptomic analysis and biochemical assays. Here, we report a chromosome-level genome assembly of L. japonica, comprising nine pseudochromosomes with a total size of 843.2 Mb. We also provide evidence for a whole-genome duplication event in the lineage leading to L. japonica, which occurred after its divergence from Dipsacales and Asterales. Moreover, gene expression analysis not only revealed correlated expression of the relevant biosynthetic genes with carotenoid accumulation, but also suggested a role for carotenoid degradation in L. japonica's dynamic flower coloration. The variation of flower color is consistent with not only the observed carotenoid accumulation pattern, but also with the release of volatile apocarotenoids that presumably serve as pollinator attractants. Beyond novel insights into the evolution and dynamics of flower coloration, the high-quality L. japonica genome sequence also provides a foundation for molecular breeding to improve desired characteristics.