Single-cell RNA sequencing illuminates the ontogeny, conservation and diversification of cartilaginous and bony fish lymphocytes.

Elucidating cellular architecture and cell-type evolution across species is central to understanding immune system function and susceptibility to disease. Adaptive immunity is a shared trait of the common ancestor of cartilaginous and bony fishes. However, evolutionary features of lymphocytes in these two jawed vertebrates remain unclear. Here, we present a single-cell RNA sequencing atlas of immune cells from cartilaginous (white-spotted bamboo shark) and bony (zebrafish and Chinese tongue sole) fishes. Cross-species comparisons show that the same cell types across different species exhibit similar transcriptional profiles. In the bamboo shark, we identify a phagocytic B cell population expressing several pattern recognition receptors, as well as a T cell sub-cluster co-expressing both T and B cell markers. In contrast to a division by function in the bony fishes, we show close linkage and poor functional specialization among lymphocytes in the cartilaginous fish. Our cross-species single-cell comparison presents a resource for uncovering the origin and evolution of the gnathostome immune system.

Discovery of an unconventional lamprey lymphocyte lineage highlights divergent features in vertebrate adaptive immune system evolution.

Lymphocyte receptors independently evolved in both jawed and jawless vertebrates with similar adaptive immune responses. However, the diversity of functional subtypes and molecular architecture in jawless vertebrate lymphocytes, comparable to jawed species, is not well defined. Here, we profile the gills, intestines, and blood of the lamprey, Lampetra morii, with single-cell RNA sequencing, using a full-length transcriptome as a reference. Our findings reveal higher tissue-specific heterogeneity among T-like cells in contrast to B-like cells. Notably, we identify a unique T-like cell subtype expressing a homolog of the nonlymphoid hematopoietic growth factor receptor, MPL-like (MPL-L). These MPL-L+ T-like cells exhibit features distinct from T cells of jawed vertebrates, particularly in their elevated expression of hematopoietic genes. We further discovered that MPL-L+ VLRA+ T-like cells are widely present in the typhlosole, gill, liver, kidney, and skin of lamprey and they proliferate in response to both a T cell mitogen and recombinant human thrombopoietin. These findings provide new insights into the adaptive immune response in jawless vertebrates, shedding new light on the evolution of adaptive immunity.

Datasets of fungal diversity and pseudo-chromosomal genomes of mangrove rhizosphere soil in China.

With climate change and anthropic influence on the coastal ecosystems, mangrove ecosystems are disappearing at an alarming rate. Accordingly, it becomes important to track, study, record and store the mangrove microbial community considering their ecological importance and potential for biotechnological applications. Here, we provide information on mangrove fungal community composition and diversity in mangrove ecosystems with different plant species and from various locations differing in relation to anthropic influences. We describe twelve newly assembled genomes, including four chromosomal-level genomes of fungal isolates from the mangrove ecosystems coupled with functional annotations. We envisage that these data will be of value for future studies including comparative genome analysis and large-scale temporal and/or spatial research to elucidate the potential mechanisms by which mangrove fungal communities assemble and evolve. We further anticipate that the genomes represent valuable resources for bioprospecting related to industrial or clinical uses.

Global citrus root microbiota unravels assembly cues and core members.

