Jeffrey L. Boore

Sequencing and comparing whole mitochondrial genomes of animals.

Comparing complete animal mitochondrial genome sequences is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. Not only are they much more informative than shorter sequences of individual genes for inferring evolutionary relatedness, but these data also provide sets of genome-level characters, such as the relative arrangements of genes, which can be especially powerful. We describe here the protocols commonly used for physically isolating mitochondrial DNA (mtDNA), for amplifying these by polymerase chain reaction (PCR) or rolling circle amplification (RCA), for cloning, sequencing, assembly, validation, and gene annotation, and for comparing both sequences and gene arrangements. On several topics, we offer general observations based on our experiences with determining and comparing complete mitochondrial DNA sequences.

Complete plastid genome sequences suggest strong selection for retention of photosynthetic genes in the parasitic plant genus Cuscuta

BackgroundPlastid genome content and protein sequence are highly conserved across land plants and their closest algal relatives. Parasitic plants, which obtain some or all of their nutrition through an attachment to a host plant, are often a striking exception. Heterotrophy can lead to relaxed constraint on some plastid genes or even total gene loss. We sequenced plastid genomes of two species in the parasitic genus Cuscuta along with a non-parasitic relative, Ipomoea purpurea, to investigate changes in the plastid genome that may result from transition to the parasitic lifestyle.ResultsAside from loss of all ndh genes, Cuscuta exaltata retains photosynthetic and photorespiratory genes that evolve under strong selective constraint. Cuscuta obtusiflora has incurred substantially more change to its plastid genome, including loss of all genes for the plastid-encoded RNA polymerase. Despite extensive change in gene content and greatly increased rate of overall nucleotide substitution, C. obtusiflora also retains all photosynthetic and photorespiratory genes with only one minor exception.ConclusionAlthough Epifagus virginiana, the only other parasitic plant with its plastid genome sequenced to date, has lost a largely overlapping set of transfer-RNA and ribosomal genes as Cuscuta, it has lost all genes related to photosynthesis and maintains a set of genes which are among the most divergent in Cuscuta. Analyses demonstrate photosynthetic genes are under the highest constraint of any genes within the plastid genomes of Cuscuta, indicating a function involving RuBisCo and electron transport through photosystems is still the primary reason for retention of the plastid genome in these species.

Complete plastid genome sequences of Drimys, Liriodendron, and Piper: implications for the phylogenetic relationships of magnoliids

BackgroundThe magnoliids with four orders, 19 families, and 8,500 species represent one of the largest clades of early diverging angiosperms. Although several recent angiosperm phylogenetic analyses supported the monophyly of magnoliids and suggested relationships among the orders, the limited number of genes examined resulted in only weak support, and these issues remain controversial. Furthermore, considerable incongruence resulted in phylogenetic reconstructions supporting three different sets of relationships among magnoliids and the two large angiosperm clades, monocots and eudicots. We sequenced the plastid genomes of three magnoliids, Drimys (Canellales), Liriodendron (Magnoliales), and Piper (Piperales), and used these data in combination with 32 other angiosperm plastid genomes to assess phylogenetic relationships among magnoliids and to examine patterns of variation of GC content.ResultsThe Drimys, Liriodendron, and Piper plastid genomes are very similar in size at 160,604, 159,886 bp, and 160,624 bp, respectively. Gene content and order are nearly identical to many other unrearranged angiosperm plastid genomes, including Calycanthus, the other published magnoliid genome. Overall GC content ranges from 34–39%, and coding regions have a substantially higher GC content than non-coding regions. Among protein-coding genes, GC content varies by codon position with 1st codon > 2nd codon > 3rd codon, and it varies by functional group with photosynthetic genes having the highest percentage and NADH genes the lowest. Phylogenetic analyses using parsimony and likelihood methods and sequences of 61 protein-coding genes provided strong support for the monophyly of magnoliids and two strongly supported groups were identified, the Canellales/Piperales and the Laurales/Magnoliales. Strong support is reported for monocots and eudicots as sister clades with magnoliids diverging before the monocot-eudicot split. The trees also provided moderate or strong support for the position of Amborella as sister to a clade including all other angiosperms.ConclusionEvolutionary comparisons of three new magnoliid plastid genome sequences, combined with other published angiosperm genomes, confirm that GC content is unevenly distributed across the genome by location, codon position, and functional group. Furthermore, phylogenetic analyses provide the strongest support so far for the hypothesis that the magnoliids are sister to a large clade that includes both monocots and eudicots.

Domestication of olive fly through a multi-regional host shift to cultivated olives: Comparative dating using complete mitochondrial genomes

