High intraspecific diversity of Restorer-of-fertility-like genes in barley

Mitochondrial atp1 mRNA knockdown by a custom-designed pentatricopeptide repeat protein alters ATP synthase

Spontaneous mutations are rare in mitochondria and the lack of mitochondrial transformation methods has hindered genetic analyses. We show that a custom-designed RNA-binding pentatricopeptide repeat (PPR) protein binds and specifically induces cleavage of ATP synthase subunit1 (atp1) mRNA in mitochondria, significantly decreasing the abundance of the Atp1 protein and the assembled F1Fo ATP synthase in Arabidopsis (Arabidopsis thaliana). The transformed plants are characterized by delayed vegetative growth and reduced fertility. Five-fold depletion of Atp1 level was accompanied by a decrease in abundance of other ATP synthase subunits and lowered ATP synthesis rate of isolated mitochondria, but no change to mitochondrial electron transport chain complexes, adenylates, or energy charge in planta. Transcripts for amino acid transport and a variety of stress response processes were differentially expressed in lines containing the PPR protein, indicating changes to achieve cellular homeostasis when ATP synthase was highly depleted. Leaves of ATP synthase-depleted lines showed higher respiratory rates and elevated steady-state levels of numerous amino acids, most notably of the serine family. The results show the value of using custom-designed PPR proteins to influence the expression of specific mitochondrial transcripts to carry out reverse genetic studies on mitochondrial gene functions and the consequences of ATP synthase depletion on cellular functions in Arabidopsis.

Comparative Genomic Analysis of Asian Cultivated Rice and Its Wild Progenitor (Oryza rufipogon) Has Revealed Evolutionary Innovation of the Pentatricopeptide Repeat Gene Family through Gene Duplication

The pentatricopeptide repeat (PPR) gene family is one of the largest gene families in land plants. However, current knowledge about the evolution of the PPR gene family remains largely limited. In this study, we performed a comparative genomic analysis of the PPR gene family in O. sativa and its wild progenitor, O. rufipogon, and outlined a comprehensive landscape of gene duplications. Our findings suggest that the majority of PPR genes originated from dispersed duplications. Although segmental duplications have only expanded approximately 11.30% and 13.57% of the PPR gene families in the O. sativa and O. rufipogon genomes, we interestingly obtained evidence that segmental duplication promotes the structural diversity of PPR genes through incomplete gene duplications. In the O. sativa and O. rufipogon genomes, 10 (~33.33%) and 22 pairs of gene duplications (~45.83%) had non-PPR paralogous genes through incomplete gene duplication. Segmental duplications leading to incomplete gene duplications might result in the acquisition of domains, thus promoting functional innovation and structural diversification of PPR genes. This study offers a unique perspective on the evolution of PPR gene structures and underscores the potential role of segmental duplications in PPR gene structural diversity.

A unique C-terminal domain contributes to the molecular function of Restorer-of-fertility proteins in plant mitochondria

Restorer-of-fertility (Rf) genes encode pentatricopeptide repeat (PPR) proteins that are targeted to mitochondria where they specifically bind to transcripts that induce cytoplasmic male sterility and repress their expression. In searching for a molecular signature unique to this class of proteins, we found that a majority of known Rf proteins have a distinct domain, which we called RfCTD (Restorer-of-fertility C-terminal domain), and its presence correlates with the ability to induce cleavage of the mitochondrial RNA target. A screen of 219 angiosperm genomes from 123 genera using a sequence profile that can quickly and accurately identify RfCTD sequences revealed considerable variation in RFL/RfCTD gene numbers across flowering plants. We observed that plant genera with bisexual flowers have significantly higher numbers of RFL genes compared to those with unisexual flowers, consistent with a role of these proteins in restoration of male fertility. We show that removing the RfCTD from the RFL protein RNA PROCESSING FACTOR 2-nad6 prevented cleavage of its RNA target, the nad6 transcript, in Arabidopsis thaliana mitochondria. We provide a simple way of identifying putative Rf candidates in genome sequences, new insights into the molecular mode of action of Rf proteins and the evolution of fertility restoration in flowering plants.

