Comprehensive analysis of metal(loid)s and associated metal(loid) resistance genes in atmospheric particulate matter.

Despite global concerns about metal(loid)s in atmospheric particulate matter (PM), the presence of metal(loid) resistance genes (MRGs) in PM remains unknown. Therefore, we conducted a comprehensive investigation of the metal(loid)s and associated MRGs in PMs in two seasons (summer and winter) in Xiamen, China. According to the geoaccumulation index (Igeo), most metal(loid)s, except for V and Mn, exhibited enrichment in PM, suggesting potential anthropogenic sources. By employing Positive Matrix Factorization (PMF) model, utilizing a dataset encompassing both total and bioaccessible metal(loid)s, along with backward trajectory simulations, traffic emissions were determined to be the primary potential contributor of metal(loid)s in summer, whereas coal combustion was observed to have a dominant contribution in winter. The major contributor to the carcinogenic risk of metal(loid)s in both summer and winter was predominantly attributed to coal combustion, which serves as the main source of bioaccessible Cr. Bacterial communities within PMs showed lower diversity and network complexity in summer than in winter, with Pseudomonadales being the dominant order. Abundant MRGs, including the As(III) S-adenosylmethionine methyltransferase gene (arsM), Cu(I)-translocating P-type ATPase gene (copA), Zn(II)/Cd(II)/Pb(II)-translocating P-type ATPase gene (zntA), and Zn(II)-translocating P-type ATPase gene (ziaA), were detected within the PMs. Seasonal variations were observed for the metal(loid) concentration, bacterial community structure, and MRG abundance. The bacterial community composition and MRG abundance within PMs were primarily influenced by temperature, rather than metal(loid)s. This research offers novel perspectives on the occurrence of metal(loid)s and MRGs in PMs, thereby contributing to the control of air pollution.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow.

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers with c . 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar.

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Scots pine - panmixia and the elusive signal of genetic adaptation.

Scots pine is the foundation species of diverse forested ecosystems across Eurasia and displays remarkable ecological breadth, occurring in environments ranging from temperate rainforests to arid tundra margins. Such expansive distributions can be favored by various demographic and adaptive processes and the interactions between them. To understand the impact of neutral and selective forces on genetic structure in Scots pine, we conducted range-wide population genetic analyses on 2321 trees from 202 populations using genotyping-by-sequencing, reconstructed the recent demography of the species and examined signals of genetic adaptation. We found a high and uniform genetic diversity across the entire range (global FST 0.048), no increased genetic load in expanding populations and minor impact of the last glacial maximum on historical population sizes. Genetic-environmental associations identified only a handful of single-nucleotide polymorphisms significantly linked to environmental gradients. The results suggest that extensive gene flow is predominantly responsible for the observed genetic patterns in Scots pine. The apparent missing signal of genetic adaptation is likely attributed to the intricate genetic architecture controlling adaptation to multi-dimensional environments. The panmixia metapopulation of Scots pine offers a good study system for further exploration into how genetic adaptation and plasticity evolve under gene flow and changing environment.

Eco-friendly approach for enhancing the floatability of non-metallic components in waste printed circuit boards: Adding gutter oil during dry grinding.

Waste printed circuit boards (WPCBs) are an attractive secondary resource that is challenging to dispose of due to its complexity. Reverse flotation is an effective method to remove non-metallic particles (NMPs) to obtain metals from WPCBs. Nevertheless, the removal of NMPs is usually inadequate in the present flotation practice. Thus, to provide a clean approach to improve the removal efficiency of NMPs, the method of adding gutter oil during dry grinding process was adopted to enhance the hydrophobic sites on the surface of NMPs to improve the floatability. The surface morphology of NMPs was analyzed by SEM, the results show that the rough morphology inhibited the adhesion of bubbles, while water occupied the cracks and pores, making it challenging for collector adsorption, which result in unstable particle-bubble adhesion. The results of FTIR indicate that both NMPs and gutter oil have -CH3, -CH2, -C = O, -C-O functional groups, which promotes the adsorption of gutter oil on the surface of NMPs. The contact angle (CA) results show that the adsorption of gutter oil on the particle surface is conducive to the formation of enhanced CA. Furthermore, the flotation enhancement effect was verified by flotation kinetic experiments. The accumulated floats yield of NMPs conditioned by gutter oil during grinding is increased from 67.05% (NMPs without conditioning) to 95.02%, and the resin recovery is increased by 31.10%. It is demonstrated that dry grinding with gutter oil can strengthen the floatability of NMPs, which provides a potential approach to increase the flotation efficiency of WPCBs.