The evolutionary history of the olive fly, Bactrocera oleae, was reconstructed in a phylogenetic and coalescent framework using full mitochondrial genome data from 21 individuals covering the entire worldwide distribution of the species. Special attention was given to reconstructing the timing of the processes under study. The early subdivision of the olive fly reflects the Quaternary differentiation between Olea europea subsp. europea in the Mediterranean area and the two lineages of Olea europea subsp. cuspidata in Africa and Asia, pointing to an early and close association between the olive fly and its host. The geographic structure and timing of olive fly differentiation in the Mediterranean indicates a clear connection with the post-glacial recolonization of wild olives in the area, and is irreconcilable with the early historical process of domestication and spread of the cultivated olive from its Levantine origin. Therefore, we suggest an early co-history of the olive fly with its wild host during the Quaternary and post-glacial periods and a multi-regional shift of olive flies to cultivated olives as these cultivars gradually replaced wild olives in historical times.

Systematics and plastid genome evolution of the cryptically photosynthetic parasitic plant genus Cuscuta (Convolvulaceae)

BackgroundThe genus Cuscuta L. (Convolvulaceae), commonly known as dodders, are epiphytic vines that invade the stems of their host with haustorial feeding structures at the points of contact. Although they lack expanded leaves, some species are noticeably chlorophyllous, especially as seedlings and in maturing fruits. Some species are reported as crop pests of worldwide distribution, whereas others are extremely rare and have local distributions and apparent niche specificity. A strong phylogenetic framework for this large genus is essential to understand the interesting ecological, morphological and molecular phenomena that occur within these parasites in an evolutionary context.ResultsHere we present a well-supported phylogeny of Cuscuta using sequences of the nuclear ribosomal internal transcribed spacer and plastid rps2, rbcL and matK from representatives across most of the taxonomic diversity of the genus. We use the phylogeny to interpret morphological and plastid genome evolution within the genus. At least three currently recognized taxonomic sections are not monophyletic and subgenus Cuscuta is unequivocally paraphyletic. Plastid genes are extremely variable with regards to evolutionary constraint, with rbcL exhibiting even higher levels of purifying selection in Cuscuta than photosynthetic relatives. Nuclear genome size is highly variable within Cuscuta, particularly within subgenus Grammica, and in some cases may indicate the existence of cryptic species in this large clade of morphologically similar species.ConclusionSome morphological characters traditionally used to define major taxonomic splits within Cuscuta are homoplastic and are of limited use in defining true evolutionary groups. Chloroplast genome evolution seems to have evolved in a punctuated fashion, with episodes of loss involving suites of genes or tRNAs followed by stabilization of gene content in major clades. Nearly all species of Cuscuta retain some photosynthetic ability, most likely for nutrient apportionment to their seeds, while complete loss of photosynthesis and possible loss of the entire chloroplast genome is limited to a single small clade of outcrossing species found primarily in western South America.

High divergence across the whole mitochondrial genome in the "pan-Antarctic" springtail Friesea grisea: Evidence for cryptic species?

Collembola are one of the few hexapod groups adapted to live in the harsh environmental conditions of Antarctic terrestrial ecosystems. Diversity is limited to a few species that can be very abundant in coastal deglaciated sites. A remarkable lack of overlap in Collembola species composition is evident between Western and Eastern Antarctica, and Friesea grisea is currently the only species whose distribution is thought to span these two main regions of the continent. However, our analysis of the complete sequences of the mitochondrial genomes from specimens obtained from each of the two regions showed unexpected genetic divergence, well above the average levels observed between populations belonging to the same species, and so indicating that these are actually separate species, despite their lack of distinguishing morphology. Detailed analysis of the two genomes showed the presence of a non-coding region observed between trnS(uga) and nad1. Other features of these mitochondrial genomes, such as base compositional bias, secondary structure features of tRNAs and the presence of regulatory elements in the control region, are described and discussed from an evolutionary standpoint.

Parallel Loss of Plastid Introns and Their Maturase in the Genus Cuscuta

