Mitochondrial genome evolution in parasitic plants.

BACKGROUND: Parasitic plants rely on their host to cover their nutritional requirements either for their entire life or a smaller part of it. Depending on the level of parasitism, a proportional reduction on the plastid genome has been found. However, knowledge on gene loss and evolution of the mitogenome of parasitic plants is only available for four hemiparasitic Viscum species (Viscaceae), which lack many of the mitochondrial genes, while the remaining genes exhibit very fast molecular evolution rates. In this study, we include another genus, Phoradendron, from the Viscaceae, as well as 10 other hemiparasitic or holoparasitic taxa from across the phylogeny of the angiosperms to investigate how fast molecular evolution works on their mitogenomes, and the extent of gene loss. RESULTS: Our observations from Viscum were replicated in Phoradendron liga, whereas the remaining parasitic plants in the study have a complete set of the core mitochondrial genes and exhibit moderate or only slightly raised substitution rates compared to most autotrophic taxa, without any statistically significant difference between the different groups (autotrophs, hemiparasites and holoparasites). Additionally, further evidence is provided for the placement of Balanophoraceae within the order Santalales, while the exact placement of Cynomoriaceae still remains elusive. CONCLUSIONS: We examine the mitochondrial gene content of 11 hemiparasitic and holoparasitic plants and confirm previous observations in Viscaceae. We show that the remaining parasitic plants do not have significantly higher substitution rates than autotrophic plants in their mitochondrial genes. We provide further evidence for the placement of Balanophoraceae in the Santalales.

Phylogeny of the Alismatales (Monocotyledons) and the relationship of Acorus (Acorales?).

A phylogenetic analysis of the early branching lineages of the monocotyledons is performed using data from two plastid genes (rbcL and matK), five mitochondrial genes (atp1, ccmB, cob, mttB and nad5) and morphology. The complete matrix includes 93 terminals representing Acorus, the 14 families currently recognized within Alismatales, and numerous lineages of monocotyledons and other angiosperms. Total evidence analysis results in an almost completely resolved strict consensus tree, but all data partitions, genomic as well as morphological, are incongruent. The effects of RNA editing and potentially processed paralogous sequences are explored and discussed. Despite a decrease in incongruence length differences after exclusion of edited sites, the major data partitions remain significantly incongruent. The 14 families of Alismatales are all found to be monophyletic, but Acorus is found to be included in Alismatales rather than being the sister group to all other monocotyledons. The placement is strongly supported by the mitochondrial data, atp1 in particular, but it cannot be explained as an artifact caused by patterns of editing or by sampling of processed paralogues.

Plastid phylogenomics and molecular evolution of Alismatales.

Past phylogenetic studies of the monocot order Alismatales left several higher-order relationships unresolved. We addressed these uncertainties using a nearly complete genus-level sampling of whole plastid genomes (gene sets representing 83 protein-coding and ribosomal genes) from members of the core alismatid families, Tofieldiaceae and additional taxa (Araceae and other angiosperms). Parsimony and likelihood analyses inferred generally highly congruent phylogenetic relationships within the order, and several alternative likelihood partitioning schemes had little impact on patterns of clade support. All families with multiple genera were resolved as monophyletic, and we inferred strong bootstrap support for most inter- and intrafamilial relationships. The precise placement of Tofieldiaceae in the order was not well supported. Although most analyses inferred Tofieldiaceae to be the sister-group of the rest of the order, one likelihood analysis indicated a contrasting Araceae-sister arrangement. Acorus (Acorales) was not supported as a member of the order. We also investigated the molecular evolution of plastid NADH dehydrogenase, a large enzymatic complex that may play a role in photooxidative stress responses. Ancestral-state reconstructions support four convergent losses of a functional NADH dehydrogenase complex in Alismatales, including a single loss in Tofieldiaceae.

The Global Genome Biodiversity Network (GGBN) Data Portal.

