Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Copy number variation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources
PubMed2004

A versatile statistical analysis algorithm to detect genome copy number variation.

We have developed a versatile statistical analysis algorithm for the detection of genomic aberrations in human cancer cell lines. The algorithm analyzes genomic data obtained from a variety of array technologies, such as oligonucleotide array, bacterial artificial chromosome array, or array-based comparative genomic hybridization, that operate by hybridizing with genomic material obtained from cancer and normal cells and allow detection of regions of the genome with altered copy number. The number of probes (i.e., resolution), the amount of uncharacterized noise per probe, and the severity of chromosomal aberrations per chromosomal region may vary with the underlying technology, biological sample, and sample preparation. Constrained by these uncertainties, our algorithm aims at robustness by using a priorless maximum a posteriori estimator and at efficiency by a dynamic programming implementation. We illustrate these characteristics of our algorithm by applying it to data obtained from representational oligonucleotide microarray analysis and array-based comparative genomic hybridization technology as well as to synthetic data obtained from an artificial model whose properties can be varied computationally. The algorithm can combine data from multiple sources and thus facilitate the discovery of genes and markers important in cancer, as well as the discovery of loci important in inherited genetic disease.

Algorithms↗
PubMed1993

Variant mapping of the Apo(B) AT rich minisatellite. Dependence on nucleotide sequence of the copy number variations. Instability of the non-canonical alleles.

Because of its variations in length, the AT rich Hyper-Variable Region (HVR) of the 3' end of the Apolipoprotein B gene is used as a polymorphic maker in genetic studies. It contains a SspI site in its repeated motif and we used this feature to precisely analyse the internal structure of the different alleles found at this locus in a Caucasian population. We performed total digestion on 194 alleles as well as Minisatellite Variant Repeat mapping (MVR mapping: partial digestion) on 54. The results show that the level of length variability (in copy number) of the 5' end of this locus is at least two times higher than that of the 3' end. This could be correlated with the difference in nucleotide sequence between the two parts of the HVR and suggests the dependence on the primary structure of the mechanism that produces length variability. A molecular model is proposed to explain this result. Moreover, the sharp analysis of the minisatellite structure by the distribution of SspI sites reveals differences between long and short alleles, indicating that in most cases, no recombination occurs between alleles of different sizes. Finally the rare alleles exhibit a non-canonical structure. These important points could explain the bimodal distribution of the frequencies of the alleles in the population.

Alleles↗
PubMed2006

The mosaic nature of intergenic 16S-23S rRNA spacer regions suggests rRNA operon copy number variation in Clostridium difficile strains.

Clostridium difficile is a major spore-forming environmental pathogen that causes serious health problems in patients undergoing antibiotic therapy. Consequently, reliable and sensitive methods for typing individual strains are required for epidemiological and environmental studies. Ribotyping is generally considered the best method, but it fails to account for sequence diversity which might exist in intergenic 16S-23S rRNA spacer regions (ISRs) within and among strains of this organism. Therefore, this study was undertaken to compare the sequence of each individual ISR in five strains of C. difficile to explore the extent of this diversity and see whether such information might provide the basis for more sensitive and discriminatory strain typing methods. After targeted PCR amplification, cloning, and sequencing, the diversity of the ISRs was used as a measure of rRNA operon copy number. In C. difficile strains 630, ATCC 43593, A, and B, 11, 11, 7, and 8 ISR length variants, respectively, were found (containing different combinations of sequence groups [i to xiii]), suggesting 11, 11, 7, and 8 rrn copies in the respective strains. Many ISRs of the same length differed markedly in their sequences, and some of these were restricted in occurrence to a single strain. Most of these ISRs did not contain any tRNA genes, and only single copies of the tRNA(Ala) gene were found in those that did. The presence of ISR sequence groups (i to xiii) varied between strains, with some found in one, two, three, four, or all five strains. We conclude that the intergenic 16S-23S rRNA spacer regions showed a high degree of diversity, not only among the rrn operons in different strains and different rrn copies in a single strain but also among ISRs of the same length. It appears that C. difficile ISRs vary more at the inter- and intragenic levels than those of other species as determined by empirical comparison of sequences. The precise characterization of these sequences has demonstrated a high level of mosaic sequence block rearrangements that are present or absent in multiple strain-variable rrn copies within and between five different strains of C. difficile.

