Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera.

PREMISE: Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants. METHODS: Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities. RESULTS: The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae. CONCLUSIONS: Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Evolution of Whirly1 in the angiosperms: sequence, splicing, and expression in a clade of early transitional mycoheterotrophic orchids.

The plastid-targeted transcription factor Whirly1 (WHY1) has been implicated in chloroplast biogenesis, plastid genome stability, and fungal defense response, which together represent characteristics of interest for the study of autotrophic losses across the angiosperms. While gene loss in the plastid and nuclear genomes has been well studied in mycoheterotrophic plants, the evolution of the molecular mechanisms impacting genome stability is completely unknown. Here, we characterize the evolution of WHY1 in four early transitional mycoheterotrophic orchid species in the genus Corallorhiza by synthesizing the results of phylogenetic, transcriptomic, and comparative genomic analyses with WHY1 genomic sequences sampled from 21 orders of angiosperms. We found an increased number of non-canonical WHY1 isoforms assembled from all but the greenest Corallorhiza species, including intron retention in some isoforms. Within Corallorhiza, phylotranscriptomic analyses revealed the presence of tissue-specific differential expression of WHY1 in only the most photosynthetically capable species and a coincident increase in the number of non-canonical WHY1 isoforms assembled from fully mycoheterotrophic species. Gene- and codon-level tests of WHY1 selective regimes did not infer significant signal of either relaxed selection or episodic diversifying selection in Corallorhiza but did so for relaxed selection in the late-stage full mycoheterotrophic orchids Epipogium aphyllum and Gastrodia elata. Additionally, nucleotide substitutions that most likely impact the function of WHY1, such as nonsense mutations, were only observed in late-stage mycoheterotrophs. We propose that our findings suggest that splicing and expression changes may precede the selective shifts we inferred for late-stage mycoheterotrophic species, which therefore does not support a primary role for WHY1 in the transition to mycoheterotrophy in the Orchidaceae. Taken together, this study provides the most comprehensive view of WHY1 evolution across the angiosperms to date.

Organellar phylogenomics at the epidendroid orchid base, with a focus on the mycoheterotrophic Wullschlaegelia.

BACKGROUND AND AIMS: Heterotrophic plants have long been a challenge for systematists, exemplified by the base of the orchid subfamily Epidendroideae, which contains numerous mycoheterotrophic species. METHODS: Here we address the utility of organellar genomes in resolving relationships at the epidendroid base, specifically employing models of heterotachy, or lineage-specific rate variation over time. We further conduct comparative analyses of plastid genome evolution in heterotrophs and structural variation in matK. KEY RESULTS: We present the first complete plastid genomes (plastomes) of Wullschlaegelia, the sole genus of the tribe Wullschlaegelieae, revealing a highly reduced genome of 37 kilobases, which retains a fraction of the genes present in related autotrophs. Plastid phylogenomic analyses recovered a strongly supported clade composed exclusively of mycoheterotrophic species with long branches. We further analyzed mitochondrial gene sets, which recovered similar relationships to those in other studies using nuclear data, but the placement of Wullschlaegelia remains uncertain. We conducted comparative plastome analyses among Wullschlaegelia and other heterotrophic orchids, revealing a suite of correlated substitutional and structural changes relative to autotrophic species. Lastly, we investigated evolutionary and structural variation in matK, which is retained in Wullschlaegelia and a few other 'late stage' heterotrophs and found evidence for structural conservation despite rapid substitution rates in both Wullschlaegelia and the leafless Gastrodia. CONCLUSIONS: Our analyses reveal the limits of what the plastid genome can tell us on orchid relationships in this part of the tree, even when applying parameter-rich heterotachy models. Our study underscores the need for increased taxon sampling across all three genomes at the epidendroid base, and illustrates the need for further research on addressing heterotachy in phylogenomic analyses.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae.

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have  been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

Morphologic and Molecular Heterogeneity of High-grade Serous Carcinoma Precursor Lesions.

