Two haplotype-resolved, gap-free genome assemblies for Actinidia latifolia and Actinidia chinensis shed light on the regulatory mechanisms of vitamin C and sucrose metabolism in kiwifruit

The telomere-to-telomere genome of Pucai () (Typha angustifolia L.): a distinctive semiaquatic vegetable with lignin and chlorophyll as quality characteristics

Pucai () (Typha angustifolia L.), within the Typha spp., is a distinctive semiaquatic vegetable. Lignin and chlorophyll are two crucial traits and quality indicators for Pucai. In this study, we assembled a 207.00-Mb high-quality gapless genome of Pucai, telomere-to-telomere (T2T) level with a contig N50 length of 13.73 Mb. The most abundant type of repetitive sequence, comprising 16.98% of the genome, is the long terminal repeat retrotransposons (LTR-RT). A total of 30 telomeres and 15 centromeric regions were predicted. Gene families related to lignin, chlorophyll biosynthesis, and disease resistance were greatly expanded, which played important roles in the adaptation of Pucai to wetlands. The slow evolution of Pucai was indicated by the σ whole-genome duplication (WGD)-associated Ks peaks from different Poales and the low activity of recent LTR-RT in Pucai. Meanwhile, we found a unique WGD event in Typhaceae. A statistical analysis and annotation of genomic variations were conducted in interspecies and intraspecies of Typha. Based on the T2T genome, we constructed lignin and chlorophyll metabolic pathways of Pucai. Subsequently, the candidate structural genes and transcription factors that regulate lignin and chlorophyll biosynthesis were identified. The T2T genomic resources will provide molecular information for lignin and chlorophyll accumulation and help to understand genome evolution in Pucai.

Construction of the super pan-genome for the genus Actinidia reveals structural variations linked to phenotypic diversity

Kiwifruits, belonging to the genus Actinidia, are acknowledged as one of the most successfully domesticated fruits in the twentieth century. Despite the rich wild resources and diverse phenotypes within this genus, insights into the genomic changes are still limited. Here, we conducted whole-genome sequencing on seven representative materials from highly diversified sections of Actinidia, leading to the assembly and annotation of 14 haplotype genomes with sizes spanning from 602.0 to 699.6 Mb. By compiling these haplotype genomes, we constructed a super pan-genome for the genus. We identified numerous structural variations (SVs, including variations in gene copy number) and highly diverged regions in these genomes. Notably, significant SV variability was observed within the intronic regions of the MED25 and TTG1 genes across different materials, suggesting their potential roles in influencing fruit size and trichome formation. Intriguingly, our findings indicated a high genetic divergence between two haplotype genomes, with one individual, tentatively named Actinidia × leiocacarpae, from sect. Leiocacarpae. This likely hybrid with a heterozygous genome exhibited notable genetic adaptations related to resistance against bacterial canker, particularly through the upregulation of the RPM1 gene, which contains a specific SV, after infection by Pseudomonas syringae pv. actinidiae. In addition, we also discussed the interlineage hybridizations and taxonomic treatments of the genus Actinidia. Overall, the comprehensive pan-genome constructed here, along with our findings, lays a foundation for examining genetic compositions and markers, particularly those related to SVs, to facilitate hybrid breeding aimed at developing desired phenotypes in kiwifruits.

Near-complete assembly and comprehensive annotation of the wheat Chinese Spring genome

A complete reference genome assembly is crucial for biological research and genetic improvement. Owing to its large size and highly repetitive nature, there are numerous gaps in the globally used wheat Chinese Spring (CS) genome assembly. In this study, we generated a 14.46 Gb near-complete assembly of the CS genome, with a contig N50 of over 266 Mb and an overall base accuracy of 99.9963%. Among the 290 gaps that remained (26, 257, and 7 gaps from the A, B, and D subgenomes, respectively), 278 were extremely high-copy tandem repeats, whereas the remaining 12 were transposable-element-associated gaps. Four chromosome assemblies were completely gap-free, including chr1D, chr3D, chr4D, and chr5D. Extensive annotation of the near-complete genome revealed 151 405 high-confidence genes, of which 59 180 were newly annotated, including 7602 newly assembled genes. Except for the centromere of chr1B, which has a gap associated with superlong GAA repeat arrays, the centromeric sequences of all of the remaining 20 chromosomes were completely assembled. Our near-complete assembly revealed that the extent of tandem repeats, such as simple-sequence repeats, was highly uneven among different subgenomes. Similarly, the repeat compositions of the centromeres also varied among the three subgenomes. With the genome sequences of all six types of seed storage proteins (SSPs) fully assembled, the expression of ω-gliadin was found to be contributed entirely by the B subgenome, whereas the expression of the other five types of SSPs was most abundant from the D subgenome. The near-complete CS genome will serve as a valuable resource for genomic and functional genomic research and breeding of wheat as well as its related species.

A telomere-to-telomere gap-free assembly integrating multi-omics uncovers the genetic mechanism of fruit quality and important agronomic trait associations in pomegranate

Pomegranate is an important perennial fruit tree distributed worldwide. Reference genomes with gaps and limit gene identification controlling important agronomic traits hinder its functional genomics and genetic improvements. Here, we reported a telomere-to-telomere (T2T) gap-free genome assembly of the distinctive cultivar 'Moshiliu'. The Moshiliu reference genome was assembled into eight chromosomes without gaps, totalling ~366.71 Mb, with 32 158 predicted protein-coding genes. All 16 telomeres and eight centromeres were characterized; combined with FISH analysis, we revealed the atypical telomere units in pomegranate as TTTTAGGG. Furthermore, a total of 16 loci associated with 15 important agronomic traits were identified based on GWAS of 146 accessions. Gene editing and biochemical experiments demonstrated that a 37.2-Kb unique chromosome translocation disrupting the coding domain sequence of PgANS was responsible for anthocyanin-less, knockout of PgANS in pomegranate exhibited a defect in anthocyanin production; a unique repeat expansion in the promoter of PgANR may affected its expression, resulting in black peel; notably, the G → A transversion located at the 166-bp coding domain of PgNST3, which caused a E56K mutation in the PgNST3 protein, closely linked with soft-seed trait. Overexpression of PgNST3A in tomato presented smaller and softer seed coats. The E56K mutation in PgNST3 protein, eliminated the binding ability of PgNST3 to the PgMYB46 promoter, which subsequently affected the thickness of the inner seed coat of soft-seeded pomegranates. Collectively, the validated gap-free genome, the identified genes controlling important traits and the CRISPR-Cas9-mediated gene knockout system all provided invaluable resources for pomegranate precise breeding.

Telomere-to-telomere genome and multi-omics analysis of Prunus avium cv. Tieton provides insights into its genomic evolution and flavonoid biosynthesis

