A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour

Why we thrive beneath a northern sky - genomic signals of selection in apple for adaptation to northern Sweden

Good understanding of the genomic regions underlying adaptation of apple to boreal climates is needed to facilitate efficient breeding of locally adapted apple cultivars. Proper infrastructure for phenotyping and evaluation is essential for identification of traits responsible for adaptation, and dissection of their genetic composition. However, such infrastructure is costly and currently not available for the boreal zone of northern Sweden. Therefore, we used historical pomological data on climate adaptation of 59 apple cultivars and whole genome sequencing to identify genomic regions that have undergone historical selection among apple cultivars recommended for cultivation in northern Sweden. We found the apple collection to be composed of two ancestral groups that are largely concordant with the grouping into 'hardy' and 'not hardy' cultivars based on the pomological literature. Using a number of genome-wide scans for signals of selection, we obtained strong evidence of positive selection at a genomic region around 29 MbHFTH1 of chromosome 1 among apple cultivars in the 'hardy' group. Using phased genotypic data from the 20 K apple Infinium® SNP array, we identified haplotypes associated with the two cultivar groups and traced transmission of these haplotypes through the pedigrees of some apple cultivars. This demonstrates that historical data from pomological literature can be analyzed by population genomic approaches as a step towards revealing the genomic control of a key property for a horticultural niche market. Such knowledge is needed to facilitate efficient breeding strategies for development of locally adapted apple cultivars in the future. The current study illustrates the response to a very strong selective pressure imposed on tree crops by climatic factors, and the importance of genetic research on this topic and feasibility of breeding efforts in the light of the ongoing climate change.

AcMYB96 promotes anthocyanin accumulation in onion (Allium cepa L) without forming the MBW complex

Anthocyanins are major flavonoid compounds with established health benefits. Although the molecular mechanisms of MYB transcription factors (TFs) within the MYB-basic helix-loop-helix (bHLH)-WD-repeat protein (MBW) complex in anthocyanin biosynthesis have been revealed, the functions of other MYB TFs that are unable to form the MBW complex in this process remain unclear. In this study, we uncovered and extensively characterized an R2R3-MYB TF in onion (Allium cepa L.), named AcMYB96, which was identified as a potential anthocyanin activator. AcMYB96 was classified into subgroup 1 of the R2R3-MYB TF family and lacked the conserved sequences required for interactions with bHLH IIIf TFs. Consistently, yeast two-hybrid assays showed that AcMYB96 did not interact with any bHLH IIIf TFs examined, including AcB2 and AtTT8. The transcription pattern of AcMYB96 correlated with the level of anthocyanin accumulation, and its role in activating anthocyanin biosynthesis was confirmed through overexpression in the epithelial cells of onion bulbs and Arabidopsis. Yeast one-hybrid, electrophoretic mobility shift, and promoter transactivation assays further demonstrated that AcMYB96 promoted anthocyanin biosynthesis by binding to the promoters of the chalcone synthase (AcCHS1), anthocyanidin synthase (AcANS), and UDP-glucose-flavonoid 3-O-glucosyltransferase (AcUFGT) genes, thereby activating their expression independent of bHLH IIIf TFs. These results demonstrate that AcMYB96 activates anthocyanin biosynthesis without forming the MBW complex, providing a theoretical foundation to further enrich the gene resources for promoting anthocyanin accumulation and breeding red onions with high anthocyanin content.

Genome-wide identification of starch phosphorylase gene family in Rosa chinensis and expression in response to abiotic stress

Chinese rose (Rosa chinensis) is an important ornamental plant, with economic, cultural, and symbolic significance. During the application of outdoor greening, adverse environments such as high temperature and drought are often encountered, which affect its application scope and ornamental quality. The starch phosphorylase (Pho) gene family participate in the synthesis and decomposition of starch, not only related to plant energy metabolism, but also plays an important role in plant stress resistance. The role of Pho in combating salinity and high temperature stress in R. chinensis remains unknown. In this work, 4 Phos from R. chinensis were detected with Pfam number of Pho (PF00343.23) and predicted by homolog-based prediction (HBP). The Phos are characterized by sequence lengths of 821 to 997 bp, and the proteins are predicted to subcellularly located in the plastid and cytoplasm. The regulatory regions of the Phos contain abundant stress and phytohormone-responsive cis-acting elements. Based on transcriptome analysis, the Phos were found to respond to abiotic stress factors such as drought, salinity, high temperature, and plant phytohormone of jasmonic acid and salicylic acid. The response of Phos to abiotic stress factors such as salinity and high temperature was confirmed by qRT-PCR analysis. To evaluate the genetic characteristics of Phos, a total of 69 Phos from 17 species were analyzed and then classified into 3 groups in phylogenetic tree. The collinearity analysis of Phos in R. chinensis and other species was conducted for the first time. This work provides a view of evolution for the Pho gene family and indicates that Phos play an important role in abiotic stress response of R. chinensis.

A haplotype-resolved genome assembly of Malus domestica 'Red Fuji'

The 'Red Fuji' apple (Malus domestica), is one of the most important and popular economic crops worldwide in the fruit industry. Using PacBio HiFi long reads and Hi-C reads, we assembled a high-quality haplotype-resolved genome of 'Red Fuji', with sizes of 668.7 and 668.8 Mb, and N50 sizes of 34.1 and 31.4 Mb. About 97.2% of sequences were anchored in 34 chromosomes. We annotated both haploid genomes, identifying a total of 95,439 protein-coding genes in the two haplotype genomes, with 98% functional annotation. The haplotype-resolved genome of 'Red Fuji' apple stands as a precise benchmark for an array of analyses, such as comparative genomics, transcriptomics, and allelic expression studies. This comprehensive resource is paramount in unraveling variations in allelic expression, advancing quality improvements, and refining breeding efforts.

A pan-TE map highlights transposable elements underlying domestication and agronomic traits in Asian rice

Transposable elements (TEs) are ubiquitous genomic components and hard to study due to being highly repetitive. Here we assembled 232 chromosome-level genomes based on long-read sequencing data. Coupling the 232 genomes with 15 existing assemblies, we developed a pan-TE map comprising both cultivated and wild Asian rice. We detected 177 084 high-quality TE variations and inferred their derived state using outgroups. We found TEs were one source of phenotypic variation during rice domestication and differentiation. We identified 1246 genes whose expression variation was associated with TEs but not single-nucleotide polymorphisms (SNPs), such as OsRbohB, and validated OsRbohB's relative expression activity using a dual-Luciferase (LUC) reporter assays system. Our pan-TE map allowed us to detect multiple novel loci associated with agronomic traits. Collectively, our findings highlight the contributions of TEs to domestication, differentiation and agronomic traits in rice, and there is massive potential for gene cloning and molecular breeding by the high-quality Asian pan-TE map we generated.

A phased genome of the highly heterozygous 'Texas' almond uncovers patterns of allele-specific expression linked to heterozygous structural variants

