the mitochondrial genomes of angiosperms exhibit a diverse array of features that contribute to increased genomic complexity , including numerous genes and introns , a large amount of intergenic and repetitive dna , and abundant cytidine - to - uridine ( c - to - u ) rna editing ( knoop 2012 ; mower et al .
for example , the 783 kb mitogenome of liriodendron tulipifera has 64 genes , 25 introns , and > 750 edit sites ( richardson et al .
2013 ) , whereas the 6.7 mb mitogenome of silene noctiflora has only 32 genes , 18 introns , and < 200 edit sites ( sloan et al . 2010b , 2012 ) .
mitochondrial genes encode for products involved either directly or indirectly in aerobic respiration , and the variation in gene content among species is caused primarily by intracellular gene transfer from the mitochondrial to the nuclear genome ( adams et al .
this process is generally assumed to be rna mediated due to the absence of introns and the conversion of most rna editing sites to their edited state in the nuclear gene copies ( nugent and palmer 1991 ; covello and gray 1992 ; wischmann and schuster 1995 ) .
plant mitochondrial introns can be classified as either group i or group ii introns based on their folded structure and splicing mechanism ( michel et al .
1994 ) , with nearly all angiosperm mitochondrial introns classified as the group ii type .
some of these introns have evolved a split arrangement , requiring a trans - splicing event to remove the fragmented intron and reconnect the independently transcribed gene halves into a continuous , functional transcript ( malek et al .
loss of introns , as well as rna edit sites , from the mitochondrial genome is often assumed to occur via an rna - mediated process termed retroprocessing , where a spliced and edited transcript is reintegrated into the genome to physically and/or functionally replace the original gene ( ran et al .
horizontal transfer can also contribute to both the gain and loss of introns in angiosperms ( vaughn et al .
the angiosperm nad1 gene epitomizes the genomic complexity of plant mitogenomes , with five exons , four introns in cis- or trans - spliced arrangements , and abundant rna editing .
in addition , embedded within the final nad1 intron [ named nad1i728 based on dombrovska and qiu ( 2004 ) intron notation ] of nearly all angiosperms is another gene ( matr ) that encodes a putative intron splicing factor termed a maturase , although recent survey sequencing has identified a few plant lineages in which matr is no longer in this position . in viscum album ,
matr was established as a freestanding gene due to loss of the host gene nad1 ( petersen et al .
the matr gene is also freestanding in several species of geranium , presumably by translocation prior to the loss of the host intron ( park et al .
matr is no longer adjacent to any nad1 exons , but it is flanked by segments of the nad1i728 intron and may still be associated with the intron through a double trans - splicing event ( guo et al .
2016 ) . in some other species of viscum ( petersen et al . 2015 ; skippington et al . 2015 ) and two species within malpighiales ( wurdack and davis 2009 ) , the matr gene appears to be missing completely from the mitochondrial genome . using the extensive genomic and transcriptomic data available for many species within geraniaceae ( weng et al .
2014 ; park et al . 2015 ; zhang et al . 2015 ; blazier et al .
2016 ) , we assessed the status of matr and nad1 in this family and present an evolutionary scenario to explain the unusual structural and functional diversity among species .
in the large majority of angiosperms , the nad1 reading frame is interrupted by four group ii introns ( fig .
1 ) . the first ( nad1i394 ) and third ( nad1i669 ) nad1 introns require trans splicing for removal whereas the second intron ( nad1i477 ) is removed by cis splicing .
the fourth intron ( nad1i728 ) , which harbors the maturase gene matr , has evolved from a cis- to trans - spliced configuration several times in angiosperms due to genomic rearrangement occurring upstream or downstream of matr ( qiu and palmer 2004 ) .
thus , nad1i728 can be found as an ancestrally cis - spliced intron as observed for melianthus villosus ( fig
. 1 ) and many other angiosperms , whereas in other species it can have a trans - spliced arrangement with matr located either in the 3 or 5 intron fragment ( qiu and palmer 2004 ) . 

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 .
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 . within geraniaceae ,
however , the structure of the nad1 gene has experienced a unique loss of complexity ( fig .
1 ) . although the first intron ( nad1i394 ) has been retained , the other three ancestral nad1 introns ( nad1i477 , nad1i669 , and nad1i728 ) are absent from all examined species . in addition , whereas most sequenced angiosperms contain > 20 sites of rna editing in this gene ( rice et al .
2013 ) as exemplified by the 28 edit sites in melianthus nad1 , nearly all of these sites have been converted to thymidines in the corresponding positions of geraniaceae nad1 genes .
the correlated loss of introns and edit sites from geraniaceae nad1 is consistent with retroprocessing activity . despite these changes ,
the geraniaceae nad1 gene is probably functional : it is full length and free of internal stop codons , and rna editing of the few remaining sites improves sequence conservation to homologous genes from non - geraniaceae species .
in addition to the presence of a putatively functional nad1 gene , most geraniaceae species contain another partial nad1 sequence at a distinct genomic position ( fig .
1 ) . unlike the putatively functional nad1 copies , however , which lack the nad1i728 intron and the associated matr gene , these partial nad1 sequences include the matr gene and parts of the nad1i728 intron ( fig .
the flanking intron sequences show clear signs of degradation based on the numerous nucleotide substitutions , insertions , and deletions that disrupt the predicted secondary structure ( fig . 2b ) .
in contrast , the matr sequence is presumably functional because it is full length , free of internal stop codons , and transcribed based on detectable reads in the rnaseq library .
the rnaseq reads also revealed 25 positions that are edited in at least one of the seven geraniales matr sequences , of which eight positions ( 32 , 43 , 326 , 1679 , 1700 , 1720 , 1756 , and 1844 ) are edited in most species ( fig .
3 ) . these data indicate that matr functions as a freestanding gene in most geraniaceae species . 

( 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 .
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 . 

( 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 .
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 .
because of the deep illumina sequencing performed here , the matr gene was easily detectable in the draft mitochondrial assemblies of most geraniaceae species .
in contrast , matr was not detected in the mitochondrial assemblies of p. x hortorum and p. citronellum ( fig .
, a homolog was recovered in the assemblies of the nuclear genome and transcriptome ( fig .
4a ) . in both pelargonium species , these nuclear matr sequences ( annotated as nmatr )
were split into two distinct genes termed nmatrt , encoding the subdomains ii to iv of the reverse transcriptase ( rt ) domain , and nmatrx , encoding subdomains v to vii of the rt domain and the entire maturase ( x ) domain .
a glutaredoxin ( grx ) gene was identified upstream of nmatrx and a copper / zinc superoxide dismutase ( cuznsod ) gene was identified downstream of nmatrt , demonstrating that both nmatr genes are located in the nuclear genome .
furthermore , both nmatr genes are transcribed , and the nuclear intron residing in the 5 utr of nmatrx is properly spliced , indicating that both nmatr genes are functionally expressed nuclear genes .
in contrast , two additional copies of nmatrt ( found in a tail - to - tail arrangement on a different scaffold ) had no transcriptional activity and may not be functional . 

( a ) the mitochondrial matr gene in the nad1 intron nad1i728 in m. villosus is shown at top , and homologous matr sequences in contigs from the pelargonium nuclear genome are shown below .
grey shading shows conserved reverse transcriptase ( rt ) and maturase ( x ) domains .
the thick blue bar following the first potential start codon of the nmatr sequences defines the location of the acquired mitochondrial targeting sequence ; targeting prediction scores for mitoprot , targetp , and predotar are shown in blue in parentheses below the target sequence .
( b ) excerpt of a phylogenetic tree of the mitochondrial matr gene of angiosperms ( thin black branches ) and the nuclear nmatr paralogs of pelargonium ( thick red branches ) .
( a ) the mitochondrial matr gene in the nad1 intron nad1i728 in m. villosus is shown at top , and homologous matr sequences in contigs from the pelargonium nuclear genome are shown below . start and stop codons are indicated by the arrow and bar , respectively .
grey shading shows conserved reverse transcriptase ( rt ) and maturase ( x ) domains .
the thick blue bar following the first potential start codon of the nmatr sequences defines the location of the acquired mitochondrial targeting sequence ; targeting prediction scores for mitoprot , targetp , and predotar are shown in blue in parentheses below the target sequence .
( b ) excerpt of a phylogenetic tree of the mitochondrial matr gene of angiosperms ( thin black branches ) and the nuclear nmatr paralogs of pelargonium ( thick red branches ) .
the full phylogenetic tree is presented in supplementary figure s1 , supplementary material online . to infer the origin of the pelargonium nmatr genes , phylogenetic analysis was used to determine their relationship to mitochondrial matr sequences from a diversity of other angiosperms ( fig .
the two nmatr sequences group together with 100% bootstrap support , and they cluster within the geraniaceae clade of mitochondrial matr genes with strong ( 95% ) bootstrap support . because p. citronellum and p. x hortorum span the full taxonomic diversity of pelargonium ( weng et al . 2012 ) , this phylogenetic result indicates that the nmatr genes were probably derived by intracellular transfer of a mitochondrial matr sequence into the nuclear genome of the common ancestor of pelargonium . the extreme sequence divergence for the pelargonium nmatr sequences probably results from the elevated substitution rates known to affect pelargonium mitochondrial genes ( parkinson et al .
2005 ; mower et al . 2007 ) and the generally higher rates of nuclear substitution relative to typical mitochondrial genomes ( wolfe et al . 1987 ; drouin et al . 2008 ) .
similar to the mitochondrial matr sequences , the nuclear nmatr genes are more conserved within the domain regions relative to the rest of the protein sequence ( supplementary fig.s2 , supplementary material online ) .
this pattern exists in both halves of the split gene , nmatrt and nmatrx , and indicates that the gene remained functional after the transfer into the nuclear genome .
many of the edited positions in the geraniaceae mitochondrial matr transcripts have been converted to a thymidine in the pelargonium nmatrt and nmatrx genes ( fig .
3 ) , suggesting that transfer to the nucleus involved a partially edited rna intermediate .
because c - to - u rna editing is not known to exist for nuclear genes we predict that the lack of conversion of some edit sites does not inhibit activity of the two nmatr gene products .
the assemblies of the nad1 gene from representative geraniaceae species revealed an unusual pattern of evolution , including ( 1 ) the first demonstrated loss of a trans - configured intron ( nad1i669 ) from a functional gene , ( 2 ) the loss of the fourth nad1 intron ( nad1i728 ) that hosts the matr gene , ( 3 ) the establishment of matr as a freestanding gene within the mitogenome , and ( 4 ) the first reported transfer of matr to the nuclear genome .
all five nad introns in trans configurations have been lost from several mistletoe mitogenomes ( petersen et al .
2015 ) , but in these cases the intron losses were due to loss of the nad genes themselves .
pcr results suggested that nad1i728 was lost from two species in malpighiales ( wurdack and davis 2009 ) , but this intron is not trans - spliced in malpighiales ( qiu and palmer 2004 ; rivarola et al .
, the absence of nad1i669 in geraniaceae represents the first reported loss of a trans - configured intron from any genome .
thus , although trans splicing appears to be a strong barrier to intron loss in eukaryotes , this result demonstrates that it is not an absolute barrier .
multiple genomic changes are necessary to explain the loss of introns from nad1 and the origin of a freestanding matr gene in geraniaceae , and one possible evolutionary scenario is shown in figure 5 . the shared loss of most , but not all , introns and edit sites in geraniaceae could be attributed to incomplete retroprocessing , in which a segment of the nad1 gene was gene converted by a transcript that was only partially edited and spliced .
most other examples of retroprocessing also involve the removal of a subset of introns and edit sites ( ran et al .
a second coexisting locus , lacking the nad1i728 intron and almost all edit sites , may have been created by a retroduplication event , involving the genomic integration of a partially spliced and edited transcript into a new genomic locus , eliminating the nad1i728 intron and additional edit sites from one of the loci .
we then propose that the nad1 locus lacking the nad1i728 intron became the functional component of the gene . at the other locus ,
once the matr gene was established as a freestanding gene , the splicing of nad1i728 became unnecessary , and the intron sequence underwent lineage - specific degradation from its 5 and 3 ends .
finally , ongoing gene conversion , which has been documented for homologous sequences in plant mitogenomes ( hao and palmer 2009 ; hao et al .
2010a ) , could explain the 100% sequence identity of the nad1 exon sequences shared between the two loci .
more generally , the newly translocated matr gene described here provides a rare glimpse into the mechanism of gene translocation . 

the effects of retroprocessing , retroduplication , gene conversion , and lineage - specific gene degradation , as described in the text , are shown .
the effects of retroprocessing , retroduplication , gene conversion , and lineage - specific gene degradation , as described in the text , are shown .
the mitochondrial - to - nuclear transfer of maturase genes has occurred multiple times during land plant evolution ( mohr and lambowitz 2003 ; guo and mower 2013 ) , but until now , there was no evidence for a nuclear transfer of matr , the sole maturase remaining in seed plant mitogenomes . in viscum
scurruloideum , matr and nearly all introns are absent from the mitogenome , but surveys of genome sequencing data did not identify a homolog in the nucleus suggesting that matr and all of its splicing targets have been lost from this species ( skippington et al .
, pcr results suggest that matr was lost from the mitogenome , but nothing is known about a potential nuclear transfer ( wurdack and davis 2009 ) .
the transfer of matr to the nucleus in pelargonium was likely facilitated by several of the unusual characteristics of this genus .
first , the heavy reduction in rna editing in this genus ( parkinson et al . 2005 ) raises the possibility of either a dna - mediated or an rna - mediated transfer event ( henze and martin 2001 ) , whereas the presence of editing sites in mitochondrial genes of most other species precludes functional transfers occurring via dna .
a similar case has been suggested for the mitochondrial rpl5 gene in grasses , in which an ancestral retroprocessing event may have facilitated multiple instances of functional dna - mediated gene transfer into the nuclear genome ( ong 2006 ) .
second , if a transfer event occurs , the very high substitution rate in pelargonium mitogenomes ( parkinson et al . 2005 ) would quickly degrade the mitochondrial copy , thus making the transferred nuclear copy essential . from a functional perspective
, the retention of a matr gene in geraniaceae after loss of its host nad1i728 intron implies that matr has functions in the mitochondrion beyond the splicing of its host intron .
one possibility is that matr facilitates the removal of other mitochondrial group ii introns , and thus is still required for one or more introns remaining in geraniaceae ( supplementary table s1 , supplementary material online ) .
a recent review of mitochondrial splicing factors suggested an association of matr to other introns ( brown et al .
2014 ) consistent with additional splicing assignments , and the plastid maturase matk was shown in vivo to have binding activity to several chloroplast introns ( zoschke et al .
it is also possible that matr performs other essential transcriptional functions , such as rna processing or stabilization .
genetic manipulation of the plant mitochondrial genome is not yet possible in plants , precluding any direct assessment of matr function .
in contrast , because nuclear transformation is possible in many plants , including pelargonium ( colling et al .
2012 ) , the nuclear location of matr in pelargonium raises the possibility of genetic and transcriptional manipulation of this gene .
furthermore , the split nature of these nmatr genes enables the independent assessment of the functions of the rt and x domains .
source of plant materials and procedures for nucleic acid extraction and illumina sequencing were described previously ( weng et al .
draft mitochondrial genomes were assembled with velvet 1.1.06 ( zerbino and birney 2008 ) using a combination of kmer ( 5191 ) and expected coverage ( 20500 ) values as previously described ( grewe et al .
the assembly with the longest average length of identifiable mitochondrial sequences ( based on blastn searches with known mitochondrial protein - coding genes from related species as query sequences ) was selected for further processing ( supplementary table s2 , supplementary material online ) .
mitochondrial gene and intron content was identified by inspection of the blastn search results ( supplementary table s1 , supplementary material online ) .
for some species , the presence of two nad1 gene sequences resulted in a failure to assemble complete loci .
these regions were manually corrected by aligning and inspecting individual sequence reads that cover the respective regions .
the assembled matr and nad1 sequences from this study were deposited in genbank ( accession numbers kx824067kx824107 ) .
a draft nuclear genome for p. citronellum was assembled using velvet with kmer ( 41 ) and expected coverage ( 20 ) values that were reduced ( relative to the mitochondrial assembly parameters ) in order to preferentially assemble the lower - depth nuclear genome .
a blastn search identified matr homologs on three contigs , with sizes of 19,729 bp ( containing nmatrx ) , 5,204 bp ( containing nmatrt ) , and 16,625 bp ( containing pseudo - nmatrt ) . in plant cells , the mitogenome is typically present at a substantially higher number of copies compared with the nuclear genome ( lamppa and bendich 1984 ; draper and hays 2000 ) .
thus , the much lower depth of coverage for these three contigs ( relative to the identified mitochondrial contigs ) provides reliable evidence that these contigs are nuclear rather than mitochondrial .
additional nuclear genes were detected in these three contigs by querying them against the non - redundant protein database using blastx , and some of the genes contain nuclear spliceosomal introns with the canonical gt and ag bases at their 5 and 3 splice sites , providing additional support that these contigs are from the nuclear genome .
mitochondrial targeting signals were predicted with the programs mitoprot ii ( claros and vincens 1996 ) , targetp 1.1 ( emanuelsson et al . 2000 ) , and predotar v1.03 ( small et al . 2004 ) . the three contig sequences were deposited in genbank ( accession numbers kx824108kx824110 ) .
rnaseq reads were mapped onto the matr and nad1 sequences and the nmatr nuclear contigs using bowtie2 2.2.6 ( langmead and salzberg 2012 ) .
mapping results were manually inspected to identify exon  intron splicing junctions . for the mitochondrial genes , edit sites
were identified by searching for c - t mismatches detectable in at least 10% of the mapped rna reads .
phylogenetic analysis of matr included sequences from 39 angiosperms ( supplementary table s3 , supplementary material online ) .
2011 ) and trimmed using gblocks ( castresana 2000 ) in codon mode with relaxed parameters ( b2 = 20 ; b4 = 5 ; b5 = 20 ) .
maximum likelihood trees were constructed with the gtr + g substitution model in raxml ( stamatakis 2006 ) .
supplementary figures s1s2 and tables s1s3 are available at genome biology and evolution online ( http://www.gbe.oxfordjournals.org/ ) .
 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 . 
Introduction
Loss of Multiple Introns from the Geraniaceae 
Establishment of 
Migration of 
Discussion
Material and Methods
Supplementary Material