Bacterial Typing Techniques↗
PubMed2006

Global variation in copy number in the human genome.

Copy number variation (CNV) of DNA sequences is functionally significant but has yet to be fully ascertained. We have constructed a first-generation CNV map of the human genome through the study of 270 individuals from four populations with ancestry in Europe, Africa or Asia (the HapMap collection). DNA from these individuals was screened for CNV using two complementary technologies: single-nucleotide polymorphism (SNP) genotyping arrays, and clone-based comparative genomic hybridization. A total of 1,447 copy number variable regions (CNVRs), which can encompass overlapping or adjacent gains or losses, covering 360 megabases (12% of the genome) were identified in these populations. These CNVRs contained hundreds of genes, disease loci, functional elements and segmental duplications. Notably, the CNVRs encompassed more nucleotide content per genome than SNPs, underscoring the importance of CNV in genetic diversity and evolution. The data obtained delineate linkage disequilibrium patterns for many CNVs, and reveal marked variation in copy number among populations. We also demonstrate the utility of this resource for genetic disease studies.

Chromosome Mapping↗
PubMed2026

Dual-dimensional profiling of host genomic variations and HPV integration in PD-L1-stratified cervical cancer via Oxford Nanopore Technology.

BACKGROUND: The integration of human papillomavirus (HPV) DNA into the host genome is a key step in the development of HPV-associated cervical cancer (CC). However, the genomic characteristics of host genomic variations and HPV integration within the context of programmed death-ligand 1 (PD-L1) expression stratification have not been systematically investigated. METHODS: Whole-genome sequencing was performed using Oxford Nanopore Technology (ONT) on six samples (three from the high PD-L1 expression group and three from the low PD-L1 expression group). The characteristics of host genomic variations under different PD-L1 expression stratifications were explored, including structural variations (SV), copy number variations (CNV), single nucleotide polymorphisms (SNP), and insertion-deletions (Indel). Subsequently, the distribution features of HPV integration sites were analyzed, different integration types were identified, and pathway analysis was conducted. RESULTS: Whole-genome SV analysis revealed that the total number of SVs and the composition of mutation types were similar between the high and low PD-L1 expression groups, with insertions (INS) and deletions (DEL) predominating in both. These variations were primarily enriched in intergenic regions and introns. In the low PD-L1 expression group, integration events were observed at multiple chromosomal loci, with the most frequent integration occurring in the KLF5 gene region on chromosome 13. No frequently integrated loci were identified in the high PD-L1 expression group. Additionally, four distinct HPV integration breakpoint patterns were preliminarily identified and analyzed. CONCLUSION: PD-L1 expression stratification did not significantly alter the overall genomic instability of the host. However, differences were observed in the distribution patterns of HPV integration sites. These findings provide new insights into the genomic heterogeneity of CC under different PD-L1 expression backgrounds and may lay the groundwork for future research exploring stratified immunotherapy based on HPV integration features.

Humans↗
PubMed2003

A sensitive method for detecting variation in copy numbers of duplicated genes.

Gene duplications are common in the vertebrate genome, and duplicated loci often show a variation in copy number that may have important phenotypic effects. Here we describe a powerful method for quantification of duplicated copies based on pyrosequencing. A reliable quantification was obtained by amplification of the duplication break-point and a corresponding nonduplicated sequence in a competitive PCR assay. A comparison with an independent method for quantification based on the Invader technology revealed an excellent correlation between the two methods. The pyrosequencing-based method was evaluated by analyzing variation in copy number at the duplicated KIT/Dominant white locus in pigs. We were able to distinguish haplotypes at this locus by combining the duplication breakpoint test with a diagnostic test for a functionally important splice mutation in the duplicated gene. An extensive allelic variation, including the presence of a new allele carrying a single KIT copy expected to encode a truncated KIT receptor, was revealed when analyzing white pigs from commercial lines.