The European sweet cherry (Prunus avium) is highly valued for its superior quality, delectable taste, and robust stress resistance, leading to its extensive cultivation in the world. However, the previous incomplete genome assemblies have impeded its evolution and genetic regulation studies. In this study, we generated a Telomere-to-Telomere gap-free genome assembly of P. avium cv. Tieton, using advanced sequencing technologies. The assembled genome comprises eight pseudochromosomes with a genome size of 342.23 Mb and a contig N50 of 40.66 Mb. Comparative genomic analysis identified several unique stress resistance-related genes, possibly associated with the species' environmental adaptation. The integrative analyses of genomics, transcriptomes and metabolomes identified some key structural genes and metabolites crucial to flavonoid biosynthesis of sweet cherry. Our analyses revealed that 85 flavonoid metabolites, which are highly differentially accumulated among five tissues (flesh, stem, leaf, bud, and seed) of cherry. Interestingly, eight abundant flavonoids (Narcissoside, Typhaneoside, Myricetin 3-0-galactoside, Diosmin, Neohesperidin, Liquiritin apioside, 5,6,7-Trimethoxyflavone and Oroxin B) were highly accumulated in cherry flesh tissues. The gene-metabolite correlation analysis revealed that seven genes (HTC8, HTC6, CYP75B1_9, CYP75B1_10, 4CL1, DFR1, and FLS1) significantly regulated flavonoid accumulation in cherry flesh. Additionally, some structural genes (4CL6, PAL3, CYP75A2, F3H1, CYP75B1_8, and CYP75B1_10) were identified in the flavonoid biosynthetic pathway and were highly expressed, aligning with high flavonoid metabolite content in cherry flesh. These identified genes and metabolites are likely pivotal in conferring sweet cherry's stress resistance and high-quality traits. These findings offer deep insights into the mechanisms of genomic evolution and flavonoid biosynthesis, which also lay a solid foundation for further function genomics studies and breeding improvement in cherry.

Chromosome-level genome assembly assisting for dissecting mechanism of anthocyanin regulation in kiwifruit (Actinidia arguta)

Actinidia arguta is a newly emerged, commercially cultivated Actinidia species. A. arguta has a beautiful appearance and is rich in anthocyanin, and is thus highly welcomed by consumers. However, the mechanism of anthocyanin regulation in A. arguta remains unclear. In this study, we assembled the nearly complete genome of the first red A. arguta cultivar, 'Tianyuanhong', with an N50 of 21 Mb. Comparative genome analysis revealed a role of the expansion/contraction of gene families in the species-specific trait formation of A. arguta. Through verification of transient overexpression and stable transformation, RNA-seq analysis revealed a key bHLH transcription factor, AaBEE1, which negatively regulates anthocyanin biosynthesis. DAP-seq analysis combined with Y1H, EMSA, Chip-qPCR and LUC suggested that AaBEE1 binds to the G-box of the AaLDOX promoter and suppresses its expression. Overall, we assembled the genome of A. arguta and clarified its AaBEE1-AaLDOX module-mediated molecular mechanism of anthocyanin regulation.

A Segregating Structural Variant Defines Novel Venom Phenotypes in the Eastern Diamondback Rattlesnake

Of all mutational mechanisms contributing to phenotypic variation, structural variants are both among the most capable of causing major effects as well as the most technically challenging to identify. Intraspecific variation in snake venoms is widely reported, and one of the most dramatic patterns described is the parallel evolution of streamlined neurotoxic rattlesnake venoms from hemorrhagic ancestors by means of deletion of snake venom metalloproteinase (SVMP) toxins and recruitment of neurotoxic dimeric phospholipase A2 (PLA2) toxins. While generating a haplotype-resolved, chromosome-level genome assembly for the eastern diamondback rattlesnake (Crotalus adamanteus), we discovered that our genome animal was heterozygous for a ∼225 Kb deletion containing six SVMP genes, paralleling one of the two steps involved in the origin of neurotoxic rattlesnake venoms. Range-wide population-genomic analysis revealed that, although this deletion is rare overall, it is the dominant homozygous genotype near the northwestern periphery of the species' range, where this species is vulnerable to extirpation. Although major SVMP deletions have been described in at least five other rattlesnake species, C. adamanteus is unique in not additionally gaining neurotoxic PLA2s. Previous work established a superficially complementary north-south gradient in myotoxin (MYO) expression based on copy number variation with high expression in the north and low in the south, yet we found that the SVMP and MYO genotypes vary independently, giving rise to an array of diverse, novel venom phenotypes across the range. Structural variation, therefore, forms the basis for the major axes of geographic venom variation for C. adamanteus.

A telomere-to-telomere genome assembly coupled with multi-omic data provides insights into the evolution of hexaploid bread wheat

The complete assembly of vast and complex plant genomes, like the hexaploid wheat genome, remains challenging. Here we present CS-IAAS, a comprehensive telomere-to-telomere (T2T) gap-free Triticum aestivum L. genome, encompassing 14.51 billion base pairs and featuring all 21 centromeres and 42 telomeres. Annotation revealed 90.8 Mb additional centromeric satellite arrays and 5,611 rDNA units. Genome-wide rearrangements, centromeric elements, transposable element expansion and segmental duplications were deciphered during tetraploidization and hexaploidization, providing a comprehensive understanding of wheat subgenome evolution. Among them, transposable element insertions during hexaploidization greatly influenced gene expression balances, thus increasing the genome plasticity of transcriptional levels. Additionally, we generated 163,329 full-length cDNA sequences and proteomic data that helped annotate 141,035 high-confidence protein-coding genes. The complete T2T reference genome (CS-IAAS), along with its transcriptome and proteome, represents a significant step in our understanding of wheat genome complexity and provides insights for future wheat research and breeding.

The haplotype-resolved genome assembly of an ancient citrus variety provides insights into the domestication history and fruit trait formation of loose-skin mandarins

BACKGROUND

Loose-skin mandarins (LSMs) are among the oldest domesticated horticultural crops, yet their domestication history and the genetic basis underlying the formation of key selected traits remain unclear.

RESULTS

We provide a chromosome-scale and haplotype-resolved assembly for the ancient Chinese citrus variety Nanfengmiju tangerine. Through the integration of 77 resequenced and 114 published citrus germplasm genomes, we categorize LSMs into 12 distinct groups based on population genomic analyses. We infer that the ancestors of modern cultivated mandarins diverged from wild mandarins in Daoxian approximately 500,000 years ago, when they entered the Yangtze and Pearl River Basins. There, they were domesticated into four ancient cultivation groups, forming the cornerstone of modern Chinese LSM cultivation. We identify selective sweeps underlying quantitative trait loci and genes related to important fruit quality traits, including sweetness and size. We reveal that the co-selection of sugar transporter and metabolism genes are associated with increased fruit sweetness. Significant alterations in the auxin and gibberellin signaling networks may contribute to the enlargement of LSM fruits. We also provide a comprehensive, high-spatiotemporal-resolution atlas of allelic gene expression during citrus fruit development. We detect 5890 allele pairs showing specific expression patterns and a significant increase in variation levels.

CONCLUSIONS

Our study provides valuable genomic resources and further revises the origin and domestication history of LSMs, offering insights for genetic improvement of citrus plants.

Genome assembly of Stewartia sinensis reveals origin and evolution of orphan genes in Theaceae

Orphan genes play crucial roles in diverse biological processes, but the evolutionary trajectories and functional divergence remain largely unexplored. The Theaceae family, including the economically and culturally important tea plant, offers a distinctive model to examine these aspects. Here, we integrated Nanopore long-read sequencing, Illumina short-read sequencing, and Hi-C methods to decode a pseudo-chromosomal genome assembly of Stewartia sinensis, from the earliest-diverging tribe of Theaceae, spanning 2.95 Gb. Comparative genomic analysis revealed the absence of recent whole-genome duplication events in the Theaceae ancestor, highlighting tandem duplications as the predominant mechanism of gene expansion. We identified 31,331 orphan genes, some of which appear to have ancient origins, suggesting early emergence with frequent gains and losses, while others seem more specific and recent. Notably, orphan genes are distinguished by shorter lengths, fewer exons and functional domains compared to genes that originate much earlier, like transcription factors. Moreover, tandem duplication contributes significantly to the adaptive evolution and characteristic diversity of Theaceae, and it is also a major mechanism driving the origination of orphan genes. This study illuminates the evolutionary dynamics of orphan genes, providing a valuable resource for understanding the origin and evolution of tea plant flavor and enhancing genetic breeding efforts.

