Chromosome-level Genome Assembly Research Articles

Chromosome-Scale Genome Assembly Provides Insights into Fresh Pine Wood Decay Strategies of the Wolfiporia hoelen

The sclerotia of Wolfiporia hoelen (Fr.) Y.C. Dai & V. Papp is an important traditional Chinese medicine with diverse pharmacological properties. This study utilized a combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing to generate a high-quality chromosome-level genome assembly of a W. hoelen strain Minling A5. There were 112 contigs in the genome, with 62.95 Mb in total length and 4.21 Mb in length for the contig N50. The average GC content was 51.89%. Based on Hi-C data, we corrected the CCS data and scaffolded them into 14 pseudo-chromosomes. The genome contained 44.37% repetitive sequences and 12,670 protein-coding genes, 86.53% (10,963) of which could be functionally annotated in at least one of the KOG, GO, Pfam, Swissprot, TrEMBL, NR, and KEGG databases. In addition, 240 transfer RNAs, 97 ribosomal RNAs, and 103 other non-coding RNAs were identified in the W. hoelen genome. A total of 755 pseudogenes were also identified, with an average length of 2665.51 bp. Further, there were 398, 100, 2837, 519, and 2068 genes annotated by CAZymes, TCDB, PHI, P450, and DFVF databases, respectively. One notable attribute of W. hoelen is its capacity to thrive in a substrate of fresh pine sawdust. Through an analysis of the growth on various pure wood sawdust culture media, we found that the growth of W. hoelen and Sparassis latifolia on pine sawdust was similar to that on broad-leaved wood sawdust, while the growth of Pleurotus ostreatus, P. eryngii, and Cyclocybe aegerita was slower than that on broad-leaved wood sawdust. By the functional annotation analysis of orthogroups in these five mushroom-forming fungi, it was determined that 645 orthogroups were specifically common in W. hoelen and S. latifolia. The genes in these specific orthogroups were significantly enriched in 12 pathways, including steroid biosynthesis, biosynthesis of antibiotics, and tyrosine metabolism. The high-quality genome and comparative genome analysis results significantly contribute to advancing our foundational knowledge of W. hoelen biology, while also offering valuable insights for the development of innovative biotechnological approaches aimed at enhancing the efficient and sustainable utilization of Pinus.

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri.

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation of orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

Chromosome-scale genome assemblies of sexually dimorphic male and female Acrossocheilus fasciatus

Acrossocheilus fasciatus is a stream-dwelling fish species of the Barbinae subfamily. It is valued for its colorfully striped appearance and delicious meat. This species is also characterized by apparent sexual dimorphism and toxic ovum. Biology and aquaculture researches of A. fasciatus are hindered by the lack of a high-quality reference genome. Here, we report chromosome-level genome assemblies of the male and female A. fasciatus. The HiFi-only genome assemblies for both female and male individuals were 899.13 Mb (N50 length of 32.58 Mb) and 885.68 Mb (N50 length of 33.06 Mb), respectively. Notably, a substantial proportion of the assembled sequences, accounting for 96.15% and 98.35% for female and male genomes, respectively, were successfully anchored onto 25 chromosomes utilizing Hi-C data. We annotated the female assembly as a reference genome and identified a total of 400.62 Mb (44.56%) repetitive sequences, 27,392 protein-coding genes, and 35,869 ncRNAs. The high-quality male and female reference genomes will provide genomic resources for developing sex-specific molecular markers, inform single-sex breeding, and elucidate genetic mechanisms of sexual dimorphism.

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis.

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Chromosome-level genome assembly of ridgetail white shrimp Exopalaemon carinicauda

Exopalaemon carinicauda, a eurythermal and euryhaline shrimp, contributes one third of the total biomass production of polyculture ponds in eastern China and is considered as a potential ideal experimental animal for research on crustaceans. We conducted a high-quality chromosome-level genome assembly of E. carinicauda combining PacBio HiFi and Hi-C sequencing data. The total assembly size was 5.86 Gb, with a contig N50 of 235.52 kb and a scaffold N50 of 138.24 Mb. Approximately 95.29% of the assembled sequences were anchored onto 45 pseudochromosomes. BUSCO analysis revealed that 92.89% of 1,013 single-copy genes were highly conserved orthologs. A total of 44, 288 protein-coding genes were predicted, of which 70.53% were functionally annotated. Given its high heterozygosity (2.62%) and large proportion of repeat sequences (71.49%), it is one of the most complex genome assemblies. This chromosome-scale genome will be a valuable resource for future molecular breeding and functional genomics research on E. carinicauda.

