Quality control (QC) and preprocessing are essential steps for sequencing data analysis to ensure the accuracy of results. However, existing tools cannot provide a satisfying solution with integrated comprehensive functions, proper architectures, and highly scalable acceleration. In this article, we demonstrate SOAPnuke as a tool with abundant functions for a "QC-Preprocess-QC" workflow and MapReduce acceleration framework. Four modules with different preprocessing functions are designed for processing datasets from genomic, small RNA, Digital Gene Expression, and metagenomic experiments, respectively. As a workflow-like tool, SOAPnuke centralizes processing functions into 1 executable and predefines their order to avoid the necessity of reformatting different files when switching tools. Furthermore, the MapReduce framework enables large scalability to distribute all the processing works to an entire compute cluster.We conducted a benchmarking where SOAPnuke and other tools are used to preprocess a ∼30× NA12878 dataset published by GIAB. The standalone operation of SOAPnuke struck a balance between resource occupancy and performance. When accelerated on 16 working nodes with MapReduce, SOAPnuke achieved ∼5.7 times the fastest speed of other tools.

Antidiabetic medication may modulate the gut microbiota and thereby alter plasma and faecal bile acid (BA) composition, which may improve metabolic health. Here we show that treatment with Acarbose, but not Glipizide, increases the ratio between primary BAs and secondary BAs and plasma levels of unconjugated BAs in treatment-naive type 2 diabetes (T2D) patients, which may beneficially affect metabolism. Acarbose increases the relative abundances of Lactobacillus and Bifidobacterium in the gut microbiota and depletes Bacteroides, thereby changing the relative abundance of microbial genes involved in BA metabolism. Treatment outcomes of Acarbose are dependent on gut microbiota compositions prior to treatment. Compared to patients with a gut microbiota dominated by Prevotella, those with a high abundance of Bacteroides exhibit more changes in plasma BAs and greater improvement in metabolic parameters after Acarbose treatment. Our work highlights the potential for stratification of T2D patients based on their gut microbiota prior to treatment.

Hypopharyngeal cancer (HPC) frequently presents at an advanced stage, resulting in poor prognosis. Although combined surgical therapy and chemoradiotherapy have improved the survival for patients with HPC over the past 3 decades, the mortality rate in late-stage diagnosis of HPC is unsatisfactory. In this study, we performed whole-exome sequencing (WES) of 23 hypopharyngeal tumor and paired adjacent normal tissue to identify novel candidate driver genes associated with hypopharyngeal carcinoma. We identified several copy number variants (CNVs) and 15 somatic mutation genes that were associated with hypopharyngeal carcinoma. Mutations in nine new genes (PRB4, NSD1, REC8, ZNF772, ZNF69, EI24, CYFIP2, NEFH, KRTAP4-5) were also indentified. PRB4 and NSD1 expression were significantly upregulated in hypopharyngeal carcinoma, which was confirmed in an independent cohort using IHC. There was a positive relationship between PRB4 and NSD1. Downregulation of PRB4 by siRNA could inhibit cell growth, colony formation and cell invasion. Notably, we here demonstrate that NSD1 could bind to the promoter regions of PRB4 and activate promoter activity by reducing the binding of H3K27me2 and increasing the binding of H3K36me2 on PRB4 promoter. In summary, we pinpoint the predominant mutations in hypopharyngeal carcinoma by WES, highlighting the substantial genetic alterations contributing to hypopharyngeal carcinoma tumorigenesis. We also indentify a novel epigenetically regulatory between PRB4 and NSD1 that contribute to hypopharyngeal carcinoma tumorigenesis. They may become potential prognostic biomarkers and therapeutic target for hypopharyngeal carcinoma treatment.

Oesophageal carcinoma is the fourth leading cause of cancer-related death in China, and more than 90% of these tumours are oesophageal squamous cell carcinoma (ESCC). Although several ESCC genomic sequencing studies have identified mutated somatic genes, the number of samples in each study was relatively small, and the molecular basis of ESCC has not been fully elucidated. Here, we performed an integrated analysis of 490 tumours by combining the genomic data from 7 previous ESCC projects. We identified 18 significantly mutated genes (SMGs). PTEN, DCDC1 and CUL3 were first reported as SMGs in ESCC. Notably, the AJUBA mutations and mutational signature4 were significantly correlated with a poorer survival in patients with ESCC. Hierarchical clustering analysis of the copy number alteration (CNA) of cancer gene census (CGC) genes in ESCC patients revealed three subtypes, and subtype3 exhibited more CNAs and marked for worse prognosis compared with subtype2. Moreover, database annotation suggested that two significantly differential CNA genes (PIK3CA and FBXW7) between subtype3 and subtype2 may serve as therapeutic drug targets. This study has extended our knowledge of the genetic basis of ESCC and shed some light into the clinical relevance, which would help improve the therapy and prognosis of ESCC patients.