Algorithms↗
PubMed1998

Simultaneous estimation of all the parameters of a stepwise mutation model.

Minisatellite and microsatellite are short tandemly repetitive sequences dispersed in eukaryotic genomes, many of which are highly polymorphic due to copy number variation of the repeats. Because mutation changes copy numbers of the repeat sequences in a generalized stepwise fashion, stepwise mutation models are widely used for studying the dynamics of these loci. We propose a minimum chi-square (MCS) method for simultaneous estimation of all the parameters in a stepwise mutation model and the ancestral allelic type of a sample. The MCS estimator requires knowing the mean number of alleles of a certain size in a sample, which can be estimated using Monte Carlo samples generated by a coalescent algorithm. The method is applied to samples of seven (CA)n repeat loci from eight human populations and one chimpanzee population. The estimated values of parameters suggest that there is a general tendency for microsatellite alleles to expand in size, because (1) each mutation has a slight tendency to cause size increase and (2) the mean size increase is larger than the mean size decrease for a mutation. Our estimates also suggest that most of these CA-repeat loci evolve according to multistep mutation models rather than single-step mutation models. We also introduced several quantities for measuring the quality of the estimation of ancestral allelic type, and it appears that the majority of the estimated ancestral allelic types are reasonably accurate. Implications of our analysis and potential extensions of the method are discussed. SINCE the discovery that a large number of loci with tandemly repeated sequences in human and many eukaryote species are highly polymorphic because of copy number variation of the repeats in different individuals (Jeffreys 1985; Litt and Luty 1989; Weber and May 1989), allele size data from such loci are rapidly becoming the dominant source of genetic markers for genome mapping, forensic testing, and population studies. Loci with repeat sequences longer than 5 bp are generally referred to as minisatellite or variable number tandem repeat loci, and those with repeat sequences between 2 to 5 bp are referred to as microsatellite or short tandem repeat loci (Tautz 1993). Because mutations change the copy number of such loci in a stepwise fashion, rapid accumulation of population samples from minisatellite and microsatellite loci has resurrected the interest of the stepwise mutation model (SMM), which was popular in the 1970s.

Chi-Square Distribution↗
PubMed2006

Increased CYP3A4 copy number in TONG/HCC cells but not in DNA from other humans.

Two recent screens for copy-number variations in the entire human genome found 12.4 gene copy number variations per person, including 2.5% of individuals with gains between 7q21.1 and 7q22.1, the chromosomal location of CYP3A4. CYP3A4 is involved in the metabolism of approximately 50% of all drugs, including many cancer chemotherapeutic agents. CYP3A4 gene copy was determined in DNA from 143 individuals: normal human livers, primary and secondary liver tumors, human hepatic cell lines, and immortalized cell lines representing eight ethnically diverse populations. CYP3A4 gene copy was normal in all but one sample, a primary human hepatocellular carcinoma cell line (TONG/HCC). Southern blots of TONG/HCC DNA revealed an approximate 10-fold increase in CYP3A and a corresponding increase in CYP3A mRNA expression and catalytic activity. Fluorescent in situ hybridization of TONG/HCC revealed specific amplification of the CYP3A4 gene on chromosome 7q21 but no amplification of the MDR1 gene that localizes 11.9 Mb upstream of CYP3A4. High resolution analysis of DNA copy number by comparative genomic hybridization confirmed amplification at 7q21.3-7q22. The amplicon spanned 1.7 Mb and contained 30 known genes, including the entire CYP3A locus. To determine whether CYP3A4 expression affected chemotherapeutic toxicity, LLC-PK1 cells were transduced with adenoviruses expressing CYP3A4 and P450 reductase. CYP3A4 conferred resistance to taxol, vinblastine and topotecan. These studies demonstrate that CYP3A4 copy number differences do not contribute to the normal variation in CYP3A4 expression. Tumors with increased CYP3A copy number (via amplification or increased chromosome 7q) would be expected to show reduced cytotoxicity to some chemotherapeutic drugs and potentially an increase in the outgrowth of drug resistant tumors.