The vast majority of traditional almond varieties are self-incompatible, and the level of variability of the species is very high, resulting in a high-heterozygosity genome. Therefore, information on the different haplotypes is particularly relevant to understand the genetic basis of trait variability in this species. However, although reference genomes for several almond varieties exist, none of them is phased and has genome information at the haplotype level. Here, we present a phased assembly of genome of the almond cv. Texas. This new assembly has 13% more assembled sequence than the previous version of the Texas genome and has an increased contiguity, in particular in repetitive regions such as the centromeres. Our analysis shows that the 'Texas' genome has a high degree of heterozygosity, both at SNPs, short indels, and structural variants level. Many of the SVs are the result of heterozygous transposable element insertions, and in many cases, they also contain genic sequences. In addition to the direct consequences of this genic variability on the presence/absence of genes, our results show that variants located close to genes are often associated with allele-specific gene expression, which highlights the importance of heterozygous SVs in almond.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

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

Introduction

Pseudocydonia sinensis, also known as Chinese quince, is a perennial shrub or small tree highly valued for its edibility and medicinal properties.

Method

This study presents the first chromosome-level genome assembly of P. sinensis, achieved using HiFi sequencing and Hi-C scaffolding technology.

Results

The assembly resulted in a high-quality genome of 576.39 Mb in size. The genome was anchored to 17 pseudo-chromosomes, with a contig N50 of 27.6 Mb and a scaffold N50 of 33.8 Mb. Comprehensive assessment using BUSCO, CEGMA and BWA tools indicates the high completeness and accuracy of the genome assembly. Our analysis identified 116 species-specific genes, 1196 expanded genes and 1109 contracted genes. Additionally, the distribution of 4DTv values suggests that the most recent duplication event occurred before the divergence of P. sinensis from both Chaenomeles pinnatifida and Pyrus pyrifolia.

Discussion

The assembly of this high-quality genome provides a valuable platform for the genetic breeding and cultivation of P. sinensis, as well as for the comparison of the genetic complexity of P. sinensis with other important crops in the Rosaceae family.

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

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

Sequence and epigenetic variations of R2R3-MYB transcription factors determine the diversity of taproot skin and flesh colors in different cultivated types of radish (Raphanus sativus L.)

KEY MESSAGE

This study found that three paralogous R2R3-MYB transcription factors exhibit functional divergence among different subspecies and cultivated types in radish. Cultivated radish taproots exhibit a wide range of color variations due to unique anthocyanin accumulation patterns in various tissues. This study investigated the universal principles of taproot color regulation that developed during domestication of different subspecies and cultivated types. The key candidate genes RsMYB1 and RsMYB2, which control anthocyanin accumulation in radish taproots, were identified using bulked segregant analysis in two genetic populations. We introduced the RsMYB1-RsF3'H-RsMYB1Met genetic model to elucidate the complex and unstable genetic regulation of taproot flesh color in Xinlimei radish. Furthermore, we analyzed the expression patterns of three R2R3-MYB transcription factors in lines with different taproot colors and investigated the relationship between RsMYB haplotypes and anthocyanin accumulation in a natural population of 56 germplasms. The results revealed that three paralogous RsMYBs underwent functional divergence during radish domestication, with RsMYB1 regulating the red flesh of Xinlimei radish, and RsMYB2 and RsMYB3 regulating the red skin of East Asian big long radish (R. sativus var. hortensis) and European small radish (R. sativus var. sativus), respectively. Moreover, RsMYB1-H1, RsMYB2-H10, and RsMYB3-H6 were identified as the primary haplotypes exerting regulatory functions on anthocyanin synthesis. These findings provide an understanding of the genetic mechanisms regulating anthocyanin synthesis in radish and offer a potential strategy for early prediction of color variations in breeding programs.

A 224-bp Indel in the Promoter of PeMYB114 Accounts for Anthocyanin Accumulation of Skin in Passion Fruit (Passiflora spp.)

Passion fruit (Passiflora spp.) is an important fruit tree in the family Passifloraceae. The color of the fruit skin, a significant agricultural trait, is determined by the content of anthocyanin in passion fruit. However, the regulatory mechanisms behind the accumulation of anthocyanin in different passion fruit skin colors remain unclear. In the study, we identified and characterized a R2R3-MYB transcription factor, PeMYB114, which functions as a transcriptional activator in anthocyanin biosynthesis. Yeast one-hybrid system and dual-luciferase analysis showed that PeMYB114 could directly activate the expression of anthocyanin structural genes (PeCHS and PeDFR). Furthermore, a natural variation in the promoter region of PeMYB114 alters its expression. PeMYB114purple accessions with the 224-bp insertion have a higher anthocyanin level than PeMYB114yellow accessions with the 224-bp deletion. The findings enhance our understanding of anthocyanin accumulation in fruits and provide genetic resources for genome design for improving passion fruit quality.

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.

Phenomics and transcriptomic profiling of fruit development in distinct apple varieties

Apple is one of the most economically important and popular temperate fruit trees. The domestication of apple has resulted in substantial phenotypic differences, particularly between wild and cultivated varieties. However, the relationship between gene expression and phenotypic variations in apple remains poorly understood. Here, we present a comprehensive dataset featuring five distinct apple varieties, including two wild varieties and three representative cultivated varieties. The dataset comprises of both phenomics data, encompassing twelve fruit quality-related traits continuously measured over two years, and transcriptomic data obtained at different developmental stages with three biological replicates. We performed basic quality control process, gene expression normalization and differential gene expression analysis to demonstrate the utility and reliability of the dataset. Our findings indicate that gene expression strongly related with phenotypic variations in apple. This dataset serves as a valuable resource, encompassing phenomics and transcriptomic data in multiple formats, thereby facilitating further exploration of the relationships between gene expression and phenotypic traits in apple.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

Haplotype-resolved genome of Prunus zhengheensis provides insight into its evolution and low temperature adaptation in apricot

Prunus zhengheensis, an extremely rare population of apricots, originated in warm South-East China and is an excellent material for genetic breeding. However, most apricots and two related species (P. sibirica, P. mandshurica) are found in the cold northern regions in China and the mechanism of their distribution is still unclear. In addition, the classification status of P. zhengheensis is controversial. Thus, we generated a high-quality haplotype-resolved genome for P. zhengheensis, exploring key genetic variations in its adaptation and the causes of phylogenetic incongruence. We found extensive phylogenetic discordances between the nuclear and organelle phylogenies of P. zhengheensis, which could be explained by incomplete lineage sorting. A 242.22-Mb pan-genome of the Armeniaca section was developed with 13 chromosomal genomes. Importantly, we identified a 566-bp insertion in the promoter of the HSFA1d gene in apricot and showed that the activity of the HSFA1d promoter increased under low temperatures. In addition, HSFA1d overexpression in Arabidopsis thaliana indicated that HSFA1d positively regulated plant growth under chilling. Therefore, we hypothesized that the insertion in the promoter of HSFA1d in apricot improved its low-temperature adaptation, allowing it to thrive in relatively cold locations. The findings help explain the weather adaptability of Armeniaca plants.

The structure, function, and evolution of plant centromeres

