GBE
Loss of a Trans-Splicing nad1 Intron from Geraniaceae and Transfer of the Maturase Gene matR to the Nucleus in Pelargonium
Felix Grewe1,2,3, Andan Zhu1,2, and Jeffrey P. Mower1,2,* 1Center for Plant Science Innovation, University of Nebraska, Lincoln, Nebraska 2Department of Agronomy and Horticulture, University of Nebraska, Lincoln, Nebraska 3Integrative Research Center, The Field Museum of Natural History, Chicago, Illinois
*Corresponding author: E-mail: [email protected]. Accepted: September 12, 2016
Abstract The mitochondrial nad1 gene of seed plants has a complex structure, including four introns in cis or trans configurations and a maturase gene (matR) hosted within the final intron. In the geranium family (Geraniaceae), however, sequencing of representative species revealed that three of the four introns, including one in a trans configuration and another that hosts matR, were lost from the nad1 gene in their common ancestor. Despite the loss of the host intron, matR has been retained as a freestanding gene in most genera of the family, indicating that this maturase has additional functions beyond the splicing of its host intron. In the common ancestor of Pelargonium, matR was transferred to the nuclear genome, where it was split into two unlinked genes that encode either its reverse transcriptase or maturase domain. Both nuclear genes are transcribed and contain predicted mitochondrial targeting signals, suggesting that they express functional proteins that are imported into mitochondria. The nuclear localization and split domain structure of matR in the Pelargonium nuclear genome offers a unique opportunity to assess the function of these two domains using transgenic approaches. Key words: Geraniaceae, intron splicing, matR maturase, nad1 gene, retroprocessing.
Introduction assumed to be RNA mediated due to the absence of introns The mitochondrial genomes of angiosperms exhibit a diverse and the conversion of most RNA editing sites to their edited array of features that contribute to increased genomic com- state in the nuclear gene copies (Nugent and Palmer 1991; plexity, including numerous genes and introns, a large amount Covello and Gray 1992; Wischmann and Schuster 1995). of intergenic and repetitive DNA, and abundant cytidine-to- Plant mitochondrial introns can be classified as either group uridine (C-to-U) RNA editing (Knoop 2012; Mower et al. I or group II introns based on their folded structure and splicing 2012). The assortment of these complex genomic features is mechanism (Michel et al. 1989; Cech et al. 1994), with nearly highly variable among individual angiosperms. For example, all angiosperm mitochondrial introns classified as the group II the 783 kb mitogenome of Liriodendron tulipifera has 64 type. Some of these introns have evolved a split arrangement, genes, 25 introns, and >750 edit sites (Richardson et al. requiring a trans-splicing event to remove the fragmented 2013), whereas the 6.7 Mb mitogenome of Silene noctiflora intron and reconnect the independently transcribed gene has only 32 genes, 18 introns, and <200 edit sites (Sloan et al. halves into a continuous, functional transcript (Malek et al. 2010b, 2012). Mitochondrial genes encode for products in- 1997; Qiu and Palmer 2004). Loss of introns, as well as RNA volved either directly or indirectly in aerobic respiration, and edit sites, from the mitochondrial genome is often assumed to the variation in gene content among species is caused primar- occur via an RNA-mediated process termed retroprocessing, ily by intracellular gene transfer from the mitochondrial to the where a spliced and edited transcript is reintegrated into the nuclear genome (Adams et al. 2002). This process is generally genome to physically and/or functionally replace the original
ß The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected]
Genome Biol. Evol. 8(10):3193–3201. doi:10.1093/gbe/evw233 Advance Access publication September 23, 2016 3193 Grewe et al. GBE
gene (Ran et al. 2010; Sloan et al. 2010b; Grewe et al. 2011). whereas most sequenced angiosperms contain >20 sites of Horizontal transfer can also contribute to both the gain and RNA editing in this gene (Rice et al. 2013) as exemplified by loss of introns in angiosperms (Vaughn et al. 1995; Sanchez- the 28 edit sites in Melianthus nad1, nearly all of these sites Puerta et al. 2008; Hepburn et al. 2012). have been converted to thymidines in the corresponding po- The angiosperm nad1 gene epitomizes the genomic com- sitions of Geraniaceae nad1 genes. The correlated loss of in- plexity of plant mitogenomes, with five exons, four introns in trons and edit sites from Geraniaceae nad1 is consistent with cis-ortrans-spliced arrangements, and abundant RNA editing. retroprocessing activity. Despite these changes, the In addition, embedded within the final nad1 intron [named Geraniaceae nad1 gene is probably functional: it is full nad1i728 based on Dombrovska and Qiu (2004) intron nota- length and free of internal stop codons, and RNA editing of tion] of nearly all angiosperms is another gene (matR)that the few remaining sites improves sequence conservation to encodes a putative intron splicing factor termed a maturase, homologous genes from non-Geraniaceae species. although recent survey sequencing has identified a few plant lineages in which matR is no longer in this position. In Viscum Establishment of matR as a album, matR was established as a freestanding gene due to Freestanding Gene in Most loss of the host gene nad1 (Petersen et al. 2015). The matR gene is also freestanding in several species of Geranium,pre- Geraniaceae Species sumably by translocation prior to the loss of the host intron In addition to the presence of a putatively functional nad1 (Park et al. 2015). In the gnetophyte Welwitschia mirabilis, gene, most Geraniaceae species contain another partial matR is no longer adjacent to any nad1 exons, but it is flanked nad1 sequence at a distinct genomic position (fig. 1). Unlike by segments of the nad1i728 intron and may still be associ- the putatively functional nad1 copies, however, which lack the ated with the intron through a double trans-splicing event nad1i728 intron and the associated matR gene, these partial (Guo et al. 2016). In some other species of Viscum (Petersen nad1 sequences include the matR gene and parts of the et al. 2015; Skippington et al. 2015) and two species within nad1i728 intron (fig. 2A). The flanking intron sequences Malpighiales (Wurdack and Davis 2009), the matR gene ap- show clear signs of degradation based on the numerous nu- pears to be missing completely from the mitochondrial cleotide substitutions, insertions, and deletions that disrupt genome. Using the extensive genomic and transcriptomic the predicted secondary structure (fig. 2B). In contrast, the data available for many species within Geraniaceae (Weng matR sequence is presumably functional because it is full et al. 2014; Park et al. 2015; Zhang et al. 2015; Blazier length, free of internal stop codons, and transcribed based et al. 2016), we assessed the status of matR and nad1 in on detectable reads in the RNAseq library. The RNAseq this family and present an evolutionary scenario to explain reads also revealed 25 positions that are edited in at least the unusual structural and functional diversity among species. one of the seven Geraniales matR sequences, of which eight positions (32, 43, 326, 1679, 1700, 1720, 1756, and 1844) Loss of Multiple Introns from the are edited in most species (fig. 3). These data indicate that Geraniaceae nad1 Gene matR functions as a freestanding gene in most Geraniaceae species. In the large majority of angiosperms, the nad1 reading frame is interrupted by four group II introns (fig. 1). The first Migration of matR into the Nuclear (nad1i394) and third (nad1i669) nad1 introns require trans splicing for removal whereas the second intron (nad1i477) is Genome of Pelargonium removed by cis splicing. The fourth intron (nad1i728), which Because of the deep Illumina sequencing performed here, the harbors the maturase gene matR, has evolved from a cis-to matR gene was easily detectable in the draft mitochondrial trans-spliced configuration several times in angiosperms due assemblies of most Geraniaceae species. In contrast, matR was to genomic rearrangement occurring upstream or down- not detected in the mitochondrial assemblies of P. x hortorum stream of matR (Qiu and Palmer 2004). Thus, nad1i728 can and P. citronellum (fig. 1). Instead, a homolog was recovered be found as an ancestrally cis-spliced intron as observed for in the assemblies of the nuclear genome and transcriptome Melianthus villosus (fig. 1) and many other angiosperms, (fig. 4A). In both Pelargonium species, these nuclear matR whereas in other species it can have a trans-spliced arrange- sequences (annotated as nmatR) were split into two distinct ment with matR located either in the 30 or 50 intron fragment genes termed nmatRT, encoding the subdomains II to IV of the (Qiu and Palmer 2004). reverse transcriptase (RT) domain, and nmatRX, encoding sub- Within Geraniaceae, however, the structure of the nad1 domains V to VII of the RT domain and the entire maturase (X) gene has experienced a unique loss of complexity (fig. 1). domain. A glutaredoxin (grx) gene was identified upstream of Although the first intron (nad1i394) has been retained, the nmatRX and a copper/zinc superoxide dismutase (CuZnSOD) other three ancestral nad1 introns (nad1i477, nad1i669, and gene was identified downstream of nmatRT, demonstrating nad1i728) are absent from all examined species. In addition, that both nmatR genes are located in the nuclear genome.
3194 Genome Biol. Evol. 8(10):3193–3201. doi:10.1093/gbe/evw233 Advance Access publication September 23, 2016 nad1 and matR Evolution in Geraniaceae GBE
Pelargonium x hortorum Pelargonium citronellum
Geranium maderense
Monsonia emarginata Geraniales Geraniaceae
Erodium moschatum
Hypseocharis bilobata
Melianthus villosus
FIG.1.—Structure of mitochondrial nad1 and matR genes of Geraniales. Mitochondrial nad1 gene sequences and fragmented copies are represented by open boxes. Homologous intergenic regions are shown by a single line. Non-homologous regions are indicated by grey arrows. Introns in cis and trans configurations are represented by full and broken ellipses, respectively; a black color indicates the presence of the maturase gene matR. Filled black and grey dots indicate the presence and absence of RNA edit sites in Geraniales, respectively.
Furthermore, both nmatR genes are transcribed, and the nu- (95%) bootstrap support. Because P. citronellum and P. x clear intron residing in the 50 UTR of nmatRX is properly hortorum span the full taxonomic diversity of Pelargonium spliced, indicating that both nmatR genes are functionally ex- (Weng et al. 2012), this phylogenetic result indicates that pressed nuclear genes. Mitochondrial targeting signals are the nmatR genes were probably derived by intracellular trans- predicted for both genes (fig. 4A), providing evidence that ferofamitochondrialmatR sequence into the nuclear their protein products are localized to mitochondria. In con- genome of the common ancestor of Pelargonium. The ex- trast, two additional copies of nmatRT (found in a tail-to-tail treme sequence divergence for the Pelargonium nmatR se- arrangement on a different scaffold) had no transcriptional quences probably results from the elevated substitution activity and may not be functional. rates known to affect Pelargonium mitochondrial genes To infer the origin of the Pelargonium nmatR genes, phy- (Parkinson et al. 2005; Moweretal.2007) and the generally logenetic analysis was used to determine their relationship to higher rates of nuclear substitution relative to typical mito- mitochondrial matR sequences from a diversity of other an- chondrial genomes (Wolfe et al. 1987; Drouin et al. 2008). giosperms (fig. 4B; supplementary fig. S1, Supplementary Similar to the mitochondrial matR sequences, the nuclear Material online). The two nmatR sequences group together nmatR genes are more conserved within the domain regions with 100% bootstrap support, and they cluster within the relative to the rest of the protein sequence (supplementary Geraniaceae clade of mitochondrial matR genes with strong fig.S2, Supplementary Material online). This pattern exists in
Genome Biol. Evol. 8(10):3193–3201. doi:10.1093/gbe/evw233 Advance Access publication September 23, 2016 3195 Grewe et al. GBE
A I II III IV V VI Melianthus exon 4 matR exon 5 Hypseocharis Monsonia Erodium Geranium
B UGCU A C EBS 1 A A G-C C-G AU C-G A GU -5 U GU C C A A C U C-G A A-U U A-U A G-C 4 AU U-A UG G-C U U G C GU IV GUG•U A-U A A C G-C C U UUCC AC G-C G UU G ••|| G matR 6 A GC C-GAGGGG U 4/4 U U U-A G U C G-C G G C AA CCU G-C G GAGA U U A-U |•| U-A AG C-G | CC AU G A-U AC CUUU G G-C CUGGGG U 10 A-U GGA | •|•• A +4 ||| GCUCUU -199/-199 A CU CCU EBS 2 U A +28/+29 CUU CU C-G A CG U• U G CC-GG G A C C-G 11 C-G A AC-G III C-G A C-G UG A I U AG A G AG CG AC-G CA U U C-G A G GC|| A G G• C • A-U A U G A C-GU C | U CG A UA A A-U A- G AA C-G G CCCG G•U A G-C UG• G- U C V A A G II U |||| AC G-C G-CU U-AC A G G GGGC G U-A G-C A-U U-A AAG U-A CAUUCA A-U U• C G-C ||||| C-GG G-C GUAAG G•U C-G G U-A G C-G C-G G-C -14 -2 G-C U-A G 78 U-A A G-C GG-C A-U• G-C U U C U A G-C C G-C ACCA G A C U A UG G Melianthus C C C C G A C C C U G C U A-U UG C UA C A Monsonia U AA A U U-A A C C U GCGGAU CAG GUGA GGUU U 6 C CGGAG GGGU Hypseocharis UGAG A A-U G •||||| ||| ||•| •|•• C ||||| |||| |||| AGGGUGU G-CA A GG UGCUUA GUC CGUU UCGG ||||||• GG CCG UC GCCUC CCCA U Erodium ACUC UCC || ||| G || CCG A CA UA G A CACG U A GGU AG A G G C UU CC GGC U ••| U Geranium U A GA U A C C G UCA A-U C-G U G u C C UU A GC-G C-G UG g GA GG G U G G G-C G-C A-UU g g U•G G• +178 a u VI G•U G-C U C G-C +4 g c U C C IBS 1 u u U C a a U a u U G UCUC UGCG u g C u C c IBS 2 u a g u
FIG.2.—Structural degradation of the nad1i728 intron but retention of matR in Geraniaceae. (A) The diagram depicts the degree of degradation of nad1i728 in selected Geraniaceae species. (B) The secondary structure of the functional cis-spliced nad1i728 intron of M. villosus (in black) was determined by a manual comparison with existing group II introns models (Michel and Ferat 1995; Qin and Pyle 1998; Toor, et al. 2001). The intron folds into a typical group II intron secondary structure, with a central core and six branched domains (I–VI). Colored intron regions represent changes in the secondary structure of the degraded introns in G. maderense (in green), E. moschatum (in magenta), H. biloba (in orange), and M. emarginata (in blue): numbers, positive numbers, and negative numbers show sequence replacements, insertions, and deletions, respectively. The position of the maturase matR gene in domain IV is indicated by a bold black arrow. Exon binding sites (EBS) within domain I of the intron sequence (capital letters) and intron binding sites (IBS) in the adjacent exon sequence (small letters) are highlighted by a thin black arrow.
3196 Genome Biol. Evol. 8(10):3193–3201. doi:10.1093/gbe/evw233 Advance Access publication September 23, 2016 nad1 and matR Evolution in Geraniaceae GBE
32 33 43 147 153 193 235 236 237 326 413 531 539 576 1344 1543 1545 1679 1700 1720 1734 1756 1758 1825 1844 Melianthus Mito E C EEC EEEEEECCEECCEEEEEEEE Hypseocharis Mito E C E CCCTE C E CACCCCCEETCE TCE Monsonia Mito E CCCCCTCCE CCCCC C C EEEC E CCE Erodium Mito EEECCCTCCE CCCCC C C EEEC E CCE Geranium Mito E C E C E CTCCE C EECC EEEEECTCCE P. x hortorum NucACACCC– – –CC– – – A C TTTTCCCCT P. citronellum Nuc T T A C C C – – – C C – – – A C TTTTCCCCT nmatRT nmatRX
FIG.3.—Nucleotide content at positions of RNA editing in matR. All positions with an edited site in at least one species are shown. Edited cytidines are marked as “E” and shaded with a black background. Thymidines are shaded grey. All other nucleotides are unshaded.
A matR (1,872 bp) in intron nad1i728 RT X Melianthus Mitogenome
1000
100
10 nmatRX
grx RT X Pelargonium Nuclear Genome 0 Mito targeting signal UTR intron (0.79, 0.65, 0.10) 1000
100
10 nmatRT RT Cu/ZnSOD 0 Mito targeting signal (0.35, 0.63, 0.79) introns 1000
100
10 500 bp RT> 0