Blotting, Southern↗
PubMed2026

Likelihood-based optimization enables accurate copy number estimation for paralogous genes using exome data.

MOTIVATION: Exome sequencing is widely used for genetic studies; however, accurate detection of copy number variants (CNV) in paralogous genes is challenging due to short-read mapping ambiguity and extensive copy-number variation. The human genome contains several hundred paralogous genes, many of which are known to harbor disease-associated CNVs. Existing exome CNV callers are primarily designed for rare CNV detection in uniquely mappable regions and are not well-suited for paralogous genes. METHODS: We describe a computational method (EdgeCopy) for copy number profiling of paralogous genes using whole-exome sequence data. EdgeCopy aggregates reads mapped to all copies of paralogous genes and relates observed read depth to copy number for multiple exome samples using an approximate composite likelihood function. The likelihood function is optimized using numerical optimization to obtain gene-level fractional copy number estimates that are discretized and refined using a Hidden Markov Model to obtain exon-level copy number estimates. RESULTS: Benchmarking of Edgecopy using experimental copy number data showed high concordance (mean = 0.973) for six disease-associated paralogous genes. We evaluated performance using whole-exome data from approximately 2400 samples across five continental populations from the 1000 Genomes Project. EdgeCopy shows robust concordance with whole-genome sequencing based estimates (0.974-0.982) across populations and 130 paralogous genes spanning a wide range of copy-number variation. In comparison, copy number analysis using a state-of-the-art exome CNV caller failed to estimate copy number for paralogous genes with very high mapping ambiguity and showed much lower concordance (0.565) for CNV events compared to EdgeCopy (0.908). AVAILABILITY: EdgeCopy is freely available at https://github.com/vibansal-lab/edgecopy.

Humans↗
PubMed1999

Genome-wide analysis of DNA copy-number changes using cDNA microarrays.

Gene amplifications and deletions frequently contribute to tumorigenesis. Characterization of these DNA copy-number changes is important for both the basic understanding of cancer and its diagnosis. Comparative genomic hybridization (CGH) was developed to survey DNA copy-number variations across a whole genome. With CGH, differentially labelled test and reference genomic DNAs are co-hybridized to normal metaphase chromosomes, and fluorescence ratios along the length of chromosomes provide a cytogenetic representation of DNA copy-number variation. CGH, however, has a limited ( approximately 20 Mb) mapping resolution, and higher-resolution techniques, such as fluorescence in situ hybridization (FISH), are prohibitively labour-intensive on a genomic scale. Array-based CGH, in which fluorescence ratios at arrayed DNA elements provide a locus-by-locus measure of DNA copy-number variation, represents another means of achieving increased mapping resolution. Published array CGH methods have relied on large genomic clone (for example BAC) array targets and have covered only a small fraction of the human genome. cDNAs representing over 30,000 radiation-hybrid (RH)-mapped human genes provide an alternative and readily available genomic resource for mapping DNA copy-number changes. Although cDNA microarrays have been used extensively to characterize variation in human gene expression, human genomic DNA is a far more complex mixture than the mRNA representation of human cells. Therefore, analysis of DNA copy-number variation using cDNA microarrays would require a sensitivity of detection an order of magnitude greater than has been routinely reported. We describe here a cDNA microarray-based CGH method, and its application to DNA copy-number variation analysis in breast cancer cell lines and tumours. Using this assay, we were able to identify gene amplifications and deletions genome-wide and with high resolution, and compare alterations in DNA copy number and gene expression.

Chromosomes, Human, Pair 17↗
PubMed1998

Trade-off between segregational stability and metabolic burden: a mathematical model of plasmid ColE1 replication control.

