A manually annotated Actinidia chinensis var. chinensis (kiwifruit) genome highlights the challenges associated with draft genomes and gene prediction in plants

Unveiling kiwifruit TCP genes: evolution, functions, and expression insights

The TEOSINTE-BRANCHED1/CYCLOIDEA/PROLEFERATING-CELL-FACTORS (TCP) gene family is a plant-specific transcriptional factor family involved in leaf morphogenesis and senescence, lateral branching, hormone crosstalk, and stress responses. To date, a systematic study on the identification and characterization of the TCP gene family in kiwifruit has not been reported. Additionally, the function of kiwifruit TCPs in regulating kiwifruit responses to the ethylene treatment and bacterial canker disease pathogen (Pseudomonas syringae pv. actinidiae, Psa) has not been investigated. Here, we identified 40 and 26 TCP genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes, respectively. The synteny analysis of AcTCPs illustrated that whole-genome duplication accounted for the expansion of the TCP family in Ac. Phylogenetic, conserved domain, and selection pressure analysis indicated that TCP family genes in Ac and Ae had undergone different evolutionary patterns after whole-genome duplication (WGD) events, causing differences in TCP gene number and distribution. Our results also suggested that protein structure and cis-element architecture in promoter regions of TCP genes have driven the function divergence of duplicated gene pairs. Three and four AcTCP genes significantly affected kiwifruit responses to the ethylene treatment and Psa invasion, respectively. Our results provided insight into general characters, evolutionary patterns, and functional diversity of kiwifruit TCPs.

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.

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.

Identification of miRNA858 long-loop precursors in seed plants

MicroRNAs (miRNAs) are a class of nonprotein-coding short transcripts that provide a layer of post-transcriptional regulation essential to many plant biological processes. MiR858, which targets the transcripts of MYB transcription factors, can affect a range of secondary metabolic processes. Although miR858 and its 187-nt precursor have been well studied in Arabidopsis (Arabidopsis thaliana), a systematic investigation of miR858 precursors and their functions across plant species is lacking due to a problem in identifying the transcripts that generate this subclass. By re-evaluating the transcript of miR858 and relaxing the length cut-off for identifying hairpins, we found in kiwifruit (Actinidia chinensis) that miR858 has long-loop hairpins (1,100 to 2,100 nt), whose intervening sequences between miRNA generating complementary sites were longer than all previously reported miRNA hairpins. Importantly, these precursors of miR858 containing long-loop hairpins (termed MIR858L) are widespread in seed plants including Arabidopsis, varying between 350 and 5,500 nt. Moreover, we showed that MIR858L has a greater impact on proanthocyanidin and flavonol levels in both Arabidopsis and kiwifruit. We suggest that an active MIR858L-MYB regulatory module appeared in the transition of early land plants to large upright flowering plants, making a key contribution to plant secondary metabolism.

Chromosome-level genome of putative autohexaploid Actinidia deliciosa provides insights into polyploidisation and evolution

Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

The kiwifruit amyloplast proteome (kfALP): a resource to better understand the mechanisms underlying amyloplast biogenesis and differentiation

The biogenesis and differentiation (B&D) of amyloplasts contributes to fruit flavor and color. Here, remodeling of starch granules, thylakoids and plastoglobules was observed during development and ripening in two kiwifruit (Actinidia spp.) cultivars - yellow-fleshed 'Hort16A' and green-fleshed 'Hayward'. A protocol was developed to purify starch-containing plastids with a high degree of intactness, and amyloplast B&D was studied using label-free-based quantitative proteomic analyses in both cultivars. Over 3000 amyloplast-localized proteins were identified, of which >98% were quantified and defined as the kfALP (kiwifruit amyloplast proteome). The kfALP data were validated by Tandem-Mass-Tag (TMT) labeled proteomics in 'Hort16A'. Analysis of the proteomic data across development and ripening revealed: 1) a conserved increase in the abundance of proteins participating in starch synthesis/degradation during both amyloplast B&D; 2) up-regulation of proteins for chlorophyll degradation and of plastoglobule-localized proteins associated with chloroplast breakdown and plastoglobule formation during amyloplast differentiation; 3) constitutive expression of proteins involved in ATP supply and protein import during amyloplast B&D. Interestingly, two different pathways of amyloplast B&D were observed in the two cultivars. In 'Hayward', significant increases in abundance of photosynthetic- and tetrapyrrole metabolism-related proteins were observed, but the opposite trend was observed in 'Hort16A'. In conclusion, analysis of the kfALP provides new insights into the potential mechanisms underlying amyloplast B&D with relevance to key fruit quality traits in contrasting kiwifruit cultivars.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

Screening and identification of photoresponse factors in kiwifruit (Actinidia arguta) development

BACKGROUND

Light is essential for kiwifruit development, in which photoresponse factors contributes greatly to the quality formation. 'Light sensitive hypocotyls, also known as light-dependent short hypocotyls' (LSH) gene family can participate in fruit development as photoresponse factor. However, the key LSH gene that determine kiwifruit development remains unclear. This study aim to screen and identify the key gene AaLSH9 in A. arguta.

MATERIALS AND METHODS

Genome-wide identification of the LSH gene family was used to analyse LSH genes in kiwifruit. Homologous cloning was used to confirm the sequence of candidate LSH genes. qRT-PCR and cluster analysis of expression pattern were used to screen the key AaLSH9 gene. Subcellular localization of AaLSH9 in tobacco leaves and overexpression of AaLSH9 in Arabidopsis thaliana hy5 mutant plants were used to define the acting place in cell and identify molecular function, respectively.

RESULTS

We identified 15 LSH genes, which were divided into two sub-families namely A and B. Domain analysis of A and B showed that they contained different domain organizations, which possibly played key roles in the evolution process. Three LSH genes, AaLSH2, AaLSH9, and AaLSH11, were successfully isolated from Actinidia arguta. The expression pattern and cluster analysis of these three AaLSH genes suggested AaLSH9 might be a key photoresponse gene participating in fruit development in A. arguta. Subcellular localization showed AaLSH9 protein was located in the nucleus. The overexpression of AaLSH9 gene in Arabidopsis thaliana hy5 mutant plants partially complemented the long hypocotyls of hy5 mutant, implying AaLSH9 played a key role as photoresponse factor in cells. In addition, the seed coat color of A. thaliana over-expressing AaLSH9 became lighter than the wide type A.thaliana. Finally, AaCOP1 was confirmed as photoresponse factor to participate in developmental process by stable transgenic A. thaliana.

CONCLUSIONS

AaLSH9 can be involved in kiwifruit (A. arguta) development as key photoresponse factor. Our results not only identified the photoresponse factors AaLSH9 and AaCOP1 but also provided insights into their key role in fruit quality improvement in the process of light response.

Agrobacterium rhizogenes-mediated marker-free transformation and gene editing system revealed that AeCBL3 mediates the formation of calcium oxalate crystal in kiwifruit

