Genome-Wide Identification of Sweet Orange WRKY Transcription Factors and Analysis of Their Expression in Response to Infection by Penicillium digitatum

The genomes of Australian wild limes

Australian wild limes occur in highly diverse range of environments and are a unique genetic resource within the genus Citrus. Here we compare the haplotype-resolved genome assemblies of six Australian native limes, including four new assemblies generated using PacBio HiFi and Hi-C sequencing data. The size of the genomes was between 315 and 391 Mb with contig N50s from 29.5 to 35 Mb. Gene completeness of the assemblies was estimated to be from 98.4 to 99.3% and the annotations from 97.7 to 98.9% based upon BUSCO, confirming the high contiguity and completeness of the assembled genomes. High collinearity was observed among the genomes and the two haplotype assemblies for each species. Gene duplication and evolutionary analysis demonstrated that the Australian citrus have undergone only one ancient whole-genome triplication event during evolution. The highest number of species-specific and expanded gene families were found in C. glauca and they were primarily enriched in purine, thiamine metabolism, amino acids and aromatic amino acids metabolism which might help C. glauca to mitigate drought, salinity, and pathogen attacks in the drier environments in which this species is found. Unique genes related to terpene biosynthesis, glutathione metabolism, and toll-like receptors in C. australasica, and starch and sucrose metabolism genes in both C. australis and C. australasica might be important candidate genes for HLB tolerance in these species. Expanded gene families were not lineage specific, however, a greater number of genes related to plant-pathogen interactions, predominantly disease resistant protein, was found in C. australasica and C. australis.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Identification of the sweet orange (Citrus sinensis) bHLH gene family and the role of CsbHLH55 and CsbHLH87 in regulating salt stress

Salt stress is one of the primary environmental stresses limiting plant growth and production and adversely affecting the growth, development, yield, and fruit quality of Citrus sinensis. bHLH (basic helix-loop-helix) genes are involved in many bioregulatory processes in plants, including growth and development, phytohormone signaling, defense responses, and biosynthesis of specific metabolites. In this study, by bioinformatics methods, 120 CsbHLHgenes were identified, and phylogenetic analysis classified them into 18 subfamilies that were unevenly distributed on nine chromosomes. The cis-acting elements of the CsbHLH genes were mainly hormone-related cis-acting elements. Seventeen CsbHLH genes exhibited significant differences in expression under salt stress. Six CsbHLH genes with significant differences in expression were randomly selected for quantitative real-time polymerase chain reaction (qRT-PCR) validation. The qRT-PCR results showed a strong correlation with the transcriptome data. Phytohormones such as jasmonic acid (JA) are essential for biotic and abiotic stress responses in plants, and CsbHLH55 and CsbHLH87 are considered candidate target genes for sweet orange MYC2 transcription factors involved in the JA signaling pathway. These genes are the main downstream effectors in the JA signaling pathway and can be activated to participate in the JA signaling pathway. Activation of the JA signaling pathway inhibits the production of reactive oxygen species and improves the salt tolerance of sweet orange plants. The CsbHLH55 and CsbHLH87 genes could be candidate genes for breeding new transgenic salt-resistant varieties of sweet orange.

Genome-wide identification of bZIP transcription factors in 12 Rosaceae species and modeling of novel mechanisms of EjbZIPs response to salt stress

In plantae, basic leucine zipper (bZIP) transcription factors (TFs) are widespread and regulate a variety of biological processes under abiotic stress. However, it has not been extensively studied in Rosaceae, and the functional effects of bZIP on Eriobotrya japonica under salt stress are still unknown. Therefore, in this study, the bZIP TF family of 12 species of Rosaceae was analyzed by bioinformatics method, and the expression profile and quantitative real-time polymerase chain reaction of E. japonica under salt stress were analyzed. The results showed that a total of 869 bZIP TFs were identified in 12 species of Rosaceae and divided into nine subfamilies. Differences in promoter cis-elements between subfamilies vary depending on their role. Species belonging to the same subfamily have a similar number of chromosomes and the number of genes contained on each chromosome. Gene duplication analysis has found segmental duplication to be a prime force in the evolution of Rosaceae species. In addition, nine EjbZIPs were significantly different, including seven up-regulated and two down-regulated in E. japonica under salt stress. Especially, EjbZIP13 was involved in the expression of SA-responsive proteins by binding to the NPR1 gene. EjbZIP27, EjbZIP30, and EjbZIP38 were highly expressed in E. japonica under salt stress, thus improving the salt tolerance capacity of the plants. These results can provide a theoretical basis for exploring the characteristics and functions of the bZIP TF family in more species and breeding salt-tolerant E. japonica varieties. It also provides a reference for resolving the response mechanism of bZIP TF in 12 Rosaceae species under salt stress.

Genome-wide identification and expression analysis of the MADS gene family in sweet orange (Citrus sinensis) infested with pathogenic bacteria

The risk of pathogenic bacterial invasion in plantations has increased dramatically due to high environmental climate change and has seriously affected sweet orange fruit quality. MADS genes allow plants to develop increased resistance, but functional genes for resistance associated with pathogen invasion have rarely been reported. MADS gene expression profiles were analyzed in sweet orange leaves and fruits infested with Lecanicillium psalliotae and Penicillium digitatum, respectively. Eighty-two MADS genes were identified from the sweet orange genome, and they were classified into five prime subfamilies concerning the Arabidopsis MADS gene family, of which the MIKC subfamily could be subdivided into 13 minor subfamilies. Protein structure analysis showed that more than 93% of the MADS protein sequences of the same subfamily between sweet orange and Arabidopsis were very similar in tertiary structure, with only CsMADS8 and AG showing significant differences. The variability of MADS genes protein structures between sweet orange and Arabidopsis subgroups was less than the variabilities of protein structures within species. Chromosomal localization and covariance analysis showed that these genes were unevenly distributed on nine chromosomes, with the most genes on chromosome 9 and the least on chromosome 2, with 36 and two, respectively. Four pairs of tandem and 28 fragmented duplicated genes in the 82 MADS gene sequences were found in sweet oranges. GO (Gene Ontology) functional enrichment and expression pattern analysis showed that the functional gene CsMADS46 was strongly downregulated of sweet orange in response to biotic stress adversity. It is also the first report that plants' MADS genes are involved in the biotic stress responses of sweet oranges. For the first time, L. psalliotae was experimentally confirmed to be the causal agent of sweet orange leaf spot disease, which provides a reference for the research and control of pathogenic L. psalliotae.

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.