Serous tubal intraepithelial carcinoma (STIC) is the fallopian tube precursor lesion for most cases of pelvic high-grade serous carcinoma (HGSC). To date, the morphologic, molecular, and clinical heterogeneity of STIC and a less atypical putative precursor lesion, termed serous tubal intraepithelial lesion, has not been well characterized. Better understanding of precursor heterogeneity could impact the clinical management of women with incidental STICs (without concurrent carcinoma) identified in cases of prophylactic or opportunistic salpingectomy. This study analyzed morphologic and molecular features of 171 STICs and 21 serous tubal intraepithelial lesions. We assessed their histologic features, Ki-67 and p53 staining patterns, and genome-wide DNA copy number alterations. We classified all precursor lesions into 2 morphologic subtypes, one with a flat surface (Flat) and the other characterized by budding, loosely adherent, or detached (BLAD) morphology. On the basis of pathology review by a panel of 8 gynecologic pathologists, we found 87 BLAD, 96 Flat, and 9 indeterminate lesions. As compared with Flat lesions, BLAD lesions were more frequently diagnostic of STIC ( P <0.0001) and were found concurrently with HGSC ( P <0.0001). BLAD morphology was also characterized by higher Ki-67 proliferation index ( P <0.0001), presence of epithelial stratification ( P <0.0001), and increased lymphocyte density ( P <0.0001). BLAD lesions also exhibited more frequent DNA copy number gain/amplification at the CCNE1 or CMYC loci canonical to HGSCs ( P <0.0001). Both BLAD morphology and STIC diagnoses are independent risk factors for an elevated Ki-67 proliferation index. No correlation was observed between BLAD and Flat lesions with respect to patient age, presence of germline BRCA1/2 mutation, or p53 staining pattern. These findings suggest that tubal precursor lesions are morphologically and molecularly heterogeneous, laying the foundation for further studies on the pathogenesis of HGSC initiation and identifying histologic features predictive of poor patient outcomes.

A plastid phylogenomic framework for the palm family (Arecaceae).

BACKGROUND: Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research. RESULTS: Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support. CONCLUSIONS: The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

A lack of population structure characterizes the invasive Lonicera japonica in West Virginia and across eastern North America.

Invasive plant species cause massive ecosystem damage globally, yet represent powerful case studies in population genetics and rapid adaptation to new habitats. The availability of digitized herbarium collections data, and the ubiquity of invasive species across the landscape make them highly accessible for studies of invasion history and population dynamics associated with their introduction, establishment, spread, and ecological interactions. Here we focus on Lonicera japonica, one of the most damaging invasive vine species in North America. We leveraged digitized collections data and contemporary field collections to reconstruct the invasion history and characterize patterns of genomic variation in the eastern USA, using a straightforward method for generating nucleotide polymorphism data and a recently published, chromosome-level genome for the species. We found an overall lack of population structure among sites in northern West Virginia, USA, as well as across sites in the central and eastern USA. Heterozygosity and population differentiation were both low based on Fst, analysis of molecular variance, principal components analysis, and cluster-based analyses. We also found evidence of high inbreeding coefficients and significant linkage disequilibrium, in line with the ability of this otherwise outcrossing, perennial species to propagate vegetatively. Our findings corroborate earlier studies based on allozyme data, and suggest that intentional, human-assisted spread explains the lack of population structure, as this species was planted for erosion control and as an ornamental, escaping cultivation repeatedly across the USA. Finally, we discuss how plant invasion genomics can be incorporated into experiential undergraduate education as a way to integrate teaching and research.

Mitochondrial genome sequencing and analysis of the invasive Microstegium vimineum: a resource for systematics, invasion history, and management.