The transformation and gene editing of the woody species kiwifruit are difficult and time-consuming. The fast and marker-free genetic modification system for kiwifruit has not been developed yet. Here, we establish a rapid and efficient marker-free transformation and gene editing system mediated by Agrobacterium rhizogenes for kiwifruit. Moreover, a removing-root-tip method was developed to significantly increase the regeneration efficiency of transgenic hairy roots. Through A. rhizogenes-mediated CRISPR/Cas9 gene editing, the editing efficiencies of CEN4 and AeCBL3 achieved 55 and 50%, respectively. And several homozygous knockout lines for both genes were obtained. Our method has been successfully applied in the transformation of two different species of kiwifruit (Actinidia chinensis 'Hongyang' and A.eriantha 'White'). Next, we used the method to study the formation of calcium oxalate (CaOx) crystals in kiwifruit. To date, little is known about how CaOx crystal is formed in plants. Our results indicated that AeCBL3 overexpression enhanced CaOx crystal formation, but its knockout via CRISPR/Cas9 significantly impaired crystal formation in kiwifruit. Together, we developed a fast maker-free transformation and highly efficient CRISPR-Cas9 gene editing system for kiwifruit. Moreover, our work revealed a novel gene mediating CaOx crystal formation and provided a clue to elaborate the underlying mechanisms.

Genome-wide identification and expression pattern analysis of the kiwifruit GRAS transcription factor family in response to salt stress

BACKGROUND

GRAS is a family of plant-specific transcription factors (TFs) that play a vital role in plant growth and development and response to adversity stress. However, systematic studies of the GRAS TF family in kiwifruit have not been reported.

RESULTS

In this study, we used a bioinformatics approach to identify eighty-six AcGRAS TFs located on twenty-six chromosomes and phylogenetic analysis classified them into ten subfamilies. It was found that the gene structure is relatively conserved for these genes and that fragmental duplication is the prime force for the evolution of AcGRAS genes. However, the promoter region of the AcGRAS genes mainly contains cis-acting elements related to hormones and environmental stresses, similar to the results of GO and KEGG enrichment analysis, suggesting that hormone signaling pathways of the AcGRAS family play a vital role in regulating plant growth and development and adversity stress. Protein interaction network analysis showed that the AcGRAS51 protein is a relational protein linking DELLA, SCR, and SHR subfamily proteins. The results demonstrated that 81 genes were expressed in kiwifruit AcGRAS under salt stress, including 17 differentially expressed genes, 13 upregulated, and four downregulated. This indicates that the upregulated AcGRAS55, AcGRAS69, AcGRAS86 and other GRAS genes can reduce the salt damage caused by kiwifruit plants by positively regulating salt stress, thus improving the salt tolerance of the plants.

CONCLUSIONS

These results provide a theoretical basis for future exploration of the characteristics and functions of more AcGRAS genes. This study provides a basis for further research on kiwifruit breeding for resistance to salt stress. RT-qPCR analysis showed that the expression of 3 AcGRAS genes was elevated under salt stress, indicating that AcGRAS exhibited a specific expression pattern under salt stress conditions.

AcHZP45 is a repressor of chlorophyll biosynthesis and activator of chlorophyll degradation in kiwifruit

The degradation of chlorophyll during fruit development is essential to reveal a more 'ripe' color that signals readiness to wild dispersers of seeds and the human consumer. Here, comparative biochemical analysis of developing fruit of Actinidia deliciosa cv. Xuxiang ('XX', green-fleshed) and Actinidia chinensis cv. Jinshi No.1 ('JS', yellow-fleshed) indicated that variation in chlorophyll content is the major contributor to differences in flesh color. Four differentially expressed candidate genes were identified: the down-regulated genes AcCRD1 and AcPOR1 involved in chlorophyll biosynthesis, and the up-regulated genes AcSGR1 and AcSGR2 driving chlorophyll degradation. Prochlorophyllide and chlorophyllide, the metabolites produced by AcCRD1 and AcPOR1, progressively reduced in 'JS', but not in 'XX', indicating that chlorophyll biosynthesis was less active in yellow-fleshed fruit. AcSGR1 and AcSGR2 were verified to be involved in chlorophyll degradation, using both transient expression in tobacco and stable overexpression in kiwifruit. Furthermore, a homeobox-leucine zipper (HD-Zip II), AcHZP45, showed significantly increased expression during 'JS' fruit ripening, which led to both repressed expression of AcCRD1 and AcPOR1 and activated expression of AcSGR1 and AcSGR2. Collectively, the present study indicated that different dynamics of chlorophyll biosynthesis and degradation coordinate the changes in chlorophyll content in kiwifruit flesh, which are orchestrated by the key transcription factor AcHZP45.

Genomic architecture of resistance to latania scale (H. lataniae) in kiwifruit (A. chinensis var. chinensis)

BACKGROUND

Latania scale (Hemiberlesia lataniae Signoret) is an armoured scale insect known to cause damage to kiwifruit plants and fruit, which ultimately reduces crop values and creates post-harvest export and quarantine issues. Resistance to H. lataniae does exist in some commercial cultivars of kiwifruit. However, some of the commercial cultivars bred in New Zealand have not inherited alleles for resistance to H. lataniae carried by their parents. To elucidate the architecture of resistance in the parents and develop molecular markers to assist breeding, these experiments analysed the inheritance of resistance to H. lataniae from families related to commercial cultivars.

RESULTS

The first experiment identified a 15.97 Mb genomic region of interest for resistance to H. lataniae in rtGBS data of 3.23 to 19.20 Mb on chromosome 10. A larger population was then QTL mapped, which confirmed the region of interest as the sole locus contributing to H. lataniae resistance. inDel markers mapping the region of low recombination under the QTL peak further narrowed the region associated with H. lataniae resistance to a 5.73 Mb region.

CONCLUSIONS

The kiwifruit populations and genomic methods used in this study identify the same non-recombinant region of chromosome 10 which confers resistance of A. chinensis var. chinensis to H. lataniae. The markers developed to target the H. lataniae resistance loci will reduce the amount of costly and time-consuming phenotyping required for breeding H. lataniae scale resistance into new kiwifruit cultivars.

Genome-wide identification of the heat shock transcription factor gene family in two kiwifruit species

High temperatures have a significant impact on plant growth and metabolism. In recent years, the fruit industry has faced a serious threat due to high-temperature stress on fruit plants caused by global warming. In the present study, we explored the molecular regulatory mechanisms that contribute to high-temperature tolerance in kiwifruit. A total of 36 Hsf genes were identified in the A. chinensis (Ac) genome, while 41 Hsf genes were found in the A. eriantha (Ae) genome. Phylogenetic analysis revealed the clustering of kiwifruit Hsfs into three distinct groups (groups A, B, and C). Synteny analysis indicated that the expansion of the Hsf gene family in the Ac and Ae genomes was primarily driven by whole genome duplication (WGD). Analysis of the gene expression profiles revealed a close relationship between the expression levels of Hsf genes and various plant tissues and stress treatments throughout fruit ripening. Subcellular localization analysis demonstrated that GFP-AcHsfA2a/AcHsfA7b and AcHsfA2a/AcHsfA7b -GFP were localized in the nucleus, while GFP-AcHsfA2a was also observed in the cytoplasm of Arabidopsis protoplasts. The results of real-time quantitative polymerase chain reaction (RT-qPCR) and dual-luciferase reporter assay revealed that the majority of Hsf genes, especially AcHsfA2a, were expressed under high-temperature conditions. In conclusion, our findings establish a theoretical foundation for analyzing the potential role of Hsfs in high-temperature stress tolerance in kiwifruit. This study also offers valuable information to aid plant breeders in the development of heat-stress-resistant plant materials.

Manual correction of genome annotation improved alternative splicing identification of Artemisia annua