Inference of the distribution of fitness effects of mutations is affected by single nucleotide polymorphism filtering methods, sample size and population structure.

The distribution of fitness effects (DFE) of new mutations has been of interest to evolutionary biologists since the concept of mutations arose. Modern population genomic data enable us to quantify the DFE empirically, but few studies have examined how data processing, sample size and cryptic population structure might affect the accuracy of DFE inference. We used simulated and empirical data (from Arabidopsis lyrata) to show the effects of missing data filtering, sample size, number of single nucleotide polymorphisms (SNPs) and population structure on the accuracy and variance of DFE estimates. Our analyses focus on three filtering methods-downsampling, imputation and subsampling-with sample sizes of 4-100 individuals. We show that (1) the choice of missing-data treatment directly affects the estimated DFE, with downsampling performing better than imputation and subsampling; (2) the estimated DFE is less reliable in small samples (<8 individuals), and becomes unpredictable with too few SNPs (<5000, the sum of 0- and 4-fold SNPs); and (3) population structure may skew the inferred DFE towards more strongly deleterious mutations. We suggest that future studies should consider downsampling for small data sets, and use samples larger than 4 (ideally larger than 8) individuals, with more than 5000 SNPs in order to improve the robustness of DFE inference and enable comparative analyses.

Genomic clines across the species boundary between a hybrid pine and its progenitor in the eastern Tibetan Plateau.

Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems.

Evolutionary origin and establishment of a dioecious diploid-tetraploid complex.

Polyploids recurrently emerge in angiosperms, but most polyploids are likely to go extinct before establishment due to minority cytotype exclusion, which may be specifically a constraint for dioecious plants. Here we test the hypothesis that a stable sex-determination system and spatial/ecological isolation facilitate the establishment of dioecious polyploids. We determined the ploidy levels of 351 individuals from 28 populations of the dioecious species Salix polyclona, and resequenced 190 individuals of S. polyclona and related taxa for genomic diversity analyses. The ploidy survey revealed a frequency 52% of tetraploids in S. polyclona, and genomic k-mer spectra analyses suggested an autopolyploid origin for them. Comparisons of diploid male and female genomes identified a female heterogametic sex-determining factor on chromosome 15, which probably also acts in the dioecious tetraploids. Phylogenetic analyses revealed two diploid clades and a separate clade/grade of tetraploids with a distinct geographic distribution confined to western and central China, where complex mountain systems create higher levels of environmental heterogeneity. Fossil-calibrated phylogenies showed that the polyploids emerged during 7.6-2.3 million years ago, and population demographic histories largely matched the geological and climatic history of the region. Our results suggest that inheritance of the sex-determining system from the diploid progenitor as intrinsic factor and spatial isolation as extrinsic factor may have facilitated the preservation and establishment of polyploid dioecious populations.

Transcriptomic and metabolomic insights into the role of the flgK gene in the pathogenicity of Pseudomonas plecoglossicida to orange-spotted grouper (Epinephelus coioides).

Pseudomonas plecoglossicida is the pathogen responsible for visceral white spot disease in large yellow croaker (Larimichthys crocea) and orange-spotted grouper (Epinephelus coioides). Previously, RNA sequencing showed that P. plecoglossicida flgK gene expression was significantly up-regulated in orange-spotted grouper spleens during infection. To explore the role of flgK in P. plecoglossicida pathogenicity, RNA interference (RNAi) was performed to silence the P. plecoglossicida flgK gene, and the mutant (flgK-RNAi strain) with the best silencing efficiency (89.40%) was chosen for further study. Results showed that flgK gene silencing significantly attenuated P. plecoglossicida motility, adhesion, and biofilm formation. Compared to those fish infected with the wild-type strain of P. plecoglossicida, orange-spotted grouper infected with the flgK-RNAi strain showed a 55% increase in the survival rate and a one-day delay in time of first death, with fewer pathogens in the spleen and fewer white spots on the spleen surface. RNAi of flgK significantly affected the transcriptome and metabolome of the spleen in infected orange-spotted grouper. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the C-type lectin receptor signaling pathway was the most significantly changed immune-related pathway and the mitogen-activated protein kinase (MAPK) signaling pathway was related to multiple immune-related pathways. Furthermore, arginine biosynthesis and glycerophospholipid metabolism were the most significantly changed metabolism-related pathways. These findings suggest that flgK is a virulence gene of P. plecoglossicida. Furthermore, flgK appears to be involved in the regulation of motility, adhesion, and biofilm formation in P. plecoglossicida, as well as in the regulation of inflammatory and immune responses of orange-spotted grouper to P. plecoglossicida infection.