A comprehensive analytical method 'Regulatome' revealed a novel pathway for aerenchyma formation under waterlogging in wheat

Waterlogging is a major abiotic stress restricting crop yield globally, and aerenchyma formation is one of the most important adaptive strategies in waterlogging-tolerant plants. However, the conservation of this process remains poorly understood, and additional pathways are yet to be identified. Here, physiological, anatomical, transcriptomic, and metabolomic analyses were conducted on wheat seedlings under normal and waterlogging conditions. Waterlogging caused growth inhibition and physiological damage, as well as induced aerenchyma formation in roots. A total of 10,346 differentially expressed genes and 3,419 differential metabolites were identified in roots. In addition to the AP2/ERF (APETALA2/ETHYLENE RESPONSIVE FACTOR) gene family, integrating analyses also revealed the role of LOB/AS2 (LATERAL ORGAN BOUNDARIES/ASYMMETRIC LEAVES2) in aerenchyma formation under waterlogging. It was revealed that the classical pathway of aerenchyma formation mediated by ethylene response, as well as synergy of calcium ion and reactive oxygen species, was deeply conserved in both monocots and eudicots during 160 million years of evolution through gene co-expression networks of cross-species. The newly introduced concept 'Regulatome' supported the classical pathway of aerenchyma formation, with a proposed model of the jasmonic acid signalling pathway involved in waterlogging, suggesting its usefulness in gene identification and function exploration. These findings provide a novel insight into the regulatory mechanisms of aerenchyma formation and breeding approaches for developing wheat cultivars with high waterlogging tolerance.

Plant sperm cell sequencing for genome phasing and determination of meiotic crossover points

Haplotype phasing represents a pivotal procedure in genome analysis, entailing the identification of specific genetic variant combinations on each chromosome. Achieving chromosome-level genome phasing constitutes a considerable challenge, particularly in organisms with large and complex genomes. To address this challenge, we have developed a robust, gamete cell-based phasing pipeline, including wet-laboratory processes for plant sperm cell isolation, short-read sequencing and a bioinformatics workflow to generate chromosome-level phasing. The bioinformatics workflow is applicable for both plant and other sperm cells, for example, those of mammals. Our pipeline ensures high-quality single-nucleotide polymorphism (SNP) calling for each sperm cell and the subsequent construction of a high-density genetic map. The genetic map facilitates accurate chromosome-level genome phasing, enables crossover event detection and could be used to correct potential assembly errors. Our bioinformatics pipeline runs on a Linux system and most of its steps can be executed in parallel, expediting the analysis process. The entire workflow can be performed over the course of 1 d. We provide a practical example from our previous research using this protocol and provide the whole bioinformatics pipeline as a Docker image to ensure its easy adaptability to other studies.

Chromatin Topological Domains Associate With the Rapid Formation of Tandem Duplicates in Plants

In eukaryotes, chromatin is compacted within nuclei under the principle of compartmentalization. On top of that, condensin II establishes eukaryotic chromosome territories, while cohesin organizes the vertebrate genome by extruding chromatin loops and forming topologically associating domains (TADs). Thus far, the formation and roles of these chromatin structures in plants remain poorly understood. This study integrates Hi-C data from diverse plant species, demonstrating that nuclear DNA content influences large-scale chromosome conformation and affects the finer details of compartmental patterns. These contrasting compartmental patterns are distinguished by gene-to-gene loops and validated through cytological observations. Additionally, a novel chromatin domain type associated with tandem duplicate gene clusters is identified. These domains are independent of H3K27me3-mediated chromatin compartmentalization and exhibit evolutionary conservation across species. Gene pairs within TAD-like domains are younger and show higher levels of coexpression. These domains potentially promote the formation of tandem duplicates, a property appears unique to the Actinidia family. Overall, this study reveals functional chromatin domains in plants and provides evidence for the role of three-dimensional chromatin architecture in gene regulation and genome evolution.

Super pan-genome reveals extensive genomic variations associated with phenotypic divergence in Actinidia

Kiwifruit is an economically and nutritionally important horticultural fruit crop worldwide. The genomic data of several kiwifruit species have been released, providing an unprecedented opportunity for pan-genome analysis to comprehensively investigate the inter- and intra-species genetic diversity and facilitate utilization for kiwifruit breeding. Here, we generated a kiwifruit super pan-genome using 15 high-quality assemblies of eight Actinidia species. For gene-based pan-genome, a total of 61,465 gene families were identified, and the softcore and dispensable genes were enriched in biological processes like response to endogenous stimulus, response to hormone and cell wall organization or biogenesis. Then, structural variations (SVs) against A. chinensis 'Donghong' were identified and then used to construct a graph-based genome. Further population-scale SVs based on resequencing data from 112 individuals of 20 species revealed extensive SVs which probably contributed to the phenotypic diversity among the Actinidia species. SV hotspot regions were found contributed to environmental adaptation. Furthermore, we systematically identified resistance gene analogs (RGAs) in the 15 assemblies and generated a pan-RGA dataset to reveal the diversity of genes potentially involved in disease resistance in Actinidia. The pan-genomic data obtained here is useful for evolutionary and functional genomic studies in Actinidia, and facilitates breeding design.

A telomere-to-telomere gap-free reference genome of Chionanthus retusus provides insights into the molecular mechanism underlying petal shape changes

Chionanthus retusus, an arbor tree of the Oleaceae family, is an ecologically and economically valuable ornamental plant for its remarkable adaptability in landscaping. During C. retusus breeding, we observed diverse floral shapes; however, no available genome for C. retusus has hindered the widespread identification of genes related to flower morphology. Thus, a de novo telomere-to-telomere (T2T) gap-free genome was generated. The assembly, incorporating high-coverage and long-read sequencing data, successfully yielded two complete haplotypes (687 and 683 Mb). The genome encompasses 42 864 predicted protein-coding genes, with all 46 telomeres and 23 centromeres in one haplotype. Whole-genome duplication analysis revealed that C. retusus underwent one fewer event of whole-genome duplication after differentiation compared to other species in the Oleaceae family. Furthermore, flower vein diversity was the main reason for the differences in floral shapes. Auxin-related genes were responsible for petal shape formation on genome-based transcriptome analysis. Specifically, the removal and retention of the first intron in CrAUX/IAA20 resulted in the production of two transcripts, and the differences in the expression levels of CrAUX/IAA20 resulted in the variations of flower veins. Compared to transcripts lacking the first intron, transcripts with intron retention caused more severe decreases in the number and length of flower veins in transgenic Arabidopsis thaliana. Our findings will deepen our understanding of flower morphology development and provide important theoretical support for the cultivation of Oleaceae.

Nanopore sequencing: flourishing in its teenage years

Over the past decade, nanopore sequencing has experienced significant advancements and changes, transitioning from an initially emerging technology to a significant instrument in the field of genomic sequencing. However, as advancements in next-generation sequencing technology persist, nanopore sequencing also improves. This paper reviews the developments, applications, and outlook on nanopore sequencing technology. Currently, nanopore sequencing supports both DNA and RNA sequencing, making it widely applicable in areas such as telomere-to-telomere (T2T) genome assembly, direct RNA sequencing (DRS), and metagenomics. The openness and versatility of nanopore sequencing have established it as a preferred option for an increasing number of research teams, signaling a transformative influence on life science research. As the nanopore sequencing technology advances, it provides a faster, more cost-effective approach with extended read lengths, demonstrating the significant potential for complex genome assembly, pathogen detection, environmental monitoring, and human disease research, offering a fresh perspective in sequencing technologies.

Genome-wide identification of SWEET gene family and the sugar transport function of three candidate genes during female flower bud induction stage of Juglans sigillata Dode