Gene annotation is essential for genome-based studies. However, algorithm-based genome annotation is difficult to fully and correctly reveal genomic information, especially for species with complex genomes. Artemisia annua L. is the only commercial resource of artemisinin production though the content of artemisinin is still to be improved. Genome-based genetic modification and breeding are useful strategies to boost artemisinin content and therefore, ensure the supply of artemisinin and reduce costs, but better gene annotation is urgently needed. In this study, we manually corrected the newly released genome annotation of A. annua using second- and third-generation transcriptome data. We found that incorrect gene information may lead to differences in structural, functional, and expression levels compared to the original expectations. We also identified alternative splicing events and found that genome annotation information impacted identifying alternative splicing genes. We further demonstrated that genome annotation information and alternative splicing could affect gene expression estimation and gene function prediction. Finally, we provided a valuable version of A. annua genome annotation and demonstrated the importance of gene annotation in future research.

Comparative transcriptomic and plastid development analysis sheds light on the differential carotenoid accumulation in kiwifruit flesh

Carotenoids are colorful lipophilic isoprenoids synthesized in all photosynthetic organisms which play roles in plant growth and development and provide numerous health benefits in the human diet (precursor of Vitamin A). The commercially popular kiwifruits are golden yellow-fleshed (Actinidia chinensis) and green fleshed (A. deliciosa) cultivars which have a high carotenoid concentration. Understanding the molecular mechanisms controlling the synthesis and sequestration of carotenoids in Actinidia species is key to increasing nutritional value of this crop via breeding. In this study we analyzed fruit with varying flesh color from three Actinidia species; orange-fleshed A. valvata (OF), yellow-fleshed A. polygama (YF) and green-fleshed A. arguta (GF). Microscopic analysis revealed that carotenoids accumulated in a crystalline form in YF and OF chromoplasts, with the size of crystals being bigger in OF compared to YF, which also contained globular substructures in the chromoplast. Metabolic profiles were investigated using ultra-performance liquid chromatography (UPLC), which showed that β-carotene was the predominant carotenoid in the OF and YF species, while lutein was the dominant carotenoid in the GF species. Global changes in gene expression were studied between OF and GF (both tetraploid) species using RNA-sequencing which showed higher expression levels of upstream carotenoid biosynthesis-related genes such as DXS, PSY, GGPPS, PDS, ZISO, and ZDS in OF species compared to GF. However, low expression of downstream pathway genes was observed in both species. Pathway regulatory genes (OR and OR-L), plastid morphology related genes (FIBRILLIN), chlorophyll degradation genes (SGR, SGR-L, RCCR, and NYC1) were upregulated in OF species compared to GF. This suggests chlorophyll degradation (primarily in the initial ripening stages) is accompanied by increased carotenoid production and localization in orange flesh tissue, a contrast from green flesh tissue. These results suggest a coordinated change in the carotenoid pathway, as well as changes in plastid type, are responsible for an orange phenotype in certain kiwifruit species.

Waterlogging Tolerance of Actinidia valvata Dunn Is Associated with High Activities of Pyruvate Decarboxylase, Alcohol Dehydrogenase and Antioxidant Enzymes

Kiwifruit (Actinidia spp.) is susceptible to waterlogging stress. Although abundant wild germplasm resources exist among Actinidia plants for improving the waterlogging tolerance of kiwifruit cultivars, the underlying mechanisms remain largely unknown. Here, a comparative study was undertaken using one wild germplasm, Maorenshen (A. valvata Dunn, MRS), and one cultivar, Miliang-1 (A. chinensis var. deliciosa (A.Chev.) A.Chev. cv. Miliang-1, ML). Under stress, the ML plantlets were seriously damaged with wilted chlorotic leaves and blackened rotten roots, whereas the symptoms of injury in the MRS plantlets were much fewer, along with higher photosynthetic rates, chlorophyll fluorescence characteristics and root activity under stress conditions. However, neither aerenchyma in the root nor adventitious roots appeared in both germplasms upon stress exposure. The activities of pyruvate decarboxylase (PDC) and alcohol dehydrogenase (ADH), as well as their transcript levels, were constitutively higher in MRS than those in ML under both normal and stress conditions. Waterlogging stress significantly enhanced the PDC and ADH enzyme activities in both germplasms, which were 60.8% and 22.4% higher in the MRS roots than those in the ML roots under waterlogging stress, respectively. Moreover, MRS displayed higher activities of antioxidant enzymes, including SOD, CAT, and APX, as well as DPPH-radical scavenging ability, and decreased H2O2 and MDA accumulation under both normal and stress conditions. Our findings suggest that the waterlogging tolerance of the wild A. valvata germplasm was associated with high PDC and ADH, as well as antioxidant ability.

The hops (Humulus lupulus) genome contains a mid-sized terpene synthase family that shows wide functional and allelic diversity

BACKGROUND

Hops (Humulus lupulus L.) are a dioecious climbing perennial, with the dried mature "cones" (strobili) of the pistillate/female inflorescences being widely used as both a bittering agent and to enhance the flavour of beer. The glandular trichomes of the bract and bracteole flowering structures of the cones produce an abundance of secondary metabolites, such as terpenoids, bitter acids and prenylated phenolics depending on plant genetics, developmental stage and environment. More knowledge is required on the functional and allelic diversity of terpene synthase (TPS) genes responsible for the biosynthesis of volatile terpenes to assist in flavour-directed hop breeding.

RESULTS

Major volatile terpene compounds were identified using gas chromatography-mass spectrometry (GC-MS) in the ripe cones of twenty-one hop cultivars grown in New Zealand. All cultivars produced the monoterpene β-myrcene and the sesquiterpenes α-humulene and β-caryophyllene, but the quantities varied broadly. Other terpenes were found in large quantities in only a smaller subset of cultivars, e.g. β-farnesene (in seven cultivars) and α-pinene (in four). In four contrasting cultivars (Wakatu™, Wai-iti™, Nelson Sauvin™, and 'Nugget'), terpene production during cone development was investigated in detail, with concentrations of some of the major terpenes increasing up to 1000-fold during development and reaching maximal levels from 50-60 days after flowering. Utilising the published H. lupulus genome, 87 putative full-length and partial terpene synthase genes were identified. Alleles corresponding to seven TPS genes were amplified from ripe cone cDNA from multiple cultivars and subsequently functionally characterised by transient expression in planta. Alleles of the previously characterised HlSTS1 produced humulene/caryophyllene as the major terpenes. HlRLS alleles produced (R)-(-)-linalool, whilst alleles of two sesquiterpene synthase genes, HlAFS1 and HlAFS2 produced α-farnesene. Alleles of HlMTS1, HlMTS2 and HlTPS1 were inactive in all the hop cultivars studied.

CONCLUSIONS

Alleles of four TPS genes were identified and shown to produce key aroma volatiles in ripe hop cones. Multiple expressed but inactive TPS alleles were also identified, suggesting that extensive loss-of-function has occurred during domestication and breeding of hops. Our results can be used to develop hop cultivars with novel/improved terpene profiles using marker-assisted breeding strategies to select for, or against, specific TPS alleles.

The terpene synthase (TPS) gene family in kiwifruit shows high functional redundancy and a subset of TPS likely fulfil overlapping functions in fruit flavour, floral bouquet and defence

