Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape.

The cistrome is the complete set of transcription factor (TF) binding sites (cis-elements) in an organism, while an epicistrome incorporates tissue-specific DNA chemical modifications and TF-specific chemical sensitivities into these binding profiles. Robust methods to construct comprehensive cistrome and epicistrome maps are critical for elucidating complex transcriptional networks that underlie growth, behavior, and disease. Here, we describe DNA affinity purification sequencing (DAP-seq), a high-throughput TF binding site discovery method that interrogates genomic DNA with in-vitro-expressed TFs. Using DAP-seq, we defined the Arabidopsis cistrome by resolving motifs and peaks for 529 TFs. Because genomic DNA used in DAP-seq retains 5-methylcytosines, we determined that >75% (248/327) of Arabidopsis TFs surveyed were methylation sensitive, a property that strongly impacts the epicistrome landscape. DAP-seq datasets also yielded insight into the biology and binding site architecture of numerous TFs, demonstrating the value of DAP-seq for cost-effective cistromic and epicistromic annotation in any organism.

Dynamics and function of DNA methylation in plants.

DNA methylation is a conserved epigenetic modification that is important for gene regulation and genome stability. Aberrant patterns of DNA methylation can lead to plant developmental abnormalities. A specific DNA methylation state is an outcome of dynamic regulation by de novo methylation, maintenance of methylation and active demethylation, which are catalysed by various enzymes that are targeted by distinct regulatory pathways. In this Review, we discuss DNA methylation in plants, including methylating and demethylating enzymes and regulatory factors, and the coordination of methylation and demethylation activities by a so-called methylstat mechanism; the functions of DNA methylation in regulating transposon silencing, gene expression and chromosome interactions; the roles of DNA methylation in plant development; and the involvement of DNA methylation in plant responses to biotic and abiotic stress conditions.

The complex polyploid genome architecture of sugarcane.

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

Brassinosteroid-dependent phosphorylation of PHOSPHATE STARVATION RESPONSE2 reduces its DNA-binding ability in rice.

Brassinosteroids (BRs) are widely used as plant growth regulators in modern agriculture. Understanding how BRs regulate nutrient signaling is crucial for reducing fertilizer usage. Here we elucidate that the central BR signaling inhibitor GSK3/SHAGGY-LIKE KINASE2 (GSK2) interacts directly with and phosphorylates PHOSPHATE STARVATION RESPONSE2 (OsPHR2), the key regulator of phosphate (Pi) signaling, to suppress its transcription factor activity in rice (Oryza sativa). We identify a critical phosphorylation site at serine residue S269 of OsPHR2 and demonstrate that phosphorylation by GSK2 or phosphor-mimic mutation of S269 substantially impairs the DNA-binding activity of OsPHR2, and thus diminishes expression of OsPHR2-induced genes and reduces Pi levels. Like BRs, Pi starvation noticeably induces GSK2 instability. We further show that this site-specific phosphorylation event is conserved in Arabidopsis (Arabidopsis thaliana), but varies among the PHR-family members, being present only in most land plants. These results unveil a distinctive post-transcriptional regulatory mechanism in Pi signaling by which BRs promote Pi acquisition, with a potential contribution to the environmental adaptability of plants during their evolution.

Reaction Mechanisms of Pol IV, RDR2, and DCL3 Drive RNA Channeling in the siRNA-Directed DNA Methylation Pathway.