Sugar Will Eventually be Exported Transporter (SWEET) transports sugar to sink organs and regulates intercellular sugar transport to provide energy for plant growth and development. In this study, twenty-two SWEET genes were identified and distributed on eleven chromosomes. Phylogenetic tree analysis showed that these genes could be divided into four subfamilies. Metabolism and transcriptome analysis showed that sucrose and fructose were accumulated in female flower buds at physiological differentiation stage (PDS). The third branch of JsSWEET1 and the fourth branch of JsSWEET9 and JsSWEET17 were highly expressed in female flower buds at undifferentiated stage (UDS) and PDS, which promoted sugar accumulation in female flower bud differentiation, so these three candidate SWEET genes were considered to be involved in the accumulation of sugar in the flower bud differentiation of Juglans sigillata. The subcellular localization showed that all three candidate genes were located on the cell membrane, and JsSWEET17 was also expressed and functioned in the vacuolar membrane. Through overexpression in callus and silencing in female flower buds at UDS and PDS, it was found that JsSWEET1 positively regulated the accumulation of sucrose in female flower buds, and JsSWEET9 and JsSWEET17 are involved in the transport and accumulation of fructose during flower bud differentiation. These results could provide a comprehensive understanding of the JsSWEETs gene family and provide theoretical guidance for further study of the function of SWEET-induced sugar accumulation in plant flower development.

Genotype-associated core bacteria enhance host resistance against kiwifruit bacterial canker

Both the phyllosphere and rhizosphere are inhabited by different kinds of microorganisms that are closely related to plant growth and health. However, it is not clear whether disease-resistant cultivars shape the microbiome to facilitate disease resistance. In this study, significant differences were found in the aboveground and belowground bacterial communities of disease-resistant and disease-susceptible cultivars grown in the same kiwifruit orchard. The phyllosphere of the resistant cultivar 'Wanjin' showed greater enrichment of Pseudomonas spp. and Sphingomonas spp. than the susceptible cultivar 'Donghong'. The rhizosphere microbes of 'Wanjin' were less affected by field location, with significantly greater bacterial abundance than those of 'Donghong' and more bacteria with potential biocontrol properties. Pseudomonas syringae pv. actinidiae (Psa) infection significantly affected the microbiome of the phyllosphere of kiwifruit plants, especially that of 'Donghong'. Resistant and susceptible kiwifruit cultivars exhibit distinct beneficial microbial recruitment strategies under Psa challenge. The phyllosphere of 'Donghong' in Jinzhai was enriched with Sphingomonas spp. and Pantoea spp. under Psa infection, while the rhizosphere of 'Wanjin' was enriched with Sphingomonas spp. and Novosphingobium spp. We further identified five key biomarkers within the microbial community associated with Psa infection. Inoculation experiments showed that Lysobacter sp. R34, Stenotrophomonas sp. R31, Pseudomonas sp. R10 and RS54, which were isolated from belowground compartments of 'Wanjin', could positively affect plant performance under Psa challenge. The combination use of Pseudomonas sp. R10 and Stenotrophomonas sp. R31 significantly improve the management of kiwifruit canker. Our findings provided novel insights into soil-microbe-plant interactions and the role of microbes in plant disease resistance and susceptibility.

Genomes of diverse Actinidia species provide insights into cis-regulatory motifs and genes associated with critical traits

BACKGROUND

Kiwifruit, belonging to the genus Actinidia, represents a unique fruit crop characterized by its modern cultivars being genetically diverse and exhibiting remarkable variations in morphological traits and adaptability to harsh environments. However, the genetic mechanisms underlying such morphological diversity remain largely elusive.

RESULTS

We report the high-quality genomes of five Actinidia species, including Actinidia longicarpa, A. macrosperma, A. polygama, A. reticulata, and A. rufa. Through comparative genomics analyses, we identified three whole genome duplication events shared by the Actinidia genus and uncovered rapidly evolving gene families implicated in the development of characteristic kiwifruit traits, including vitamin C (VC) content and fruit hairiness. A range of structural variations were identified, potentially contributing to the phenotypic diversity in kiwifruit. Notably, phylogenomic analyses revealed 76 cis-regulatory elements within the Actinidia genus, predominantly associated with stress responses, metabolic processes, and development. Among these, five motifs did not exhibit similarity to known plant motifs, suggesting the presence of possible novel cis-regulatory elements in kiwifruit. Construction of a pan-genome encompassing the nine Actinidia species facilitated the identification of gene DTZ79_23g14810 specific to species exhibiting extraordinarily high VC content. Expression of DTZ79_23g14810 is significantly correlated with the dynamics of VC concentration, and its overexpression in the transgenic roots of kiwifruit plants resulted in increased VC content.

CONCLUSIONS

Collectively, the genomes and pan-genome of diverse Actinidia species not only enhance our understanding of fruit development but also provide a valuable genomic resource for facilitating the genome-based breeding of kiwifruit.

Genome-wide identification of kiwifruit K+ channel Shaker family members and their response to low-K+ stress

BACKGROUND

'Hongyang' kiwifruit (Actinidia chinensis cv 'Hongyang') is a high-quality variety of A. chinensis with the advantages of high yield, early ripening, and high stress tolerance. Studies have confirmed that the Shaker K+ genes family is involved in plant uptake and distribution of potassium (K+).

RESULTS

Twenty-eight Shaker genes were identified and analyzed from the 'Hongyang' kiwifruit (A. chinensis cv 'Hongyang') genome. Subcellular localization results showed that more than one-third of the AcShaker genes were on the cell membrane. Phylogenetic analysis indicated that the AcShaker genes were divided into six subfamilies (I-VI). Conservative model, gene structure, and structural domain analyses showed that AcShaker genes of the same subfamily have similar sequence features and structure. The promoter cis-elements of the AcShaker genes were classified into hormone-associated cis-elements and environmentally stress-associated cis-elements. The results of chromosomal localization and duplicated gene analysis demonstrated that AcShaker genes were distributed on 18 chromosomes, and segmental duplication was the prime mode of gene duplication in the AcShaker family. GO enrichment analysis manifested that the ion-conducting pathway of the AcShaker family plays a crucial role in regulating plant growth and development and adversity stress. Compared with the transcriptome data of the control group, all AcShaker genes were expressed under low-K+stress, and the expression differences of three genes (AcShaker15, AcShaker17, and AcShaker22) were highly significant. The qRT-PCR results showed a high correlation with the transcriptome data, which indicated that these three differentially expressed genes could regulate low-K+ stress and reduce K+ damage in kiwifruit plants, thus improving the resistance to low-K+ stress. Comparison between the salt stress and control transcriptomic data revealed that the expression of AcShaker11 and AcShaker18 genes was significantly different and lower under salt stress, indicating that both genes could be involved in salt stress resistance in kiwifruit.

CONCLUSION

The results showed that 28 AcShaker genes were identified and all expressed under low K+ stress, among which AcShaker22 was differentially and significantly upregulated. The AcShaker22 gene can be used as a candidate gene to cultivate new varieties of kiwifruit resistant to low K+ and provide a reference for exploring more properties and functions of the AcShaker genes.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

Changes in fruit quality properties and phytochemical substances of kiwifruit (Actinidia deliciosa) grown in different agro-ecological conditions during cold storage

BACKGROUND

The changes in the physical structures of the products are the first things that consumers pay attention to. Therefore, it is essential and significant importance to take measures to improve the storage conditions of products and to minimize quality losses. The main objective of the study was to evaluate the effects of agro-ecological conditions on bioactive compounds and fruit quality of kiwifruit during cold storage. The 'Hayward' kiwifruit cultivar grown in Ordu, Giresun, Samsun, Rize, and Yalova provinces of Türkiye were kept at 0 ± 0.5 °C and relative humidity of 90 ± 5% for 150 d.