Repetitive Elements, Sequence Turnover and Cyto-Nuclear Gene Transfer in Gymnosperm Mitogenomes.

Among the three genomes in plant cells, the mitochondrial genome (mitogenome) is the least studied due to complex recombination and intergenomic transfer. In gymnosperms only ∼20 mitogenomes have been released thus far, which hinders a systematic investigation into the tempo and mode of mitochondrial DNA evolution in seed plants. Here, we report the complete mitogenome sequence of Platycladus orientalis (Cupressaceae). This mitogenome is assembled as two circular-mapping chromosomes with a size of ∼2.6 Mb and which contains 32 protein-coding genes, three rRNA and seven tRNA genes, and 1,068 RNA editing sites. Repetitive sequences, including dispersed repeats, transposable elements (TEs), and tandem repeats, made up 23% of the genome. Comparative analyses with 17 other mitogenomes representing the five gymnosperm lineages revealed a 30-fold difference in genome size, 80-fold in repetitive content, and 230-fold in substitution rate. We found dispersed repeats are highly associated with mitogenome expansion (r = 0.99), and most of them were accumulated during recent duplication events. Syntenic blocks and shared sequences between mitogenomes decay rapidly with divergence time (r = 0.53), with the exceptions of Ginkgo and Cycads which retained conserved genome structure over long evolutionary time. Our phylogenetic analysis supports a sister group relationship of Cupressophytes and Gnetophytes; both groups are unique in that they lost 8-12 protein-coding genes, of which 4-7 intact genes are likely transferred to nucleus. These two clades also show accelerated and highly variable substitution rates relative to other gymnosperms. Our study highlights the dynamic and enigmatic evolution of gymnosperm mitogenomes.

UV-B and UV-C radiation trigger both common and distinctive signal perceptions and transmissions in Pinus tabuliformis Carr.

In plants, ultraviolet (UV)-light is an important driver for growth and natural distribution, and is also a valuable tool for manipulating productivity as well as biotic interactions. Understanding of plant responses to different UV radiation is sparse, especially from a systems biology perspective and particularly for conifers. Here, we evaluated the physiological and transcriptomic responses to the short-term application of high-irradiance UV-B and UV-C waves on Pinus tabuliformis Carr., a major conifer in Northern China. By undertaking time-ordered gene coexpression network analyses and network comparisons incorporating physiological traits and gene expression variation, we uncovered communalities but also differences in P. tabuliformis responses to UV-B and UV-C. Both types of spectral bands caused a significant inhibition of photosynthesis, and conversely, the improvement of antioxidant capacity, flavonoid production and signaling pathways related to stress resistance, indicating a clear switch from predominantly primary metabolism to enhanced defensive metabolism in pine. We isolated distinct subnetworks for photoreceptor-mediated signal transduction, maximum quantum efficiency of photosystem II (Fv/Fm) regulation and flavonoid biosynthesis in response to UV-B and UV-C radiation. From these subnetworks, we further identified phototropins as potentially important elements in both UV-B and UV-C signaling and, for the first time, suggesting peptide hormones to be involved in promoting flavonoid biosynthesis against UV-B, while these hormones seem not to be implicated in the defense against UV-C exposure. The present study employed an effective strategy for disentangling the complex physiological and genetic regulatory mechanisms in a nonmodel plant species, and thus, provides a suitable reference for future functional evaluations and artificial UV-light mediated growing strategies in plant production.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance.

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