Premise of the Research: Plants remain underrepresented among species with sequenced mitochondrial genomes (mitogenomes), due to the difficulty in assembly with short-read technology. Invasive species lag behind crops and other economically important species in this respect, representing a lack of tools for management and land conservation efforts. Methodology: The mitogenome of Microstegium vimineum, one of the most damaging invasive plant species in North America, was sequenced and analyzed using long-read data, providing a resource for biologists and managers. We conducted analyses of genome content, phylogenomic analyses among grasses and relatives based on mitochondrial coding regions, and an analysis of mitochondrial single nucleotide polymorphism in this invasive grass species. Pivotal Results: The assembly is 478,010 bp in length and characterized by two large, inverted repeats, and a large, direct repeat. However, the genome could not be circularized, arguing against a "master circle" structure. Long-read assemblies with data subsets revealed several alternative genomic conformations, predominantly associated with large repeats. Plastid-like sequences comprise 2.4% of the genome, with further evidence of Class I and Class II transposable element-like sequences. Phylogenetic analysis placed M. vimineum with other Microstegium species, excluding M. nudum, but with weak support. Analysis of polymorphic sites across 112 accessions of M. vimineum from the native and invasive ranges revealed a complex invasion history. Conclusions: We present an in-depth analysis of mitogenome structure, content, phylogenetic relationships, and range-wide genomic variation in M. vimineum's invasive US range. The mitogenome of M. vimineum is typical of other andropogonoid grasses, yet mitochondrial sequence variation across the invasive and native ranges is extensive. Our findings suggest multiple introductions to the US over the last century, with subsequent spread, secondary contact, long-distance dispersal, and possibly post-invasion selection on awn phenotypes. Efforts to produce genomic resources for invasive species, including sequenced mitochondrial genomes, will continue to provide tools for their effective management, and to help predict and prevent future invasions.

Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages.

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes.

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

Lineage and role in integrative taxonomy of a heterotrophic orchid complex.

Lineage-based species definitions applying coalescent approaches to species delimitation have become increasingly popular. Yet, the application of these methods and the recognition of lineage-only definitions have recently been questioned. Species delimitation criteria that explicitly consider both lineages and evidence for ecological role shifts provide an opportunity to incorporate ecologically meaningful data from multiple sources in studies of species boundaries. Here, such criteria were applied to a problematic group of mycoheterotrophic orchids, the Corallorhiza striata complex, analysing genomic, morphological, phenological, reproductive-mode, niche, and fungal host data. A recently developed method for generating genomic polymorphism data-ISSRseq-demonstrates evidence for four distinct lineages, including a previously unidentified lineage in the Coast Ranges and Cascades of California and Oregon, USA. There is divergence in morphology, phenology, reproductive mode, and fungal associates among the four lineages. Integrative analyses, conducted in population assignment and redundancy analysis frameworks, provide evidence of distinct genomic lineages and a similar pattern of divergence in the extended data, albeit with weaker signal. However, none of the extended data sets fully satisfy the condition of a significant role shift, which requires evidence of fixed differences. The four lineages identified in the current study are recognized at the level of variety, short of comprising different species. This study represents the most comprehensive application of lineage + role to date and illustrates the advantages of such an approach.

The ecology of palm genomes: repeat-associated genome size expansion is constrained by aridity.

Genome size varies 2400-fold across plants, influencing their evolution through changes in cell size and cell division rates which impact plants' environmental stress tolerance. Repetitive element expansion explains much genome size diversity, and the processes structuring repeat 'communities' are analogous to those structuring ecological communities. However, which environmental stressors influence repeat community dynamics has not yet been examined from an ecological perspective. We measured genome size and leveraged climatic data for 91% of genera within the ecologically diverse palm family (Arecaceae). We then generated genomic repeat profiles for 141 palm species, and analysed repeats using phylogenetically informed linear models to explore relationships between repeat dynamics and environmental factors. We show that palm genome size and repeat 'community' composition are best explained by aridity. Specifically, Ty3-gypsy and TIR elements were more abundant in palm species from wetter environments, which generally had larger genomes, suggesting amplification. By contrast, Ty1-copia and LINE elements were more abundant in drier environments. Our results suggest that water stress inhibits repeat expansion through selection on upper genome size limits. However, elements that may associate with stress-response genes (e.g. Ty1-copia) have amplified in arid-adapted palm species. Overall, we provide novel evidence of climate influencing the assembly of repeat 'communities'.

Digitized collections elucidate invasion history and patterns of awn polymorphism in Microstegium vimineum.

