The jujube (Ziziphus jujuba Mill.), a member of family Rhamnaceae, is a major dry fruit and a traditional herbal medicine for more than one billion people. Here we present a high-quality sequence for the complex jujube genome, the first genome sequence of Rhamnaceae, using an integrated strategy. The final assembly spans 437.65 Mb (98.6% of the estimated) with 321.45 Mb anchored to the 12 pseudo-chromosomes and contains 32,808 genes. The jujube genome has undergone frequent inter-chromosome fusions and segmental duplications, but no recent whole-genome duplication. Further analyses of the jujube-specific genes and transcriptome data from 15 tissues reveal the molecular mechanisms underlying some specific properties of the jujube. Its high vitamin C content can be attributed to a unique high level expression of genes involved in both biosynthesis and regeneration. Our study provides insights into jujube-specific biology and valuable genomic resources for the improvement of Rhamnaceae plants and other fruit trees.

Integrative analysis of multi-omics layers at single cell level is critical for accurate dissection of cell-to-cell variation within certain cell populations. Here we report scCAT-seq, a technique for simultaneously assaying chromatin accessibility and the transcriptome within the same single cell. We show that the combined single cell signatures enable accurate construction of regulatory relationships between cis-regulatory elements and the target genes at single-cell resolution, providing a new dimension of features that helps direct discovery of regulatory patterns specific to distinct cell identities. Moreover, we generate the first single cell integrated map of chromatin accessibility and transcriptome in early embryos and demonstrate the robustness of scCAT-seq in the precise dissection of master transcription factors in cells of distinct states. The ability to obtain these two layers of omics data will help provide more accurate definitions of "single cell state" and enable the deconvolution of regulatory heterogeneity from complex cell populations.

Polyploidization has provided much genetic variation for plant adaptive evolution, but the mechanisms by which the molecular evolution of polyploid genomes establishes genetic architecture underlying species differentiation are unclear. Brassica is an ideal model to increase knowledge of polyploid evolution. Here we describe a draft genome sequence of Brassica oleracea, comparing it with that of its sister species B. rapa to reveal numerous chromosome rearrangements and asymmetrical gene loss in duplicated genomic blocks, asymmetrical amplification of transposable elements, differential gene co-retention for specific pathways and variation in gene expression, including alternative splicing, among a large number of paralogous and orthologous genes. Genes related to the production of anticancer phytochemicals and morphological variations illustrate consequences of genome duplication and gene divergence, imparting biochemical and morphological variation to B. oleracea. This study provides insights into Brassica genome evolution and will underpin research into the many important crops in this genus.

Hepatitis B virus (HBV) can integrate into the human genome, contributing to genomic instability and hepatocarcinogenesis. Here by conducting high-throughput viral integration detection and RNA sequencing, we identify 4,225 HBV integration events in tumour and adjacent non-tumour samples from 426 patients with HCC. We show that HBV is prone to integrate into rare fragile sites and functional genomic regions including CpG islands. We observe a distinct pattern in the preferential sites of HBV integration between tumour and non-tumour tissues. HBV insertional sites are significantly enriched in the proximity of telomeres in tumours. Recurrent HBV target genes are identified with few that overlap. The overall HBV integration frequency is much higher in tumour genomes of males than in females, with a significant enrichment of integration into chromosome 17. Furthermore, a cirrhosis-dependent HBV integration pattern is observed, affecting distinct targeted genes. Our data suggest that HBV integration has a high potential to drive oncogenic transformation.

Prunus mume (mei), which was domesticated in China more than 3,000 years ago as ornamental plant and fruit, is one of the first genomes among Prunus subfamilies of Rosaceae been sequenced. Here, we assemble a 280M genome by combining 101-fold next-generation sequencing and optical mapping data. We further anchor 83.9% of scaffolds to eight chromosomes with genetic map constructed by restriction-site-associated DNA sequencing. Combining P. mume genome with available data, we succeed in reconstructing nine ancestral chromosomes of Rosaceae family, as well as depicting chromosome fusion, fission and duplication history in three major subfamilies. We sequence the transcriptome of various tissues and perform genome-wide analysis to reveal the characteristics of P. mume, including its regulation of early blooming in endodormancy, immune response against bacterial infection and biosynthesis of flower scent. The P. mume genome sequence adds to our understanding of Rosaceae evolution and provides important data for improvement of fruit trees.

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.

