Explosive invasion of plant mitochondria by a group I intron

Analysis of the Aging-Related AP2/ERF Transcription Factor Gene Family in Osmanthus fragrans

Ethylene-Responsive Factor (ERF) is a key element found in the middle and lower reaches of the ethylene signal transduction pathway. It is widely distributed in plants and plays important roles in plant growth and development, hormone signal transduction, and various stress processes. Although there is research on AP/ERF family members, research on AP2/ERF in Osmanthus fragrans is lacking. Thus, in this work, AP2/ERF in O. fragrans was extensively and comprehensively analyzed. A total of 298 genes encoding OfAP2/ERF proteins with complete AP2/ERF domains were identified. Based on the number of AP2/ERF domains and the similarity among amino acid sequences between AP2/ERF proteins from A. thaliana and O. fragrans, the 298 putative OfAP2/ERF proteins were divided into four different families, including AP2 (45), ERF (247), RAV (5), and SOLOIST (1). In addition, the exon-intron structure characteristics of these putative OfAP2/ERF genes and the conserved protein motifs of their encoded OfAP2/ERF proteins were analyzed, and the results were found to be consistent with those of the population classification. A tissue-specific analysis showed the spatiotemporal expression of OfAP2/ERF in the stems and leaves of O. fragrans at different developmental stages. Specifically, 21 genes were not expressed in any tissue, while high levels of expression were found for 25 OfAP2/ERF genes in several tissues, 60 genes in the roots, 34 genes in the stems, 37 genes in young leaves, 34 genes in old leaves, 32 genes in the early flowering stage, 18 genes in the full flowering stage, and 37 genes in the late flowering stage. Quantitative RT-PCR experiments showed that OfERF110a and OfERF110b had the highest expression levels at the full-bloom stage (S4), and this gradually decreased with the senescence of petals. The expression of OfERF119c decreased first and then increased, while the expression levels of OfERF4c and OfERF5a increased constantly. This indicated that these genes may play roles in flower senescence and the ethylene response. In the subsequent subcellular localization experiments, we found that ERF1-4 was localized in the nucleus, indicating that it was expressed in the nucleus. In yeast self-activation experiments, we found that OfERF112, OfERF228, and OfERF23 had self-activation activity. Overall, these results suggest that OfERFs may have the function of regulating petal senescence in O. fragrans.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

Evolution of cox2 introns in angiosperm mitochondria and efficient splicing of an elongated cox2i691 intron

In angiosperms, the mitochondrial cox2 gene harbors up to two introns, commonly referred to as cox2i373 and cox2i691. We studied the cox2 from 222 fully-sequenced mitogenomes from 30 angiosperm orders and analyzed the evolution of their introns. Unlike cox2i373, cox2i691 shows a distribution among plants that is shaped by frequent intron loss events driven by localized retroprocessing. In addition, cox2i691 exhibits sporadic elongations, frequently in domain IV of introns. Such elongations are poorly related to repeat content and two of them showed the presence of LINE transposons, suggesting that increasing intron size is very likely due to nuclear intracelular DNA transfer followed by incorporation into the mitochondrial DNA. Surprisingly, we found that cox2i691 is erroneously annotated as absent in 30 mitogenomes deposited in public databases. Although each of the cox2 introns is ∼1.5 kb in length, a cox2i691 of 4.2 kb has been reported in Acacia ligulata (Fabaceae). It is still unclear whether its unusual length is due to a trans-splicing arrangement or the loss of functionality of the interrupted cox2. Through analyzing short-read RNA sequencing of Acacia with a multi-step computational strategy, we found that the Acacia cox2 is functional and its long intron is spliced in cis in a very efficient manner despite its length.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

Rickettsial DNA and a trans-splicing rRNA group I intron in the unorthodox mitogenome of the fern Haplopteris ensiformis

Plant mitochondrial genomes can be complex owing to highly recombinant structures, lack of gene syntenies, heavy RNA editing and invasion of chloroplast, nuclear or even foreign DNA by horizontal gene transfer (HGT). Leptosporangiate ferns remained the last major plant clade without an assembled mitogenome, likely owing to a demanding combination of the above. We here present both organelle genomes now for Haplopteris ensiformis. More than 1,400 events of C-to-U RNA editing and over 500 events of reverse U-to-C edits affect its organelle transcriptomes. The Haplopteris mtDNA is gene-rich, lacking only the ccm gene suite present in ancestral land plant mitogenomes, but is highly unorthodox, indicating extraordinary recombinogenic activity. Although eleven group II introns known in disrupted trans-splicing states in seed plants exist in conventional cis-arrangements, a particularly complex structure is found for the mitochondrial rrnL gene, which is split into two parts needing reassembly on RNA level by a trans-splicing group I intron. Aside from ca. 80 chloroplast DNA inserts that complicated the mitogenome assembly, the Haplopteris mtDNA features as an idiosyncrasy 30 variably degenerated protein coding regions from Rickettiales bacteria indicative of heavy bacterial HGT on top of tRNA genes of chlamydial origin.

Super-Mitobarcoding in Plant Species Identification? It Can Work! The Case of Leafy Liverworts Belonging to the Genus Calypogeia

Molecular identification of species is especially important where traditional taxonomic methods fail. The genus Calypogeia belongs to one of the tricky taxons. The simple morphology of these species and a tendency towards environmental plasticity make them complicated in identification. The finding of the universal single-locus DNA barcode in plants seems to be 'the Holy Grail'; therefore, researchers are increasingly looking for multiloci DNA barcodes or super-barcoding. Since the mitochondrial genome has low sequence variation in plants, species delimitation is usually based on the chloroplast genome. Unexpectedly, our research shows that super-mitobarcoding can also work! However, our outcomes showed that a single method of molecular species delimitation should be avoided. Moreover, it is recommended to interpret the results of molecular species delimitation alongside other types of evidence, such as ecology, population genetics or comparative morphology. Here, we also presented genetic data supporting the view that C. suecica is not a homogeneous species.

Generation of a mitochondrial protein compendium in Dictyostelium discoideum