Volatile terpenes are important compounds that influence fruit flavour and aroma of kiwifruit. Terpenes in plants also impact on the floral bouquet and defence against pests and pathogens in leaves and fruit. To better understand the overlapping roles that terpenes may fulfil in plants, a systematic gene, chemical and biochemical analysis of terpenes and terpene synthases (TPS) was undertaken in Red5 kiwifruit (Actinidia spp.). Analysis of the Red5 genome shows it contains only 22 TPS gene models, of which fifteen encode full-length TPS. Thirteen TPS can account for the major terpene volatiles produced in different tissues of Red5 kiwifruit and in response to different stimuli. The small Red5 TPS family displays surprisingly high functional redundancy with five TPS producing linalool/nerolidol. Treatment of leaves and fruit with methyl jasmonate enhanced expression of a subset of defence-related TPS genes and stimulated the release of terpenes. Six TPS genes were induced upon herbivory of leaves by the economically important insect pest Ctenopseustis obliquana (brown-headed leaf roller) and emission, but not accumulation, of (E)- and (Z)-nerolidol was strongly linked to herbivory. Our results provide a framework to understand the overlapping biological and ecological roles of terpenes in Actinidia and other horticultural crops.

Identification and Characterization of Long Non-Coding RNAs: Implicating Insights into Their Regulatory Role in Kiwifruit Ripening and Softening during Low-Temperature Storage

Long non-coding RNAs (lncRNAs) are crucial players regulating many biological processes in plants. However, limited knowledge is available regarding their roles in kiwifruit ripening and softening. In this study, using lncRNA-seq technology, 591 differentially expressed (DE) lncRNAs (DELs) and 3107 DE genes (DEGs) were identified from kiwifruit stored at 4 °C for 1, 2, and 3 weeks in comparison with non-treated control fruits. Of note, 645 DEGs were predicted to be targets of DELs (DEGTLs), including some DE protein-coding genes (such as β-amylase and pectinesterase). DEGTL-based GO enrichment analysis revealed that these genes were significantly enriched in cell wall modification and pectinesterase activity in 1 W vs. CK and 3 W vs. CK, which might be closely related to the fruit softening during low-temperature storage. Moreover, KEGG enrichment analysis revealed that DEGTLs were significantly associated with starch and sucrose metabolism. Our study revealed that lncRNAs play critical regulatory roles in kiwifruit ripening and softening under low-temperature storage, mainly by mediating the expression of starch and sucrose metabolism and cell wall modification related genes.

Telomere-to-telomere and haplotype-resolved genome of the kiwifruit Actinidia eriantha

Actinidia eriantha is a characteristic fruit tree featuring with great potential for its abundant vitamin C and strong disease resistance. It has been used in a wide range of breeding programs and functional genomics studies. Previously published genome assemblies of A. eriantha are quite fragmented and not highly contiguous. Using multiple sequencing strategies, we get the haplotype-resolved and gap-free genomes of an elite breeding line "Midao 31" (MD), termed MDHAPA and MDHAPB. The new assemblies anchored to 29 pseudochromosome pairs with a length of 619.3 Mb and 611.7 Mb, as well as resolved 27 and 28 gap-close chromosomes in a telomere-to-telomere (T2T) manner. Based on the haplotype-resolved genome, we found that most alleles experienced purifying selection and coordinately expressed. Owing to the high continuity of assemblies, we defined the centromeric regions of A. eriantha, and identified the major repeating monomer, which is designated as Ae-CEN153. This resource lays a solid foundation for further functional genomics study and horticultural traits improvement in kiwifruit.

Chromosome-scale genome assembly of a natural diploid kiwifruit (Actinidia chinensis var. deliciosa)

The most commercialized kiwifruit, Actinidia chinensis var. deliciosa (Acd), is an allohexaploid (2n = 6x = 174), making high-quality assemblage genome challenging. We previously discovered a rare naturally occurring diploid Acd plant. Here, chromosome-level de novo genome assembly for this diploid Acd was reported, reaching approximately 621.98 Mb in length with contig and scaffold N50 values of 10.08 and 21.09 Mb, respectively, 99.66% of the bases anchored to 29 pseudochromosomes, and 38,990 protein-coding genes and 42.29% repetitive elements annotated. The divergence time of A. chinensis cv. 'Red5' and 'Hongyang' (11.1-27.7 mya) was more recent compared with the divergence time of them and Acd (19.9-41.2 mya), with the divergence time of A. eriantha cv. 'White' being the earliest (22.9-45.7 mya) among that of the four Actinidia species. The 4DTv distance distribution highlighted three recent whole-genome duplication events in Acd. This is the first high-quality diploid Acd genome, which lays an important foundation for not only kiwifruit functional genomics studies but also further elucidating genome evolution of allohexaploid Acd.

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

Kiwifruit is a recently domesticated horticultural fruit crop with substantial economic and nutritional value, especially because of the high content of vitamin C in its fruit. In this study, we de novo assembled two telomere-to-telomere kiwifruit genomes from Actinidia chinensis var. 'Donghong' (DH) and Actinidia latifolia 'Kuoye' (KY), with total lengths of 608 327 852 and 640 561 626 bp for 29 chromosomes, respectively. With a burst of structural variants involving inversion, translocations, and duplications within 8.39 million years, the metabolite content of DH and KY exhibited differences in saccharides, lignans, and vitamins. A regulatory ERF098 transcription factor family has expanded in KY and Actinidia eriantha, both of which have ultra-high vitamin C content. With each assembly phased into two complete haplotypes, we identified allelic variations between two sets of haplotypes, leading to protein sequence variations in 26 494 and 27 773 gene loci and allele-specific expression of 4687 and 12 238 homozygous gene pairs. Synchronized metabolome and transcriptome changes during DH fruit development revealed the same dynamic patterns in expression levels and metabolite contents; free fatty acids and flavonols accumulated in the early stages, but sugar substances and amino acids accumulated in the late stages. The AcSWEET9b gene that exhibits allelic dominance was further identified to positively correlate with high sucrose content in fruit. Compared with wild varieties and other Actinidia species, AcSWEET9b promoters were selected in red-flesh kiwifruits that have increased fruit sucrose content, providing a possible explanation on why red-flesh kiwifruits are sweeter. Collectively, these two gap-free kiwifruit genomes provide a valuable genetic resource for investigating domestication mechanisms and genome-based breeding of kiwifruit.

Two genes, ANS and UFGT2, from Vaccinium spp. are key steps for modulating anthocyanin production

Anthocyanins are a major group of red to blue spectrum plant pigments with many consumer health benefits. Anthocyanins are derived from the flavonoid pathway and diversified by glycosylation and methylation, involving the concerted action of specific enzymes. Blueberry and bilberry (Vaccinium spp.) are regarded as 'superfruits' owing to their high content of flavonoids, especially anthocyanins. While ripening-related anthocyanin production in bilberry (V. myrtillus) and blueberry (V. corymbosum) is regulated by the transcriptional activator MYBA1, the role of specific structural genes in determining the concentration and composition of anthocyanins has not been functionally elucidated. We isolated three candidate genes, CHALCONE SYNTHASE (VmCHS1), ANTHOCYANIDIN SYNTHASE (VmANS) and UDP-GLUCOSE : FLAVONOID-3-O-GLYCOSYLTRANSFERASE (VcUFGT2), from Vaccinium, which were predominantly expressed in pigmented fruit skin tissue and showed high homology between bilberry and blueberry. Agrobacterium-mediated transient expression of Nicotiana benthamiana showed that overexpression of VcMYBA1 in combination with VmANS significantly increased anthocyanin concentration (3-fold). Overexpression of VmCHS1 showed no effect above that induced by VcMYBA1, while VcUFGT2 modulated anthocyanin composition to produce delphinidin-3-galactosylrhamnoside, not naturally produced in tobacco. In strawberry (Fragaria × ananassa), combined transient overexpression of VcUFGT2 with a FLAVONOID 3´,5´-HYDROXYLASE from kiwifruit (Actinidia melanandra) modulated the anthocyanin profile to include galactosides and arabinosides of delphinidin and cyanidin, major anthocyanins in blueberry and bilberry. These findings provide insight into the role of the final steps of biosynthesis in modulating anthocyanin production in Vaccinium and may contribute to the targeted breeding of new cultivars with improved nutritional properties.