Introduction: Citrus is one of the most important fruit crops worldwide, and the root-associated microbiota can have a profound impact on tree health and growth. Methods: In a collaborative effort, the International Citrus Microbiome Consortium investigated the global citrus root microbiota with samples collected from nine citrus-producing countries across six continents. We analyzed 16S rDNA and ITS2 amplicon sequencing data to identify predominant prokaryotic and fungal taxa in citrus root samples. Comparative analyses were conducted between root-associated microbial communities and those from the corresponding rhizosphere and bulk soil samples. Additionally, genotype-based group-wise comparisons were performed to assess the impact of citrus genotype on root microbiota composition. Results: Ten predominant prokaryotic phyla, containing nine bacterial phyla including Proteobacteria, Actinobacteria, Acidobacteria, and Bacteroidetes and one archaeal phylum (Thaumarchaeota), and multiple fungal phyla including Ascomycota and Basidiomycota were identified in the citrus root samples. Compared with the microbial communities from the corresponding rhizosphere and bulk soil samples from the same trees, the prokaryotic and fungal communities in the roots exhibited lower diversity and complexity but greater modularity compared to those in the rhizosphere. In total, 30 root-enriched and 150 root-depleted genera in bacterial community were identified, whereas 21 fungal genera were enriched, and 147 fungal genera were depleted in the root niche compared with the rhizosphere. The citrus genotype significantly affected the root prokaryotic and fungal communities. In addition, we have identified the core root prokaryotic genera comprising Acidibacter, Allorhizobium, Bradyrhizobium, Chitinophaga, Cupriavidus, Devosia, Dongia, Niastella, Pseudomonas, Sphingobium, Steroidobacter and Streptomyces, and the core fungal genera including Acrocalymma, Cladosporium, Fusarium, Gibberella, Mortierella, Neocosmospora and Volutella. The potential functions of these core genera of root microbiota were predicted. Conclusion: Overall, this study provides new insights into the assembly of microbial communities and identifies core members of citrus root microbiota across a wide geographic range. The findings offer valuable information for manipulating root microbiota to enhance plant growth and health.

Horizontal gene transfer and symbiotic microorganisms regulate the adaptive evolution of intertidal algae, Porphyra sense lato.

Intertidal algae may adapt to environmental challenges by acquiring genes from other organisms and relying on symbiotic microorganisms. Here, we obtained a symbiont-free and chromosome-level genome of Pyropia haitanensis (47.2 Mb), a type of intertidal algae, by using multiple symbiont screening methods. We identified 286 horizontal gene transfer (HGT) genes, 251 of which harbored transposable elements (TEs), reflecting the importance of TEs for facilitating the transfer of genes into P. haitanensis. Notably, the bulked segregant analysis revealed that two HGT genes, sirohydrochlorin ferrochelatase and peptide-methionine (R)-S-oxide reductase, play a significant role in the adaptation of P. haitanensis to heat stress. Besides, we found Pseudomonas, Actinobacteria, and Bacteroidetes are the major taxa among the symbiotic bacteria of P. haitanensis (nearly 50% of the HGT gene donors). Among of them, a heat-tolerant actinobacterial strain (Saccharothrix sp.) was isolated and revealed to be associated with the heat tolerance of P. haitanensis through its regulatory effects on the genes involved in proline synthesis (proC), redox homeostasis (ggt), and protein folding (HSP20). These findings contribute to our understanding of the adaptive evolution of intertidal algae, expanding our knowledge of the HGT genes and symbiotic microorganisms to enhance their resilience and survival in challenging intertidal environments.

The Chromosome-Scale Genome of Chitala ornata Illuminates the Evolution of Early Teleosts.

Teleosts are the most prolific vertebrates, occupying the vast majority of aquatic environments, and their pectoral fins have undergone remarkable physiological transformations throughout their evolution. Studying early teleost fishes, such as those belonging to the Osteoglossiformes order, could offer crucial insights into the adaptive evolution of pectoral fins within this group. In this study, we have assembled a chromosomal-level genome for the Clown featherback (Chitala ornata), achieving the highest quality genome assembly for Osteoglossiformes to date, with a contig N50 of 32.78 Mb and a scaffold N50 of 40.73 Mb. By combining phylogenetic analysis, we determined that the Clown featherback diverged approximately 202 to 203 million years ago (Ma), aligning with continental separation events. Our analysis revealed the intriguing discovery that a unique deletion of regulatory elements is adjacent to the Gli3 gene, specifically in teleosts. This deletion might be tied to the specialized adaptation of their pectoral fins. Furthermore, our findings indicate that specific contractions and expansions of transposable elements (TEs) in teleosts, including the Clown featherback, could be connected to their adaptive evolution. In essence, this study not only provides a high-quality genomic resource for Osteoglossiformes but also sheds light on the evolutionary trajectory of early teleosts.

Sirenian genomes illuminate the evolution of fully aquatic species within the mammalian superorder afrotheria.

