PubMed HealthSearch

SEARCH · PubMed Health

Results for “Purifying selection”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 37 records · Page 2Linked to original sources

Are rare variants responsible for susceptibility to complex diseases?

Little is known about the nature of genetic variation underlying complex diseases in humans. One popular view proposes that mapping efforts should focus on identification of susceptibility mutations that are relatively old and at high frequency. It is generally assumed-at least for modeling purposes-that selection against complex disease mutations is so weak that it can be ignored. In this article, I propose an explicit model for the evolution of complex disease loci, incorporating mutation, random genetic drift, and the possibility of purifying selection against susceptibility mutations. I show that, for the most plausible range of mutation rates, neutral susceptibility alleles are unlikely to be at intermediate frequencies and contribute little to the overall genetic variance for the disease. Instead, it seems likely that the bulk of genetic variance underlying diseases is due to loci where susceptibility mutations are mildly deleterious and where there is a high overall mutation rate to the susceptible class. At such loci, the total frequency of susceptibility mutations may be quite high, but there is likely to be extensive allelic heterogeneity at many of these loci. I discuss some practical implications of these results for gene mapping efforts.

Alleles

Cell-mediated cytotoxicity. Characterization of the effector cells.

Isolated human mononuclear cells were fractionated according to their membrane characteristics or physical properties. Adherent cells were depleted by filtration through glass columns; phagocytic cells were removed by iron treatment and cell subpopulations capable of forming rosettes with sheep erythrocytes (E), erythrocyte-antibody-complement (EAC) and chicken erythrocyte-antibody complexes (CEA) were separated by centrifugation of Ficoll-Hypaque gradients. The functional activity of the cell subpopulations obtained was assayed by testing PHA-induced cytoxicity (PIC), antibody-dependent cytoxicity (ADCC) and blast transformation by PHA. The results of this study demonstrate that: (1) cells reacting in PIC and ADCC assays are different, adherent and phagocytic cells being necessary for full expression of PIC and not for ADCC; (2) PHA induces direct blast transformation of purified E-RFC in the absence of PIC cytotoxic cells; (3) cell populations specifically enriched in E or EAC rosette-forming cells are not cytotoxic neither in the PHA nor in antibody mediated cytotoxic assays; (4) cells participating in ADCC can be selectively purified by centrifugation of CEA rosettes.

Antibodies

Genetic and environmental interactions outweigh mitonuclear coevolution for complex traits in Drosophila.

The interdependent relationship between mitochondrial and nuclear genomes is a powerful model for understanding how epistasis shapes the architecture and evolution of complex traits. Once considered a neutral marker, mitochondrial DNA variation is now recognized as critical to phenotypic evolution because of its epistatic interactions and history of coevolution with the nuclear genome. A central challenge in evolutionary genetics is to quantify the relative importance of stabilizing and directional selection shaping complex trait distributions within and among species. Both can act on interacting and/or co-evolving genes contributing to quantitative traits, but resolving their relative roles is complicated by the complex architecture of most traits. Here, we use a panel of 90 Drosophila mitonuclear genotypes to quantify the relative contributions of mitochondrial, nuclear, and environmental variation and their interactions to four metabolically demanding complex traits. We sample both within-species and between-species mitochondrial variation and observe stronger interaction effects attributable to within-species variation, consistent with stabilizing selection maintaining mitonuclear function. Additionally, culturing the flies on a mitochondrial Complex I inhibitor, rotenone, reveals significant genotype x environment (G×E and G×G×E) interaction effects, providing insight into how genetic variation can be maintained across changing environments. Our results have broader implications in medicine, where mitochondrial DNA donors with longer purifying selection histories may be safer for mitochondrial replacement therapies.

Journal Article

Sexual selection, genomic evolution and population fitness in Drosophila pseudoobscura.

Sexual selection shapes the genome in unique ways. It is also likely to have significant fitness consequences, such as purging deleterious mutations from the genome or conversely maintaining genetic load in a population via sexual conflict. Here, we examined what the influence of sexual selection has on genomic variation potentially underlying population fitness using experimentally evolved Drosophila pseudoobscura populations. Sexual selection was manipulated by keeping replicate lines in elevated polyandry or strict monogamy for approximately 200 generations followed by individual-based sequencing. Using pi (π), fixation index (Fst)and recombination rate measures, we confirmed signatures of selection were not dispersed but mainly localized to the third and X chromosome. Overall mutational load was similar between lines but our analysis of the distribution of fitness effects revealed considerable variation between lines and chromosomes. Furthermore, we found that the distribution of transposable elements differs between the lines, with a higher load in monogamous lines. Our results suggest that complex interactions between purifying selection and sexual conflict are shaping the genome, particularly on chromosome 3 and the sex chromosome; sexual selection influences divergence across chromosomes but in a more complex way than proposed by simple 'purging' of deleterious loci.