Centromeres are essential regions of eukaryotic chromosomes responsible for the formation of kinetochore complexes, which connect to spindle microtubules during cell division. Notably, although centromeres maintain a conserved function in chromosome segregation, the underlying DNA sequences are diverse both within and between species and are predominantly repetitive in nature. The repeat content of centromeres includes high-copy tandem repeats (satellites), and/or specific families of transposons. The functional region of the centromere is defined by loading of a specific histone 3 variant (CENH3), which nucleates the kinetochore and shows dynamic regulation. In many plants, the centromeres are composed of satellite repeat arrays that are densely DNA methylated and invaded by centrophilic retrotransposons. In some cases, the retrotransposons become the sites of CENH3 loading. We review the structure of plant centromeres, including monocentric, holocentric, and metapolycentric architectures, which vary in the number and distribution of kinetochore attachment sites along chromosomes. We discuss how variation in CENH3 loading can drive genome elimination during early cell divisions of plant embryogenesis. We review how epigenetic state may influence centromere identity and discuss evolutionary models that seek to explain the paradoxically rapid change of centromere sequences observed across species, including the potential roles of recombination. We outline putative modes of selection that could act within the centromeres, as well as the role of repeats in driving cycles of centromere evolution. Although our primary focus is on plant genomes, we draw comparisons with animal and fungal centromeres to derive a eukaryote-wide perspective of centromere structure and function.

The genome of Haberlea rhodopensis provides insights into the mechanisms for tolerance to multiple extreme environments

Haberlea rhodopensis, a resurrection species, is the only plant known to be able to survive multiple extreme environments, including desiccation, freezing temperatures, and long-term darkness. However, the molecular mechanisms underlying tolerance to these stresses are poorly studied. Here, we present a high-quality genome of Haberlea and found that ~ 23.55% of the 44,306 genes are orphan. Comparative genomics analysis identified 89 significantly expanded gene families, of which 25 were specific to Haberlea. Moreover, we demonstrated that Haberlea preserves its resurrection potential even in prolonged complete darkness. Transcriptome profiling of plants subjected to desiccation, darkness, and low temperatures revealed both common and specific footprints of these stresses, and their combinations. For example, PROTEIN PHOSPHATASE 2C (PP2C) genes were substantially induced in all stress combinations, while PHYTOCHROME INTERACTING FACTOR 1 (PIF1) and GROWTH RESPONSE FACTOR 4 (GRF4) were induced only in darkness. Additionally, 733 genes with unknown functions and three genes encoding transcription factors specific to Haberlea were specifically induced/repressed upon combination of stresses, rendering them attractive targets for future functional studies. The study provides a comprehensive understanding of the genomic architecture and reports details of the mechanisms of multi-stress tolerance of this resurrection species that will aid in developing strategies that allow crops to survive extreme and multiple abiotic stresses.

Near-gapless and haplotype-resolved apple genomes provide insights into the genetic basis of rootstock-induced dwarfing

Dwarfing rootstocks have transformed the production of cultivated apples; however, the genetic basis of rootstock-induced dwarfing remains largely unclear. We have assembled chromosome-level, near-gapless and haplotype-resolved genomes for the popular dwarfing rootstock 'M9', the semi-vigorous rootstock 'MM106' and 'Fuji', one of the most commonly grown apple cultivars. The apple orthologue of auxin response factor 3 (MdARF3) is in the Dw1 region of 'M9', the major locus for rootstock-induced dwarfing. Comparing 'M9' and 'MM106' genomes revealed a 9,723-bp allele-specific long terminal repeat retrotransposon/gypsy insertion, DwTE, located upstream of MdARF3. DwTE is cosegregated with the dwarfing trait in two segregating populations, suggesting its prospective utility in future dwarfing rootstock breeding. In addition, our pipeline discovered mobile mRNAs that may contribute to the development of dwarfed scion architecture. Our research provides valuable genomic resources and applicable methodology, which have the potential to accelerate breeding dwarfing rootstocks for apple and other perennial woody fruit trees.

Chromosomal level genome assemblies of two Malus crabapple cultivars Flame and Royalty

Malus hybrid 'Flame' and Malus hybrid 'Royalty' are representative ornamental crabapples, rich in flavonoids and serving as the preferred materials for studying the coloration mechanism. We generated two sets of high-quality chromosome-level and haplotype-resolved genome of 'Flame' with sizes of 688.2 Mb and 675.7 Mb, and those of 'Royalty' with sizes of 674.1 Mb and 663.6 Mb, all anchored to 17 chromosomes and with a high BUSCO completeness score nearly 99.0%. A total of 47,833 and 47,307 protein-coding genes were annotated in the two haplotype genomes of 'Flame', and the numbers of 'Royalty' were 46,305 and 46,920 individually. The assembled high-quality genomes offer new resources for studying the origin and adaptive evolution of crabapples and the molecular basis of the accumulation of flavonoids and anthocyanins, facilitating molecular breeding of Malus plants.

Genome-Wide Association Study and Prediction of Tassel Weight of Tropical Maize Germplasm in Multi-Parent Population

Tassel weight (TW) is a crucial agronomic trait that significantly affects pollen supply and grain yield development in maize breeding. To improve maize yield and develop new varieties, a comprehensive understanding of the genetic mechanisms underlying tassel weight is essential. In this study, tropical maize inbred lines, namely CML312, CML373, CML444, and YML46, were selected as female parents and crossed with the elite maize inbred line Ye107, which served as the common male parent, to develop a multi-parent population comprising four F8 recombinant inbred line (RIL) subpopulations. Using 6616 high-quality single nucleotide polymorphism (SNP) markers, we conducted genome-wide association analysis (GWAS) and genomic selection (GS) on 642 F8 RILs in four subpopulations across three different environments. Through GWAS, we identified 16 SNPs that were significantly associated with TW, encompassing two stable loci expressed across multiple environments. Furthermore, within the candidate regions of these SNPs, we discovered four novel candidate genes related to TW, namely Zm00001d044362, Zm00001d011048, Zm00001d011049, and Zm00001d031173 distributed on chromosomes 1, 3, and 8, which have not been previously reported. These genes are involved in processes such as signal transduction, growth and development, protein splicing, and pollen development, all of which play crucial roles in inflorescence meristem development, directly affecting TW. The co-localized SNP, S8_137379725, on chromosome 8 was situated within a 16.569 kb long terminal repeat retrotransposon (LTR-RT), located 22.819 kb upstream and 26.428 kb downstream of the candidate genes (Zm00001d011048 and Zm00001d011049). When comparing three distinct GS models, the BayesB model demonstrated the highest accuracy in predicting TW. This study establishes the theoretical foundation for future research into the genetic mechanisms underlying maize TW and the efficient breeding of high-yielding varieties with desired tassel weight through GS.

Genome-wide identification, classification and expression analysis of MYB gene family in coconut (Cocos nucifera L.)