Chromosome-scale assembly and evolution of the tetraploid Salvia splendens (Lamiaceae) genome.

Polyploidization plays a key role in plant evolution, but the forces driving the fate of homoeologs in polyploid genomes, i.e., paralogs resulting from a whole-genome duplication (WGD) event, remain to be elucidated. Here, we present a chromosome-scale genome assembly of tetraploid scarlet sage (Salvia splendens), one of the most diverse ornamental plants. We found evidence for three WGD events following an older WGD event shared by most eudicots (the γ event). A comprehensive, spatiotemporal, genome-wide analysis of homoeologs from the most recent WGD unveiled expression asymmetries, which could be associated with genomic rearrangements, transposable element proximity discrepancies, coding sequence variation, selection pressure, and transcription factor binding site differences. The observed differences between homoeologs may reflect the first step toward sub- and/or neofunctionalization. This assembly provides a powerful tool for understanding WGD and gene and genome evolution and is useful in developing functional genomics and genetic engineering strategies for scarlet sage and other Lamiaceae species.

Haplotype-resolved genome assembly and allele-specific gene expression in cultivated ginger.

Ginger (Zingiber officinale) is one of the most valued spice plants worldwide; it is prized for its culinary and folk medicinal applications and is therefore of high economic and cultural importance. Here, we present a haplotype-resolved, chromosome-scale assembly for diploid ginger anchored to 11 pseudochromosome pairs with a total length of 3.1 Gb. Remarkable structural variation was identified between haplotypes, and two inversions larger than 15 Mb on chromosome 4 may be associated with ginger infertility. We performed a comprehensive, spatiotemporal, genome-wide analysis of allelic expression patterns, revealing that most alleles are coordinately expressed. The alleles that exhibited the largest differences in expression showed closer proximity to transposable elements, greater coding sequence divergence, more relaxed selection pressure, and more transcription factor binding site differences. We also predicted the transcription factors potentially regulating 6-gingerol biosynthesis. Our allele-aware assembly provides a powerful platform for future functional genomics, molecular breeding, and genome editing in ginger.

Landscape Setting Drives the Microbial Eukaryotic Community Structure in Four Swedish Mountain Lakes over the Holocene.

On the annual and interannual scales, lake microbial communities are known to be heavily influenced by environmental conditions both in the lake and in its terrestrial surroundings. However, the influence of landscape setting and environmental change on shaping these communities over a longer (millennial) timescale is rarely studied. Here, we applied an 18S metabarcoding approach to DNA preserved in Holocene sediment records from two pairs of co-located Swedish mountain lakes. Our data revealed that the microbial eukaryotic communities were strongly influenced by catchment characteristics rather than location. More precisely, the microbial communities from the two bedrock lakes were largely dominated by unclassified Alveolata, while the peatland lakes showed a more diverse microbial community, with Ciliophora, Chlorophyta and Chytrids among the more predominant groups. Furthermore, for the two bedrock-dominated lakes-where the oldest DNA samples are dated to only a few hundred years after the lake formation-certain Alveolata, Chlorophytes, Stramenopiles and Rhizaria taxa were found prevalent throughout all the sediment profiles. Our work highlights the importance of species sorting due to landscape setting and the persistence of microbial eukaryotic diversity over millennial timescales in shaping modern lake microbial communities.

Divergent patterns between phenotypic and genetic variation in Scots pine.

In boreal forests, autumn frost tolerance in seedlings is a critical fitness component because it determines survival rates during regeneration. To understand the forces that drive local adaptation in this trait, we conducted freezing tests in a common garden setting for 54 Pinus sylvestris (Scots pine) populations (>5000 seedlings) collected across Scandinavia into western Russia, and genotyped 24 of these populations (>900 seedlings) at >10 000 SNPs. Variation in cold hardiness among populations, as measured by QST , was above 80% and followed a distinct cline along latitude and longitude, demonstrating significant adaptation to climate at origin. In contrast, the genetic differentiation was very weak (mean FST 0.37%). Despite even allele frequency distribution in the vast majority of SNPs among all populations, a few rare alleles appeared at very high or at fixation in marginal populations restricted to northwestern Fennoscandia. Genotype-environment associations showed that climate variables explained 2.9% of the genetic differentiation, while genotype-phenotype associations revealed a high marker-estimated heritability of frost hardiness of 0.56, but identified no major loci. Very extensive gene flow, strong local adaptation, and signals of complex demographic history across markers are interesting topics of forthcoming studies on this species to better clarify signatures of selection and demography.

