A comprehensive evaluation of the potential of three next-generation short-read-based plant pan-genome construction strategies for the identification of novel non-reference sequence.

Pan-genome studies are important for understanding plant evolution and guiding the breeding of crops by containing all genomic diversity of a certain species. Three short-read-based strategies for plant pan-genome construction include iterative individual, iteration pooling, and map-to-pan. Their performance is very different under various conditions, while comprehensive evaluations have yet to be conducted nowadays. Here, we evaluate the performance of these three pan-genome construction strategies for plants under different sequencing depths and sample sizes. Also, we indicate the influence of length and repeat content percentage of novel sequences on three pan-genome construction strategies. Besides, we compare the computational resource consumption among the three strategies. Our findings indicate that map-to-pan has the greatest recall but the lowest precision. In contrast, both two iterative strategies have superior precision but lower recall. Factors of sample numbers, novel sequence length, and the percentage of novel sequences' repeat content adversely affect the performance of all three strategies. Increased sequencing depth improves map-to-pan's performance, while not affecting the other two iterative strategies. For computational resource consumption, map-to-pan demands considerably more than the other two iterative strategies. Overall, the iterative strategy, especially the iterative pooling strategy, is optimal when the sequencing depth is less than 20X. Map-to-pan is preferable when the sequencing depth exceeds 20X despite its higher computational resource consumption.

Aluminum ion chemistry of Na4Fe3(PO4)2(P2O7) for all-climate full Na-ion battery.

Na4Fe3(PO4)2(P2O7) (NFPP) is currently drawing increased attention as a sodium-ion batteries (SIBs) cathode due to the cost-effective and NASICON-type structure features. Owing to the sluggish electron and Na+ conductivities, however, its real implementation is impeded by the grievous capacity decay and inferior rate capability. Herein, multivalent cation substituted microporous Na3.9Fe2.9Al0.1(PO4)2(P2O7) (NFAPP) with wide operation-temperature is elaborately designed through regulating structure/interface coupled electron/ion transport. Greatly, the derived Na vacancy and charge rearrangement induced by trivalent Al3+ substitution lower the ions diffusion barriers, thereby endowing faster electron transport and Na+ mobility. More importantly, the existing Al-O-P bonds strengthen the local environment and alleviate the volume vibration during (de)sodiation, enabling highly reversible valence variation and structural evolution. As a result, remarkable cyclability (over 10,000 loops), ultrafast rate capability (200 C), and exceptional all-climate stability (-40-60 °C) in half/full cells are demonstrated. Given this, the rational work might provide an actionable strategy to promote the electrochemical property of NFPP, thus unveiling the great application prospect of sodium iron mixed phosphate materials.

CCLHunter: An efficient toolkit for cancer cell line authentication.

Cancer cell lines are essential in cancer research, yet accurate authentication of these cell lines can be challenging, particularly for consanguineous cell lines with close genetic similarities. We introduce a new Cancer Cell Line Hunter (CCLHunter) method to tackle this challenge. This approach utilizes the information of single nucleotide polymorphisms, expression profiles, and kindred topology to authenticate 1389 human cancer cell lines accurately. CCLHunter can precisely and efficiently authenticate cell lines from consanguineous lineages and those derived from other tissues of the same individual. Our evaluation results indicate that CCLHunter has a complete accuracy rate of 93.27%, with an accuracy of 89.28% even for consanguineous cell lines, outperforming existing methods. Additionally, we provide convenient access to CCLHunter through standalone software and a web server at https://ngdc.cncb.ac.cn/cclhunter.

Transcriptome-wide association studies: recent advances in methods, applications and available databases.

Genome-wide association study has identified fruitful variants impacting heritable traits. Nevertheless, identifying critical genes underlying those significant variants has been a great task. Transcriptome-wide association study (TWAS) is an instrumental post-analysis to detect significant gene-trait associations focusing on modeling transcription-level regulations, which has made numerous progresses in recent years. Leveraging from expression quantitative loci (eQTL) regulation information, TWAS has advantages in detecting functioning genes regulated by disease-associated variants, thus providing insight into mechanisms of diseases and other phenotypes. Considering its vast potential, this review article comprehensively summarizes TWAS, including the methodology, applications and available resources.

Gastrodin ameliorates cognitive dysfunction in diabetes by inhibiting PAK2 phosphorylation.

