Kiwifruit Information Resource (KIR): a comparative platform for kiwifruit genomics

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.

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.

Novel insight into ferroptosis-related genes, molecular subtypes, and immune characteristics in intracranial aneurysms

OBJECTIVES

This study aimed to identify the role of ferroptosis in intracranial aneurysm (IA).

METHODS

GSE122897, GSE75436, GSE15629, and GSE75434 datasets were downloaded from the Gene Expression Omnibus database. The differentially expressed ferroptosis-related genes (DEFRGs) were selected to construct a diagnostic model integrating with machine learning. Then, a consensus clustering algorithm was performed to classify IA patients into distinct ferroptosis-related clusters. Functional analyses, including GO, KEGG, GSVA, and GSEA analyses, were conducted to elucidate the underlying mechanisms. ssGSEA and xCell algorithms were performed to uncover the immune characteristics.

RESULTS

We identified 28 DEFRGs between IAs and controls from the GSE122897 dataset. GO and KEGG results showed that these genes were enriched in cytokine activity, ferroptosis, and the IL-17 signaling pathway. Immune analysis showed that the IAs had higher levels of immune infiltration. A four FRGs model (MT3, CDKN1A, ZEP69B, and ABCC1) was established and validated with great IA diagnostic ability. We divided the IA samples into two clusters and found that cluster 2 had a higher proportion of rupture and immune infiltration. We identified 10 ferroptosis phenotypes-related markers in IAs.

CONCLUSION

Ferroptosis and the immune microenvironment are closely associated with IAs, providing a basis for understanding the IA development.

Metabolic modeling-based drug repurposing in Glioblastoma

The manifestation of intra- and inter-tumor heterogeneity hinders the development of ubiquitous cancer treatments, thus requiring a tailored therapy for each cancer type. Specifically, the reprogramming of cellular metabolism has been identified as a source of potential drug targets. Drug discovery is a long and resource-demanding process aiming at identifying and testing compounds early in the drug development pipeline. While drug repurposing efforts (i.e., inspecting readily available approved drugs) can be supported by a mechanistic rationale, strategies to further reduce and prioritize the list of potential candidates are still needed to facilitate feasible studies. Although a variety of ‘omics’ data are widely gathered, a standard integration method with modeling approaches is lacking. For instance, flux balance analysis is a metabolic modeling technique that mainly relies on the stoichiometry of the metabolic network. However, exploring the network’s topology typically neglects biologically relevant information. Here we introduce Transcriptomics-Informed Stoichiometric Modelling And Network analysis (TISMAN) in a recombinant innovation manner, allowing identification and validation of genes as targets for drug repurposing using glioblastoma as an exemplar.

A Comparative Analysis of the Microbiome of Kiwifruit at Harvest Under Open-Field and Rain-Shelter Cultivation Systems

The composition of microbial communities can directly affect fruit quality, health status, and storability. The present study characterized the epiphytes and endophytes of “Hongyang” and “Cuiyu” kiwifruit at harvest under grown under open-field (OF) and rain-shelter (RS) cultivation systems. Disease incidence in kiwifruit was significantly lower (p < 0.05) under the RS system than it was under the OF system. High-throughput sequencing [16S V3-V4 ribosomal region and the fungal internal transcribed spacer (ITS2)] was conducted to compare the composition of the epiphytic and endophytic microbial community of kiwifruit under the two cultivation systems. Results indicated that the abundance of Actinobacteria, Bacteroidetes, Enterobacteriales, Acetobacterales, Sphingomonas, Pseudomonas, and Sphingobacterium was higher under the RS system, relative to the OF system, while the abundance of Capnodiales, Hypocreales, Vishniacozyma, and Plectosphaerella was also higher under the RS system. Some of these bacterial and fungal taxa have been reported to as act as biocontrol agents and reduce disease incidence. Notably, the α-diversity of the epiphytic bacterial and fungal communities on kiwifruit was higher under RS cultivation. In summary, RS cultivation reduced natural disease incidence in kiwifruit, which may be partially attributed to differences in the structure and composition of the microbial community present in and on kiwifruit.