The social ameba Dictyostelium discoideum has emerged as a powerful model to study mitochondrial genetics and bioenergetics. However, a comprehensive inventory of mitochondrial proteins that is critical to understanding mitochondrial processes has yet to be curated. Here, we utilized high-throughput multiplexed protein quantitation and homology analyses to generate a high-confidence mitochondrial protein compendium consisting of 936 proteins. Our proteomic approach, which utilizes mass spectrometry in combination with mathematical modeling, was validated through mitochondrial targeting sequence prediction and live-cell imaging. Our final compendium consists of 936 proteins. Nearly, a third of D. discoideum mitochondrial proteins do not have homologs in humans, budding yeasts, or an ancestral alphaproteobacteria. Additionally, we leverage our compendium to highlight the complexity of metabolic reprogramming during starvation-induced development. Our compendium lays a foundation to investigate mitochondrial processes that are unique in ameba and to understand the functions of conserved mitochondrial proteins in D. discoideum.

Shifts from cis-to trans-splicing of five mitochondrial introns in Tolypanthus maclurei

Shifts from cis-to trans-splicing of mitochondrial introns tend to correlate with relative genome rearrangement rates during vascular plant evolution, as is particularly apparent in some lineages of gymnosperms. However, although many angiosperms have also relatively high mitogenomic rearrangement rates, very few cis-to trans-splicing shifts except for five trans-spliced introns shared in seed plants have been reported. In this study, we sequenced and characterized the mitogenome of Tolypanthus maclurei, a hemiparasitic plant from the family Loranthaceae (Santalales). The mitogenome was assembled into a circular chromosome of 256,961 bp long, relatively small compared with its relatives from Santalales. It possessed a gene content of typical angiosperm mitogenomes, including 33 protein-coding genes, three rRNA genes and ten tRNA genes. Plastid-derived DNA fragments took up 9.1% of the mitogenome. The mitogenome contained one group I intron (cox1i729) and 23 group II introns. We found shifts from cis-to trans-splicing of five additional introns in its mitogenome, of which two are specific in T. maclurei. Moreover, atp1 is a chimeric gene and phylogenetic analysis indicated that a 356 bp region near the 3′ end of atp1 of T. maclurei was acquired from Lamiales via horizontal gene transfer. Our results suggest that shifts to trans-splicing of mitochondrial introns may not be uncommon among angiosperms.

Plastid phylogenomics of Pleurothallidinae (Orchidaceae): Conservative plastomes, new variable markers, and comparative analyses of plastid, nuclear, and mitochondrial data

We present the first comparative plastome study of Pleurothallidinae with analyses of structural and molecular characteristics and identification of the ten most-variable regions to be incorporated in future phylogenetic studies. We sequenced complete plastomes of eight species in the subtribe and compared phylogenetic results of these to parallel analyses of their nuclear ribosomal DNA operon (26S, 18S, and 5.8S plus associated spacers) and partial mitochondrial genome sequences (29–38 genes and partial introns). These plastomes have the typical quadripartite structure for which gene content is similar to those of other orchids, with variation only in the composition of the ndh genes. The independent loss of ndh genes had an impact on which genes border the inverted repeats and thus the size of the small single-copy region, leading to variation in overall plastome length. Analyses of 68 coding sequences indicated the same pattern of codon usage as in other orchids, and 13 protein-coding genes under positive selection were detected. Also, we identified 62 polymorphic microsatellite loci and ten highly variable regions, for which we designed primers. Phylogenomic analyses showed that the top ten mutational hotspots represent well the phylogenetic relationships found with whole plastome sequences. However, strongly supported incongruence was observed among plastid, nuclear ribosomal DNA operon, and mitochondrial DNA trees, indicating possible occurrence of incomplete lineage sorting and/or introgressive hybridization. Despite the incongruence, the mtDNA tree retrieved some clades found in other analyses. These results, together with performance in recent studies, support a future role for mitochondrial markers in Pleurothallidinae phylogenetics.

Horizontal Gene Transfers in Plants

In plants, as in all eukaryotes, the vertical transmission of genetic information through reproduction ensures the maintenance of the integrity of species. However, many reports over the past few years have clearly shown that horizontal gene transfers, referred to as HGTs (the interspecific transmission of genetic information across reproductive barriers) are very common in nature and concern all living organisms including plants. The advent of next-generation sequencing technologies (NGS) has opened new perspectives for the study of HGTs through comparative genomic approaches. In this review, we provide an up-to-date view of our current knowledge of HGTs in plants.

Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant

Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.

SPInDel Analysis of the Non-Coding Regions of cpDNA as a More Useful Tool for the Identification of Rye (Poaceae: Secale) Species

Secale is a small but very diverse genus from the tribe Triticeae (family Poaceae), which includes annual, perennial, self-pollinating and open-pollinating, cultivated, weedy and wild species of various phenotypes. Despite its high economic importance, classification of this genus, comprising 3–8 species, is inconsistent. This has resulted in significantly reduced progress in the breeding of rye which could be enriched with functional traits derived from wild rye species. Our previous research has suggested the utility of non-coding sequences of chloroplast and mitochondrial DNA in studies on closely related species of the genus Secale. Here we applied the SPInDel (Species Identification by Insertions/Deletions) approach, which targets hypervariable genomic regions containing multiple insertions/deletions (indels) and exhibiting extensive length variability. We analysed a total of 140 and 210 non-coding sequences from cpDNA and mtDNA, respectively. The resulting data highlight regions which may represent useful molecular markers with respect to closely related species of the genus Secale, however, we found the chloroplast genome to be more informative. These molecular markers include non-coding regions of chloroplast DNA: atpB-rbcL and trnT-trnL and non-coding regions of mitochondrial DNA: nad1B-nad1C and rrn5/rrn18. Our results demonstrate the utility of the SPInDel concept for the characterisation of Secale species.

Controlling Geminiviruses before Transmission: Prospects