Chromosome-level genome assembly of Solanum pimpinellifolium

Solanum pimpinellifolium, the closest wild relative of the domesticated tomato, has high potential for use in breeding programs aimed at developing multi-pathogen resistance and quality improvement. We generated a chromosome-level genome assembly of S. pimpinellifolium LA1589, with a size of 833 Mb and a contig N50 of 31 Mb. We anchored 98.80% of the contigs into 12 pseudo-chromosomes, and identified 74.47% of the sequences as repetitive sequences. The genome evaluation revealed BUSCO and LAI score of 98.3% and 14.49, respectively, indicating high quality of this assembly. A total of 41,449 protein-coding genes were predicted in the genome, of which 89.17% were functionally annotated. This high-quality genome assembly serves as a valuable resource for accelerating the biological discovery and molecular breeding of this important horticultural crop.

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution.

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

High-density linkage to physical mapping in a unique Tall × Dwarf coconut (Cocos nucifera L.) outbred F2 uncovers a major QTL for flowering time colocalized with the FLOWERING LOCUS T (FT).

The coconut tree crop (Cocos nucifera L.) provides vital resources for millions of people worldwide. Coconut germplasm is largely classified into 'Tall' (Typica) and 'Dwarf' (Nana) types. While Tall coconuts are outcrossing, stress tolerant, and late flowering, Dwarf coconuts are inbred and flower early with a high rate of bunch emission. Precocity determines the earlier production of a plantation and facilitates management and harvest. A unique outbred F2 population was used, generated by intercrossing F1 hybrids between Brazilian Green Dwarf from Jiqui (BGDJ) and West African Tall (WAT) cultivars. Single-nucleotide polymorphism (SNP) markers fixed for alternative alleles in the two varieties, segregating in an F2 configuration, were used to build a high-density linkage map with ~3,000 SNPs, anchored to the existing chromosome-level genome assemblies, and a quantitative trait locus (QTL) mapping analysis was carried out. The linkage map established the chromosome numbering correspondence between the two reference genome versions and the relationship between recombination rate, physical distance, and gene density in the coconut genomes. Leveraging the strong segregation for precocity inherited from the Dwarf cultivar in the F2, a major effect QTL with incomplete dominance was mapped for flowering time. FLOWERING LOCUS T (FT) gene homologs of coconut previously described as putatively involved in flowering time by alternative splice variant analysis were colocalized within a ~200-kb window of the major effect QTL [logarithm of the odds (LOD) = 11.86]. Our work provides strong phenotype-based evidence for the role of the FT locus as the putative underlying functional variant for the flowering time difference between Dwarf and Tall coconuts. Major effect QTLs were also detected for developmental traits of the palm, plausibly suggesting pleiotropism of the FT locus for other precocity traits. Haplotypes of the two SNPs flanking the flowering time QTL inherited from the Dwarf parent BGDJ caused a reduction in the time to flower of approximately 400 days. These SNPs could be used for high-throughput marker-assisted selection of early-flowering and higher-productivity recombinant lines, providing innovative genetic material to the coconut industry.

A chromosome-level genome assembly of Chinese quince (Pseudocydonia sinensis).

Pseudocydonia sinensis, also known as Chinese quince, is a perennial shrub or small tree highly valued for its edibility and medicinal properties. This study presents the first chromosome-level genome assembly of P. sinensis, achieved using HiFi sequencing and Hi-C scaffolding technology. The assembly resulted in a high-quality genome of 576.39 Mb in size. The genome was anchored to 17 pseudo-chromosomes, with a contig N50 of 27.6 Mb and a scaffold N50 of 33.8 Mb. Comprehensive assessment using BUSCO, CEGMA and BWA tools indicates the high completeness and accuracy of the genome assembly. Our analysis identified 116 species-specific genes, 1196 expanded genes and 1109 contracted genes. Additionally, the distribution of 4DTv values suggests that the most recent duplication event occurred before the divergence of P. sinensis from both Chaenomeles pinnatifida and Pyrus pyrifolia. The assembly of this high-quality genome provides a valuable platform for the genetic breeding and cultivation of P. sinensis, as well as for the comparison of the genetic complexity of P. sinensis with other important crops in the Rosaceae family.

Chromosome-level genome assembly and annotation of the cold-water species Ophiura sarsii