PREMISE: Digitized collections can help illuminate the mechanisms behind the establishment and spread of invasive plants. These databases provide a record of traits in space and time that allows for investigation of abiotic and biotic factors that influence invasive species. METHODS: Over 1100 digitized herbarium records were examined to investigate the invasion history and trait variation of Microstegium vimineum. Presence-absence of awns was investigated to quantify geographic patterns of this polymorphic trait, which serves several functions in grasses, including diaspore burial and dispersal to germination sites. Floret traits were further quantified, and genomic analyses of contemporary samples were conducted to investigate the history of M. vimineum's introduction and spread into North America. RESULTS: Herbarium records revealed similar patterns of awn polymorphism in native and invaded ranges of M. vimineum, with awned forms predominating at higher latitudes and awnless forms at lower latitudes. Herbarium records and genomic data suggested initial introduction and spread of the awnless form in the southeastern United States, followed by a putative secondary invasion and spread of the awned form from eastern Pennsylvania. Awned forms have longer florets, and floret size varies significantly with latitude. There is evidence of a transition zone with short-awned specimens at mid-latitudes. Genomic analyses revealed two distinct clusters corresponding to awnless and awned forms, with evidence of admixture. CONCLUSIONS: Our results demonstrate the power of herbarium data to elucidate the invasion history of a problematic weed in North America and, together with genomic data, reveal a possible key trait in introduction success: presence or absence of an awn.

Chromosome Level Genome Assembly and Annotation of Highly Invasive Japanese Stiltgrass (Microstegium vimineum).

The invasive Japanese stiltgrass (Microstegium vimineum) affects a wide range of ecosystems and threatens biodiversity across the eastern USA. However, the mechanisms underlying rapid adaptation, plasticity, and epigenetics in the invasive range are largely unknown. We present a chromosome-level assembly for M. vimineum to investigate genome dynamics, evolution, adaptation, and the genomics of phenotypic plasticity. We generated a 1.12-Gb genome with scaffold N50 length of 53.44 Mb respectively, taking a de novo assembly approach that combined PacBio and Dovetail Genomics Omni-C sequencing. The assembly contains 23 pseudochromosomes, representing 99.96% of the genome. BUSCO assessment indicated that 80.3% of Poales gene groups are present in the assembly. The genome is predicted to contain 39,604 protein-coding genes, of which 26,288 are functionally annotated. Furthermore, 66.68% of the genome is repetitive, of which unclassified (35.63%) and long-terminal repeat (LTR) retrotransposons (26.90%) are predominant. Similar to other grasses, Gypsy (41.07%) and Copia (32%) are the most abundant LTR-retrotransposon families. The majority of LTR-retrotransposons are derived from a significant expansion in the past 1-2 Myr, suggesting the presence of relatively young LTR-retrotransposon lineages. We find corroborating evidence from Ks plots for a stiltgrass-specific duplication event, distinct from the more ancient grass-specific duplication event. The assembly and annotation of M. vimineum will serve as an essential genomic resource facilitating studies of the invasion process, the history and consequences of polyploidy in grasses, and provides a crucial tool for natural resource managers.

Automated outbreak detection of hospital-associated pathogens: Value to infection prevention programs.

OBJECTIVE: To assess the utility of an automated, statistically-based outbreak detection system to identify clusters of hospital-acquired microorganisms. DESIGN: Multicenter retrospective cohort study. SETTING: The study included 43 hospitals using a common infection prevention surveillance system. METHODS: A space-time permutation scan statistic was applied to hospital microbiology, admission, discharge, and transfer data to identify clustering of microorganisms within hospital locations and services. Infection preventionists were asked to rate the importance of each cluster. A convenience sample of 10 hospitals also provided information about clusters previously identified through their usual surveillance methods. RESULTS: We identified 230 clusters in 43 hospitals involving Gram-positive and -negative bacteria and fungi. Half of the clusters progressed after initial detection, suggesting that early detection could trigger interventions to curtail further spread. Infection preventionists reported that they would have wanted to be alerted about 81% of these clusters. Factors associated with clusters judged to be moderately or highly concerning included high statistical significance, large size, and clusters involving Clostridioides difficile or multidrug-resistant organisms. Based on comparison data provided by the convenience sample of hospitals, only 9 (18%) of 51 clusters detected by usual surveillance met statistical significance, and of the 70 clusters not previously detected, 58 (83%) involved organisms not routinely targeted by the hospitals' surveillance programs. All infection prevention programs felt that an automated outbreak detection tool would improve their ability to detect outbreaks and streamline their work. CONCLUSIONS: Automated, statistically-based outbreak detection can increase the consistency, scope, and comprehensiveness of detecting hospital-associated transmission.

