Chromosome-level genome assembly of the glass catfish(Kryptopterus vitreolus)reveals molecular clues to its transparent phenotype

Glass catfish(Kryptopterus vitreolus)are notable in the aquarium trade for their highly transparent body pattern.This transparency is due to the loss of most reflective iridophores and light-absorbing melanophores in the main body,although certain black and silver pigments remain in the face and head.To date,however,the molecular mechanisms underlying this transparent phenotype remain largely unknown.To explore the genetic basis of this transparency,we constructed a chromosome-level haplotypic genome assembly for the glass catfish,encompassing 32 chromosomes and 23344 protein-coding genes,using PacBio and Hi-C sequencing technologies and standard assembly and annotation pipelines.Analysis revealed a premature stop codon in the putative albinism-related tyrp1b gene,encoding tyrosinase-related protein 1,rendering it a nonfunctional pseudogene.Notably,a synteny comparison with over 30 other fish species identified the loss of the endothelin-3(edn3b)gene in the glass catfish genome.To investigate the role of edn3b,we generated edn3b^(−/−)mutant zebrafish,which exhibited a remarkable reduction in black pigments in body surface stripes compared to wild-type zebrafish.These findings indicate that edn3b loss contributes to the transparent phenotype of the glass catfish.Our high-quality chromosome-scale genome assembly and identification of key genes provide important molecular insights into the transparent phenotype of glass catfish.These findings not only enhance our understanding of the molecular mechanisms underlying transparency in glass catfish,but also offer a valuable genetic resource for further research on pigmentation in various animal species.

Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs

Background Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones.Nevertheless,the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear.Results Here,we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical,high-altitude,and frigid regions.Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds.Analysis of selection signals revealed regions under selection for adaptation in tropical(55.5 Mb),high-altitude(43.6 Mb),and frigid(17.72 Mb)regions.The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments,including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments.Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation.These pathways included blood circulation,protein degradation,and inflammation for adaptation to tropical environments;heart and lung development,hypoxia response,and DNA damage repair for high-altitude adaptation;and thermogenesis,cold-induced vasodilation(CIVD),and the cell cycle for adaptation to frigid environments.By examining the chromatin state of the selection signatures,we identified the lung and ileum as two candidate functional tissues for environmental adaptation.Finally,we identified a mutation(chr1:G246,175,129A)in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments.Conclusions In this study,we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical,high-altitude,and frigid environments.Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

The Genetic Association between CDKN1A and Heart Failure:Genome-Wide Exploration of m6A-SNPs and Mendelian Randomization

Objective N6-methyladenosine(m^(6)A)is a common epigenetic modification in eukaryotes.In this study,we explore the potential impact of m^(6)A-associated single nucleotide polymorphisms(m^(6)A-SNPs)on heart failure(HF).Methods Data from genome-wide association studies(GWAS)investigating HF in humans and from m^(6)A-SNPs datasets were used to identify HF-associated m^(6)A-SNPs.Their functions were explored using expression quantitative trait locus(eQTL),gene expression,and gene enrichment analyses.Mediation protein quantitative trait locus(pQTL)-Mendelian randomization(MR)was used to investigate the potential mechanism between critical protein levels and risk factors for HF.Results We screened 44 HF-associated m^(6)A-SNPs,including 10 m^(6)A-SNPs that showed eQTL signals and differential expressions in HF.The SNP rs1801270 in CDKN1A showed the strongest association with HF(P=7.75×10^(−6)).Additionally,MR verified the genetic association between the CDKN1A protein and HF,as well as the mediating effect of blood pressure(BP)in this pathway.Higher circulating level of CDKN1A was associated with a lower risk of HF(odds ratio[OR]=0.82,95%confidence interval[CI]:0.69 to 0.99).The proportions of hypertension,systolic BP,and diastolic BP were 48.10%,28.94%,and 18.02%,respectively.Associations of PDIA6(P=1.30×10^(−2))and SMAD3(P=4.80×10^(−2))with HF were also detected.Conclusion Multiple HF-related m^(6)A-SNPs were identified in this study.Genetic associations of CDKN1A and other proteins with HF and its risk factors were demonstrated,providing new ideas for further exploration of the molecular mechanisms of HF.

Genome-wide association study identifies 12 new genetic loci associated with growth traits in pigs

