AppleMDO: A Multi-Dimensional Omics Database for Apple Co-Expression Networks and Chromatin States

Immunoproteomic profile of Malus domestica in Mexican pediatric patients. Evidence of new allergen prospects

Background: Apple (Malus domestica) is a fruit commonly associated with allergic oral symptoms in the Mexican pediatric population; however, knowledge of its allergenic proteins is limited. This information is crucial as sensitization frequencies to specific allergens can vary among different populations. The main allergic symptomatology before apple ingestion derives from primary sensitizations induced by pollen, promoting cross-reactivity with the main allergenic protein of apple. Therefore, this study aims to identify new potential sensitizing proteins to apple using immunoproteomic techniques. Methods: We collected serum samples from 14 pediatric patients with confirmed immunoglobulin E (IgE)-mediated apple allergy and used these samples to assess immunoreactivity to apple protein extracts through 2D-western blot assays. The spots corresponding to the 2D-SDS-PAGE were analyzed using nanoLC-MS/MS. Results: We identified 11 non-redundant proteins, including Mal d 2 and Mal d 1, the latter showing a high frequency of sensitization (79%) in our patients, and being considered the main apple allergenic protein. The remaining identified proteins have not been previously described as apple allergens in the International Union of Immunological Societies databases. However, three of these may be categorized as pan-allergens. Conclusions: This study shows evidence that the repertoire of apple allergens in the Mexican population could differ from those reported internationally, highlighting the importance of studies in different countries to improve the certainty of allergy diagnosis and allow the implementation of precision medicine.

HortGenome Search Engine, a universal genomic search engine for horticultural crops