RESULTS

The kiwifruit obtained from the provinces of Yalova, Ordu, and Giresun experienced the least weight loss during cold storage. Kiwifruit from Samsun and Yalova provinces had the lowest fruit firmness, while those from Giresun had the highest on 150th d. The changes were observed in the skin and flesh colors of the kiwifruit belonging to all cultivation areas. The amount of vitamin C increased throughout the study in all ecological conditions, but the Yalova province's kiwifruit was found to have the highest levels. Additionally, in all ecologies, kiwifruit showed an increase in antioxidant activity, total phenolics, and total flavonoids, all known to have beneficial effects on human health. The total antioxidant activity and total phenolics were highest in the kiwifruit of Yalova province, but the total flavonoids were found in the kiwifruit of Rize and Ordu provinces.

CONCLUSION

The study's results revealed that kiwifruit's bioactive compounds and quality parameters may vary depending on the cultivation area. Additionally, it can be stated that Yalova province kiwifruit experiences the least amount of postharvest quality losses.

Comparative transcriptome and metabolome reveal the role of acidic electrolyzed oxidizing water in improving postharvest disease resistance of longan fruit

Acidic electrolyzed oxidizing water (AEOW) was applied to suppress disease development and maintain good quality of fresh fruit. However, the involvement of AEOW in improving disease resistance of fresh longan remains unknown. Here, transcriptomic and metabolic analyses were performed to compare non-treated and AEOW-treated longan during storage. The transcriptome analysis showed AEOW-induced genes associated with phenylpropanoid and flavonoid biosynthesis. The metabolome analysis found the contents of coumarin, phenolic acid, and tannin maintained higher levels in AEOW-treated longan than non-treated longan. Moreover, the weighted correlation network analysis (WGCNA) was performed to identify hub genes, and a gene-metabolite correlation network associated with AEOW-improved disease resistance in longan was constructed by the co-analysis of transcriptomics and metabolomics. These findings identified a series of important genes and metabolites involving in AEOW-induced disease resistance of longan fruit, expanding our knowledges on fruit disease resistance and quality maintenance at the transcript and metabolic levels.

Divergent Roles of the Auxin Response Factors in Lemongrass (Cymbopogon flexuosus (Nees ex Steud.) W. Watson) during Plant Growth

Auxin Response Factors (ARFs) make up a plant-specific transcription factor family that mainly couples perception of the phytohormone, auxin, and gene expression programs and plays an important and multi-faceted role during plant growth and development. Lemongrass (Cymbopogon flexuosus) is a representative Cymbopogon species widely used in gardening, beverages, fragrances, traditional medicine, and heavy metal phytoremediation. Biomass yield is an important trait for several agro-economic purposes of lemongrass, such as landscaping, essential oil production, and phytoremediation. Therefore, we performed gene mining of CfARFs and identified 26 and 27 CfARF-encoding genes in each of the haplotype genomes of lemongrass, respectively. Phylogenetic and domain architecture analyses showed that CfARFs can be divided into four groups, among which groups 1, 2, and 3 correspond to activator, repressor, and ETTN-like ARFs, respectively. To identify the CfARFs that may play major roles during the growth of lemongrass plants, RNA-seq was performed on three tissues (leaf, stem, and root) and four developmental stages (3-leaf, 4-leaf, 5-leaf. and mature stages). The expression profiling of CfARFs identified several highly expressed activator and repressor CfARFs and three CfARFs (CfARF3, 18, and 35) with gradually increased levels during leaf growth. Haplotype-resolved transcriptome analysis revealed that biallelic expression dominance is frequent among CfARFs and contributes to their gene expression patterns. In addition, co-expression network analysis identified the modules enriched with CfARFs. By establishing orthologous relationships among CfARFs, sorghum ARFs, and maize ARFs, we showed that CfARFs were mainly expanded by whole-genome duplications, and that the duplicated CfARFs might have been divergent due to differential expression and variations in domains and motifs. Our work provides a detailed catalog of CfARFs in lemongrass, representing a first step toward characterizing CfARF functions, and may be useful in molecular breeding to enhance lemongrass plant growth.

A telomere-to-telomere gap-free reference genome assembly of avocado provides useful resources for identifying genes related to fatty acid biosynthesis and disease resistance

Avocado (Persea americana Mill.) is an economically valuable plant because of the high fatty acid content and unique flavor of its fruits. Its fatty acid content, especially the relatively high unsaturated fatty acid content, provides significant health benefits. We herein present a telomere-to-telomere gapless genome assembly (841.6 Mb) of West Indian avocado. The genome contains 40 629 predicted protein-coding genes. Repeat sequences account for 57.9% of the genome. Notably, all telomeres, centromeres, and a nucleolar organizing region are included in this genome. Fragments from these three regions were observed via fluorescence in situ hybridization. We identified 376 potential disease resistance-related nucleotide-binding leucine-rich repeat genes. These genes, which are typically clustered on chromosomes, may be derived from gene duplication events. Five NLR genes (Pa11g0262, Pa02g4855, Pa07g3139, Pa07g0383, and Pa02g3196) were highly expressed in leaves, stems, and fruits, indicating they may be involved in avocado disease responses in multiple tissues. We also identified 128 genes associated with fatty acid biosynthesis and analyzed their expression patterns in leaves, stems, and fruits. Pa02g0113, which encodes one of 11 stearoyl-acyl carrier protein desaturases mediating C18 unsaturated fatty acid synthesis, was more highly expressed in the leaves than in the stems and fruits. These findings provide valuable insights that enhance our understanding of fatty acid biosynthesis in avocado.

Graph-Based Pangenome of Actinidia chinensis Reveals Structural Variations Mediating Fruit Degreening

Fruit ripening is associated with the degreening process (loss of chlorophyll) that occurs in most fruit species. Kiwifruit is one of the special species whose fruits may maintain green flesh by accumulating a large amount of chlorophyll even after ripening. However, little is known about the genetic variations related to the fruit degreening process. Here, a graph-based kiwifruit pangenome by analyzing 14 chromosome-scale haplotype-resolved genome assemblies from seven representative cultivars or lines in Actinidia chinensis is built. A total of 49,770 non-redundant gene families are identified, with core genes constituting 46.6%, and dispensable genes constituting 53.4%. A total of 84,591 non-redundant structural variations (SVs) are identified. The pangenome graph integrating both reference genome sequences and variant information facilitates the identification of SVs related to fruit color. The SV in the promoter of the AcBCM gene determines its high expression in the late developmental stage of fruits, which causes chlorophyll accumulation in the green-flesh fruits by post-translationally regulating AcSGR2, a key enzyme of chlorophyll catabolism. Taken together, a high-quality pangenome is constructed, unraveled numerous genetic variations, and identified a novel SV mediating fruit coloration and fruit quality, providing valuable information for further investigating genome evolution and domestication, QTL genes function, and genomics-assisted breeding.

Genome assembly of autotetraploid Actinidia arguta highlights adaptive evolution and enables dissection of important economic traits