Rhododendron delavayi Franch. is globally famous as an ornamental plant. Its distribution in southwest China covers several different habitats and environments. However, not much research had been conducted on Rhododendron spp. at the molecular level, which hinders understanding of its evolution, speciation, and synthesis of secondary metabolites, as well as its wide adaptability to different environments. Here, we report the genome assembly and gene annotation of R. delavayi var. delavayi (the second genome sequenced in the Ericaceae), which will facilitate the study of the family. The genome assembly will have further applications in genome-assisted cultivar breeding. The final size of the assembled R. delavayi var. delavayi genome (695.09 Mb) was close to the 697.94 Mb, estimated by k-mer analysis. A total of 336.83 gigabases (Gb) of raw Illumina HiSeq 2000 reads were generated from 9 libraries (with insert sizes ranging from 170 bp to 40 kb), achieving a raw sequencing depth of ×482.6. After quality filtering, 246.06 Gb of clean reads were obtained, giving ×352.55 coverage depth. Assembly using Platanus gave a total scaffold length of 695.09 Mb, with a contig N50 of 61.8 kb and a scaffold N50 of 637.83 kb. Gene prediction resulted in the annotation of 32 938 protein-coding genes. The genome completeness was evaluated by CEGMA and BUSCO and reached 95.97% and 92.8%, respectively. The gene annotation completeness was also evaluated by CEGMA and BUSCO and reached 97.01% and 87.4%, respectively. Genome annotation revealed that 51.77% of the R. delavayi genome is composed of transposable elements, and 37.48% of long terminal repeat elements (LTRs). The de novo assembled genome of R. delavayi var. delavayi (hereinafter referred to as R. delavayi) is the second genomic resource of the family Ericaceae and will provide a valuable resource for research on future comparative genomic studies in Rhododendron species. The availability of the R. delavayi genome sequence will hopefully provide a tool for scientists to tackle open questions regarding molecular mechanisms underlying environmental interactions in the genus Rhododendron, more accurately understand the evolutionary processes and systematics of the genus, facilitate the identification of genes encoding pharmaceutically important compounds, and accelerate molecular breeding to release elite varieties.

Achieving a near complete understanding of how the genome of an individual affects the phenotypes of that individual requires deciphering the order of variations along homologous chromosomes in species with diploid genomes. However, true diploid assembly of long-range haplotypes remains challenging.

Cultivated rice varieties are all diploid, and polyploidization of rice has long been desired because of its advantages in genome buffering, vigorousness, and environmental robustness. However, a workable route remains elusive. Here, we describe a practical strategy, namely de novo domestication of wild allotetraploid rice. By screening allotetraploid wild rice inventory, we identified one genotype of Oryza alta (CCDD), polyploid rice 1 (PPR1), and established two important resources for its de novo domestication: (1) an efficient tissue culture, transformation, and genome editing system and (2) a high-quality genome assembly discriminated into two subgenomes of 12 chromosomes apiece. With these resources, we show that six agronomically important traits could be rapidly improved by editing O. alta homologs of the genes controlling these traits in diploid rice. Our results demonstrate the possibility that de novo domesticated allotetraploid rice can be developed into a new staple cereal to strengthen world food security.

Rich fossil evidence suggests that many traits and functions related to terrestrial evolution were present long before the ancestor of lobe- and ray-finned fishes. Here, we present genome sequences of the bichir, paddlefish, bowfin, and alligator gar, covering all major early divergent lineages of ray-finned fishes. Our analyses show that these species exhibit many mosaic genomic features of lobe- and ray-finned fishes. In particular, many regulatory elements for limb development are present in these fishes, supporting the hypothesis that the relevant ancestral regulation networks emerged before the origin of tetrapods. Transcriptome analyses confirm the homology between the lung and swim bladder and reveal the presence of functional lung-related genes in early ray-finned fishes. Furthermore, we functionally validate the essential role of a jawed vertebrate highly conserved element for cardiovascular development. Our results imply the ancestors of jawed vertebrates already had the potential gene networks for cardio-respiratory systems supporting air breathing.

