Genome-wide analyses of genes encoding FK506-binding proteins reveal their involvement in abiotic stress responses in apple

Pan-genome survey of Septoria pistaciarum, causal agent of Septoria leaf spot of pistachios, across three Aegean sub-regions of Greece

Septoria pistaciarum, a causal agent of Septoria leaf spot disease of pistachio, is a fungal pathogen that causes substantial losses in the cultivation, worldwide. This study describes the first pan-genome-based survey of this phytopathogen-comprising a total of 27 isolates, with 9 isolates each from 3 regional units of Greece (Pieria, Larissa and Fthiotida). The reference isolate (SPF8) assembled into a total of 43.1 Mb, with 38.6% contained within AT-rich regions of approximately 37.5% G:C. The genomes of the 27 isolates exhibited on average 42% gene-coding and 20% repetitive regions. The genomes of isolates from the southern Fthiotida region appeared to more diverged from each other than the other regions based on SNP-derived trees, and also contained isolates similar to both the Pieria and Larissa regions. In contrast, isolates of the Pieria and Larissa were less diverse and distinct from one another. Asexual reproduction appeared to be typical, with no MAT1-2 locus detected in any isolate. Genome-based prediction of infection mode indicated hemibiotrophic and saprotrophic adaptations, consistent with its long latent phase. Gene prediction and orthology clustering generated a pan-genome-wide gene set of 21,174 loci. A total of 59 ortholog groups were predicted to contain candidate effector proteins, with 36 (61%) of these either having homologs to known effectors from other species or could be assigned predicted functions from matches to conserved domains. Overall, effector prediction suggests that S. pistaciarum employs a combination of defensive effectors with roles in suppression of host defenses, and offensive effectors with a range of cytotoxic activities. Some effector-like ortholog groups presented as divergent versions of the same protein, suggesting region-specific adaptations may have occurred. These findings provide insights and future research directions in uncovering the pathogenesis and population dynamics of S. pistaciarum toward the efficient management of Septoria leaf spot of pistachio.

High Arctic seawater and coastal soil microbiome co-occurrence and composition structure and their potential hydrocarbon biodegradation

The accelerated decline in Arctic sea-ice cover and duration is enabling the opening of Arctic marine passages and improving access to natural resources. The increasing accessibility to navigation and resource exploration and production brings risks of accidental hydrocarbon releases into Arctic waters, posing a major threat to Arctic marine ecosystems where oil may persist for many years, especially in beach sediment. The composition and response of the microbial community to oil contamination on Arctic beaches remain poorly understood. To address this, we analyzed microbial community structure and identified hydrocarbon degradation genes among the Northwest Passage intertidal beach sediments and shoreline seawater from five high Arctic beaches. Our results from 16S/18S rRNA genes, long-read metagenomes, and metagenome-assembled genomes reveal the composition and metabolic capabilities of the hydrocarbon microbial degrader community, as well as tight cross-habitat and cross-kingdom interactions dominated by lineages that are common and often dominant in the polar coastal habitat, but distinct from petroleum hydrocarbon-contaminated sites. In the polar beach sediment habitats, Granulosicoccus sp. and Cyclocasticus sp. were major potential hydrocarbon-degraders, and our metagenomes revealed a small proportion of microalgae and algal viruses possessing key hydrocarbon biodegradative genes. This research demonstrates that Arctic beach sediment and marine microbial communities possess the ability for hydrocarbon natural attenuation. The findings provide new insights into the viral and microalgal communities possessing hydrocarbon degradation genes and might represent an important contribution to the removal of hydrocarbons under harsh environmental conditions in a pristine, cold, and oil-free environment that is threatened by oil spills.

Genome-wide association study and pathway analysis to decipher loci associated with Fusarium ear rot resistance in tropical maize germplasm