Molecular genetics of cytoplasmic male sterility and restorer-of-fertility for the fine tuning of pollen production in crops

Cytoplasmic male sterility (CMS) is an increasingly important issue within the context of hybrid seed production. Its genetic framework is simple: S-cytoplasm for male sterility induction and dominant allele of the restorer-of-fertility gene (Rf) for suppression of S. However, breeders sometimes encounter a phenotype of CMS plants too complex to be explained via this simple model. The molecular basis of CMS provides clue to the mechanisms that underlie the expression of CMS. Mitochondria have been associated with S, and several unique ORFs to S-mitochondria are thought to be responsible for the induction of male sterility in various crops. Their functions are still the subject of debate, but they have been hypothesized to emit elements that trigger sterility. Rf suppresses the action of S by various mechanisms. Some Rfs, including those that encode the pentatricopeptide repeat (PPR) protein and other proteins, are now considered members of unique gene families that are specific to certain lineages. Additionally, they are thought to be complex loci in which several genes in a haplotype simultaneously counteract an S-cytoplasm and differences in the suite of genes in a haplotype can lead to multiple allelism including strong and weak Rf at phenotypic level. The stability of CMS is influenced by factors such as the environment, cytoplasm, and genetic background; the interaction of these factors is also important. In contrast, unstable CMS becomes inducible CMS if its expression can be controlled. CMS becomes environmentally sensitive in a genotype-dependent manner, suggesting the feasibility of controlling the expression of CMS.

Triticeae genome sequences reveal huge expansions of gene families implicated in fertility restoration

Breakthroughs in assembly of whole-genome sequencing and targeted sequence capture data have accelerated comparative genomics analyses in cereals with big and complex genomes such as wheat. This newly acquired information has revealed unexpected expansions in two large gene families linked to restoration of fertility in species that exhibit cytoplasmic male sterility. Extreme levels of copy-number and structural variation detected within and between species illustrate the genetic diversity among the family members and reveal the evolutionary mechanisms at work. This new knowledge will greatly facilitate the development of hybrid production strategies in wheat and related species.

A super PPR cluster for restoring fertility revealed by genetic mapping, homocap-seq and de novo assembly in cotton

Rf candidate genes were related to the super D05_PPR-cluster and verified to be individually nonfunctional. Restorer of fertility (Rf) genes of cytoplasmic male sterility (CMS) is commonly found to be PPR (pentatricopeptide repeat) genes, which are mostly located in a cluster of PPR genes with high similarity. Here, Homocap-seq was applied to analyze PPR clusters in 'three lines,' and we found broad variations within the D05_PPR-cluster in a restorer line and deduced that the D05_PPR-cluster was associated with fertility restoration. Genetic mapping of Rf and Homocap-seq analysis of three genotypes in the F₂ population validated that the D05_PPR-cluster was the origin of Rf. Three Rf candidates were cloned that were the most actively expressed genes in the D05_PPR-cluster in the restorer line as revealed by their high-depth amplicons. However, further transgenic experiments showed that none of the candidates could restore fertility of the CMS line independently. Then, the members of the brand-new super D05_PPR-cluster in the restorer line, containing 14 full-length PPRs and at least 13 PPR homologous sequences, were identified by long-read resequencing, which validated the effectiveness of variation and expression prediction of Homocap-seq. Additionally, we found that several PPR duplications, including 2 of the 3 Rf candidates, had undergone site-specific selection as potentially important anther development-associated genes. Finally, we proposed that multiple PPRs were coordinately responsible for the fertility restoration of the CMS line.

Small RNAs in Plant Immunity and Virulence of Filamentous Pathogens

Gene silencing guided by small RNAs governs a broad range of cellular processes in eukaryotes. Small RNAs are important components of plant immunity because they contribute to pathogen-triggered transcription reprogramming and directly target pathogen RNAs. Recent research suggests that silencing of pathogen genes by plant small RNAs occurs not only during viral infection but also in nonviral pathogens through a process termed host-induced gene silencing, which involves trans-species small RNA trafficking. Similarly, small RNAs are also produced by eukaryotic pathogens and regulate virulence. This review summarizes the small RNA pathways in both plants and filamentous pathogens, including fungi and oomycetes, and discusses their role in host-pathogen interactions. We highlight secondarysmall interfering RNAs of plants as regulators of immune receptor gene expression and executors of host-induced gene silencing in invading pathogens. The current status and prospects of trans-species gene silencing at the host-pathogen interface are discussed.