Whitefly (Bemisia tabaci)-transmitted Geminiviruses cause serious diseases of crop plants in tropical and sub-tropical regions. Plants, animals, and their microbial symbionts have evolved complex ways to interact with each other that impact their life cycles. Blocking virus transmission by altering the biology of vector species, such as the whitefly, can be a potential approach to manage these devastating diseases. Virus transmission by insect vectors to plant hosts often involves bacterial endosymbionts. Molecular chaperonins of bacterial endosymbionts bind with virus particles and have a key role in the transmission of Geminiviruses. Hence, devising new approaches to obstruct virus transmission by manipulating bacterial endosymbionts before infection opens new avenues for viral disease control. The exploitation of bacterial endosymbiont within the insect vector would disrupt interactions among viruses, insects, and their bacterial endosymbionts. The study of this cooperating web could potentially decrease virus transmission and possibly represent an effective solution to control viral diseases in crop plants.

The insertion of a mitochondrial selfish element into the nuclear genome and its consequences

Homing endonucleases (HE) are enzymes capable of cutting DNA at highly specific target sequences, the repair of the generated double-strand break resulting in the insertion of the HE-encoding gene ("homing" mechanism). HEs are present in all three domains of life and viruses; in eukaryotes, they are mostly found in the genomes of mitochondria and chloroplasts, as well as nuclear ribosomal RNAs. We here report the case of a HE that accidentally integrated into a telomeric region of the nuclear genome of the fungal maize pathogen Ustilago maydis. We show that the gene has a mitochondrial origin, but its original copy is absent from the U. maydis mitochondrial genome, suggesting a subsequent loss or a horizontal transfer from a different species. The telomeric HE underwent mutations in its active site and lost its original start codon. A potential other start codon was retained downstream, but we did not detect any significant transcription of the newly created open reading frame, suggesting that the inserted gene is not functional. Besides, the insertion site is located in a putative RecQ helicase gene, truncating the C-terminal domain of the protein. The truncated helicase is expressed during infection of the host, together with other homologous telomeric helicases. This unusual mutational event altered two genes: The integrated HE gene subsequently lost its homing activity, while its insertion created a truncated version of an existing gene, possibly altering its function. As the insertion is absent in other field isolates, suggesting that it is recent, the U. maydis 521 reference strain offers a snapshot of this singular mutational event.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Mitochondrial genome of the nonphotosynthetic mycoheterotrophic plant Hypopitys monotropa, its structure, gene expression and RNA editing

Heterotrophic plants—plants that have lost the ability to photosynthesize—are characterized by a number of changes at all levels of organization. Heterotrophic plants are divided into two large categories—parasitic and mycoheterotrophic (MHT). The question of to what extent such changes are similar in these two categories is still open. The plastid genomes of nonphotosynthetic plants are well characterized, and they exhibit similar patterns of reduction in the two groups. In contrast, little is known about the mitochondrial genomes of MHT plants. We report the structure of the mitochondrial genome of Hypopitys monotropa, a MHT member of Ericaceae, and the expression of its genes. In contrast to its highly reduced plastid genome, the mitochondrial genome of H. monotropa is larger than that of its photosynthetic relative Vaccinium macrocarpon, and its complete size is ~810 Kb. We observed an unusually long repeat-rich structure of the genome that suggests the existence of linear fragments. Despite this unique feature, the gene content of the H. monotropa mitogenome is typical of flowering plants. No acceleration of substitution rates is observed in mitochondrial genes, in contrast to previous observations in parasitic non-photosynthetic plants. Transcriptome sequencing revealed the trans-splicing of several genes and RNA editing in 33 of 38 genes. Notably, we did not find any traces of horizontal gene transfer from fungi, in contrast to plant parasites, which extensively integrate genetic material from their hosts.

Variation in protein gene and intron content among land plant mitogenomes

The functional content of the mitochondrial genome (mitogenome) is highly diverse across eukaryotes. Among land plants, our understanding of the variation in mitochondrial gene and intron content is improving from concerted efforts to densely sample mitogenomes from diverse land plants. Here I review the current state of knowledge regarding the diversity in content of protein genes and introns in the mitogenomes of all major land plant lineages. Mitochondrial protein gene content is largely conserved among mosses and liverworts, but it varies substantially among and within other land plant lineages due to convergent losses of genes encoding ribosomal proteins and, to a lesser extent, genes for proteins involved in cytochrome c maturation and oxidative phosphorylation. Mitochondrial intron content is fairly stable within each major land plant lineage, but highly variable among lineages, resulting from occasional gains and many convergent losses over time. Trans-splicing has evolved dozens of times in various vascular plant lineages, particularly those with relatively higher rates of mitogenomic rearrangement. Across eukaryotes, mitochondrial protein gene and intron content has been shaped massive convergent evolution.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

Mitochondrial Genome of Fagopyrum esculentum and the Genetic Diversity of Extranuclear Genomes in Buckwheat

Fagopyrum esculentum (common buckwheat) is an important agricultural non-cereal grain plant. Despite extensive genetic studies, the information on its mitochondrial genome is still lacking. Using long reads generated by single-molecule real-time technology coupled with circular consensus sequencing (CCS) protocol, we assembled the buckwheat mitochondrial genome and detected that its prevalent form consists of 10 circular chromosomes with a total length of 404 Kb. In order to confirm the presence of a multipartite structure, we developed a new targeted assembly tool capable of processing long reads. The mitogenome contains all genes typical for plant mitochondrial genomes and long inserts of plastid origin (~6.4% of the total mitogenome length). Using this new information, we characterized the genetic diversity of mitochondrial and plastid genomes in 11 buckwheat cultivars compared with the ancestral subspecies, F. esculentum ssp. ancestrale. We found it to be surprisingly low within cultivars: Only three to six variations in the mitogenome and one to two in the plastid genome. In contrast, the divergence with F. esculentum ssp. ancestrale is much higher: 220 positions differ in the mitochondrial genome and 159 in the plastid genome. The SNPs in the plastid genome are enriched in non-synonymous substitutions, in particular in the genes involved in photosynthesis: psbA, psbC, and psbH. This presumably reflects the selection for the increased photosynthesis efficiency as a part of the buckwheat breeding program.