The cold-water species Ophiura sarsii, a brittle star, is a key echinoderm in the Arctic continental shelf region, highly sensitive to climate change. However, the absence of a high-quality genome has hindered a thorough understanding of its adaptive evolution. In this study, we reported the first chromosome-level genome assembly of O. sarsii. The genome assembly totalled 1.57 Gb, encompassing 19 chromosomes with a GC content of 37.11% and a scaffold N50 length of 78.03 Mb. The Benchmarking Universal Single-Copy Orthologs (BUSCO) assessment yielded a completeness estimate of 93.5% for this assembly. We predicted a total of 27,099 protein-coding genes, with 25,079 functionally annotated. The genome was comprised of 58.09% transposable elements. This chromosome-level genome of O. sarsii contributes to our understanding of the origin and evolution of marine organisms.

The chromosome-level genomes of the herbal magnoliids Warburgia ugandensis and Saururus chinensis

Warburgia ugandensis and Saururus chinensis are two of the most important medicinal plants in magnoliids and are widely utilized in traditional Kenya and Chinese medicine, respectively. The absence of higher-quality reference genomes has hindered research on the medicinal compound biosynthesis mechanisms of these plants. We report the chromosome-level genome assemblies of W. ugandensis and S. chinensis, and generated 1.13 Gb and 0.53 Gb genomes from 74 and 27 scaffolds, respectively, using BGI-DIPSEQ, Nanopore, and Hi-C sequencing. The scaffold N50 lengths were 82.97 Mb and 48.53 Mb, and the assemblies were anchored to 14 and 11 chromosomes of W. ugandensis and S. chinensis, respectively. In total, 24,739 and 20,561 genes were annotated, and 98.5% and 98% of the BUSCO genes were fully represented, respectively. The chromosome-level genomes of W. ugandensis and S. chinensis will be valuable resources for understanding the genetics of these medicinal plants, studying the evolution of magnoliids and angiosperms and conserving plant genetic resources.

The chromosome-level genome assembly of the dwarfing apple interstock Malus hybrid ‘SH6’

Malus hybrid ‘SH6’ (M. honanensis × M. domestica)is a commonly used apple interstock in China, known for its excellent dwarfing characteristics and cold tolerance. In this study, a combined strategy utilizing PacBio HiFi, Hi-C and parental resequencing data were employed to assemble two haploid genomes for ‘SH6’. After chromosome anchoring, the final hapH genome size was 596.63 Mb, with a contig N50 of 34.38 Mb. The hapR genome was 649.37 Mb, with a contig N50 of 36.84 Mb. Further analysis predicted that repeated sequences made up 59.69% and 62.52% of the entire genome, respectively. Gene annotations revealed 45,435 genes for hapH and 48,261 genes for hapR. Combined with genomic synteny we suggest that the hapR genome originates from its maternal parent M. domestica cv. Ralls Janet, while the hapH genome comes from its paternal parent, M. honanensis. The assembled genome significantly contributes to the discovery of genes associated with apple dwarfing and the molecular mechanisms governing them.

Aroma formation in single- and double-petal jasmines (Jasminum sambac) during flowering via volatile metabolome and transcriptome

Jasminum sambac (L.) Aiton exists in single-petal (SP), double-petal (DP), and multi-petal (MP) phenotypes, which have a sweet floral scent and release aroma when the flowers open. A total of 103 volatiles were discovered from SP and DP jasmines. 16 volatiles were identified as their key odorants. Caryophyllene, geraniol, linalool, α-farnesene, β-ocimene, methyl anthranilate, methyl salicylate, benzeneacetaldehyde and hexanal may form the overall aroma of those two jasmines during flowering. High-resolution transcriptional profiles of those two jasmines at four flowering stages were generated from our recently reported chromosome-level genome assembly, from which we identified 35 differentially expressed genes (DEGs) involved in terpenoid metabolic pathways and 10 DEGs involved in the phenylpropane/benzene metabolic pathway, which may shape the characteristic aroma of those two jasmines at the molecular level. Through co-expression network analysis of key odorants, we found that bZIP8, bHLH6, MYB4, and AP2/ERF3 may contribute to the regulation of aroma volatiles in those two jasmines. These results offer an important theoretical foundation for studying the regulation of aroma formation for SP and DP jasmines.

An improved chromosome-level genome assembly and annotation of Echeneis naucrates