New PCR-specific markers for pollen fertility restoration QRfp-4R in rye (Secale cereale L.) with Pampa sterilizing cytoplasm

Pampa cytoplasmic male sterility phenomenon is used extensively in the rye hybrid breeding programs. It relies on sterilizing action of the cytoplasm resulting in non-viable pollen of female lines. The sterilizing effect is problematic for reversion, and efficient restores are needed. The most promising QTL is located on chromosome 4R, but other chromosomes may also code the trait. Advanced recombinant inbred lines formed bi-parental mapping population genotyped with DArTseq markers. Genetic mapping allowed the seven linkage groups to construct with numerous markers and represent all rye chromosomes. Single marker analysis and composite interval mapping were conducted to identify markers linked to the pollen fertility. Association mapping was used to detect additional markers associated with the trait. A highly significant QTL (QRfp-4R) that explained 42.3% of the phenotypic variation was mapped to the distal part of the long arm of the 4R chromosome. The markers localized in the QRfp-4R region achieve R² association values up to 0.59. The homology of the 43 marker sequences to the loci responsible for fertility restoration in other species and transcription termination factor (mTERF) linked to Rf genes was established. Ten markers were successfully converted into PCR-specific conditions, and their segregation pattern was identical to that of unconverted DArTs.

A mitochondria-localized pentatricopeptide repeat protein is required to restore hau cytoplasmic male sterility in Brassica napus

A mitochondria-localized pentatricopeptide repeat protein was identified by positional cloning and transferred into the hau CMS line, where it successfully restored fertility Cytoplasmic male sterility (CMS) is a maternally inherited trait that can be controlled by restorer-of-fertility (Rf) genes present in the nucleus. The hau CMS was identified as a new form of CMS associated with the mitochondrial transcript orf288; however, a lack of a restorer gene has limited its utilization in Brassica crops. Here, the combination of Brassica 60 K array with bulk segregant analysis and map-based cloning was used to delimit the Rfh locus to an 82.2-kb region on chromosome A09. A candidate gene encoding a mitochondria-localized pentatricopeptide repeat (PPR) protein was identified and transferred into the hau CMS line, where it successfully restored the fertility of the hau CMS plants. Furthermore, the expression analysis showed that Rfh was highly expressed in the flower buds, and the sequence analysis results implied that functional divergence between RFH and rfh could be due to 59 amino acid residue differences in the deduced protein sequences. In addition, a co-separated molecular marker was developed based on the divergent sequences between the dominant and recessive alleles. These results will help enable the heterosis of Brassica crops in the future.

The genetic basis of cytoplasmic male sterility and fertility restoration in wheat

Hybrid wheat varieties give higher yields than conventional lines but are difficult to produce due to a lack of effective control of male fertility in breeding lines. One promising system involves the Rf1 and Rf3 genes that restore fertility of wheat plants carrying Triticum timopheevii-type cytoplasmic male sterility (T-CMS). Here, by genetic mapping and comparative sequence analyses, we identify Rf1 and Rf3 candidates that can restore normal pollen production in transgenic wheat plants carrying T-CMS. We show that Rf1 and Rf3 bind to the mitochondrial orf279 transcript and induce cleavage, preventing expression of the CMS trait. The identification of restorer genes in wheat is an important step towards the development of hybrid wheat varieties based on a CMS-Rf system. The characterisation of their mode of action brings insights into the molecular basis of CMS and fertility restoration in plants.

A Lineage-Specific Paralog of Oma1 Evolved into a Gene Family from Which a Suppressor of Male Sterility-Inducing Mitochondria Emerged in Plants