Convergent horizontal gene transfer and cross-talk of mobile nucleic acids in parasitic plants

Horizontal gene transfer (HGT), the movement and genomic integration of DNA across species boundaries, is commonly associated with bacteria and other microorganisms, but functional HGT (fHGT) is increasingly being recognized in heterotrophic parasitic plants that obtain their nutrients and water from their host plants through direct haustorial feeding. Here, in the holoparasitic stem parasite Cuscuta, we identify 108 transcribed and probably functional HGT events in Cuscuta campestris and related species, plus 42 additional regions with host-derived transposon, pseudogene and non-coding sequences. Surprisingly, 18 Cuscuta fHGTs were acquired from the same gene families by independent HGT events in Orobanchaceae parasites, and the majority are highly expressed in the haustorial feeding structures in both lineages. Convergent retention and expression of HGT sequences suggests an adaptive role for specific additional genes in parasite biology. Between 16 and 20 of the transcribed HGT events are inferred as ancestral in Cuscuta based on transcriptome sequences from species across the phylogenetic range of the genus, implicating fHGT in the successful radiation of Cuscuta parasites. Genome sequencing of C. campestris supports transfer of genomic DNA-rather than retroprocessed RNA-as the mechanism of fHGT. Many of the C. campestris genes horizontally acquired are also frequent sources of 24-nucleotide small RNAs that are typically associated with RNA-directed DNA methylation. One HGT encoding a leucine-rich repeat protein kinase overlaps with a microRNA that has been shown to regulate host gene expression, suggesting that HGT-derived parasite small RNAs may function in the parasite-host interaction. This study enriches our understanding of HGT by describing a parasite-host system with unprecedented gene exchange that points to convergent evolution of HGT events and the functional importance of horizontally transferred coding and non-coding sequences.

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha-proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein-coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome-wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad-scale genome streamlining via loss of protein-coding genes as well as noncoding, parasitic DNA elements.

In Silico Analyses of Burial Codon Bias Among the Species of Dipterocarpaceae Through Molecular and Phylogenetic Data

Introduction:

DNA barcode, a molecular marker, is used to distinguish among the closely related species, and it can be applied across a broad range of taxa to understand ecology and evolution. MaturaseK gene (matK) and rubisco bisphosphate carboxylase/oxygenase form I gene (rbcL) of the chloroplast are highly conserved in a plant system, which are used as core barcode. This present endeavor entails the comprehensive examination of the under threat plant species based on success of discrimination on DNA barcode under selection pressure.

Result:

The family Dipterocarpaceae comprising of 15 genera is under threat due to some factors, namely, deforestation, habitat alteration, poor seed, pollen dispersal, etc. Species of this family was grouped into 6 clusters for matK and 5 clusters and 2 sub-clusters for rbcL in the phylogenetic tree by using neighbor-joining method. Cluster I to cluster VI of matK and cluster I to cluster V of rbcL genes were analyzed by various codon and substitution bias tools. Mutational pressure guided the codon bias which was favored by the avoidance of higher GC content and significant negative correlation between GC12 and GC3 (in sub-cluster I of cluster I [0.03 < P], cluster I [0.00001 < P], and cluster II [0.01 < P] of rbcL, and cluster IV [0.013 < P] of matK). After refining the results, it could be speculated that the lower null expectation values (R = 0.5 or <0.5) were less divergent from the evolutionary perspective. Apart from that, the higher null expectation values (R = >0.85) also showed the same result, which possibly could be due to the negative impact of very high and low transition rate than transversion.

Conclusion:

Through the analysis of inter-generic, inter/intra-specific variation and phylogenetic data, it was found that both selection and mutation played an important role in synonymous codon choice in these genes, but they acted inconsistently on the genes, both matK and rbcL. In vitro stable proteins of both matK and rbcL were selected through natural selection rather than mutational selection. matK gene had higher individual discrimination and barcode success compared with rbcL. These discriminatory approaches may describe the problem related to the extinction of plant species. Hence, it becomes very imperative to identify and detect the under threat plant species in advance.

Horizontally-acquired genetic elements in the mitochondrial genome of a centrohelid Marophrys sp. SRT127

Mitochondrial genomes exhibit diverse features among eukaryotes in the aspect of gene content, genome structure, and the mobile genetic elements such as introns and plasmids. Although the number of published mitochondrial genomes is increasing at tremendous speed, those of several lineages remain unexplored. Here, we sequenced the complete mitochondrial genome of a unicellular heterotrophic eukaryote, Marophrys sp. SRT127 belonging to the Centroheliozoa, as the first report on this lineage. The circular-mapped mitochondrial genome, which is 113,062 bp in length, encodes 69 genes typically found in mitochondrial genomes. In addition, the Marophrys mitochondrial genome contains 19 group I introns. Of these, 11 introns have genes for homing endonuclease (HE) and phylogenetic analyses of HEs have shown that at least five Marophrys HEs are related to those in green algal plastid genomes, suggesting intron transfer between the Marophrys mitochondrion and green algal plastids. We also discovered a putative mitochondrial plasmid in linear form. Two genes encoded in the circular-mapped mitochondrial genome were found to share significant similarities to those in the linear plasmid, suggesting that the plasmid was integrated into the mitochondrial genome. These findings expand our knowledge on the diversity and evolution of the mobile genetic elements in mitochondrial genomes.

Genetic diversity and relationship between cultivated, weedy and wild rye species as revealed by chloroplast and mitochondrial DNA non-coding regions analysis