MYB transcription factors regulate the growth, development, and secondary metabolism of plant species. To investigate the origin of color variations in coconut pericarp, we identified and analyzed the MYB gene family present in coconut. According to the sequence of MYB genes in Arabidopsis thaliana, homologous MYB gene sequences were found in the whole genome database of coconut, the conserved sequence motifs within MYB proteins were analyzed by Motif Elicitation (MEME) tool, and the sequences without conservative structure were eliminated. Additionally, we employed RNA-seq technology to generate gene expression signatures of the R2R3-MYB genes across distinctive coconut parts exhibiting diverse colors. To validate these profiles, we conducted quantitative PCR (qPCR). Through comprehensive genome-wide screening, we successfully identified a collection of 179 MYB genes in coconut. Subsequent phylogenetic analysis categorized these 179 coconut MYB genes into 4-subfamilies: 124 R2R3-MYB, 4 3R-MYB types, 4 4R-MYB type, and 47 unknown types. Furthermore, these genes were further divided into 34 subgroups, with 28 of these subgroups successfully classified into known subfamilies found in Arabidopsis thaliana. By mapping the CnMYB genes onto the 16 chromosomes of the coconut genome, we unveiled a collinearity association between them. Moreover, a preservation of gene structure and motif distribution was observed across the CnMYB genes. Our research encompassed a thorough investigation of the R2R3-MYB genes present in the coconut genome, including the chromosomal localization, gene assembly, conserved regions, phylogenetic associations, and promoter cis-acting elements of the studied genes. Our findings revealed a collection of 12 R2R3-MYB candidate genes, namely CnMYB8, CnMYB15, CnMYB27, CnMYB28, CnMYB61, CnMYB63, CnMYB68, CnMYB94, CnMYB101, CnMYB150, CnMYB153, and CnMYB164. These genes showed differential expressions in diverse tissues and developmental stages of four coconut species, such as CnMYB68, CnMYB101, and CnMYB28 exhibited high expression in majority of tissues and coconut species, while CnMYB94 and CnMYB164 showed lower expression. These findings shed light on the crucial functional divergence of CnMYB genes across various coconut tissues, suggesting these genes as promising candidate genes for facilitating color development in this important crop.

Effect of structural variation in the promoter region of RsMYB1.1 on the skin color of radish taproot

Accumulation of anthocyanins in the taproot of radish is an agronomic trait beneficial for human health. Several genetic loci are related to a red skin or flesh color of radish, however, the functional divergence of candidate genes between non-red and red radishes has not been investigated. Here, we report that a novel genetic locus on the R2 chromosome, where RsMYB1.1 is located, is associated with the red color of the skin of radish taproot. A genome-wide association study (GWAS) of 66 non-red-skinned (nR) and 34 red-skinned (R) radish accessions identified three nonsynonymous single nucleotide polymorphisms (SNPs) in the third exon of RsMYB1.1. Although the genotypes of SNP loci differed between the nR and R radishes, no functional difference in the RsMYB1.1 proteins of nR and R radishes in their physical interaction with RsTT8 was detected by yeast-two hybrid assay or in anthocyanin accumulation in tobacco and radish leaves coexpressing RsMYB1.1 and RsTT8. By contrast, insertion- or deletion-based GWAS revealed that one large AT-rich low-complexity sequence of 1.3-2 kb was inserted in the promoter region of RsMYB1.1 in the nR radishes (RsMYB1.1nR), whereas the R radishes had no such insertion; this represents a presence/absence variation (PAV). This insertion sequence (RsIS) was radish specific and distributed among the nine chromosomes of Raphanus genomes. Despite the extremely low transcription level of RsMYB1.1nR in the nR radishes, the inactive RsMYB1.1nR promoter could be functionally restored by deletion of the RsIS. The results of a transient expression assay using radish root sections suggested that the RsIS negatively regulates the expression of RsMYB1.1nR, resulting in the downregulation of anthocyanin biosynthesis genes, including RsCHS, RsDFR, and RsANS, in the nR radishes. This work provides the first evidence of the involvement of PAV in an agronomic trait of radish.

Gapless genome assembly and epigenetic profiles reveal gene regulation of whole-genome triplication in lettuce

BACKGROUND

Lettuce, an important member of the Asteraceae family, is a globally cultivated cash vegetable crop. With a highly complex genome (∼2.5 Gb; 2n = 18) rich in repeat sequences, current lettuce reference genomes exhibit thousands of gaps, impeding a comprehensive understanding of the lettuce genome.

FINDINGS

Here, we present a near-complete gapless reference genome for cutting lettuce with high transformability, using long-read PacBio HiFi and Nanopore sequencing data. In comparison to stem lettuce genome, we identify 127,681 structural variations (SVs, present in 0.41 Gb of sequence), reflecting the divergence of leafy and stem lettuce. Interestingly, these SVs are related to transposons and DNA methylation states. Furthermore, we identify 4,612 whole-genome triplication genes exhibiting high expression levels associated with low DNA methylation levels and high N6-methyladenosine RNA modifications. DNA methylation changes are also associated with activation of genes involved in callus formation.

CONCLUSIONS

Our gapless lettuce genome assembly, an unprecedented achievement in the Asteraceae family, establishes a solid foundation for functional genomics, epigenomics, and crop breeding and sheds new light on understanding the complexity of gene regulation associated with the dynamics of DNA and RNA epigenetics in genome evolution.

Genome variation and LTR-RT analyses of an ancient peach landrace reveal mechanism of blood-flesh fruit color formation and fruit maturity date advancement

Peach (Prunus persica) landrace has typical regional characteristics, strong environmental adaptability, and contains many valuable genes that provide the foundation for breeding excellent varieties. Therefore, it is necessary to assemble the genomes of specific landraces to facilitate the localization and utilization of these genes. Here, we de novo assembled a high-quality genome from an ancient blood-fleshed Chinese landrace Tianjin ShuiMi (TJSM) that originated from the China North Plain. The assembled genome size was 243.5 Mb with a contig N50 of 23.7 Mb and a scaffold N50 of 28.6 Mb. Compared with the reported peach genomes, our assembled TJSM genome had the largest number of specific structural variants (SVs) and long terminal repeat-retrotransposons (LTR-RTs). Among the LTR-RTs with the potential to regulate their host genes, we identified a 6688 bp LTR-RT (named it blood TE) in the promoter of NAC transcription factor-encoding PpBL, a gene regulating peach blood-flesh formation. The blood TE was not only co-separated with the blood-flesh phenotype but also associated with fruit maturity date advancement and different intensities of blood-flesh color formation. Our findings provide new insights into the mechanism underlying the development of the blood-flesh color and determination of fruit maturity date and highlight the potential of the TJSM genome to mine more variations related to agronomic traits in peach fruit.

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.

Chromosome-level phased genome assembly of "Antonovka" identified candidate apple scab-resistance genes highly homologous to HcrVf2 and HcrVf1 on linkage group 1