Cytoplasmic male sterility (MS) in plants is caused by MS-inducing mitochondria, which have emerged frequently during plant evolution. Nuclear restorer-of-fertility (Rf)genes can suppress their cognate MS-inducing mitochondria. Whereas many Rfs encode a class of RNA-binding protein, the sugar beet (Caryophyllales) Rf encodes a protein resembling Oma1, which is involved in the quality control of mitochondria. In this study, we investigated the molecular evolution of Oma1 homologs in plants. We analyzed 37 plant genomes and concluded that a single copy is the ancestral state in Caryophyllales. Among the sugar beet Oma1 homologs, the orthologous copy is located in a syntenic region that is preserved in Arabidopsis thaliana. The sugar beet Rf is a complex locus consisting of a small Oma1 homolog family (RF-Oma1 family) unique to sugar beet. The gene arrangement in the vicinity of the locus is seen in some but not all Caryophyllalean plants and is absent from Ar. thaliana. This suggests a segmental duplication rather than a whole-genome duplication as the mechanism of RF-Oma1 evolution. Of thirty-seven positively selected codons in RF-Oma1, twenty-six of these sites are located in predicted transmembrane helices. Phylogenetic network analysis indicated that homologous recombination among the RF-Oma1 members played an important role to generate protein activity related to suppression. Together, our data illustrate how an evolutionarily young Rf has emerged from a lineage-specific paralog. Interestingly, several evolutionary features are shared with the RNA-binding protein type Rfs. Hence, the evolution of the sugar beet Rf is representative of Rf evolution in general.

Genome-wide identification of the restorer-of-fertility-like (RFL) gene family in Brassica napus and expression analysis in Shaan2A cytoplasmic male sterility

Background

Cytoplasmic male sterility (CMS) is very important in hybrid breeding. The restorer-of-fertility (Rf) nuclear genes rescue the sterile phenotype. Most of the Rf genes encode pentatricopeptide repeat (PPR) proteins.

Results

We investigated the restorer-of-fertility-like (RFL) gene family in Brassica napus. A total of 53 BnRFL genes were identified. While most of the BnRFL genes were distributed on 10 of the 19 chromosomes, gene clusters were identified on chromosomes A9 and C8. The number of PPR motifs in the BnRFL proteins varied from 2 to 19, and the majority of BnRFL proteins harbored more than 10 PPR motifs. An interaction network analysis was performed to predict the interacting partners of RFL proteins. Tissue-specific expression and RNA-seq analyses between the restorer line KC01 and the sterile line Shaan2A indicated that BnRFL1, BnRFL5, BnRFL6, BnRFL8, BnRFL11, BnRFL13 and BnRFL42 located in gene clusters on chromosomes A9 and C8 were highly expressed in KC01.

Conclusions

In the present study, identification and gene expression analysis of RFL gene family in the CMS system were conducted, and seven BnRFL genes were identified as candidates for the restorer genes in Shaan2A CMS. Taken together, this method might provide new insight into the study of Rf genes in other CMS systems.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07163-z.

The molecular basis for allelic differences suggests Restorer-of-fertility 1 is a complex locus in sugar beet (Beta vulgaris L.)

BACKGROUND

Cytoplasmic male sterility (CMS) is a widely used trait for hybrid seed production in many crops. Sugar beet CMS is associated with a unique mitochondrial protein named preSATP6 that forms a 250-kDa complex. Restorer-of-fertility 1 (Rf1) is a nuclear gene that suppresses CMS and is, hence, one of the targets of sugar beet breeding. Rf1 has dominant, semi-dominant and recessive alleles, suggesting that it may be a multi-allelic locus; however, the molecular basis for differences in genetic action is obscure. Molecular cloning of Rf1 revealed a gene (orf20) whose protein products produced in transgenics can bind with preSATP6 to generate a novel 200-kDa complex. The complex is also detected in fertility-restored anthers concomitant with a decrease in the amount of the 250-kDa complex. Molecular diversity of the Rf1 locus involves organizational diversity of a gene cluster composed of orf20-like genes (RF-Oma1s). We examined the possibility that members of the clustered RF-Oma1 in this locus could be associated with fertility restoration.

RESULTS