Diabetes is associated with higher prevalence of cognitive dysfunction, while the underlying mechanism is still elusive. In this study, we aim to explore the potential mechanism of diabetes-induced cognitive dysfunction and assess the therapeutic effects of Gastrodin on cognitive dysfunction. Diabetes was induced by a single injection of streptozotocin. The Morris Water Maze Test was employed to assess the functions of spatial learning and memory. Transcriptome was used to identify the potential factors involved. Western blot and immunofluorescence were applied to detect the protein expression. Our results have shown that spatial learning was impaired in diabetic rats, coupled with damaged hippocampal pyramidal neurons. Gastrodin intervention ameliorated the spatial learning impairments and neuronal damages. Transcriptomics analysis identified differential expression genes critical for diabetes-induced hippocampal damage and Gastrodin treatment, which were further confirmed by qPCR and western blot. Moreover, p21 activated kinase 2 (PAK2) was found to be important for diabetes-induced hippocampal injury and its inhibitor could promote the survival of primary hippocampal neurons. It suggested that PAK2 pathway may be involved in cognitive dysfunction in diabetes and could be a therapeutic target for Gastrodin intervention.

Elevated ZNF704 expression is associated with poor prognosis of uveal melanoma and promotes cancer cell growth by regulating AKT/mTOR signaling.

BACKGROUND: Uveal melanoma (UM) is the most common intraocular malignancy in adults, with a poor survival prognosis. To date, limited understanding of UM's molecular mechanisms constitutes an obstacle to developing effective therapy. In this study, we examined key regulators mediating UM progression and their clinical relevance. METHODS: Transcriptomics of UM patients and cells were analyzed via RNA sequencing and bioinformatic analysis. Zinc finger protein 704 (ZNF704) was identified as prognosis-related biomarker for UM based on clinical characteristics and RNA-seq data from The Cancer Genome Atlas (TCGA). Gene expression was knocked down by specific shRNAs/siRNAs and overexpressed by transfection with plasmids inserted with investigated gene cDNA. Cell proliferation, viability and invasion abilities were determined by CCK8, colony formation and transwell assays, respectively. For cell cycle and apoptosis, cells were PI or PI/Annexin V-APC stained and analyzed by flow cytometry. Standard immunoblotting and quantitative RT-PCR were employed to assess the mRNA and protein abundance. To determine tumor growth in vivo, 4-week-old BALB/c-nu immune-deficient nude mice were inoculated with tumor cells. RESULTS: Analysis of differential expressed genes (DEGs) and survival analysis identified ZNF704 as a novel biomarker of UM. Prognostic analysis indicated ZNF704 as an independent predictor of UM overall survival. Expression of ZNF704 is elevated in UM tissues relative to adjacent normal choroid tissues. Knockdown of ZNF704 suppressed the growth and migration of UM cells and vice versa. In addition, expression of ZNF704 arrest UM cells at G0/G1 phase and inhibit cell apoptosis. RNA sequencing analysis indicated that SORBS3 were dysregulated after ZNF704 downregulation. Gene Set Enrichment Analysis (GSEA) revealed that upon ZNF704 knowndown, genes related with PI3K/AKT/mTOR, EMT and metastasis are enriched. Mechanistically, ZNF704 activates AKT/mTOR/glycolysis signaling pathway in UM cells. Moreover, expression of SORBS3 is downregulated by ZNF704 and knockdown of SORBS3 restored tumor cell viability in ZNF704 silenced cells. CONCLUSIONS: ZNF704 predicts poor prognosis of UM and exhibit pro-oncogenic effect in UM progression in vivo and in vitro, mediated through AKT/mTOR signaling pathway and suppression of SORBS3 expression.

TWAS Atlas: a curated knowledgebase of transcriptome-wide association studies.