Identification of GH17 gene family in Vitis vinifera and expression analysis of GH17 under various adversities

Glycoside hydrolase (GH, EC 3.2.1) is a group of enzymes that hydrolyzes glycosidic bonds and play a role in the hydrolysis and synthesis of sugars in living organisms. Vitis vinifera is an important fruit crop and it harbors GH17 gene family however, their function in grapes has not been systematically investigated. In this study, a total of 870 GH17 genes were identified from 14 plant species and their structural domain, sequence alignment, phylogenetic tree, collinear analysis, with the expression profiles of VviGH17 gene family was performed. The promoter analysis of VviGH17 gene showed the presence of cis-acting elements, which are responsive to plant growth and development. In addition, elements for plant hormones were found that are triggered in response to abiotic/biological stress. Transcriptomic data led to the identification of several VviGH17 genes, which are associated with bud dormancy and in response to abiotic stress. Transcript analysis was carried out for some of the selected VviGH17 genes RT-qPCR. VviGH17-16 and VviGH17-30 genes were differentially expressed during bud dormancy, fruit development and different abiotic stresses. Moreover, VviGH17-37 and VviGH17-44 were differentially expressed at fruit development, in response to abiotic stress. In addition, subcellular localization predicts that the VviGH17-16, VviGH17-30, and VviGH17-37 genes were located in the cell membrane, while VviGH17-44 gene was located in the vacuole. In conclusion, our study led to the identification of several GH17s and their probable role in development and stress.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01014-1.

Comparative Analysis of Physicochemical Characteristics, Nutritional and Functional Components and Antioxidant Capacity of Fifteen Kiwifruit (Actinidia) Cultivars-Comparative Analysis of Fifteen Kiwifruit (Actinidia) Cultivars

Physicochemical characteristics, nutritional and functional components, and the antioxidant capacity of 15 kinds of domestic and imported kiwifruit in China were studied. Kiwifruit was classified according to flesh color or species, and the differences were analyzed and compared. Results demonstrated Ruiyu had the highest sugar-acid ratio, and Hongshi No.2 was an excellent cultivar with strong antioxidant capacity. TPC (total polyphenol content) and AAC (ascorbic acid content) showed a significant positive correlation. TPC was the greatest antioxidant contributor in the DPPH and FRAP assays. The sugar-acid ratio and TFC (total flavonoids content) in red-fleshed kiwifruit were significantly higher than those in yellow-fleshed and green-fleshed ones. The composition of free amino acids had a tendency to distinguish A. deliciosa and A. chinensis, but this needs further verification. In addition, the contents of mineral elements, folic acid and L-5-methyltetrahydrofolate were also analyzed. Generally, kiwifruit contains comprehensive nutrients and has strong antioxidant capacity. Cultivar is one of the main factors affecting nutritional and functional properties and antioxidant capacity.

Kiwifruit Genome Database (KGD): a comprehensive resource for kiwifruit genomics

