PILER: identification and classification of genomic repeats

Haplotype-resolved genome assembly and genome-wide association study identifies the candidate gene closely related to sugar content and tuber yield in Solanum tuberosum

As an important noncereal food crop grown worldwide, the genetic improvement of potato in tuber yield and quality is largely constrained due to the lacking of a high-quality reference genome and understanding of the regulatory mechanism underlying the formation of superior alleles. Here, a chromosome-scale haplotype-resolved genome assembled from an anther-cultured progeny of 'Ningshu 15', a tetraploid variety featured by its high starch content and drought resistance was presented. The assembled genome size was 1.653 Gb, with a contig N50 of approximately 1.4 Mb and a scaffold N50 of 61 Mb. The long terminal repeat assembly index score of the two identified haplotypes of 'Ningshu 15' was 11.62 and 11.94, respectively. Comparative genomic analysis revealed that positive selection occurred in gene families related to starch, sucrose, fructose and mannose metabolism, and carotenoid biosynthesis. Further genome-wide association study in 141 accessions identified a total number of 53 quantitative trait loci related to fructose, glucose, and sucrose content. Among them, a tonoplast sugar transporter encoding gene, StTST2, closely associated with glucose content was identified. Constitutive expression of StTST2 in potato and Arabidopsis increased the photosynthetic rate, chlorophyll and sugar content, biomass tuber and seed production in transgenic plants. In addition, co-immunoprecipitation assays demonstrated that StTST2 directly interacted with SUT2. Our study provides a high-quality genome assembly and new genetic locus of potato for molecular breeding.

Chromosome-level genome assembly of scalloped spiny lobster Panulirus homarus homarus

Lobsters, aquatic organisms of significant economic value, hold an important position in the global aquaculture and fisheries industries. However, due to overfishing and ecological change, the populations of certain lobster species have declined dramatically, prompting conservation efforts in various countries. However, limited genomics research has restricted our capacity to conserve and exploit lobster germplasm resources. Here, we present a chromosome-level reference genome for Panulirus homarus homarus constructed using PacBio long-read sequencing and Hi-C data. The genome assembly size was 2.61 Gb, with a contig N50 of 5.43 Mb, and a scaffold N50 of 36.69 Mb. The assembled sequences were anchored to 73 chromosomes, covering 96.05% of the total genome. A total of 25,580 protein-coding genes were predicted, and 99.98% of the genes were functionally annotated using various protein databases. The high-quality genome assembly provides a valuable resource for studying the biology and evolutionary history of P. h. homarus, and could facilitate sustainable resource management, aquaculture, and conservation of the species.

New insights into the cold tolerance of upland switchgrass by integrating a haplotype-resolved genome and multi-omics analysis

BACKGROUND

Switchgrass (Panicum virgatum L.) is a bioenergy and forage crop. Upland switchgrass exhibits superior cold tolerance compared to the lowland ecotype, but the underlying molecular mechanisms remain unclear.

RESULTS

Here, we present a high-quality haplotype-resolved genome of the upland ecotype "Jingji31." We then conduct multi-omics analysis to explore the mechanism underlying its cold tolerance. By comparative transcriptome analysis of the upland and lowland ecotypes, we identify many genes with ecotype-specific differential expression, particularly members of the cold-responsive (COR) gene family, under cold stress. Notably, AFB1, ATL80, HOS10, and STRS2 gene families show opposite expression changes between the two ecotypes. Based on the haplotype-resolved genome of "Jingji31," we detect more cold-induced allele-specific expression genes in the upland ecotype than in the lowland ecotype, and these genes are significantly enriched in the COR gene family. By genome-wide association study, we detect an association signal related to the overwintering rate, which overlaps with a selective sweep region containing a cytochrome P450 gene highly expressed under cold stress. Heterologous overexpression of this gene in rice alleviates leaf chlorosis and wilting under cold stress. We also verify that expression of this gene is suppressed by a structural variation in the promoter region.

CONCLUSIONS

Based on the high-quality haplotype-resolved genome and multi-omics analysis of upland switchgrass, we characterize candidate genes responsible for cold tolerance. This study advances our understanding of plant cold tolerance, which provides crop breeding for improved cold tolerance.

Chromosome-level genomes of Arctic and Antarctic mosses: Aulacomnium turgidum and Polytrichastrum alpinum

Bryophytes play a crucial role in the ecosystems of polar regions. These simple plants are among the predominant vegetation types in both Arctic and Antarctic landscapes, where they contribute significantly to biodiversity and ecological stability. Here, we report the chromosome-level genomes of two polar moss species, the Arctic Aulacomnium turgidum and Antarctic Polytrichastrum alpinum. Utilizing a combination of Illumina short reads, Nanopore long reads, and Hi-C data, we assembled genomes of 277.84 Mb for A. turgidum and 498.33 Mb for P. alpinum, respectively. These assemblies were anchored to 11 chromosomes for A. turgidum and 8 chromosomes for P. alpinum. Both species exhibited a sex chromosome with distinct genomic characteristics. Gene annotations revealed 25,999 protein-coding genes in A. turgidum and 28,070 in P. alpinum. The high completeness of the gene space was validated via BUSCO, achieving impressive scores of 98.2% and 98.0%. These high-quality genomes provide critical resources for studying the adaptive evolution and stress tolerance mechanisms of mosses in extreme polar environments.

Comparative genomics provides insights into chromosomal evolution and immunological adaptation in horseshoe bats

Horseshoe bats are natural hosts of zoonotic viruses, yet the genetic basis of their antiviral immunity is poorly understood. Here we generated two new chromosomal-level genome assemblies for horseshoe bat species (Rhinolophus) and three close relatives, and show that, during their diversification, horseshoe bats underwent extensive chromosomal rearrangements and gene expansions linked to segmental duplications. These expansions have generated new adaptive variations in type I interferons and the interferon-stimulated gene ANXA2R, which potentially enhance antiviral states, as suggested by our functional assays. Genome-wide selection screens, including of candidate introgressed regions, uncover numerous putative molecular adaptations linked to immunity, including in viral receptors. By expanding taxon coverage to ten horseshoe bat species, we identify new variants of the SARS-CoV-2 receptor ACE2, and report convergent functionally important residues that could explain wider patterns of susceptibility across mammals. We conclude that horseshoe bats have numerous signatures of adaptation, including some potentially related to immune response to viruses, in genomic regions with diverse and multiscale mutational changes.

Chromosome-level genome of the parthenogenetic booklouse Liposcelis bostrychophila reveals high heterozygosity and a nonhomologous chromosome

Psocodean species are emerging as significant sanitary and stored-product pests, posing threats to human health and global food security. Out of an estimated 10 000 species, the whole genome sequences of only 4 species have been published. Genomic resources are crucial for establishing effective pest control and enhancing our understanding of the evolution of psocodean species. In this study, we employed Illumina and PacBio sequencing along with Hi-C scaffolding techniques to generate a chromosome-level genome assembly for the parthenogenetic booklouse Liposcelis bostrychophila. The assembled genome of this booklouse measures 291.67 Mb in length and comprises 9 chromosomes. Notably, the genome of L. bostrychophila exhibits a high level of heterozygosity and features a distinctive nonhomologous chromosome. This heterozygous characteristic of the parthenogenetic booklouse genome may arise from high mutation rates, based on genomic variations analysis across multiple generations. Our analysis revealed significantly expanded gene families, primarily associated with the detoxification and feeding habits of L. bostrychophila. These include integument esterases (ESTs), ATP-binding cassette (ABC) transporter genes and gustatory receptors (GRs). The high-quality genome sequence of L. bostrychophila provides valuable resources for further study on the molecular mechanisms of stress resistance. It enables researchers to identify crucial functional genes and facilitates research on the population genetics, evolution and phylogeny of booklice.