Actinidia arguta, the most widely distributed Actinidia species and the second cultivated species in the genus, can be distinguished from the currently cultivated Actinidia chinensis on the basis of its small and smooth fruit, rapid softening, and excellent cold tolerance. Adaptive evolution of tetraploid Actinidia species and the genetic basis of their important agronomic traits are still unclear. Here, we generated a chromosome-scale genome assembly of an autotetraploid male A. arguta accession. The genome assembly was 2.77 Gb in length with a contig N50 of 9.97 Mb and was anchored onto 116 pseudo-chromosomes. Resequencing and clustering of 101 geographically representative accessions showed that they could be divided into two geographic groups, Southern and Northern, which first diverged 12.9 million years ago. A. arguta underwent two prominent expansions and one demographic bottleneck from the mid-Pleistocene climate transition to the late Pleistocene. Population genomics studies using paleoclimate data enabled us to discern the evolution of the species' adaptation to different historical environments. Three genes (AaCEL1, AaPME1, and AaDOF1) related to flesh softening were identified by multi-omics analysis, and their ability to accelerate flesh softening was verified through transient expression assays. A set of genes that characteristically regulate sexual dimorphism located on the sex chromosome (Chr3) or autosomal chromosomes showed biased expression during stamen or carpel development. This chromosome-level assembly of the autotetraploid A. arguta genome and the genes related to important agronomic traits will facilitate future functional genomics research and improvement of A. arguta.

Chromosome-scale genome assembly of oil-tea tree Camellia crapnelliana

Camellia crapnelliana Tutch., belonging to the Theaceae family, is an excellent landscape tree species with high ornamental values. It is particularly an important woody oil-bearing plant species with high ecological, economic, and medicinal values. Here, we first report the chromosome-scale reference genome of C. crapnelliana with integrated technologies of SMRT, Hi-C and Illumina sequencing platforms. The genome assembly had a total length of ~2.94 Gb with contig N50 of ~67.5 Mb, and ~96.34% of contigs were assigned to 15 chromosomes. In total, we predicted 37,390 protein-coding genes, ~99.00% of which could be functionally annotated. The chromosome-scale genome of C. crapnelliana will become valuable resources for understanding the genetic basis of the fatty acid biosynthesis, and greatly facilitate the exploration and conservation of C. crapnelliana.

Haplotype-resolved chromosomal-level genome assembly reveals regulatory variations in mulberry fruit anthocyanin content

Understanding the intricate regulatory mechanisms underlying the anthocyanin content (AC) in fruits and vegetables is crucial for advanced biotechnological customization. In this study, we generated high-quality haplotype-resolved genome assemblies for two mulberry cultivars: the high-AC 'Zhongsang5801' (ZS5801) and the low-AC 'Zhenzhubai' (ZZB). Additionally, we conducted a comprehensive analysis of genes associated with AC production. Through genome-wide association studies (GWAS) on 112 mulberry fruits, we identified MaVHAG3, which encodes a vacuolar-type H+-ATPase G3 subunit, as a key gene linked to purple pigmentation. To gain deeper insights into the genetic and molecular processes underlying high AC, we compared the genomes of ZS5801 and ZZB, along with fruit transcriptome data across five developmental stages, and quantified the accumulation of metabolic substances. Compared to ZZB, ZS5801 exhibited significantly more differentially expressed genes (DEGs) related to anthocyanin metabolism and higher levels of anthocyanins and flavonoids. Comparative analyses revealed expansions and contractions in the flavonol synthase (FLS) and dihydroflavonol 4-reductase (DFR) genes, resulting in altered carbon flow. Co-expression analysis demonstrated that ZS5801 displayed more significant alterations in genes involved in late-stage AC regulation compared to ZZB, particularly during the phase stage. In summary, our findings provide valuable insights into the regulation of mulberry fruit AC, offering genetic resources to enhance cultivars with higher AC traits.

High-integrity Pueraria montana var. lobata genome and population analysis revealed the genetic diversity of Pueraria genus

Pueraria montana var. lobata (P. lobata) is a traditional medicinal plant belonging to the Pueraria genus of Fabaceae family. Pueraria montana var. thomsonii (P. thomsonii) and Pueraria montana var. montana (P. montana) are its related species. However, evolutionary history of the Pueraria genus is still largely unknown. Here, a high-integrity, chromosome-level genome of P. lobata and an improved genome of P. thomsonii were reported. It found evidence for an ancient whole-genome triplication and a recent whole-genome duplication shared with Fabaceae in three Pueraria species. Population genomics of 121 Pueraria accessions demonstrated that P. lobata populations had substantially higher genetic diversity, and P. thomsonii was probably derived from P. lobata by domestication as a subspecies. Selection sweep analysis identified candidate genes in P. thomsonii populations associated with the synthesis of auxin and gibberellin, which potentially play a role in the expansion and starch accumulation of tubers in P. thomsonii. Overall, the findings provide new insights into the evolutionary and domestication history of the Pueraria genome and offer a valuable genomic resource for the genetic improvement of these species.

Haplotype-resolved gapless genome and chromosome segment substitution lines facilitate gene identification in wild rice

The abundant genetic variation harbored by wild rice (Oryza rufipogon) has provided a reservoir of useful genes for rice breeding. However, the genome of wild rice has not yet been comprehensively assessed. Here, we report the haplotype-resolved gapless genome assembly and annotation of wild rice Y476. In addition, we develop two sets of chromosome segment substitution lines (CSSLs) using Y476 as the donor parent and cultivated rice as the recurrent parents. By analyzing the gapless reference genome and CSSL population, we identify 254 QTLs associated with agronomic traits, biotic and abiotic stresses. We clone a receptor-like kinase gene associated with rice blast resistance and confirm its wild rice allele improves rice blast resistance. Collectively, our study provides a haplotype-resolved gapless reference genome and demonstrates a highly efficient platform for gene identification from wild rice.

Haplotype-resolved T2T reference genomes for wild and domesticated accessions shed new insights into the domestication of jujube

Chinese jujube (Ziziphus jujuba Mill.) is one of the most important deciduous tree fruits in China, with substantial economic and nutritional value. Jujube was domesticated from its wild progenitor, wild jujube (Z. jujuba var. spinosa), and both have high medicinal value. Here we report the 767.81- and 759.24-Mb haplotype-resolved assemblies of a dry-eating 'Junzao' jujube (JZ) and a wild jujube accession (SZ), using a combination of multiple sequencing strategies. Each assembly yielded two complete haplotype-resolved genomes at the telomere-to-telomere (T2T) level, and ~81.60 and 69.07 Mb of structural variations were found between the two haplotypes within JZ and SZ, respectively. Comparative genomic analysis revealed a large inversion on each of chromosomes 3 and 4 between JZ and SZ, and numerous genes were affected by structural variations, some of which were associated with starch and sucrose metabolism. A large-scale population analysis of 672 accessions revealed that wild jujube originated from the lower reaches of the Yellow River and was initially domesticated at local sites. It spread widely and was then independently domesticated at the Shanxi-Shaanxi Gorge of the middle Yellow River. In addition, we identified some new selection signals regions on genomes, which are involved in the tissue development, pollination, and other aspects of jujube tree morphology and fertilization domestication. In conclusion, our study provides high-quality reference genomes of jujube and wild jujube and new insights into the domestication history of jujube.

Integrated analysis of metabolome, transcriptome, and bioclimatic factors of Acer truncatum seeds reveals key candidate genes related to unsaturated fatty acid biosynthesis, and potentially optimal production area

BACKGROUND

Lipids found in plant seeds are essential for controlling seed dormancy, dispersal, and defenses against biotic and abiotic stress. Additionally, these lipids provide nutrition and energy and are therefore important to the human diet as edible oils. Acer truncatum, which belongs to the Aceaceae family, is widely cultivated around the world for its ornamental value. Further because its seed oil is rich in unsaturated fatty acids (UFAs)- i.e. α-linolenic acid (ALA) and nervonic acid (NA)- and because it has been validated as a new food resource in China, the importance of A. truncatum has greatly risen. However, it remains unknown how UFAs are biosynthesized during the growth season, to what extent environmental factors impact their content, and what areas are potentially optimal for their production.