Animals

SARS-CoV-2 intra-host variation shows evidence of transmission and convergent evolution in a university surveillance cohort.

Monitoring and understanding the transmission and evolution of SARS-CoV-2 remains a significant public health priority. Within-host genetic variation provides insight into viral evolution during infection and may help infer transmission events. In this study, we analysed intra-host variation in SARS-CoV-2 genome sequences from Boston University's testing mandate. Focusing on intra-host single nucleotide variants (iSNVs), we inferred transmission events and assessed the selective forces shaping within-host viral evolution. To minimize false-positive iSNVs resulting from systematic biases, we implemented stringent data filtering and developed a heuristic to exclude contamination-derived artefacts arising from batched sequencing. We find that intra-host variation is limited and infrequently transmitted during acute infections, suggesting that shared iSNVs serve as highly specific but insensitive markers of transmission. We also observed incomplete purifying selection shaping within-host diversity, with the loci most affected changing among variants of concern. Finally, we identified a highly recurrent iSNV (G11083T) which may represent a site of positive selection. Our results highlight that within-host variation provides insight into within-host pathogen evolution, in spite of its limited use in genomic epidemiology.

SARS-CoV-2

Transgenerational continuity: Persistence as a dimension of inheritance and evolution.

Transgenerational continuity (TC) describes the persistence of inherited molecular architectures across generations. Progress in identity-by-descent (IBD) detection, recombination dynamics, and epigenetic research highlights the growing need for a more comprehensive model of inheritance. This theoretical framework synthesizes evidence from genomics, population studies, and epigenetics to outline how inherited molecular architectures, which are transmitted through IBD, together with heritable epigenetic modifications, can preserve ancestral information across generations. IBD captures genomic continuity across three nested scales, where recent familial segments link close relatives, population-level haplotypes are shared across cohorts, and archaic fragments from Neanderthal and Denisovan admixture persist as molecular fossils of ancient lineages. Although recombination and selection reshape these regions, their persistence across time scales highlights the evolutionary durability of genomic continuity. Epigenetic memory reflects regulatory persistence, whereby molecular modifications can preserve functional states across cell divisions and sometimes across generations. Together with familial and population-level IBD persistence and the long-term retention of introgressed haplotypes, these findings demonstrate that inherited molecular architectures can persist across multiple timescales. Evolutionary processes shape this persistence. Purifying selection preferentially removes deleterious inherited variants, whereas positive selection can favor the persistence of functionally relevant genomic architectures. From this perspective, evolutionary dynamics arise not only from the generation of variation, but also from the differential persistence of inherited molecular architectures through selection. Transgenerational continuity therefore provides a conceptual framework in which persistence serves as an explanatory dimension of inheritance and evolution that complements variation and explains the persistence of biological identity across generations and evolutionary time.

Biological identity

Properties of immunoaffinity purified 106-kDa Ca2+ release channels from the skeletal sarcoplasmic reticulum.