Echeneis naucrates, as known as live sharksucker, is famous for the behavior of attaching to hosts using a highly modified dorsal fin with oval-shaped sucking disc. Here, we generated an improved high-quality chromosome-level genome assembly of E. naucrates using Illumina short reads, PacBio long reads and Hi-C data. Our assembled genome spans 572.85 Mb with a contig N50 of 23.19 Mb and is positioned to 24 pseudo-chromosomes. Additionally, at least one telomere was identified for 23 out of 24 chromosomes. Furthermore, we identified a total of 22,161 protein-coding genes, of which 21,402 genes (96.9%) were annotated successfully with functions. The combination of ab initio predictions and Repbase-based searches revealed that 15.57% of the assembled E. naucrates genome was identified as repetitive sequences. The completeness of the genome assembly and the gene annotation were estimated to be 97.5% and 95.4% with BUSCO analyses. This work enhances the utility of the live sharksucker genome and provides a valuable groundwork for the future study of genomics, biology and adaptive evolution in this species.

Gene gain and loss from the Asian corn borer W chromosome.

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function. Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad. Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

Systematic characterization of gene families and functional analysis of PvRAS3 and PvRAS4 involved in rosmarinic acid biosynthesis in Prunella vulgaris.

Prunella vulgaris is an important material for Chinese medicines with rosmarinic acid (RA) as its index component. Based on the chromosome-level genome assembly we obtained recently, 51 RA biosynthesis-related genes were identified. Sequence feature, gene expression pattern and phylogenetic relationship analyses showed that 17 of them could be involved in RA biosynthesis. In vitro enzymatic assay showed that PvRAS3 catalyzed the condensation of p-coumaroyl-CoA and caffeoyl-CoA with pHPL and DHPL. Its affinity toward p-coumaroyl-CoA was higher than caffeoyl-CoA. PvRAS4 catalyzed the condensation of p-coumaroyl-CoA with pHPL and DHPL. Its affinity toward p-coumaroyl-CoA was lower than PvRAS3. UPLC and LC-MS/MS analyses showed the existence of RA, 4-coumaroyl-3',4'-dihydroxyphenyllactic acid, 4-coumaroyl-4'-hydroxyphenyllactic acid and caffeoyl-4'-hydroxyphenyllactic acid in P. vulgaris. Generation and analysis of pvras3 homozygous mutants showed significant decrease of RA, 4-coumaroyl-3',4'-dihydroxyphenyllactic acid, 4-coumaroyl-4'-hydroxyphenyllactic acid and caffeoyl-4'-hydroxyphenyllactic acid and significant increase of DHPL and pHPL. It suggests that PvRAS3 is the main enzyme catalyzing the condensation of acyl donors and acceptors during RA biosynthesis. The role of PvRAS4 appears minor. The results provide significant information for quality control of P. vulgaris medicinal materials.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae).

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Chromosome-level genome assembly of the yeast Lodderomyces beijingensis reveals the genetic nature of metabolic adaptations and identifies subtelomeres as hotspots for amplification of mating type loci.

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterized its genome sequence.

The Spartina alterniflora genome sequence provides insights into the salt-tolerance mechanisms of exo-recretohalophytes.

Spartina alterniflora is an exo-recretohalophyte Poaceae species that is able to grow well in seashore, but the genomic basis underlying its adaptation to salt tolerance remains unknown. Here, we report a high-quality, chromosome-level genome assembly of S. alterniflora constructed through PacBio HiFi sequencing, combined with high-throughput chromosome conformation capture (Hi-C) technology and Illumina-based transcriptomic analyses. The final 1.58 Gb genome assembly has a contig N50 size of 46.74 Mb. Phylogenetic analysis suggests that S. alterniflora diverged from Zoysia japonica approximately 21.72 million years ago (MYA). Moreover, whole-genome duplication (WGD) events in S. alterniflora appear to have expanded gene families and transcription factors relevant to salt tolerance and adaptation to saline environments. Comparative genomics analyses identified numerous species-specific genes, significantly expanded genes and positively selected genes that are enriched for 'ion transport' and 'response to salt stress'. RNA-seq analysis identified several ion transporter genes including the high-affinity K+ transporters (HKTs), SaHKT1;2, SaHKT1;3 and SaHKT1;8, and high copy number of Salt Overly Sensitive (SOS) up-regulated under high salt conditions, and the overexpression of SaHKT2;4 in Arabidopsis thaliana conferred salt tolerance to the plant, suggesting specialized roles for S. alterniflora to adapt to saline environments. Integrated metabolomics and transcriptomics analyses revealed that salt stress activate glutathione metabolism, with differential expressions of several genes such as γ-ECS, GSH-S, GPX, GST and PCS in the glutathione metabolism. This study suggests several adaptive mechanisms that could contribute our understanding of evolutional basis of the halophyte.