Breeding for host resistance is the most efficient and environmentally safe method to curb the spread of fusarium ear rot (FER). However, conventional breeding for resistance to FER is hampered by the complex polygenic nature of this trait, which is highly influenced by environmental conditions. This study aimed to identify genomic regions, single nucleotide polymorphisms (SNPs), and putative candidate genes associated with FER resistance as well as candidate metabolic pathways and pathway genes involved in it. A panel of 151 tropical inbred maize lines were used to assess the genetic architecture of FER resistance over two seasons. During the study period, seven SNPs associated with FER resistance were identified on chromosomes 1, 2, 4, 5, and 9, accounting for 4-11% of the phenotypic variance. These significant markers were annotated into four genes. Seven significant metabolic pathways involved in FER resistance were identified using the Pathway Association Study Tool, the most significant being the superpathway of the glyoxylate cycle. Overall, this study confirmed that resistance to FER is indeed a complex mechanism controlled by several small to medium-effect loci. Our findings may contribute to fast-tracking the efforts to develop disease-resistant maize lines through marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10722-023-01793-4.

Genome-Wide Identification and Expression Analysis of the Cucumber FKBP Gene Family in Response to Abiotic and Biotic Stresses

The FKBP (FK506-binding protein) gene family is an important member of the PPlase protease family and plays a vital role during the processes of plant growth and development. However, no studies of the FKBP gene family have been reported in cucumber. In this study, 19 FKBP genes were identified in cucumber, which were located on chromosomes 1, 3, 4, 6, and 7. Phylogenetic analysis divided the cucumber FKBP genes into three subgroups. The FKBP genes in the same subgroup exhibited similar structures and conserved motifs. The cis-acting elements analysis revealed that the promoters of cucumber FKBP genes contained hormone-, stress-, and development-related cis-acting elements. Synteny analysis of the FKBP genes among cucumber, Arabidopsis, and rice showed that 12 kinds of syntenic relationships were detected between cucumber and Arabidopsis FKBP genes, and 3 kinds of syntenic relationships were observed between cucumber and rice FKBP genes. The tissue-specific expression analysis showed that some FKBP genes were expressed in all tissues, while others were only highly expressed in part of the 10 types of tissues. The expression profile analysis of cucumber FKBP genes under 13 types of stresses showed that the CsaV3_1G007080 gene was differentially expressed under abiotic stresses (high temperature, NaCl, silicon, and photoperiod) and biotic stresses (downy mildew, green mottle mosaic virus, Fusarium wilt, phytophthora capsica, angular leaf spot, and root-knot nematode), which indicated that the CsaV3_1G007080 gene plays an important role in the growth and development of cucumber. The interaction protein analysis showed that most of the proteins in the FKBP gene family interacted with each other. The results of this study will lay the foundation for further research on the molecular biological functions of the cucumber FKBP gene family.

Genome-wide characterization of FK506-binding proteins, parvulins and phospho-tyrosyl phosphatase activators in wheat and their regulation by heat stress

Peptidyl-prolyl cis-trans isomerases (PPIases) are ubiquitous proteins which are essential for cis-trans isomerisation of peptide bonds preceding the proline residue. PPIases are categorized into four sub-families viz., cyclophilins, FK506-binding proteins (FKBPs), parvulins and protein phosphatase 2A phosphatase activators (PTPAs). Apart from catalysing the cis-trans isomerization, these proteins have also been implicated in diverse cellular functions. Though PPIases have been identified in several important crop plants, information on these proteins, except cyclophilins, is scanty in wheat. In order to understand the role of these genes in wheat, we carried out genome-wide identification using computational approaches. The present study resulted in identification of 71 FKBP (TaFKBP) 12 parvulin (TaPar) and 3 PTPA (TaPTPA) genes in hexaploid wheat genome, which are distributed on different chromosomes with uneven gene densities. The TaFKBP and TaPar proteins, besides PPIase domain, also contain additional domains, indicating functional diversification. In silico prediction also revealed that TaFKBPs are localized to ER, nucleus, chloroplast and cytoplasm, while the TaPars are confined to cytoplasm and nucleus. The TaPTPAs, on the contrary, appear to be present only in the cytoplasm. Evolutionary studies predicted that most of the TaFKBP, TaPar and TaPTPA genes in hexaploid wheat have been derived from their progenitor species, with some events of loss or gain. Syntenic analysis revealed the presence of many collinear blocks of TaFKBP genes in wheat and its sub-genome donors. qRT-PCR analysis demonstrated that expression of TaFKBP and TaPar genes is regulated differentially by heat stress, suggesting their likely involvement in thermotolerance. The findings of this study will provide basis for further functional characterization of these genes and their likely applications in crop improvement.

Genome-Wide Identification and Analysis of FKBP Gene Family in Wheat (Triticum asetivum)