In eukaryotes with multiple small RNA pathways, the mechanisms that channel RNAs within specific pathways are unclear. Here, we reveal the reactions that account for channeling in the small interfering RNA (siRNA) biogenesis phase of the Arabidopsis RNA-directed DNA methylation pathway. The process begins with template DNA transcription by NUCLEAR RNA POLYMERASE IV (Pol IV), whose atypical termination mechanism, induced by nontemplate DNA base-pairing, channels transcripts to the associated RNA-dependent RNA polymerase RDR2. RDR2 converts Pol IV transcripts into double-stranded RNAs and then typically adds an extra untemplated 3' terminal nucleotide to the second strands. The dicer endonuclease DCL3 cuts resulting duplexes to generate 24- and 23-nt siRNAs. The 23-nt RNAs bear the untemplated terminal nucleotide of the RDR2 strand and are underrepresented among ARGONAUTE4-associated siRNAs. Collectively, our results provide mechanistic insights into Pol IV termination, Pol IV-RDR2 coupling, and RNA channeling, from template DNA transcription to siRNA strand discrimination.

Genome copy number predicts extreme evolutionary rate variation in plant mitochondrial DNA.

Nuclear and organellar genomes can evolve at vastly different rates despite occupying the same cell. In most bilaterian animals, mitochondrial DNA (mtDNA) evolves faster than nuclear DNA, whereas this trend is generally reversed in plants. However, in some exceptional angiosperm clades, mtDNA substitution rates have increased up to 5,000-fold compared with closely related lineages. The mechanisms responsible for this acceleration are generally unknown. Because plants rely on homologous recombination to repair mtDNA damage, we hypothesized that mtDNA copy numbers may predict evolutionary rates, as lower copy numbers may provide fewer templates for such repair mechanisms. In support of this hypothesis, we found that copy number explains 47% of the variation in synonymous substitution rates of mtDNA across 60 diverse seed plant species representing ~300 million years of evolution. Copy number was also negatively correlated with mitogenome size, which may be a cause or consequence of mutation rate variation. Both relationships were unique to mtDNA and not observed in plastid DNA. These results suggest that homologous recombinational repair plays a role in driving mtDNA substitution rates in plants and may explain variation in mtDNA evolution more broadly across eukaryotes. Our findings also contribute to broader questions about the relationships between mutation rates, genome size, selection efficiency, and the drift-barrier hypothesis.

Transgenerational increases in DNA methylation in Arabidopsis plants defective in active DNA demethylation.

Spontaneous gain or loss of DNA methylation occurs in plant and animal genomes, and DNA methylation changes can lead to meiotically stable epialleles that generate heritable phenotypic diversity. However, it is unclear whether transgenerational epigenetic stability may be regulated by any cellular factors. Here, we examined spontaneously occurring variations in DNA methylation in wild-type and ros1 mutant Arabidopsis plants that were propagated for ten generations from single-seed descent. We found that the ros1 mutant, which is defective in active DNA demethylation, showed an increased transgenerational epimutation rate. The ros1 mutation led to more spontaneously gained methylation than lost methylation at individual cytosines, compared to the wild type which had similar numbers of spontaneously gained and lost methylation cytosines. Consistently, transgenerational differentially methylated regions were also biased toward hypermethylation in the ros1 mutant. Our results reveal a genetic contribution of the ROS1 DNA demethylase to transgenerational epigenetic stability and suggest that ROS1 may have an unexpected surveillance function in preventing transgenerational DNA methylation increases.

Simultaneous profiling of chromatin accessibility and DNA methylation in complete plant genomes using long-read sequencing.

Epigenetic regulations, including chromatin accessibility, nucleosome positioning and DNA methylation intricately shape genome function. However, current chromatin profiling techniques relying on short-read sequencing technologies fail to characterise highly repetitive genomic regions and cannot detect multiple chromatin features simultaneously. Here, we performed Simultaneous Accessibility and DNA Methylation Sequencing (SAM-seq) of purified plant nuclei. Thanks to the use of long-read nanopore sequencing, SAM-seq enables high-resolution profiling of m6A-tagged chromatin accessibility together with endogenous cytosine methylation in plants. Analysis of naked genomic DNA revealed significant sequence preference biases of m6A-MTases, controllable through a normalisation step. By applying SAM-seq to Arabidopsis and maize nuclei we obtained fine-grained accessibility and DNA methylation landscapes genome-wide. We uncovered crosstalk between chromatin accessibility and DNA methylation within nucleosomes of genes, TEs, and centromeric repeats. SAM-seq also detects DNA footprints over cis-regulatory regions. Furthermore, using the single-molecule information provided by SAM-seq we identified extensive cellular heterogeneity at chromatin domains with antagonistic chromatin marks, suggesting that bivalency reflects cell-specific regulations. SAM-seq is a powerful approach to simultaneously study multiple epigenetic features over unique and repetitive sequences, opening new opportunities for the investigation of epigenetic mechanisms.

