Nanopore sequencing has been widely used for the reconstruction of microbial genomes. Owing to higher error rates, errors on the genome are corrected via neural networks trained by Nanopore reads. However, the systematic errors usually remain uncorrected. This paper designs a model that is trained by homologous sequences for the correction of Nanopore systematic errors. The developed program, Homopolish, outperforms Medaka and HELEN in bacteria, viruses, fungi, and metagenomic datasets. When combined with Medaka/HELEN, the genome quality can exceed Q50 on R9.4 flow cells. We show that Nanopore-only sequencing can produce high-quality microbial genomes sufficient for downstream analysis.

Key innovations have facilitated novel niche utilization, such as the movement of the algal predecessors of land plants into terrestrial habitats where drastic fluctuations in light intensity, ultraviolet radiation and water limitation required a number of adaptations. The NDH (NADH dehydrogenase-like) complex of Viridiplantae plastids participates in adapting the photosynthetic response to environmental stress, suggesting its involvement in the transition to terrestrial habitats. Although relatively rare, the loss or pseudogenization of plastid NDH genes is widely distributed across diverse lineages of photoautotrophic seed plants and mutants/transgenics lacking NDH function demonstrate little difference from wild type under non-stressed conditions. This study analyzes large transcriptomic and genomic datasets to evaluate the persistence and loss of NDH expression across plants. Nuclear expression profiles showed accretion of the NDH gene complement at key transitions in land plant evolution, such as the transition to land and at the base of the angiosperm lineage. While detection of transcripts for a selection of non-NDH, photosynthesis related proteins was independent of the state of NDH, coordinate, lineage-specific loss of plastid NDH genes and expression of nuclear-encoded NDH subunits was documented in Pinaceae, gnetophytes, Orchidaceae and Geraniales confirming the independent and complete loss of NDH in these diverse seed plant taxa. The broad phylogenetic distribution of NDH loss and the subtle phenotypes of mutants suggest that the NDH complex is of limited biological significance in contemporary plants. While NDH activity appears dispensable under favorable conditions, there were likely sufficiently frequent episodes of abiotic stress affecting terrestrial habitats to allow the retention of NDH activity. These findings reveal genetic factors influencing plant/environment interactions in a changing climate through 450 million years of land plant evolution.

Long-read phasing has been used for reconstructing diploid genomes, improving variant calling and resolving microbial strains in metagenomics. However, the phasing blocks of existing methods are broken by large Structural Variations (SVs), and the efficiency is unsatisfactory for population-scale phasing. This article presents a novel algorithm, LongPhase, which can simultaneously phase single nucleotide polymorphisms (SNPs) and SVs of a human genome in 10–20 min, 10× faster than the state-of-the-art WhatsHap, HapCUT2 and Margin. In particular, co-phasing SNPs and SVs produces much larger haplotype blocks (N50 = 25 Mbp) than those of existing methods (N50 = 10–15 Mbp). We show that LongPhase combined with Nanopore ultra-long reads is a cost-effective and highly contiguous solution, which can produce between one and 26 blocks per chromosome arm without the need for additional trios, chromosome-conformation and strand-seq data. LongPhase is freely available at https://github.com/twolinin/LongPhase/. Supplementary data are available at Bioinformatics online.

The bamboo Bambusa edulis has a long juvenile phase in situ, but can be induced to flower during in vitro tissue culture, providing a readily available source of material for studies on reproductive biology and flowering. In this report, in vitro-derived reproductive and vegetative materials of B. edulis were harvested and used to generate transcriptome databases by use of two sequencing platforms: Illumina and 454. Combination of the two datasets resulted in high transcriptome quality and increased length of the sequence reads. In plants, many MADS genes control flower development, and the ABCDE model has been developed to explain how the genes function together to create the different whorls within a flower.

The accuracy of Oxford Nanopore Technology (ONT) sequencing has significantly improved thanks to new flowcells, sequencing kits, and basecalling algorithms. However, novel modification types untrained in the basecalling models can seriously reduce the quality. Here we reports a set of ONT-sequenced genomes with unexpected low quality due to novel modification types. Demodification by whole-genome amplification significantly improved the quality but lost the epigenome. We also developed a reference-based method, Modpolish, for correcting modification-mediated errors while retaining the epigenome when a sufficient number of closely-related genomes is publicly available (default: top 20 genomes with at least 95% identity). Modpolish not only significantly improved the quality of in-house sequenced genomes but also public datasets sequenced by R9.4 and R10.4 (simplex). Our results suggested that novel modifications are prone to ONT systematic errors. Nevertheless, these errors are correctable by nucleotide demodification or Modpolish without prior knowledge of modifications.

Motivation: Integrated analysis of single nucleotide polymorphisms (SNPs) and structure variations showed that the extent of linkage disequilibrium is common across different types of genetic variants. A subset of SNPs (called tag SNPs) is sufficient for capturing alleles of bi-allelic and even multi-allelic variants. However, accuracy and power of tag SNPs are affected by several factors, including genotyping failure, errors and tagging bias of certain alleles. In addition, different sets of tag SNPs should be selected for fulfilling requirements of various genotyping platforms and projects. Results: This study formulates the problem of selecting tag SNPs into a four-objective optimization problem that minimizes the total amount of tag SNPs, maximizes tolerance for missing data, enlarges and balances detection power of each allele class. To resolve this problem, we propose evolutionary algorithms incorporated with greedy initialization to find non-dominated solutions considering all objectives simultaneously. This method provides users with great flexibility to extract different sets of tag SNPs for different platforms and scenarios (e.g. up to 100 tags and 10% missing rate). Compared to conventional methods, our method explores larger search space and requires shorter convergence time. Experimental results revealed strong and weak conflicts among these objectives. In particular, a small number of additional tag SNPs can provide sufficient tolerance and balanced power given the low missing and error rates of today's genotyping platforms. Availability: The software is freely available at Bioinformatics online and http://cilab.cs.ccu.edu.tw/service_dl.html Contact:ckting@cs.ccu.edu.tw; ythuang@cs.ccu.edu.tw