FK506-binding protein (FKBP) genes have been found to play vital roles in plant development and abiotic stress responses. However, limited information is available about this gene family in wheat (Triticum aestivum L.). In this study, a total of 64 FKBP genes were identified in wheat via a genome-wide analysis involving a homologous search of the latest wheat genome data, which was unevenly distributed in 21 chromosomes, encoded 152 to 649 amino acids with molecular weights ranging from 16 kDa to 72 kDa, and was localized in the chloroplast, cytoplasm, nucleus, mitochondria, peroxisome and endoplasmic reticulum. Based on sequence alignment and phylogenetic analysis, 64 TaFKBPs were divided into four different groups or subfamilies, providing evidence of an evolutionary relationship with Aegilops tauschii, Brachypodium distachyon, Triticum dicoccoides, Arabidopsis thaliana and Oryza sativa. Hormone-related, abiotic stress-related and development-related cis-elements were preferentially presented in promoters of TaFKBPs. The expression levels of TaFKBP genes were investigated using transcriptome data from the WheatExp database, which exhibited tissue-specific expression patterns. Moreover, TaFKBPs responded to drought and heat stress, and nine of them were randomly selected for validation by qRT-PCR. Yeast cells expressing TaFKBP19-2B-2 or TaFKBP18-6B showed increased influence on drought stress, indicating their negative roles in drought tolerance. Collectively, our results provide valuable information about the FKBP gene family in wheat and contribute to further characterization of FKBPs during plant development and abiotic stress responses, especially in drought stress.

MicroRNAs expression profiles in early responses to different levels of water deficit in Setaria viridis

Water deficit is a major constraint for crops of economic importance in almost all agricultural regions. However, plants have an active defense system to adapt to these adverse conditions, acting in the reprogramming of gene expression responsible for encoding microRNAs (miRNAs). These miRNAs promote the regulation to the target gene expression by the post-transcriptional (PTGS) and transcriptional gene silencing (TGS), modulating several pathways including defense response to water deficit. The broader knowledge of the miRNA expression profile and its regulatory networks in response to water deficit can provide evidence for the development of new biotechnological tools for genetic improvement of several important crops. In this study, we used Setaria viridis accession A10.1 as a C4 model plant to widely investigate the miRNA expression profile in early responses to different levels of water deficit. Ecophysiological studies in Setaria viridis under water deficit and after rewatering demonstrated a drought tolerant accession, capable of a rapid recovery from the stress. Deep small RNA sequencing and degradome studies were performed in plants submitted to drought to identify differentially expressed miRNA genes and their predicted targets, using in silico analysis. Our findings showed that several miRNAs were differentially modulated in response to distinctive levels of water deficit and after rewatering. The predicted mRNA targets mainly corresponded to genes related to cell wall remodeling, antioxidant system and drought-related transcription factors, indicating that these genes are rapidly regulated in early responses to drought stress. The implications of these modulations are extensively discussed, and higher-effect miRNAs are suggested as major players for potential use in genetic engineering to improve drought tolerance in economically important crops, such as sugarcane, maize, and sorghum.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01226-z.

MdWRKY30, a group IIa WRKY gene from apple, confers tolerance to salinity and osmotic stresses in transgenic apple callus and Arabidopsis seedlings

Abiotic stresses threaten the productivity and quality of economically important perennial fruit crops such as apple (Malus × domestica Borkh.). WRKY transcription factors play various roles in plant responses to abiotic stress, but little is known regarding WRKY genes in apple. Here, we carried out functional characterization of an apple Group IIa WRKY gene (MdWRKY30). qRT-PCR analysis found that MdWRKY30 expression was induced by salt and drought stress. A subcellular localization assay showed that MdWRKY30 is localized to the nucleus. A transactivation assay found that MdWRKY30 has no transcriptional activation activity. A Y2H assay indicated that MdWRKY26, MdWRKY28, and MdWRKY30 interact with each other to form heterodimers and homodimers. Transgenic analysis revealed that the overexpression of MdWRKY30 in Arabidopsis enhanced salt and osmotic tolerance in the seedling stage, as well as during the seed germination and greening cotyledon stages. MdWRKY30 overexpression enhanced tolerance to salt and osmotic stresses in transgenic apple callus through transcriptional regulation of stress-related genes. Together, our results demonstrate that MdWRKY30 is an important regulator of salinity and osmotic stress tolerance in apple.