The clade-specific target recognition mechanisms of plant RISCs.

Eukaryotic Argonaut proteins (AGOs) assemble RNA-induced silencing complexes (RISCs) with guide RNAs that allow binding to complementary RNA sequences and subsequent silencing of target genes. The model plant Arabidopsis thaliana encodes 10 different AGOs, categorized into three distinct clades based on amino acid sequence similarity. While clade 1 and 2 RISCs are known for their roles in post-transcriptional gene silencing, and clade 3 RISCs are associated with transcriptional gene silencing in the nucleus, the specific mechanisms of how RISCs from each clade recognize their targets remain unclear. In this study, I conducted quantitative binding analyses between RISCs and target nucleic acids with mismatches at various positions, unveiling distinct target binding characteristics unique to each clade. Clade 1 and 2 RISCs require base pairing not only in the seed region but also in the 3' supplementary region for stable target RNA binding, with clade 1 exhibiting a higher stringency. Conversely, clade 3 RISCs tolerate dinucleotide mismatches beyond the seed region. Strikingly, they bind to DNA targets with an affinity equal to or surpassing that of RNA, like prokaryotic AGO complexes. These insights challenge existing views on plant RNA silencing and open avenues for exploring new functions of eukaryotic AGOs.

Diversity of plant DNA in stool is linked to dietary quality, age, and household income.

Eating a varied diet is a central tenet of good nutrition. Here, we develop a molecular tool to quantify human dietary plant diversity by applying DNA metabarcoding with the chloroplast trnL-P6 marker to 1,029 fecal samples from 324 participants across two interventional feeding studies and three observational cohorts. The number of plant taxa per sample (plant metabarcoding richness or pMR) correlated with recorded intakes in interventional diets and with indices calculated from a food frequency questionnaire in typical diets (ρ = 0.40 to 0.63). In adolescents unable to collect validated dietary survey data, trnL metabarcoding detected 111 plant taxa, with 86 consumed by more than one individual and four (wheat, chocolate, corn, and potato family) consumed by >70% of individuals. Adolescent pMR was associated with age and household income, replicating prior epidemiologic findings. Overall, trnL metabarcoding promises an objective and accurate measure of the number and types of plants consumed that is applicable to diverse human populations.

Differential contributions of double-strand break repair pathways to DNA rearrangements following the irradiation of Arabidopsis seeds and seedlings with ion beams.

DNA rearrangements, including inversions, translocations, and large insertions/deletions (indels), are crucial for crop evolution, domestication, and improvement. The rearrangements are frequently induced by ion beams via the mis-repair of DNA double-strand breaks (DSBs). Unfortunately, how ion beam-induced DSBs are repaired has not been comprehensively analyzed and the mechanisms underlying DNA rearrangements remain unclear. In this study, clonal sectors originating from single mutated cells in carbon ion-irradiated plants were used for whole-genome sequencing analyses after Arabidopsis seeds and seedlings were irradiated. Comparative analyses of the induced mutations (e.g., size and frequency of indels and microhomology at the junctions of the rearrangements) in the irradiated materials suggested that the broken/rejoined DSB ends were more extensively processed in seedlings than in seeds. A mutation to canonical non-homologous end-joining (c-NHEJ), which is a DSB repair pathway with minimal processing of DSB ends, increased the sensitivity to ion beams more in the seeds than in the seedlings, which was consistent with the junction analysis results, indicative of the minor contribution of c-NHEJ to the carbon ion-induced DSB repair in seedlings. Considering the characteristics of the large templated insertions in irradiated seedlings, ion-beam-induced DSBs in seedlings are likely repaired primarily by a polymerase theta-mediated pathway. Polymerase theta-deficient seedlings were more sensitive to ion beams than the c-NHEJ-deficient seedlings, consistent with this hypothesis. This study revealed the key characteristics of ion beam-induced DSBs and the associated repair mechanisms related to the physiological status of the irradiated materials, with implications for elucidating the occurrence and induction of rearrangements.