Pituitary tumors are frequently associated with mutations in the AIP gene and are sometimes associated with hypersecretion of growth hormone. It is unclear whether other factors besides an enlarged pituitary contribute to the hypersecretion. In a genetic screen for suppressors of reduced neurotransmitter release, we identified a mutation in Caenorhabditis elegans AIPR-1 (AIP-related-1), which causes profound increases in evoked and spontaneous neurotransmitter release, a high frequency of spontaneous calcium transients in motor neurons and an enlarged readily releasable pool of vesicles. Calcium bursts and hypersecretion are reversed by mutations in the ryanodine receptor but not in the voltage-gated calcium channel, indicating that these phenotypes are caused by a leaky ryanodine receptor. AIPR-1 is physically associated with the ryanodine receptor at synapses. Finally, the phenotypes in aipr-1 mutants can be rescued by presynaptic expression of mouse AIP, demonstrating that a conserved function of AIP proteins is to inhibit calcium release from ryanodine receptors.

The outbreak of coronavirus disease 2019 (COVID-19) is a global health emergency. Various omics results have been reported for COVID-19, but the molecular hallmarks of COVID-19, especially in those patients without comorbidities, have not been fully investigated. Here we collect blood samples from 231 COVID-19 patients, prefiltered to exclude those with selected comorbidities, yet with symptoms ranging from asymptomatic to critically ill. Using integrative analysis of genomic, transcriptomic, proteomic, metabolomic and lipidomic profiles, we report a trans-omics landscape for COVID-19. Our analyses find neutrophils heterogeneity between asymptomatic and critically ill patients. Meanwhile, neutrophils over-activation, arginine depletion and tryptophan metabolites accumulation correlate with T cell dysfunction in critical patients. Our multi-omics data and characterization of peripheral blood from COVID-19 patients may thus help provide clues regarding pathophysiology of and potential therapeutic strategies for COVID-19.

Mei (Prunus mume) is an ornamental woody plant that has been domesticated in East Asia for thousands of years. High diversity in floral traits, along with its recent genome sequence, makes mei an ideal model system for studying the evolution of woody plants. Here, we investigate the genetic architecture of floral traits in mei and its domestication history by sampling and resequencing a total of 351 samples including 348 mei accessions and three other Prunus species at an average sequencing depth of 19.3×. Highly-admixed population structure and introgression from Prunus species are identified in mei accessions. Through a genome-wide association study (GWAS), we identify significant quantitative traits locus (QTLs) and genomic regions where several genes, such as MYB108, are positively associated with petal color, stigma color, calyx color, and bud color. Results from this study shed light on the genetic basis of domestication in flowering plants, particularly woody plants.

Calorie restriction (CR) and fasting are common approaches to weight reduction, but the maintenance is difficult after resuming food consumption. Meanwhile, the gut microbiome associated with energy harvest alters dramatically in response to nutrient deprivation. Here, we reported that CR and high-fat diet (HFD) both remodeled the gut microbiota with similar microbial composition, Parabacteroides distasonis was most significantly decreased after CR or HFD. CR altered microbiota and reprogramed metabolism, resulting in a distinct serum bile acid profile characterized by depleting the proportion of non-12α-hydroxylated bile acids, ursodeoxycholic acid and lithocholic acid. Downregulation of UCP1 expression in brown adipose tissue and decreased serum GLP-1 were observed in the weight-rebound mice. Moreover, treatment with Parabacteroides distasonis or non-12α-hydroxylated bile acids ameliorated weight regain via increased thermogenesis. Our results highlighted the gut microbiota-bile acid crosstalk in rebound weight gain and Parabacteroides distasonis as a potential probiotic to prevent rapid post-CR weight gain.

Pigs are valuable large animal models for biomedical and genetic research, but insights into the tissue- and cell-type-specific transcriptome and heterogeneity remain limited. By leveraging single-cell RNA sequencing, we generate a multiple-organ single-cell transcriptomic map containing over 200,000 pig cells from 20 tissues/organs. We comprehensively characterize the heterogeneity of cells in tissues and identify 234 cell clusters, representing 58 major cell types. In-depth integrative analysis of endothelial cells reveals a high degree of heterogeneity. We identify several functionally distinct endothelial cell phenotypes, including an endothelial to mesenchymal transition subtype in adipose tissues. Intercellular communication analysis predicts tissue- and cell type-specific crosstalk between endothelial cells and other cell types through the VEGF, PDGF, TGF-β, and BMP pathways. Regulon analysis of single-cell transcriptome of microglia in pig and 12 other species further identifies MEF2C as an evolutionally conserved regulon in the microglia. Our work describes the landscape of single-cell transcriptomes within diverse pig organs and identifies the heterogeneity of endothelial cells and evolutionally conserved regulon in microglia.