UV-B-induced molecular mechanisms of stress physiology responses in the major northern Chinese conifer Pinus tabuliformis Carr.

During their lifetimes, plants are exposed to different abiotic stress factors eliciting various physiological responses and triggering important defense processes. For UV-B radiation responses in forest trees, the genetics and molecular regulation remain to be elucidated. Here, we exposed Pinus tabuliformis Carr., a major conifer from northern China, to short-term high-intensity UV-B and employed a systems biology approach to characterize the early physiological processes and the hierarchical gene regulation, which revealed a temporal transition from primary to secondary metabolism, the buildup of enhanced antioxidant capacity and stress-signaling activation. Our findings showed that photosynthesis and biosynthesis of photosynthetic pigments were inhibited, while flavonoids and their related derivates biosynthesis, as well as glutathione and glutathione S-transferase mediated antioxidant processes, were enhanced. Likewise, stress-related phytohormones (jasmonic acid, salicylic acid and ethylene), kinase and reactive oxygen species signal transduction pathways were activated. Biological processes regulated by auxin and karrikin were, for the first time, found to be involved in plant defense against UV-B by promoting the biosynthesis of flavonoids and the improvement of antioxidant capacity in our research system. Our work evaluated the physiological and transcriptome perturbations in a conifer's response to UV-B, and generally, highlighted the necessity of a systems biology approach in addressing plant stress biology.

Chromosome-level genome assembly of a parent species of widely cultivated azaleas.

Azaleas (Ericaceae) comprise one of the most diverse ornamental plants, renowned for their cultural and economic importance. We present a chromosome-scale genome assembly for Rhododendron simsii, the primary ancestor of azalea cultivars. Genome analyses unveil the remnants of an ancient whole-genome duplication preceding the radiation of most Ericaceae, likely contributing to the genomic architecture of flowering time. Small-scale gene duplications contribute to the expansion of gene families involved in azalea pigment biosynthesis. We reconstruct entire metabolic pathways for anthocyanins and carotenoids and their potential regulatory networks by detailed analysis of time-ordered gene co-expression networks. MYB, bHLH, and WD40 transcription factors may collectively regulate anthocyanin accumulation in R. simsii, particularly at the initial stages of flower coloration, and with WRKY transcription factors controlling progressive flower coloring at later stages. This work provides a cornerstone for understanding the underlying genetics governing flower timing and coloration and could accelerate selective breeding in azalea.

Microhomologies Are Associated with Tandem Duplications and Structural Variation in Plant Mitochondrial Genomes.

Short tandem repeats (STRs) contribute to structural variation in plant mitochondrial genomes, but the mechanisms underlying their formation and expansion are unclear. In this study, we detected high polymorphism in the nad7-1 region of the Pinus tabuliformis mitogenome caused by the rapid accumulation of STRs and rearrangements over a few million years ago. The STRs in nad7-1 have a 7-bp microhomology (TAG7) flanking the repeat array. We then scanned the mitogenomes of 136 seed plants to understand the role of microhomology in the formation of STR and mitogenome evolution. A total of 13,170 STRs were identified, and almost half of them were associated with microhomologies. A substantial amount (1,197) of microhomologies was long enough to mediate structural variation, and the length of microhomology is positively correlated with the length of tandem repeat unit. These results suggest that microhomology may be involved in the formation of tandem repeat via microhomology-mediated pathway, and the formation of longer duplicates required greater length of microhomology. We examined the abundance of these 1,197 microhomologies, and found 75% of them were enriched in the plant mitogenomes. Further analyses of the 400 prevalent microhomologies revealed that 175 of them showed differential enrichment between angiosperms and gymnosperms and 186 differed between angiosperms and conifers, indicating lineage-specific usage and expansion of microhomologies. Our study sheds light on the sources of structural variation in plant mitochondrial genomes and highlights the importance of microhomology in mitochondrial genome evolution.