Au nanoclusters, decorated with graphene quantum dots (GQDs), were obtained through photocatalytic reduction of AuCl43- by UV irradiation, and then cytarabine (Cyt) was loaded to the Au/GQDs via charge-dipole interactions. Mercaptopropionic acid (MPA) was anchored to the Cyt-loaded Au/GQDs through the formation of Au-S bond, which was further encapsulated by polyethyleneimine (PEI) via charge-dipole interactions. The delivery of Cyt from the quaternary complex (Au/GQDs/MPA/PEI) is pH-sensitive and can be modulated by near-infrared (NIR) irradiation. The results of cell viability test indicate that the developed nanoplatform can be used for chemo-photothermal combination therapy of cancer cells, and the efficacy of chemo-photothermal combination therapy is significantly higher than that of the single mode of photothermal therapy (PTT) or chemotherapy.

Anglerfishes are a highly diverse group of species with unique characteristics. Here, we report the first chromosome-level genome of a species in the order Lophiiformes, the yellow goosefish (Lophius litulon), obtained by whole genome shotgun sequencing and high-throughput chromatin conformation capture. Approximately 97.20% of the assembly spanning 709.23 Mb could be anchored to 23 chromosomes with a contig N50 of 164.91 kb. The BUSCO value was 95.4%, suggesting that the quality of the assembly was high. A comparative gene family analysis identified expanded and contracted gene families, and these may be associated with adaptation to the benthic environment and the lack of scales in the species. A majority of positively selected genes were related to metabolic processes, suggesting that digestive and metabolic system evolution expanded the diversity of yellow goosefish prey. Our study provides a valuable genetic resource for understanding the mechanisms underlying the unique features of the yellow goosefish and for investigating anglerfish evolution.

Lung metastasis is the major cause of breast cancer-related mortality. The neutrophil-associated inflammatory microenvironment aids tumor cells in metastatic colonization in lungs. Here, we show that tumor-secreted protease cathepsin C (CTSC) promotes breast-to-lung metastasis by regulating recruitment of neutrophils and formation of neutrophil extracellular traps (NETs). CTSC enzymatically activates neutrophil membrane-bound proteinase 3 (PR3) to facilitate interleukin-1β (IL-1β) processing and nuclear factor κB activation, thus upregulating IL-6 and CCL3 for neutrophil recruitment. In addition, the CTSC-PR3-IL-1β axis induces neutrophil reactive oxygen species production and formation of NETs, which degrade thrombospondin-1 and support metastatic growth of cancer cells in the lungs. CTSC expression and secretion are associated with NET formation and lung metastasis in human breast tumors. Importantly, targeting CTSC with compound AZD7986 effectively suppresses lung metastasis of breast cancer in a mouse model. Overall, our findings reveal a mechanism of how tumor cells regulate neutrophils in metastatic niches and support CTSC-targeting approaches for cancer treatment.

There is a dearth of information on the occurrence and risks of antibiotics in the urban rivers from plateau areas. This study investigated 83 antibiotics in water and sediments of an urban river and effluents of sewage treatment plants (E-STPs) in Xining, Qinghai (northeastern Tibetan Plateau). Fifty-three antibiotics were detected, and the concentrations of individual antibiotics varied in the range of undetected (ND)-552 ng/L in water, ND-164 ng/g in sediments, and ND-3821 ng/L in E-STPs. Seasonal differences of antibiotic concentrations were significant for water samples (p < 0.05) but insignificant for sediments (p > 0.05). In urban area, E-STP is the main source of antibiotics in the river, while runoff from manured cropland contributes partially to antibiotics in the river in the suburban area. The antibiotic compositions in water were different from those in sediments, but were similar to those in E-STPs. Notably, because of strong solar radiation and long sunshine hours in the plateau area, low levels of quinolones, which are sensitive to photolysis, were observed in river water. Moreover, norfloxacin and enrofloxacin, observed in urban river from other regions of China, were not detected in the Huangshui River water. The occurrence of ofloxacin, erythromycin, roxithromycin, clarithromycin, and trimethoprim in E-STPs may induce a possible risk to antibiotic resistance evolution. Trimethoprim, anhydroerythromycin, sulfamethoxazole, sulfapyridine, and clindamycin in river water could pose low to medium risks to aquatic organisms. Further investigation on the occurrence and distribution of antibiotic resistance genes in the Huangshui River is urgently needed.

