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 approximately 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 more than 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, non-plant organisms.
Alignment of orchid HGT extents
Alignment of orchid mitochondrial sequences representing the extent of the fungal HGT region for each taxon, as identified via BLAST.
Orchid_HGT_extents_alignment_FINAL.phy
Fungal and orchid atp6 alignment
Alignment of fungal and orchid atp6 sequences. Genbank accession number are appended to names of fungal sequences.
atp6_orchid_fungi_alignment_curated.phy
Ustilago maydis coxI and Dendrobium catenatum fungal coxI alignment
Alignment of coxI of Ustilago maydis and the fungal-derived coxI remnant of the Dendrobium catenatum mitogenome.
dendrobium_ustilago_cox1_alignment.fasta
Alignment of Erycina pusilla mitogenomic BACs
Alignment of four Genbank-archived BACs representing mitogenomic sequences of Erycina pusilla. Each sequence is named using its GenBank accession number.
Erycina_mito_HGT_BAC_clones alignment.phy
Hexalectris brevicaulis vs. NCBI Basidiomycete RefSeq Genomes BLAST Result
BLAST result summary of Hexalectris brevicaulis fungal-derived sequence extent queried against the NCBI RefSeq Genomes database. Search conducted 13 March 2019.
NCBI Blast_Hexa_brev_extent_vs_Basidiomycota_RefSeq_Genomes.pdf
Corallorhiza striata var. striata 187MT Kraken2 results
Kraken2 results of 100bp segments of long range PCR of mitogenomic sequence of Corallorhiza striata var. striata (187MT) containing fungal-derived sequence.
C_striata_187MT_100bp_interval_IDs.txt
Alignment of tRNA region common across orchid fungal-derived regions
Alignment of fungal-derived sequence containing the same tRNA genes across orchid taxa sampled.
HGT_core_tRNAs_aligned_genera.phy