Chromosome-scale assembly of the Xenocypris davidi using PacBio HiFi reads and Hi-C technologies

Xenocypris davidi is a benthic fish species widely distributed in the water systems south of the Yellow River in China, playing a significant role in aquatic ecosystems. Despite its ecological and economic importance, genomic resources for X. davidi are limited, hindering a comprehensive understanding of its evolutionary adaptations and genetic improvements. This study presents the first chromosome-level genome assembly of X. davidi, utilizing PacBio long-reads, Illumina short reads, and Hi-C sequencing data. The genome assembly spans 1.05 Gb with a scaffold N50 length of 33.99 Mb, and 95.12% of the genome sequence was successfully anchored onto 24 pseudochromosomes. We identified 27,360 protein-coding genes, of which 26,672 were functionally annotated. This genome sequence provides a valuable resource for exploring the molecular basis of agronomic traits in X. davidi and will facilitate its genetic enhancement.

Chromosome-level genome assembly of Eimeria tenella at the single-oocyst level

BACKGROUND

Eimeria are obligate protozoan parasites, and more than 1,500 species have been reported. However, Eimeria genomes lag behind many other eukaryotes since obtaining many oocysts is difficult due to a lack of sustainable in vitro culture, highly repetitive sequences, and mixed species infections. To address this challenge, we used whole-genome amplification of a single oocyst followed by long-read sequencing and obtained a chromosome-level genome of Eimeria tenella.

RESULTS

The assembled genome was 52.13 Mb long, encompassing 15 chromosomes and 46.94% repeat sequences. In total, 7,296 protein-coding genes were predicted, exhibiting high completeness, with 92.00% single-copy BUSCO genes. To the best of our knowledge, this is the first chromosome-level assembly of E. tenella using a combination of single-oocyst whole-genome amplification and long-read sequencing. Comparative genomic and transcriptome analyses confirmed evolutionary relationship and supported estimates of divergence time of apicomplexan parasites and identified AP2 and Myb gene families that may play indispensable roles in regulating the growth and development of E. tenella.

CONCLUSION

This high-quality genome assembly and the established sequencing strategy provide valuable community resources for comparative genomic and evolutionary analyses of the Eimeria clade. Additionally, our study also provides a valuable resource for exploring the roles of AP2 and Myb transcription factor genes in regulating the development of Eimeria parasites.

Decoding Pecan's Fungal Foe: A Genomic Insight into Colletotrichum plurivorum Isolate W-6

Pecan (Carya illinoinensis) is a world-renowned nut crop that is highly favored by consumers for its high content of healthy nutrients. For a long time, anthracnose has severely threatened the yield and quality of pecan, causing significant economic losses to the global pecan industry. Here, we report the 54.57-Mb gapless chromosome-level assembly of the pathogenic ascomycetes Colletotrichum plurivorum isolate W-6 from pecan plantations in Southeast China. Six of 12 chromosomes contain, at least, telomeric repeats (CCCTAA)n or (TTAGGG)n at one end. A total of 14,343 protein-coding genes were predicted. Pathogenicity- and virulence-related annotations revealed 137 to 4558 genes associated with the TCDB, PHI, Cyt_P450, DFVF, effector, and secretome databases, respectively. A comparative analysis of isolate W-6, together with 51 other Colletotrichum strains, reveled 13 genes unique to the Orchidearum complex to which isolate W-6 belongs, highlighting the major facilitator superfamily transporters. The detailed analyses of MFS transporters associated with secondary metabolite gene clusters in isolate W-6 led to the identification and protein structure analyses of two key virulence factor candidates in DHA1 subclass, prlG and azaK, which were reported as efflux transporters of antibiotics in other pathogenic fungi. The assembly and further functional investigation of two pathogenic genes identified here potentially provide important resources for better understanding the biology and lifestyle of Colletotrichum and pave the way for designing more efficient strategies to control anthracnose in pecan plantations.

Characterization and Genome Analysis of Fusarium oxysporum Provides Insights into the Pathogenic Mechanisms of the Pokkah Boeng Disease in China

Pokkah Boeng Disease (PBD) is a severe and devastating disease that causes significant damage and yield losses in China. The pathogenic fungus Fusarium oxysporum is responsible for the rapid onset of top rot symptoms in sugarcane. In this study, we selected a representative strain, BS2-6, to perform morphological observations of colonies and determine pathogenicity. We examined the effects of BS2-6 infestation on the ultrastructure of sugarcane leaves. Moreover, we sequenced the whole genome of BS2-6 and examined the effects of various nitrogen sources and chemical reagents on its growth and pathogenicity. Our results indicate that sugarcane leaves inoculated with BS2-6 quickly succumb to heart leaf and growing rot. Ultrastructural analysis revealed that the surface tissues of the diseased leaves were destroyed with mycelium, and conidia blocked leaf stomata, which ultimately led to the degradation of leaf tissues. Ammoniacal nitrogen significantly promoted mycelial growth, pigment secretion, and the expression of genes related to secondary metabolite synthesis, thereby accelerating the development of PBD. In addition, we found that carbendazim effectively inhibited the growth of BS2-6 at various concentrations. These findings provide important insights for the effective prevention and control of PBD during sugarcane production.

Genomic Analysis of Penicillium griseofulvum CF3 Reveals Potential for Plant Growth Promotion and Disease Resistance

Penicillium griseofulvum CF3 is a fungus isolated from healthy strawberry soil, with the potential to promote the growth of plants and enhance their resistance to diseases. However, the genome sequence of P. griseofulvum CF3 remains unclear. Therefore, we performed the whole-genome CCS sequencing of P. griseofulvum CF3 using the PacBio Sequel II platform. The assembled genome comprised 104 contigs, with a total length of 37,564,657 bp, encoding 13,252 protein-coding genes. Comprehensive functional annotation was performed using various BLAST databases, including the non-redundant (Nr) protein sequence database, Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG), EuKaryotic Orthologous Groups (KOG), and the Carbohydrate-Active enZymes (CAZy) database, to identify and predict protein-coding genes, tRNAs, and rRNAs. The Antibiotics and Secondary Metabolites Analysis Shell (Antismash) analysis identified 50 biosynthetic gene clusters involved in secondary metabolite production within the P. griseofulvum CF3 genome. The whole-genome sequencing of P. griseofulvum CF3 helps us to understand its potential mechanisms in promoting plant growth and enhancing disease resistance, paving the way for the application of the CF3 strain in sustainable crop production.

Chromosome-level genome assembly of tetraploid Chinese cherry (Prunus pseudocerasus)

Chinese cherry belongs to the family Rosaceae, genus Prunus, and has high nutritional and economic value. 'Duiying' is a Chinese cherry variety local to Beijing, and has better performance than sweet cherry in terms of disease resistance. However, disease resistance resources of 'Duiying' have not been fully exploited partially due to the lack of a high-quality genome. In this study, we report a high-quality chromosome-scale genome assembly for Chinese cherry 'Duiying', by combining PacBio HiFi, Bionano and Hi-C sequencing data. The assembled genome has a size of 1035.19 Mb, with a scaffold N50 of 28.99 Mb, and 978.61 Mb (94.54%) assembled into 32 pseudochromosomes. A total of 547.16 Mb (52.86%) sequences were identified as repetitive sequences, and 114,451 protein-coding genes were annotated. Moreover, a total of 1635 microRNA (miRNA), 6637 transfer RNA (tRNA), 38,258 ribosomal RNA (rRNA), and 169 small nuclear RNAs (snRNA) genes were identified. The genome assembly presented here provides valuable genomic resources to enhance our understanding of genetic and molecular basis of Chinese cherry.