Plastid genome content and arrangement are highly conserved across most land plants and their closest relatives, streptophyte algae, with nearly all plastid introns having invaded the genome in their common ancestor at least 450 million years ago. One such intron, within the transfer RNA trnK-UUU, contains a large open reading frame that encodes a presumed intron maturase, matK. This gene is missing from the plastid genomes of two species in the parasitic plant genus Cuscuta but is found in all other published land plant and streptophyte algal plastid genomes, including that of the nonphotosynthetic angiosperm Epifagus virginiana and two other species of Cuscuta. By examining matK and plastid intron distribution in Cuscuta, we add support to the hypothesis that its normal role is in splicing seven of the eight group IIA introns in the genome. We also analyze matK nucleotide sequences from Cuscuta species and relatives that retain matK to test whether changes in selective pressure in the maturase are associated with intron deletion. Stepwise loss of most group IIA introns from the plastid genome results in substantial change in selective pressure within the hypothetical RNA-binding domain of matK in both Cuscuta and Epifagus, either through evolution from a generalist to a specialist intron splicer or due to loss of a particular intron responsible for most of the constraint on the binding region. The possibility of intron-specific specialization in the X-domain is implicated by evidence of positive selection on the lineage leading to C. nitida in association with the loss of six of seven introns putatively spliced by matK. Moreover, transfer RNA gene deletion facilitated by parasitism combined with an unusually high rate of intron loss from remaining functional plastid genes created a unique circumstance on the lineage leading to Cuscuta subgenus Grammica that allowed elimination of matK in the most species-rich lineage of Cuscuta.

Analysis of the complete mitochondrial genome sequences of the soybean rust pathogens Phakopsora pachyrhizi and P. meibomiae

The mitochondrial (mt) genomes of two soybean rust pathogens, Phakopsora pachyrhizi and P. meibomiae, have been sequenced. The mt genome of P. pachyrhizi is a circular 31 825-bp molecule with a mean GC content of 34.6%, while P. meibomiae possesses a 32 520-bp circular molecule with a mean GC content of 34.9%. Both mt genomes contain the genes encoding ATP synthase subunits 6, 8 and 9 (atp6, atp8 and atp9), cytochrome oxidase subunits I, II and III (cox1, cox2 and cox3), apocytochrome b (cob), reduced nicotinamide adenine dinucleotide ubiquinone oxidoreductase subunits (nad1, nad2, nad3, nad4, nad4L, nad5 and nad6), the large and small mt ribosomal RNA genes, one ORF coding for a ribosomal protein (rps3), and a set of 24 tRNA genes that recognize codons for all amino acids. The order of the protein-coding genes and tRNA is identical in the two Phakopsora species, and all genes are transcribed from the same DNA strand clockwise. Introns were identified in the cox1, cob and rnl genes of both species, with three of the introns having ORFs with motifs similar to the LAGLIDADG endonucleases of other fungi. Phylogenetic analysis of the 14 shared protein-coding genes agrees with commonly accepted fungal taxonomy.

Development of simple sequence repeat markers for the soybean rust fungus, Phakopsora pachyrhizi

Twenty‐four simple sequence repeat markers were developed for Phakopsora pachyrhizi, a fungal pathogen of soybean (Glycine max) and other legumes. All 24 of the loci were evaluated on 28 isolates of P. pachyrhizi. Twenty‐one loci were polymorphic, with allelic diversity ranging from two to eight alleles, and null alleles were observed for eight of the 24 loci. A preliminary screen with the closely related species, P. meibomiae, indicated that these primer pairs are specific to P. pachyrhizi.

DETECTING EVOLUTIONARY TRANSFER OF GENES USING PhIGs1

Organisms have acquired plastids by convoluted paths that have provided multiple opportunities for gene transfer into a host nucleus from intracellular organisms, including the cyanobacterial ancestor of plastids, the proteobacterial ancestor of mitochondria, and both green and red algae whose engulfment has led to secondary acquisition of plastids. These gene movements are most accurately demonstrated by building phylogenetic trees that identify the evolutionary origin of each gene, and one effective tool for this is “PhIGs” (Phylogenetically Inferred Groups; http://PhIGs.org), a set of databases and computer tools with a Web interface for whole‐genome evolutionary analysis. PhIGs takes as input gene sets of completely sequenced genomes, builds clusters of genes using a novel, graph‐based approach, and reconstructs the evolutionary relationships among all gene families. The user can view and download the sequence alignments, compare intron‐exon structures, and follow links to functional genomic databases. Currently, PhIGs contains 652,756 genes from 45 genomes grouped into 61,059 gene families. Graphical displays show the relative positions of these genes among genomes. PhIGs has been used to detect the evolutionary transfer of hundreds of genes from cyanobacteria and red algae into oömycete nuclear genomes, revealing that even though they have no plastids, their ancestors did, having secondarily acquired them from an intracellular red alga. A great number of genomes are soon to become available that are relevant to our broader understanding of the movement of genes among intracellular compartments after engulfing other organisms, and PhIGs will be an effective tool to interpret these gene movements.