Interaction of AcMADS68 with transcription factors regulates anthocyanin biosynthesis in red-fleshed kiwifruit

In red-fleshed kiwifruit, anthocyanin pigmentation is a crucial commercial trait. The MYB-bHLH-WD40 (MBW) complex and other transcription factors regulate its accumulation. Herein, a new SEP gene, AcMADS68, was identified as a regulatory candidate for anthocyanin biosynthesis in the kiwifruit by transcriptome data and bioinformatic analyses. AcMADS68 alone could not induce the accumulation of anthocyanin both in Actinidia arguta fruit and tobacco leaves. However, in combination with AcMYBF110, AcMYB123, and AcbHLH1, AcMADS68 co-overexpression increased anthocyanin biosynthesis, whereas its silencing reduced anthocyanin accumulation. The results of the dual-luciferase reporter, firefly luciferase complementation, yeast two-hybrid and co-immunoprecipitation assays showed that AcMADS68 could interact with both AcMYBF110 and AcMYB123 but not with AcbHLH1, thereby co-regulating anthocyanin biosynthesis by promoting the activation of the target genes, including AcANS, AcF3GT1, and AcGST1. Moreover, AcMADS68 also could activate the promoter of AcbHLH1 surported by dual-luciferase reporter and yeast one-hybrid assays, thereby further amplifying the regulation signals from the MBW complex, thus resulting in enhanced anthocyanin accumulation in the kiwifruit. These findings may facilitate better elucidation of various regulatory mechanisms underlying anthocyanin accumulation and contribute to the quality enhancement of red-fleshed kiwifruit.

Development of a 135K SNP genotyping array for Actinidia arguta and its applications for genetic mapping and QTL analysis in kiwifruit

Kiwifruit (Actinidia spp) is a woody, perennial and deciduous vine. In this genus, there are multiple ploidy levels but the main cultivated cultivars are polyploid. Despite the availability of many genomic resources in kiwifruit, SNP genotyping is still a challenge given these different levels of polyploidy. Recent advances in SNP array technologies have offered a high-throughput genotyping platform for genome-wide DNA polymorphisms. In this study, we developed a high-density SNP genotyping array to facilitate genetic studies and breeding applications in kiwifruit. SNP discovery was performed by genome-wide DNA sequencing of 40 kiwifruit genotypes. The identified SNPs were stringently filtered for sequence quality, predicted conversion performance and distribution over the available Actinidia chinensis genome. A total of 134 729 unique SNPs were put on the array. The array was evaluated by genotyping 400 kiwifruit individuals. We performed a multidimensional scaling analysis to assess the diversity of kiwifruit germplasm, showing that the array was effective to distinguish kiwifruit accessions. Using a tetraploid F1 population, we constructed an integrated linkage map covering 3060.9 cM across 29 linkage groups and performed QTL analysis for the sex locus that has been identified on Linkage Group 3 (LG3) in Actinidia arguta. Finally, our dataset presented evidence of tetrasomic inheritance with partial preferential pairing in A. arguta. In conclusion, we developed and evaluated a 135K SNP genotyping array for kiwifruit. It has the advantage of a comprehensive design that can be an effective tool in genetic studies and breeding applications in this high-value crop.

Telomere-to-telomere and gap-free reference genome assembly of the kiwifruit Actinidia chinensis

Kiwifruit is an economically and nutritionally important fruit crop with extremely high contents of vitamin C. However, the previously released versions of kiwifruit genomes all have a mass of unanchored or missing regions. Here, we report a highly continuous and completely gap-free reference genome of Actinidia chinensis cv. 'Hongyang', named Hongyang v4.0, which is the first to achieve two de novo haploid-resolved haplotypes, HY4P and HY4A. HY4P and HY4A have a total length of 606.1 and 599.6 Mb, respectively, with almost the entire telomeres and centromeres assembled in each haplotype. In comparison with Hongyang v3.0, the integrity and contiguity of Hongyang v4.0 is markedly improved by filling all unclosed gaps and correcting some misoriented regions, resulting in ~38.6-39.5 Mb extra sequences, which might affect 4263 and 4244 protein-coding genes in HY4P and HY4A, respectively. Furthermore, our gap-free genome assembly provides the first clue for inspecting the structure and function of centromeres. Globally, centromeric regions are characterized by higher-order repeats that mainly consist of a 153-bp conserved centromere-specific monomer (Ach-CEN153) with different copy numbers among chromosomes. Functional enrichment analysis of the genes located within centromeric regions demonstrates that chromosome centromeres may not only play physical roles for linking a pair of sister chromatids, but also have genetic features for participation in the regulation of cell division. The availability of the telomere-to-telomere and gap-free Hongyang v4.0 reference genome lays a solid foundation not only for illustrating genome structure and functional genomics studies but also for facilitating kiwifruit breeding and improvement.

Genome-Wide Identification of Kiwifruit SGR Family Members and Functional Characterization of SGR2 Protein for Chlorophyll Degradation

The STAY-GREEN (SGR) proteins play an important role in chlorophyll (Chl) degradation and are closely related to plant photosynthesis. However, the availability of inadequate studies on SGR motivated us to conduct a comprehensive study on the identification and functional dissection of SGR superfamily members in kiwifruit. Here, we identified five SGR genes for each of the kiwifruit species [Actinidia chinensis (Ac) and Actinidia eriantha (Ae)]. The phylogenetic analysis showed that the kiwifruit SGR superfamily members were divided into two subfamilies the SGR subfamily and the SGRL subfamily. The results of transcriptome data and RT-qPCR showed that the expression of the kiwifruit SGRs was closely related to light and plant developmental stages (regulated by plant growth regulators), which were further supported by the presence of light and the plant hormone-responsive cis-regulatory element in the promoter region. The subcellular localization analysis of the AcSGR2 protein confirmed its localization in the chloroplast. The Fv/Fm, SPAD value, and Chl contents were decreased in overexpressed AcSGR2, but varied in different cultivars of A. chinensis. The sequence analysis showed significant differences within AcSGR2 proteins. Our findings provide valuable insights into the characteristics and evolutionary patterns of SGR genes in kiwifruit, and shall assist kiwifruit breeders to enhance cultivar development.

Genome-wide analysis of MADS-box gene family in kiwifruit (Actinidia chinensis var. chinensis) and their potential role in floral sex differentiation