Growth traits are among the most important economic traits in pigs and are regulated by polygenes with complex regulatory mechanisms.As the major indicators of growth performance,the backfat thickness(BFT),loin eye area(LEA),and days to 100 kg(D100)traits are commonly used to the genetics improvement in pigs.However,the available genetic markers for these traits are limited.To uncover novel loci and candidate genes associated with growth performance,we collected the phenotypic information of BFT,LEA,and D100 in 1,186 pigs and genotyped all these individuals using the Neogen GGP porcine 80K BeadChip.We performed a genome-wide association study(GWAS)using 4 statistical models,including mixed linear models(MLM),fixed and random model circulating probability unification(FarmCPU),settlement of MLM under progressively exclusive relationships(SUPER),Bayesian-information and linkage-disequilibrium Iteratively nested keyway(Blink),and identified 5,3,and 6 high-confidence single nucleotide polymorphisms(SNPs)associated with BFT,LEA,and D100,respectively.Variant annotation and quantitative trait locus(QTL)mapping analysis suggested that6 genes(SKAP2,SATB1,PDE7B,PPP1R16B,WNT3,and WNT9B)were potentially associated with growth performance in pigs.Transcriptome analysis suggested that the expression of Src Kinase Associated Phosphoprotein 2(SKAP2)was higher in prenatal muscles than in postnatal muscles,and the expression of Phosphodiesterase 7B(PDE7B)continuously increased during the prenatal stages and gradually decreased after birth,implying their potential roles in prenatal skeletal muscle development.Overall,this study provides new candidate loci and genes for the genetic improvement of pigs.

Improved genome annotation of Brassica oleracea highlights the importance of alternative splicing

Brassica oleracea has been developed into many important crops,including cabbage,kale,cauliflower,broccoli and so on.The genome and gene annotation of cabbage(cultivar JZS),a representative morphotype of B.oleracea,has been widely used as a common reference in biological research.Although its genome assembly has been updated twice,the current gene annotation still lacks information on untranslated regions(UTRs)and alternative splicing(AS).Here,we constructed a high-quality gene annotation(JZSv3)using a full-length transcriptome acquired by nanopore sequencing,yielding a total of 59452 genes and 75684 transcripts.Additionally,we re-analyzed the previously reported transcriptome data related to the development of different tissues and cold response using JZSv3 as a reference,and found that 3843 out of 11908 differentially expressed genes(DEGs)underwent AS during the development of different tissues and 309 out of 903 cold-related genes underwent AS in response to cold stress.Meanwhile,we also identified many AS genes,including BolLHCB5 and BolHSP70,that displayed distinct expression patterns within variant transcripts of the same gene,highlighting the importance of JZSv3 as a pivotal reference for AS analysis.Overall,JZSv3 provides a valuable resource for exploring gene function,especially for obtaining a deeper understanding of AS regulation mechanisms.

Haplotype-resolved and chromosome-scale genomes provide insights into co-adaptation between the Amur tiger and Amur leopard

DEAR EDITOR,Big cats,such as Amur tigers(Panthera tigris altaica)and Amur leopards(P.pardus orientalis),are apex predator and have evolved specialized traits for hunting and carnivory(Moya et al.,2022),thus playing a crucial role in maintaining biodiversity and ecosystem integrity by regulating prey-predator dynamics.However,human-induced pressures,habitat fragmentation,and environmental alterations have restricted these species in small and isolated populations.Currently,all extant big cats are categorized as endangered or threatened according to their conservation status.Amur tigers and Amur leopards share overlapping geographic ranges,habitats,and certain prey species in the forests of Northeast Asia(Jiang et al.,2015).To reduce interspecies conflict,these carnivores exhibit differentiated dietary and temporal niches.Amur tigers predominantly prey on large ungulates,while Amur leopards hunt small to medium-sized animals(Sugimoto et al.,2016).Additionally,they occupy different temporal niches,with tigers being active at night and leopards more active during the day.Despite spatial and temporal niche partitioning,interspecific competition between these two species is inevitable.Tigers,benefiting from their greater size,have a competitive advantage over leopards,which can manifest in occasional leopard predation by tigers and declines in leopard populations with increasing tiger density(Jiang et al.,2015).Tigers also displace leopards from marginal habitats in nature reserves where they coexist.