Apple scab, a fungal disease caused by Venturia inaequalis, leads to losses in both yield and fruit quality of apples (Malus domestica Borkh.). Most commercial apple cultivars, including those containing the well-characterized Rvi6-scab-resistance locus on linkage group (LG) 1, are susceptible to scab. HcrVf2 and HcrVf1 are considered the main paralogs of the Rvi6 locus. The major apple scab-resistance loci Vhc1 in "Honeycrisp" and Rvi17 in "Antonovka," were identified in close proximity to HcrVf2. In this study, we used long-read sequencing and in silico gene sequence characterization to identify candidate resistance genes homologous to HcrVf2 and HcrVf1 in Honeycrisp and Antonovka. Previously published chromosome-scale phased assembly of Honeycrisp and a newly assembled phased genome of Antonovka 172670-B were used to identify HcrVf2 and HcrVf1 homologs spanning Vhc1 and Rvi17 loci. In combination with 8 available Malus assemblies, 43 and 46 DNA sequences highly homologous to HcrVf2 and HcrVf1, respectively, were identified on LG 1 and 6, with identity and coverage ranging between 87-95 and 81-95%, respectively. Among these homologs, 2 candidate genes in Antonovka and Honeycrisp haplome A are located in close physical proximity to the scab-resistance marker Ch-Vf1 on LG 1. They showed the highest identity and coverage (95%) of HcrVf2 and only minor changes in the protein motifs. They were identical by state between each other, but not with HcrVf2. This study offers novel genomic resources and insights into the Vhc1 and Rvi17 loci on LG 1 and identifies candidate genes for further resistance characterization.

Dihydrochalcone glycoside biosynthesis in Malus is regulated by two MYB-like transcription factors and is required for seed development

Dihydrochalcones (DHCs) including phlorizin (phloretin 2'-O-glucoside) and its positional isomer trilobatin (phloretin 4'-O-glucoside) are the most abundant phenylpropanoids in apple (Malus spp.). Transcriptional regulation of DHC production is poorly understood despite their importance in insect- and pathogen-plant interactions in human physiology research and in pharmaceuticals. In this study, segregation in hybrid populations and bulked segregant analysis showed that the synthesis of phlorizin and trilobatin in Malus leaves are both single-gene-controlled traits. Promoter sequences of PGT1 and PGT2, two glycosyltransferase genes involved in DHC glycoside synthesis, were shown to discriminate Malus with different DHC glycoside patterns. Differential PGT1 and PGT2 promoter activities determined DHC glycoside accumulation patterns between genotypes. Two transcription factors containing MYB-like DNA-binding domains were then shown to control DHC glycoside patterns in different tissues, with PRR2L mainly expressed in leaf, fruit, flower, stem, and seed while MYB8L mainly expressed in stem and root. Further hybridizations between specific genotypes demonstrated an absolute requirement for DHC glycoside production in Malus during seed development which explains why no Malus spp. with a null DHC chemotype have been reported.

De novo reconstruction of satellite repeat units from sequence data

Satellite DNA are long tandemly repeating sequences in a genome and may be organized as high-order repeats (HORs). They are enriched in centromeres and are challenging to assemble. Existing algorithms for identifying satellite repeats either require the complete assembly of satellites or only work for simple repeat structures without HORs. Here we describe Satellite Repeat Finder (SRF), a new algorithm for reconstructing satellite repeat units and HORs from accurate reads or assemblies without prior knowledge on repeat structures. Applying SRF to real sequence data, we show that SRF could reconstruct known satellites in human and well-studied model organisms. We also find satellite repeats are pervasive in various other species, accounting for up to 12% of their genome contents but are often underrepresented in assemblies. With the rapid progress in genome sequencing, SRF will help the annotation of new genomes and the study of satellite DNA evolution even if such repeats are not fully assembled.

Comparative Analysis of Transposable Elements in Strawberry Genomes of Different Ploidy Levels

Transposable elements (TEs) make up a large portion of plant genomes and play a vital role in genome structure, function, and evolution. Cultivated strawberry (Fragaria x ananassa) is one of the most important fruit crops, and its octoploid genome was formed through several rounds of genome duplications from diploid ancestors. Here, we built a pan-genome TE library for the Fragaria genus using ten published strawberry genomes at different ploidy levels, including seven diploids, one tetraploid, and two octoploids, and performed comparative analysis of TE content in these genomes. The TEs comprise 51.83% (F. viridis) to 60.07% (F. nilgerrensis) of the genomes. Long terminal repeat retrotransposons (LTR-RTs) are the predominant TE type in the Fragaria genomes (20.16% to 34.94%), particularly in F. iinumae (34.94%). Estimating TE content and LTR-RT insertion times revealed that species-specific TEs have shaped each strawberry genome. Additionally, the copy number of different LTR-RT families inserted in the last one million years reflects the genetic distance between Fragaria species. Comparing cultivated strawberry subgenomes to extant diploid ancestors showed that F. vesca and F. iinumae are likely the diploid ancestors of the cultivated strawberry, but not F. viridis. These findings provide new insights into the TE variations in the strawberry genomes and their roles in strawberry genome evolution.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.

A haplotype resolved chromosome-scale assembly of North American wild apple Malus fusca and comparative genomics of the fire blight Mfu10 locus

The Pacific crabapple (Malus fusca) is a wild relative of the commercial apple (Malus × domestica). With a range extending from Alaska to Northern California, M. fusca is extremely hardy and disease resistant. The species represents an untapped genetic resource for the development of new apple cultivars with enhanced stress resistance. However, gene discovery and utilization of M. fusca have been hampered by the lack of genomic resources. Here, we present a high‐quality, haplotype‐resolved, chromosome‐scale genome assembly and annotation for M. fusca. The genome was assembled using high‐fidelity long‐reads and scaffolded using genetic maps and high‐throughput chromatin conformation capture sequencing, resulting in one of the most contiguous apple genomes to date. We annotated the genome using public transcriptomic data from the same species taken from diverse plant structures and developmental stages. Using this assembly, we explored haplotypic structural variation within the genome of M. fusca, identifying thousands of large variants. We further showed high sequence co‐linearity with other domesticated and wild Malus species. Finally, we resolve a known quantitative trait locus associated with resistance to fire blight (Erwinia amylovora). Insights gained from the assembly of a reference‐quality genome of this hardy wild apple relative will be invaluable as a tool to facilitate DNA‐informed introgression breeding.

MdAIL5 overexpression promotes apple adventitious shoot regeneration by regulating hormone signaling and activating the expression of shoot development-related genes

Adventitious shoot (AS) regeneration is a significant factor in the genetic transformation of horticultural plants. It is also a noteworthy approach to their vegetative propagation. AS regeneration remains highly dependent on the genotype or maturity of explants. We here found that the AS regeneration abilities of apple leaves were positively correlated with MdAIL5 expression. MdAIL5 overexpression dramatically increased AS regeneration efficiency. Notably, MdAIL5 overexpression could restore the AS formation ability of explants to a certain extent, which was lost with an increase in maturity. Endogenous hormone detection revealed that MdAIL5 overexpression changed the contents of auxin, cytokinin (CK), and other hormones in apple leaves. Transcriptome analysis revealed that many genes related to auxin, CK, and brassinolide signaling pathways were significantly and differentially expressed between MdAIL5-overexpressing transgenic apple and wild-type apple plants. Yeast one-hybrid assays, the electrophoretic mobility shift assay, and the dual-luciferase reporter assay revealed that MdAIL5 directly binds to MdARF9 and MdHB14 promoters and positively affects their expression. We here established a model of MdAIL5 regulating AS formation, which acts as a theoretical basis for facilitating genotype- or explant maturity-independent AS regeneration in the future.