Kiwifruit (Actinidia chinensis Planch.) is a functionally dioecious plant, which displays diverse morphology in male and female flowers. MADS-box is an ancient and huge gene family that plays a key role in plant floral organ differentiation. In this study, we have identified 89 MADS-box genes from A. chinensis Red 5 genome. These genes are distributed on 26 chromosomes and are classified into type I (21 genes) and type II (68 genes). Overall, type II AcMADS-box genes have more complex structures than type I with more exons, protein domains, and motifs, indicating that type II genes may have more diverse functions. Gene duplication analysis showed that most collinearity occurred in type II AcMADS-box genes, which was consistent with a large number of type II genes. Analysis of cis-acting elements in promoters showed that AcMADS-box genes are mainly associated with light and phytohormone responsiveness. The expression profile of AcMADS-box genes in different tissues showed that most genes were highly expressed in flowers. Further, the qRT-PCR analysis of the floral organ ABCDE model-related genes in male and female flowers revealed that AcMADS4, AcMADS56, and AcMADS70 were significantly expressed in female flowers. It indicated that those genes may play an important role in the sex differentiation of kiwifruit. This work provided a comprehensive analysis of the AcMADS-box genes and may help facilitate our understanding of the sex differentiation regulatory mechanism in kiwifruit.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

CRISPR-Cas9-mediated mutagenesis of kiwifruit BFT genes results in an evergrowing but not early flowering phenotype

Phosphatidylethanolamine-binding protein (PEBP) genes regulate flowering and architecture in many plant species. Here, we study kiwifruit (Actinidia chinensis, Ac) PEBP genes with homology to BROTHER OF FT AND TFL1 (BFT). CRISPR-Cas9 was used to target AcBFT genes in wild-type and fast-flowering kiwifruit backgrounds. The editing construct was designed to preferentially target AcBFT2, whose expression is elevated in dormant buds. Acbft lines displayed an evergrowing phenotype and increased branching, while control plants established winter dormancy. The evergrowing phenotype, encompassing delayed budset and advanced budbreak after defoliation, was identified in multiple independent lines with edits in both alleles of AcBFT2. RNA-seq analyses conducted using buds from gene-edited and control lines indicated that Acbft evergrowing plants had a transcriptome similar to that of actively growing wild-type plants, rather than dormant controls. Mutations in both alleles of AcBFT2 did not promote flowering in wild-type or affect flowering time, morphology and fertility in fast-flowering transgenic kiwifruit. In summary, editing of AcBFT2 has the potential to reduce plant dormancy with no adverse effect on flowering, giving rise to cultivars better suited for a changing climate.

A Comparative Genomics Approach for Analysis of Complete Mitogenomes of Five Actinidiaceae Plants

Actinidiaceae, an economically important plant family, includes the Actinidia, Clematoclethra and Saurauia genus. Kiwifruit, with remarkably high vitamin C content, is an endemic species widely distributed in China with high economic value. Although many Actinidiaceae chloroplast genomes have been reported, few complete mitogenomes of Actinidiaceae have been studied. Here, complete circular mitogenomes of the four kiwifruit species and Saurauia tristyla were assembled. Codon usage, sequence repeats, RNA editing, gene transfers, selective pressure, and phylogenetic relationships in the four kiwifruit species and S. tristyla were comparatively analyzed. This research will contribute to the study of phylogenetic relationships within Actiniaceae and molecular barcoding in kiwifruit.

Alcohol acyl transferase genes at a high-flavor intensity locus contribute to ester biosynthesis in kiwifruit

Volatile esters are key compounds contributing to flavor intensity in commonly consumed fruits including apple (Malus domestica), strawberry (Fragaria spp.), and banana (Musa sapientum). In kiwifruit (Actinidia spp.), ethyl butanoate and other esters have been proposed to contribute fruity, sweet notes to commercial cultivars. Here, we investigated the genetic basis for ester production in Actinidia in an A. chinensis mapping population (AcMPO). A major quantitative trait loci for the production of multiple esters was identified at the high-flavor intensity (HiFI) locus on chromosome 20. This locus co-located with eight tandemly arrayed alcohol acyl transferase genes in the Red5 genome that were expressed in a ripening-specific fashion that corresponded with ester production. Biochemical characterization suggested two genes at the HiFI locus, alcohol acyl transferase 16-b/c (AT16-MPb/c), probably contributed most to the production of ethyl butanoate. A third gene, AT16-MPa, probably contributed more to hexyl butanoate and butyl hexanoate production, two esters that segregated in AcMPO. Sensory analysis of AcMPO indicated that fruit from segregating lines with high ester concentrations were more commonly described as being "fruity" as opposed to "beany". The downregulation of AT16-MPa-c by RNAi reduced ester production in ripe "Hort16A" fruit by >90%. Gas chromatography-olfactometry indicated the loss of the major "fruity" notes contributed by ethyl butanoate. A comparison of unimproved Actinidia germplasm with those of commercial cultivars indicated that the selection of fruit with high concentrations of alkyl esters (but not green note aldehydes) was probably an important selection trait in kiwifruit cultivation. Understanding ester production at the HiFI locus is a critical step toward maintaining and improving flavor intensity in kiwifruit.

Genomic and Bioinformatic Resources for Perennial Fruit Species

In the post-genomic era, data management and development of bioinformatic tools are critical for the adequate exploitation of genomics data. In this review, we address the actual situation for the subset of crops represented by the perennial fruit species. The agronomical singularity of these species compared to plant and crop model species provides significant challenges on the implementation of good practices generally not addressed in other species. Studies are usually performed over several years in non-controlled environments, usage of rootstock is common, and breeders heavily rely on vegetative propagation. A reference genome is now available for all the major species as well as many members of the economically important genera for breeding purposes. Development of pangenome for these species is beginning to gain momentum which will require a substantial effort in term of bioinformatic tool development. The available tools for genome annotation and functional analysis will also be presented.

Genome-Wide Expression Profiling Analysis of Kiwifruit GolS and RFS Genes and Identification of AcRFS4 Function in Raffinose Accumulation

The raffinose synthetase (RFS) and galactinol synthase (GolS) are two critical enzymes for raffinose biosynthesis, which play an important role in modulating plant growth and in response to a variety of biotic or abiotic stresses. Here, we comprehensively analyzed the RFS and GolS gene families and their involvement in abiotic and biotic stresses responses at the genome-wide scale in kiwifruit. A total of 22 GolS and 24 RFS genes were identified in Actinidia chinensis and Actinidia eriantha genomes. Phylogenetic analysis showed that the GolS and RFS genes were clustered into four and six groups, respectively. Transcriptomic analysis revealed that abiotic stresses strongly induced some crucial genes members including AcGolS1/2/4/8 and AcRFS2/4/8/11 and their expression levels were further confirmed by qRT-PCR. The GUS staining of AcRFS4Pro::GUS transgenic plants revealed that the transcriptionlevel of AcRFS4 was significantly increased by salt stress. Overexpression of AcRFS4 in Arabidopsis demonstrated that this gene enhanced the raffinose accumulation and the tolerance to salt stress. The co-expression networks analysis of hub transcription factors targeting key AcRFS4 genes indicated that there was a strong correlation between AcNAC30 and AcRFS4 expression under salt stress. Furthermore, the yeast one-hybrid assays showed that AcNAC30 could bind the AcRFS4 promoter directly. These results may provide insights into the evolutionary and functional mechanisms of GolS and RFS genes in kiwifruit.

The red flesh of kiwifruit is differentially controlled by specific activation-repression systems