Methods for sensitive and high-throughput evaluation of CRISPR RNA-guided nucleases (RGNs) off-targets (OTs) are essential for advancing RGN-based gene therapies. Here we report SURRO-seq for simultaneously evaluating thousands of therapeutic RGN OTs in cells. SURRO-seq captures RGN-induced indels in cells by pooled lentiviral OTs libraries and deep sequencing, an approach comparable and complementary to OTs detection by T7 endonuclease 1, GUIDE-seq, and CIRCLE-seq. Application of SURRO-seq to 8150 OTs from 110 therapeutic RGNs identifies significantly detectable indels in 783 OTs, of which 37 OTs are found in cancer genes and 23 OTs are further validated in five human cell lines by targeted amplicon sequencing. Finally, SURRO-seq reveals that thermodynamically stable wobble base pair (rG•dT) and free binding energy strongly affect RGN specificity. Our study emphasizes the necessity of thoroughly evaluating therapeutic RGN OTs to minimize inevitable off-target effects.

The grand challenge in the development of atomically dispersed metallic catalysts is their low metal-atom loading density, uncontrollable localization and ambiguous interactions with supports, posing difficulty in maximizing their catalytic performance. Here, we achieve an interface catalyst consisting of atomic cobalt array covalently bound to distorted 1T MoS2 nanosheets (SA Co-D 1T MoS2). The phase of MoS2 transforming from 2H to D-1T, induced by strain from lattice mismatch and formation of Co-S covalent bond between Co and MoS2 during the assembly, is found to be essential to form the highly active single-atom array catalyst. SA Co-D 1T MoS2 achieves Pt-like activity toward HER and high long-term stability. Active-site blocking experiment together with density functional theory (DFT) calculations reveal that the superior catalytic behaviour is associated with an ensemble effect via the synergy of Co adatom and S of the D-1T MoS2 support by tuning hydrogen binding mode at the interface.

Itaconate is a newly discovered endogenous metabolite promoting an anti-inflammatory program during innate immune response, but the precise mechanisms underlying its effect remains poorly understood owing primarily to the limitations of available itaconate-monitoring techniques. Here, we develop and validate a genetically encoded fluorescent itaconate biosensor, BioITA, for directly monitoring itaconate dynamics in subcellular compartments of living macrophages. Utilizing BioITA, we monitor the itaconate dynamics in response to lipopolysaccharide (LPS) stimulation in the context of modulating itaconate transportation and metabolism. Moreover, we show that STING activation induces itaconate production, and injection of AAVs expressing cytosolic BioITA into mice allows directly reporting elevation of itaconate level in activated macrophages derived from LPS-injected mice. Thus, BioITA enables subcellular resolution imaging of itaconate in living macrophages.

Proper threat-reward decision-making is critical to animal survival. Emerging evidence indicates that the motor system may participate in decision-making but the neural circuit and molecular bases for these functions are little known. We found in C. elegans that GABAergic motor neurons (D-MNs) bias toward the reward behavior in threat-reward decision-making by retrogradely inhibiting a pair of premotor command interneurons, AVA, that control cholinergic motor neurons in the avoidance neural circuit. This function of D-MNs is mediated by a specific ionotropic GABA receptor (UNC-49) in AVA, and depends on electrical coupling between the two AVA interneurons. Our results suggest that AVA are hub neurons where sensory inputs from threat and reward sensory modalities and motor information from D-MNs are integrated. This study demonstrates at single-neuron resolution how motor neurons may help shape threat-reward choice behaviors through interacting with other neurons.