A model of ColE1 copy number control has been developed where molecular details of replication are connected both to segregational stability and metabolic burden. Efficient replication control reduces copy number variation and increases segregational stability for a given average copy number. Copy number variation is predicted to depend on the type of inhibition mechanism as well as RNA I and RNA II turnover rate constants. It is shown that when both RNA I and RNA II transcription frequencies and the rate constant for degradation of free RNA I are very large, a hyperbolic inhibition mechanism must compensate with a 1.4 times greater average copy number to obtain the same segregational stability as an exponential inhibition mechanism. How sensitively the replication frequency responds to changes in RNA I concentration depends on the type of inhibition mechanism and the number of attempts to form an RNA II replication primer per plasmid and cell cycle. If RNA I is too stable, it will not follow changes in plasmid concentration closely, and when the transcription frequency for RNA I is only slightly higher than for RNA II, RNA I concentration becomes randomized. In both these cases, the proportionality between the single cell RNA I and plasmid concentrations is lost and this impairs copy number control. Thresholds in the rate for degradation of free RNA I as well as in RNA I and RNA II transcription frequencies have been computed, where an increase in these rate constants has a negligible effect on segregational stability but a corresponding decrease leads to segregational disaster. This indicates that there exists a well defined optimal set of rate constants where the regulation system works well without excessive metabolic load. A number of new experiments are suggested to address features of particular importance for the evolution of ColE1 copy number control.

DNA Replication↗
PubMed2002

Intratumoral regional variations in copy number of the chromosomal part revealed by microdissection and combined ploidy and comparative genomic hybridization analyses in esophageal squamous cell carcinoma.

Intratumoral regional variations in the copy number of chromosomal material were analyzed to demonstrate the time sequence of chromosomal changes in progression of individual squamous cell carcinoma of the esophagus. We applied combined DNA ploidy and comparative genomic hybridization (CGH) analyses to multiple DNA samples extracted from microdissected, formalin-fixed, paraffin-embedded tissues, and amplified and labeled according to degenerate oligonucleotide-primed polymerase chain reaction. We examined two cases: one with a deep invasive tumor and the other with a superficial spreading tumor. We found that each sample had unique aberrations in addition to the ones common to all or some of the samples in a tumor. Based on previous studies (Okada et al., Cancer Genet Cytogenet 2000;118:99-107), we classified significant shifts of the green to red (G/R) ratio into small and large ones, which were within and beyond the range of 0.65 to 1.35, respectively. Most of the large-shift aberrations were found to be common to all or some of the samples in each case. These were thought to represent earlier events in the DNA-diploid stage, while small shifts may possibly reflect one-copy changes after tetraploidization or chromosomal instability. Based on the breakpoints and on the absolute copy numbers of altered chromosomal parts inferred from DNA ploidy and the shift size of the G/R ratio, we reconstructed the sequence of accumulation and divergence of chromosomal alterations as a dendrogram in each case. This method of temporal analysis may enable us to extract important early events from numerous aberrations screened by CGH in individual tumors.

Anatomy, Regional↗
PubMed1998

Requirements for rapid plasmid ColE1 copy number adjustments: a mathematical model of inhibition modes and RNA turnover rates.

The random distribution of ColE1 plasmids between the daughter cells at cell division introduces large copy number variations. Statistic variation associated with limited copy number in single cells also causes fluctuations to emerge spontaneously during the cell cycle. Efficient replication control out of steady state is therefore important to tame such stochastic effects of small numbers. In the present model, the dynamic features of copy number control are divided into two parts: first, how sharply the replication frequency per plasmid responds to changes in the concentration of the plasmid-coded inhibitor, RNA I, and second, how tightly RNA I and plasmid concentrations are coupled. Single (hyperbolic)- and multiple (exponential)-step inhibition mechanisms are compared out of steady state and it is shown how the response in replication frequency depends on the mode of inhibition. For both mechanisms, sensitivity of inhibition is "bought" at the expense of a rapid turnover of a replication preprimer, RNA II. Conventional, single-step, inhibition kinetics gives a sloppy replication control even at high RNA II turnover rates, whereas multiple-step inhibition has the potential of working with unlimited precision. When plasmid concentration changes rapidly, RNA I must be degraded rapidly to be "up to date" with the change. Adjustment to steady state is drastically impaired when the turnover rate constants of RNA I decrease below certain thresholds, but is basically unaffected for a corresponding increase. Several features of copy number control that are shown to be crucial for the understanding of ColE1-type plasmids still remain to be experimentally characterized. It is shown how steady-state properties reflect dynamics at the heart of regulation and therefore can be used to discriminate between fundamentally different copy number control mechanisms. The experimental tests of the predictions made require carefully planned assays, and some suggestions for suitable experiments arise naturally from the present work. It is also discussed how the presence of the Rom protein may affect dynamic qualities of copy number control.