Zea mays(L.) P1 locus for cob glume color identified as a post-domestication selection target with an effect on temperate maize genomes

Artificial selection during domestication and post-domestication improvement results in loss of genetic diversity near target loci. However, the genetic locus associated with cob glume color and the nature of the genomic pattern surrounding it was elusive and the selection effect in that region was not clear. An association mapping panel consisting of 283 diverse modern temperate maize elite lines was genotyped by a chip containing over 55,000 evenly distributed SNPs. Ten-fold resequencing at the target region on 40 of the panel lines and 47 tropical lines was also undertaken. A genome-wide association study(GWAS) for cob glume color confirmed the P1 locus, which is located on the short arm of chromosome 1, with a-log10 P value for surrounding SNPs higher than the Bonferroni threshold(α/n, α < 0.001) when a mixed linear model(MLM) was implemented. A total of 26 markers were identified in a 0.78 Mb region surrounding the P1 locus, including 0.73 Mb and 0.05 Mb upstream and downstream of the P1 gene, respectively. A clear linkage disequilibrium(LD) block was found and LD decayed very rapidly with increasing physical distance surrounding the P1 locus. The estimates of π and Tajima's D were significantly(P < 0.001) lower at both ends compared to the locus. Upon comparison of temperate and tropical lines at much finer resolution by resequencing(180-fold finer than chip SNPs), a more structured LD block pattern was found among the 40 resequenced temperate lines. All evidence indicates that the P1 locus in temperate maize has not undergone neutral evolution but has been subjected to artificial selection during post-domestication selection or improvement. The information and analytical results generated in this study provide insights as to how breeding efforts have affected genome evolution in crop plants.

Genome engineering using the CRISPR/Cas system

Recently, an epoch-making genome engineering technology using clustered regularly at interspaced short palindromic repeats(CRISPR) and CRISPR associated(Cas) nucleases, was developed. Previous technologies for genome manipulation require the time-consuming design and construction of genome-engineered nucleases for each target and have, therefore, not been widely used in mouse research where standard techniques based on homologous recombination are commonly used. The CRISPR/Cas system only requires the design of sequences complementary to a target locus, making this technology fast and straightforward. In addition, CRISPR/Cas can be used to generate mice carrying mutations in multiple genes in a single step, an achievement not possible using other methods. Here, we review the uses of this technology in genetic analysis and manipulation, including achievements made possible to date and the prospects for future therapeutic applications.

Complete genome sequence of the rifamycin SV-producing Amycolatopsis mediterranei U32 revealed its genetic characteristics in phylogeny and metabolism

Amycolatopsis mediterranei is used for industry-scale production of rifamycin, which plays a vital role in antimyco- bacterial therapy. As the first sequenced genome of the genus Amycolatopsis, the chromosome of strain U32 comprising 10 236 715 base pairs, is one of the largest prokaryotic genomes ever sequenced so far. Unlike the linear topology found in streptomycetes, this chromosome is circular, particularly similar to that of Saccharopolyspora erythraea and Nocardia farcinica, representing their close relationship in phylogeny and taxonomy. Although the predicted 9 228 protein-coding genes in the A. mediterranei genome shared the greatest number of orthologs with those of S. erythraea, it was unexpectedly followed by Streptomyces coelicolor rather than N. farcinica, indicating the distinct metabolic characteristics evolved via adaptation to diverse ecological niches. Besides a core region analogous to that common in streptomycetes, a novel ＇quasicore＇ with typical core characteristics is defined within the non-core region, where 21 out of the total 26 gene clusters for secondary metabolite production are located. The rifamycin biosynthesis gene cluster located in the core encodes a cytochrome P450 enzyme essential for the conversion of rifamycin SV to B, revealed by comparing to the highly homologous cluster of the rifamycin B-producing strain S699 and further confirmed by genetic complementation. The genomic information of A. mediterranei demonstrates a metabolic network orchestrated not only for extensive utilization of various carbon sources and inorganic nitrogen compounds but also for effective funneling of metabolic intermediates into the secondary antibiotic synthesis process under the control of a seemingly complex regulatory mechanism.

Metagenome of microorganisms associated with the toxic Cyanobacteria Microcystis aeruginosa analyzed using the 454 sequencing platform