Single cell approaches have increased our knowledge about the cell type composition of the non-human primate (NHP), but a detailed characterization of area-specific regulatory features remains outstanding. We generated single-cell transcriptomic and chromatin accessibility (single-cell ATAC) data of 358,237 cells from prefrontal cortex (PFC), primary motor cortex (M1) and primary visual cortex (V1) of adult female cynomolgus monkey brain, and integrated this dataset with Stereo-seq (spatial enhanced resolution omics-sequencing) of the corresponding cortical areas to assign topographic information to molecular states. We identified area-specific chromatin accessible sites and their targeted genes, including the cell type-specific transcriptional regulatory network associated with excitatory neurons heterogeneity. We reveal calcium ion transport and axon guidance genes related to specialized functions of PFC and M1, identified the similarities and differences between adult macaque and human oligodendrocyte trajectories, and mapped the genetic variants and gene perturbations of human diseases to NHP cortical cells. This resource establishes a transcriptomic and chromatin accessibility combinatory regulatory landscape at a single-cell and spatially resolved resolution in NHP cortex.

Culture-independent metagenomic studies have revolutionized our understanding of the gut microbiota. However, the lack of full genomes from cultured species is still a limitation for in-depth studies of the gut microbiota. Here we present a substantially expanded version of our Cultivated Genome Reference (CGR), termed CGR2, providing 3324 high-quality draft genomes from isolates selected from a large-scale cultivation of bacterial isolates from fecal samples of healthy Chinese individuals. The CGR2 classifies 527 species (179 previously unidentified species) from 8 phyla, and uncovers a genomic and functional diversity of Collinsella aerofaciens. The CGR2 genomes match 126 metagenome-assembled genomes without cultured representatives in the Unified Human Gastrointestinal Genome (UHGG) collection and harbor 3767 unidentified secondary metabolite biosynthetic gene clusters, providing a source of natural compounds with pharmaceutical potentials. We uncover accurate phage-bacterium linkages providing information on the evolutionary characteristics of interaction between bacteriophages and bacteria at the strain level.

Human cancer cell lines have long served as tools for cancer research and drug discovery, but the presence and the source of intra-cell-line heterogeneity remain elusive. Here, we perform single-cell RNA-sequencing and ATAC-sequencing on 42 and 39 human cell lines, respectively, to illustrate both transcriptomic and epigenetic heterogeneity within individual cell lines. Our data reveal that transcriptomic heterogeneity is frequently observed in cancer cell lines of different tissue origins, often driven by multiple common transcriptional programs. Copy number variation, as well as epigenetic variation and extrachromosomal DNA distribution all contribute to the detected intra-cell-line heterogeneity. Using hypoxia treatment as an example, we demonstrate that transcriptomic heterogeneity could be reshaped by environmental stress. Overall, our study performs single-cell multi-omics of commonly used human cancer cell lines and offers mechanistic insights into the intra-cell-line heterogeneity and its dynamics, which would serve as an important resource for future cancer cell line-based studies.

Using a whole-genome-sequencing approach to explore germplasm resources can serve as an important strategy for crop improvement, especially in investigating wild accessions that may contain useful genetic resources that have been lost during the domestication process. Here we sequence and assemble a draft genome of wild soybean and construct a recombinant inbred population for genotyping-by-sequencing and phenotypic analyses to identify multiple QTLs relevant to traits of interest in agriculture. We use a combination of de novo sequencing data from this work and our previous germplasm re-sequencing data to identify a novel ion transporter gene, GmCHX1, and relate its sequence alterations to salt tolerance. Rapid gain-of-function tests show the protective effects of GmCHX1 towards salt stress. This combination of whole-genome de novo sequencing, high-density-marker QTL mapping by re-sequencing and functional analyses can serve as an effective strategy to unveil novel genomic information in wild soybean to facilitate crop improvement.

Mudskippers are amphibious fishes that have developed morphological and physiological adaptations to match their unique lifestyles. Here we perform whole-genome sequencing of four representative mudskippers to elucidate the molecular mechanisms underlying these adaptations. We discover an expansion of innate immune system genes in the mudskippers that may provide defence against terrestrial pathogens. Several genes of the ammonia excretion pathway in the gills have experienced positive selection, suggesting their important roles in mudskippers' tolerance to environmental ammonia. Some vision-related genes are differentially lost or mutated, illustrating genomic changes associated with aerial vision. Transcriptomic analyses of mudskippers exposed to air highlight regulatory pathways that are up- or down-regulated in response to hypoxia. The present study provides a valuable resource for understanding the molecular mechanisms underlying water-to-land transition of vertebrates.