Cloning, sequencing, and expression analysis of 32 NAC transcription factors (MdNAC) in apple

BACKGROUND

NAC transcription factors play important roles in the regulation of plant growth, development, abiotic and biotic stress responses. The transcriptional level of MdNACs in different tissues and under various biotic and abiotic stress treatments was determined to provide a solid foundation for studying the function of MdNACs.

METHODS

Thirty-two full-length cDNA sequences of Md NACs were isolated by homologous comparison and RT-PCR confirmation, and the obtained cDNA sequences and the deduced amino acid sequences were analyzed with bioinformatics methods. The prediction of subcellular locations of MdNAC proteins was performed using CELLO v.2.5, PSORT, and SoftBerry ProtComp 9.0. Expression levels of MdNACs were detected in 16 different tissues using an array. Expression patterns of MdNACs were detected in response to Alternaria alternata apple pathotype (AAAP) infection using RNA-seq, and the expression of MdNACs was analyzed under NaCl and mannitol treatments using RT-qPCR.

RESULTS

The sequencing results produced 32 cDNAs (designated as MdNAC24-39, MdNAC54-65, and MdNAC67-70 with GenBank accession No. MG099861-MG099876, MG099891-MG099902, and MG099904-MG099907, respectively). Phylogenetic analysis revealed that MdNAC34 belonged to the ATAF group, MdNAC63 belonged to the AtNAC3 group, MdNAC24, MdNAC26-30, MdNAC32-33, MdNAC35, MdNAC37-39, MdNAC56-57, MdNAC59-62, MdNAC64-65, and MdNAC67-70 belonged to the NAM group, and MdNAC25, MdNAC36, MdNAC54-55, and MdNAC58 belonged to the VND group. Predictions of subcellular localization showed that MdNAC24-27, MdNAC29-30, MdNAC33-37, MdNAC39, MdNAC54-65, and MdNAC67-70 proteins were located in the nucleus, MdNAC28 proteins were located in the cytoplasm, MdNAC31-32 proteins were located in the nucleus and cytoplasm, and MdNAC38 proteins were located in the nucleus and plasma membrane. Array results indicated that 32 MdNACs were expressed in all examined tissues at various expression levels. RNA-seq results showed that expression levels of MdNAC26-28, MdNAC33-34, MdNAC60, MdNAC62-65, and MdNAC68 were induced, but MdNAC24, MdNAC32, and MdNAC58 were down-regulated in response to AAAP infection. Under salt treatment, MdNAC24, MdNAC27, MdNAC29, MdNAC34, MdNAC37, MdNAC39, MdNAC54, MdNAC59, and MdNAC63 transcription levels were induced. Under mannitol treatment, MdNAC32 and MdNAC54 transcription levels were induced, but MdNAC24, MdNAC28, MdNAC30, MdNAC33, MdNAC35, MdNAC37, MdNAC55, MdNAC56, MdNAC58, and MdNAC59 were down-regulated. Taken together, the results indicated that the cloned MdNAC genes were expressed constitutively in all examined tissues. These genes were up-regulated or down-regulated in response to AAAP infection and to salt or mannitol, which suggested they may be involved in the regulation of growth, development, and stress response in apple.

Isolation, sequencing, and expression analysis of 30 AP2/ERF transcription factors in apple

BACKGROUND

AP2/ERF transcription factors are involved in the regulation of plant growth, development, and stress responses. Our research objective was to characterize novel apple (Malus × domestica Borkh.) genes encoding AP2/ERF transcription factors involved in regulation of plant growth, development, and stress response. The transcriptional level of apple AP2/ERF genes in different tissues and under various biotic and abiotic stress was determined to provide valuable insights into the function of AP2/ERF transcription factors in apple.

METHODS

Thirty full-length cDNA sequences of apple AP2/ERF genes were isolated from 'Zihong Fuji' apple (Malus × domestica cv. Zihong Fuji) via homologous comparison and RT-PCR confirmation, and the obtained cDNA sequences and the deduced amino acid sequences were analyzed with bioinformatics methods. Expression levels of apple AP2/ERF genes were detected in 16 different tissues using a known array. Expression patterns of apple AP2/ERF genes were detected in response to Alternaria alternata apple pathotype (AAAP) infection using RNA-seq with existing data, and the expression of apple AP2/ERF genes was analyzed under NaCl and mannitol treatments using qRT-PCR.