How diverse a monocentric chromosome can be? Repeatome and centromeric organization of Juncus effusus (Juncaceae).

Juncus is the largest genus of Juncaceae and was considered holocentric for a long time. Recent findings, however, indicated that 11 species from different clades of the genus have monocentric chromosomes. Thus, the Juncus centromere organization and evolution need to be reassessed. We aimed to investigate the major repetitive DNA sequences of two accessions of Juncus effusus and its centromeric structure by employing whole-genome analyses, fluorescent in situ hybridization, CENH3 immunodetection, and chromatin immunoprecipitation sequencing. We showed that the repetitive fraction of the small J. effusus genome (~270 Mbp/1C) is mainly composed of Class I and Class II transposable elements (TEs) and satellite DNAs. Three identified satellite DNA families were mainly (peri)centromeric, with two being associated with the centromeric protein CENH3, but not strictly centromeric. Two types of centromere organization were discerned in J. effusus: type 1 was characterized by a single CENH3 domain enriched with JefSAT1-155 or JefSAT2-180, whereas type 2 showed multiple CENH3 domains interrupted by other satellites, TEs or genes. Furthermore, while type 1 centromeres showed a higher degree of satellite identity along the array, type 2 centromeres had less homogenized arrays along the multiple CENH3 domains per chromosome. Although the analyses confirmed the monocentric organization of J. effusus chromosomes, our data indicate a more dynamic arrangement of J. effusus centromeres than observed for other plant species, suggesting it may constitute a transient state between mono- and holocentricity.

The interplay of homology-directed repair pathways in the repair of zebularine-induced DNA-protein crosslinks in Arabidopsis.

DNA-protein crosslinks (DPCs) are highly toxic DNA lesions represented by proteins covalently bound to the DNA. Persisting DPCs interfere with fundamental genetic processes such as DNA replication and transcription. Cytidine analog zebularine (ZEB) has been shown to crosslink DNA METHYLTRANSFERASE1 (MET1). Recently, we uncovered a critical role of the SMC5/6-mediated SUMOylation in the repair of DPCs. In an ongoing genetic screen, we identified two additional candidates, HYPERSENSITIVE TO ZEBULARINE 2 and 3, that were mapped to REGULATOR OF TELOMERE ELONGATION 1 (RTEL1) and polymerase TEBICHI (TEB), respectively. By monitoring the growth of hze2 and hze3 plants in response to zebularine, we show the importance of homologous recombination (HR) factor RTEL1 and microhomology-mediated end-joining (MMEJ) polymerase TEB in the repair of MET1-DPCs. Moreover, genetic interaction and sensitivity assays showed the interdependency of SMC5/6 complex, HR, and MMEJ in the homology-directed repair of MET1-DPCs in Arabidopsis. Altogether, we provide evidence that MET1-DPC repair in plants is more complex than originally expected.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping.

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Genome-wide characterization of single-stranded DNA in rice.