Rvi4 and Rvi15 are the same apple scab resistance genes

The apple (Malus x domestica) scab (Venturia inaequalis) resistance genes Rvi4 and Rvi15 were mapped to a similar region on the top of linkage group 2 and both resistance genes elicit the same type of resistance reaction, i.e., a hypersensitive response; hence, it is suspected that the two genes may be the same. As the two resistance genes Rvi4 and Rvi15 are currently used in apple breeding, it is important to clarify whether the two resistance genes are the same or not. Several approaches were used to make this determination. First, the pedigree of the genotype GMAL 2473, the source of Rvi15, was reconstructed. GMAL 2473 was found to be an F1 of 'Russian seedling', the genotype, which is known to also be the source of Rvi4. Next, it was further demonstrated that 'Regia', a cultivar known to carry Rvi4 (and Rvi2), carries the same gene (Vr2-C), which was demonstrated to be the gene inducing Rvi15 resistance. Finally, it was shown that transgenic lines carrying Vr2-C are compatible with race 4 apple scab isolates. Taken all together, these results definitively demonstrate that Rvi4 and Rvi15 are the same resistance gene. For future studies, we suggest referring to this resistance with the first name that was assigned to this gene, namely Rvi4.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01421-0.

Time course of changes in the transcriptome during russet induction in apple fruit

BACKGROUND

Russeting is a major problem in many fruit crops. Russeting is caused by environmental factors such as wounding or moisture exposure of the fruit surface. Despite extensive research, the molecular sequence that triggers russet initiation remains unclear. Here, we present high-resolution transcriptomic data by controlled russet induction at very early stages of fruit development. During Phase I, a patch of the fruit surface is exposed to surface moisture. For Phase II, moisture exposure is terminated, and the formerly exposed surface remains dry. We targeted differentially expressed transcripts as soon as 24 h after russet induction.

RESULTS

During moisture exposure (Phase I) of 'Pinova' apple, transcripts associated with the cell cycle, cell wall, and cuticle synthesis (SHN3) decrease, while those related to abiotic stress increase. NAC35 and MYB17 were the earliest induced genes during Phase I. They are therefore linked to the initial processes of cuticle microcracking. After moisture removal (Phase II), the expression of genes related to meristematic activity increased (WOX4 within 24 h, MYB84 within 48 h). Genes related to lignin synthesis (MYB52) and suberin synthesis (MYB93, WRKY56) were upregulated within 3 d after moisture removal. WOX4 and AP2B3 are the earliest differentially expressed genes induced in Phase II. They are therefore linked to early events in periderm formation. The expression profiles were consistent between two different seasons and mirrored differences in russet susceptibility in a comparison of cultivars. Furthermore, expression profiles during Phase II of moisture induction were largely identical to those following wounding.

CONCLUSIONS

The combination of a unique controlled russet induction technique with high-resolution transcriptomic data allowed for the very first time to analyse the formation of cuticular microcracks and periderm in apple fruit immediately after the onset of triggering factors. This data provides valuable insights into the spatial-temporal dynamics of russeting, including the synthesis of cuticles, dedifferentiation of cells, and impregnation of cell walls with suberin and lignin.

Stone cell formation in the pedicel of pears and apples

MAIN CONCLUSION

For the first time, stone cells in pear and apple pedicel were studied. The lignification of the pedicel outer part was correlated with flesh, and the secondary cell wall biosynthesis genes were activated. Fruit pedicels act as bridges between the fruit and the shoot. They have secondary thickened cell walls that presumably function in mechanical support, water and nutrient transport. Stone cells are cells with a secondary cell wall thickening. In pears, yet not in apples, the stone cells affect the flesh texture. There have been few reports on stone cell formation in pear and apple pedicels; therefore, we studied these cells for the first time. The apple pedicel had few stone cells in the cortex. The formation of stone cells in pear continued until seven weeks after flowering (WAF), and the density was significantly higher than in apple. The stone cell formation degree (SFD) of pear was 3.6-7.1 times higher than that of apple. Total lignin and lignin non-condensed structure (G and S units) content in the pear pedicle outer part was 1.5-2.7 times higher than that of the apple at harvest. The SFD of the pedicel outer part had a positive correlation with the G and S units content of the flesh. The total lignin and G and S units content between flesh and the pedicel outer part were positively correlated. Correlation analysis revealed a positive relationship between fruit and pedicel formation of the stone cells. The WGCNA showed that NST3 was linked to NAC028, MYB46, CESA, POD, LAC, and VSR6. These genes were highly expressed in the outer part of the pear pedicel, while they were suppressed in that issue of the apple at 4 WAF.

Transcription factor NnMYB5 controls petal color by regulating GLUTATHIONE S-TRANSFERASE2 in Nelumbo nucifera

Lotus (Nelumbo spp.) is an important aquatic ornamental genus in the family Nelumbonaceae comprising only 2 species: Nelumbo lutea with yellow flowers and Nelumbo nucifera with red or white flowers. The petal color variations between these 2 species have previously been associated with the potential activities of FLAVONOL SYNTHASE (FLS) and MYB5. However, the underlying genetic mechanisms of flower color divergence within the N. nucifera species remain unclear. Here, quantitative trait locus mapping led to the identification of MYB5, a candidate gene controlling petal color in N. nucifera. Genotyping of 213 natural lotus accessions revealed an 80 kb presence/absence variant (PAV) of the NnMYB5 gene that is associated with petal color variation. Transcriptome analysis, dual-luciferase, and yeast 1-hybrid assays showed that NnMYB5 could directly activate the anthocyanin transporter gene GLUTATHIONE S-TRANSFERASE2 (NnGST2). Heterologous expression of NnGST2 in Arabidopsis (Arabidopsis thaliana) and its overexpression in lotus petals induced anthocyanin accumulation. Deletion of the 80 kb PAV within NnMYB5 inactivated NnGST2 expression and blocked anthocyanin accumulation in white N. nucifera petals. In contrast, the anthocyanin deficiency of N. lutea occurred due to pseudogenized NlMYB5 alleles. Our results establish a regulatory link between NnMYB5 and NnGST2 in petal anthocyanin accumulation and demonstrate the independent mechanisms controlling flower coloration in Nelumbo.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Low temperature-induced regulatory network rewiring via WRKY regulators during banana peel browning

Banana (Musa spp.) fruits, as typical tropical fruits, are cold sensitive, and lower temperatures can disrupt cellular compartmentalization and lead to severe browning. How tropical fruits respond to low temperature compared to the cold response mechanisms of model plants remains unknown. Here, we systematically characterized the changes in chromatin accessibility, histone modifications, distal cis-regulatory elements, transcription factor binding, and gene expression levels in banana peels in response to low temperature. Dynamic patterns of cold-induced transcripts were generally accompanied by concordant chromatin accessibility and histone modification changes. These upregulated genes were enriched for WRKY binding sites in their promoters and/or active enhancers. Compared to banana peel at room temperature, large amounts of banana WRKYs were specifically induced by cold and mediated enhancer-promoter interactions regulating critical browning pathways, including phospholipid degradation, oxidation, and cold tolerance. This hypothesis was supported by DNA affinity purification sequencing, luciferase reporter assays, and transient expression assay. Together, our findings highlight widespread transcriptional reprogramming via WRKYs during banana peel browning at low temperature and provide an extensive resource for studying gene regulation in tropical plants in response to cold stress, as well as potential targets for improving cold tolerance and shelf life of tropical fruits.

Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis

Somatic embryogenesis (SE), like zygotic embryo development, is a progressive process. Early SE is the beginning of a switch from a somatic to an embryogenic state and is an important stage for initiating chromatin reprogramming of SE. Previous studies suggest that changes in chromatin accessibility occur during early SE, although information on the 3D structure of chromatin is not yet available. Here, we present a chromosome-level genome assembly of longan (Dimocarpus longan) using PacBio combined with high-through chromosome conformation capture scaffolding, which resulted in a 446 Mb genome assembly anchored onto 15 scaffolds. During early SE, chromatin was concentrated and then decondensed, and a large number of long terminal repeat retrotransposons (LTR-RTs) were enriched in the local chromatin interaction region, suggesting LTR-RTs were involved in chromatin reorganization. Early SE was accompanied by the transformation from A to B compartments, and the interactions between B compartments were enhanced. Results from chromatin accessibility, monomethylation of histone H3 at lysine 4 (H3K4me1) modification, and transcription analyses further revealed a gene regulatory network for cell wall thickening during SE. Particularly, we found that the H3K4me1 differential peak binding motif showed abnormal activation of ethylene response factor transcription factors and participation in SE. The chromosome-level genomic and multiomics analyses revealed the 3D conformation of chromatin during early SE, providing insight into the molecular mechanisms underlying cell wall thickening and the potential regulatory networks of TFs during early SE in D. longan. These results provide additional clues for revealing the molecular mechanisms of plant SE.

Non-precious metal-based heterostructure catalysts for hydrogen evolution reaction: mechanisms, design principles, and future prospects

As a highly promising clean energy source to replace fossil fuels in the 21st century, hydrogen energy has garnered considerable attention, with water electrolysis emerging as a key hydrogen production technology. The development of highly active and stable non-precious metal-based catalysts for the hydrogen evolution reaction (HER) is crucial for achieving efficient and low-cost hydrogen production through electrolysis. Recently, heterostructure composite catalysts comprising two or more non-precious metals have demonstrated outstanding catalytic performance. First, we introduced the basic mechanism of the HER and, based on the reported HER theory, discussed the essence of constructing heterostructures to improve the catalytic activity of non-noble metal-based catalysts, that is, the coupling effect between components effectively regulates the electronic structure and the position of d-band centers. Then three catalytic effects of non-precious metal-based heterogeneous catalysts are described: synergistic effect, electron transfer effect and support effect. Lastly, we emphasized the potential of non-precious metal-based heterogeneous catalysts to replace precious metal-based catalysts, and summarized the future prospects and challenges.

The Synergistic Effects of Environmental and Genetic Factors on the Regulation of Anthocyanin Accumulation in Plant Tissues

Anthocyanin accumulation is responsible for the coloration of apple fruit, and their accumulation depends on the expression of anthocyanin biosynthesis-related genes. Light is an environmental stimulus that induces fruit color by regulating genes involved in the anthocyanin biosynthesis pathway. In this study, the roles of light and genetic factors on fruit coloration and anthocyanin accumulation in apple fruit were investigated. Three genes in the anthocyanin biosynthesis pathway, MdCHS, MdANS, and MdUFGT1, were synthesized and cloned into a viral-based expression vector system for transient expression in 'Ruby S' apple fruits. Apple fruits were agroinfiltrated with expression vectors harboring MdCHS, MdANS, and MdUFGT1. Agroinfiltrated apple fruits were then either kept in the dark (bagged fruits) or exposed to light (exposed fruits). The agroinfiltrated fruits showed significantly different coloration patterns, transcript expression levels, and anthocyanin accumulation compared to the control fruits. Moreover, these parameters were higher in exposed fruits than in bagged fruits. For stable expression, MdCHS was introduced into a binary vector under the control of the rice α-amylase 3D (RAmy3D) promoter. The ectopic overexpression of MdCHS in transgenic rice calli showed a high accumulation of anthocyanin content. Taken together, our findings suggest that light, together with the overexpression of anthocyanin biosynthesis genes, induced the coloration and accumulation of anthocyanin content in apple fruits by upregulating the expression of the genes involved in the anthocyanin biosynthesis pathway.

Transcriptomic approach to uncover dynamic events in the development of mid-season sunburn in apple fruit

Apples grown in high heat, high light, and low humidity environments are at risk for sun injury disorders like sunburn and associated crop losses. Understanding the physiological and molecular mechanisms underlying sunburn will support improvement of mitigation strategies and breeding for more resilient varieties. Numerous studies have highlighted key biochemical processes involved in sun injury, such as the phenylpropanoid and reactive oxygen species (ROS) pathways, demonstrating both enzyme activities and expression of related genes in response to sunburn conditions. Most previous studies have focused on at-harvest activity of a small number of genes in response to heat stress. Thus, it remains unclear how stress events earlier in the season affect physiology and gene expression. Here, we applied heat stress to mid-season apples in the field and collected tissue along a time course-24, 48, and 72 h following a heat stimulus-to investigate dynamic gene expression changes using a transcriptomic lens. We found a relatively small number of differentially expressed genes (DEGs) and enriched functional terms in response to heat treatments. Only a few of these belonged to pathways previously described to be involved in sunburn, such as the AsA-GSH pathway, while most DEGs had not yet been implicated in sunburn or heat stress in pome fruit.

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

Haplotype-resolved genomes of wild octoploid progenitors illuminate genomic diversifications from wild relatives to cultivated strawberry

Strawberry is an emerging model for studying polyploid genome evolution and rapid domestication of fruit crops. Here we report haplotype-resolved genomes of two wild octoploids (Fragaria chiloensis and Fragaria virginiana), the progenitor species of cultivated strawberry. Substantial variation is identified between species and between haplotypes. We redefine the four subgenomes and track the genetic contributions of diploid species by additional sequencing of the diploid F. nipponica genome. We provide multiple lines of evidence that F. vesca and F. iinumae, rather than other described extant species, are the closest living relatives of these wild and cultivated octoploids. In response to coexistence with quadruplicate gene copies, the octoploid strawberries have experienced subgenome dominance, homoeologous exchanges and coordinated expression of homoeologous genes. However, some homoeologues have substantially altered expression bias after speciation and during domestication. These findings enhance our understanding of the origin, genome evolution and domestication of strawberries.

Small RNA and Degradome Sequencing Reveal Roles of miRNAs in the Petal Color Fading of Malus Crabapple