Transcriptome-wide association studies (TWASs), as a practical and prevalent approach for detecting the associations between genetically regulated genes and traits, are now leading to a better understanding of the complex mechanisms of genetic variants in regulating various diseases and traits. Despite the ever-increasing TWAS outputs, there is still a lack of databases curating massive public TWAS information and knowledge. To fill this gap, here we present TWAS Atlas (https://ngdc.cncb.ac.cn/twas/), an integrated knowledgebase of TWAS findings manually curated from extensive literature. In the current implementation, TWAS Atlas collects 401,266 high-quality human gene-trait associations from 200 publications, covering 22,247 genes and 257 traits across 135 tissue types. In particular, an interactive knowledge graph of the collected gene-trait associations is constructed together with single nucleotide polymorphism (SNP)-gene associations to build up comprehensive regulatory networks at multi-omics levels. In addition, TWAS Atlas, as a user-friendly web interface, efficiently enables users to browse, search and download all association information, relevant research metadata and annotation information of interest. Taken together, TWAS Atlas is of great value for promoting the utility and availability of TWAS results in explaining the complex genetic basis as well as providing new insights for human health and disease research.

TSomVar: a tumor-only somatic and germline variant identification method with random forest.

Somatic variants act as critical players during cancer occurrence and development. Thus, an accurate and robust method to identify them is the foundation of cutting-edge cancer genome research. However, due to low accessibility and high individual-/sample-specificity of the somatic variants in tumor samples, the detection is, to date, still crammed with challenges, particularly when lacking paired normal samples as control. To solve this burning issue, we developed a tumor-only somatic and germline variant identification method (TSomVar) using the random forest algorithm established on sample-specific variant datasets derived from genotype imputation, reads-mapping level annotation and functional annotation. We trained TSomVar by using genomic variant datasets of three major cancer types: colorectal cancer, hepatocellular carcinoma and skin cutaneous melanoma. Compared with existing tumor-only somatic variant identification tools, TSomVar shows excellent performances in somatic variant detection with higher accuracy and better capability of recalling for test datasets from colorectal cancer and skin cutaneous melanoma. In addition, TSomVar is equipped with the competence of accurately identifying germline variants in tumor samples. Taken together, TSomVar will undoubtedly facilitate and revolutionize somatic variant explorations in cancer research.

CancerSCEM: a database of single-cell expression map across various human cancers.

With the proliferating studies of human cancers by single-cell RNA sequencing technique (scRNA-seq), cellular heterogeneity, immune landscape and pathogenesis within diverse cancers have been uncovered successively. The exponential explosion of massive cancer scRNA-seq datasets in the past decade are calling for a burning demand to be integrated and processed for essential investigations in tumor microenvironment of various cancer types. To fill this gap, we developed a database of Cancer Single-cell Expression Map (CancerSCEM, https://ngdc.cncb.ac.cn/cancerscem), particularly focusing on a variety of human cancers. To date, CancerSCE version 1.0 consists of 208 cancer samples across 28 studies and 20 human cancer types. A series of uniformly and multiscale analyses for each sample were performed, including accurate cell type annotation, functional gene expressions, cell interaction network, survival analysis and etc. Plus, we visualized CancerSCEM as a user-friendly web interface for users to browse, search, online analyze and download all the metadata as well as analytical results. More importantly and unprecedentedly, the newly-constructed comprehensive online analyzing platform in CancerSCEM integrates seven analyze functions, where investigators can interactively perform cancer scRNA-seq analyses. In all, CancerSCEM paves an informative and practical way to facilitate human cancer studies, and also provides insights into clinical therapy assessments.

Enhancing Surface Chemical Stability of LiMn2 O4 Cathode by Strontium Enrichment at Grain Boundaries.

LiMn2 O4 (LMO) cathodes suffer from limited cycle life, resulting from Mn dissolution and side reactions between electrode and electrolyte. In this study, Sr-modified LMO is prepared by using a simple strategy. The nature and position of large-radius Sr ions are investigated, alongside their influence on the structural stability of the bulk. SrMnO3 (SMO) is found to be enriched at grain boundaries of LMO, with Mn-O-Sr bonds forming at the SMO/LMO interface. Furthermore, stable SMO alleviates the migration of Mn ions in LMO associated with structural integrity and suppresses side reactions between the electrode and electrolyte. The modified LMO cathodes maintain their structural integrity and display improved rate performance and cycling stability under harsh conditions. Remarkably, the discharge capacity of a Sr-modified LMO||Li half-cell maintains 94.8 % at 25 °C and 79.6 % at 55 °C after 500 cycles. Consequently, enrichment of strontium at grain boundaries presents a promising strategy for developing cathodes for long-term use.

RefRGim: an intelligent reference panel reconstruction method for genotype imputation with convolutional neural networks.

Genotype imputation is a statistical method for estimating missing genotypes from a denser haplotype reference panel. Existing methods usually performed well on common variants, but they may not be ideal for low-frequency and rare variants. Previous studies showed that the population similarity between study and reference panels is one of the key factors influencing the imputation accuracy. Here, we developed an imputation reference panel reconstruction method (RefRGim) using convolutional neural networks (CNNs), which can generate a study-specified reference panel for each input data based on the genetic similarity of individuals from current study and references. The CNNs were pretrained with single nucleotide polymorphism data from the 1000 Genomes Project. Our evaluations showed that genotype imputation with RefRGim can achieve higher accuracies than original reference panel, especially for low-frequency and rare variants. RefRGim will serve as an efficient reference panel reconstruction method for genotype imputation. RefRGim is freely available via GitHub: https://github.com/shishuo16/RefRGim.

Functionalized ionic liquids based on vegetable oils for rare earth elements recovery.

Functionalized ionic liquids (FILs) based on vegetable oils have been directly synthesized and used for the first time to extract rare earth elements (REEs). Carbon dioxide gas is introduced to successfully strip REE-loaded FILs in the presence of water. The novel extraction process reveals some advantages of accessibility, biocompatibility and sustainability as well as cost efficiency.

The Global Landscape of SARS-CoV-2 Genomes, Variants, and Haplotypes in 2019nCoVR.

On January 22, 2020, China National Center for Bioinformation (CNCB) released the 2019 Novel Coronavirus Resource (2019nCoVR), an open-access information resource for the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). 2019nCoVR features a comprehensive integration of sequence and clinical information for all publicly available SARS-CoV-2 isolates, which are manually curated with value-added annotations and quality evaluated by an automated in-house pipeline. Of particular note, 2019nCoVR offers systematic analyses to generate a dynamic landscape of SARS-CoV-2 genomic variations at a global scale. It provides all identified variants and their detailed statistics for each virus isolate, and congregates the quality score, functional annotation, and population frequency for each variant. Spatiotemporal change for each variant can be visualized and historical viral haplotype network maps for the course of the outbreak are also generated based on all complete and high-quality genomes available. Moreover, 2019nCoVR provides a full collection of SARS-CoV-2 relevant literature on the coronavirus disease 2019 (COVID-19), including published papers from PubMed as well as preprints from services such as bioRxiv and medRxiv through Europe PMC. Furthermore, by linking with relevant databases in CNCB, 2019nCoVR offers data submission services for raw sequence reads and assembled genomes, and data sharing with NCBI. Collectively, SARS-CoV-2 is updated daily to collect the latest information on genome sequences, variants, haplotypes, and literature for a timely reflection, making 2019nCoVR a valuable resource for the global research community. 2019nCoVR is accessible at https://bigd.big.ac.cn/ncov/.

RETRACTED: CGVD: a genomic variation database for Chinese populations.

Precision medicine calls upon deeper coverage of population-based sequencing and thorough gene-content and phenotype-based analysis, which lead to a population-associated genomic variation map or database. The Chinese Genomic Variation Database (CGVD; https://bigd.big.ac.cn/cgvd/) is such a database that has combined 48.30 million (M) SNVs and 5.77 M small indels, identified from 991 Chinese individuals of the Chinese Academy of Sciences Precision Medicine Initiative Project (CASPMI) and 301 Chinese individuals of the 1000 Genomes Project (1KGP). The CASPMI project includes whole-genome sequencing data (WGS, 25-30×) from ∼1000 healthy individuals of the CASPMI cohort. To facilitate the usage of such variations for pharmacogenomics studies, star-allele frequencies of the drug-related genes in the CASPMI and 1KGP populations are calculated and provided in CGVD. As one of the important database resources in BIG Data Center, CGVD will continue to collect more genomic variations and to curate structural and functional annotations to support population-based healthcare projects and studies in China and worldwide.

Whole Genome Analyses of Chinese Population and De Novo Assembly of A Northern Han Genome.

To unravel the genetic mechanisms of disease and physiological traits, it requires comprehensive sequencing analysis of large sample size in Chinese populations. Here, we report the primary results of the Chinese Academy of Sciences Precision Medicine Initiative (CASPMI) project launched by the Chinese Academy of Sciences, including the de novo assembly of a northern Han reference genome (NH1.0) and whole genome analyses of 597 healthy people coming from most areas in China. Given the two existing reference genomes for Han Chinese (YH and HX1) were both from the south, we constructed NH1.0, a new reference genome from a northern individual, by combining the sequencing strategies of PacBio, 10× Genomics, and Bionano mapping. Using this integrated approach, we obtained an N50 scaffold size of 46.63 Mb for the NH1.0 genome and performed a comparative genome analysis of NH1.0 with YH and HX1. In order to generate a genomic variation map of Chinese populations, we performed the whole-genome sequencing of 597 participants and identified 24.85 million (M) single nucleotide variants (SNVs), 3.85 M small indels, and 106,382 structural variations. In the association analysis with collected phenotypes, we found that the T allele of rs1549293 in KAT8 significantly correlated with the waist circumference in northern Han males. Moreover, significant genetic diversity in MTHFR, TCN2, FADS1, and FADS2, which associate with circulating folate, vitamin B12, or lipid metabolism, was observed between northerners and southerners. Especially, for the homocysteine-increasing allele of rs1801133 (MTHFR 677T), we hypothesize that there exists a "comfort" zone for a high frequency of 677T between latitudes of 35-45 degree North. Taken together, our results provide a high-quality northern Han reference genome and novel population-specific data sets of genetic variants for use in the personalized and precision medicine.

Complete Sequence and Analysis of Coconut Palm (Cocos nucifera) Mitochondrial Genome.

Coconut (Cocos nucifera L.), a member of the palm family (Arecaceae), is one of the most economically important crops in tropics, serving as an important source of food, drink, fuel, medicine, and construction material. Here we report an assembly of the coconut (C. nucifera, Oman local Tall cultivar) mitochondrial (mt) genome based on next-generation sequencing data. This genome, 678,653bp in length and 45.5% in GC content, encodes 72 proteins, 9 pseudogenes, 23 tRNAs, and 3 ribosomal RNAs. Within the assembly, we find that the chloroplast (cp) derived regions account for 5.07% of the total assembly length, including 13 proteins, 2 pseudogenes, and 11 tRNAs. The mt genome has a relatively large fraction of repeat content (17.26%), including both forward (tandem) and inverted (palindromic) repeats. Sequence variation analysis shows that the Ti/Tv ratio of the mt genome is lower as compared to that of the nuclear genome and neutral expectation. By combining public RNA-Seq data for coconut, we identify 734 RNA editing sites supported by at least two datasets. In summary, our data provides the second complete mt genome sequence in the family Arecaceae, essential for further investigations on mitochondrial biology of seed plants.

Identification and analysis of mouse non-coding RNA using transcriptome data.

Transcripts are expressed spatially and temporally and they are very complicated, precise and specific; however, most studies are focused on protein-coding related genes. Recently, massively parallel cDNA sequencing (RNA-seq) has emerged to be a new and promising tool for transcriptome research, and numbers of non-coding RNAs, especially lincRNAs, have been widely identified and well characterized as important regulators of diverse biological processes. In this study, we used ultra-deep RNA-seq data from 15 mouse tissues to study the diversity and dynamic of non-coding RNAs in mouse. Using our own criteria, we identified totally 16,249 non-coding genes (21,569 non-coding RNAs) in mouse. We annotated these non-coding RNAs by diverse properties and found non-coding RNAs are generally shorter, have fewer exons, express in lower level and are more strikingly tissue-specific compared with protein-coding genes. Moreover, these non-coding RNAs show significant enrichment with transcriptional initiation and elongation signals including histone modifications (H3K4me3, H3K27me3 and H3K36me3), RNAPII binding sites and CAGE tags. The gene set enrichment analysis (GSEA) result revealed several sets of lincRNAs associated with diverse biological processes such as immune effector process, muscle development and sexual reproduction. Taken together, this study provides a more comprehensive annotation of mouse non-coding RNAs and gives an opportunity for future functional and evolutionary study of mouse non-coding RNAs.

Identification and analysis of house-keeping and tissue-specific genes based on RNA-seq data sets across 15 mouse tissues.

Recently, RNA-seq has become widely used technology for transcriptome profiling due to its single-base accuracy and high-throughput speciality. In this study, we applied a computational approach on an integrated RNA-seq dataset across 15 normal mouse tissues, and consequently assigned 8408 house-keeping (HK) genes and 2581 tissue-specific (TS) genes among UCSC RefGene annotation. Apart from some basic genomic features, we also performed expression, function and pathway analysis with clustering, DAVID and Ingenuity Pathway Analysis, indicating the physiological connections (tissues) and diverse biological roles of HK genes (fundamental processes) and TS genes (tissue-corresponding processes). Moreover, we used RT-PCR method to test 18 candidate HK genes and finally identified a novel list of highly stable internal control genes: Ywhae, Ddb 1, Eif4h, etc. In summary, this study provides a new HK gene and TS gene resource for further genetic and evolution research and helps us better understand morphogenesis and biological diversity in mouse.