Kiwifruit (Actinidia spp.) plants produce economically important fruits containing abundant, balanced phytonutrients with extraordinarily high vitamin C contents. Since the release of the first kiwifruit reference genome sequence in 2013, large volumes of genome and transcriptome data have been rapidly accumulated for a handful of kiwifruit species. To efficiently store, analyze, integrate, and disseminate these large-scale datasets to the research community, we constructed the Kiwifruit Genome Database (KGD; http://kiwifruitgenome.org/). The database currently contains all publicly available genome and gene sequences, gene annotations, biochemical pathways, transcriptome profiles derived from public RNA-Seq datasets, and comparative genomic analysis results such as syntenic blocks and homologous gene pairs between different kiwifruit genome assemblies. A set of user-friendly query interfaces, analysis tools and visualization modules have been implemented in KGD to facilitate translational and applied research in kiwifruit, which include JBrowse, a popular genome browser, and the NCBI BLAST sequence search tool. Other notable tools developed within KGD include a genome synteny viewer and tools for differential gene expression analysis as well as gene ontology (GO) term and pathway enrichment analysis.

Full-length transcriptome sequencing provides insights into the evolution of apocarotenoid biosynthesis in Crocus sativus

Crocus sativus, containing remarkably amounts of crocin, picrocrocin and safranal, is the source of saffron with tremendous medicinal, economic and cultural importance. Here, we present a high-quality full-length transcriptome of the sterile triploid C. sativus, using the PacBio SMRT sequencing technology. This yields 31,755 high-confidence predictions of protein-coding genes, with 50.1% forming paralogous gene pairs. Analysis on distribution of Ks values suggests that the current genome of C. sativus is probably a product resulting from at least two rounds of whole-genome duplication (WGD) events occurred at ~28 and ~114 million years ago (Mya), respectively. We provide evidence demonstrating that the recent β WGD event confers a major impact on family expansion of secondary metabolite genes, possibly leading to an enhanced accumulation of three distinct compounds: crocin, picrocrocin and safranal. Phylogenetic analysis unravels that the founding member (CCD2) of CCD enzymes necessary for the biosynthesis of apocarotenoids in C. sativus might be evolved from the CCD1 family via the β WGD event. Based on the gene expression profiling, CCD2 is found to be expressed at an extremely high level in the stigma. These findings may shed lights on further genomic refinement of the characteristic biosynthesis pathways and promote germplasm utilization for the improvement of saffron quality.

A MYB/bHLH complex regulates tissue-specific anthocyanin biosynthesis in the inner pericarp of red-centered kiwifruit Actinidia chinensis cv. Hongyang

Many Actinidia cultivars are characterized by anthocyanin accumulation, specifically in the inner pericarp, but the underlying regulatory mechanism remains elusive. Here we report two interacting transcription factors, AcMYB123 and AcbHLH42, that regulate tissue-specific anthocyanin biosynthesis in the inner pericarp of Actinidia chinensis cv. Hongyang. Through transcriptome profiling analysis we identified five MYB and three bHLH transcription factors that were upregulated in the inner pericarp. We show that the combinatorial action of two of them, AcMYB123 and AcbHLH42, is required for activating promoters of AcANS and AcF3GT1 that encode the dedicated enzymes for anthocyanin biosynthesis. The presence of anthocyanin in the inner pericarp appears to be tightly associated with elevated expression of AcMYB123 and AcbHLH42. RNA interference repression of AcMYB123, AcbHLH42, AcF3GT1 and AcANS in 'Hongyang' fruits resulted in significantly reduced anthocyanin biosynthesis. Using both transient assays in Nicotiana tabacum leaves or Actinidia arguta fruits and stable transformation in Arabidopsis, we demonstrate that co-expression of AcMYB123 and AcbHLH42 is a prerequisite for anthocyanin production by activating transcription of AcF3GT1 and AcANS or the homologous genes. Phylogenetic analysis suggests that AcMYB123 or AcbHLH42 are closely related to TT2 or TT8, respectively, which determines proanthocyanidin biosynthesis in Arabidopsis, and to anthocyanin regulators in monocots rather than regulators in dicots. All these experimental results suggest that AcMYB123 and AcbHLH42 are the components involved in spatiotemporal regulation of anthocyanin biosynthesis specifically in the inner pericarp of kiwifruit.

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

Background

Most published genome sequences are drafts, and most are dominated by computational gene prediction. Draft genomes typically incorporate considerable sequence data that are not assigned to chromosomes, and predicted genes without quality confidence measures. The current Actinidia chinensis (kiwifruit) ‘Hongyang’ draft genome has 164 Mb of sequences unassigned to pseudo-chromosomes, and omissions have been identified in the gene models.

Results

A second genome of an A. chinensis (genotype Red5) was fully sequenced. This new sequence resulted in a 554.0 Mb assembly with all but 6 Mb assigned to pseudo-chromosomes. Pseudo-chromosomal comparisons showed a considerable number of translocation events have occurred following a whole genome duplication (WGD) event some consistent with centromeric Robertsonian-like translocations. RNA sequencing data from 12 tissues and ab initio analysis informed a genome-wide manual annotation, using the WebApollo tool. In total, 33,044 gene loci represented by 33,123 isoforms were identified, named and tagged for quality of evidential support. Of these 3114 (9.4%) were identical to a protein within ‘Hongyang’ The Kiwifruit Information Resource (KIR v2). Some proportion of the differences will be varietal polymorphisms. However, as most computationally predicted Red5 models required manual re-annotation this proportion is expected to be small. The quality of the new gene models was tested by fully sequencing 550 cloned ‘Hort16A’ cDNAs and comparing with the predicted protein models for Red5 and both the original ‘Hongyang’ assembly and the revised annotation from KIR v2. Only 48.9% and 63.5% of the cDNAs had a match with 90% identity or better to the original and revised ‘Hongyang’ annotation, respectively, compared with 90.9% to the Red5 models.

Conclusions

Our study highlights the need to take a cautious approach to draft genomes and computationally predicted genes. Our use of the manual annotation tool WebApollo facilitated manual checking and correction of gene models enabling improvement of computational prediction. This utility was especially relevant for certain types of gene families such as the EXPANSIN like genes. Finally, this high quality gene set will supply the kiwifruit and general plant community with a new tool for genomics and other comparative analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4656-3) contains supplementary material, which is available to authorized users.