In this study, the 454 pyrosequencing technology was used to analyze the DNA of the Microcystis aeruginosa symbiosis system from cyanobacterial algal blooms in Taihu Lake, China. We generated 183 228 reads with an average length of 248 bp. Running the 454 assembly algorithm over our sequences yielded 22 239 significant contigs. After excluding the M. aeruginosa sequences, we obtained 1 322 assembled contigs longer than 1 000 bp. Taxonomic analysis indicated that four kingdoms were represented in the community: Archaea (n = 9; 0.01%), Bacteria (n = 98 921; 99.6%), Eukaryota (n = 373; 3.7%), and Viruses (n = 18; 0.02%). The bacterial sequences were predominantly Alphaproteobacteria (n = 41 805; 83.3%), Betaproteobacteria (n = 5 254; 10.5%) and Gammaproteobacteria (n = 1 180; 2.4%). Gene annotations and assignment of COG (clusters of orthologous groups) functional categories indicate that a large number of the predicted genes are involved in metabolic, genetic, and environmental information processes. Our results demonstrate the extraordinary diversity of a microbial community in an ectosymbiotic system and further establish the tremendous utility of pyrosequencing.

Genetic dissection of rice appearance quality and cooked rice elongation by genome-wide association study

Appearance and cooked rice elongation are key quality traits of rice. Although some QTL for these traits have been identified, understanding of the genetic relationship between them remains limited. In the present study, large phenotypic variation was observed in 760 accessions from the 3 K Rice Genomes Project for both appearance quality and cooked rice elongation. Most component traits of appearance quality and cooked rice elongation showed significant pairwise correlations, but a low correlation was found between appearance quality and cooked rice elongation. A genome-wide association study identified 74 QTL distributed on all 12 chromosomes for grain length, grain width, length to width ratio, degree of endosperm with chalkiness, rice elongation difference, and elongation index. Thirteen regions containing QTL stably expressed in multiple environments and/or exerting pleiotropic effects on multiple traits were detected. By gene-based association analysis and haplotype analysis, 46 candidate genes, including five cloned genes, and 49 favorable alleles were identified for these 13 QTL. The effect of the candidate gene Wx on rice elongation difference was validated by a transgenic strategy. These results shed light on the genetic bases of appearance quality and cooked rice elongation and provide gene resources for improving rice quality by molecular breeding.

Sequencing the Cotton Genomes-Gossypium spp.

The genomes of most major crops,including cotton,will be fully sequenced in the next few years.Cotton is unusual,although not unique,in that we will need to sequence not only

Whole Genome Analysis Reveals New Insights into Macrolide Resistance in Mycoplasma pneumoniae

Objective Mutations in 23 S rRNA gene are known to be associated with macrolide resistance in Mycoplasma pneumoniae（M. pneumoniae）. However, these mutations alone do not fully explain the high resistance rates in Asia. The aim of this study was to investigate other possible mutations involved in macrolide resistance in M. pneumoniae. Methods The whole genomes of 10 clinical isolates of M. pneumoniae with macrolide resistance were sequenced by Illumina Hi Seq2000 platform. The role of the macrolide-specific efflux transporter was assessed by efflux-pump inhibition assays with reserpine and carbonyl cyanide m-chlorophenyl-hydrazone（CCCP）. Results A total of 56 single nucleotide polymorphisms（SNPs） were identified in 10 clinical isolates in comparison to the reference strains M129 and FH. Strikingly, 4 of 30 SNPs causing non-synonymous mutations were clustered in macrolide-specific efflux system gene mac B encoding macrolide-specific efflux pump protein of the ATP-binding cassette transporter family. In assays of the minimal inhibitory concentrations（MIC） of macrolide antibiotics in the presence of the efflux pump inhibitors caused a significant decrease of MICs, even under detectable levels in some strains. Conclusion Our study suggests that macrolide efflux pump may contribute to macrolide resistance in M. pneumoniae in addition to the common point mutations in 23 S r RNA gene.

Endogenous viral elements in algal genomes