Six yet uncharacterized RF-Oma1s from dominant and recessive alleles were examined to determine whether they could generate the 200-kDa complex. Analyses of transgenic calli revealed that three RF-Oma1s from a dominant allele could generate the 200-kDa complex, suggesting that clustered RF-Oma1s in the dominant allele can participate in fertility restoration. None of the three copies from two recessive alleles was 200-kDa generative. The absence of this ability was confirmed by analyzing mitochondrial complexes in anthers of plants having these recessive alleles. Together with our previous data, we designed a set of PCR primers specific to the 200-kDa generative RF-Oma1s. The amount of mRNA measured by this primer set inversely correlated with the amount of the 250-kDa complex in anthers and positively correlated with the strength of the Rf1 alleles.

CONCLUSIONS

Fertility restoration by sugar beet Rf1 can involve multiple RF-Oma1s clustered in the locus, implying that stacking 200-kDa generative copies in the locus strengthens the efficacy, whereas the absence of 200-kDa generative copies in the locus makes the allele recessive irrespective of the copy number. We propose that sugar beet Rf1 is a complex locus.

What Does the Molecular Genetics of Different Types of Restorer-of-Fertility Genes Imply?

Cytoplasmic male sterility (CMS) is a widely used trait for hybrid seed production. Although male sterility is caused by S cytoplasm (male-sterility inducing mitochondria), the action of S cytoplasm is suppressed by restorer-of-fertility (Rf), a nuclear gene. Hence, the genetics of Rf has attained particular interest among plant breeders. The genetic model posits Rf diversity in which an Rf specifically suppresses the cognate S cytoplasm. Molecular analysis of Rf loci in plants has identified various genes; however, pentatricopeptide repeat (PPR) protein (a specific type of RNA-binding protein) is so prominent as the Rf-gene product that Rfs have been categorized into two classes, PPR and non-PPR. In contrast, several shared features between PPR- and some non-PPR Rfs are apparent, suggesting the possibility of another grouping. Our present focus is to group Rfs by molecular genetic classes other than the presence of PPRs. We propose three categories that define partially overlapping groups of Rfs: association with post-transcriptional regulation of mitochondrial gene expression, resistance gene-like copy number variation at the locus, and lack of a direct link to S-orf (a mitochondrial ORF associated with CMS). These groups appear to reflect their own evolutionary background and their mechanism of conferring S cytoplasm specificity.

Fine mapping of restorer-of-fertility gene based on high-density genetic mapping and collinearity analysis in pepper (Capsicum annuum L.)

The pepper restorer-of-fertility (CaRf) gene was fine mapped based on conjoint analysis of recombinants and collinearity between the two pepper reference genomes. Capana06g003028, encoding an Rf-like PPR protein, was proposed as the strongest candidate for pepper CaRf based on sequence comparison and expression analysis. The cytoplasmic male sterility (CMS)/restorer-of-fertility (Rf) system not only provides an excellent model to dissect genetic interactions between the mitochondria and nucleus but also plays a vital role in high-efficiency hybrid seed production in crops including pepper (Capsicum spp.). Although two important CMS candidate genes, orf507 and Ψatp6-2, have previously been suggested, the pepper Rf gene (CaRf) has not yet been isolated. In this study, a high-density genetic map comprising 7566 SNP markers in 1944 bins was first constructed with the array genotyping results from 317 F₂ individuals. Based on this map, the CaRf gene was preliminarily mapped to a region of 1.15 Mb in length at the end of chromosome P6. Then, by means of a conjoint analysis of recombinants and collinearity between the two pepper reference genomes, an important candidate interval with 270.10 kb in length was delimited for CaRf. Finally, Capana06g003028, which encodes an Rf-like PPR protein, was proposed as the strongest candidate for CaRf based on sequence analysis and expression characteristics in sterile and fertile plants. The high-density genetic map and fine mapping results provided here lay a foundation for the application of molecular breeding, as well as cloning and functional analysis of CaRf, in pepper.

The Expansion and Diversification of Pentatricopeptide Repeat RNA-Editing Factors in Plants