The genus Secale is small but very diverse. Despite the high economic importance, phylogenetic relationships of rye species have not been fully determined, and they are extremely important for the process of breeding of new cultivars that can be enriched with functional traits derived from wild rye species. The study analyzed the degree of relationship of 35 accessions of the genus Secale, representing 13 most often distinguished species and subspecies, originating from various seed collections in the world, based on the analysis of non-coding regions of the chloroplast (cpDNA) and mitochondrial genome (mtDNA), widely used in phylogenetic and population plant studies, because of a higher rate of evolution than the coding regions. There was no clear genetic structure between different species and subspecies, which may indicated the introgression between these taxa. The obtained data confirmed that S. vavilovii was very similar to S. cereale, which confirmed the assumption that they might share a common ancestor. The results also confirmed the divergence of S. sylvestre from other species and subspecies of rye. Areas that may be useful molecular markers in studies on closely related species of the genus Secale were also indicated.

Efficient Identification of Pulsatilla (Ranunculaceae) Using DNA Barcodes and Micro-Morphological Characters

Pulsatilla (Ranunculaceae) comprises about 40 species, many of which have horticultural and/or medicinal importance. However, the recognition and identification of wild Pulsatilla species is difficult due to the presence of complex morphological characters. DNA barcoding is a powerful molecular tool capable of rapidly and accurately distinguishing between species. Here, we assessed the effectiveness of four commonly used DNA barcoding loci-rbcL (R), trnH-psbA ( T ), matK (M), and ITS (I)-to identify species of Pulsatilla from a comprehensive sampling group. Among the four barcoding single loci, the nuclear ITS marker showed the highest interspecific distances and the highest rate of correct identification. Among the eleven combinations, the chloroplast multi-locus R+T and R+M+T combinations were found to have the best species discrimination rate, followed by R+M. Overall, we propose that the R+M+T combination and the ITS marker on its own are, respectively, the best multi- and single-locus barcodes for discriminating among species of Pulsatilla. The phylogenetic analysis was able to distinguish species of Pulsatilla to the subgenus level, but the analysis also showed relatively low species resolution. This may be caused by incomplete lineage sorting and/or hybridization events in the evolutionary history of the genus, or by the resolution limit of the candidate barcodes. We also investigated the leaf epidermis of eight representative species using scanning electronic microscopy. The resulting micro-morphological characters were valuable for identification of related species. Using additional genome fragments, or even whole chloroplast genomes combined with micro-morphological data may permit even higher resolution of species in Pulsatilla.

Recurrent Loss, Horizontal Transfer, and the Obscure Origins of Mitochondrial Introns in Diatoms (Bacillariophyta)

We sequenced mitochondrial genomes from five diverse diatoms (Toxarium undulatum, Psammoneis japonica, Eunotia naegelii, Cylindrotheca closterium, and Nitzschia sp.), chosen to fill important phylogenetic gaps and help us characterize broadscale patterns of mitochondrial genome evolution in diatoms. Although gene content was strongly conserved, intron content varied widely across species. The vast majority of introns were of group II type and were located in the cox1 or rnl genes. Although recurrent intron loss appears to be the principal underlying cause of the sporadic distributions of mitochondrial introns across diatoms, phylogenetic analyses showed that intron distributions superficially consistent with a recurrent-loss model were sometimes more complicated, implicating horizontal transfer as a likely mechanism of intron acquisition as well. It was not clear, however, whether diatoms were the donors or recipients of horizontally transferred introns, highlighting a general challenge in resolving the evolutionary histories of many diatom mitochondrial introns. Although some of these histories may become clearer as more genomes are sampled, high rates of intron loss suggest that the origins of many diatom mitochondrial introns are likely to remain unclear.

Origin and evolution of chloroplast group I introns in lichen algae

The history of group I introns is characterized by repeated horizontal transfers, even among phylogenetically distant species. The symbiogenetic thalli of lichens are good candidates for the horizontal transfer of genetic material among distantly related organisms, such as fungi and green algae. The main goal of this study was to determine whether there were different trends in intron distribution and properties among Chlorophyte algae based on their phylogenetic relationships and living conditions. Therefore, we investigated the occurrence, distribution and properties of group I introns within the chloroplast LSU rDNA in 87 Chlorophyte algae including lichen and free-living Trebouxiophyceae compared to free-living non-Trebouxiophyceae species. Overall, our findings showed that there was high diversity of group I introns and homing endonucleases (HEs) between Trebouxiophyceae and non-Trebouxiophyceae Chlorophyte algae, with divergence in their distribution patterns, frequencies and properties. However, the differences between lichen Trebouxiophyceae and free-living Trebouxiophyceae were smaller. An exception was the cL2449 intron, which was closely related to ω elements in yeasts. Such introns seem to occur more frequently in lichen Trebouxiophyceae compared to free-living Trebouxiophyceae. Our data suggest that lichenization and maintenance of lichen symbiosis for millions of years of evolution may have facilitated horizontal transfers of specific introns/HEs between symbionts. The data also suggest that sequencing of more chloroplast genes harboring group I introns in diverse algal groups may help us to understand the group I intron/HE transmission process within these organisms.

Hybrids and horizontal transfer: introgression allows adaptive allele discovery

Evolution has devised countless remarkable solutions to diverse challenges. Understanding the mechanistic basis of these solutions provides insights into how biological systems can be subtly tweaked without maladaptive consequences. The knowledge gained from illuminating these mechanisms is equally important to our understanding of fundamental evolutionary mechanisms as it is to our hopes of developing truly rational plant breeding and synthetic biology. In particular, modern population genomic approaches are proving very powerful in the detection of candidate alleles for mediating consequential adaptations that can be tested functionally. Especially striking are signals gained from contexts involving genetic transfers between populations, closely related species, or indeed between kingdoms. Here we discuss two major classes of these scenarios, adaptive introgression and horizontal gene flow, illustrating discoveries made across kingdoms.

Unparalleled replacement of native mitochondrial genes by foreign homologs in a holoparasitic plant