Endogenous viral elements （EVEs） are host-genomic fragments originated from viral genomes. They have been found universally in animal and plant genomes. Here we carried out a systematic screening and analy-sis of EVEs in algal genomes and found that EVEs commonly exist in algal genomes. We classified the EVE fragments into three categories according to the length of EVE fragments. Due to the probability of sequence similarity by chance, we ignored the potential function of medium-length EVE fragments. However, long-length EVE fragments probably had capability to encode protein domains or even entire proteins, and some short-length EVE fragments had high similarity with host＇s siRNA sequences and possibly served functions of small RNAs. Therefore, short and long EVE fragments might provide regulomic and proteomic novelty to the host＇s metabolism and adaptation. We also found several EVE fragments shared by more than 3 algal genomes. By phylogenetic analysis of the shared EVEs and their corresponding species, we found that the integration of viral fragments into host genomes was an ancient event, possibly before the divergence of Chlorophytes and Ochrophytes. Our findings show that there is a frequent genetic flow from viruses to algal genomes. Moreover, study on algal EVEs shed light on the virus-host interaction in large timescale and could also help us understand the balance of marine ecosystems.

Human myoblast genome therapy

Human Myoblast Genome Therapy (HMGT) is a platform technology of cell transplantation, nuclear transfer, and tissue engineering. Unlike stem cells, myoblasts are differentiated, immature cells destined to become muscles. Myoblasts cultured from satellite cells of adult muscle biopsies survive, develop, and function to revitalize degenerative muscles upon transplantation. Injection injury activates regeneration of host myofibers that fuse with the engrafted myoblasts, sharing their nuclei in a common gene pool of the syncytium. Thus, through nuclear transfer and complementation, the normal human genome can be transferred into muscles of patients with genetic disorders to achieve phenotype repair or disease prevention. Myoblasts are safe and efficient gene transfer vehicles endogenous to muscles that constitute 50% of body weight. Results of over 280 HMGT procedures on Duchenne Muscular Dystrophy (DMD) subjects in the past 15 years demonstrated absolute safety. Myoblast-injected DMD muscles showed improved histology. Strength increase at 18 months post-operatively averaged 123%. In another application of HMGT on ischemic cardiomyopathy, the first human myoblast transfer into porcine myocardium revealed that it was safe and effective. Clinical trials on approximately 220 severe cardiomyopathy patients in 15 countries showed a <10% mortality. Most subjects received autologous cells implanted on the epicardial surface during coronory artery bypass graft, or injected on the endomyocardial surface percutaneously through guiding catheters. Significant increases in left ventricular ejection fraction, wall thickness, and wall motion have been reported, with reduction in perfusion defective areas, angina, and shortness of breath. As a new modality of treatment for disease in the skeletal muscle or myocardium, HMGT emerged as safe and effective. Large randomized multi-center trials are under way to confirm these preliminary results. The future of HMGT is bright and exciting.

Genome assembly of the maize inbred line A188 provides a new reference genome for functional genomics

The current assembled maize genomes cannot represent the broad genetic diversity of maize germplasms.Acquiring more genome sequences is critical for constructing a pan-genome and elucidating the linkage between genotype and phenotype in maize.Here we describe the genome sequence and annotation of A188,a maize inbred line with high phenotypic variation relative to other lines,acquired by single-molecule sequencing and optical genome mapping.We assembled a 2210-Mb genome with a scaffold N50 size of 11.61 million bases(Mb),compared to 9.73 Mb for B73 and 10.2 Mb for Mo17.Based on the B73_Ref Gen_V4 genome,295 scaffolds(2084.35 Mb,94.30%of the final genome assembly)were anchored and oriented on ten chromosomes.Comparative analysis revealed that~30%of the predicted A188 genes showed large structural divergence from B73,Mo17,and W22 genomes,which causes high protein divergence and may lead to phenotypic variation among the four inbred lines.As a line with high embryonic callus(EC)induction capacity,A188 provides a convenient tool for elucidating the molecular mechanism underlying the formation of EC in maize.Combining our new A188 genome with previously reported QTL and RNA sequencing data revealed eight genes with large structural variation and two differentially expressed genes playing potential roles in maize EC induction.

First Principle Material Genome Approach for All Solid-State Batteries

Due to ever-increasing concern about safety issues in using alkali metal ionic batteries, all solid-state batteries (ASSBs) have attracted tremendous attention. The foundation to enable high-performance ASSBs lies in delivering ultra-fast ionic conductors that are compatible with both alkali anodes and high-voltage cathodes. Such a challenging task cannot be fulfilled, without solid understanding covering materials stability and properties, interfacial reactions, structural integrity, and electrochemical windows. Here in this work, we will review recent advances on fundamental modeling in the framework of material genome initiative based on the density functional theory (DFT), focusing on solid alkali batteries. Efforts are made in offering a dependable road chart to formulate competitive materials and construct "better" batteries.