Single-stranded DNA (ssDNA) is essential for various DNA-templated processes in both eukaryotes and prokaryotes. However, comprehensive characterizations of ssDNA still lag in plants compared to nonplant systems. Here, we conducted in situ S1-sequencing, with starting gDNA ranging from 5 µg to 250 ng, followed by comprehensive characterizations of ssDNA in rice (Oryza sativa L.). We found that ssDNA loci were substantially associated with a subset of non-B DNA structures and functional genomic loci. Subtypes of ssDNA loci had distinct epigenetic features. Importantly, ssDNA may act alone or partly coordinate with non-B DNA structures, functional genomic loci, or epigenetic marks to actively or repressively modulate gene transcription, which is genomic region dependent and associated with the distinct accumulation of RNA Pol II. Moreover, distinct types of ssDNA had differential impacts on the activities and evolution of transposable elements (TEs) (especially common or conserved TEs) in the rice genome. Our study showcases an antibody-independent technique for characterizing non-B DNA structures or functional genomic loci in plants. It lays the groundwork and fills a crucial gap for further exploration of ssDNA, non-B DNA structures, or functional genomic loci, thereby advancing our understanding of their biology in plants.

Population epigenetics: DNA methylation in the plant omics era.

DNA methylation plays an important role in many biological processes. The mechanisms underlying the establishment and maintenance of DNA methylation are well understood thanks to decades of research using DNA methylation mutants, primarily in Arabidopsis (Arabidopsis thaliana) accession Col-0. Recent genome-wide association studies (GWASs) using the methylomes of natural accessions have uncovered a complex and distinct genetic basis of variation in DNA methylation at the population level. Sequencing following bisulfite treatment has served as an excellent method for quantifying DNA methylation. Unlike studies focusing on specific accessions with reference genomes, population-scale methylome research often requires an additional round of sequencing beyond obtaining genome assemblies or genetic variations from whole-genome sequencing data, which can be cost prohibitive. Here, we provide an overview of recently developed bisulfite-free methods for quantifying methylation and cost-effective approaches for the simultaneous detection of genetic and epigenetic information. We also discuss the plasticity of DNA methylation in a specific Arabidopsis accession, the contribution of DNA methylation to plant adaptation, and the genetic determinants of variation in DNA methylation in natural populations. The recently developed technology and knowledge will greatly benefit future studies in population epigenomes.

The Arabidopsis ATR-SOG1 signaling module regulates pleiotropic developmental adjustments in response to 3'-blocked DNA repair intermediates.