Multiomics analysis provides insights into musk secretion in muskrat and musk deer

BACKGROUND

Musk, secreted by the musk gland of adult male musk-secreting mammals, holds significant pharmaceutical and cosmetic potential. However, understanding the molecular mechanisms of musk secretion remains limited, largely due to the lack of comprehensive multiomics analyses and available platforms for relevant species, such as muskrat (Ondatra zibethicus Linnaeus) and Chinese forest musk deer (Moschus berezovskii Flerov).

RESULTS

We generated chromosome-level genome assemblies for the 2 species of muskrat (Ondatra zibethicus Linnaeus) and musk deer (Moschus berezovskii Flerov), along with 168 transcriptomes from various muskrat tissues. Comparative analysis with 11 other vertebrate genomes revealed genes and amino acid sites with signs of adaptive convergent evolution, primarily linked to lipid metabolism, cell cycle regulation, protein binding, and immunity. Single-cell RNA sequencing in muskrat musk glands identified increased acinar/glandular epithelial cells during secretion, highlighting the role of lipometabolism in gland development and evolution. Additionally, we developed MuskDB (http://muskdb.cn/home/), a freely accessible multiomics database platform for musk-secreting mammals.

CONCLUSIONS

The study concludes that the evolution of musk secretion in muskrats and musk deer is likely driven by lipid metabolism and cell specialization. This underscores the complexity of the musk gland and calls for further investigation into musk secretion-specific genetic variants.

A multiomics investigation into the evolution and specialized metabolisms of three Toxicodendron cultivars

Toxicodendron species are economically and medicinally important trees because of their rich sources of natural products. We present three chromosome-level genome assemblies of Toxicodendron vernicifluum 'Dali', Toxicodendron succedaneum 'Vietnam', and T. succedaneum 'Japan', which display diverse production capacities of specialized metabolites. Genome synteny and structural variation analyses revealed large genomic differences between the two species (T. vernicifluum and T. succedaneum) but fewer differences between the two cultivars within the species. Despite no occurrence of recent whole-genome duplications, Toxicodendron showed evidence of local duplications. The genomic modules with high expression of genes encoding metabolic flux regulators and chalcone synthase-like enzymes were identified via multiomics analyses, which may be responsible for the greater urushiol accumulation in T. vernicifluum 'Dali' than in other Toxicodendron species. In addition, our analyses revealed the regulatory patterns of lipid metabolism in T. succedaneum 'Japan', which differ from those of other Toxicodendron species and may contribute to its high lipid accumulation. Furthermore, we identified the key regulatory elements of lipid metabolism at each developmental stage, which could aid in molecular breeding to improve the production of urushiol and lipids in Toxicodendron species. In summary, this study provides new insights into the genomic underpinnings of the evolution and diversity of specialized metabolic pathways in three Toxicodendron cultivars through multiomics studies.

Chromosome-level genome assembly of banaba (Lagerstroemia speciosa L.)

Lagerstroemia speciosa is a famous medicinal and ornamental plant of the Lagerstroemia genus, with large and gorgeous flowers and a flowering period of 3 to 5 months. L. speciosa extracts have been used for many years in folk medicine to treat diabetes. Here, we used PacBio and Hi-C sequencing technologies to obtain a high-quality whole genome map of L. speciosa at chromosome level. The assembled genome is 306.76 Mb, with a scaffold N50 of 13.03 Mb. 98.75% of contigs were anchored to 24 pseudochromosomes, and 38.58% of the contigs (118.11 Mb) were identified as repeats. 93.54% of the 29,351 protein-coding genes were annotated. In addition, 146 miRNAs, 511 tRNAs, 2,733 rRNAs and 679 snRNAs were annotated. This high-quality genome assembly provides a valuable resource for understanding the species evolution in the Lythraceae family and promoting the study of important traits of L. speciosa.

The male and female genomes of golden pompano (Trachinotus ovatus) provide insights into the sex chromosome evolution and rapid growth

INTRODUCTION

Golden pompano (Trachinotus ovatus) is economically significant important for offshore cage aquaculture in China and Southeast Asian countries. Lack of high-quality genomic data and accurate gene annotations greatly restricts its genetic breeding progress.

OBJECTIVES

To decode the mechanisms of sex determination and rapid growth in golden pompano and facilitate the sex- and growth-aimed genetic breeding.

METHODS

Genome assemblies of male and female golden pompano were generated using Illumina, PacBio, BioNano, genetic maps and Hi-C sequencing data. Genomic comparisons, whole genome re-sequencing of 202 F1 individuals, QTL mapping and gonadal transcriptomes were used to analyze the sex determining region, sex chromosome evolution, SNP loci, and growth candidate genes. Zebrafish model was used to investigate the functions of growth candidate gene.

RESULTS

Female (644.45 Mb) and male (652.12 Mb) genomes of golden pompano were assembled and annotated at the chromosome level. Both genomes are highly conserved and no new or highly differentiated sex chromosomes occur. A 3.5 Mb sex determining region on LG15 was identified, where Hsd17b1, Micall2 and Lmx1a were putative candidates for sex determination. Three SNP loci significantly linked to growth were pinpointed, and a growth-linked gene gpsstr1 was identified by locus BSNP1369 (G → C, 17489695, Chr23). Loss of sstr1a (homologue of gpsstr1) in zebrafish caused growth retardation.

CONCLUSION

This study provides insights into sex chromosome evolution, sex determination and rapid growth of golden pompano.

Acquisition and evolution of the neurotoxin domoic acid biosynthesis gene cluster in Pseudo-nitzschia species

Of the hitherto over 60 taxonomically identified species in the genus of Pseudo-nitzschia, 26 have been confirmed to be toxigenic. Nevertheless, the acquisition and evolution of the toxin biosynthesis (dab) genes by this extensive group of Pseudo-nitzschia species remains unclear. Through constructing chromosome-level genomes of three Pseudo-nitzschia species and draft genomes of ten additional Pseudo-nitzschia species, putative genomic integration sites for the dab genes in Pseudo-nitzschia species were explored. A putative breakpoint was observed in syntenic regions in the dab gene cluster-lacking Pseudo-nitzschia species, suggesting potential independent losses of dab genes. The breakpoints between this pair of conserved genes were also identified in some dab genes-possessing Pseudo-nitzschia species, suggesting that the dab gene clusters transposed to other loci after the initial integration. A "single acquisition, multiple independent losses (SAMIL)" model is proposed to explain the acquisition and evolution of the dab gene cluster in Pseudo-nitzschia species.

Whole-Genome Sequencing of the Entomopathogenic Fungus Fusarium solani KMZW-1 and Its Efficacy Against Bactrocera dorsalis

Fusarium solani KMZW-1 is recognized for its potential as a biocontrol agent against agricultural and forestry pests, particularly due to its compatibility with integrated pest management (IPM) strategies. This study aimed to investigate the complete genome of F. solani KMZW-1 and assess its pathogenicity against Bactrocera dorsalis. Whole-genome sequencing revealed a genome size of 47,239,278 bp, comprising 27 contigs, with a GC content of 51.16% and fungus identified as F. solani KMZW-1. The genome completeness was assessed as 97.93% using BUSCO analysis, the DFVF sequence identifier was Fusarium 0G092560.1, and AntiSMASH analysis identified 35 gene clusters associated with secondary metabolite biosynthesis, providing insights into the genetic basis of its pathogenic mechanisms and biocontrol potential. Comparative genomic analysis found 269 unique genes for F. solani KMZW-1, and collinearity analysis exhibited a high degree of synteny with Fusarium solani-melongenae. The pathogenicity of F. solani KMZW-1 was assessed using concentrations ranging from 1 × 104 to 1 × 1011 conidia/mL. Higher concentrations (1 × 1010 to 1 × 1011 conidia/mL) resulted in significantly increased cumulative mortality rates of B. dorsalis adults compared to the control group. Notably, the pathogenicity was higher in male adults than in females. Probit analysis yielded LC50 (50% lethal concentration) values of 5.662 for female and 4.486 for male B. dorsalis adults. In summary, F. solani, KMZW-1 exhibits strong insecticidal activity against B. dorsalis and shows potential as a biocontrol agent with IPM strategies. These findings provide robust genomic evidence supporting the use of F. solani KMZW-1 in managing against B. dorsalis populations.

Comparative analysis of chromosome-level genomes provides insights into chromosomal evolution in Chiroptera

Chiroptera (bats) presents a fascinating model due to its remarkable variation in chromosome numbers, which range from 14 to 62. This astonishing diversity makes bats an excellent subject for studying chromosome evolution. The black-bearded tomb bat (Taphozous melanopogon) occupies a pivotal phylogenetic position within Chiroptera, emphasizing its crucial role in the systematic examination of bat chromosome evolution. In this study, we present the first chromosome-level genome of T. melanopogon within the family Emballonuridae. Together with previously published genomes, we construct a strongly supported phylogenetic tree of bats, which supports that Emballonuridae forms a basal group within Yangochiroptera. Furthermore, we reconstruct ancestral karyotypes at key nodes along the bat phylogeny and conduct a synteny analysis among the genomes of 12 bat species. Our findings identified evolutionary breakpoint regions (EBRs) that are of particular interest. Notably, some bat genomes exhibit an enrichment of genes related to host defense against microbial pathogens within EBRs. Remarkably, one species possesses multiple copies of some β-defensin genes, while six other species have experienced the loss of some β-defensin genes due to EBRs. Furthermore, some olfactory receptor genes are located in EBRs of 12 species, 4 of which have a significant enrichment in sensory perception of smell. Together, our comparative genomic analysis underscores the potential link between chromosome rearrangements and the adaptation of bats to defend against microbial pathogens.

Chromosome-level genome assembly of the cashmere goat

The goat, an early domesticated ruminant, is a reliable source of cashmere, meat and milk in global agricultural production. Despite this, the genome of cashmere-rich goats has yet to be characterized. Here, we assembled the nearly complete genome of a cashmere goat from a highly economically valuable Inner Mongolian Cashmere buck, utilizing a combination of PacBio HiFi, ONT ultra-long reads, and Hi-C technologies. The size of this genome is 2.76 Gb, with a contig N50 of 95.22 Mb. All assembled sequences were anchored onto 29 autosomes and both sex chromosomes, with only two gaps present on the X chromosome. We identified 1,333.29 Mb (48.26%) of repetitive sequences and predicted 22,480 protein-coding genes. Assembly quality assessment of the genome demonstrated that our assembled cashmere goat genome surpasses the continuity, completeness, and accuracy of other published goat genomes. Taken together, we provided the first cashmere goat assembly, bridging the gap in the genome of important economic breeds of domestic goats, and providing a valuable reference resource for goat genetics and genome research.

Two high quality chromosome-scale genome assemblies of female and male silver pomfret (Pampus argenteus)

Pampus argenteus is a highly commercial marine fish whose population is declining sharply. Here, we generated a female P. argenteus genome, spanning 536.33 Mb with contig N50 of 1.79 Mb; 24070 genes (99.50% of 24,182) were functionally annotated. To improve quality of it, we assembled a 553.79 Mb genome of male fish with contig N50 of 24.75 Mb through HiFi and ultra-long ONT sequence technologies; 550.82 Mb were anchored onto 24 gap-free chromosomes; 22,892 genes (98.1% of 23,346) were functionally annotated; the QV value was 51.55 with 98.9% of BUSCO and 99.39% coverage of Illumina reads. Finally, we compared this genome with previous published one, revealing 37,301 SVs. 52.82 Mb and 18.05 Mb SDs were characterized in our and published assemblies, respectively, and 48.96 Mb PURs were constructed. Thus, this genome assembly exhibits excellent completeness, continuity and accuracy comparing to the published one, which can be current preferred reference genome. Overall, these works help aquaculture and wild resources recovery of P. argenteus and provide a valuable genetic resource for study.

Transposable elements in Drosophila montana from harsh cold environments

BACKGROUND

Substantial discoveries during the past century have revealed that transposable elements (TEs) can play a crucial role in genome evolution by affecting gene expression and inducing genetic rearrangements, among other molecular and structural effects. Yet, our knowledge on the role of TEs in adaptation to extreme climates is still at its infancy. The availability of long-read sequencing has opened up the possibility to identify and study potential functional effects of TEs with higher precision. In this work, we used Drosophila montana as a model for cold-adapted organisms to study the association between TEs and adaptation to harsh climates.

RESULTS

Using the PacBio long-read sequencing technique, we de novo identified and manually curated TE sequences in five Drosophila montana genomes from eco-geographically distinct populations. We identified 489 new TE consensus sequences which represented 92% of the total TE consensus in D. montana. Overall, 11-13% of the D. montana genome is occupied by TEs, which as expected are non-randomly distributed across the genome. We identified five potentially active TE families, most of them from the retrotransposon class of TEs. Additionally, we found TEs present in the five analyzed genomes that were located nearby previously identified cold tolerant genes. Some of these TEs contain promoter elements and transcription binding sites. Finally, we detected TEs nearby fixed and polymorphic inversion breakpoints.

CONCLUSIONS

Our research revealed a significant number of newly identified TE consensus sequences in the genome of D. montana, suggesting that non-model species should be studied to get a comprehensive view of the TE repertoire in Drosophila species and beyond. Genome annotations with the new D. montana library allowed us to identify TEs located nearby cold tolerant genes, and present at high population frequencies, that contain regulatory regions and are thus good candidates to play a role in D. montana cold stress response. Finally, our annotations also allow us to identify for the first time TEs present in the breakpoints of three D. montana inversions.

High-quality genome of a modern soybean cultivar and resequencing of 547 accessions provide insights into the role of structural variation

Soybean provides protein, oil and multiple health-related compounds. Understanding the effects of structural variations (SVs) on economic traits in modern breeding is important for soybean improvement. Here we assembled the high-quality genome of modern cultivar Nongdadou2 (NDD2) and identified 25,814 SV-gene pairs compared to 29 reported genomes, with 13 NDD2-private SVs validated in 547 deep-resequencing (average = 18.05-fold) accessions, which advances our understanding of genomic variation biology. We found some insertions/deletions involved in seed protein and weight formation, an inversion related to adaptation to drought and a large intertranslocation implicated in a key divergence event in soybean. Of 749,714 SVs from 547 accessions, 6,013 were significantly associated with 22 yield-related and seed-quality-related traits determined in ten location × year environments. We uncovered 1,761 associated SVs that hit genes or regulatory regions, with 12 in GmMQT influencing oil and isoflavone contents. Our work provides resources and insights into SV roles in soybean improvement.

High-quality chromosome-level genomic insights into molecular adaptation to low-temperature stress in Madhuca longifolia in southern subtropical China

BACKGROUND

Madhuca longifolia, the energy-producing and medicinal tropical tree originally from southern India, faces difficulties in adapting to the low temperatures of late autumn and early winter in subtropical southern China, impacting its usability. Therefore, understanding the molecular mechanisms controlling the ability of this species to adapt to environmental challenges is essential for optimising horticulture efforts. Accordingly, this study aimed to elucidate the molecular responses of M. longifolia to low-temperature stress through genomic and transcriptomic analyses to inform strategies for its effective cultivation and utilisation in colder climates.

RESULTS

Herein, the high-quality reference genome and genomic assembly for M. longifolia are presented for the first time. Using Illumina sequencing, Hi-C technology, and PacBio HiFi sequencing, we assembled a chromosome-level genome approximately 737.92 Mb in size, investigated its genomic features, and conducted an evolutionary analysis of the genus Madhuca. Additionally, using transcriptome sequencing, we identified 17,941 differentially expressed genes related to low-temperature response. Through bioinformatics analysis of the WRKY gene family, 15 genes crucial for M. longifolia low-temperature resistance were identified.

CONCLUSIONS

This research not only lays the groundwork for the successful ecological adaptation and cultivation of M. longifolia in China's southern subtropical regions but also offers valuable insights for the genetic enhancement of cold tolerance in tropical species, contributing to their sustainable horticulture and broader industrial, medicinal, and agricultural use.

Genome assembly of Stephania longa provides insight into cepharanthine biosynthesis

Introduction

Stephania longa, a medicinal plant renowned for producing cepharanthine, has gained significance due to the compound's notable antiviral properties against SARS-CoV-2. However, a comprehensive genetic understanding of S. longa has been lacking. This study aimed to develop a high-quality, chromosome-level genome assembly to uncover the genetic intricacies and evolutionary narrative of this species. By integrating genomic data with metabolomic and transcriptomic analyses, we sought to identify key genes involved in cepharanthine biosynthesis.

Methods

We employed a multi-faceted approach comprising genome assembly, phylogenetic analysis, gene family dynamics investigation, metabolomic profiling, and gene expression analysis across various tissues of S. longa. This integrated strategy enabled the identification of key genes involved in cepharanthine biosynthesis and elucidated the species' evolutionary history.

Results

Our phylogenetic analysis clarified the placement of the genus Stephania within the Ranunculales order and revealed its notably high mutation rate. We identified gene family expansions and signs of positive selection likely contributing to Stephania's unique metabolic capabilities. Metabolomic profiling uncovered complex regulatory mechanisms orchestrating the biosynthesis and distribution of cepharanthine and related metabolites. Through the integration of genomic, transcriptomic, and metabolomic data, we identified genes with expression patterns and evolutionary trajectories suggesting pivotal roles in cepharanthine biosynthesis, including those involved in crucial biosynthetic steps.

Discussion

This comprehensive study, integrating genomic, metabolomic, and transcriptomic approaches, provides valuable insights into S. longa's biosynthetic potential. It not only enhances our understanding of the species but also establishes a foundation for future investigations into the biosynthesis and therapeutic exploitation of cepharanthine and related alkaloids.

A Genomic Sequence Resource of Diaporthe mahothocarpus GZU-Y2 Causing Leaf Spot Blight in Camellia oleifera

Diaporthe mahothocarpus GZU-Y2, a new pathogen responsible for leaf spot blight disease, leads to significant damage and economic losses in some Camellia oleifera plantations. The current study annotated the genome of the D. mahothocarpus isolate GZU-Y2 to advance our knowledge of the pathogen and facilitate improved disease management of leaf spot blight. The initial Pacbio-Illumina hybrid draft genome for GZU-Y2 resulted in a high-quality assembly with 62 contigs, characterized by an N50 length of 7.07 Mb. The complete genome of isolate GZU-Y2 was 58.97 Mbp, with a GC content of 50.65%. Importantly, the assembly exhibits remarkable integrity, with 97.93% of complete BUSCO validating genome completeness. The prediction results showed that a total of 15,918 protein-coding genes were annotated using multiple bioinformatics databases. The genome assembly and annotation resource reported here will be useful for the further study of fungal infection mechanisms and pathogen-host interaction.

New insights into decoding the lifestyle of endophytic Fusarium lateritium Fl617 via comparing genomes

Fungal-plant interactions have persisted for 460 million years, and almost all terrestrial plants on Earth have endophytic fungi. However, the mechanism of symbiosis between endophytic fungi and host plants has been inconclusive. In this dissertation, we used a strain of endophytic Fusarium lateritium (Fl617), which was found in the previous stage to promote disease resistance in tomato, and selected the pathogenic Fusarium oxysporum Fo4287 and endophytic Fusarium oxysporum Fo47, which are in the same host and the closest relatives of Fl617, to carry out a comparative genomics analysis of the three systems and to provide a new perspective for the elucidation of the special lifestyle of the fungal endophytes. We found that endophytic F. lateritium has a smaller genome, fewer clusters and genes associated with pathogenicity, and fewer plant cell wall degrading enzymes (PCWDEs). There were also relatively fewer secondary metabolisms and typical Fusarium spp. toxins, and a lack of the key Fusarium spp. pathogenicity factor, secreted in xylem (SIX), but the endophytic fungi may be more sophisticated in their regulation of the colonization process. It is hypothesized that the endophytic fungi may have maintained their symbiosis with plants due to the relatively homogeneous microenvironment in plants for a long period of time, considering only plant interactions and discarding the relevant pathogenicity factors, and that their endophytic evolutionary tendency may tend to be genome streamlining and to enhance the fineness of the regulation of plant interactions, thus maintaining their symbiotic status with plants.

Chromosome-level genome sequencing and multi-omics of the Hungarian White Goose (Anser anser domesticus) reveals novel miRNA-mRNA regulation mechanism of waterfowl feather follicle development

The Hungarian White Goose (Anser anser domesticus) is an excellent European goose breed, with high feather and meat production. Despite its importance in the poultry industry, no available genome assembly information has been published. This study aimed to present Chromosome-level and functional genome sequencing of the Hungarian White Goose. The results showed that the genome assembly has a total length of 1115.82 Mb, 39 pairs of chromosomes, 92.98% of the BUSCO index, and contig N50 and scaffold N50 were up to 2.32 Mb and 60.69 Mb, respectively. Annotation of the genome assembly revealed 19550 genes, 286 miRNAs, etc. We identified 235 expanded and 1,167 contracted gene families in this breed compared with the other 16 species. We performed a positive selection analysis between this breed and four species of Anatidae to uncover the genetic information underlying feather follicle development. Further, we detected the function of miR-199-x, miR-143-y, and miR-23-z on goose embryonic skin fibroblast. In summary, we have successfully generated a highly complete genome sequence of the Hungarian white goose, which will provide a great resource to improve our understanding of gene functions and enhance the studies on feather follicle development at the genomic level.

Single-cell mitochondrial sequencing reveals low-frequency mitochondrial mutations in naturally aging mice

Mitochondria play a crucial role in numerous biological processes; however, limited methods and research have focused on revealing mitochondrial heterogeneity at the single-cell level. In this study, we optimized the DNBelab C4 single-cell ATAC (assay for transposase-accessible chromatin) sequencing workflow for single-cell mitochondrial sequencing (C4_mtscATAC-seq). We validated the effectiveness of our C4_mtscATAC-seq protocol by sequencing the HEK-293T cell line with two biological replicates, successfully capturing both mitochondrial content (~68% of total sequencing data) and open chromatin status simultaneously. Subsequently, we applied C4_mtscATAC-seq to investigate two mouse tissues, spleen and bone marrow, obtained from two mice aged 2 months and two mice aged 23 months. Our findings revealed higher mitochondrial DNA (mtDNA) content in young tissues compared to more variable mitochondrial content in aged tissues, consistent with higher activity scores of nuclear genes associated with mitochondrial replication and transcription in young tissues. We detected a total of 22, 15, and 21 mtDNA mutations in the young spleen, aged spleen, and bone marrow, respectively, with most variant allele frequencies (VAF) below 1%. Moreover, we observed a higher number of mtDNA mutations with higher VAF in aged tissues compared to young tissues. Importantly, we identified three mtDNA variations (m.9821A>T, m.15219T>C, and m.15984C>T) with the highest VAF in both aged spleen and aged bone marrow. By comparing cells with and without these mtDNA variations, we analyzed differential open chromatin status to identify potential genes associated with these mtDNA variations, including transcription factors such as KLF15 and NRF1. Our study presents an alternative single-cell mitochondrial sequencing method and provides crude insights into age-related single-cell mitochondrial variations.

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

Genomic and single-cell analyses reveal genetic signatures of swimming pattern and diapause strategy in jellyfish

Jellyfish exhibit innovative swimming patterns that contribute to exploring the origins of animal locomotion. However, the genetic and cellular basis of these patterns remains unclear. Herein, we generated chromosome-level genome assemblies of two jellyfish species, Turritopsis rubra and Aurelia coerulea, which exhibit straight and free-swimming patterns, respectively. We observe positive selection of numerous genes involved in statolith formation, hair cell ciliogenesis, ciliary motility, and motor neuron function. The lineage-specific absence of otolith morphogenesis- and ciliary movement-related genes in T. rubra may be associated with homeostatic structural statocyst loss and straight swimming pattern. Notably, single-cell transcriptomic analyses covering key developmental stages reveal the enrichment of diapause-related genes in the cyst during reverse development, suggesting that the sustained diapause state favours the development of new polyps under favourable conditions. This study highlights the complex relationship between genetics, locomotion patterns and survival strategies in jellyfish, thereby providing valuable insights into the evolutionary lineages of movement and adaptation in the animal kingdom.

HiTE: a fast and accurate dynamic boundary adjustment approach for full-length transposable element detection and annotation

Recent advancements in genome assembly have greatly improved the prospects for comprehensive annotation of Transposable Elements (TEs). However, existing methods for TE annotation using genome assemblies suffer from limited accuracy and robustness, requiring extensive manual editing. In addition, the currently available gold-standard TE databases are not comprehensive, even for extensively studied species, highlighting the critical need for an automated TE detection method to supplement existing repositories. In this study, we introduce HiTE, a fast and accurate dynamic boundary adjustment approach designed to detect full-length TEs. The experimental results demonstrate that HiTE outperforms RepeatModeler2, the state-of-the-art tool, across various species. Furthermore, HiTE has identified numerous novel transposons with well-defined structures containing protein-coding domains, some of which are directly inserted within crucial genes, leading to direct alterations in gene expression. A Nextflow version of HiTE is also available, with enhanced parallelism, reproducibility, and portability.

Chromosome-level genome of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus

Recent conservation efforts to protect rare and endangered aquatic species have intensified. Nevertheless, the ornate spiny lobster (Panulirus ornatus), which is prevalent in the Indo-Pacific waters, has been largely ignored. In the absence of a detailed genomic reference, the conservation and population genetics of this crustacean are poorly understood. Here, We assembled a comprehensive chromosome-level genome for P. ornatus. This genome-among the most detailed for lobsters-spans 2.65 Gb with a contig N50 of 51.05 Mb, and 99.11% of the sequences with incorporated to 73 chromosomes. The ornate spiny lobster genome comprises 65.67% repeat sequences and 22,752 protein-coding genes with 99.20% of the genes functionally annotated. The assembly of the P. ornatus genome provides valuable insights into comparative crustacean genomics and endangered species conservation, and lays the groundwork for future research on the speciation, ecology, and evolution of the ornate spiny lobster.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Chromosome-level Alstonia scholaris genome unveils evolutionary insights into biosynthesis of monoterpenoid indole alkaloids

Alstonia scholaris of the Apocynaceae family is a medicinal plant with a rich source of bioactive monoterpenoid indole alkaloids (MIAs), which possess anti-cancer activity like vinca alkaloids. To gain genomic insights into MIA biosynthesis, we assembled a high-quality chromosome-level genome for A. scholaris using nanopore and Hi-C data. The 444.95 Mb genome contained 35,488 protein-coding genes. A total of 20 chromosomes were assembled with a scaffold N50 of 21.75 Mb. The genome contained a cluster of strictosidine synthases and tryptophan decarboxylases with synteny to other species and a saccharide-terpene cluster involved in the monoterpenoid biosynthesis pathway of the MIA upstream pathway. The multi-omics data of A. scholaris provide a valuable resource for understanding the evolutionary origins of MIAs and for discovering biosynthetic pathways and synthetic biology efforts for producing pharmaceutically useful alkaloids.

Genome-Wide Identification and Characterization of Homeobox Transcription Factors in Phoma sorghina var. saccharum Causing Sugarcane Twisted Leaf Disease

A homeobox transcription factor is a conserved transcription factor, ubiquitous in eukaryotes, that regulates the tissue formation of structure, cell differentiation, proliferation, and cancer. This study identified the homeobox transcription factor family and its distribution in Phoma sorghina var. saccharum at the whole genome level. It elucidated the gene structures and evolutionary characteristics of this family. Additionally, knockout experiments were carried out and the preliminary function of these transcription factors was studied. Through bioinformatics approaches, nine homeobox transcription factors (PsHOX1-PsHOX9) were identified in P. sorghina var. saccharum, and these contained HOX-conserved domains and helix-turn-helix secondary structures. Nine homeobox gene deletion mutants were obtained using the homologous recombinant gene knockout technique. Protoplast transformation was mediated by polyethylene glycol (PEG) and the transformants were identified using PCR. The knockouts of PsHOX1, PsHOX2, PsHOX3, PsHOX4, PsHOX6, PsHOX8, and PsHOX9 genes resulted in a smaller growth diameter in P. sorghina var. saccharum. In contrast, the knockouts of the PsHOX3, PsHOX6, and PsHOX9 genes inhibited the formation of conidia and led to a significant decrease in the pathogenicity. This study's results will provide insights for understanding the growth and development of P. sorghina var. saccharum. The pathogenic mechanism of the affected sugarcane will provide an essential theoretical basis for preventing and controlling sugarcane twisted leaf disease.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

Chromosome-scale genome assembly and annotation of Cotoneaster glaucophyllus

Cotoneaster glaucophyllus is a semi-evergreen plant that blossoms in late summer, producing dense, attractive, fragrant white flowers with significant ornamental and ecological value. Here, a chromosome-scale genome assembly was obtained by integrating PacBio and Illumina sequencing data with the aid of Hi-C technology. The genome assembly was 563.3 Mb in length, with contig N50 and scaffold N50 values of ~6 Mb and ~31 Mb, respectively. Most (95.59%) of the sequences were anchored onto 17 pseudochromosomes (538.4 Mb). We predicted 35,856 protein-coding genes, 1,401 miRNAs, 655 tRNAs, 425 rRNAs, and 795 snRNAs. The functions of 34,967 genes (97.52%) were predicted. The availability of this chromosome-level genome will provide valuable resources for molecular studies of this species, facilitating future research on speciation, functional genomics, and comparative genomics within the Rosaceae family.

Chromosome-level genome assembly of marine diatom Skeletonema tropicum

Skeletonema tropicum is a marine diatom of the genus Skeletonema that also includes many well-known species including S. marinoi. S. tropicum is a high temperature preferring species thriving in tropical ocean regions or temperate ocean regions during summer-autumn. However, mechanisms of ecological adaptation of S. tropicum remain poorly understood due partially to the lack of a high-quality whole genome assembly. Here, we report the first high-quality chromosome-scale genome assembly for S. tropicum, using cutting-edge technologies including PacBio single molecular sequencing and high-throughput chromatin conformation capture. The assembled genome has a size of 78.78 Mb with a scaffold N50 of 3.17 Mb, anchored to 23 pseudo-chromosomes. In total, 20,613 protein-coding genes were predicted, of which 17,757 (86.14%) genes were functionally annotated. Collinearity analysis of the genomes of S. tropicum and S. marinoi revealed that these two genomes were highly homologous. This chromosome-level genome assembly of S. tropicum provides a valuable genomic platform for comparative analysis of mechanisms of ecological adaption.

Genome assembly of Melilotus officinalis provides a new reference genome for functional genomics

BACKGROUND

Sweet yellow clover (Melilotus officinalis) is a diploid plant (2n = 16) that is native to Europe. It is an excellent legume forage. It can both fix nitrogen and serve as a medicine. A genome assembly of Melilotus officinalis that was collected from Best corporation in Beijing is available based on Nanopore sequencing. The genome of Melilotus officinalis was sequenced, assembled, and annotated.

RESULTS

The latest PacBio third generation HiFi assembly and sequencing strategies were used to produce a Melilotus officinalis genome assembly size of 1,066 Mbp, contig N50 = 5 Mbp, scaffold N50 = 130 Mbp, and complete benchmarking universal single-copy orthologs (BUSCOs) = 96.4%. This annotation produced 47,873 high-confidence gene models, which will substantially aid in our research on molecular breeding. A collinear analysis showed that Melilotus officinalis and Medicago truncatula shared conserved synteny. The expansion and contraction of gene families showed that Melilotus officinalis expanded by 565 gene families and shrank by 56 gene families. The contacted gene families were associated with response to stimulus, nucleotide binding, and small molecule binding. Thus, it is related to a family of genes associated with peptidase activity, which could lead to better stress tolerance in plants.

CONCLUSIONS

In this study, the latest PacBio technology was used to assemble and sequence the genome of the Melilotus officinalis and annotate its protein-coding genes. These results will expand the genomic resources available for Melilotus officinalis and should assist in subsequent research on sweet yellow clover plants.

The genome of Stephania japonica provides insights into the biosynthesis of cepharanthine

Stephania japonica is an early-diverging eudicotyledon plant with high levels of cepharanthine, proven to be effective in curing coronavirus infections. Here, we report a high-quality S. japonica genome. The genome size is 688.52 Mb, and 97.37% sequences anchor to 11 chromosomes. The genome comprises 67.46% repetitive sequences and 21,036 genes. It is closely related to two Ranunculaceae species, which diverged from their common ancestor 55.90-71.02 million years ago (Mya) with a whole-genome duplication 85.59-96.75 Mya. We further reconstruct ancestral karyotype of Ranunculales. Several cepharanthine biosynthesis genes are identified and verified by western blot. Two genes (Sja03G0243 and Sja03G0241) exhibit catalytic activity as shown by liquid chromatography-mass spectrometry. Then, cepharanthine biosynthesis genes, transcription factors, and CYP450 family genes are used to construct a comprehensive network. Finally, we construct an early-diverging eudicotyledonous genome resources (EEGR) database. As the first genome of the Menispermaceae family to be released, this study provides rich resources for genomic studies.

Chromosome-level genome assembly of Platycarya strobilacea

Platycarya strobilacea belongs to the walnut family (Juglandaceae), is commonly known as species endemic to East Asia, and is an ecologically important, wind pollinated, woody deciduous tree. To facilitate this ancient tree for the ecological value and conservation of this ancient tree, we report a new high-quality genome assembly of P. strobilacea. The genome size was 677.30 Mb, with a scaffold N50 size of 45,791,698 bp, and 98.43% of the assembly was anchored to 15 chromosomes. We annotated 32,246 protein-coding genes in the genome, of which 96.30% were functionally annotated in six databases. This new high-quality assembly of P. strobilacea provide valuable resource for the phylogenetic and evolutionary analysis of the walnut family and angiosperm.

Genomics insights into flowering and floral pattern formation: regional duplication and seasonal pattern of gene expression in Camellia

BACKGROUND

The formation and domestication of ornamental traits are influenced by various aspects, such as the recognition of esthetic values and cultural traditions. Camellia japonica is widely appreciated and domesticated around the world mainly due to its rich variations in ornamental traits. Ornamental camellias have a diverse range of resources, including different bud variations from Camellia spp. as well as inter- and intra- specific hybridization. Despite research on the formation of ornamental traits, a basic understanding of their genetics and genomics is still lacking.

RESULTS

Here, we report the chromosomal-level reference genome of C. japonica through combining multiple DNA-sequencing technologies and obtain a high-density genetic linkage map of 4255 markers by sequencing 98 interspecific F1 hybrids between C. japonica and C. chekiangoleosa. We identify two whole-genome duplication events in C. japonica: one is a shared ancient γ event, and the other is revealed to be specific to genus Camellia. Based on the micro-collinearity analysis, we find large-scale segmental duplication of chromosome 8, resulting to two copies of the AGAMOUS loci, which may play a key role in the domestication of floral shapes. To explore the regulatory mechanisms of seasonal flowering, we have analyzed year-round gene expression patterns of C. japonica and C. azalea-a sister plant of continuous flowering that has been widely used for cross breeding. Through comparative analyses of gene co-expression networks and annual gene expression patterns, we show that annual expression rhythms of some important regulators of seasonal growth and development, including GIGANTEA and CONSTANS of the photoperiod pathway, have been disrupted in C. azalea. Furthermore, we reveal that the distinctive expression patterns of FLOWERING LOCUS T can be correlated with the seasonal activities of flowering and flushing. We demonstrate that the regulatory module involved in GIGANTEA, CONSTANS, and FLOWERING LOCUS T is central to achieve seasonality.

CONCLUSIONS

Through the genomic and comparative genomics characterizations of ornamental Camellia spp., we propose that duplication of chromosomal segments as well as the establishment of gene expression patterns has played a key role in the formation of ornamental traits (e.g., flower shape, flowering time). This work provides a valuable genomic platform for understanding the molecular basis of ornamental traits.

Genome-wide transcriptome profiling reveals molecular response pathways of Trichoderma harzianum in response to salt stress

Trichoderma harzianum exhibits a strong biological control effect on many important plant pathogens, such as Fusarium oxysporum, Botrytis cinerea, and Meloidogyne. However, its biocontrol effectiveness is weakened or reduced under salt stress. The aim of this study was to investigate the molecular response of T. harzianum to salt stress at the whole-genome level. Here, we present a 44.47 Mb near-complete genome assembly of the T. harzianum qt40003 strain for the first time, which was assembled de novo with 7.59 Gb Nanopore sequencing long reads (~170-fold) and 5.2 Gb Illumina short reads (~116-fold). The assembled qt40003 genome contains 12 contigs, with a contig N50 of 4.81 Mb, in which four of the 12 contigs were entirely reconstructed in a single chromosome from telomere to telomere. The qt40003 genome contains 4.27 Mb of repeat sequences and 12,238 protein-coding genes with a BUSCO completeness of 97.5%, indicating the high accuracy and completeness of our gene annotations. Genome-wide transcriptomic analysis was used to investigate gene expression changes related to salt stress in qt40003 at 0, 2% (T2), and 4% (T4) sodium chloride concentrations. A total of 2,937 and 3,527 differentially expressed genes (DEGs) were obtained under T2 and T4 conditions, respectively. GO enrichment analysis showed that the T2-treatment DEGs were highly enriched in detoxification (p < 0.001), while the T4 DEGs were mainly enriched in cell components, mostly in cellular detoxification, cell surface, and cell wall. KEGG metabolic pathway analysis showed that 91 and 173 DEGs were significantly enriched in the T2 and T4 treatments, respectively (p < 0.01), mainly in the glutathione metabolism pathway. We further experimentally analyzed the differentially expressed glutathione transferase genes in the glutathione metabolic pathway, most of which were downregulated (13/15). In addition, we screened 13 genes related to active oxygen clearance, including six upregulated and seven downregulated genes, alongside five fungal hydrophobic proteins, of which two genes were highly expressed. Our study provides high-quality genome information for the use of T. harzianum for biological control and offers significant insights into the molecular responses of T. harzianum under salt-stress conditions.

Chromosome-level genome assembly of the Asian citrus psyllid, Diaphorina citri

Diaphorina citri is a global citrus pest. As a vector insect, it can transmit the causative agents of citrus huanglongbing, causing irreversible losses to the citrus industry. The acquisition of genomic information can provide a molecular genetic basis for effective control of D. citri. Here, the DNBSEQ™ , Oxford Nanopore Technologies, and Hi-C technologies are applied to generate a high-quality chromosome-level genome of D. citri. The genome size of D. citri was 523.78 Mb with a scaffold N50 of 47.05 Mb distributed on 13 chromosomes. A total of 250.64 Mb (47.85%) repeat sequences and 24 048 protein-coding genes were predicted. Genome resequencing of female and male individuals indicated that the sex chromosome system of D. citri is XO. Phylogenetic analysis demonstrated that D. citri and Pachypsylla venusta, which separated from their most recent common ancestor about 336.62 million years ago, were the most closely related. Additionally, we identified genes potentially involved in detoxification metabolism, pathogen transmission, and honeydew secretion for further investigation. The high-quality genome provides an important reference for developing effective management strategies of D. citri.

Near telomere-to-telomere genome of the model plant Physcomitrium patens

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Genome Sequencing and Analysis of Nigrospora oryzae, a Rice Leaf Disease Fungus

Nigrospora oryzae is one of several fungal pathogens known to cause brown streaks, leaf spots, and latent infections in rice. In this study, the entire 42.09-Mb genome of N. oryzae was sequenced at a depth of 169× using the Oxford Nanopore Technologies platform. The draft genome sequence was comprised of 26 scaffolds, possessed an average GC content of 58.83%, and contained a total of 10,688 protein-coding genes. Analysis of the complete genome sequence revealed that CAZyme-encoding genes account for 6.11% of all identified genes and that numerous transcription factors (TFs) associated with diverse biological processes belong predominantly to Zn-clus (22.20%) and C2H2 (10.59%) fungal TF classes. In addition, genes encoding 126 transport proteins and 3307 pathogen-host interaction proteins were identified. Comparative analysis of the previously reported N. oryzae reference strain GZL1 genome and the genome of a representative strain ZQ1 obtained here revealed 9722 colinear genes. Collectively, these findings provide valuable insights into N. oryzae genetic mechanisms and phenotypic characteristics.

Improved chromosome-level genome assembly of Indian sandalwood (Santalum album)

Santalum album is a well-known aromatic and medicinal plant that is highly valued for the essential oil (EO) extracted from its heartwood. In this study, we present a high-quality chromosome-level genome assembly of S. album after integrating PacBio Sequel, Illumina HiSeq paired-end and high-throughput chromosome conformation capture sequencing technologies. The assembled genome size is 207.39 M with a contig N50 of 7.33 M and scaffold N50 size of 18.31 M. Compared with three previously published sandalwood genomes, the N50 length of the genome assembly was longer. In total, 94.26% of the assembly was assigned to 10 pseudo-chromosomes, and the anchor rate far exceeded that of a recently released value. BUSCO analysis yielded a completeness score of 94.91%. In addition, we predicted 23,283 protein-coding genes, 89.68% of which were functionally annotated. This high-quality genome will provide a foundation for sandalwood functional genomics studies, and also for elucidating the genetic basis of EO biosynthesis in S. album.

The draft genome of the Temminck's tragopan (Tragopan temminckii) with evolutionary implications

BACKGROUND

High-quality genome data of birds play a significant role in the systematic study of their origin and adaptive evolution. The Temminck's tragopan (Tragopan temminckii) (Galliformes, Phasianidae), a larger pheasant, is one of the most abundant and widely distributed species of the genus Tragopan, and was defined as class II of the list of national key protected wild animals in China. The absence of a sequenced genome has restricted previous evolutionary trait studies of this taxa.

RESULTS

The whole genome of the Temminck's tragopan was sequenced using Illumina and PacBio platform, and then de novo assembled and annotated. The genome size was 1.06 Gb, with a contig N50 of 4.17 Mb. A total of 117.22 Mb (11.00%) repeat sequences were identified. 16,414 genes were predicted using three methods, with 16,099 (98.08%) annotated as functional genes based on five databases. In addition, comparative genome analyses were conducted across 12 Galliformes species. The results indicated that T. temminckii was the first species to branch off from the clade containing Lophura nycthemera, Phasianus colchicus, Chrysolophus pictus, Syrmaticus mikado, Perdix hodgsoniae, and Meleagris gallopavo, with a corresponding divergence time of 31.43 million years ago (MYA). Expanded gene families associated with immune response and energy metabolism were identified. Genes and pathways associated with plumage color and feather development, immune response, and energy metabolism were found in the list of positively selected genes (PSGs).

CONCLUSIONS

A genome draft of the Temminck's tragopan was reported, genome feature and comparative genome analysis were described, and genes and pathways related to plumage color and feather development, immune response, and energy metabolism were identified. The genomic data of the Temminck's tragopan considerably contribute to the genome evolution and phylogeny of the genus Tragopan and the whole Galliformes species underlying ecological adaptation strategies.

Genome and haplotype provide insights into the population differentiation and breeding improvement of Gossypium barbadense

INTRODUCTION

Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare.

OBJECTIVES

We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton.

METHODS

Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3-79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton.

RESULTS

One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian.

CONCLUSION

Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.