RESULTS

In this study, transcriptome and metabolome of A. truncatum seeds at three representative developmental stages was used to find the accumulation patterns of all major FAs. Cumulatively, 966 metabolites and 87,343 unigenes were detected; the differential expressed unigenes and metabolites were compared between stages as follows: stage 1 vs. 2, stage 1 vs. 3, and stage 2 vs. 3 seeds, respectively. Moreover, 13 fatty acid desaturases (FADs) and 20 β-ketoacyl-CoA synthases (KCSs) were identified, among which the expression level of FAD3 (Cluster-7222.41455) and KCS20 (Cluster-7222.40643) were consistent with the metabolic results of ALA and NA, respectively. Upon analysis of the geographical origin-affected diversity from 17 various locations, we found significant variation in phenotypes and UFA content. Notably, in this study we found that 7 bioclimatic variables showed considerable influence on FAs contents in A. truncatum seeds oil, suggesting their significance as critical environmental parameters. Ultimately, we developed a model for potentially ecological suitable regions in China.

CONCLUSION

This study provides a comprehensive understanding of the relationship between metabolome and transcriptome in A. truncatum at various developmental stages of seeds and a new strategy to enhance seed FA content, especially ALA and NA. This is particularly significant in meeting the increasing demands for high-quality edible oil for human consumption. The study offers a scientific basis for A. truncatum's novel utilization as a woody vegetable oil rather than an ornamental plant, potentially expanding its cultivation worldwide.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

A haplotype-resolved gap-free genome assembly provides novel insight into monoterpenoid diversification in Mentha suaveolens 'Variegata'

Mentha is a commonly used spice worldwide, which possesses medicinal properties and fragrance. These characteristics are conferred, at least partially, by essential oils such as menthol. In this study, a gap-free assembly with a genome size of 414.3 Mb and 31,251 coding genes was obtained for Mentha suaveolens 'Variegata'. Based on its high heterozygosity (1.5%), two complete haplotypic assemblies were resolved, with genome sizes of 401.9 and 405.7 Mb, respectively. The telomeres and centromeres of each haplotype were almost fully annotated. In addition, we detected a total of 41,135 structural variations. Enrichment analysis demonstrated that genes involved in terpenoid biosynthesis were affected by these structural variations. Analysis of volatile metabolites showed that M. suaveolens mainly produces piperitenone oxide rather than menthol. We identified three genes in the M. suaveolens genome which encode isopiperitenone reductase (ISPR), a key rate-limiting enzyme in menthol biosynthesis. However, the transcription levels of ISPR were low. Given that other terpenoid biosynthesis genes were expressed, M. suaveolens ISPRs may account for the accumulation of piperitenone oxide in this species. The findings of this study may provide a valuable resource for improving the detection rate and accuracy of genetic variants, thereby enhancing our understanding of their impact on gene function and expression. Moreover, our haplotype-resolved gap-free genome assembly offers novel insights into molecular marker-assisted breeding of Mentha.

Haplotype-resolved genome assembly provides insights into evolutionary history of the Actinidia arguta tetraploid

Actinidia arguta, known as hardy kiwifruit, is a widely cultivated species with distinct botanical characteristics such as small and smooth-fruited, rich in beneficial nutrients, rapid softening and tolerant to extremely low temperatures. It contains the most diverse ploidy types, including diploid, tetraploid, hexaploid, octoploid, and decaploid. Here we report a haplotype-resolved tetraploid genome (A. arguta cv. 'Longcheng No.2') containing four haplotypes, each with 40,859, 41,377, 39,833 and 39,222 protein-coding genes. We described the phased genome structure, synteny, and evolutionary analyses to identify and date possible WGD events. Ks calculations for both allelic and paralogous genes pairs throughout the assembled haplotypic individuals showed its tetraploidization is estimated to have formed ~ 1.03 Mya following Ad-α event occurred ~ 18.7 Mya. Detailed annotations of NBS-LRRs or CBFs highlight the importance of genetic variations coming about after polyploidization in underpinning ability of immune responses or environmental adaptability. WGCNA analysis of postharvest quality indicators in combination with transcriptome revealed several transcription factors were involved in regulating ripening kiwi berry texture. Taking together, the assembly of an A. arguta tetraploid genome provides valuable resources in deciphering complex genome structure and facilitating functional genomics studies and genetic improvement for kiwifruit and other crops.

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.

AcMYB10 Involved in Anthocyanin Regulation of 'Hongyang' Kiwifruit Induced via Fruit Bagging and High-Postharvest-Temperature Treatments

Light and temperature are key factors influencing the accumulation of anthocyanin in fruit crops. To assess the effects of fruit bagging during development and high post-ripening temperature on 'Hongyang' kiwifruit, we compared the pigmentation phenotypes and expression levels of anthocyanin-related genes between bagged and unbagged treatments, and between 25 °C and 37 °C postharvest storage temperatures. Both the bagging and 25 °C treatments showed better pigmentation phenotypes with higher anthocyanin concentrations. The results of the qRT-PCR analysis revealed that the gene expression levels of LDOX (leucoanthocyanidin dioxygenase), F3GT (UDP-flavonoid 3-O-glycosyltransferase ), AcMYB10, and AcbHLH42 were strongly correlated and upregulated by both the bagging treatment and 25 °C storage. The results of bimolecular fluorescence complementation and luciferase complementation imaging assays indicated an interaction between AcMYB10 and AcbHLH42 in plant cells, whereas the results of a yeast one-hybrid assay further demonstrated that AcMYB10 activated the promoters of AcLODX and AcF3GT. These results strongly suggest that enhanced anthocyanin synthesis is caused by the promoted expression of AcLODX and AcF3GT, regulated by the complex formed by AcMYB10-AcbHLH42.

The telomere-to-telomere (T2T) genome provides insights into the evolution of specialized centromere sequences in sandalwood

BACKGROUND

Sandalwood, a prized hemiparasitic plant, is highly sought in the commercial market because of its aromatic core materia. The structure and stability of the genome are instrumental in the rapid adaptation of parasitic plants to their surroundings. However, there is a conspicuous lack of research on the genomic-level adaptive evolution of sandalwood.

RESULTS

In this study, we assembled a gap-free telomere-to-telomere (T2T) reference genome for Santalum album using PacBio HiFi, Hi-C, and ultra-long ONT data. The T2T reference genome (Sal_t2t) encompassed annotations of 24,171 genes and 25.34% repetitive sequences, in addition to all 10 centromeres and 20 telomeres across the 10 chromosomes. The results revealed that the 3 distinct parasitic species of Santalales had diverse centromeric compositions. The Copia-type long terminal repeat transposon emerged as the most significant in the S. album genome, constituting the primary sequence of the centromere and influencing gene expression. Third, in sandalwood, the presence of Copia affected the size of the centromeres and, consequently, the genome size. Identification of the sandalwood T2T genome in this study also enabled the identification of more precise organelle transfer fragments.

CONCLUSIONS

Our research provides a sandalwood T2T genome, laying the groundwork for future investigations on the evolution of energy organs in parasitic plants. Moreover, it offers novel insights into the function and evolution of centromeres, as well as the mechanisms of adaptation and parasitism.

The haplotype-resolved telomere-to-telomere carnation (Dianthus caryophyllus) genome reveals the correlation between genome architecture and gene expression