Sirenians of the superorder Afrotheria were the first mammals to transition from land to water and are the only herbivorous marine mammals. Here, we generated a chromosome-level dugong (Dugong dugon) genome. A comparison of our assembly with other afrotherian genomes reveals possible molecular adaptations to aquatic life by sirenians, including a shift in daily activity patterns (circadian clock) and tolerance to a high-iodine plant diet mediated through changes in the iodide transporter NIS (SLC5A5) and its co-transporters. Functional in vitro assays confirm that sirenian amino acid substitutions alter the properties of the circadian clock protein PER2 and NIS. Sirenians show evidence of convergent regression of integumentary system (skin and its appendages) genes with cetaceans. Our analysis also uncovers gene losses that may be maladaptive in a modern environment, including a candidate gene (KCNK18) for sirenian cold stress syndrome likely lost during their evolutionary shift in daily activity patterns. Genomes from nine Australian locations and the functionally extinct Okinawan population confirm and date a genetic break ~10.7 thousand years ago on the Australian east coast and provide evidence of an associated ecotype, and highlight the need for whole-genome resequencing data from dugong populations worldwide for conservation and genetic management.

Exploring the cobia (Rachycentron canadum) genome: unveiling putative male heterogametic regions and identification of sex-specific markers.

BACKGROUND: Cobia (Rachycentron canadum) is the only member of the Rachycentridae family and exhibits considerable sexual dimorphism in growth rate. Sex determination in teleosts has been a long-standing basic biological question, and the molecular mechanisms of sex determination/differentiation in cobia are completely unknown. RESULTS: Here, we reported 2 high-quality, chromosome-level annotated male and female cobia genomes with assembly sizes of 586.51 Mb (contig/scaffold N50: 86.0 kb/24.3 Mb) and 583.88 Mb (79.9 kb/22.5 Mb), respectively. Synteny inference among perciform genomes revealed that cobia and the remora Echeneis naucrates were sister groups. Further, whole-genome resequencing of 31 males and 60 females, genome-wide association study, and sequencing depth analysis identified 3 short male-specific regions within a 10.7-kb continuous genomic region on male chromosome 18, which hinted at an undifferentiated sex chromosome system with a putative XX/XY mode of sex determination in cobia. Importantly, the only 2 genes within/between the male-specific regions, epoxide hydrolase 1 (ephx1, renamed cephx1y) and transcription factor 24 (tcf24, renamed ctcf24y), showed testis-specific/biased gene expression, whereas their counterparts cephx1x and ctf24x, located in female chromosome 18, were similarly expressed in both sexes. In addition, male-specific PCR targeting the cephx1y gene revealed that this genomic feature is conserved in cobia populations from Panama, Brazil, Australia, and Japan. CONCLUSION: The first comprehensive genomic survey presented here is a valuable resource for future studies on cobia population structure and dynamics, conservation, and evolutionary history. Furthermore, it establishes evidence of putative male heterogametic regions with 2 genes playing a potential role in the sex determination of the species, and it provides further support for the rapid evolution of sex-determining mechanisms in teleost fish.

Functional vertical connectivity of microbial communities in the ocean.

Sinking particles are a critical conduit for the transport of surface microbes to the ocean's interior. Vertical connectivity of phylogenetic composition has been shown; however, the functional vertical connectivity of microbial communities has not yet been explored in detail. We investigated protein and taxa profiles of both free-living and particle-attached microbial communities from the surface to 3000 m depth using a combined metaproteomic and 16S rRNA amplicon sequencing approach. A clear compositional and functional vertical connectivity of microbial communities was observed throughout the water column with Oceanospirillales, Alteromonadales, and Rhodobacterales as key taxa. The surface-derived particle-associated microbes increased the expression of proteins involved in basic metabolism, organic matter processing, and environmental stress response in deep waters. This study highlights the functional vertical connectivity between surface and deep-sea microbial communities via sinking particles and reveals that a considerable proportion of the deep-sea microbes might originate from surface waters and have a major impact on the biogeochemical cycles in the deep sea.

A new invertebrate NPY-like polypeptide, ZoaNPY, from the Zoanthus sociatus, as a novel ligand of human NPY Y2 receptor rescues vascular insufficiency via PLC/PKC and Src- FAK-dependent signaling pathways.

