BGI Research, Sanya

The past two decades has witnessed a remarkable increase in the number of microbial genomes retrieved from marine systems1,2. However, it has remained challenging to translate this marine genomic diversity into biotechnological and biomedical applications3,4. Here we recovered 43,191 bacterial and archaeal genomes from publicly available marine metagenomes, encompassing a wide range of diversity with 138 distinct phyla, redefining the upper limit of marine bacterial genome size and revealing complex trade-offs between the occurrence of CRISPR–Cas systems and antibiotic resistance genes. In silico bioprospecting of these marine genomes led to the discovery of a novel CRISPR–Cas9 system, ten antimicrobial peptides, and three enzymes that degrade polyethylene terephthalate. In vitro experiments confirmed their effectiveness and efficacy. This work provides evidence that global-scale sequencing initiatives advance our understanding of how microbial diversity has evolved in the oceans and is maintained, and demonstrates how such initiatives can be sustainably exploited to advance biotechnology and biomedicine. Analysis of 43,191 genomes obtained from publicly available marine bacterial and archaeal metagenome data provides insights into marine bacterial evolution, CRISPR–Cas defence and antibiotic resistance genes, and demonstrates the potential of marine metagenomes for biotechnological applications.

SUMMARY: Large-scale comparative genomic studies have provided important insights into species evolution and diversity, but also lead to a great challenge to visualize. Quick catching or presenting key information hidden in the vast amount of genomic data and relationships among multiple genomes requires an efficient visualization tool. However, current tools for such visualization remain inflexible in layout and/or require advanced computation skills, especially for visualization of genome-based synteny. Here, we developed an easy-to-use and flexible layout tool, NGenomeSyn [multiple (N) Genome Synteny], for publication-ready visualization of syntenic relationships of the whole genome or local region and genomic features (e.g. repeats, structural variations, genes) across multiple genomes with a high customization. NGenomeSyn provides an easy way for its users to visualize a large amount of data with a rich layout by simply adjusting options for moving, scaling, and rotation of target genomes. Moreover, NGenomeSyn could be applied on the visualization of relationships on non-genomic data with similar input formats. AVAILABILITY AND IMPLEMENTATION: NGenomeSyn is freely available at GitHub (https://github.com/hewm2008/NGenomeSyn) and Zenodo (https://doi.org/10.5281/zenodo.7645148).

Quantifying spatiotemporal dynamics during embryogenesis is crucial for understanding congenital diseases. We developed Spateo (https://github.com/aristoteleo/spateo-release), a 3D spatiotemporal modeling framework, and applied it to a 3D mouse embryogenesis atlas at E9.5 and E11.5, capturing eight million cells. Spateo enables scalable, partial, non-rigid alignment, multi-slice refinement, and mesh correction to create molecular holograms of whole embryos. It introduces digitization methods to uncover multi-level biology from subcellular to whole organ, identifying expression gradients along orthogonal axes of emergent 3D structures, e.g., secondary organizers such as midbrain-hindbrain boundary (MHB). Spateo further jointly models intercellular and intracellular interaction to dissect signaling landscapes in 3D structures, including the zona limitans intrathalamica (ZLI). Lastly, Spateo introduces "morphometric vector fields" of cell migration and integrates spatial differential geometry to unveil molecular programs underlying asymmetrical murine heart organogenesis and others, bridging macroscopic changes with molecular dynamics. Thus, Spateo enables the study of organ ecology at a molecular level in 3D space over time.

Systematic exploration of the hadal zone, Earth's deepest oceanic realm, has historically faced technical limitations. Here, we collected 1,648 sediment samples at 6-11 km in the Mariana Trench, Yap Trench, and Philippine Basin for the Mariana Trench Environment and Ecology Research (MEER) project. Metagenomic and 16S rRNA gene amplicon sequencing generated the 92-Tbp MEER dataset, comprising 7,564 species (89.4% unreported), indicating high taxonomic novelty. Unlike in reported environments, neutral drift played a minimal role, while homogeneous selection (HoS, 50.5%) and dispersal limitation (DL, 43.8%) emerged as dominant ecological drivers. HoS favored streamlined genomes with key functions for hadal adaptation, e.g., aromatic compound utilization (oligotrophic adaptation) and antioxidation (high-pressure adaptation). Conversely, DL promoted versatile metabolism with larger genomes. These findings indicated that environmental factors drive the high taxonomic novelty in the hadal zone, advancing our understanding of the ecological mechanisms governing microbial ecosystems in such an extreme oceanic environment.

Global warming poses new challenges for insects to adapt to higher temperatures. Trehalose is the main blood sugar in insects and plays an important role in energy metabolism and stress resistance. The transmembrane transport of trehalose mainly depends on the trehalose transporter (TRET1). Plutella xylostella (L.) is a worldwide agricultural pest; however, the effects of the trehalose transport mechanism and trehalose distribution in tissues on the development, reproduction and temperature adaptation of P. xylostella have yet to be reported. In this study, PxTret1-like was cloned and analyzed regarding its expression pattern. It was found that the expression of PxTret1-like was affected by ambient temperature. The knockout mutation of PxTret1-like was generated using a CRISPR/Cas9 system by targeted knockout. The trehalose content and trehalase activity of mutant P. xylostella increased at different developmental stages. The trehalose content increased in the fat body of the fourth-instar P. xylostella, and decreased in the hemolymph, and there was no significant change in glucose in the fat body and hemolymph. Mutant strains of P. xylostella showed a significantly reduced survival rate, fecundity and ability to withstand extreme temperatures. The results showed that PxTret1-like could affect the development, reproduction and temperature adaptability of P. xylostella by regulating the trehalose content in the fat body and hemolymph.

induced enhanced neuronal activation in the hippocampal CA1 subregion. These findings indicate that gut microbial characteristics that are unique to patients with C-ADHD may be masked when considering a more heterogeneous group of patients. We link the gut microbiota to brain function in an ADHD animal model, suggesting the relevance of testing a potential bacteria-based intervention for some aspects of ADHD.

Obesity, a pervasive global public health challenge, impacts human well-being, and there is a growing interest in the potential of herbs for their anti-obesity effects. However, the specific herbs that possess these properties are not well understood. This study aims to identify and characterize herbs that promote weight loss and to elucidate the mechanisms involved. We employed a high-fat diet (HFD)-induced obesity mouse model to assess the efficacy of various herbs in weight reduction and to explore how these herbs interact with the gut microbiota. By leveraging the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP), we identified active ingredients present in the herbs and analyzed their regulatory effects on the gut microbiota. Among the 65 herbs, inulin, mori follum, corn silk, maydis seigma, crataegus, taraxacum, portulacae herba, rhizoma gastrodiae, rose rugosae flos, thalluslaminariae were found to significantly reduce body weight and restore the gut microbiota dysbiosis associated with HFD-induced obesity. We found that <i>Blautia</i>, which was enriched by the anti-obesity herbs, was inversely related to triglyceride (TG) levels. Conversely, <i>Desulfovibrio</i>, reduced by the anti-obesity herbs, was a positively linked to weight and weight gain. Furthermore, we discovered that arachidonic acid, shared by three anti-obesity herbs, was associated with beneficial bacteria such as <i>Blautia</i>. Our study identified 10 herbs with anti-obesity properties, detailing their active constituents and their influence on the gut microbiota. Our finding not only identified nearly 10 herbs that significantly reduce body weight, but also contribute to improving the understanding of the regulatory network between herbs and the gut microbiota.

Principal component analysis (PCA) is an important and widely used unsupervised learning method that determines population structure based on genetic variation. Genome sequencing of thousands of individuals usually generate tens of millions of SNPs, making it challenging for PCA analysis and interpretation. Here we present VCF2PCACluster, a simple, fast and memory-efficient tool for Kinship estimation, PCA and clustering analysis, and visualization based on VCF formatted SNPs. We implemented five Kinship estimation methods and three clustering methods for its users to choose from. Moreover, unlike other PCA tools, VCF2PCACluster possesses a clustering function based on PCA result, which enabling users to automatically and clearly know about population structure. We demonstrated the same accuracy but a higher performance of this tool in performing PCA analysis on tens of millions of SNPs compared to another popular PLINK2 software, especially in peak memory usage that is independent of the number of SNPs in VCF2PCACluster.

Abstract Fulminant myocarditis (FM) is a severe inflammatory condition of the myocardium that often results in sudden death, particularly in young individuals. In this study, we employed single-nucleus and spatial transcriptomics to perform a comprehensive analysis of coxsackievirus B3 (CVB3)-induced FM in A/J mice, spanning seven distinct time points pre- and post-treatment. Our findings reveal that mesothelial cells play a critical role in the early stage of myocarditis by acting as primary targets for CVB3 infection. This triggers the activation of macrophages, initiating a cascade of inflammation. Subsequently, pro-inflammatory Inflammatory_Mac and T cells infiltrate the myocardium, driving tissue damage. We also identified Cd8 + effector T cells as key mediators of cardiomyocyte injury. These cells release cytotoxic molecules, particularly IFN-γ, which modulates the expression of Spi1 , a factor implicated in exacerbating cardiomyocyte death and amplifying disease progression. Therapeutic interventions targeting the IFN-γ/ Spi1 axis demonstrated significant efficacy in FM models. Notably, intravenous immunoglobulin (IVIG) treatment reduced mortality, suppressed viral proliferation, and mitigated the hyperinflammatory state of FM. IVIG therapy also downregulated IFN-γ and Spi1 expression, underscoring its immunomodulatory and therapeutic potential. This comprehensive spatiotemporal transcriptomic analysis provides profound insights into the pathogenesis of FM and highlights actionable therapeutic targets, paving the way for more effective management strategies for this life-threatening condition.

Adult zebrafish robustly regenerate injured hearts through a complex orchestration of molecular and cellular activities. However, this remarkable process, which is largely non-existent in humans, remains incompletely understood. Here, we utilize integrated spatial transcriptomics (Stereo-seq) and single-cell RNA-sequencing (scRNA-seq) to generate a spatially-resolved molecular and cellular atlas of regenerating zebrafish heart across eight stages. We characterize the cascade of cardiomyocyte cell states responsible for producing regenerated myocardium and explore a potential role for tpm4a in cardiomyocyte re-differentiation. Moreover, we uncover the activation of ifrd1 and atp6ap2 genes as a unique feature of regenerative hearts. Lastly, we reconstruct a 4D “virtual regenerating heart” comprising 569,896 cells/spots derived from 36 scRNA-seq libraries and 224 Stereo-seq slices. Our comprehensive atlas serves as a valuable resource to the cardiovascular and regeneration scientific communities and their ongoing efforts to understand the molecular and cellular mechanisms underlying vertebrate heart regeneration. Adult zebrafish possess the ability to regenerate injured hearts through a complex orchestration of molecular and cellular activities. Here, the authors present a single cell and spatially-resolved atlas of the regenerating zebrafish heart, which serves as a valuable resource for research on heart regeneration.

The gut microbiota influences host immunity and metabolism, and changes in its composition and function have been implicated in several non-communicable diseases. Here, comparing germ-free (GF) and specific pathogen-free (SPF) mice using spatial transcriptomics, single-cell RNA sequencing, and targeted bile acid metabolomics across multiple organs, we systematically assessed how the gut microbiota's absence affected organ morphology, immune homeostasis, bile acid, and lipid metabolism. Through integrated analysis, we detect marked aberration in B, myeloid, and T/natural killer cells, altered mucosal zonation and nutrient uptake, and significant shifts in bile acid profiles in feces, liver, and circulation, with the alternate synthesis pathway predominant in GF mice and pronounced changes in bile acid enterohepatic circulation. Particularly, autophagy-driven lipid droplet breakdown in ileum epithelium and the liver's zinc finger and BTB domain-containing protein (ZBTB20)-Lipoprotein lipase (LPL) (ZBTB20-LPL) axis are key to plasma lipid homeostasis in GF mice. Our results unveil the complexity of microbiota-host interactions in the crosstalk between commensal gut bacteria and the host.

Spatial-transcriptomics technologies have demonstrated exceptional performance in characterizing brain and visceral organ tissues, as well as brain and retinal organoids. However, it has not yet been proven whether spatial transcriptomics can effectively characterize primary tissue-derived organoids, as the standardized tissue sectioning or slicing methods are not applicable for such organoids. Herein, we present a technique, lamination-based organoid spatially resolved transcriptomics (LOSRT), for organoid-spatially resolved transcriptomics based on organoid lamination. Primary mouse lung and liver-derived organoids were used in this study. The organoids were formulated using the droplet-engineering method and laminated using a homemade device with weight compression. This technique preserved most cells in individual organoids while maintaining delicate epithelium structures in laminated domains that can be recognized through visual segmentation. The mouse lung and liver organoids were resolved comprising various cell types, including alveolar cells, damage-associated transient progenitor cells, basal cells, macrophages, endothelial cells, fibroblasts, hepatocytes, and hepatic stellate cells. The distribution and count of cells were confirmed using immunohistology and identified with spatial transcriptomic features. This study reports an automated and integrated spatial transcriptomics method for primary organoids. It has the potential to standardize and rapidly characterize primary tissue-derived organoids.

Most grouper species are functional protogynous hermaphrodites, but the genetic basis and the molecular mechanisms underlying the regulation of this unique reproductive strategy remain enigmatic. In this study, we report a high-quality chromosome-level genome assembly of the representative orange-spotted grouper (Epinephelus coioides). No duplication or deletion of sex differentiation-related genes was found in the genome, suggesting that sex development in this grouper may be related to changes in regulatory sequences or environmental factors. Transcriptomic analyses showed that aromatase and retinoic acid are probably critical to promoting ovarian fate determination, and follicle-stimulating hormone triggers the female-to-male sex change. Socially controlled sex-change studies revealed that, in sex-changing fish, the brain's response to the social environment may be mediated by activation of the phototransduction cascade and the melatonin synthesis pathway. In summary, our genomic and experimental results provide novel insights into the molecular mechanisms of sex differentiation and sex change in the protogynous groupers.

Marine environments contain diverse halogenated organic compounds (HOCs), both anthropogenic and natural, nourishing a group of versatile organohalide-respiring bacteria (OHRB). Here, we identified a novel OHRB (Peptococcaceae DCH) with conserved motifs but phylogenetically diverse reductive dehalogenase catalytic subunit (RdhAs) from marine enrichment culture. Further analyses clearly demonstrate the horizontal gene transfer of rdhA s among marine OHRB. Moreover, 2,4,6-trichlorophenol (TCP) was dechlorinated to 2,4-dichlorophenol and terminated at 4-chlorophenol in culture. Dendrosporobacter and Methanosarcina were the two dominant genera, and the constructed and verified metabolic pathways clearly demonstrated that the former provided various substrates for other microbes, while the latter drew nutrients, but might provide little benefit to microbial dehalogenation. Furthermore, Dendrosporobacter could readily adapt to TCP, and sporulation-related proteins of Dendrosporobacter were significantly upregulated in TCP-free controls, whereas other microbes ( e.g., Methanosarcina and Aminivibrio ) became more active, providing insights into how HOCs shape microbial communities. Additionally, sulfate could affect the dechlorination of Peptococcaceae DCH, but not debromination. Considering their electron accessibility and energy generation, the results clearly demonstrate that bromophenols are more suitable than chlorophenols for the enrichment of OHRB in marine environments. This study will greatly enhance our understanding of marine OHRB ( rdhA s), auxiliary microbes, and microbial HOC adaptive mechanisms.

Breast cancers present intricate microenvironments comprising heterotypic cellular interactions, yet a comprehensive spatial map remained to be established. Here, we employed the DNA nanoball-based genome-wide in situ sequencing (Stereo-seq) to visualize the geospatial architecture of 30 primary breast tumors and metastatic lymph nodes across different molecular subtypes. This unprecedented high-resolution atlas unveils the fine structure of the tumor vasculature, highlighting heterogeneity in phenotype, spatial distribution, and intercellular communication within both endothelial and perivascular cells. In particular, venular smooth muscle cells are identified as the primary source of CCL21/CCL19 within the microenvironment. In collaboration with ACKR1-positive endothelial cells, they create a chemokine-rich venular niche to synergistically promote lymphocyte extravasation into tumors. High venule density predicts increased immune infiltration and improved clinical outcomes. This study provides a detailed spatial landscape of human breast cancer, offering key insights into the venular regulation of tumor immune infiltration. This study utilizes Stereo-seq to comprehensively map the geospatial architecture of breast cancer, revealing venular niches that promote immune infiltration and better clinical outcomes, offering new insights into tumor microenvironment regulation.

Summary Climate change‐induced drought is a major threat to agriculture. C 4 crops have a higher water use efficiency (WUE) and better adaptability to drought than C 3 crops due to their smaller stomatal morphology and faster response. However, our understanding of stomatal behaviours in both C 3 and C 4 Poaceae crops is limited by knowledge gaps in physical traits of guard cell (GC) and subsidiary cell (SC). We employed infrared gas exchange analysis and a stomatal assay to explore the relationship between GC/SC sizes and stomatal kinetics across diverse drought conditions in two C 3 (wheat and barley) and three C 4 (maize, sorghum and foxtail millet) upland Poaceae crops. Through statistical analyses, we proposed a GCSC‐ τ model to demonstrate how morphological differences affect stomatal kinetics in C 4 Poaceae crops. Our findings reveal that morphological variations specifically correlate with stomatal kinetics in C 4 Poaceae crops, but not in C 3 ones. Subsequent modelling and experimental validation provide further evidence that GC/SC sizes significantly impact stomatal kinetics, which affects stomatal responses to different drought conditions and thereby WUE in C 4 Poaceae crops. These findings emphasize the crucial advantage of GC/SC morphological characteristics and stomatal kinetics for the drought adaptability of C 4 Poaceae crops, highlighting their potential as future climate‐resilient crops.

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Plant evolution is driven by key innovations of functional traits that enables their survivals in diverse ecological environments. However, plant adaptive evolution from land to atmospheric niches remains poorly understood. In this study, we use the epiphytic Tillandsioideae subfamily of Bromeliaceae as model plants to explore their origin, evolution and diversification. We provide a comprehensive phylogenetic tree based on nuclear transcriptomic sequences, indicating that core tillandsioids originated approximately 11.3 million years ago in the Andes. The geological uplift of the Andes drives the divergence of tillandsioids into tank-forming and atmospheric types. Our genomic and transcriptomic analyses reveal gene variations and losses associated with adaptive traits such as impounding tanks and absorptive trichomes. Furthermore, we uncover specific nitrogen-fixing bacterial communities in the phyllosphere of tillandsioids as potential source of nitrogen acquisition. Collectively, our study provides integrative multi-omics insights into the adaptive evolution of tillandsioids in response to elevated aerial habitats. The mechanisms by which tillandsioids adapt to elevated aerial habitats remain largely unexplored. Here, the authors report their evolution and link life history, diversification, comparative genomic, and functional changes to processes underlying the unique biology of air plants.

DNA-segment copy number variations (DSCNVs), such as deletions and duplications, are important sources of genomic structural variation. However, the types and sizes of DSCNVs, as well as their genome-wide distribution and potential functions, are poorly understood in wheat. Here, we identified 198 985 DSCNVs by investigating 491 genomes of common wheat and found that they account for 20% of the entire genome. Interestingly, approximately 38% of genes are linked to DSCNVs. The number of DSCNVs within each accession ranges from 47 366 to 96 342, and their total sizes vary from 421.3 to 1267.9 Mb. We found that 957 and 1304 DSCNVs have been favored by breeders in China and the United States, respectively. By conducting DSCNV-based genome-wide association studies for the principal components of plant developmental and yield-component traits, we identified 34 loci as directly or indirectly involved in controlling the formation of multiple traits. Notably, a newly discovered DSCNV covering TaFT-D1 is significantly associated with flowering time and other agronomic traits. Overall, our findings highlight the potential of DSCNVs to drive fundamental discoveries in plant science. The comprehensive DSCNV map and the DSCNV-associated genes will also facilitate future research efforts to improve wheat yield, quality, and adaptation.

authors share their thoughts on three questions: how do you use AI tools ethically during manuscript writing? What benefits and risks do you believe AI tools will bring to scientific writing? Do you have any recommendations for better policies regulating AI tools' use in scientific writing?