Horizontal gene transfer (HGT) among flowering plant mitochondria occurs frequently and, in most cases, leads to nonfunctional transgenes in the recipient genome. Parasitic plants are particularly prone to this phenomenon, but their mitochondrial genomes (mtDNA) have been largely unexplored. We undertook a large-scale mitochondrial genomic study of the holoparasitic plant Lophophytum mirabile (Balanophoraceae). Comprehensive phylogenetic analyses were performed to address the frequency, origin, and impact of HGT. The sequencing of the complete mtDNA of L. mirabile revealed the unprecedented acquisition of host-derived mitochondrial genes, representing 80% of the protein-coding gene content. All but two of these foreign genes replaced the native homologs and are probably functional in energy metabolism. The genome consists of 54 circular-mapping chromosomes, 25 of which carry no intact genes. The likely functional replacement of up to 26 genes in L. mirabile represents a stunning example of the potential effect of rampant HGT on plant mitochondria. The use of host-derived genes may have a positive effect on the host-parasite relationship, but could also be the result of nonadaptive forces.

Evolution of group I introns in Porifera: new evidence for intron mobility and implications for DNA barcoding

BACKGROUND

Mitochondrial introns intermit coding regions of genes and feature characteristic secondary structures and splicing mechanisms. In metazoans, mitochondrial introns have only been detected in sponges, cnidarians, placozoans and one annelid species. Within demosponges, group I and group II introns are present in six families. Based on different insertion sites within the cox1 gene and secondary structures, four types of group I and two types of group II introns are known, which can harbor up to three encoding homing endonuclease genes (HEG) of the LAGLIDADG family (group I) and/or reverse transcriptase (group II). However, only little is known about sponge intron mobility, transmission, and origin due to the lack of a comprehensive dataset. We analyzed the largest dataset on sponge mitochondrial group I introns to date: 95 specimens, from 11 different sponge genera which provided novel insights into the evolution of group I introns.

RESULTS

For the first time group I introns were detected in four genera of the sponge family Scleritodermidae (Scleritoderma, Microscleroderma, Aciculites, Setidium). We demonstrated that group I introns in sponges aggregate in the most conserved regions of cox1. We showed that co-occurrence of two introns in cox1 is unique among metazoans, but not uncommon in sponges. However, this combination always associates an active intron with a degenerating one. Earlier hypotheses of HGT were confirmed and for the first time VGT and secondary losses of introns conclusively demonstrated.

CONCLUSION

This study validates the subclass Spirophorina (Tetractinellida) as an intron hotspot in sponges. Our analyses confirm that most sponge group I introns probably originated from fungi. DNA barcoding is discussed and the application of alternative primers suggested.

Mobile Introns Shape the Genetic Diversity of Their Host Genes

Self-splicing introns populate several highly conserved protein-coding genes in fungal and plant mitochondria. In fungi, many of these introns have retained their ability to spread to intron-free target sites, often assisted by intron-encoded endonucleases that initiate the homing process. Here, leveraging population genomic data from Saccharomyces cerevisiae, Schizosaccharomyces pombe, and Lachancea kluyveri, we expose nonrandom patterns of genetic diversity in exons that border self-splicing introns. In particular, we show that, in all three species, the density of single nucleotide polymorphisms increases as one approaches a mobile intron. Through multiple lines of evidence, we rule out relaxed purifying selection as the cause of uneven nucleotide diversity. Instead, our findings implicate intron mobility as a direct driver of host gene diversity. We discuss two mechanistic scenarios that are consistent with the data: either endonuclease activity and subsequent error-prone repair have left a mutational footprint on the insertion environment of mobile introns or nonrandom patterns of genetic diversity are caused by exonic coconversion, which occurs when introns spread to empty target sites via homologous recombination. Importantly, however, we show that exonic coconversion can only explain diversity gradients near intron-exon boundaries if the conversion template comes from outside the population. In other words, there must be pervasive and ongoing horizontal gene transfer of self-splicing introns into extant fungal populations.

Complete mitochondrial genomes from the ferns Ophioglossum californicum and Psilotum nudum are highly repetitive with the largest organellar introns

Currently, complete mitochondrial genomes (mitogenomes) are available from all major land plant lineages except ferns. Sequencing of fern mitogenomes could shed light on the major evolutionary transitions that established mitogenomic diversity among extant lineages. In this study, we generated complete mitogenomes from the adder's tongue fern (Ophioglossum californicum) and the whisk fern (Psilotum nudum). The Psilotum mitogenome (628 kb) contains a rich complement of genes and introns, some of which are the largest of any green plant organellar genome. In the Ophioglossum mitogenome (372 kb), gene and intron content is slightly reduced, including the loss of all four mitochondrial ccm genes. Transcripts of nuclear Ccm genes also were not detected, suggesting loss of the entire mitochondrial cytochrome c maturation pathway from Ophioglossum. Both fern mitogenomes are highly repetitive, yet they show extremely low levels of active recombination. Transcriptomic sequencing uncovered ˜1000 sites of C-to-U RNA editing in both species, plus a small number (< 60) of U-to-C edit sites. Overall, the first mitochondrial genomes of ferns show a mix of features shared with lycophytes and/or seed plants and several novel genomic features, enabling a robust reconstruction of the mitogenome in the common ancestor of vascular plants.

Horizontal gene transfer is more frequent with increased heterotrophy and contributes to parasite adaptation

Horizontal gene transfer (HGT) is the transfer of genetic material across species boundaries and has been a driving force in prokaryotic evolution. HGT involving eukaryotes appears to be much less frequent, and the functional implications of HGT in eukaryotes are poorly understood. We test the hypothesis that parasitic plants, because of their intimate feeding contacts with host plant tissues, are especially prone to horizontal gene acquisition. We sought evidence of HGTs in transcriptomes of three parasitic members of Orobanchaceae, a plant family containing species spanning the full spectrum of parasitic capabilities, plus the free-living Lindenbergia Following initial phylogenetic detection and an extensive validation procedure, 52 high-confidence horizontal transfer events were detected, often from lineages of known host plants and with an increasing number of HGT events in species with the greatest parasitic dependence. Analyses of intron sequences in putative donor and recipient lineages provide evidence for integration of genomic fragments far more often than retro-processed RNA sequences. Purifying selection predominates in functionally transferred sequences, with a small fraction of adaptively evolving sites. HGT-acquired genes are preferentially expressed in the haustorium-the organ of parasitic plants-and are strongly biased in predicted gene functions, suggesting that expression products of horizontally acquired genes are contributing to the unique adaptive feeding structure of parasitic plants.