Anthocyanins are visual cues for pollination and seed dispersal. Fruit containing anthocyanins also appeals to consumers due to its appearance and health benefits. In kiwifruit (Actinidia spp.) studies have identified at least two MYB activators of anthocyanin, but their functions in fruit and the mechanisms by which they act are not fully understood. Here, transcriptome and small RNA high-throughput sequencing were used to comprehensively identify contributors to anthocyanin accumulation in kiwifruit. Stable overexpression in vines showed that both 35S::MYB10 and MYB110 can upregulate anthocyanin biosynthesis in Actinidia chinensis fruit, and that MYB10 overexpression resulted in anthocyanin accumulation which was limited to the inner pericarp, suggesting that repressive mechanisms underlie anthocyanin biosynthesis in this species. Furthermore, motifs in the C-terminal region of MYB10/110 were shown to be responsible for the strength of activation of the anthocyanic response. Transient assays showed that both MYB10 and MYB110 were not directly cleaved by miRNAs, but that miR828 and its phased small RNA AcTAS4-D4(-) efficiently targeted MYB110. Other miRNAs were identified, which were differentially expressed between the inner and outer pericarp, and cleavage of SPL13, ARF16, SCL6 and F-box1, all of which are repressors of MYB10, was observed. We conclude that it is the differential expression and subsequent repression of MYB activators that is responsible for variation in anthocyanin accumulation in kiwifruit species.

Two plant NLR proteins confer strain-specific resistance conditioned by an effector from Pseudomonas syringae pv. actinidiae

Pseudomonas syringae pv. actinidiae (Psa) causes bacterial canker, a devastating disease threatening the Actinidia fruit industry. In a search for non-host resistance genes against Psa, we found that the nucleotide-binding leucine-rich repeat receptor (NLR) protein ZAR1 from both Arabidopsis and Nicotiana benthamiana (Nb) recognizes HopZ5 and triggers cell death. The recognition requires ZED1 in Arabidopsis and JIM2 in Nb plants, which are members of the ZRK pseudokinases and known components of the ZAR1 resistosome. Surprisingly, Arabidopsis ZAR1 and RPM1, another NLR known to recognize HopZ5, confer disease resistance to HopZ5 in a strain-specific manner. Thus, ZAR1, but not RPM1, is solely required for resistance to P. s. maculicola ES4326 (Psm) carrying hopZ5, whereas RPM1 is primarily required for resistance to P. s. tomato DC3000 (Pst) carrying hopZ5. Furthermore, the ZAR1-mediated resistance to Psm hopZ5 in Arabidopsis is insensitive to SOBER1, which encodes a deacetylase known to suppress the RPM1-mediated resistance to Pst hopZ5. In addition, hopZ5 enhances P. syringae virulence in the absence of ZAR1 or RPM1, and that SOBER1 abolishes such virulence function. Together the study suggests that ZAR1 may be used for improving Psa resistance in Actinidia and uncovers previously unknown complexity of effector-triggered immunity and effector-triggered virulence.

Genome-Wide Identification of the LHC Gene Family in Kiwifruit and Regulatory Role of AcLhcb3.1/3.2 for Chlorophyll a Content

Light-harvesting chlorophyll a/b-binding (LHC) protein is a superfamily that plays a vital role in photosynthesis. However, the reported knowledge of LHCs in kiwifruit is inadequate and poorly understood. In this study, we identified 42 and 45 LHC genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes. Phylogenetic analysis showed that the kiwifruit LHCs of both species were grouped into four subfamilies (Lhc, Lil, PsbS, and FCII). Expression profiles and qRT-PCR results revealed expression levels of LHC genes closely related to the light, temperature fluctuations, color changes during fruit ripening, and kiwifruit responses to Pseudomonas syringae pv. actinidiae (Psa). Subcellular localization analysis showed that AcLhcb1.5/3.1/3.2 were localized in the chloroplast while transient overexpression of AcLhcb3.1/3.2 in tobacco leaves confirmed a significantly increased content of chlorophyll a. Our findings provide evidence of the characters and evolution patterns of kiwifruit LHCs genes in kiwifruit and verify the AcLhcb3.1/3.2 genes controlling the chlorophyll a content.

Kiwifruit MYBS1-like and GBF3 transcription factors influence l-ascorbic acid biosynthesis by activating transcription of GDP-L-galactose phosphorylase 3

Plant-derived Vitamin C (l-ascorbic acid (AsA)) is crucial for human health and wellbeing and thus increasing AsA content is of interest to plant breeders. In plants GDP-l-galactose phosphorylase (GGP) is a key biosynthetic control step and here evidence is presented for two new transcriptional activators of GGP. AsA measurement, transcriptomics, transient expression, hormone application, gene editing, yeast 1/2-hybrid, and electromobility shift assay (EMSA) methods were used to identify two positively regulating transcription factors. AceGGP3 was identified as the most highly expressed GGP in Actinidia eriantha fruit, which has high fruit AsA. A gene encoding a 1R-subtype myeloblastosis (MYB) protein, AceMYBS1, was found to bind the AceGGP3 promoter and activate its expression. Overexpression and gene-editing show AceMYBS1 effectively increases AsA accumulation. The bZIP transcription factor AceGBF3 (a G-box binding factor), also was shown to increase AsA content, and was confirmed to interact with AceMYBS1. Co-expression experiments showed that AceMYBS1 and AceGBF3 additively promoted AceGGP3 expression. Furthermore, AceMYBS1, but not GBF3, was repressed by abscisic acid, resulting in reduced AceGGP3 expression and accumulation of AsA. This study sheds new light on the roles of MYBS1 homologues and ABA in modulating AsA synthesis, and adds to the understanding of mechanisms underlying AsA accumulation.

The genome sequencing and comparative analysis of a wild kiwifruit Actinidia eriantha

The current kiwifruit industry is mainly based on the cultivars derived from the species Actinidia chinensis (Ac) which may bring risks such as canker disease. Introgression of desired traits from wild relatives is an important method for improving kiwifruit cultivars. Actinidia eriantha (Ae) is a particularly important taxon used for hybridization or introgressive breeding of new kiwifruit cultivars because of its valued species-specific traits. Here, we assembled a chromosome-scale high-quality genome of a Ae sample which was directly collected from its wild populations. Our analysis revealed that 41.3% of the genome consists of repetitive elements, comparable to the percentage in Ac and Ae cultivar "White" genomes. The genomic structural variation, including the presence/absence-variation (PAV) of genes, is distinct between Ae and Ac, despite both sharing the same two kiwifruit-specific whole genome duplication (WGD) events. This suggests that a post-WGD divergence mechanism occurred during their evolution. We further investigated genes involved in ascorbic acid biosynthesis and disease-resistance of Ae, and we found introgressive genome could contribute to the complex relationship between Ae and other representative kiwifruit taxa. Collectively, the Ae genome offers valuable genetic resource to accelerate kiwifruit breeding applications.

Genome-wide identification and characterization of the TIFY gene family in kiwifruit

Background

The TIFY gene family is a group of plant-specific transcription factors involved in regulation of plant growth and development and a variety of stress responses. However, the TIFY family has not yet been well characterized in kiwifruit, a popular fruit with important nutritional and economic value.

Results