Horticultural crops comprising fruit, vegetable, ornamental, beverage, medicinal and aromatic plants play essential roles in food security and human health, as well as landscaping. With the advances of sequencing technologies, genomes for hundreds of horticultural crops have been deciphered in recent years, providing a basis for understanding gene functions and regulatory networks and for the improvement of horticultural crops. However, these valuable genomic data are scattered in warehouses with various complex searching and displaying strategies, which increases learning and usage costs and makes comparative and functional genomic analyses across different horticultural crops very challenging. To this end, we have developed a lightweight universal search engine, HortGenome Search Engine (HSE; http://hort.moilab.net), which allows for the querying of genes, functional annotations, protein domains, homologs, and other gene-related functional information of more than 500 horticultural crops. In addition, four commonly used tools, including 'BLAST', 'Batch Query', 'Enrichment analysis', and 'Synteny Viewer' have been developed for efficient mining and analysis of these genomic data.

Genome-wide identification and expression analysis of protein arginine methyltransferase and JmjC domain-containing family in apple

Histone methylation is an important type of histone modification that regulates gene expression in plants. In this study, we identified 14 arginine methylation-related genes (Protein Arginine Methyltransferase, MdPRMT) and 32 demethylation-related genes (JmjC Domain-Containing Family, MdJMJ) in apple. Furthermore, we investigated the phylogenetic relationship, chromosome distribution, gene structure, motif analysis, promoter sequence analysis, and expression patterns of MdPRMT and MdJMJ genes. Homology analysis showed a high degree of conservation and homology between PRMT and JMJ genes in Arabidopsis and apple. We identified the types of duplicated genes in the MdJMJ and MdPRMT gene families, found a large number of whole-genome duplicates (WGD) gene pairs and a small number of tandem duplicates (TD) pairs, transposed duplication (TRD) gene pairs as well as proximal duplicates (PD) pairs, and discussed the possible evolutionary pathways of the gene families from the perspective of duplicated genes. Homology analysis showed a high degree of conservation and homology between PRMT and JMJ genes in Arabidopsis and apple. In addition, the promoter regions of MdPRMT and MdJMJ contain numerous cis-acting elements involved in plant growth and development, hormone response, and stress responses. Based on the transcriptional profiles of MdPRMT and MdJMJ in different tissues and developmental stages, it was found that MdPRMT and MdJMJ may play multiple roles in apple growth and development, for example, MdJMJ21 may be involved in the regulation of apple endosperm formation. MdPRMT and MdJMJ exhibit different expression patterns in response to hormone signaling in apple, MdJMJ3, MdJMJ18, MdJMJ30, MdPRMT2, MdPRMT13, and MdPRMT14 may play roles in apple response to drought stress, while the expression of MdJMJ13, MdPRMT3, MdPRMT4, and MdPRMT6 is affected by cold stress. Our study provides a foundation for determining the functional roles of MdPRMT and MdJMJ genes in apple.

RoseAP: an analytical platform for gene function of Rosa rugosa

Rosa rugosa, a perennial shrub belonging to family Rosaceae, is a well-known ornamental plant. Its petals contain an abundance of essential oils and anthocyanins with enormous economic and health benefits when used as edible or cosmetic ingredients. The whole genome of R. rugosa was sequenced in 2021, which provided opportunities and challenges for gene regulation. However, many gene functions remain unknown. Therefore, an analytical platform named RoseAP (http://www.gzybioinformatics.cn/RoseAP/index.php) for the functional analysis of R. rugosa genes was constructed. It improved the gene annotation rate by integrating and analyzing genomic and transcriptomic datasets. First, 38,815 genes, covering 97.76% of the coding genes, were annotated functionally and structurally using a variety of algorithms and rules. Second, a total of 33 transcriptome samples were integrated, including 23 samples from our lab and 10 samples from the SRA database. A co-expression network containing approximately 29,657 positive or negative gene pairs, covering 74.7% of the coding genes, was constructed based on PCC and MR algorithms. Network analysis revealed that the DFR function was closely related to anthocyanin metabolism. It demonstrated the reliability of the network. Several SAUR genes of R. rugosa shared similar expression patterns. RoseAP was used to determine the sequence, structure, functional annotation, expression profile, regulatory network, and functional modules at the transcriptional and protein levels by inputting gene IDs. In addition, auxiliary analytical tools, including BLAST, gene set enrichment, orthologue conversion, gene sequence extraction, gene expression value extraction, and JBrowse, were utilized. Regular updates to RoseAP are expected to facilitate mining of gene function and promote genetic improvement in R. rugosa.

Genome-wide analysis of MdPLATZ genes and their expression during axillary bud outgrowth in apple (Malus domestica Borkh.)

BACKGROUND

Branching is a plastic character that affects plant architecture and spatial structure. The trait is controlled by a variety of plant hormones through coordination with environmental signals. Plant AT-rich sequence and zinc-binding protein (PLATZ) is a transcription factor that plays an important role in plant growth and development. However, systematic research on the role of the PLATZ family in apple branching has not been conducted previously.

RESULTS

In this study, a total of 17 PLATZ genes were identified and characterized from the apple genome. The 83 PLATZ proteins from apple, tomato, Arabidopsis, rice, and maize were classified into three groups based on the topological structure of the phylogenetic tree. The phylogenetic relationships, conserved motifs, gene structure, regulatory cis-acting elements, and microRNAs of the MdPLATZ family members were predicted. Expression analysis revealed that MdPLATZ genes exhibited distinct expression patterns in different tissues. The expression patterns of the MdPLATZ genes were systematically investigated in response to treatments that impact apple branching [thidazuron (TDZ) and decapitation]. The expression of MdPLATZ1, 6, 7, 8, 9, 15, and 16 was regulated during axillary bud outgrowth based on RNA-sequencing data obtained from apple axillary buds treated by decapitation or exogenous TDZ application. Quantitative real-time PCR analysis showed that MdPLATZ6 was strongly downregulated in response to the TDZ and decapitation treatments, however, MdPLATZ15 was significantly upregulated in response to TDZ, but exhibited little response to decapitation. Furthermore, the co-expression network showed that PLATZ might be involved in shoot branching by regulating branching-related genes or mediating cytokinin or auxin pathway.

CONCLUSION

The results provide valuable information for further functional investigation of MdPLATZ genes in the control of axillary bud outgrowth in apple.

Genome-wide analysis of the MdZR gene family revealed MdZR2.2-induced salt and drought stress tolerance in apple rootstock

The DNL-type zinc finger protein constitutes a zinc ribbon protein (ZR) family, which belongs to a branch of zinc finger protein and plays an essential role in response to abiotic stress. Here, we identified six apple (Malus domestica) MdZR genes. Based on their phylogenetic relationship and gene structure, the MdZR genes were divided into three categories, including MdZR1, MdZR2, and MdZR3. Subcellular results showed that the MdZRs are located on the nuclear and membrane. The transcriptome data showed that MdZR2.2 is expressed in various tissues. The expression analysis results showed that MdZR2.2 was significantly upregulated under salt and drought treatments. Thus, we selected MdZR2.2 for further research. Overexpression of MdZR2.2 in apple callus improved their tolerance to drought and salt stress and ability to scavenge reactive oxygen species (ROS). In contrast, transgenic apple roots with silenced MdZR2.2 grew more poorly than the wild type when subjected to salt and drought stress, which reduced their ability to scavenge ROS. To our knowledge, this is the first study to analyze the MdZR protein family. This study identified a gene that responds to drought and salt stress. Our findings lay a foundation for a comprehensive analysis of the MdZR family members.

Kinase MxMPK4-1 and calmodulin-binding protein MxIQM3 enhance apple root acidification during Fe deficiency

Iron (Fe) deficiency is a long-standing issue in plant mineral nutrition. Ca2+ signals and the mitogen-activated protein kinase (MAPK) cascade are frequently activated in parallel to perceive external cues, but their interplay under Fe deficiency stress remains largely unclear. Here, the kinase MxMPK4-1, which is induced during the response to Fe deficiency stress in apple rootstock Malus xiaojinensis, cooperates with IQ-motif containing protein3 (MxIQM3). MxIQM3 gene expression, protein abundance, and phosphorylation level increased under Fe deficiency stress. The overexpression of MxIQM3 in apple calli and rootstocks mitigated the Fe deficiency phenotype and improved stress tolerance, whereas RNA interference or silencing of MxIQM3 in apple calli and rootstocks, respectively, worsened the phenotype and reduced tolerance to Fe deficiency. MxMPK4-1 interacted with MxIQM3 and subsequently phosphorylated MxIQM3 at Ser393, and co-expression of MxMPK4-1 and MxIQM3 in apple calli and rootstocks enhanced Fe deficiency responses. Furthermore, MxIQM3 interacted with the central-loop region of the plasma membrane (PM) H+-ATPase MxHA2. Phospho-mimicking mutation of MxIQM3 at Ser393 inhibited binding to MxHA2, but phospho-abolishing mutation promoted interaction with both the central-loop and C terminus of MxHA2, demonstrating phosphorylation of MxIQM3 caused dissociation from MxHA2 and therefore increased H+ secretion. Moreover, Ca2+/MxCAM7 (Calmodulin7) regulated the MxMPK4-1-MxIQM3 module in response to Fe deficiency stress. Overall, our results demonstrate that MxMPK4-1-MxIQM3 forms a functional complex and positively regulates PM H+-ATPase activity in Fe deficiency responses, revealing a versatile mechanism of Ca2+/MxCAM7 signaling and MAPK cascade under Fe deficiency stress.

Identification of PLATZ genes in Malus and expression characteristics of MdPLATZs in response to drought and ABA stresses

Plant AT-rich sequences and zinc-binding proteins (PLATZ) play crucial roles in response to environmental stresses. Nevertheless, PLATZ gene family has not been systemically studied in Rosaceae species, such as in apple, pear, peach, or strawberry. In this study, a total of 134 PLATZ proteins were identified from nine Rosaceae genomes and were classified into seven phylogenetic groups. Subsequently, the chromosomal localization, duplication, and collinearity relationship for apple PLATZ genes were investigated, and segmental duplication is a major driving-force in the expansion of PLATZ in Malus. Expression profiles analysis showed that PLATZs had distinct expression patterns in different tissues, and multiple genes were significantly changed after drought and ABA treatments. Furthermore, the co-expression network combined with RNA-seq data showed that PLATZ might be involved in drought stress by regulating ABA signaling pathway. In summary, this study is the first in-depth and systematic identification of PLATZ gene family in Rosaceae species, especially for apple, and provided specific PLATZ gene resource for further functional research in response to abiotic stress.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

MbHY5-MbYSL7 mediates chlorophyll synthesis and iron transport under iron deficiency in Malus baccata

Iron (Fe) plays an important role in cellular respiration and catalytic reactions of metalloproteins in plants and animals. Plants maintain iron homeostasis through absorption, translocation, storage, and compartmentalization of iron via a cooperative regulative network. Here, we showed different physiological characteristics in the leaves and roots of Malus baccata under Fe sufficiency and Fe deficiency conditions and propose that MbHY5 (elongated hypocotyl 5), an important transcription factor for its function in photomorphogenesis, participated in Fe deficiency response in both the leaves and roots of M. baccata. The gene co-expression network showed that MbHY5 was involved in the regulation of chlorophyll synthesis and Fe transport pathway under Fe-limiting conditions. Specifically, we found that Fe deficiency induced the expression of MbYSL7 in root, which was positively regulated by MbHY5. Overexpressing or silencing MbYSL7 influenced the expression of MbHY5 in M. baccata.

Analyzing lignin biosynthesis pathways in rattan using improved co-expression networks of NACs and MYBs

BACKGROUND

The rattan is a valuable plant resource with multiple applications in tropical forests. Calamus simplicifolius and Daemonorops jenkinsiana are the two most representative rattan species, supplying over 95% of the raw materials for the rattan industry. Hence, the wood properties of both rattans have always attracted researchers' attention.

RESULTS

We re-annotated the genomes, obtained 81 RNA-Seq datasets, and developed an improved pipeline to increase the reliability of co-expression networks of both rattans. Based on the data and pipeline, co-expression relationships were detected in 11 NACs, 49 MYBs, and 86 lignin biosynthesis genes in C. simplicifolius and four NACs, 59 MYBs, and 76 lignin biosynthesis genes in D. jenkinsiana, respectively. Among these co-expression pairs, several genes had a close relationship to the development of wood properties. Additionally, we detected the enzyme gene on the lignin biosynthesis pathway was regulated by either NAC or MYB, while LACCASES was regulated by both NAC and MYB. For D. jenkinsiana, the lignin biosynthesis regulatory network was characterized by positive regulation, and MYB possible negatively regulate non-expressed lignin biosynthesis genes in stem tissues. For C. simplicifolius, NAC may positively regulate highly expressed genes and negatively regulate non-expressed lignin biosynthesis genes in stem tissues. Furthermore, we established core regulatory networks of NAC and MYB for both rattans.

CONCLUSIONS

This work improved the accuracy of rattan gene annotation by integrating an efficient co-expression network analysis pipeline, enhancing gene coverage and accuracy of the constructed network, and facilitating an understanding of co-expression relationships among NAC, MYB, and lignin biosynthesis genes in rattan and other plants.

A Genome-Wide Analysis of StTGA Genes Reveals the Critical Role in Enhanced Bacterial Wilt Tolerance in Potato During Ralstonia solanacearum Infection

TGA is one of the members of TGACG sequence-specific binding protein family, which plays a crucial role in the regulated course of hormone synthesis as a stress-responsive transcription factor (TF). Little is known, however, about its implication in response to bacterial wilt disease in potato (Solanum tuberosum) caused by Ralstonia solanacearum. Here, we performed an in silico identification and analysis of the members of the TGA family based on the whole genome data of potato. In total, 42 StTGAs were predicted to be distributed on four chromosomes in potato genome. Phylogenetic analysis showed that the proteins of StTGAs could be divided into six sub-families. We found that many of these genes have more than one exon according to the conserved motif and gene structure analysis. The heat map inferred that StTGAs are generally expressed in different tissues which are at different stages of development. Genomic collinear analysis showed that there are homologous relationships among potato, tomato, pepper, Arabidopsis, and tobacco TGA genes. Cis-element in silico analysis predicted that there may be many cis-acting elements related to abiotic and biotic stress upstream of StTGA promoter including plant hormone response elements. A representative member StTGA39 was selected to investigate the potential function of the StTGA genes for further analysis. Quantitative real-time polymerase chain reaction (qRT-PCR) assays indicated that the expression of the StTGAs was significantly induced by R. solanacearum infection and upregulated by exogenous salicylic acid (SA), abscisic acid (ABA), gibberellin 3 (GA₃), and methyl jasmonate (MeJA). The results of yeast one-hybrid (Y1H) assay showed that StTGA39 regulates S. tuberosum BRI1-associated receptor kinase 1 (StBAK1) expression. Thus, our study provides a theoretical basis for further research of the molecular mechanism of the StTGA gene of potato tolerance to bacterial wilt.

Genome-wide identification and comparative analysis of genes encoding AAPs in apple (Malus × domestica Borkh.)

Amino acid permeases (AAPs) play important roles in plant amino acid transport and nitrogen metabolism. In this study, we carried a comprehensive analysis for apple genes encoding AAPs using bioinformatics and molecular biology. Eleven MdAAPs were identified by a genome-wide search and comparative genomic analysis revealed relatively conserved gene composition, transmembrane characteristics, and protein structures. Phylogenetic tree construction and analysis of the conserved motifs of MdAAPs and AtAAPs showed that AAPs can be classified into three groups (I, II, and III). We compared the promoters of the identified genes and did gene functional annotation and qRT-PCR and found a relationship between apple AAPs and nitrogen deficiency. The expression profile data implied that MdAAPs exhibit diversified distributions and functions in different tissues.

Gene Co-Expression Network Tools and Databases for Crop Improvement

Transcriptomics has significantly grown as a functional genomics tool for understanding the expression of biological systems. The generated transcriptomics data can be utilised to produce a gene co-expression network that is one of the essential downstream omics data analyses. To date, several gene co-expression network databases that store correlation values, expression profiles, gene names and gene descriptions have been developed. Although these resources remain scattered across the Internet, such databases complement each other and support efficient growth in the functional genomics area. This review presents the features and the most recent gene co-expression network databases in crops and summarises the present status of the tools that are widely used for constructing the gene co-expression network. The highlights of gene co-expression network databases and the tools presented here will pave the way for a robust interpretation of biologically relevant information. With this effort, the researcher would be able to explore and utilise gene co-expression network databases for crops improvement.

Comprehensive Genome-Wide Identification and Transcript Profiling of GABA Pathway Gene Family in Apple (Malus domestica)

γ-Aminobutyric Acid (GABA), a four-carbon non-protein amino acid, is a significant component of the free amino acid pool in most prokaryotic and eukaryotic organisms. GABA is involved in pH regulation, maintaining C/N balance, plant development and defence, as well as a compatible osmolyte and an alternative pathway for glutamate utilization via anion flux. Glutamate decarboxylase (GAD, EC 4.1.1.15) and GABA transaminase (GABA-T, EC 2.6.1.19) are two key enzymes involved in the synthesis and metabolism of GABA. Recently, GABA transporters (GATs), protein and aluminium-activated malate transporter (ALMT) proteins which function as GABA receptors, have been shown to be involved in GABA regulation. However, there is no report on the characterization of apple GABA pathway genes. In this study, we performed a genome-wide analysis and expression profiling of the GABA pathway gene family in the apple genome. A total of 24 genes were identified including five GAD genes (namely MdGAD 1–5), two GABA-T genes (namely MdGABA-T 1,2), 10 GAT genes (namely GAT 1–10) and seven ALMT genes (namely MdALMT1–7). These genes were randomly distributed on 12 chromosomes. Phylogenetic analyses grouped GABA shunt genes into three clusters—cluster I, cluster II, and cluster III—which had three, four, and five genes, respectively. The expression profile analysis revealed significant MdGAD4 expression levels in both fruit and flower organs, except pollen. However, there were no significant differences in the expression of other GABA shunt genes in different tissues. This work provides the first characterization of the GABA shunt gene family in apple and suggests their importance in apple response to abiotic stress. These results can serve as a guide for future studies on the understanding and functional characterization of these gene families.

Starch branching enzymes as putative determinants of postharvest quality in horticultural crops

Starch branching enzymes (SBEs) are key determinants of the structure and amount of the starch in plant organs, and as such, they have the capacity to influence plant growth, developmental, and fitness processes, and in addition, the industrial end-use of starch. However, little is known about the role of SBEs in determining starch structure-function relations in economically important horticultural crops such as fruit and leafy greens, many of which accumulate starch transiently. Further, a full understanding of the biological function of these types of starches is lacking. Because of this gap in knowledge, this minireview aims to provide an overview of SBEs in horticultural crops, to investigate the potential role of starch in determining postharvest quality. A systematic examination of SBE sequences in 43 diverse horticultural species, identified SBE1, 2 and 3 isoforms in all species examined except apple, olive, and Brassicaceae, which lacked SBE1, but had a duplicated SBE2. Among our findings after a comprehensive and critical review of published data, was that as apple, banana, and tomato fruits ripens, the ratio of the highly digestible amylopectin component of starch increases relative to the more digestion-resistant amylose fraction, with parallel increases in SBE2 transcription, fruit sugar content, and decreases in starch. It is tempting to speculate that during the ripening of these fruit when starch degradation occurs, there are rearrangements made to the structure of starch possibly via branching enzymes to increase starch digestibility to sugars. We propose that based on the known action of SBEs, and these observations, SBEs may affect produce quality, and shelf-life directly through starch accumulation, and indirectly, by altering sugar availability. Further studies where SBE activity is fine-tuned in these crops, can enrich our understanding of the role of starch across species and may improve horticulture postharvest quality.

An apple (Malus domestica) AP2/ERF transcription factor modulates carotenoid accumulation

Color is an important trait for horticultural crops. Carotenoids are one of the main pigments for coloration and have important implications for photosynthesis in plants and benefits for human health. Here, we identified an APETALA2 (AP2)/ETHYLENE RESPONSE FACTOR (ERF) transcription factor named MdAP2-34 in apple (Malus domestica Borkh.). MdAP2-34 expression exhibited a close correlation with carotenoid content in 'Benin Shogun' and 'Yanfu 3' fruit flesh. MdAP2-34 promotes carotenoid accumulation in MdAP2-34-OVX transgenic apple calli and fruits by participating in the carotenoid biosynthesis pathway. The major carotenoid contents of phytoene and β-carotene were much higher in overexpressing MdAP2-34 transgenic calli and fruit skin, yet the predominant compound of lutein showed no obvious difference, indicating that MdAP2-34 regulates phytoene and β-carotene accumulation but not lutein. MdPSY2-1 (phytoene synthase 2) is a major gene in the carotenoid biosynthesis pathway in apple fruit, and the MdPSY2-1 gene is directly bound and transcriptionally activated by MdAP2-34. In addition, overexpressing MdPSY2-1 in apple calli mainly increases phytoene and total carotenoid contents. Our findings will advance and extend our understanding of the complex molecular mechanisms of carotenoid biosynthesis in apple, and this research is valuable for accelerating the apple breeding process.

Insights into angiosperm evolution, floral development and chemical biosynthesis from the Aristolochia fimbriata genome

Aristolochia, a genus in the magnoliid order Piperales, has been famous for centuries for its highly specialized flowers and wide medicinal applications. Here, we present a new, high-quality genome sequence of Aristolochia fimbriata, a species that, similar to Amborella trichopoda, lacks further whole-genome duplications since the origin of extant angiosperms. As such, the A. fimbriata genome is an excellent reference for inferences of angiosperm genome evolution, enabling detection of two novel whole-genome duplications in Piperales and dating of previously reported whole-genome duplications in other magnoliids. Genomic comparisons between A. fimbriata and other angiosperms facilitated the identification of ancient genomic rearrangements suggesting the placement of magnoliids as sister to monocots, whereas phylogenetic inferences based on sequence data we compiled yielded ambiguous relationships. By identifying associated homologues and investigating their evolutionary histories and expression patterns, we revealed highly conserved floral developmental genes and their distinct downstream regulatory network that may contribute to the complex flower morphology in A. fimbriata. Finally, we elucidated the genetic basis underlying the biosynthesis of terpenoids and aristolochic acids in A. fimbriata.

Genome-wide identification, characterization and expression analysis of MATE family genes in apple (Malus × domestica Borkh)

BACKGROUND

As an important group of the multidrug efflux transporter family, the multidrug and toxic compound extrusion (MATE) family has a wide range of functions and is distributed in all kingdoms of living organisms. However, only two MATE genes in apple have been analyzed and genome-wide comprehensive analysis of MATE family is needed.

RESULTS

In this study, a total of 66 MATE (MdMATE) candidates encoding putative MATE transporters were identified in the apple genome. These MdMATE genes were classified into four groups by phylogenetic analysis with MATE genes in Arabidopsis. Synteny analysis reveals that whole genome duplication (WGD) and segmental duplication events played a major role in the expansion of MATE gene family in apple. MdMATE genes show diverse expression patterns in different tissues/organs and developmental stages. Analysis of cis-regulatory elements in MdMATE promoter regions indicates that the function of MdMATE genes is mainly related to stress response. Besides, the changes of gene expression levels upon different pathogen infections reveal that MdMATE genes are involved in biotic stress response.

CONCLUSIONS

In this work, we systematically identified MdMATE genes in apple genome using a set of bioinformatics approaches. Our comprehensive analysis provided valuable resources for improving disease resistance in apple and further functional characterization of MATE genes in other species.

LjaFGD: Lonicera japonica functional genomics database

Lonicera japonica Thunb., a traditional Chinese herb, has been used for treating human diseases for thousands of years. Recently, the genome of L. japonica has been decoded, providing valuable information for research into gene function. However, no comprehensive database for gene functional analysis and mining is available for L. japonica. We therefore constructed LjaFGD (www.gzybioinformatics.cn/LjaFGD and bioinformatics.cau.edu.cn/LjaFGD), a database for analyzing and comparing gene function in L. japonica. We constructed a gene co-expression network based on 77 RNA-seq samples, and then annotated genes of L. japonica by alignment against protein sequences from public databases. We also introduced several tools for gene functional analysis, including Blast, motif analysis, gene set enrichment analysis, heatmap analysis, and JBrowse. Our co-expression network revealed that MYB and WRKY transcription factor family genes were co-expressed with genes encoding key enzymes in the biosynthesis of chlorogenic acid and luteolin in L. japonica. We used flavonol synthase 1 (LjFLS1) as an example to show the reliability and applicability of our database. LjaFGD and its various associated tools will provide researchers with an accessible platform for retrieving functional information on L. japonica genes to further biological discovery.

Genome-wide analysis of SET-domain group histone methyltransferases in apple reveals their role in development and stress responses

Background

Histone lysine methylation plays an important role in plant development and stress responses by activating or repressing gene expression. Histone lysine methylation is catalyzed by a class of SET-domain group proteins (SDGs). Although an increasing number of studies have shown that SDGs play important regulatory roles in development and stress responses, the functions of SDGs in apple remain unclear.

Results

A total of 67 SDG members were identified in the Malus×domestica genome. Syntenic analysis revealed that most of the MdSDG duplicated gene pairs were associated with a recent genome-wide duplication event of the apple genome. These 67 MdSDG members were grouped into six classes based on sequence similarity and the findings of previous studies. The domain organization of each MdSDG class was characterized by specific patterns, which was consistent with the classification results. The tissue-specific expression patterns of MdSDGs among the 72 apple tissues in the different apple developmental stages were characterized to provide insight into their potential functions in development. The expression profiles of MdSDGs were also investigated in fruit development, the breaking of bud dormancy, and responses to abiotic and biotic stress; the results indicated that MdSDGs might play a regulatory role in development and stress responses. The subcellular localization and putative interaction network of MdSDG proteins were also analyzed.

Conclusions

This work presents a fundamental comprehensive analysis of SDG histone methyltransferases in apple and provides a basis for future studies of MdSDGs involved in apple development and stress responses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07596-0.

An Atlas of Genomic Resources for Studying Rosaceae Fruits and Ornamentals

Rosaceae, a large plant family of more than 3,000 species, consists of many economically important fruit and ornamental crops, including peach, apple, strawberry, raspberry, cherry, and rose. These horticultural crops are not only important economic drivers in many regions of the world, but also major sources of human nutrition. Additionally, due to the diversity of fruit types in Rosaceae, this plant family offers excellent opportunities for investigations into fleshy fruit diversity, evolution, and development. With the development of high-throughput sequencing technologies and computational tools, an increasing number of high-quality genomes and transcriptomes of Rosaceae species have become available and will greatly facilitate Rosaceae research and breeding. This review summarizes major genomic resources and genome research progress in Rosaceae, highlights important databases, and suggests areas for further improvement. The availability of these big data resources will greatly accelerate research progress and enhance the agricultural productivity of Rosaceae.

Intricate genetic variation networks control the adventitious root growth angle in apple

Background

The root growth angle (RGA) typically determines plant rooting depth, which is significant for plant anchorage and abiotic stress tolerance. Several quantitative trait loci (QTLs) for RGA have been identified in crops. However, the underlying mechanisms of the RGA remain poorly understood, especially in apple rootstocks. The objective of this study was to identify QTLs, validate genetic variation networks, and develop molecular markers for the RGA in apple rootstock.

Results

Bulked segregant analysis by sequencing (BSA-seq) identified 25 QTLs for RGA using 1955 hybrids of the apple rootstock cultivars ‘Baleng Crab’ (Malus robusta Rehd., large RGA) and ‘M9’ (M. pumila Mill., small RGA). With RNA sequencing (RNA-seq) and parental resequencing, six major functional genes were identified and constituted two genetic variation networks for the RGA. Two single nucleotide polymorphisms (SNPs) of the MdLAZY1 promoter damaged the binding sites of MdDREB2A and MdHSFB3, while one SNP of MdDREB2A and MdIAA1 affected the interactions of MdDREB2A/MdHSFB3 and MdIAA1/MdLAZY1, respectively. A SNP within the MdNPR5 promoter damaged the interaction between MdNPR5 and MdLBD41, while one SNP of MdLBD41 interrupted the MdLBD41/MdbHLH48 interaction that affected the binding ability of MdLBD41 on the MdNPR5 promoter. Twenty six SNP markers were designed on candidate genes in each QTL interval, and the marker effects varied from 0.22°-26.11°.

Conclusions

Six diagnostic markers, SNP592, G122, b13, Z312, S1272, and S1288, were used to identify two intricate genetic variation networks that control the RGA and may provide new insights into the accuracy of the molecular markers. The QTLs and SNP markers can potentially be used to select deep-rooted apple rootstocks.

Genomics-assisted prediction of salt and alkali tolerances and functional marker development in apple rootstocks

Background

Saline, alkaline, and saline-alkaline stress severely affect plant growth and development. The tolerance of plants to these stressors has long been important breeding objectives, especially for woody perennials like apple. The aims of this study were to identify quantitative trait loci (QTLs) and to develop genomics-assisted prediction models for salt, alkali, and salt-alkali tolerance in apple rootstock.

Results

A total of 3258 hybrids derived from the apple rootstock cultivars ‘Baleng Crab’ (Malus robusta Rehd., tolerant) × ‘M9’ (M. pumila Mill., sensitive) were used to identify 17, 13, and two QTLs for injury indices of salt, alkali, and salt–alkali stress via bulked segregant analysis. The genotype effects of single nucleotide polymorphism (SNP) markers designed on candidate genes in each QTL interval were estimated. The genomic predicted value of an individual hybrid was calculated by adding the sum of all marker genotype effects to the mean phenotype value of the population. The prediction accuracy was 0.6569, 0.6695, and 0.5834 for injury indices of salt, alkali, and salt–alkali stress, respectively. SNP182G on MdRGLG3, which changes a leucine to an arginine at the vWFA-domain, conferred tolerance to salt, alkali, and salt-alkali stress. SNP761A on MdKCAB, affecting the Kv_beta domain that cooperated with the linked allelic variation SNP11, contributed to salt, alkali, and salt–alkali tolerance in apple rootstock.

Conclusions

The genomics-assisted prediction models can potentially be used in breeding saline, alkaline, and saline-alkaline tolerant apple rootstocks. The QTLs and the functional markers may provide insight for future studies into the genetic variation of plant abiotic stress tolerance.

Coexpression Analysis Reveals Dynamic Modules Regulating the Growth and Development of Cirri in the Rattans (Calamus simplicifolius and Daemonorops jenkinsiana)

Rattan is regarded as one of the major non-timber forest products, second only to wood and bamboo, worldwide. Although the published genomes of Calamus simplicifolius and Daemonorops jenkinsiana have facilitated genome-wide gene functional analyses, coexpression networks (CENs) provide more comprehensive and complete annotations of gene function at the transcriptome level. Thus, we analyzed the CENs of the two rattans, C. simplicifolius and D. jenkinsiana, by integrating the genome sequences and analyzing in-house transcriptome data from different development stages of their cirri using a well-developed strategy. A total of 3,504 and 3,027 functional modules were identified in C. simplicifolius and D. jenkinsiana, respectively, based on a combination of CENs, gene family classification, and function enrichment tools. These modules covered the major developmental processes, including photosynthesis, lignin biosynthesis, flavonoid biosynthesis, and phenylpropanoid biosynthesis. Reference annotations were refined using CENs and functional modules. Moreover, we obtained novel insights into the regulation of cirrus growth and development in rattans. Furthermore, Rattan-NET (http://rattan.bamboogdb.org/), an online database with analysis tools for gene set enrichment analysis, module enrichment, network comparison analysis, and cis-element analysis, was constructed for the easy analysis of gene function and regulation modules involved in the growth and development of cirri in rattans.