Limited mitogenomic degradation in response to a parasitic lifestyle in Orobanchaceae

In parasitic plants, the reduction in plastid genome (plastome) size and content is driven predominantly by the loss of photosynthetic genes. The first completed mitochondrial genomes (mitogenomes) from parasitic mistletoes also exhibit significant degradation, but the generality of this observation for other parasitic plants is unclear. We sequenced the complete mitogenome and plastome of the hemiparasite Castilleja paramensis (Orobanchaceae) and compared them with additional holoparasitic, hemiparasitic and nonparasitic species from Orobanchaceae. Comparative mitogenomic analysis revealed minimal gene loss among the seven Orobanchaceae species, indicating the retention of typical mitochondrial function among Orobanchaceae species. Phylogenetic analysis demonstrated that the mobile cox1 intron was acquired vertically from a nonparasitic ancestor, arguing against a role for Orobanchaceae parasites in the horizontal acquisition or distribution of this intron. The C. paramensis plastome has retained nearly all genes except for the recent pseudogenization of four subunits of the NAD(P)H dehydrogenase complex, indicating a very early stage of plastome degradation. These results lend support to the notion that loss of ndh gene function is the first step of plastome degradation in the transition to a parasitic lifestyle.

Monilophyte mitochondrial rps1 genes carry a unique group II intron that likely originated from an ancient paralog in rpl2

Intron patterns in plant mitochondrial genomes differ significantly between the major land plant clades. We here report on a new, clade-specific group II intron in the rps1 gene of monilophytes (ferns). This intron, rps1i25g2, is strikingly similar to rpl2i846g2 previously identified in the mitochondrial rpl2 gene of seed plants, ferns, and the lycophyte Phlegmariurus squarrosus Although mitochondrial ribosomal protein genes are frequently subject to endosymbiotic gene transfer among plants, we could retrieve the mitochondrial rps1 gene in a taxonomically wide sampling of 44 monilophyte taxa including basal lineages such as the Ophioglossales, Psilotales, and Marattiales with the only exception being the Equisetales (horsetails). Introns rps1i25g2 and rpl2i846g2 were likewise consistently present with only two exceptions: Intron rps1i25g2 is lost in the genus Ophioglossum and intron rpl2i846g2 is lost in Equisetum bogotense Both intron sequences are moderately affected by RNA editing. The unprecedented primary and secondary structure similarity of rps1i25g2 and rpl2i846g2 suggests an ancient retrotransposition event copying rpl2i846g2 into rps1, for which we suggest a model. Our phylogenetic analysis adding the new rps1 locus to a previous data set is fully congruent with recent insights on monilophyte phylogeny and further supports a sister relationship of Gleicheniales and Hymenophyllales.

Efficient Identification of the Forest Tree Species in Aceraceae Using DNA Barcodes

Aceraceae is a large forest tree family that comprises many economically and ecologically important species. However, because interspecific and/or intraspecific morphological variations result from frequent interspecific hybridization and introgression, it is challenging for non-taxonomists to accurately recognize and identify the tree species in Aceraceae based on a traditional approach. DNA barcoding is a powerful tool that has been proposed to accurately distinguish between species. In this study, we assessed the effectiveness of three core standard markers (matK, rbcL and ITS) plus the chloroplast locus trnS-trnG as Aceraceae barcodes. A total of 231 sequences representing 85 species in this forest family were collected. Of these four barcode markers, the discrimination power was highest for the ITS (I) region (50%) and was progressively reduced in the other three chloroplast barcodes matK (M), trnS-trnG (T) and rbcL (R); the discrimination efficiency of the ITS marker was also greater than any two-locus combination of chloroplast barcodes. However, the combinations of ITS plus single or combined chloroplast barcodes could improve species resolution significantly; T+I (90.5% resolution) and R+M+T+I (90.5% resolution) differentiated the highest portion of species in Aceraceae. Our current results show that the nuclear ITS fragment represents a more promising DNA barcode marker than the maternally inherited chloroplast barcodes. The most efficient and economical method to identify tree species in Aceraceae among single or combined DNA barcodes is the combination of T+I (90.5% resolution).

Dynamic evolution of Geranium mitochondrial genomes through multiple horizontal and intracellular gene transfers

The exchange of genetic material between cellular organelles through intracellular gene transfer (IGT) or between species by horizontal gene transfer (HGT) has played an important role in plant mitochondrial genome evolution. The mitochondrial genomes of Geraniaceae display a number of unusual phenomena including highly accelerated rates of synonymous substitutions, extensive gene loss and reduction in RNA editing. Mitochondrial DNA sequences assembled for 17 species of Geranium revealed substantial reduction in gene and intron content relative to the ancestor of the Geranium lineage. Comparative analyses of nuclear transcriptome data suggest that a number of these sequences have been functionally relocated to the nucleus via IGT. Evidence for rampant HGT was detected in several Geranium species containing foreign organellar DNA from diverse eudicots, including many transfers from parasitic plants. One lineage has experienced multiple, independent HGT episodes, many of which occurred within the past 5.5 Myr. Both duplicative and recapture HGT were documented in Geranium lineages. The mitochondrial genome of Geranium brycei contains at least four independent HGT tracts that are absent in its nearest relative. Furthermore, G. brycei mitochondria carry two copies of the cox1 gene that differ in intron content, providing insight into contrasting hypotheses on cox1 intron evolution.

Miniaturized mitogenome of the parasitic plant Viscum scurruloideum is extremely divergent and dynamic and has lost all nad genes