A total of 27 and 21 TIFY genes were identified in the genomes of Actinidia eriantha and A. chinensis, respectively. Phylogenetic analyses showed that kiwifruit TIFY genes could be classified into four major groups, JAZ, ZML, TIFY and PPD, and the JAZ group could be further clustered into six subgroups (JAZ I to JAZ VI). Members within the same group or subgroup have similar exon-intron structures and conserved motif compositions. The kiwifruit TIFY genes are unevenly distributed on the chromosomes, and the segmental duplication events played a vital role in the expansion of the TIFY genes in kiwifruit. Syntenic analyses of TIFY genes between kiwifruit and other five plant species (including Arabidopsis thaliana, Camellia sinensis, Oryza sativa, Solanum lycopersicum and Vitis vinifera) and between the two kiwifruit species provided valuable clues for understanding the potential evolution of the kiwifruit TIFY family. Molecular evolutionary analysis showed that the evolution of kiwifruit TIFY genes was primarily constrained by intense purifying selection. Promoter cis-element analysis showed that most kiwifruit TIFY genes possess multiple cis-elements related to stress-response, phytohormone signal transduction and plant growth and development. The expression pattern analyses indicated that TIFY genes might play a role in different kiwifruit tissues, including fruit at specific development stages. In addition, several TIFY genes with high expression levels during Psa (Pseudomonas syringae pv. actinidiae) infection were identified, suggesting a role in the process of Pas infection.

Conclusions

In this study, the kiwifruit TIFY genes were identified from two assembled kiwifruit genomes. In addition, their basic physiochemical properties, chromosomal localization, phylogeny, gene structures and conserved motifs, synteny analyses, promoter cis-elements and expression patters were systematically examined. The results laid a foundation for further understanding the function of TIFY genes in kiwifruit, and provided a new potential approach for the prevention and treatment of Psa infection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08398-8.

Patterns of genomic divergence in sympatric and allopatric speciation of three Mihoutao (Actinidia) species

Isolation by geographic distance is involved in the formation of potential genomic islands and the divergence of genomes, which often result in speciation. The mechanisms of sympatric and allopatric speciation associated with geographic distance remain a topic of interest to evolutionary biologists. Here, we examined genomic divergence in three Actinidia species from large-scale sympatric and allopatric regions. Genome sequence data revealed that hexaploid Actinidia deliciosa originated from Actinidia chinensis and supported the speciation-with-gene-flow model in sympatric regions. The common ancestor of Actinidia setosa and A. deliciosa migrated from the mainland to the Taiwan Island ~2.91 Mya and formed A. setosa ~0.92 Mya, and the speciation of A. setosa is consistent with the divergence-after-speciation model with selective sweeps. Geographic isolation resulted in population contraction and accelerated the process of lineage sorting and speciation due to natural selection. Genomic islands contained genes associated with organ development, local adaptation, and stress resistance, indicating selective sweeps on a specific set of traits. Our results highlight the patterns of genomic divergence in sympatric and allopatric speciation, with the mediation of geographic isolation in the formation of genomic islands during Actinidia speciation.

A MADS-box gene with similarity to FLC is induced by cold and correlated with epigenetic changes to control budbreak in kiwifruit

Temperate perennials require exposure to chilling temperatures to resume growth in the following spring. Growth and dormancy cycles are controlled by complex genetic regulatory networks and are governed by epigenetic mechanisms, but the specific genes and mechanisms remain poorly understood. To understand how seasonal changes and chilling regulate dormancy and growth in the woody perennial vine kiwifruit (Ac, Actinidia chinensis), a transcriptome study of kiwifruit buds in the field and controlled conditions was performed. A MADS-box gene with homology to Arabidopsis FLOWERING LOCUS C (FLC) was identified and characterized. Elevated expression of AcFLC-like (AcFLCL) was detected during bud dormancy and chilling. A long noncoding (lnc) antisense transcript with an expression pattern opposite to AcFLCL and shorter sense noncoding RNAs were identified. Chilling induced an increase in trimethylation of lysine-4 of histone H3 (H3K4me3) in the 5' end of the gene, indicating multiple layers of epigenetic regulation in response to cold. Overexpression of AcFLCL in kiwifruit gave rise to plants with earlier budbreak, whilst gene editing using CRISPR-Cas9 resulted in transgenic lines with substantially delayed budbreak, suggesting a role in activation of growth. These results have implications for the future management and breeding of perennials for resilience to changing climate.

Characterization and Identification of a Ripening-Related Gene AaPG18 in Actinidia arguta

Actinidia arguta (A. arguta) is a kind of climacteric fruit that quickly softens and limits fruit shelf-life and commercial value. Therefore, it is of great significance to develop kiwifruit genotypes with an extended shelf-life of fruit. However, the ripening and softening mechanisms remain unclear in A. arguta. Here, we demonstrated that a key polygalacturonase (PG)-encoding gene AaPG18 was involved in A. arguta ripening through the degradation of the cell wall. Fruits were harvested at three developmental stages (S1, S2, and S3) for high-throughput transcriptome sequencing, based on which two candidate transcripts c109562_g1 and c111961_g1 were screened. The genome-wide identification of the PG gene family assigned c109562_g1 and c111961_g1 to correspond to AaPG4 and AaPG18, respectively. The expression profiles of candidate genes at six preharvest stages of fruit showed significantly higher expression levels of AaPG18 than AaPG4, indicating AaPG18 might be a key gene during fruit ripening processes. The subcellular localization displayed AaPG18 was located at the cytoplasmic membrane. The transient overexpression of AaPG18 in strawberry and the following morphological observation suggested AaPG18 played a key role in maintaining the stability of cell morphology. The homologous transient transformation in A. arguta “RB-4” proved the crucial function of AaPG18 in fruit ripening processes by causing the rapid redness of the fruit, which was an indicator of fruit maturity. All in all, our results identified AaPG18 as a key candidate gene involved in cell wall degeneration, which provides a basis for the subsequent exploration of the molecular mechanisms underlying the ripening and softening of A. arguta fruit.

GreeNC 2.0: a comprehensive database of plant long non-coding RNAs

The Green Non-Coding Database (GreeNC) is one of the reference databases for the study of plant long non-coding RNAs (lncRNAs). Here we present our most recent update where 16 species have been updated, while 78 species have been added, resulting in the annotation of more than 495 000 lncRNAs. Moreover, sequence clustering was applied providing information about sequence conservation and gene families. The current version of the database is available at: http://greenc.sequentiabiotech.com/wiki2/Main_Page.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

Transcription factors AcERF74/75 respond to waterlogging stress and trigger alcoholic fermentation-related genes in kiwifruit

Kiwifruit plants have a fleshy, shallow root system which is sensitive to waterlogging stress, which results in a decrease in crop yield or even plants death. Although the waterlogging stress responses in kiwifruit have attracted much attention, the underlying molecular mechanism remains unclear. In this study, waterlogging led to drastic inhibition of root growth of 'Donghong' kiwifruit (Actinidia chinensis) plants grown in vitro, which was accompanied by significant elevation of endogenous acetaldehyde and ethanol contents. RNA-seq of roots of plants waterlogged for 0, 1 and 2 days revealed that a total of 149 genes were up- or down-regulated, including seven biosynthetic genes related to the glycolysis/gluconeogenesis pathway and 10 transcription factors. Analyses with real-time PCR, dual-luciferase assays and EMSA demonstrated that AcERF74 and AcERF75, two members of the ERF-VII subfamily, directly upregulated AcADH1 (alcohol dehydrogenase). Moreover, the overexpression of AcERF74/75 in transgenic calli resulted in dramatic increase of endogenous ethanol contents through the triggering of AcADH1 and AcADH2 expression. Although the AcPDC2 (pyruvate decarboxylase) expression was also enhanced in transgenic lines, the endogenous acetaldehyde contents showed no significant changes. These results illustrated that AcERF74/75 are two transcriptional activators on alcoholic fermentation related genes and are responsive to waterlogging stress in kiwifruit.