- Research article
- Open Access
Identification, characterization and functional analysis of grape (Vitis vinifera L.) mitochondrial transcription termination factor (mTERF) genes in responding to biotic stress and exogenous phytohormone
BMC Genomics volume 22, Article number: 136 (2021)
Mitochondrial transcription termination factor (mTERF) is a large gene family which plays a significant role during plant growth under various environmental stresses. However, knowledge of mTERF genes in grapevine (Vitis L.) is limited.
In this research, a comprehensive analysis of grape mTERF (VvmTERF) genes, including chromosome locations, phylogeny, protein motifs, gene structures, gene duplications, synteny analysis and expression profiles, was conducted. As a result, a total of 25 mTERF genes were identified from the grape genome, which are distributed on 13 chromosomes with diverse densities and segmental duplication events. The grape mTERF gene family is classified into nine clades based on phylogenetic analysis and structural characteristics. These VvmTERF genes showed differential expression patterns in response to multiple phytohormone treatments and biotic stresses, including treatments with abscisic acid and methyl jasmonate, and inoculation of Plasmopara viticola and Erysiphe necator.
These research findings, as the first of its kind in grapevine, will provide useful information for future development of new stress tolerant grape cultivars through genetic manipulation of VvmTERF genes.
In eukaryotes, genetic information is not only stored in the nucleus, but also in organelle genomes such as mitochondria and chloroplasts. However, these organelles’ gene pool has dramatically reduced during their evolution, which is due to the loss of their genes, and continuous transfer of organelle-nuclear genes [1,2,3]. In living organisms, the organelle gene expression system largely depends on nuclear-coding proteins, which include RNA polymerase, sigma factor, as well as specific RNA maturation factors [4,5,6,7,8]. Meanwhile, some organelle protein families including, PPRs, HAT, OPRs and mTERFs which have similar modular structures consisting of repetitive helical motifs also play an important role in their gene expression mechanism [4, 9].
Mitochondrial transcription termination factor (mTERF) genes comprise a large family which plays an essential role in the regulation of mitochondrial gene transcription . MTERF proteins possess a unique repetition of 30 amino acids residue, which enables them to recognize and bind to specific sites on mitochondrial genome known as typical mTERF motif . In human for instance, mTERFs comprises a proline at position 8, 11, 18 and 25. Therefore, the motif is conserved for leucine or hydrophobic amino acids, indicating that there are at least three leucine motifs in the mTERFs . Former research indicated that mTERF proteins could have multiple biological functions of intracellular regulation. For instance, human mTERF1, with 342 amino acids in length, can bind to 28 nucleotide sequences downstream of the 3’end of 16SrRNA, leading to the termination of mitochondrial gene transcription [13, 14]. The mTERF1 protein possess regulatory function for transcriptional initiation of mitochondrial rDNA and mitochondrial DNA replication . In addition, mTERF2 protein showed a significant downregulation of mitochondrial transcription level in vitro, suggesting that mTERF2 protein may affect mitochondrial transcription by binding with regulatory activators of mtDNA transcription initiation .
In recent years, plant mTERF genes and their roles in mitochondrial gene expression regulation have received a good deal of attention . Bioinformatic analysis shows that these mTERF genes are a large and complex family existing in metazoans and plants . There are at least identified 35 mTERFs in Arabidopsis thaliana, mainly located in mitochondria or chloroplasts, participating in abiotic stresses [19, 20]. For instance, the seed germination rate of mterf1 (soldat 10) mutant was considerably lower than wild type under the same condition . Over expression of the AtmTERF5 (MDA1) gene affects the germination rate of transgenic lines under simulated drought stress as higher germination rate was observed under mannitol treatment . Besides, A. thaliana mterf9 mutant was insensitive to ABA treatment. Under the treatment of NaCl and ABA, the root growth retardation of mTERF9 mutants displayed the phenotype of short root and lighter fresh weight compared to the wild type . Furthermore, previous studies in maize showed that ZmmTERF4 protein can co-immunoprecipitated with multiple chloroplast introns leading to the disruption of splicing in Zm-mterf4 mutants, indicating its key role in meditating the communication between organelle and the nucleus . The evidence expands the functional knowledge of the mTERF family.
As a large economic worth fruit crop , grapevine is an important candidate for identifying stress resistance genes to leading to better grape quality. At present, the basic structure and preliminary functions of mTERF family proteins have been continuously explored, but their detailed functions and regulation mechanisms under different stresses still remain unknown. This study introduces the members of the grape mTERF gene family (VvmTERF) and determine their potentiality in stress resistance, aiming to afford an essential information of the grape mTERF gene family and providing a resource for functional research in grape biology study.
Identification of mTERF genes in grape genome
MTERF genes in the grape genome were identified by BLASTP with HMMER 3.0  searching key domain mTERF PFAM file (PF02636) and previous reports [20, 27]. A total of 25 grape mTERF genes were identified, which were named as VvmTERF1-VvmTERF25 according to sequence of their chromosomal locations (Table 1). A high conserved mTERF domain was found in all the VvmTERF proteins.
Phylogenetic analysis and classification of grape mTERF genes
In order to evaluate the evolutionary relationship of VvmTERF gene family, a total of 91 mTERF genes from Arabidopsis (35), maize (31) and grape (25) genomes were collected for a phylogenetic tree construction using MEGA5.0 software (Fig. 1 and Figure S1). Detailed sequence information of Arabidopsis and maize mTERF genes were obtained from a previous study . The tree topology result demonstrated that nine groups (Clade I–IX) were classified according to homologous genes of maize and Arabidopsis. Of the 25 VvmTERF genes, Clade VII contained 7 genes, the most among all the clades, while other clades had 1 to 5 members, respectively. One grape mTERF gene, VvmTERF24, belonged to Clade I where 2 members were identified in Arabidopsis and in maize [20, 28]. It is worth noting that the well functional characterized mTERF genes from Arabidopsis, such as SOLDAT10 (AtmTERF1, AT2G03050), BSM/RUG2 (AtmTERF4, AT4G02990), and SHOT1 (AtmTERF18, AT3G60400) were distributed in group II, IV and VI, respectively. Meanwhile, a certain of grape mTERF genes belong to these groups, indicated their close evolutionary relationships with Arabidopsis mTERF genes from the same group.
Exon–intron structure analysis of VvmTERF genes
Structure analysis on the exon and intron boundaries of the VvmTERF genes will provide important clues as they played significant roles in evolution of various gene families. The number of exons per grape mTERF gene ranged from 1 to 22 (Fig. 2). Among them, VvmTERF20 had the highest number of exons of 22, followed by VvmTERF14 (10), VvmTERF16 (7), VvmTERF18 (6), VvmTERF9 (6), VvmTERF24 (6) and VvmTERF4 (6), while VvmTERF3, VvmTERF8, VvmTERF11–13 and VvmTERF21 had only one exon each. These results indicated that during the long evolution of VvmTERF gene family, both exon loss and gain have occurred, which might lead to diversified function among the other closely related mTERF genes. In clade IV, for instance, the number of exons was quite large, ranging from three to ten, while the genes in clade I and IX had a relatively smaller number, ranging from one to six exons. It is interesting that VvmTERF7, 8, 11 and 12 demonstrated similar exon/intron structures and came from the same clade while most VvmTERF genes showed distinct structures. This difference in exon/intron patterns might be resulted from a series of gene replication events.
Conserved motifs and subcellular localization analysis of VvmTERF
Searching for putative conserved motifs in grape mTERF proteins analysis was conducted via Pfam  and SMART  databases. In order to predict the potential motifs in the putative grape mTERF gene family gene sequences, the MEME (Multiple Em for Motif Elicitation) program  was used and 15 mTERF motifs in grape were identified and clustered (Fig. 2 and Table 2) using the ClustalW 2.0 program . Among all, class VII sequences had more than 10 mTERF motifs, and clade IX mTERF sequences showed 5–8 mTERF motifs. Identified in human mTERF proteins previously , conversed mTERF motifs containing repeats of leucine zipper-like heptad X3LX3 structure was also found in grape mTERF motifs (Table 2), suggesting that fundamental structures and functions of mTERF proteins in Vitis might be similar to human mTERF proteins.
Aiming to find predicted motifs shared among related proteins within the grape mTERF gene family, the MEME database program  was performed. As shown in Fig. 2, a total of 15 motifs were discovered in these 25 proteins. Among them, motifs 2 and 8 were found in most grape mTERF proteins. Motif sequences comparison with PFAM mTERF domain alignment revealed that motifs 1, 4 and 5 partly covered the PFAM mTERF domain (PF02536), and motif 5 belonged to specific organelle-targeting mTERF proteins, such as the group IV grape mTERF proteins (Fig. 2). It is highly probabe that group-specific motifs lead to characteristic functions in various life activities.
In plants, the subcellular localization of a protein is closely related to its biological function. Table 3 indicates the predicted cellular location of VvmTERF proteins for future functional research. Based on protein sequence, subcellular localization prediction demonstrated that there are 12 VvmTERFs associated with chloroplasts or mitochondria, which may imply that functions of VvmTERF proteins are related to these organelles.
Synteny analysis of VvmTERF and AtmTERF genes
Arabidopsis is a well-studied model species which can provide available genomic information to a less-studied species through genomic comparison method [33, 34]. As showed in Fig. 3, a large-scale syntenies study containing 6 pairs of grape and Arabidopsis mTERF genes were recognized. Grape orthologues including VvmTERF2, VvmTERF6, VvmTERF13, VvmTERF15, VvmTERF24 and VvmTERF25 displayed synteny location with Arabidopsis mTERF genes AtmTERF6, AtmTERF4, AtmTERF19, AtmTERF10, AtmTERF9 and AtmTERF17, respectively (Table S2). The number of synteny results indicated that several mTERF genes might arise before the divergence of Arabidopsis and grape lineages, and also suggested that partial deletion of the grape genes might occur in specific syntenic locations during genome evolution.
Cis-element analysis of grape mTERF gene promoters
To understand the possible regulatory mechanism of VvmTERF genes in multiple stress responses and functions in chloroplast and mitochondrion, a 2-kb sequence upstream of the precited transcription start site (TSS) of each VvmTERF gene was analyzed by the PlantCARE database. Meanwhile, Actin1 was chosen in grape genome as the housekeeping gene (Fig. 4). The sequences of VvmTERF gene promoters were found to contain various hormone regulation-related cis-elements such as those responsive to auxin, MeJA (Methyl Jasmonate), gibberellin, abscisic acid and salicylic acid. In addition, various defense and stress-related elements were also observed. These elements included light and wound responsive elements, osmotic stress-related elements, and low temperature and drought responsive elements.
Analysis of expression profiles among the grape mTERF genes in different tissues and organs
To discover the potential function of VvmTERF proteins during different stages of grape development, the tissue/organ-specific gene expression profiles of VvmTERF were analyzed in the V. vinifera cv. Corvina global gene expression atlas from the GEO DataSet (GSE36128). This dataset contained expression information of 54 sample tissues and organs in different developmental stages acquired by microarray database (Fig. 5). The results showed that some VvmTERF genes such as VvmTERF6, 9, 11 and 23 displayed similar expression patterns in different tissues and organs, while other VvmTERF genes like VvmTERF1, 3, 10 and 16 demonstrated tissue/organ-specific expression patterns, suggesting multiple roles played by these VvmTERF genes family in grapevine.
Expression patterns of VvmTERF genes under different exogenous hormone treatments
To explore potential stress-related genes characterized in this research, plant signaling and regulatory hormones including, ABA, MeJA, SA and Eth were used for exogenous treatment . Interestingly, almost all these VvmTERF gene expressions were influenced by exogenous hormone treatments (Fig. 6 and Figure S2). For instance, after the ABA treatment, a total of 13 VvmTERF genes displayed multiple degrees of up regulation while 8 genes were down regulated. MeJA treatment led to the expression increase of 17 VvmTERF genes and decrease on 7 genes. However, the expression patterns under SA and Eth treatments were different from those regulated by ABA and MeJA as more down regulated genes were observed. A total of 5 VvmTERF genes were up regulated and 12 were down regulated by SA, while 7 were up regulated and 14 were down regulated by exogenous Eth hormone treatment. According to the semi-quantitative RT-PCR result, VvmTERF2, VvmTERF6, VvmTERF16, VvmTERF22 and VvmTERF23, which were downregulated by the Eth treatment, displayed significant upregulation under MeJA treatment, indicating an existence of different regulatory networks among these phytohormones. Meanwhile, the existence number of responsive cis-elements of the promoter region may also play a role in the gene expression regulation.
Expression profiles of VvmTERF genes in response to biotic infections
In order to adapt to changing environments, the ability to tolerate diverse array of stresses is an essential evolutionary trait in plant kingdom. Identification and functional analysis of genes involved in biological signal transduction pathways is of great significance in providing a fundamental information for plant development and stress responses. To investigate their role in responding to biotic stress, express analysis of the 25 VvmTERF genes were conducted in potted ‘Thompson Seedless’ grapevines in greenhouse after inoculating with powdery mildew (PM) and downy mildew (DM) pathogens. As shown in Fig. 6, most VvmTERF genes demonstrated a tendency of downward expression after the inoculation. For instance, the expression of clade VII genes-VvmTERF5, 7–12-decreased in both E. necator and P. viticola treatments, while VvmTERF7 and VvmTERF10 genes have slightly decreased after P. viticola inoculation (Fig. 6, Fig. S3 and S4). Besides, the expression level of VvmTERF6, VvmTERF14 and VvmTERF19 held steady under both biotic treatments. On the other hand, VvmTERF11, VvmTERF17 and VvmTERF21 displayed an increasing trend in both PM and DM treatments in comparison with the control. Based on semi quantitative RT-PCR analysis, three grape mTERF genes (VvmTERF2, VvmTERF4 and VvmTERF20) were chosen for further detailed analysis using real-time qPCR (Figs. 7 and 8). The qPCR results were consistent with the those obtained by semi quantitative RT-PCR.
Widely identified in metazoans and plants, mitochondrial transcription termination factors (mTERFs) can regulate the expression of organelle genes at different levels [36, 37]. Previous research showed that mTERF plays a significant regulatory role in the evolution of mitochondrial genes, leading to a fundamental role in mitochondrial functionality, biological evolution, gene diagnosis and treatment . In plants, the expression of mitochondrial genes are fundamental for various plant biological functions . To fully explore the functions of grape mTERF genes, it is essential to identify and characterize mTERF genes in grape genome. In the current study, 25 grape mTERF genes were identified and their structures, evolutionary traits and expression patterns in responding to biotic stresses and hormone treatments were also analyzed.
Identification, annotation and evolution of VvmTERF genes
In this study we investigated mTERF genes in grapevine and found 25 genes belonging to mTERF family that can be mapped onto the sequenced grape genome. The number of grape mTERF genes was less than that of Arabidopsis, despite the fact that grape has a much larger genome which would suggest the un-sequenced genomic gaps or mis-annotated genes of grape genome.
According to phylogenetic analysis, mTERF genes of Arabidopsis have been classified into eight groups . In this study, a constructed phylogenetic tree which gathered the mTERF proteins from Arabidopsis, maize and grape had close topological framework to the tree constructed. Based on phylogenetic classification, grape mTERF genes were classified into nine groups. The number of clades I and VII of VvmTERF genes are small that could be due to a different pattern of duplication events. Furthermore, most of the VvmTERF genes were closely related to AtmTERF genes, which is in accordance with the fact that both grape and Arabidopsis are eudicots and exiting an appearance of close evolutional distance. As a result of highly conserved features, those mTERF genes which contain the same subclass may display similar functions. Multiple Arabidopsis mTERF genes functions have been tested, for instance, AtmTERF5, 9, 10 and 11 have functions on the resistance to salt and osmotic stress, and AtmTERF5, 9 and 10 also play roles in responding to ABA regulation . Although the Arabidopsis mTERF genes can provide predicted characterization for VvmTERF genes, the functional analysis of VvmTERF genes homologs still need more detailed experimental demonstration.
Expansion and synteny analysis of VvmTERF gene family
During evolution, tandem, segmental and whole genome duplications have been commonly found in many organisms [34, 41]. In our study, based on the chromosome locations, motifs and sequences, we concluded that some of the VvmTERF genes such as VvmTERF3–4 and VvmTERF8–11 might arise by tandem duplications. This leads to rearrangement and extension of the mTERF gene family, reported in other plant species [20, 28, 42], such as Arabidopsis genome . Some of the Arabidopsis mTERF genes were considered to be generated by tandem duplications and one block duplication [19, 20]. On the other hand, it has been demonstrated that grape has gone through whole-genome duplication events distinctly . Therefore, tandem and segmental duplications could probably contribute for most gene extensions, although there are different opinions on the exact nature and timing of these events [44, 45]. Similarly, tandem and segmental duplications have probably played a key role for grape mTERF gene expansions and their structural and functional diversity. Therefore, according to the comparison with respective orthologs of the model plant Arabidopsis, the putative functions of grape mTERF genes can be speculated. The current work analyzed the tandem duplication events of the 25 grape mTERF genes on the 13 grape chromosomes based on the research techniques of Holub , within 200 kb length on all chromosomes containing more than two genes that will be deemed to regard as a tandem duplication event.
In order to research a less-studied species, we often use genomic comparison method which could effectively transfer genomic knowledge obtained from a well-studied model species (e.g. Arabidopsis) to a less studied organism [33, 47]. In this research, as seen in Fig. 3, synteny analysis of the grape and Arabidopsis genomes identified six pairs of mTERF genes (VvmTERF2-AtmTERF6, VvmTERF6-AtmTERF4, VvmTERF13-AtmTERF19, VvmTERF15-AtmTERF10, VvmTERF24-AtmTERF9 and VvmTERF25-AtmTERF17) located in syntenic genomic regions (Fig. 3). Accompanied by selected genes loss, Arabidopsis and grape genomes have also gone through multiple and crucial chromosomal rearrangements and fusion processes during their evolution, which results in introduction of genes mismatches on chromosomes. In this case, we can deduce that the mTERF genes of grape and Arabidopsis in the same linear region may have a common ancestor. The first identified mTERF gene in Arabidopsis is AtmTERF1 (SOLDAT10), which is mainly involved in fluorescent phenotype and O2− signaling cell death . Furthermore, the AtmTERF4 (BSM/RUG2) gene is crucial for plant development. The rug2–1 and bsm mutant are unable to grow compared with the wild type plant. Sequence analysis revealed that VvmTERF6 was homologous to AtmTERF4 which might imply that the VvmTERF6 may have similar function in regulating plant development.
VvmTERF genes response to hormone treatments and biotic stresses
In a previous study, Linder  firstly described the mTERF gene family in plants leading to number of studies on identification of this family in various plant species. Identification and functional analysis of mTERF gene family in maize and pepper are among the best examples [28, 42]. However, information about mTERF functions in plants is still rather limited and needs further investigations. In this report, we analyzed the expression patterns of VvmTERF after exposing to various biotic and abiotic stresses including pathogens and phytohormones. Under these different treatments, the VvmTERF genes showed various expression patterns, implying their crucial role in plant growth and response to environmental stresses.
Among the phytohormones, it has been reported that ABA is widely involved in various biological function in response to biotic and abiotic stress [48,49,50], while plant hormone Eth, SA and MeJA have synergistic effects on biological stress signals after pathogen infection . Our results revealed that the grape mTERF genes responded to different plant hormones, which were consistent with former studies in other species such as maize and capsicum mTERF genes [28, 42].
A total of 35 mTERF genes were identified in A. thaliana, among which 6 mTERF genes were discovered specifically functional. Gene SOLDAT10 and SHOT1 for instance can respond to abiotic stresses [21, 52], as well as gene TWIRT1 has meristem function  and gene RUG2 is associated with leaf morphology . VvmTERF6 is a homolog of AtmTERF4 that has functions of organelles development and photoautotrophic growth. In Fig. 5, expression alteration of VvmTERF6 gene in the development stages of grape showed distinct upregulated pattern, indicating potential function during growth periods. It is interesting to note that 11 light responsive cis-elements were found from the promoter region of VvmTERF6 gene in Fig. 4, suggesting a role in grape growth and development. Named MDA1, AtmTERF5 played a role in rooting process of model plant Arabidopsis and showed responses to abiotic treatment. On account of reducing sensitivity to hormone abscisic acid, mutant mda1 seedlings are exhibited insensitiveness to osmotic and salt stresses, while grown Arabidopsis mda1 plants demonstrate reduced tolerance to cold, salinity or ABA treatment . Hence, since this gene family has expanded in plants, scientists found that this family is a suitable candidate for many mutant-generating strategies in order to determine potent phenotypic, or even showed embryo-lethal feature. In addition to the results acquired from the mutant analysis, previously published mTERF genes expression data which had also clearly showed the potential role for mTERFs in plant stress response. In the report, ZmmTERF12 and ZmmTERF28 were down- and up-regulated by ABA hormone treatment, respectively . From the same clade with ZmmTERF28, VvmTERF25 also displayed quick response to ABA after 0.5 h treatment. It was found that abscisic acid responsive element appeared in the promoter region of VvmTERF25, implying the reason that ABA lead to a positive regulation of VvmTERF25. Taken together, we analyzed the response of VvmTERF genes to various plant hormones and found potential key targets for enhancing grape stress resistance, which provided basic information for future studies on function of VvmTERF genes and its related signal transduction network.
Previous research demonstrated that some of the plant mTERF genes such as SOLDAT10 and SHOT1 participate in ROS scavenging process . Given this, the role of mTERF genes in response to biotic stresses needs more information. Further we conducted the related experiment here and the results confirmed our hypothesis before. In this study, we found that some grape mTERF genes responded to powdery or downy mildew treatments. For example, VvmTERF22 was up regulated by PM at 12 h, and VvmTERF21 positively responded to downy mildew infection. We detected several VvmTERF genes expression levels which are actively responsive to PM or DM treatments, indicating these genes might display significant characters in the grape protection, but further research is needful to demonstrate that they participate in biotic stress responses in grapevine.
The family of transcriptional activators encoded by the mTERF genes is widely found in plants and animals. Although significant progresses have been made in identifying mTERF genes in model plant species, little information of mTERF genes has been known in fruit crops. In this study, we identified a total of 25 VvmTERF genes in the grape genome, and also investigated their structural, phylogenetical and syntenic features. Through comparative analysis of homology between grape and Arabidopsis, it is found that several mTERF genes of grape and Arabidopsis are located in the homologous regions, indicating that they may present close evolutionary relationship. Expression analysis of VvmTERF genes showed that multiple genes could respond to different biological stresses and hormone treatments. Results from this study have paved the way for future research to investigate roles of VvmTERF genes on disease resistance of grapevines.
Identification and annotation of grape mTERF genes
Conserved mTERF domains were first used to detect grape genes in mTERF HMM (Hidden Markov Model) file (PF02536) from the Pfam database  using the HMMER 3.0 package . The domains were further used as the query to search the GenBank nonredundant protein and the Grape Genome Database (http://www.rosaceae.org/projects/grape_genome) using the BLAST program. All mTERF proteins with an E value < 0.01 were collected and the domains were manually checked in each candidate VvmTERF gene.
Multiple sequence alignment, phylogenetic analysis and classification
A total of 25 VvmTERF genes containing mTERF core domains were identified. The CLUSTALX software was then used for multiple sequence alignment analysis including grape mTERF genes and those from Arabidopsis (AtmTERF) and maize (ZmmTERF) . Based on the neighbor-joining method and maximum likelihood method for 1000 bootstrap replicates, a phylogenetic tree was constructed using MEGA5.0 software . The VvmTERF genes were classified into clades ground on multiple sequence alignments with those AtmTERF and ZmmTERF genes.
MEME motifs, conserved sequences and subcellular localization analysis of grape mTERF proteins
The identification of known conserved motifs in grape mTERF proteins was conducted by BLAST against the SMART  and Pfam  database searching. The potential motifs in the putative mTERF family gene sequences were predicted by Multiple Em for Motif Elicitation (MEME) software  with the parameters as follow: the optimum width of every single motif distributed between 6 to 50, and the maximum number of motifs to find was 15. After that, the collection and cluster of mTERF motifs from grape mTERF proteins were conducted using the ClustalW 2.0 program software , and graphical representation of amino acid residues was arranged by TBtools .
Exon–intron structure analysis, synteny analysis and gene duplication
According to alignments of grape mTERF gene coding sequences and their respective full-length sequences, the exon-intron structure was determined on Grape Genome Browser: http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/. And the online program Gene Structure Display Server (GSDS: http://gsds.cbi.pku.edu.ch)  was carried out to obtain relative diagrams. Besides, the definition of mTERF genes with tandem duplication events were adjacent homologous genes on a single chromosome, while gene duplication events among diverse chromosomes were defined as segmental duplications . The specific physical location of each VvmTERF gene on its individual chromosome determines whether it is considered in a tandem duplication event. Therefore a synteny analysis map of grape mTERF gene was constructed via the syntenic blocks, and a further synteny analysis between grape and AtmTERF genes was acquired from the Plant Genome Duplication Database . The generation of related diagrams were illustrated using Circos version website (http://circos.ca/).
Cis-element analysis of grape mTERF genes promoter
The 2000 bps upstream promoter sequence of each grape mTERF gene coding regions were obtained from the grape genome database (https://wwwdev.genoscope.cns.fr/vitis). PlantCARE online analysis program (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) was used to find the predicted cis-element.
Expression profiles of grape mTERF gene family in different tissues and organs
The expression profile of grape mTERF gene family was confirmed in a V. vinifera cv ‘Corvina’ (clone48) gene expression atlas of various organs at different developmental stages. Microarray data were collected from the NCBI gene expression omnibus (GEO) datasets under the series entry GSE36128 (http://www.ncbi.nlm.nih.gov/geo/) . The mean expression value of grape mTERF genes in all tissues and organs were analyzed and detailed displayed by Multiple Experiment Viewer software (MeV) . Measured using RNA-Seq data, the expression patterns of VvmTERF genes in various berry developmental stages were gained from gene expression omnibus (GEO) database of NCBI (GSE77218) .
Plant materials and stress treatments
To validate the expression regulation of grape mTERF genes under abiotic and biotic stresses, grape leaves and organs were obtained from V. vinifera ‘Thompson Seedless’ grape grown in a greenhouse. When the shoots of the grapevines reached 30 cm long with fully expanded young leaves, the plants were subjected to hormone treatment. Hormone treatment was carried out on grape leaves under similar growth condition, spraying with 100 μM salicylic acid (SA), 300 μM abscisic acid (ABA), 50 μM methyl jasmonate (MeJA), and 0.5 g/l ethylene (Eth), respectively . Leaves from the treated vines were collected at 0.5, 1, 3, 6, 12, 24, and 48 h post treatment. Grape leaves sprayed with sterile water were used as the negative control.
In term of biotic stress, powdery mildew (Erysiphe necator) and downy mildew (Plasmopara viticola) pathogens were used to inoculate young leaves of V. vinifera ‘Thompson Seedless’ following previous protocol [64, 65]. Leaves were sampled at 6, 12, 24, 48, 72, 96, and 120 h post inoculation and untreated leaves were collected as the negative control. At each time point of all treatments, nine leaves from three separate plants were homogenized, and the treatments were conducted three times independently. These grape leaves were frozen in liquid nitrogen and stored at − 80 °C for further use.
Semiquantitative PCR and real-time quantitative PCR analysis
The total RNA was extracted as described by Zhang et al. . The genomic DNA was digested using RNase-free DNase I kit (OMEGA Bio Inc., USA). The grape Actin1 gene (GenBank Accession number AY680701) and EF1-α (GenBank Accession number EC931777) gene were chosen as housekeeping genes and amplified with the primers showing in Table S1, which also includes Gene-specific primers for the 25 VvmTERF genes. For the semi-quantitative reverse transcription-PCR experiment, the volume of reaction system was 20 μL which includes 1 μL cDNA template, 1 μL forward and reverse gene-specific primers (10 μM), 10 μL PCR Master Mix (Qingke Biotech Co. Ltd., Shanghai, China) and 8 μL sterile water, the specific proportion and program were set according to the PCR Master Mix instruction book. Each PCR reaction was conducted in duplicate. The Gene Tools software was used for quality control of semi-quantitative PCR results, further log-transformed values of the relative expression patterns of VvmTERF genes under various phytohormone and biotic stresses treatment were used to perform hierarchical cluster using Mev software .
Quantitative real-time PCR analysis was performed with an IQ5 real-time PCR instrument (Bio-Rad, Hercules, CA, USA). All reactions were performed in triplicate with a reaction system of 20 μL including 1 μL specific primers (10 μM), 1 μL cDNA, 10 μL SYBR green (Yeasen Biotech Co Ltd., Shanghai, China), and 8 μL sterile water, the specific proportion was on the instruction book as well as the PCR parameters. The grape Actin1 gene (GenBank Accession number AY680701) was chosen as the housekeeping gene. The expression levels of grape mTERF genes were analyzed using IQ5 software with the normalized expression method. The t-test was conducted using the SPSS software (SPSS 17.0, Chicago, IL, USA, *P < 0.05, **P < 0.01).
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional files.
Mitochondrial transcription termination factor
Vitis vinifera mitochondrial transcription termination factor
Pentatricopeptide repeat proteins
Half atetratricopeptide proteins
Octotricopeptide repeat proteins
Reverse transcription quantitative PCR
Kleine T, Voigt C, Leister D. Plastid signalling to the nucleus: messengers still lost in the mists? Trends Genet. 2009;25(4):185–90.
Martin W, Rujan T, Richly E, Hansen A, Cornelsen S, Lins T, Leister D, Stoebe B, Hasegawa M, Penny D. Evolutionary analysis of Arabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus. P Natl Acad Sci USA. 2002;99(19):12246–51.
Timmis JN, Ayliffe MA, Huang CY, Martin W. Endosymbiotic gene transfer: organelle genomes forge eukaryotic chromosomes. Nat Rev Genet. 2004;5(2):123–U116.
Hammani K, Barkan A. An mTERF domain protein functions in group II intron splicing in maize chloroplasts. Nucleic Acids Res. 2014;42(8):5033–42.
Barkan A. Studying the structure and processing of chloroplast transcripts. Methods Mol Biol. 2011;774:183–97.
Rackham O, Mercer TR, Filipovska A. The human mitochondrial transcriptome and the RNA-binding proteins that regulate its expression. Wiley Interdiscip Rev RNA. 2012;3(5):675–95.
Stoppel R, Meurer J. The cutting crew - ribonucleases are key players in the control of plastid gene expression. J Exp Bot. 2012;63(4):1663–73.
Yagi Y, Shiina T. Recent advances in the study of chloroplast gene expression and its evolution. Front Plant Sci. 2014;5:61.
Kleine T, Leister D. Emerging functions of mammalian and plant mTERFs. Biochim Biophys Acta. 2015;1847(9):786–97.
Linder T, Park CB, Asin-Cayuela J, Pellegrini M, Larsson NG, Falkenberg M, Samuelsson T, Gustafsson CM. A family of putative transcription termination factors shared amongst metazoans and plants. Curr Genet. 2005;48(4):265–9.
Roberti M, Bruni F, Polosa PL, Manzari C, Gadaleta MN, Cantatore P. MTERF3, the most conserved member of the mTERF-family, is a modular factor involved in mitochondrial protein synthesis. Bba-Bioenergetics. 2006;1757(9–10):1199–206.
Roberti M, Polosa PL, Bruni F, Manzari C, Deceglie S, Gadaleta MN, Cantatore P. The MTERF family proteins: mitochondrial transcription regulators and beyond. Bba-Bioenergetics. 2009;1787(5):303–11.
Kruse B, Narasimhan N, Attardi G. Termination of transcription in human mitochondria: identification and purification of a DNA binding protein factor that promotes termination. Cell. 1989;58(2):391–7.
FernandezSilva P, MartinezAzorin F, Micol V, Attardi G. The human mitochondrial transcription termination factor (mTERF) is a multizipper protein but binds to DNA as a monomer, with evidence pointing to intramolecular leucine zipper interactions. EMBO J. 1997;16(5):1066–79.
Martin M, Cho JY, Cesare AJ, Griffith JD, Attardi G. Termination factor-mediated DNA loop between termination and initiation sites drives mitochondrial rRNA synthesis. Cell. 2005;123(7):1227–40.
Wenz T, Luca C, Torraco A, Moraes CT. mTERF2 regulates oxidative phosphorylation by modulating mtDNA transcription. Cell Metab. 2009;9(6):499–511.
Quesada V. The roles of mitochondrial transcription termination factors (MTERFs) in plants. Physiol Plantarum. 2016;157(3):389–99.
Robles P, Micol JL, Quesada V. Unveiling plant mTERF functions. Mol Plant. 2012;5(2):294–6.
Babiychuk E, Vandepoele K, Wissing J, Garcia-Diaz M, De Rycke R, Akbari H, Joubes J, Beeckman T, Jansch L, Frentzen M, et al. Plastid gene expression and plant development require a plastidic protein of the mitochondrial transcription termination factor family. P Natl Acad Sci USA. 2011;108(16):6674–9.
Kleine T. Arabidopsis thaliana mTERF proteins: evolution and functional classification. Front Plant Sci. 2012;3:233.
Meskauskiene R, Wursch M, Laloi C, Vidi PA, Coll NS, Kessler F, Baruah A, Kim C, Apel K. A mutation in the Arabidopsis mTERF-related plastid protein SOLDAT10 activates retrograde signaling and suppresses 1O(2)-induced cell death. Plant J. 2009;60(3):399–410.
Robles P, Micol JL, Quesada V. Arabidopsis MDA1, a nuclear-encoded protein, functions in chloroplast development and abiotic stress responses. Plos One. 2012;7(8):e42924.
Robles P, Micol JL, Quesada V. Mutations in the plant-conserved MTERF9 alter chloroplast gene expression, development and tolerance to abiotic stress in Arabidopsis thaliana. Physiol Plantarum. 2015;154(2):297–313.
Kamel H. Alice BJNAR. An mTERF domain protein functions in group II intron splicing in maize chloroplasts. 2014;8:8.
Kreamer RJA. US table grape exports scoring big in world markets; 1995.
Eddy SR. Hidden Markov models. Curr Opin Struct Biol. 1996;6(3):361–5.
Inal B, Ilhan E, Büyük İ, Altıntaş S. Transcriptome wide characterization of water deficit responsive grape mTERF transcription. J Plant Biochem Biot. 2020;29(1):102-13.
Zhao YX, Cai MJ, Zhang XB, Li YR, Zhang JH, Zhao HL, Kong F, Zheng YL, Qiu FZ. Genome-wide identification, evolution and expression analysis of mTERF gene family in Maize. Plos One. 2014;9(4):e94126.
Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40(D1):D290–301.
Letunic I, Doerks T, Bork P. SMART 7: recent updates to the protein domain annotation resource. Nucleic Acids Res. 2012;40(D1):D302–5.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren JY, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:W202–8.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.
Lyons E, Pedersen B, Kane J, Alam M, Ming R, Tang HB, Wang XY, Bowers J, Paterson A, Lisch D, et al. Finding and comparing Syntenic regions among Arabidopsis and the Outgroups papaya, poplar, and grape: CoGe with Rosids. Plant Physiol. 2008;148(4):1772–81.
Guo CL, Guo RR, Xu XZ, Gao M, Li XQ, Song JY, Zheng Y, Wang XP. Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family. J Exp Bot. 2014;65(6):1513–28.
Fujita M, Fujita Y, Noutoshi Y, Takahashi F, Narusaka Y, Yamaguchi-Shinozaki K, KJCoipb S. Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Curr Opin Plant Biol. 2006;9(4):436–42.
Nedialkov YA, Opron K, Assaf F, Artsimovitch I, Kireeva ML, Kashlev M, Cukier RI, Nudler E, Burton ZF. The RNA polymerase bridge helix YFI motif in catalysis, fidelity and translocation. Bba-Gene Regul Mech. 2013;1829(2):187–98.
Leister D, Kleine T. Extending the repertoire of mTERF proteins with functions in Organellar gene expression. Mol Plant. 2020;13(6):817–9.
DiMaggio C, Markenson D, Loo GT, Redlener I. The willingness of US emergency medical technicians to respond to terrorist incidents. Biosecur Bioterror. 2005;3(4):331–7.
Karaca M, Ince AG, Ay ST, Turgut K, Onus AN. PCR-RFLP and DAMD-PCR genotyping for Salvia species. J Sci Food Agr. 2010;88(14):2508–16.
Xu DR, Leister D, Kleine T. Arabidopsis thaliana mTERF10 and mTERF11, but Not mTERF12, are involved in the response to salt stress. Front Plant Sci. 2017;8:1213.
Xu G, Guo C, Shan H, Kong H. Divergence of duplicate genes in exon-intron structure. Proc Natl Acad Sci U S A. 2012;109(4):1187–92.
Tang B, Xie L, Yi T, Lv J, Yang H, Cheng X, Liu F, Zou X. Genome-wide identification and characterization of the mitochondrial transcription termination factors (mTERFs) in Capsicum annuum L. Int J Mol Sci. 2019;21(1):269.
Maere S, De Bodt S, Raes J, Casneuf T, Van Montagu M, Kuiper M, Van de Peer Y. Modeling gene and genome duplications in eukaryotes. Proc Natl Acad Sci U S A. 2005;102(15):5454–9.
Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449(7161):463–U465.
Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D, Pindo M, Fitzgerald LM, Vezzulli S, Reid J, et al. A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS One. 2007;2(12):e1326.
Holub EB. The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet. 2001;2(7):516–27.
Zhang YC, Mao LY, Wang H, Brocker C, Yin XJ, Vasiliou V, Fei ZJ, Wang XP. Genome-wide identification and Analysis of Grape Aldehyde Dehydrogenase (ALDH) Gene Superfamily. Plos One. 2012;7(2):e32153.
Mauch-Mani B, Mauch F. The role of abscisic acid in plant-pathogen interactions. Curr Opin Plant Biol. 2005;8(4):409–14.
Davies WJ, Jones HG. Abscisic acid physiology and biochemistry. Abscisic acid: physiology and biochemistry. Bios Sci. 1992;16(2):99.
Jameson PE, Murray BG. Abscisic acid physiology and biochemistry. New Zeal J Bot. 1992;30(3):369-71.
Glazebrook J. Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu Rev Phytopathol. 2005;43:205–27.
Kim M, Lee U, Small I, Cc d F-S, Vierling E. Mutations in an Arabidopsis mitochondrial transcription termination factor-related protein enhance Thermotolerance in the absence of the major molecular chaperone HSP101. Plant Cell. 2012;24(8):3349–65.
Mokry M, Nijman IJ, van Dijken A, Benjamins R, Heidstra R, Scheres B, Cuppen E. Identification of factors required for meristem function in Arabidopsis using a novel next generation sequencing fast forward genetics approach. BMC Genomics. 2011;12.
Quesada V, Sarmiento-Manus R, Gonzalez-Bayon R, Hricova A, Perez-Marcos R, Gracia-Martinez E, Medina-Ruiz L, Leyva-Diaz E, Ponce MR, Micol JL. Arabidopsis RUGOSA2 encodes an mTERF family member required for mitochondrion, chloroplast and leaf development. Plant J. 2011;68(4):738–53.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, SJMb K. Evolution. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.
Chen C, Chen H, He Y, RJB X. TBtools, a toolkit for biologists integrating various biological data handling tools with a user-friendly interface, vol. 289660; 2018.
Yu CS, Chen YC, Lu CH, Hwang JK. Prediction of protein subcellular localization. Proteins-Structure Function and Bioinformatics. 2006;64(3):643–51.
Hu B, Jin JP, Guo AY, Zhang H, Luo JC, Gao G. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2015;31(8):1296–7.
Liu Y, Jiang H, Chen W, Qian Y, Ma Q, Cheng B, Zhu SJPGR. Genome-wide analysis of the auxin response factor (ARF) gene family in maize (Zea mays) 2011; 63(3):225–234.
Tang H, Wang X, Bowers JE, Ming R, Alam M, AHJGr P. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008;18(12):1944–54.
Fasoli M, Dal Santo S, Zenoni S, Tornielli GB, Farina L, Zamboni A, Porceddu A, Venturini L, Bicego M, Murino V, et al. The grapevine expression atlas reveals a deep Transcriptome shift driving the entire plant into a maturation program. Plant Cell. 2012;24(9):3489–505.
Saeed AI, Hagabati NK, Braisted JC, Liang W, Sharov V, Howe EA, Li JW, Thiagarajan M, White JA, Quackenbush J. TM4 microarray software suite. Method Enzymol. 2006;411:134.
Shangguan L, Mu Q, Fang X, Zhang K, JJPO F. RNA-Sequencing Reveals Biological Networks during Table Grapevine ('Fujiminori') Fruit Development. PLoS One. 2017;12(1):e0170571.
Guan X, Zhao H, Xu Y, YJP W. Transient expression of glyoxal oxidase from the Chinese wild grape Vitis pseudoreticulata can suppress powdery mildew in a susceptible genotype. Protoplasma. 2011;248(2):415–23.
Li X, Wu J, Yin L, Zhang Y, Qu J, Lu JJPP. Biochemistry. Comparative transcriptome analysis reveals defense-related genes and pathways against downy mildew in Vitis amurensis grapevine, vol. 95; 2015. p. 1–14.
The authors thank all editors and reviewers for their comments on this manuscript.
This work was supported by Shanghai Agriculture Applied Technology Development Program (Grant No. 2018-02-08-00-08-F01552); National Key R & D Program of China (2018 YFD100300); the National Natural Science Foundation of China (No. 31801830 and 31701775); the Shanghai Jiao Tong University New Scholar Start-Up Fund (AF1500068); Shanghai Municipal Commission for Science and Technology (18391900400); China Agriculture Research System (CARS-29-yc-2).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phylogenetic analysis among the grape, Arabidopsis and maize mTERF proteins. The unrooted tree was constructed using MEGA5.0 software by Maximum Likelihood method for 1000 bootstrap replicates. Three dot colors mean different species. Yellow, green and red represent maize, Arabidopsis and grape, respectively.
Expression patterns of 25 VvmTERF genes under hormone Eth, MeJA, SA and ABA treatments analyzed by semi-quantitative RT-PCR. Actin1 and EF-1α (GenBank Accession number AY680701 and EC931777) were used as internal reference genes. The upper and lower bands indicate treatment and control, respectively.
Expression patterns of 25 VvmTERF genes after inoculation of downy mildew and powdery mildew treatments analyzed by semi-quantitative RT-PCR. Actin1 and EF-1α (GenBank Accession number AY680701 and EC931777) were used as internal reference genes. The upper and lower bands indicate treatment and control, respectively.
Expression profiles of 25 VvmTERF genes under exogenous hormone and biotic treatments. Numbers in boxes represent different expression levels. 0 indicates no change, number less than zero means down-regulated expression and more than zero means up-regulated expression.
Primer sequences used for semi-quantitative RT-PCR and quantitative real-time PCR of the 25 grape mTERF genes. Primer 5.0 software were used to designed specific primers.
Synteny blocks of VvmTERF genes between the grape and Arabidopsis genomes.
About this article
Cite this article
Yin, X., Gao, Y., Song, S. et al. Identification, characterization and functional analysis of grape (Vitis vinifera L.) mitochondrial transcription termination factor (mTERF) genes in responding to biotic stress and exogenous phytohormone. BMC Genomics 22, 136 (2021). https://doi.org/10.1186/s12864-021-07446-z
- Bioinformatics analysis
- Expression profile analysis
- Grapevine (V. vinifera L.)
- mTERF family