Our recent multi-omics studies have revealed rich sources of novel bioactive proteins and polypeptides from marine organisms including cnidarians. In the present study, we initially conducted a transcriptomic analysis to review the composition profile of polypeptides from Zoanthus sociatus. Then, a newly discovered NPY-like polypeptide-ZoaNPY was selected for further in silico structural, binding and virtually pharmacological studies. To evaluate the pro-angiogenic effects of ZoaNPY, we employed an in vitro HUVECs model and an in vivo zebrafish model. Our results indicate that ZoaNPY, at 1-100 pmol, enhances cell survival, migration and tube formation in the endothelial cells. Besides, treatment with ZoaNPY could restore a chemically-induced vascular insufficiency in zebrafish embryos. Western blot results demonstrated the application of ZoaNPY could increase the phosphorylation of proteins related to angiogenesis signaling including PKC, PLC, FAK, Src, Akt, mTOR, MEK, and ERK1/2. Furthermore, through molecular docking and surface plasmon resonance (SPR) verification, ZoaNPY was shown to directly and physically interact with NPY Y2 receptor. In view of this, all evidence showed that the pro-angiogenic effects of ZoaNPY involve the activation of NPY Y2 receptor, thereby activating the Akt/mTOR, PLC/PKC, ERK/MEK and Src- FAK-dependent signaling pathways. Furthermore, in an excision wound model, the treatment with ZoaNPY was shown to accelerate the wound healing process in mice. Our findings provide new insights into the discovery and development of novel pro-angiogenic drugs derived from NPY-like polypeptides in the future.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution.

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

Spatially resolved single-cell atlas of ascidian endostyle provides insight into the origin of vertebrate pharyngeal organs.

The pharyngeal endoderm, an innovation of deuterostome ancestors, contributes to pharyngeal development by influencing the patterning and differentiation of pharyngeal structures in vertebrates; however, the evolutionary origin of the pharyngeal organs in vertebrates is largely unknown. The endostyle, a distinct pharyngeal organ exclusively present in basal chordates, represents a good model for understanding pharyngeal organ origins. Using Stereo-seq and single-cell RNA sequencing, we constructed aspatially resolved single-cell atlas for the endostyle of the ascidian Styela clava. We determined the cell composition of the hemolymphoid region, which illuminates a mixed ancestral structure for the blood and lymphoid system. In addition, we discovered a cluster of hair cell-like cells in zone 3, which has transcriptomic similarity with the hair cells of the vertebrate acoustico-lateralis system. These findings reshape our understanding of the pharynx of the basal chordate and provide insights into the evolutionary origin of multiplexed pharyngeal organs.

In situ community transcriptomics illuminates CO2-fixation potentials and supporting roles of phagotrophy and proton pump in plankton in a subtropical marginal sea.

Lineage-wise physiological activities of plankton communities in the ocean are important but challenging to characterize. Here, we conducted whole-assemblage metatranscriptomic profiling at continental shelf and slope sites in the South China Sea to investigate carbon fixation potential in different lineages. RuBisCO expression, the proxy of Calvin carbon fixation (CCF) potential, was mainly contributed by Bacillariophyta, Chlorophyta, Cyanobacteria, and Haptophyta, which was differentially affected by environmental factors among lineages. CCF potential exhibited positive or negative correlations with phagotrophy gene expression, suggesting phagotrophy possibly enhances or complements CCF. Our data also reveal significant non-Calvin carbon fixation (NCF) potential, as indicated by the active expression of genes in all five currently recognized NCF pathways, mainly contributed by Flavobacteriales, Alteromonadales, and Oceanospirillales. Furthermore, in Flavobacteriales, Alteromonadales, Pelagibacterales, and Rhodobacterales, NCF potential was positively correlated with proton-pump rhodopsin (PPR) expression, suggesting that NCF might be energetically supported by PPR. The novel insights into the lineage-differential potential of carbon fixation, widespread mixotrophy, and PPR as an energy source for NCF lay a methodological and informational foundation for further research to understand carbon fixation and the trophic landscape in the ocean.IMPORTANCEMarine plankton plays an important role in global carbon cycling and climate regulation. Phytoplankton and cyanobacteria fix CO2 to produce organic compounds using solar energy and mainly by the Calvin cycle, whereas autotrophic bacteria and archaea may fix CO2 by non-Calvin cycle carbon fixation pathways. How active individual lineages are in carbon fixation and mixotrophy, and what energy source bacteria may employ in non-Calvin carbon fixation, in a natural plankton assemblage are poorly understood and underexplored. Using metatranscriptomics, we studied carbon fixation in marine plankton with lineage resolution in tropical marginal shelf and slope areas. Based on the sequencing results, we characterized the carbon fixation potential of different lineages and assessed Calvin- and non-Calvin- carbon fixation activities and energy sources. Data revealed a high number of unigenes (4.4 million), lineage-dependent differential potentials of Calvin carbon fixation and responses to environmental conditions, major contributors of non-Calvin carbon fixation, and their potential energy source.