The Global Genome Biodiversity Network (GGBN) was formed in 2011 with the principal aim of making high-quality well-documented and vouchered collections that store DNA or tissue samples of biodiversity, discoverable for research through a networked community of biodiversity repositories. This is achieved through the GGBN Data Portal (http://data.ggbn.org), which links globally distributed databases and bridges the gap between biodiversity repositories, sequence databases and research results. Advances in DNA extraction techniques combined with next-generation sequencing technologies provide new tools for genome sequencing. Many ambitious genome sequencing projects with the potential to revolutionize biodiversity research consider access to adequate samples to be a major bottleneck in their workflow. This is linked not only to accelerating biodiversity loss and demands to improve conservation efforts but also to a lack of standardized methods for providing access to genomic samples. Biodiversity biobank-holding institutions urgently need to set a standard of collaboration towards excellence in collections stewardship, information access and sharing and responsible and ethical use of such collections. GGBN meets these needs by enabling and supporting accessibility and the efficient coordinated expansion of biodiversity biobanks worldwide.

Evolution of Asparagus L. (Asparagaceae): Out-of-South-Africa and multiple origins of sexual dimorphism.

In the most comprehensive study to date we explored the phylogeny and evolution of the genus Asparagus, with emphasis on the southern African species. We included 211 accessions, representing 77 (92%) of the southern African, 6 (17%) of the tropical African, 10 (56%) of the strictly European and 6 (9%) of the Eurasian species. We analyzed DNA sequences from three plastid regions (trnH-psbA, trnD-T, ndhF) and from the nuclear region phytochrome C (PHYC) with parsimony and maximum likelihood methods, and recovered a monophyletic Asparagus. The phylogeny conflicts with all previous infra-generic classifications. It has many strongly supported clades, corroborated by morphological characters, which may provide a basis for a revised taxonomy. Additionally, the phylogeny indicates that many of the current species delimitations are problematic. Using biogeographic analyses that account for phylogenetic uncertainty (S-DIVA) and take into account relative branch lengths (Lagrange) we confirm the origin of Asparagus in southern Africa, and find no evidence that the dispersal of Asparagus follow the Rand flora pattern. We find that all truly dioecious species of Asparagus share a common origin, but that sexual dimorphism has arisen independently several times.

Phylogeny of the Liliales (Monocotyledons) with special emphasis on data partition congruence and RNA editing.

A phylogenetic analysis of the monocot order Liliales was performed using sequence data from three mitochondrial (atp1, cob, nad5) and two plastid genes (rbcL, ndhF). The complete data matrix includes 46 terminals representing all 10 families currently included in Liliales. The two major partitions, mitochondrial and plastid data, were congruent, and parsimony analysis resulted in 50 equally parsimonious trees and a well resolved consensus tree confirming monophyly of all families. Mitochondrial genes are known to include RNA edited sites, and in some cases unprocessed genes are replaced by retro-processed gene copies, that is processed paralogs. To test the effects on phylogeny reconstruction of predicted edited sites and potentially unintentionally sampled processed paralogs, a number of analyses were performed using subsets of the complete data matrix. In general, predicted edited sites were more homoplasious than the other characters and increased incongruence among most data partitions. The predicted edited sites have a non-random phylogenetic signal in conflict with the signal of the non-edited sites. The potentially misleading signal was caused partially by the apparent presence of processed paralogs in Galanthus (Amaryllidaceae), part of the outgroup, but also by a deviating evolutionary pattern of predicted edited sites in Liliaceae compared with the remainder of the Liliales. Despite the problems that processed paralogs may cause, we argue that they should not a priori be excluded from phylogenetic analysis.

Genes and processed paralogs co-exist in plant mitochondria.

RNA-mediated gene duplication has been proposed to create processed paralogs in the plant mitochondrial genome. A processed paralog may retain signatures left by the maturation process of its RNA precursor, such as intron removal and no need of RNA editing. Whereas it is well documented that an RNA intermediary is involved in the transfer of mitochondrial genes to the nucleus, no direct evidence exists for insertion of processed paralogs in the mitochondria (i.e., processed and un-processed genes have never been found simultaneously in the mitochondrial genome). In this study, we sequenced a region of the mitochondrial gene nad1, and identified a number of taxa were two different copies of the region co-occur in the mitochondria. The two nad1 paralogs differed in their (a) presence or absence of a group II intron, and (b) number of edited sites. Thus, this work provides the first evidence of co-existence of processed paralogs and their precursors within the plant mitochondrial genome. In addition, mapping the presence/absence of the paralogs provides indirect evidence of RNA-mediated gene duplication as an essential process shaping the mitochondrial genome in plants.

Phylogeny of the Asparagales based on three plastid and two mitochondrial genes.

PREMISE OF THE STUDY: The Asparagales, with ca. 40% of all monocotyledons, include a host of commercially important ornamentals in families such as Orchidaceae, Alliaceae, and Iridaceae, and several important crop species in genera such as Allium, Aloe, Asparagus, Crocus, and Vanilla. Though the order is well defined, the number of recognized families, their circumscription, and relationships are somewhat controversial. METHODS: Phylogenetic analyses of Asparagales were based on parsimony and maximum likelihood using nucleotide sequence variation in three plastid genes (matK, ndhF, and rbcL) and two mitochondrial genes (atp1 and cob). Branch support was assessed using both jackknife analysis implementing strict-consensus (SC) and bootstrap analysis implementing frequency-within-replicates (FWR). The contribution of edited sites in the mitochondrial genes to topology and branch support was investigated. KEY RESULTS: The topologies recovered largely agree with previous results, though some clades remain poorly resolved (e.g., Ruscaceae). When the edited sites were included in the analysis, the plastid and mitochondrial genes were highly incongruent. However, when the edited sites were removed, the two partitions became congruent. CONCLUSIONS: Some deeper nodes in the Asparagales tree remain poorly resolved or unresolved as do the relationships of certain monogeneric families (e.g., Aphyllanthaceae, Ixioliriaceae, Doryanthaceae), whereas support for many families increases. However, the increased support is dominated by plastid data, and the potential influence of mitochondrial and biparentially inherited single or low-copy nuclear genes should be investigated.

Genes from oxidative phosphorylation complexes II-V and two dual-function subunits of complex I are transcribed in Viscum album despite absence of the entire mitochondrial holo-complex I.

Mistletoes (Viscum) and close relatives are unique among flowering plants in having a drastically altered electron transport chain. Lack of complex I genes has previously been reported for the mitochondrial genome, and here we report an almost complete absence of nuclear-encoded complex I genes in the transcriptome of Viscum album. Compared to Arabidopsis with approximately 40 nuclear complex I genes, we recover only transcripts of two dual-function genes: gamma carbonic anhydrase and L-galactono-1,4-lactone dehydrogenase. The complement of genes belonging to complexes II-V of the oxidative phosphorylation pathway appears to be in accordance with other vascular plants. Additionally, transcripts encoding alternative NAD(P)H dehydrogenases and alternative oxidase were found. Despite sequence divergence, structural modeling suggests that the encoded proteins are structurally conserved. Complex I loss is a special feature in Viscum species and relatives, as all other parasitic flowering plants investigated to date seem to have a complete OXPHOS system. Hence, Viscum offers a unique system for specifically investigating molecular consequences of complex I absence, such as the role of complex I subunits involved in secondary functions.

The genomic basis of the plant island syndrome in Darwin's giant daisies.

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the 'plant island syndrome', include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin's giant daisies.

Phylogenetic relationships in the millipede family Julidae.

A phylogenetic analysis of 40 species (22 genera) of the Palaearctic millipede family Julidae was made based on partial sequences of the mitochondrial 16S rRNA (16S) gene and the nuclear 28S rRNA (28S) gene, respectively. The two data sets (16S rDNA and 28S rDNA) were analysed individually and in combination using direct optimization as implemented in POY. The 16S rDNA and the 28S rDNA sequences vary from 410 to 449 bp and from 467 to 525 bp in length, respectively. All searches were performed under six different gap opening costs, an extension gap cost of 1, and a substitution cost of 2. Incongruence length difference values were used to select the preferred tree. The order Julida was recovered as monophyletic under all weight sets. The family Julidae was recovered as monophyletic except under one weight set where the genus Nepalmatoiulus is sister to all other Julida. Within Julidae, a clade of Paectophyllini + Calyptophyllini is sister to all others on the preferred tree but this relationship is not robust. A hitherto unrecognized clade of (South) east Asian genera (Anaulaciulus and Nepalmatoiulus) was recovered under five weight sets. Another "new" robust clade (Oncoiulini + Schizophyllini) is congruent with a hitherto unrecognized complex morphological character. Further clades recovered within the Julidae partly conflict with the accepted classification, which is only to a limited extent based on phylogenetic arguments.

When is enough, enough in phylogenetics? A case in point from Hordeum (Poaceae).

Direct optimization was used to reconstruct the phylogeny of the 26 diploid taxa included in the genus Hordeum. The total data set was composed of 16 nucleotide sequence regions from the nuclear as well as the plastid genome. The nine nuclear regions were from single-copy, protein coding genes located on six of the seven chromosome pairs in the diploid H. vulgare genome. The seven plastid regions comprise protein coding genes as well as intergenic regions. Studies of character congruence between data partitions showed no correlation between chromosomal location and congruence among the nuclear sequences and a level of congruence among the plastid sequences comparable with the level among the nuclear sequences. Combined analysis of all data resolved the phylogeny completely with most clades being robust and well supported. However, due to incongruence among data partitions some relationships are still and likely to remain ambiguously inferred. Rather than adding still more genes to the phylogenetic analyses, patterns of incongruence may be better explored by adding data from multiple specimens per taxon. For some species relationships the plastid data appear positively misleading, emphasizing the need for caution if plastid data are the only or dominant type of data used for phylogenetic reconstruction and subsequent re-classification. © The Willi Hennig Society 2011.

Are substitution rates and RNA editing correlated?

BACKGROUND: RNA editing is a post-transcriptional process that, in seed plants, involves a cytosine to uracil change in messenger RNA, causing the translated protein to differ from that predicted by the DNA sequence. RNA editing occurs extensively in plant mitochondria, but large differences in editing frequencies are found in some groups. The underlying processes responsible for the distribution of edited sites are largely unknown, but gene function, substitution rate, and gene conversion have been proposed to influence editing frequencies. RESULTS: We studied five mitochondrial genes in the monocot order Alismatales, all showing marked differences in editing frequencies among taxa. A general tendency to lose edited sites was observed in all taxa, but this tendency was particularly strong in two clades, with most of the edited sites lost in parallel in two different areas of the phylogeny. This pattern is observed in at least four of the five genes analyzed. Except in the groups that show an unusually low editing frequency, the rate of C-to-T changes in edited sites was not significantly higher that in non-edited 3rd codon positions. This may indicate that selection is not actively removing edited sites in nine of the 12 families of the core Alismatales. In all genes but ccmB, a significant correlation was found between frequency of change in edited sites and synonymous substitution rate. In general, taxa with higher substitution rates tend to have fewer edited sites, as indicated by the phylogenetically independent correlation analyses. The elimination of edited sites in groups that lack or have reduced levels of editing could be a result of gene conversion involving a cDNA copy (retroprocessing). If so, this phenomenon could be relatively common in the Alismatales, and may have affected some groups recurrently. Indirect evidence of retroprocessing without a necessary correlation with substitution rate was found mostly in families Alismataceae and Hydrocharitaceae (e.g., groups that suffered a rapid elimination of all their edited sites, without a change in substitution rate). CONCLUSIONS: The effects of substitution rate, selection, and/or gene conversion on the dynamics of edited sites in plant mitochondria remain poorly understood. Although we found an inverse correlation between substitution rate and editing frequency, this correlation is partially obscured by gene retroprocessing in lineages that have lost most of their edited sites. The presence of processed paralogs in plant mitochondria deserves further study, since most evidence of their occurrence is circumstantial.

Alignment, clade robustness and fungal phylogenetics- Crepidotaceae and sister families revisited.

The circumscription and delimitation of the agaric family Crepidotaceae and related families in the Agaricales have been analysed using nLSU sequence data and direct optimization. The influence of alignment parameters on clade robustness has been explored, and many clades are found to have a low robustness. Nevertheless, based on clades with a high or moderately high robustness to increased gap costs, a few taxonomic changes are suggested: Crepidotaceae has been amended to include Crepidotus, Episphaeria, Simocybe, Pleuroflammula and Inocybe, the latter genus adding the homoplasious character mycorrhization to the family, Pellidiscus is included in Crepidotus, a position supported by micromorphology, two other genera formerly included in the family, Chromocyphella and Phaeosolenia, are excluded, and a new family, Chromocyphellaceae fam. nov., is erected for them based on fruit body shape, microcharacters, and sequences. Monophyly of Tubariaceae as recently circumscribed was not recovered, but tentatively we accept the genera Tubaria and Phaeomarasmius as members of the family. © The Willi Hennig Society 2009.

Stowaway MITEs in Hordeum (Poaceae): evolutionary history, ancestral elements and classification.

Given the lack of direct observational data relating to transposition of Stowaway miniature inverted repeat transposable elements, phylogenetic methods may provide a means of generating data that adds to our knowledge of these elements. In a phylogenetic framework the evolutionary history of homologous elements may be traced, and the nucleotide sequence of elements at or close to the time of insertion can be reconstructed. Based on a phylogeny of the diploid species of the genus Hordeum we explore evolutionary aspects of four non-homologous groups of Stowaway elements inserted into three nuclear genes: nucellin, xylose isomerase, and barley leucine zipper 1. The data illustrate how elements starting from a high degree of sequence similarity between terminal inverted repeat regions gradually degrade, and confirm previous notions about preferential insertion at particular TA target sites. It is shown how creation of consensus sequences as estimates of ancestral elements may be positively misleading. The Stowaway family of transposable elements is often further divided into subfamilies based on sequence similarity between elements. Sequence similarity data from the elements discovered in the xylose isomerase gene, and other elements found through BLAST searches in GenBank, reveal inconsistency in the rules used for classification. In order to reflect natural groups, a classification of transposable elements must be based on phylogenetic evidence rather than raw similarity. © The Willi Hennig Society 2009.

Genome Reports: Contracted Genes and Dwarfed Plastome in Mycoheterotrophic Sciaphila thaidanica (Triuridaceae, Pandanales).

With a reduced need for photosynthesis, the plastome of parasitic and mycoheterotrophic plants degrades. In the tiny, fully mycoheterotrophic plant Sciaphila thaidanica, we find one of the smallest plastomes yet encountered. Its size is just 12,780 bp and it contains only 20 potentially functional housekeeping genes. Thus S. thaidanica fits the proposed model of gene loss in achlorophyllous plants. The most astonishing feature of the plastome is its extremely compact nature, with more than half of the genes having overlapping reading frames. Additionally, intergenic sequences have been reduced to a bare minimum, and the retained genes have been reduced in length both compared with the orthologous genes in another mycoheterotrophic species of Sciaphila and in the autotrophic relative Carludovica.

Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes.

The order Alismatales is a hotspot for evolution of plant mitochondrial genomes characterized by remarkable differences in genome size, substitution rates, RNA editing, retrotranscription, gene loss and intron loss. Here we have sequenced the complete mitogenomes of Zostera marina and Stratiotes aloides, which together with previously sequenced mitogenomes from Butomus and Spirodela, provide new evolutionary evidence of genome size reduction, gene loss and transfer to the nucleus. The Zostera mitogenome includes a large portion of DNA transferred from the plastome, yet it is the smallest known mitogenome from a non-parasitic plant. Using a broad sample of the Alismatales, the evolutionary history of ribosomal protein gene loss is analyzed. In Zostera almost all ribosomal protein genes are lost from the mitogenome, but only some can be found in the nucleus.

Localized Retroprocessing as a Model of Intron Loss in the Plant Mitochondrial Genome.

Loss of introns in plant mitochondrial genes is commonly explained by retroprocessing. Under this model, an mRNA is reverse transcribed and integrated back into the genome, simultaneously affecting the contents of introns and edited sites. To evaluate the extent to which retroprocessing explains intron loss, we analyzed patterns of intron content and predicted RNA editing for whole mitochondrial genomes of 30 species in the monocot order Alismatales. In this group, we found an unusually high degree of variation in the intron content, even expanding the hitherto known variation among angiosperms. Some species have lost some two-third of the cis-spliced introns. We found a strong correlation between intron content and editing frequency, and detected 27 events in which intron loss is consistent with the presence of nucleotides in an edited state, supporting retroprocessing. However, we also detected seven cases of intron loss not readily being explained by retroprocession. Our analyses are also not consistent with the entire length of a fully processed cDNA copy being integrated into the genome, but instead indicate that retroprocessing usually occurs for only part of the gene. In some cases, several rounds of retroprocessing may explain intron loss in genes completely devoid of introns. A number of taxa retroprocessing seem to be very common and a possibly ongoing process. It affects the entire mitochondrial genome.