Despite the enormous diversity among parasitic angiosperms in form and structure, life-history strategies, and plastid genomes, little is known about the diversity of their mitogenomes. We report the sequence of the wonderfully bizarre mitogenome of the hemiparasitic aerial mistletoe Viscum scurruloideum. This genome is only 66 kb in size, making it the smallest known angiosperm mitogenome by a factor of more than three and the smallest land plant mitogenome. Accompanying this size reduction is exceptional reduction of gene content. Much of this reduction arises from the unexpected loss of respiratory complex I (NADH dehydrogenase), universally present in all 300+ other angiosperms examined, where it is encoded by nine mitochondrial and many nuclear nad genes. Loss of complex I in a multicellular organism is unprecedented. We explore the potential relationship between this loss in Viscum and its parasitic lifestyle. Despite its small size, the Viscum mitogenome is unusually rich in recombinationally active repeats, possessing unparalleled levels of predicted sublimons resulting from recombination across short repeats. Many mitochondrial gene products exhibit extraordinary levels of divergence in Viscum, indicative of highly relaxed if not positive selection. In addition, all Viscum mitochondrial protein genes have experienced a dramatic acceleration in synonymous substitution rates, consistent with the hypothesis of genomic streamlining in response to a high mutation rate but completely opposite to the pattern seen for the high-rate but enormous mitogenomes of Silene. In sum, the Viscum mitogenome possesses a unique constellation of extremely unusual features, a subset of which may be related to its parasitic lifestyle.

The SLOW GROWTH3 Pentatricopeptide Repeat Protein Is Required for the Splicing of Mitochondrial NADH Dehydrogenase Subunit7 Intron 2 in Arabidopsis

Mitochondria play an important role in maintaining metabolic and energy homeostasis in the cell. In plants, impairment in mitochondrial functions usually has detrimental effects on growth and development. To study genes that are important for plant growth, we have isolated a collection of slow growth (slo) mutants in Arabidopsis (Arabidopsis thaliana). One of the slo mutants, slo3, has a significant reduction in mitochondrial complex I activity. The slo3 mutant has a four-nucleotide deletion in At3g61360 that encodes a pentatricopeptide repeat (PPR) protein. The SLO3 protein contains nine classic PPR domains belonging to the P subfamily. The small deletion in the slo3 mutant changes the reading frame and creates a premature stop codon in the first PPR domain. We demonstrated that the SLO3-GFP is localized to the mitochondrion. Further analysis of mitochondrial RNA metabolism revealed that the slo3 mutant was defective in splicing of NADH dehydrogenase subunit7 (nad7) intron 2. This specific splicing defect led to a dramatic reduction in complex I activity in the mutant as revealed by blue native gel analysis. Complementation of slo3 by 35S:SLO3 or 35S:SLO3-GFP restored the splicing of nad7 intron 2, the complex I activity, and the growth defects of the mutant. Together, these results indicate that the SLO3 PPR protein is a splicing factor of nad7 intron 2 in Arabidopsis mitochondria.

Homologous recombination and retention of a single form of most genes shape the highly chimeric mitochondrial genome of a cybrid plant

The structure and evolution of angiosperm mitochondrial genomes are driven by extremely high rates of recombination and rearrangement. An excellent experimental system for studying these events is offered by cybrid plants, in which parental mitochondria usually fuse and their genomes recombine. Little is known about the extent, nature and consequences of mitochondrial recombination in these plants. We conducted the first study in which the organellar genomes of a cybrid - between Nicotiana tabacum and Hyoscyamus niger - were sequenced and compared to those of its parents. This cybrid mitochondrial genome is highly recombinant, reflecting at least 30 crossovers and five gene conversions between its parental genomes. It is also surprisingly large (41% and 64% larger than the parental genomes), yet contains single alleles for 90% of mitochondrial genes. Recombination produced a remarkably chimeric cybrid mitochondrial genome and occurred entirely via homologous mechanisms involving the double-strand break repair and/or break-induced replication pathways. Retention of a single form of most genes could be advantageous to minimize intracellular incompatibilities and/or reflect neutral forces that preferentially eliminate duplicated regions. We discuss the relevance of these findings to the surprisingly frequent occurrence of horizontal gene - and genome - transfer in angiosperm mitochondrial DNAs.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative “parasitism genes.” Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.

Massively parallel determination and modeling of endonuclease substrate specificity

We describe the identification and characterization of novel homing endonucleases using genome database mining to identify putative target sites, followed by high throughput activity screening in a bacterial selection system. We characterized the substrate specificity and kinetics of these endonucleases by monitoring DNA cleavage events with deep sequencing. The endonuclease specificities revealed by these experiments can be partially recapitulated using 3D structure-based computational models. Analysis of these models together with genome sequence data provide insights into how alternative endonuclease specificities were generated during natural evolution.

The chloroplast genome of Hyoscyamus niger and a phylogenetic study of the tribe Hyoscyameae (Solanaceae)

The tribe Hyoscyameae (Solanaceae) is restricted to Eurasia and includes the genera Archihyoscyamus, Anisodus, Atropa, Atropanthe, Hyoscyamus, Physochlaina, Przewalskia and Scopolia. Even though the monophyly of Hyoscyameae is strongly supported, the relationships of the taxa within the tribe remain unclear. Chloroplast markers have been widely used to elucidate plant relationships at low taxonomic levels. Identification of variable chloroplast intergenic regions has been developed based on comparative genomics of chloroplast genomes, but these regions have a narrow phylogenetic utility. In this study, we present the chloroplast genome sequence of Hyoscyamus niger and make comparisons to other solanaceous plastid genomes in terms of gene order, gene and intron content, editing sites, origins of replication, repeats, and hypothetical open reading frames. We developed and sequenced three variable plastid markers from eight species to elucidate relationships within the tribe Hyoscyameae. The presence of a horizontally transferred intron in the mitochondrial cox1 gene of some species of the tribe is considered here a likely synapomorphy uniting five genera of the Hyoscyameae. Alternatively, the cox1 intron could be a homoplasious character acquired twice within the tribe. A homoplasious inversion in the intergenic plastid spacer trnC-psbM was recognized as a source of bias and removed from the data set used in the phylogenetic analyses. Almost 12 kb of plastid sequence data were not sufficient to completely resolve relationships among genera of Hyoscyameae but some clades were identified. Two alternative hypotheses of the evolution of the genera within the tribe are proposed.