Plastid genomes of the North American Rhus integrifolia-ovata complex and phylogenomic implications of inverted repeat structural evolution in Rhus L.

Plastid genomes (plastomes) represent rich sources of information for phylogenomics, from higher-level studies to below the species level. The genus Rhus (sumac) has received a significant amount of study from phylogenetic and biogeographic perspectives, but genomic studies in this genus are lacking. Rhus integrifolia and R. ovata are two shrubby species of high ecological importance in the southwestern USA and Mexico, where they occupy coastal scrub and chaparral habitats. They hybridize frequently, representing a fascinating system in which to investigate the opposing effects of hybridization and divergent selection, yet are poorly characterized from a genomic perspective. In this study, complete plastid genomes were sequenced for one accession of R. integrifolia and one each of R. ovata from California and Arizona. Sequence variation among these three accessions was characterized, and PCR primers potentially useful in phylogeographic studies were designed. Phylogenomic analyses were conducted based on a robustly supported phylogenetic framework based on 52 complete plastomes across the order Sapindales. Repeat content, rather than the size of the inverted repeat, had a stronger relative association with total plastome length across Sapindales when analyzed with phylogenetic least squares regression. Variation at the inverted repeat boundary within Rhus was striking, resulting in major shifts and independent gene losses. Specifically, rps19 was lost independently in the R. integrifolia-ovata complex and in R. chinensis, with a further loss of rps22 and a major contraction of the inverted repeat in two accessions of the latter. Rhus represents a promising novel system to study plastome structural variation of photosynthetic angiosperms at and below the species level.

Ancient Mitochondrial Gene Transfer between Fungi and the Orchids.

The mitochondrial genomes (mitogenomes) of plants are known to incorporate and accumulate DNA from intra- and extracellular donors. Despite the intimate relationships formed between flowing plants (angiosperms) and fungi, lengthy fungal-like sequence has not been identified in angiosperm mitogenomes to date. Here, we present multiple lines of evidence documenting horizontal gene transfer (HGT) between the mitogenomes of fungi and the ancestors of the orchids, plants that are obligate parasites of fungi during their early development. We show that the ancestor of the orchids acquired an ∼270-bp fungal mitogenomic region containing three transfer RNA genes. We propose that the short HGT was later replaced by a second HGT event transferring >8 kb and 14 genes from a fungal mitogenome to that of the ancestor of the largest orchid subfamily, Epidendroideae. Our results represent the first evidence of genomic-scale HGT between fungal and angiosperm mitogenomes and demonstrate that the length intergenic spacer regions of angiosperm mitogenomes can effectively fossilize the genomic remains of ancient, nonplant organisms.

Unprecedented Parallel Photosynthetic Losses in a Heterotrophic Orchid Genus.

Heterotrophic plants are evolutionary experiments in genomic, morphological, and physiological change. Yet, genomic sampling gaps exist among independently derived heterotrophic lineages, leaving unanswered questions about the process of genome modification. Here, we have sequenced complete plastid genomes for all species of the leafless orchid genus Hexalectris, including multiple individuals for most, and leafy relatives Basiphyllaea and Bletia. Our objectives are to determine the number of independent losses of photosynthesis and to test hypotheses on the process of genome degradation as a result of relaxed selection. We demonstrate four to five independent losses of photosynthesis in Hexalectris based on degradation of the photosynthetic apparatus, with all but two species displaying evidence of losses, and variation in gene loss extending below the species level. Degradation in the atp complex is advanced in Hexalectris warnockii, whereas only minimal degradation (i.e., physical loss) has occurred among some "housekeeping" genes. We find genomic rearrangements, shifts in Inverted Repeat boundaries including complete loss in one accession of H. arizonica, and correlations among substitutional and genomic attributes. Our unprecedented finding of multiple, independent transitions to a fully mycoheterotrophic lifestyle in a single genus reveals that the number of such transitions among land plants is likely underestimated. This study underscores the importance of dense taxon sampling, which is highly informative for advancing models of genome evolution in heterotrophs. Mycoheterotrophs such as Hexalectris provide forward-genetic opportunities to study the consequences of radical genome evolution beyond what is possible with mutational studies in model organisms alone.