Genomics of post-bottleneck recovery in the northern elephant seal.

Populations and species are threatened by human pressure, but their fate is variable. Some depleted populations, such as that of the northern elephant seal (Mirounga angustirostris), recover rapidly even when the surviving population was small. The northern elephant seal was hunted extensively and taken by collectors between the early 1800s and 1892, suffering an extreme population bottleneck as a consequence. Recovery was rapid and now there are over 200,000 individuals. We sequenced 260 modern and 8 historical northern elephant seal nuclear genomes to assess the impact of the population bottleneck on individual northern elephant seals and to better understand their recovery. Here we show that inbreeding, an increase in the frequency of alleles compromised by lost function, and allele frequency distortion, reduced the fitness of breeding males and females, as well as the performance of adult females on foraging migrations. We provide a detailed investigation of the impact of a severe bottleneck on fitness at the genomic level and report on the role of specific gene systems.

Single-cell transcriptome landscape of zebrafish liver reveals hepatocytes and immune cell interactions in understanding nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is becoming the most common chronic liver disease in the world. Immunity is the major contributing factor in NAFLD; however, the interaction of immune cells and hepatocytes in disease progression has not been fully elucidated. As a popular species for studying NAFLD, zebrafish, whose liver is a complex immune system mediated by immune cells and non-immune cells in maintaining immune tolerance and homeostasis. Understanding the cellular composition and immune environment of zebrafish liver is of great significance for its application in NAFLD. Here, we established a liver atlas that consists of 10 cell types using single-cell RNA sequencing (scRNA-seq). By examining the heterogeneity of hepatocytes and analyzing the expression of NAFLD-associated genes in the specific cluster, we provide a potential target cell model to study NAFLD. Additionally, our analysis identified two subtypes of distinct resident macrophages with inflammatory and non-inflammatory functions and characterized the successive stepwise development of T cell subclusters in the liver. Importantly, we uncovered the possible regulation of macrophages and T cells on target cells of fatty liver by analyzing the cellular interaction between hepatocytes and immune cells. Our data provide valuable information for an in-depth study of immune cells targeting hepatocytes to regulate the immune balance in NAFLD.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis).

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes.

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

A single-cell atlas of West African lungfish respiratory system reveals evolutionary adaptations to terrestrialization.

The six species of lungfish possess both lungs and gills and are the closest extant relatives of tetrapods. Here, we report a single-cell transcriptome atlas of the West African lungfish (Protopterus annectens). This species manifests the most extreme form of terrestrialization, a life history strategy to survive dry periods that can last for years, characterized by dormancy and reversible adaptive changes of the gills and lungs. Our atlas highlights the cell type diversity of the West African lungfish, including gene expression consistent with phenotype changes of terrestrialization. Comparison with terrestrial tetrapods and ray-finned fishes reveals broad homology between the swim bladder and lung cell types as well as shared and idiosyncratic changes of the external gills of the West African lungfish and the internal gills of Atlantic salmon. The single-cell atlas presented here provides a valuable resource for further exploration of the respiratory system evolution in vertebrates and the diversity of lungfish terrestrialization.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits.

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

The complete and fully-phased diploid genome of a male Han Chinese.

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.