Base excision repair and active DNA demethylation produce repair intermediates with DNA molecules blocked at the 3'-OH end by an aldehyde or phosphate group. However, both the physiological consequences of these accumulated single-strand DNAs break with 3'-blocked ends (DNA 3'-blocks) and the signaling pathways responding to unrepaired DNA 3'-blocks remain unclear in plants. Here, we investigated the effects of DNA 3'-blocks on plant development using the zinc finger DNA 3'-phosphoesterase (zdp) AP endonuclease2 (ape2) double mutant, in which 3'-blocking residues are poorly repaired. The accumulation of DNA 3'-blocked triggered diverse developmental defects that were dependent on the ATM and RAD3-related (ATR)-suppressor of gamma response 1 (SOG1) signaling module. SOG1 mutation rescued the developmental defects of zdp ape2 leaves by preventing cell endoreplication and promoting cell proliferation. However, SOG1 mutation caused intensive meristematic cell death in the radicle of zdp ape2 following germination, resulting in rapid termination of radicle growth. Notably, mutating FORMAMIDOPYRIMIDINE DNA GLYCOSYLASE (FPG) in zdp ape2 sog1 partially recovered its radicle growth, demonstrating that DNA 3'-blocks generated by FPG caused the meristematic defects. Surprisingly, despite lacking a functional radicle, zdp ape2 sog1 mutants compensated the lack of root growth by generating anchor roots having low levels of DNA damage response. Our results reveal dual roles of SOG1 in regulating root establishment when seeds germinate with excess DNA 3'-blocks.

Evolutionarily Conserved Principles Predict 3D Chromatin Organization.

Topologically associating domains (TADs), CTCF loop domains, and A/B compartments have been identified as important structural and functional components of 3D chromatin organization, yet the relationship between these features is not well understood. Using high-resolution Hi-C and HiChIP, we show that Drosophila chromatin is organized into domains we term compartmental domains that correspond precisely with A/B compartments at high resolution. We find that transcriptional state is a major predictor of Hi-C contact maps in several eukaryotes tested, including C. elegans and A. thaliana. Architectural proteins insulate compartmental domains by reducing interaction frequencies between neighboring regions in Drosophila, but CTCF loops do not play a distinct role in this organism. In mammals, compartmental domains exist alongside CTCF loop domains to form topological domains. The results suggest that compartmental domains are responsible for domain structure in all eukaryotes, with CTCF playing an important role in domain formation in mammals.

Sorting of mitochondrial and plastid heteroplasmy in Arabidopsis is extremely rapid and depends on MSH1 activity.

The fate of new mitochondrial and plastid mutations depends on their ability to persist and spread among the numerous organellar genome copies within a cell (heteroplasmy). The extent to which heteroplasmies are transmitted across generations or eliminated through genetic bottlenecks is not well understood in plants, in part because their low mutation rates make these variants so infrequent. Disruption of MutS Homolog 1 (MSH1), a gene involved in plant organellar DNA repair, results in numerous de novo point mutations, which we used to quantitatively track the inheritance of single nucleotide variants in mitochondrial and plastid genomes in Arabidopsis. We found that heteroplasmic sorting (the fixation or loss of a variant) was rapid for both organelles, greatly exceeding rates observed in animals. In msh1 mutants, plastid variants sorted faster than those in mitochondria and were typically fixed or lost within a single generation. Effective transmission bottleneck sizes (N) for plastids and mitochondria were N ∼ 1 and 4, respectively. Restoring MSH1 function further increased the rate of heteroplasmic sorting in mitochondria (N ∼ 1.3), potentially because of its hypothesized role in promoting gene conversion as a mechanism of DNA repair, which is expected to homogenize genome copies within a cell. Heteroplasmic sorting also favored GC base pairs. Therefore, recombinational repair and gene conversion in plant organellar genomes can potentially accelerate the elimination of heteroplasmies and bias the outcome of this sorting process.

Integration of large and diverse angiosperm DNA fragments into Asian Gnetum mitogenomes.

BACKGROUND: Horizontal gene transfer (HGT) events have rarely been reported in gymnosperms. Gnetum is a gymnosperm genus comprising 25‒35 species sympatric with angiosperms in West African, South American, and Southeast Asian rainforests. Only a single acquisition of an angiosperm mitochondrial intron has been documented to date in Asian Gnetum mitogenomes. We wanted to develop a more comprehensive understanding of frequency and fragment length distribution of such events as well as their evolutionary history in this genus. RESULTS: We sequenced and assembled mitogenomes from five Asian Gnetum species. These genomes vary remarkably in size and foreign DNA content. We identified 15 mitochondrion-derived and five plastid-derived (MTPT) foreign genes. Our phylogenetic analyses strongly indicate that these foreign genes were transferred from diverse eudicots-mostly from the Rubiaceae genus Coptosapelta and ten genera of Malpighiales. This indicates that Asian Gnetum has experienced multiple independent HGT events. Patterns of sequence evolution strongly suggest DNA-mediated transfer between mitochondria as the primary mechanism giving rise to these HGT events. Most Asian Gnetum species are lianas and often entwined with sympatric angiosperms. We therefore propose that close apposition of Gnetum and angiosperm stems presents opportunities for interspecific cell-to-cell contact through friction and wounding, leading to HGT. CONCLUSIONS: Our study reveals that multiple HGT events have resulted in massive amounts of angiosperm mitochondrial DNA integrated into Asian Gnetum mitogenomes. Gnetum and its neighboring angiosperms are often entwined with each other, possibly accounting for frequent HGT between these two phylogenetically remote lineages.