The RNA-binding pentatricopeptide repeat (PPR) family comprises hundreds to thousands of genes in most plants, but only a few dozen in algae, indicating massive gene expansions during land plant evolution. The nature and timing of these expansions has not been well defined due to the sparse sequence data available from early-diverging land plant lineages. In this study, we exploit the comprehensive OneKP datasets of over 1000 transcriptomes from diverse plants and algae toward establishing a clear picture of the evolution of this massive gene family, focusing on the proteins typically associated with RNA editing, which show the most spectacular variation in numbers and domain composition across the plant kingdom. We characterize over 2 250 000 PPR motifs in over 400 000 proteins. In lycophytes, polypod ferns, and hornworts, nearly 10% of expressed protein-coding genes encode putative PPR editing factors, whereas they are absent from algae and complex-thalloid liverworts. We show that rather than a single expansion, most land plant lineages with high numbers of editing factors have continued to generate novel sequence diversity. We identify sequence variations that imply functional differences between PPR proteins in seed plants versus non-seed plants and variations we propose to be linked to seed-plant-specific editing co-factors. Finally, using the sequence variations across the datasets, we develop a structural model of the catalytic DYW domain associated with C-to-U editing and identify a clade of unique DYW variants that are strong candidates as U-to-C RNA-editing factors, given their phylogenetic distribution and sequence characteristics.

Transcriptomics, chromosome engineering and mapping identify a restorer-of-fertility region in the CMS wheat system msH1

An original RNA-seq mapping strategy, validated with chromosome engineering and physical mapping, identifies candidate genes for fertility restoration in the 6H^chS chromosome of Hordeum chilense in the wheat msH1 system. Cytoplasmic male sterility (CMS) is a valuable trait for hybrid seed production. The msH1 CMS system in common wheat results from the incompatibility between the nuclear genome of wheat and the cytoplasm of the wild barley Hordeum chilense. This work aims to identify H. chilense candidate genes for fertility restoration in the msH1 system with a multidisciplinary strategy based on chromosome engineering, differential expression analysis and genome mapping. Alloplasmic isogenic wheat lines differing for fertility, associated with the presence of an acrocentric chromosome H^chac resulting from the rearrangement of the short arms of H. chilense chromosomes 1H^ch and 6H^ch, were used for transcriptome sequencing. Two novel RNA-seq mapping approaches were designed and compared to identify differentially expressed genes of H. chilense associated with male fertility restoration. Minichromosomes (H^chmi), new smaller reorganizations of the H^chac also restoring fertility, were obtained and used to validate the candidate genes. This strategy was successful identifying a putative restorer-of-fertility region on 6H^chS, with six candidate genes, including the ortholog of the barley restorer gene Rfm1. Additionally, transcriptomics gave preliminary insights on sterility and restoration networks showing the importance of energy supply, stress, protein metabolism and RNA processing.

Nucleotide sequence polymorphism in the RFL-PPR genes of potato

Cytoplasmic male sterility (CMS) is widely used for hybrid seed production in cultivated Solanaceae species. However, there is very limited information about CMS-Rf genetic systems in potato (Solanum tuberosum). Studying the CMS-Rf systems in potato is both of theoretical and practical significance due to the emergence of a new revolutionary strategy of reinventing potato as adiploid inbred line-based crop to develop F₁ hybrid seed potato breeding (Lindhout et al. 2011; Jansky et al. 2016). To search for potato Rf gene candidates, the comparative genetic approach was applied. Based on similarity to petunia Rf-PPR592 gene, 38 fragments were identified in five loci of the whole-genome nucleotide sequence of the accession DM 1-3 516 R44 S. tuberosum Phureja group (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The putative encoded mitochondrial proteins have 589-597 amino acid residues, similarto RF-PPR proteins of petunia and chili pepper and contain 14 or 15 PPR motifs. Primers have been developed flanking the most variable 782-865 bp regions of the selected loci, and polymorphism of the cloned fragments has been investigated in a subset of nine potato genotypes. The amplified fragments included seven or eight PPR motifs and lacked introns. The SNP frequencies ranged from 7.0 to 19.8% depending on the locus, while the ratio of nonsynonymous to synonymous substitutions varied between 0.9 and 2.1. Positions 1, 3 and 6 were the most variable in the studied PPR motifs. Our results demonstrated that the analysed sequences belong to the RFL-PPR gene subfamily and may be considered as Rf gene candidates in potato.