The sulfhydryl-gated 106-kDa Ca(2+)-release channel (SG-106) was purified by biotin-avidin chromatography from skeletal sarcoplasmic reticulum (SR) vesicles and used as an antigen to raise polyclonal antibodies. Western blots showed that the antisera crossreacted with the antigenic SG-106 and not with SR Ca2+, Mg(2+)-ATPase or with junctional foot proteins (JFPs) (Zaidi et al., 1989, J. Biol. Chem. 264(36), 21, 725-21, 736; 21, 737-21, 747). Polyclonal antibody-affinity columns were used to selectively purify SG-106-kDa proteins which, upon incorporation in planar bilayers, revealed the presence of a cationic channels with properties similar to "native" Ca(2+)-release channels obtained through the fusion of SR vesicles with planar bilayers. In agreement with measurements of Ca2+ release from SR vesicles, sulfhydryl oxidizing and reducing agents (i.e., 2,2'-dithiodipyridine and dithiothreitol) respectively increased and decreased the open-time probability of 106-kDa Ca(2+)-release channels. In contrast with reports on JFPs, ryanodine at 0.5-1 nM increased the open-time probability and at 2-10 nM locked 106-kDa Ca(2+)-release channels in a closed state rather than an open subconductance state. The SG-106 was activated by millimolar ATP, inhibited by millimolar Mg2+, and blocked by micromolar ruthenium red. Adriamycin (2-10 microM) caused a transient activation of SG-106 Ca(2+)-release channels, followed by closure in about 5 min, and intermittent activation to a subconductance state. Polyclonal antibodies used to purify the SG-106 also activated the channel when added to the cis side but not the trans side of the bilayer. Thus, SG-106 channels possess features that are similar to "native" SR Ca(2+)-release channels, are immunologically distinct from JFPs, and interact in seconds with nanomolar ryanodine in planar bilayers.

Animals

Characterization of West Nile virus Koutango lineage from phlebotomine sandflies in Kenya.

The West Nile virus (WNV), primarily transmitted by mosquitoes, is one of the most widespread flaviviruses globally, with past outbreaks occurring in the USA and Europe. Recent studies in parts of Africa, including Kenya, have identified the West Nile virus Koutango lineage (WN-KOUTV) among phlebotomine sandfly populations, however, our understanding of this virus remains limited. This study aimed to characterize WN-KOUTV from phlebotomine sandflies. Sandflies were sampled between 12th -16th March 2021 and 16th -20th March 2023 from six villages each in Baringo and Isiolo Counties, using CDC light traps. Female sandflies were taxonomically identified and pooled based on genus and site of collection. Virus isolation was performed in Vero cells. Viral genomes were determined using next-generation sequencing. Phylogenetic and molecular clock analyses were done to decipher the virus's evolutionary relationships. Comparative analyses of amino acid sequences were performed to determine variations. Protein modeling in Pymol was conducted to elucidate variations in key protein regions. Evolutionary pressure analysis investigated the selection pressures on the virus. In vitro experiments were done to investigate the virus growth kinetics in mammalian Vero E6 and mosquito C6/36 cells. We report the isolation of WN-KOUTV from Salabani in Baringo and Aremet in Isiolo, Kenya. The isolated WN-KOUTVs clustered with previously identified WN-KOUTV strains. Comparative analysis revealed a unique amino acid at NS5 653. The WN-KOUTV lineage as a whole is under purifying selective pressure, with diversifying pressure acting at site NS3 267. The current WN-KOUTV replicated in Vero E6 and C6/36 cells comparable to West Nile virus Lineage 1a, isolated from mosquitoes. Subsequent isolations of WN-KOUTV in phlebotomine sandflies suggest potential vectors, however, vector competence studies would confirm this. Replication in mammalian and insect cell lines suggests there may exist a vector/host relationship. We speculate the close genetic relationship of WN-KOUTV strains from East and West Africa may potentially be enabled by bird migratory routes between the two regions. If proven, this could point to a potential future pandemic pathway for this virus.

Animals

Evolutionary consequences of domestication on the selective effects of new amino acid changing mutations in canids.

The domestication of wild canids led to dogs no longer living in the wild but instead residing alongside humans. Extreme changes in behavior and diet associated with domestication may have led to the relaxation of the selective pressure on traits that may be less important in the domesticated context. Thus, here we hypothesize that strongly deleterious mutations may have become less deleterious in domesticated populations. We test this hypothesis by estimating the distribution of fitness effects (DFE) for new amino acid changing mutations using whole-genome sequence data from 24 gray wolves and 61 breed dogs. We find that the DFE is strikingly similar across canids, with 26-28% of new amino acid changing mutations being neutral/nearly neutral (|s| < 1e-5), and 41-48% under strong purifying selection (|s| > 1e-2). Our results are robust to different model assumptions suggesting that the DFE is stable across short evolutionary timescales, even in the face of putative drastic changes in the selective pressure caused by artificial selection during domestication and breed formation. On par with previous works describing DFE evolution, our data indicate that the DFE of amino acid changing mutations depends more strongly on genome structure and organismal characteristics, and less so on shifting selective pressures or environmental factors. Given the constant DFE and previous data showing that genetic variants that differentiate wolf and dog populations are enriched in regulatory elements, we speculate that domestication may have had a larger impact on regulatory variation than on amino acid changing mutations.

Journal Article

Selection profiles in RNA viruses reflect the characteristics of viruses more than individual proteins.

Proteins that are exposed on the surface of a virus are frequently subject to strong selection to escape from neutralizing antibodies. To investigate whether surface-exposed (SE) and non-exposed (NE) proteins encoded by RNA viruses exhibit different patterns of evolution under selection, we analyzed 244 protein-coding genes from 28 species of RNA viruses representing 15 taxonomic families. First, we show that gene-wide rates of non-synonymous (dN) and synonymous (dS) substitutions do not differentiate between SE and NE proteins. To incorporate variation in substitution rates among codon sites, we inferred the posterior distribution over a fixed grid of dN and dS rates for each alignment. This 'evolutionary fingerprint' provides a common framework for comparing the selection profiles of non-homologous genes. Next, we computed the Wasserstein distance for every pair of fingerprints, which is analogous to amount of work required to reshape one distribution to another. After compensating for differences in genetic variation among alignments, we found a small but significant difference between the fingerprints of SE and NE proteins (PERMANOVA, P&#x2009;=&#x2009;0.03). However, we observed larger and more significant effects of whether the virus is enveloped (P&#x2009;<&#x2009;10-5) and the interaction between these factors (P=6.9&#xd7;10-4). The latter effects were driven by high levels of purifying selection in capsid proteins of Picornaviruses. Furthermore, greater amounts of variation in fingerprints were explained by significant differences among virus families and modes of transmission (P&#x2009;<&#x2009;10-5). These results imply the pattern of selection on a virus protein is shaped more by characteristics of the virus than the protein itself.

RNA Viruses

Widespread horizontal transfer and strong selection enhance microbial adaptation in Antarctic soils.

Terrestrial Antarctica harbors compositionally diverse and functionally distinct microbial life. Yet the eco-evolutionary processes underlying adaptation to Antarctica's polyextreme conditions remain largely unknown. Here, we address how horizontal gene transfer (HGT) and de novo mutations influence microbial adaptation in 16 Antarctic soils using combined short- and long-read datasets. Phylogenetic reconciliation and mobile genetic element analysis of 676 metagenome-assembled genomes show frequent HGT across communities. While transferred genes span diverse functional categories, those involved in energy metabolism are exchanged at higher frequency. Genes for aerotrophy, i.e. the consumption of atmospheric trace gases to provide energy, carbon, and hydration, are among the most frequently disseminated. Approximately a quarter of carbon monoxide dehydrogenases and [NiFe]-hydrogenases are predicted to be horizontally acquired and are often associated with mobile genetic elements. Analysis of polymorphisms suggests widespread purifying selection, particularly for aerotrophy genes, providing further evidence that aerotrophy is critical for microbial survival in Antarctica. Genetic variation in hydrogenases is tightly associated with predicted protein structures, with intense selection acting on critical sites preserving stability and function. Together, these findings show that previously unrecognized eco-evolutionary dynamics shape the composition and function of Antarctic microbial communities, and confirm aerotrophy is a strongly selected and horizontally disseminated trait.

Antarctic Regions

Calcium-binding modulator protein from the unfertilized egg of the sea urchin Arbacia punctulata.

We have purified and partly characterized a calcium-binding protein from the unfertilized egg of the sea urchin Arbacia punctulata. This protein closely resembles the calcium-binding modulator protein of bovine brain in its molecular weight, electrophoretic mobility, amino acid analysis, and peptide map. It activates bovine brain phosphodiesterase in the presence of calcium but has no effect on the phosphodiesterase of the Arbacia egg. Densitometric scanning of acrylamide gels of arbacia egg homogenates shows the modulator protein to represent 0.1% of the total protein of the egg. At 10(-4) M free calcium, the protein binds four calcium ions per 17,000-dalton molecule. We have used a column of rabbit skeletal muscle troponin-I covalently coupled to Sepharose 4B as an affinity column to selectively purify the Arbacia egg calcium-binding protein. This column has also been used to purify bovine brain modulator protein and may prove of general use in isolating similar proteins from other sources. The technique may be particularly helpful when only small quantities of starting material are available.

3',5'-Cyclic-AMP Phosphodiesterases

Evidence of genome-wide relaxed selection on mildly deleterious mutations in an ancient subterranean catfish.

About one hundred subterranean catfish species have been described, resulting from repeated colonization of cave environments by multiple surface lineages. Most cave-dwelling species are found in the Americas, in particular in South America, but a few species also live in Central and North America. Despite the availability of high-quality genome assemblies for two cave species, the Mexican blind catfish Prietella phreatophila and the Colombian blind catfish Trichomycterus rosablanca, genomic approaches to investigate genetic changes associated with subterranean life or to estimate cave colonization times remain largely unexplored. To fill this gap, we additionally sequenced the genomes of four blind and depigmented subterranean catfishes from Peru (three Trichomycterus and one Astroblepus), as well as the genomes of four close surface relatives. We first extracted a large set of light-related genes, such as phototransduction and crystallin genes, and found contrasting decays of these sequences in different cave species, from 1% of pseudogenes in T. rosablanca to 48% in P. phreatophila. Two independent molecular dating methods gave congruent ages, indicating that these catfishes colonized subterranean habitats at different times, ranging from Early Pliocene to Late Pleistocene, supporting the hypothesis that surface catfishes repeatedly and rapidly adapted to subterranean habitats. The oldest cave species, P. phreatophila, appears to have been thriving in the dark for over 3.5 million years. Moreover, a genome-wide analysis of protein-coding genes suggests weaker purifying selection on mildly deleterious mutations in this cavefish than in other catfish lineages, likely reflecting a long-term small effective population size.

cavefishes

Diversification of Cellulose Synthase (CESA) Genes in Mosses Suggests Both Ancient and Recent Gene duplications.

Cellulose is an important polysaccharide that constitutes all plant cell walls, giving them strength and stability. The plant cellulose synthase (CESA) gene family, which encodes the catalytic subunits of cellulose synthesis complexes (CSCs), has diversified independently in several plant lineages, providing an interesting model for understanding selection for gene duplication. Here we quantified the presence of CESA genes across mosses to understand how the process of gene family diversification occurred in this group and how it parallels diversification in other groups. We first examined the CESA gene family in eight species of mosses across seven families for which whole genome assemblies were available. We then identified CESA genes from additional species, for which only short-read sequence data was available, by using BLAST searches and targeted gene assemblies. We validated this approach by comparing the assembled paralogs from the short-read data to the genes identified from whole genome assemblies in the eight reference species. This approach allowed us to identify paralogs directly from short-read data and greatly expand our sample set. Results from the combined empirical data support the hypothesis that CESA genes diversified within the moss lineage at least as early as the mesozoic period, during or possibly even prior to the onset of moss diversification, but also continue to diversify within modern species. In addition, we found evidence for purifying selection as the dominant force shaping these genes and observed that different lineages experienced different levels of evolutionary constraint. Lastly, our approach to assemble paralogs has the potential to allow researchers to improve analyses of gene duplication events.

Physcomitrium patens

Transposable Element Dynamics Drive the Genomic Evolution and Phenotypic Diversification of Allotetraploid Common Carp.

An important question in evolutionary biology is how polyploidization generates raw material for phenotypic diversification. Transposable elements (TEs) represent an underestimated source of genetic variation in eukaryotic genomes. By integrating 516 whole-genome resequencing datasets and 236 transcriptomes from common carp (Cyprinus carpio), a representative allotetraploid fish, we constructed the first population-scale landscape of TE insertions in teleosts. TE insertions are widespread in the carp genome and preferentially associated with stress-responsive genes, with DNA transposons as major contributors. Relaxed purifying selection and TE burst events coexist, generating abundant variation for subsequent subspecies differentiation. Compared with a closely related diploid species, carp exhibits more exonic TE insertions and shorter TE-gene distances, and multiple TE superfamilies expanded during tetraploidization. Genome-wide association analyses uncovered intragenic TE variants underlying domesticated traits missed by SNPs, including DNA transposon deletions associated with scale reduction and altered body shape. Notably, lighter-colored individuals harbor homozygous deletions of LTR and DNA transposons within mdfic2, whose knockout in zebrafish reduces pigmentation. Most trait-associated variants reflect lineage-specific loss of ancient TE insertions rather than recent transposition. Overall, these findings highlight the distinct role of TEs in polyploid genome evolution and phenotypic diversification, providing new insights into TE dynamics in vertebrates.

allotetraploidization

Impact of High-Titer Convalescent Plasma on Clinical and Virologic Outcomes Among Veterans Hospitalized With SARS-CoV-2 Infection: VA CoronavirUs Research and Efficacy Studies-1 (VA CURES-1).

In the initial absence of proven therapies, empirical COVID-19 convalescent plasma (CCP) was rapidly introduced for individuals hospitalized for COVID-19. Seventy-five participants were randomized from November 2020 to June 2021 in a double-blind, multi-site, placebo-controlled, randomized trial (VA CURES-1) evaluating the impact of CCP vs. saline in Veterans hospitalized with COVID-19 with hypoxemia. The composite primary outcome was acute hypoxemic respiratory failure or all-cause death by Day 29. We analyzed clinical outcomes, nasal viral RNA, plasma cytokines and viral evolution over time. Among 40 participants receiving saline and 35 receiving CCP with high neutralizing titers (median 1:1420), the percent reaching the primary outcome was similar (10%), as were time to clinical recovery and to nasal viral clearance. By whole genome sequencing, viral molecular complexity evolved pre- to posttreatment more frequently in recipients of saline vs. CCP (4 of 7 (57.1%) vs. 1 of 4 (25%), respectively), based on numbers of mixed allele positions. Numbers of amino acid-changing, non-synonymous mutations in the spike protein were greater in saline vs. CCP recipients. Both outcomes suggested purifying selection (reduced overall viral infection complexity) following CCP. In conclusion, convalescent plasma showed no significant clinical impact but may influence SARS-CoV-2 complexity. Trial Registration: ClinicalTrials.gov Identifier: NCT04539275.

Aged

Dynamic and non-additive gene regulation shapes maize responses to simultaneous salt and cold stress.

Salt and cold stresses often occur together in nature and severely impact crop productivity, yet their transcriptional regulation remains poorly understood. Here, we conducted a time-series transcriptomic analysis of maize under salt, cold, and their combination at 0, 6, 12, and 24&#xa0;h. Differential expression analysis revealed dynamic, condition-specific gene responses grouped into eight distinct temporal patterns. Promoter motif analysis of genes within each pattern identified 5-39 significantly enriched motifs, with over 40% lacking known counterparts, suggesting the involvement of previously uncharacterized cis-regulatory elements in stress-responsive transcriptional regulation. By comparing combined stress responses to the sum of single-stress effects, we found that about 74% of DEGs showed non-additive patterns, suggesting that combined stress triggers a distinct transcriptional program. Evolutionary analysis showed that additive DEGs tend to be more recently evolved, subject to weaker purifying selection, and enriched in transposed duplications, contrasting with the stronger constraint observed in non-additive DEGs. WGCNA identified 24 co-expression modules, among which 65 hub DEGs were detected in modules significantly correlated with specific stress conditions. Furthermore, we reconstructed 228, 20, and 200 sequential transcription factor cascades spanning 6&#xa0;h, 12&#xa0;h, and 24&#xa0;h under cold, salt, and combined stress, respectively, with no cascade shared across all three conditions. Together, these results reveal that maize responses to combined salt and cold stress are largely non-additive and temporally dynamic, with distinct evolutionary patterns underlying different response types, offering insights and candidate regulators for enhancing crop stress resilience.

Zea mays

Assembly and Characterization of the First Complete Mitochondrial Genome of Tussilago farfara L.: Insights into Biological Functions and Phylogenetic Relationships within the Asteraceae Family.

Tussilago farfara L., a member of the Asteraceae family, is an economically valuable species due to its edible and medicinal properties. To elucidate the structural characteristics, genetic mechanisms, and evolutionary pathways of the organelle genomes of T. farfara, we sequenced, assembled, and annotated its mitochondrial genome for the first time. The complete mitochondrial genome of T. farfara spans 306,024&#xa0;bp and contains 33 mitochondrial protein-coding genes (PCGs), 3 rRNAs, and 22 tRNAs. Analysis of the nucleotide substitution rate and genetic diversity revealed that most mitochondrial genome genes may have undergone purifying selection, indicating a slow evolutionary rate and a relatively conserved genomic structure. We further identified 13 fragments of chloroplast-derived DNA integrated into the mitochondrial genome, evidencing intracellular gene transfer. Collinearity analysis showed that Arctium lappa shares the most extensive mitochondrial homologous sequences and the highest sequence similarity with T. farfara. Phylogenetic analysis based on the mitochondrial genome helped to clarify the evolutionary and taxonomic position of T. farfara within the Asteraceae family. The mitochondrial genome sequence of T. farfara provides a valuable genomic resource for species identification and for evolutionary studies within the Asteraceae family.

Genome, Mitochondrial