NextPolish2: A Repeat-aware Polishing Tool for Genomes Assembled Using HiFi Long Reads.

The high-fidelity (HiFi) long-read sequencing technology developed by PacBio has greatly improved the base-level accuracy of genome assemblies. However, these assemblies still contain base-level errors, particularly within the error-prone regions of HiFi long reads. Existing genome polishing tools usually introduce overcorrections and haplotype switch errors when correcting errors in genomes assembled from HiFi long reads. Here, we describe an upgraded genome polishing tool - NextPolish2, which can fix base errors remaining in those "highly accurate" genomes assembled from HiFi long reads without introducing excessive overcorrections and haplotype switch errors. We believe that NextPolish2 has a great significance to further improve the accuracy of telomere-to-telomere (T2T) genomes. NextPolish2 is freely available at https://github.com/Nextomics/NextPolish2.

Plant-on-chip: Core morphogenesis processes in the tiny plant Wolffia australiana.

A plant can be thought of as a colony comprising numerous growth buds, each developing to its own rhythm. Such lack of synchrony impedes efforts to describe core principles of plant morphogenesis, dissect the underlying mechanisms, and identify regulators. Here, we use the minimalist known angiosperm to overcome this challenge and provide a model system for plant morphogenesis. We present a detailed morphological description of the monocot Wolffia australiana, as well as high-quality genome information. Further, we developed the plant-on-chip culture system and demonstrate the application of advanced technologies such as single-nucleus RNA-sequencing, protein structure prediction, and gene editing. We provide proof-of-concept examples that illustrate how W. australiana can decipher the core regulatory mechanisms of plant morphogenesis.

Efficient assembly of nanopore reads via highly accurate and intact error correction.

Long nanopore reads are advantageous in de novo genome assembly. However, nanopore reads usually have broad error distribution and high-error-rate subsequences. Existing error correction tools cannot correct nanopore reads efficiently and effectively. Most methods trim high-error-rate subsequences during error correction, which reduces both the length of the reads and contiguity of the final assembly. Here, we develop an error correction, and de novo assembly tool designed to overcome complex errors in nanopore reads. We propose an adaptive read selection and two-step progressive method to quickly correct nanopore reads to high accuracy. We introduce a two-stage assembler to utilize the full length of nanopore reads. Our tool achieves superior performance in both error correction and de novo assembling nanopore reads. It requires only 8122 hours to assemble a 35X coverage human genome and achieves a 2.47-fold improvement in NG50. Furthermore, our assembly of the human WERI cell line shows an NG50 of 22 Mbp. The high-quality assembly of nanopore reads can significantly reduce false positives in structure variation detection.

Genome and single-cell RNA-sequencing of the earthworm Eisenia andrei identifies cellular mechanisms underlying regeneration.

The earthworm is particularly fascinating to biologists because of its strong regenerative capacity. However, many aspects of its regeneration in nature remain elusive. Here we report chromosome-level genome, large-scale transcriptome and single-cell RNA-sequencing data during earthworm (Eisenia andrei) regeneration. We observe expansion of LINE2 transposable elements and gene families functionally related to regeneration (for example, EGFR, epidermal growth factor receptor) particularly for genes exhibiting differential expression during earthworm regeneration. Temporal gene expression trajectories identify transcriptional regulatory factors that are potentially crucial for initiating cell proliferation and differentiation during regeneration. Furthermore, early growth response genes related to regeneration are transcriptionally activated in both the earthworm and planarian. Meanwhile, single-cell RNA-sequencing provides insight into the regenerative process at a cellular level and finds that the largest proportion of cells present during regeneration are stem cells.

DeepSignal: detecting DNA methylation state from Nanopore sequencing reads using deep-learning.

MOTIVATION: The Oxford Nanopore sequencing enables to directly detect methylation states of bases in DNA from reads without extra laboratory techniques. Novel computational methods are required to improve the accuracy and robustness of DNA methylation state prediction using Nanopore reads. RESULTS: In this study, we develop DeepSignal, a deep learning method to detect DNA methylation states from Nanopore sequencing reads. Testing on Nanopore reads of Homo sapiens (H. sapiens), Escherichia coli (E. coli) and pUC19 shows that DeepSignal can achieve higher performance at both read level and genome level on detecting 6 mA and 5mC methylation states comparing to previous hidden Markov model (HMM) based methods. DeepSignal achieves similar performance cross different DNA methylation bases, different DNA methylation motifs and both singleton and mixed DNA CpG. Moreover, DeepSignal requires much lower coverage than those required by HMM and statistics based methods. DeepSignal can achieve 90% above accuracy for detecting 5mC and 6 mA using only 2× coverage of reads. Furthermore, for DNA CpG methylation state prediction, DeepSignal achieves 90% correlation with bisulfite sequencing using just 20× coverage of reads, which is much better than HMM based methods. Especially, DeepSignal can predict methylation states of 5% more DNA CpGs that previously cannot be predicted by bisulfite sequencing. DeepSignal can be a robust and accurate method for detecting methylation states of DNA bases. AVAILABILITY AND IMPLEMENTATION: DeepSignal is publicly available at https://github.com/bioinfomaticsCSU/deepsignal. SUPPLEMENTARY INFORMATION: Supplementary data are available at bioinformatics online.

N6-Methyladenine DNA Modification in the Human Genome.