Two of the most economically important plants in the Artocarpus genus are jackfruit (A. heterophyllus Lam.) and breadfruit (A. altilis (Parkinson) Fosberg). Both species are long-lived trees that have been cultivated for thousands of years in their native regions. Today they are grown throughout tropical to subtropical areas as an important source of starch and other valuable nutrients. There are hundreds of breadfruit varieties that are native to Oceania, of which the most commonly distributed types are seedless triploids. Jackfruit is likely native to the Western Ghats of India and produces one of the largest tree-borne fruit structures (reaching up to 45 kg). To-date, there is limited genomic information for these two economically important species. Here, we generated 273 Gb and 227 Gb of raw data from jackfruit and breadfruit, respectively. The high-quality reads from jackfruit were assembled into 162,440 scaffolds totaling 982 Mb with 35,858 genes. Similarly, the breadfruit reads were assembled into 180,971 scaffolds totaling 833 Mb with 34,010 genes. A total of 2822 and 2034 expanded gene families were found in jackfruit and breadfruit, respectively, enriched in pathways including starch and sucrose metabolism, photosynthesis, and others. The copy number of several starch synthesis-related genes were found to be increased in jackfruit and breadfruit compared to closely-related species, and the tissue-specific expression might imply their sugar-rich and starch-rich characteristics. Overall, the publication of high-quality genomes for jackfruit and breadfruit provides information about their specific composition and the underlying genes involved in sugar and starch metabolism.

In plants, parasitism triggers the reductive evolution of plastid genomes (plastomes). To disentangle the molecular evolutionary associations between feeding on other plants below- or aboveground and general transitions from facultative to obligate parasitism, we analyzed 34 complete plastomes of autotrophic, root- and stem-feeding hemiparasitic, and holoparasitic Santalales. We observed inexplicable losses of housekeeping genes and tRNAs in hemiparasites and dramatic genomic reconfiguration in holoparasitic Balanophoraceae, whose plastomes have exceptionally low GC contents. Genomic changes are related primarily to the evolution of hemi- or holoparasitism, whereas the transition from a root- to a stem-feeding mode plays no major role. In contrast, the rate of molecular evolution accelerates in a stepwise manner from autotrophs to root- and then stem-feeding parasites. Already the ancestral transition to root-parasitism coincides with a relaxation of selection in plastomes. Another significant selectional shift in plastid genes occurs as stem-feeders evolve, suggesting that this derived form coincides with trophic specialization despite the retention of photosynthetic capacity. Parasitic Santalales fill a gap in our understanding of parasitism-associated plastome degeneration. We reveal that lifestyle-genome associations unfold interdependently over trophic specialization and feeding mode transitions, where holoparasitic Balanophoraceae provide a system for exploring the functional realms of plastomes.

Mangroves are a group of plant species that occupy the coastal intertidal zone and are major components of this ecologically important ecosystem. Mangroves belong to about twenty diverse families. Here, we sequenced and assembled chloroplast genomes of 14 mangrove species from eight families spanning five rosid orders and one asterid order: Fabales (Pongamia pinnata), Lamiales (Avicennia marina), Malpighiales (Excoecaria agallocha, Bruguiera sexangula, Kandelia obovata, Rhizophora stylosa, and Ceriops tagal), Malvales (Hibiscus tiliaceus, Heritiera littoralis, and Thespesia populnea), Myrtales (Laguncularia racemosa, Sonneratia ovata, and Pemphis acidula), and Sapindales (Xylocarpus moluccensis). These chloroplast genomes range from 149 kb to 168 kb in length. A conserved structure of two inverted repeats (IRa and IRb, ~25.8 kb), one large single-copy region (LSC, ~89.0 kb), and one short single-copy region (SSC, ~18.9 kb) as well as ~130 genes (85 protein-coding, 37 tRNAs, and 8 rRNAs) was observed. We found the lowest divergence in the IR regions among the four regions. We also identified simple sequence repeats (SSRs), which were found to be variable in numbers. Most chloroplast genes are highly conserved, with only four genes under positive selection or relaxed pressure. Combined with publicly available chloroplast genomes, we carried out phylogenetic analysis and confirmed the previously reported phylogeny within rosids, including the positioning of obscure families in Malpighiales. Our study reports 14 mangrove chloroplast genomes and illustrates their genome features and evolution.