The Malus crabapple is an important woody ornamental plant. The fading of petals during its development significantly affects their ornamental value. Petal color is related to anthocyanin content and miRNAs play an important role in the post-transcriptional regulation of anthocyanin synthesis. However, the mechanisms underlying miRNA regulation of petal fading have rarely been studied. Transcriptome and small RNA sequencing of petals from the blooming phases of Malus. 'Indian Summer' varieties S1 (small bud), S2 (initial-flowering), and S3 (late-flowering) allowed us to identify 230 known miRNAs and 17 novel miRNAs, including 52 differentially expressed miRNAs which targeted 494 genes and formed 823 miRNA-target pairs. Based on the target gene annotation results, miRNA-target pairs were screened that may be involved in the fading process of Malus crabapple petals through three different pathways: anthocyanin synthesis, transport, and degradation, involving mcr-miR858-MYB1\MYB5 and mcr-miR396-McCHI inhibiting anthocyanin synthesis; mcr-miR167, mcr-miR390, mcr-miR535, and mcr-miR858 inhibiting anthocyanin transport from the cytoplasm to the vacuole by targeting ABC transporter genes (ABCB, ABCC, ABCD, and ABCG); and mcr-miR398 targeting the superoxide dismutase genes (CZSOD2 and CCS) to accelerate anthocyanin degradation. These findings offer a novel approach to understanding the mechanism of petal fading and serve as a reference for other plants with floral fading.

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

Efficient homology-based annotation of transposable elements using minimizers

Premise

Transposable elements (TEs) make up more than half of the genomes of complex plant species and can modulate the expression of neighboring genes, producing significant variability of agronomically relevant traits. The availability of long-read sequencing technologies allows the building of genome assemblies for plant species with large and complex genomes. Unfortunately, TE annotation currently represents a bottleneck in the annotation of genome assemblies.

Methods and Results

We present a new functionality of the Next-Generation Sequencing Experience Platform (NGSEP) to perform efficient homology-based TE annotation. Sequences in a reference library are treated as long reads and mapped to an input genome assembly. A hierarchical annotation is then assigned by homology using the annotation of the reference library. We tested the performance of our algorithm on genome assemblies of different plant species, including Arabidopsis thaliana, Oryza sativa, Coffea humblotiana, and Triticum aestivum (bread wheat). Our algorithm outperforms traditional homology-based annotation tools in speed by a factor of three to >20, reducing the annotation time of the T. aestivum genome from months to hours, and recovering up to 80% of TEs annotated with RepeatMasker with a precision of up to 0.95.

Conclusions

NGSEP allows rapid analysis of TEs, especially in very large and TE-rich plant genomes.

Telomere-to-telomere carrot (Daucus carota) genome assembly reveals carotenoid characteristics

Carrot (Daucus carota) is an Apiaceae plant with multi-colored fleshy roots that provides a model system for carotenoid research. In this study, we assembled a 430.40 Mb high-quality gapless genome to the telomere-to-telomere (T2T) level of "Kurodagosun" carrot. In total, 36 268 genes were identified and 34 961 of them were functionally annotated. The proportion of repeat sequences in the genome was 55.3%, mainly long terminal repeats. Depending on the coverage of the repeats, 14 telomeres and 9 centromeric regions on the chromosomes were predicted. A phylogenetic analysis showed that carrots evolved early in the family Apiaceae. Based on the T2T genome, we reconstructed the carotenoid metabolic pathway and identified the structural genes that regulate carotenoid biosynthesis. Among the 65 genes that were screened, 9 were newly identified. Additionally, some gene sequences overlapped with transposons, suggesting replication and functional differentiation of carotenoid-related genes during carrot evolution. Given that some gene copies were barely expressed during development, they might be functionally redundant. Comparison of 24 cytochrome P450 genes associated with carotenoid biosynthesis revealed the tandem or proximal duplication resulting in expansion of CYP gene family. These results provided molecular information for carrot carotenoid accumulation and contributed to a new genetic resource.

Chromosome-level genome assembly of Murraya paniculata sheds light on biosynthesis of floral volatiles

BACKGROUND

Murraya paniculata (L.) Jack, commonly called orange jessamine in the family Rutaceae, is an important ornamental plant in tropical and subtropical regions which is famous for its strong fragrance. Although genome assemblies have been reported for many Rutaceae species, mainly in the genus Citrus, full genomic information has not been reported for M. paniculata, which is a prerequisite for in-depth genetic studies on Murraya and manipulation using genetic engineering techniques. Here, we report a high-quality chromosome-level genome assembly of M. paniculata and aim to provide insights on the molecular mechanisms of flower volatile biosynthesis.

RESULTS

The genome assembly with a contig N50 of 18.25 Mb consists of 9 pseudomolecules and has a total length of 216.86 Mb. Phylogenetic analysis revealed that M. paniculata diverged from the common ancestor approximately 25 million years ago and has not undergone any species-specific whole genome duplication events. Genome structural annotation and comparative genomics analysis revealed that there are obvious differences in transposon contents among the genomes of M. paniculata and Citrus species, especially in the upstream regions of genes. Research on the flower volatiles of M. paniculata and C. maxima at three flowering stages revealed significant differences in volatile composition with the flowers of C. maxima lacking benzaldehyde and phenylacetaldehyde. Notably, there are transposons inserted in the upstream region of the phenylacetaldehyde synthase (PAAS) genes Cg1g029630 and Cg1g029640 in C. maxima, but not in the upstream region of three PAAS genes Me2G_2379, Me2G_2381, and Me2G_2382 in M. paniculata. Our results indicated that compared to the low expression levels of PAAS genes in C. maxima, the higher expression levels of the three PAAS genes in M. paniculata are the main factor affecting the phenylacetaldehyde biosynthesis and causing the content difference of phenylacetaldehyde. The phenylacetaldehyde synthetic activities of the enzymes encoded by M. paniculata PAAS genes were validated by in vitro analyses.

CONCLUSIONS

Our study provides useful genomic resources of M. paniculata for further research on Rutaceae plants, identifies new PAAS genes, and provides insights into how transposons contribute to variations in flower volatiles among Murraya and Citrus plants.

Chromosome-scale assembly and quantitative trait locus mapping for major economic traits of the Culter alburnus genome using Illumina and PacBio sequencing with Hi-C mapping information

Topmouth culter (Culter alburnus) is an economically important freshwater fish with high nutritional value. However, its potential genetic advantages have not been fully exploited. Therefore, we aimed to determine the genome sequence of C. alburnus and examine quantitative trait loci (QTLs) related to major economic traits. The results showed that 24 pseudochromosomes were anchored by 914.74 Mb of the C. alburnus genome sequence. De novo sequencing identified 31,279 protein-coding genes with an average length of 8507 bp and average coding sequ ence of 1115 bp. In addition, a high-density genetic linkage map consisting of 24 linkage groups was constructed based on 353,532 high-quality single nucleotide polymorphisms and 4,710 bin markers. A total of 28 QTLs corresponding to 11 genes, 26 QTLs corresponding to 11 genes, and 12 QTLs corresponding to 5 genes were identified for sex, intermuscular spine number and body weight traits, respectively. In this study, we assembled an accurate and nearly complete genome of C. alburnus by combining Illumina, PacBio, and high-throughput Chromosome conformation capture (Hi-C) technologies. In addition, we identified QTLs that explained variances in intermuscular spine number, body weight, and sex differences in C. alburnus. These genetic markers or candidate genes associated with growth traits provide a basis for marker-assisted selection in C. alburnus.