Bacteriocin Plasmids↗
PubMed1999

Molecular characterization of ribosomal gene variation within and among NORs segregating in specialized populations of chicken.

The molecular organization of the 18S, 5.8S, and 28S ribosomal RNA gene repeat units, located at the single nucleolus organizer region (NOR) locus in the chicken, was investigated in genetically distinct populations of research and commercial chickens. Substantial gene repeat variation within and among NORs was documented. Intact ribosomal gene repeat size ranged from 11 kb to over 50 kb. Unique combinations of ribosomal genes, of different size, were specific to particular populations. It was determined that the basis for the ribosomal gene repeat size variation was intergenic spacer (IGS) length heterogeneity. Interestingly, in different populations, the location of the variation that contributes to length heterogeneity was specific to particular IGS subregions. In addition to IGS variation, an inbred line of Red Jungle Fowl exhibited coding region variation. Ribosomal gene copy number variation was also studied, and line averages ranged from 279 to 368. Average rDNA array size (a function of copy number and gene repeat length) was calculated for each of the populations and found to vary over a range of two megabases, from 5 to 7 Mb.

Animals↗
PubMed1998

Extraordinary ribosomal spacer length heterogeneity in a neotyphodium endophyte hybrid: implications for concerted evolution.

An extraordinary level of length heterogeneity was found in the ribosomal DNA (rDNA) of an asexual hybrid Neotyphodium grass endophyte, isolate Lp1. This hybrid Neotyphodium endophyte is an interspecific hybrid between two grass endophytes, Neotyphodium lolii, and a sexual form, Epichlöe typhina, and the length heterogeneity was not found in either of these progenitor species. The length heterogeneity in the hybrid is localized to the intergenic spacer (IGS) and is the result of copy-number variation of a tandemly repeated subrepeat class within the IGS, the 111-/119-bp subrepeats. Copy number variation of this subrepeat class appears to be a consequence of mitotic unequal crossing over that occurs between these subrepeats. This implies that unequal crossing over plays a role in the concerted evolution of the whole rDNA. Changes in the pattern of IGS length variants occurred in just two rounds of single-spore purification. Analysis of the IGS length heterogeneity revealed features that are unexpected in a simple model of unequal crossing over. Potential refinements of the molecular details of unequal crossing over are presented, and we also discuss evidence for a combination of homogenization mechanisms that drive the concerted evolution of the Lp1 rDNA.

Acremonium↗
PubMed2002

DNA copy number analysis by MAPH: molecular diagnostic applications.

DNA copy number variation is an important cause of genetic disease. There are several techniques available to detect copy number changes of various sizes, each with their limitations in resolution and cost. Here we outline the development of multiplex amplifiable probe hybridization (MAPH) into a high-throughput diagnostic technique for detecting copy number variation of almost any size. Its application in testing for genetic mutations causing diseases, such as familial breast cancer, Charcot-Marie-Tooth disease Type 1A, Duchenne/Becker muscular dystrophy and familial colorectal cancer is described, as well as its use in identifying chromosomal changes in some individuals with mental retardation. The analysis of the data produced by MAPH is also considered, along with its potential for automation and development of microarray-based MAPH.

Animals↗
PubMed2025

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera.

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Genome, Mitochondrial↗