Rationale: KRAS is one of the most frequently mutated oncogenes in cancers. The protein's picomolar affinity for GTP/GDP and smooth protein structure resulting in the absence of known allosteric regulatory sites makes its genomic-level activating mutations a difficult but attractive target. Methods: Two CRISPR systems, genome-editing CRISPR/SpCas9 and transcription-regulating dCas9-KRAB, were developed to deplete the KRAS G12S mutant allele or repress its transcription, respectively, with the goal of treating KRAS-driven cancers. Results: SpCas9 and dCas9-KRAB systems with a sgRNA targeting the mutant allele blocked the expression of the mutant KRAS gene, leading to an inhibition of cancer cell proliferation. Local adenoviral injections using SpCas9 and dCas9-KRAB systems suppressed tumor growth in vivo. The gene-depletion system (SpCas9) performed more effectively than the transcription-suppressing system (dCas9-KRAB) on tumor inhibition. Application of both Cas9 systems to wild-type KRAS tumors did not affect cell proliferation. Furthermore, through bioinformatic analysis of 31555 SNP mutations of the top 20 cancer driver genes, the data showed that our mutant-specific editing strategy could be extended to a reference list of oncogenic mutations with high editing potentials. This pipeline could be applied to analyze the distribution of PAM sequences and survey the best alternative targets for gene editing. Conclusion: We successfully developed both gene-depletion and transcription-suppressing systems to specifically target an oncogenic KRAS mutant allele that led to significant tumor regression. These findings show the potential of CRISPR-based strategies for the treatment of tumors with driver gene mutations.

The diamondback moth, Plutella xylostella is a cosmopolitan pest that has evolved resistance to all classes of insecticide, and costs the world economy an estimated US $4-5 billion annually. We analyse patterns of variation among 532 P. xylostella genomes, representing a worldwide sample of 114 populations. We find evidence that suggests South America is the geographical area of origin of this species, challenging earlier hypotheses of an Old-World origin. Our analysis indicates that Plutella xylostella has experienced three major expansions across the world, mainly facilitated by European colonization and global trade. We identify genomic signatures of selection in genes related to metabolic and signaling pathways that could be evidence of environmental adaptation. This evolutionary history of P. xylostella provides insights into transoceanic movements that have enabled it to become a worldwide pest.

Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.

Carbon catabolite repression 4 (CCR4) is a conserved mRNA deadenylase regulating posttranscriptional gene expression. However, regulation of CCR4 in virus infections is less understood. Here, we characterized a pro-viral role of CCR4 in replication of a plant cytorhabdovirus, Barley yellow striate mosaic virus (BYSMV). The barley (Hordeum vulgare) CCR4 protein (HvCCR4) was identified to interact with the BYSMV phosphoprotein (P). The BYSMV P protein recruited HvCCR4 from processing bodies (PBs) into viroplasm-like bodies. Overexpression of HvCCR4 promoted BYSMV replication in plants. Conversely, knockdown of the small brown planthopper CCR4 inhibited viral accumulation in the insect vector. Biochemistry experiments revealed that HvCCR4 was recruited into N-RNA complexes by the BYSMV P protein and triggered turnover of N-bound cellular mRNAs, thereby releasing RNA-free N protein to bind viral genomic RNA for optimal viral replication. Our results demonstrate that the co-opted CCR4-mediated RNA decay facilitates cytorhabdovirus replication in plants and insects.

Epstein-Barr virus (EBV) is a well-recognized oncogenic virus that can induce host cell metabolic reprogramming and tumorigenesis by targeting vital metabolic enzymes or regulators. This study aims to explore the role of wild-type isocitrate dehydrogenase 2 (IDH2) in metabolic reprogramming and tumorigenesis induced by EBV-encoded latent membrane protein 1 (LMP1).