Carnation (Dianthus caryophyllus) is one of the most valuable commercial flowers, due to its richness of color and form, and its excellent storage and vase life. The diverse demands of the market require faster breeding in carnations. A full understanding of carnations is therefore required to guide the direction of breeding. Hence, we assembled the haplotype-resolved gap-free carnation genome of the variety 'Baltico', which is the most common white standard variety worldwide. Based on high-depth HiFi, ultra-long nanopore, and Hi-C sequencing data, we assembled the telomere-to-telomere (T2T) genome to be 564 479 117 and 568 266 215 bp for the two haplotypes Hap1 and Hap2, respectively. This T2T genome exhibited great improvement in genome assembly and annotation results compared with the former version. The improvements were seen when different approaches to evaluation were used. Our T2T genome first informs the analysis of the telomere and centromere region, enabling us to speculate about specific centromere characteristics that cannot be identified by high-order repeats in carnations. We analyzed allele-specific expression in three tissues and the relationship between genome architecture and gene expression in the haplotypes. This demonstrated that the length of the genes, coding sequences, and introns, the exon numbers and the transposable element insertions correlate with gene expression ratios and levels. The insertions of transposable elements repress expression in gene regulatory networks in carnation. This gap-free finished T2T carnation genome provides a valuable resource to illustrate the genome characteristics and for functional genomics analysis in further studies and molecular breeding.

Determination of the effects of pre-harvest bagging treatment on kiwifruit appearance and quality via transcriptome and metabolome analyses

Bagging is an effective cultivation strategy to produce attractive and pollution-free kiwifruit. However, the effect and metabolic regulatory mechanism of bagging treatment on kiwifruit quality remain unclear. In this study, transcriptome and metabolome analyses were conducted to determine the regulatory network of the differential metabolites and genes after bagging. Using outer and inner yellow single-layer fruit bags, we found that bagging treatment improved the appearance of kiwifruit, increased the soluble solid content (SSC) and carotenoid and anthocyanin levels, and decreased the chlorophyll levels. We also identified 41 differentially expressed metabolites and 897 differentially expressed genes (DEGs) between the bagged and control 'Hongyang' fruit. Transcriptome and metabolome analyses revealed that the increase in SSC after bagging treatment was mainly due to the increase in D-glucosamine metabolite levels and eight DEGs involved in amino sugar and nucleotide sugar metabolic pathways. A decrease in glutamyl-tRNA reductase may be the main reason for the decrease in chlorophyll. Downregulation of lycopene epsilon cyclase and 9-cis-epoxycarotenoid dioxygenase increased carotenoid levels. Additionally, an increase in the levels of the taxifolin-3'-O-glucoside metabolite, flavonoid 3'-monooxygenase, and some transcription factors led to the increase in anthocyanin levels. This study provides novel insights into the effects of bagging on the appearance and internal quality of kiwifruit and enriches our theoretical knowledge on the regulation of color pigment synthesis in kiwifruit.

Telomere-to-telomere pear (Pyrus pyrifolia) reference genome reveals segmental and whole genome duplication driving genome evolution

Previously released pear genomes contain a plethora of gaps and unanchored genetic regions. Here, we report a telomere-to-telomere (T2T) gap-free genome for the red-skinned pear, 'Yunhong No. 1' (YH1; Pyrus pyrifolia), which is mainly cultivated in Yunnan Province (southwest China), the pear's primary region of origin. The YH1 genome is 501.20 Mb long with a contig N50 length of 29.26 Mb. All 17 chromosomes were assembled to the T2T level with 34 characterized telomeres. The 17 centromeres were predicted and mainly consist of centromeric-specific monomers (CEN198) and long terminal repeat (LTR) Gypsy elements (≥74.73%). By filling all unclosed gaps, the integrity of YH1 is markedly improved over previous P. pyrifolia genomes ('Cuiguan' and 'Nijisseiki'). A total of 1531 segmental duplication (SD) driven duplicated genes were identified and enriched in stress response pathways. Intrachromosomal SDs drove the expansion of disease resistance genes, suggesting the potential of SDs in adaptive pear evolution. A large proportion of duplicated gene pairs exhibit dosage effects or sub-/neo-functionalization, which may affect agronomic traits like stone cell content, sugar content, and fruit skin russet. Furthermore, as core regulators of anthocyanin biosynthesis, we found that MYB10 and MYB114 underwent various gene duplication events. Multiple copies of MYB10 and MYB114 displayed obvious dosage effects, indicating role differentiation in the formation of red-skinned pear fruit. In summary, the T2T gap-free pear genome provides invaluable resources for genome evolution and functional genomics.

Genome of Phyllanthus emblica: the medicinal plant Amla with super antioxidant properties

Phyllanthus emblica or Indian gooseberry, commonly known as amla, is an important medicinal horticultural plant used in traditional and modern medicines. It bears stone fruits with immense antioxidant properties due to being one of the richest natural sources of vitamin C and numerous flavonoids. This study presents the first genome sequencing of this species performed using 10x Genomics and Oxford Nanopore Technology. The draft genome assembly was 519 Mbp in size and consisted of 4,384 contigs, N50 of 597 Kbp, 98.4% BUSCO score, and 37,858 coding sequences. This study also reports the genome-wide phylogeny of this species with 26 other plant species that resolved the phylogenetic position of P. emblica. The presence of three ascorbate biosynthesis pathways including L-galactose, galacturonate, and myo-inositol pathways was confirmed in this genome. A comprehensive comparative evolutionary genomic analysis including gene family expansion/contraction and identification of multiple signatures of adaptive evolution provided evolutionary insights into ascorbate and flavonoid biosynthesis pathways and stone fruit formation through lignin biosynthesis. The availability of this genome will be beneficial for its horticultural, medicinal, dietary, and cosmetic applications and will also help in comparative genomics analysis studies.

The allotetraploid horseradish genome provides insights into subgenome diversification and formation of critical traits

Polyploidization can provide a wealth of genetic variation for adaptive evolution and speciation, but understanding the mechanisms of subgenome evolution as well as its dynamics and ultimate consequences remains elusive. Here, we report the telomere-to-telomere (T2T) gap-free reference genome of allotetraploid horseradish (Armoracia rusticana) sequenced using a comprehensive strategy. The (epi)genomic architecture and 3D chromatin structure of the A and B subgenomes differ significantly, suggesting that both the dynamics of the dominant long terminal repeat retrotransposons and DNA methylation have played critical roles in subgenome diversification. Investigation of the genetic basis of biosynthesis of glucosinolates (GSLs) and horseradish peroxidases reveals both the important role of polyploidization and subgenome differentiation in shaping the key traits. Continuous duplication and divergence of essential genes of GSL biosynthesis (e.g., FMO_GS-OX, IGMT, and GH1 gene family) contribute to the broad GSL profile in horseradish. Overall, the T2T assembly of the allotetraploid horseradish genome expands our understanding of polyploid genome evolution and provides a fundamental genetic resource for breeding and genetic improvement of horseradish.

High-Quality Assembly and Comparative Analysis of Actinidia latifolia and A. valvata Mitogenomes

Kiwifruit (Actinidia) has been recently domesticated as a horticultural crop with remarkably economic and nutritional value. In this study, by combining sequence datasets from Oxford Nanopore long-reads and Illumina short-reads, we de novo assembled two mitogenomes of Actinidia latifolia and A. valvata, respectively. The results indicated that the A. latifolia mitogenome has a single, circular, 825,163 bp molecule while the A. valvata mitogenome possesses two distinct circular molecules, 781,709 and 301,558 bp, respectively. We characterized the genome structure, repeated sequences, DNA transfers, and dN/dS selections. The phylogenetic analyses showed that A. valvata and A. arguta, or A. latifolia and A. eriantha, were clustered together, respectively. This study provides valuable sequence resources for evolutionary study and molecular breeding in kiwifruit.