RESULTS

The sequencing results produced 30 cDNAs (designated as MdERF3-8, MdERF11, MdERF16-19, MdERF22-28, MdERF31-35, MdERF39, MdAP2D60, MdAP2D62-65, and MdRAV2). Phylogenetic analysis revealed that MdERF11/16, MdERF33/35, MdERF34/39, and MdERF18/23 belonged to groups A-2, A-4, A-5, and A-6 of the DREB subfamily, respectively; MdERF31, MdERF19, MdERF4/25/28/32, MdERF24, MdERF5/6/27, and MdERF3/7/8/17/22/26 belonged to groups B-1, B-2, B-3, B-4, B-5, and B-6 of the ERF subfamily, respectively; MdAP2D60 and MdAP2D62/63/64/65 belonged to the AP2 subfamily; and MdRAV2 belonged to the RAV subfamily. Array results indicated that 30 apple AP2/ERF genes were expressed in all examined tissues to different degrees. RNA-seq results using previously reported data showed that many members of the apple ERF and DREB subfamilies were induced by Alternaria alternate apple pathotype (AAAP) infection. Under salt treatment, many members in the apple ERF and DREB subfamilies were transcriptionally up or down-regulated. Under mannitol treatment, many members of the apple ERF, DREB, and AP2 subfamilies were induced at the transcriptional level. Taken together, the results indicated that the cloned apple AP2/ERF genes were expressed in all examined tissues. These genes were up-regulated or down-regulated in response to AAAP infection and to salt or mannitol treatment, which suggested they may be involved in regulating growth, development, and stress response in apple.

Genome-wide identification, expression, and interaction analysis for ovate family proteins in peach

Ovate family proteins (OFPs), which are involved in aspects of plant development and growth, is a class of plant-specific transcription factors. Although OFPs have been reported in some species, little is known about their evolution, structure, fruit developmental expression, and interactions among OFP members in peach (Prunus persica). In this study, 15 peach OFP (PpOFP) genes were identified. Phylogenetic analysis showed that 716 OFPs from 32 species were divided into 15 subgroups; 10 subgroups (Ia, Ib, Ic, Id, Ie, If, Ig, Ih, Ii, and Ij) were composed of dicotyledonous plants only and five (IIa, IIb, IIc, IId, and IIe) were composed of monocotyledonous plants only. Structure analysis showed that the OFP genes in monocotyledonous and dicotyledonous plants had no introns. Chromosomal localization analysis showed that 15 PpOFP genes were unevenly mapped on seven chromosomes. Furthermore, eight of the 15 PpOFP genes were cloned successfully by the RT-PCR method. To some extent, eight PpOFPs were expressed in all the detected peach tissues. In addition, analysis by Y2H and BiFC technologies indicated that both PpOFP4 and PpOFP5 formed homodimers with themselves, and PpOFP5 interacted with PpOFP7 or PpOFP8 to form heterodimers. These results serve as the theoretical basis for the analysis of the biological function and regulation of peach OFP transcription factors in the growth, development and other conditions, as well as evolution studies of OFP transcription factors in higher plants.

TOR inhibitors: from mammalian outcomes to pharmacogenetics in plants and algae

Target of rapamycin (TOR) is a conserved eukaryotic phosphatidylinositol 3-kinase-related kinase that regulates growth and metabolism in response to environment in plants and algae. The study of the plant and algal TOR pathway has largely depended on TOR inhibitors first developed for non-photosynthetic eukaryotes. In animals and yeast, fundamental work on the TOR pathway has benefited from the allosteric TOR inhibitor rapamycin and more recently from ATP-competitive TOR inhibitors (asTORis) that circumvent the limitations of rapamycin. The asTORis, developed for medical application, inhibit TOR complex 1 (TORC1) more efficiently than rapamycin and also inhibit rapamycin-resistant TORCs. This review presents knowledge on TOR inhibitors from the mammalian field and underlines important considerations for plant and algal biologists. It discusses the use of rapamycin and asTORis in plants and algae and concludes with guidelines for physiological studies and genetic screens with TOR inhibitors.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.