Low-depth whole genome sequencing reveals copy number variations associated with higher pathologic grading and more aggressive subtypes of lung non-mucinous adenocarcinoma

Objective:Histology grade,subtypes and TNM stage of lung adenocarcinomas are useful predictors of prognosis and survival.The aim of the study was to investigate the relationship between chromosomal instability,morphological subtypes and the grading system used in lung non-mucinous adenocarcinoma(LNMA).Methods:We developed a whole genome copy number variation(WGCNV)scoring system and applied next generation sequencing to evaluate CNVs present in 91 LNMA tumor samples.Results:Higher histological grades,aggressive subtypes and more advanced TNM staging were associated with an increased WGCNV score,particularly in CNV regions enriched for tumor suppressor genes and oncogenes.In addition,we demonstrate that 24-chromosome CNV profiling can be performed reliably from specific cell types(<100 cells)isolated by sample laser capture microdissection.Conclusions:Our findings suggest that the WGCNV scoring system we developed may have potential value as an adjunct test for predicting the prognosis of patients diagnosed with LNMA.

Research on Copyright Protection Method of Material Genome Engineering Data Based on Zero-Watermarking

In order to effectively solve the problem of copyright protection of materials genome engineering data,this paper proposes a method for copyright protection of materials genome engineering data based on zero-watermarking technology.First,the important attribute values are selected from the materials genome engineering database;then,use the method of remainder to group the selected attribute values and extract eigenvalues;then,the eigenvalues sequence is obtained by the majority election method;finally,XOR the sequence with the actual copyright information to obtain the watermarking information and store it in the third-party authentication center.When a copyright dispute requires copyright authentication for the database to be detected.First,the zero-watermarking construction algorithm is used to obtain an eigenvalues sequence;then,this sequence is XORed with the watermarking information stored in the third-party authentication center to obtain copyright information to-be-detected.Finally,the ownership is determined by calculating the similarity between copyright information to-be-detected and copyright information that has practical significance.The experimental result shows that the zero-watermarking method proposed in this paper can effectively resist various common attacks,and can well achieve the copyright protection of material genome engineering database.

Nuclear and chloroplast genome diversity revealed by low-coverage whole-genome shotgun sequence in 44 Brassica oleracea breeding lines

Whole-genome shogun sequence(WGS)data generated by next-generation sequencing(NGS)platforms are a valuable resource for crop improvement.We produced 5–6×WGS coverage of 44 Brassica oleracea breeding lines representing seven subspecies/morphotypes:cabbage,broccoli,cauliflower,kailan,kale,Brussels sprout,and kohlrabi to systematically evaluate the nuclear and chloroplast(Cp)diversity in the 44 B.oleracea breeding lines.We then exploited the impact of low-coverage NGS by evaluating nuclear genome diversity and assembly,annotation of complete chloroplast(Cp)genomes and 45 S nuclear ribosomal DNA(45S nrDNA)sequences,and copy number variation for major repeats.Nuclear genome diversity analysis has revealed a total of 496463 SNPs and 37493 indels in the nuclear genome across the 44 accessions.Interestingly,some SNPs showed subspecies enrichment at certain chromosomal regions.The assembly of complete Cp genomes contained 153361–153372 bp with 37 variants including SNPs and indels.The 45S nrDNA transcription unit was 5802 bp long with a total of 31 SNPs from the 44 lines.The phylogenetic tree inferred from the nuclear and Cp genomes coincided and clustered broccoli,cauliflower,and kailan in one group and cabbage,Brussels sprout,kale,and kohlrabi in another group.The morphotypes diverged during the last 0.17 million years.The Cp genome diversity reflected the unique cytoplasm of each subspecies,and revealed that the cytoplasm of many breeding lines was replaced and intermingled via inter-subspecies crosses during the breeding process instead.The polymorphic Cp markers provide a classification system for the cytoplasm types in B.oleracea.Furthermore,copy numbers of major transposable elements(TEs)showed high diversity among the 44 accessions,indicating that many TEs have become active recently.Overall,we demonstrated a comprehensive utilization of low-coverage NGS data and might shed light on the genetic diversity and evolution of diverse B.oleracea morphotypes.