DNA N6-methyladenine (6mA) modification is the most prevalent DNA modification in prokaryotes, but whether it exists in human cells and whether it plays a role in human diseases remain enigmatic. Here, we showed that 6mA is extensively present in the human genome, and we cataloged 881,240 6mA sites accounting for ∼0.051% of the total adenines. [G/C]AGG[C/T] was the most significantly associated motif with 6mA modification. 6mA sites were enriched in the coding regions and mark actively transcribed genes in human cells. DNA 6mA and N6-demethyladenine modification in the human genome were mediated by methyltransferase N6AMT1 and demethylase ALKBH1, respectively. The abundance of 6mA was significantly lower in cancers, accompanied by decreased N6AMT1 and increased ALKBH1 levels, and downregulation of 6mA modification levels promoted tumorigenesis. Collectively, our results demonstrate that DNA 6mA modification is extensively present in human cells and the decrease of genomic DNA 6mA promotes human tumorigenesis.

[Speciation distribution and risk assessment of heavy metals in sediments in suburban outfall of industrial oasis region].

The speciation distribution and potential environmental risk of Cu, Zn, Cd, Pb, Ni, Cr in sediments in suburban outfall of industrial oasis region, Baiyin City were studied by Tessier sequential extraction method, correlation analysis, enrichment factor (EF) and risk assessment code (RAC). The results indicated that, the average concentration of heavy metals in sediments samples of Dongdagou and Xidagou streams exceeded the background levels in Gansu Province, with Cr, as an exception. The enrichment factor suggested that the enrichment of Cd in sediments of Dongdagou and Xidagou streams were very serious, which posed a strong pollution level. Furthermore, the analysis of chemical speciation indicated that Zn, Ni, and Cr in sediments of Dongdagou stream were mainly dominated by the fraction of residual, the existence of Cu and Cd was mainly in organic forms, while Pb was composed mostly by its Fe-Mn oxides fraction; Pb in sediments of Xidagou stream existed by Fe-Mn oxides fraction, other metals mainly appeared with the residual fraction. The risk assessment code (RAC) showed that the risk level of heavy metals in sediments of Dongdagou stream descended in the order: Ni > Cd > Pb > Cu > Zn > Cr, Ni posed a highest risk level. The risk level of heavy metals in sediments of Xidagou stream decreased in the order: Pb > Cd > Ni > Cu > Zn > Cr, Cd, Pb and Ni posed a highest risk level.

The protein level of hypoxia-inducible factor-1alpha is increased in the plateau pika (Ochotona curzoniae) inhabiting high altitudes.

The plateau pika (Ochotona curzoniae) is a high hypoxia-tolerant species living only at 3,000-5,000 m above sea-level on the Qinghai-Tibetan plateau. Hypoxia-inducible factor-1 (HIF-1) is a key transcription factor that regulates a variety of cellular and systemic adaptations to hypoxia. To investigate how the plateau pika adapts to a high-altitude hypoxic environment at the molecular level, we examined the expression pattern of the HIF-1alpha protein in the pika by Western blot and immunohistochemical analysis. We found that HIF-1alpha protein is expressed at a significantly high level in the pika, which is higher in most tissues (particularly in the lung, liver, spleen and kidney) of the plateau pika than that of mice living at sea-level. Importantly, we found that the protein levels of HIF-1alpha in the lung, liver, spleen and kidney of the pika were increased with increased habitat altitudes. We observed that the plateau pika HIF-1alpha localized to the nucleus of cells by an immunostaining analysis, and enhanced HRE-driven gene expression by luciferase reporter assays. Our study suggests that the HIF-1alpha protein levels are related to the adaptation of the plateau pika to the high-altitude hypoxic environment.

Natural selection and adaptive evolution of leptin in the ochotona family driven by the cold environmental stress.

BACKGROUND: Environmental stress can accelerate the evolutionary rate of specific stress-response proteins and create new functions specialized for different environments, enhancing an organism's fitness to stressful environments. Pikas (order Lagomorpha), endemic, non-hibernating mammals in the modern Holarctic Region, live in cold regions at either high altitudes or high latitudes and have a maximum distribution of species diversification confined to the Qinghai-Tibet Plateau. Variations in energy metabolism are remarkable for them living in cold environments. Leptin, an adipocyte-derived hormone, plays important roles in energy homeostasis. METHODOLOGY/PRINCIPAL FINDINGS: To examine the extent of leptin variations within the Ochotona family, we cloned the entire coding sequence of pika leptin from 6 species in two regions (Qinghai-Tibet Plateau and Inner Mongolia steppe in China) and the leptin sequences of plateau pikas (O. curzonia) from different altitudes on Qinghai-Tibet Plateau. We carried out both DNA and amino acid sequence analyses in molecular evolution and compared modeled spatial structures. Our results show that positive selection (PS) acts on pika leptin, while nine PS sites located within the functionally significant segment 85-119 of leptin and one unique motif appeared only in pika lineages-the ATP synthase alpha and beta subunit signature site. To reveal the environmental factors affecting sequence evolution of pika leptin, relative rate test was performed in pikas from different altitudes. Stepwise multiple regression shows that temperature is significantly and negatively correlated with the rates of non-synonymous substitution (Ka) and amino acid substitution (Aa), whereas altitude does not significantly affect synonymous substitution (Ks), Ka and Aa. CONCLUSIONS/SIGNIFICANCE: Our findings support the viewpoint that adaptive evolution may occur in pika leptin, which may play important roles in pikas' ecological adaptation to extreme environmental stress. We speculate that cold, and probably not hypoxia, may be the primary environmental factor for driving adaptive evolution of pika leptin.