PCPPI: a comprehensive database for the prediction of Penicillium-crop protein-protein interactions

Penicillium expansum , the causal agent of blue mold, is one of the most prevalent post-harvest pathogens, infecting a wide range of crops after harvest. In response, crops have evolved various defense systems to protect themselves against this and other pathogens. Penicillium -crop interaction is a multifaceted process and mediated by pathogen- and host-derived proteins. Identification and characterization of the inter-species protein-protein interactions (PPIs) are fundamental to elucidating the molecular mechanisms underlying infection processes between P. expansum and plant crops. Here, we have developed PCPPI, the Penicillium -Crop Protein-Protein Interactions database, which is constructed based on the experimentally determined orthologous interactions in pathogen-plant systems and available domain-domain interactions (DDIs) in each PPI. Thus far, it stores information on 9911 proteins, 439 904 interactions and seven host species, including apple, kiwifruit, maize, pear, rice, strawberry and tomato. Further analysis through the gene ontology (GO) annotation indicated that proteins with more interacting partners tend to execute the essential function. Significantly, semantic statistics of the GO terms also provided strong support for the accuracy of our predicted interactions in PCPPI. We believe that all the PCPPI datasets are helpful to facilitate the study of pathogen-crop interactions and freely available to the research community.

Database URL

: http://bdg.hfut.edu.cn/pcppi/index.html.

Identification of Conserved and Novel MicroRNAs in Blueberry

MicroRNAs (miRNAs) are a class of small endogenous RNAs that play important regulatory roles in cells by negatively affecting gene expression at both transcriptional and post-transcriptional levels. There have been extensive studies aiming to identify miRNAs and to elucidate their functions in various plant species. In the present study, we employed the high-throughput sequencing technology to profile miRNAs in blueberry fruits. A total of 9,992,446 small RNA tags with sizes ranged from 18 to 30 nt were obtained, indicating that blueberry fruits have a large and diverse small RNA population. Bioinformatic analysis identified 412 conserved miRNAs belonging to 29 families, and 35 predicted novel miRNAs that are likely to be unique to blueberries. Among them, expression profiles of five conserved miRNAs were validated by stem loop qRT-PCR. Furthermore, the potential target genes of conserved and novel miRNAs were predicted and subjected to Gene Ontology (GO) annotation. Enrichment analysis of the GO-represented biological processes and molecular functions revealed that these target genes were potentially involved in a wide range of metabolic pathways and developmental processes. Particularly, anthocyanin biosynthesis has been predicted to be directly or indirectly regulated by diverse miRNA families. This study is the first report on genome-wide miRNA profile analysis in blueberry and it provides a useful resource for further elucidation of the functional roles of miRNAs during fruit development and ripening.