SEARCH · Search PubMed

Results for “evolutionary dynamics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 667 records · Page 37Linked to original sources

PubMed2026

Potential survival strategies of novel comammox and nitrite-oxidizing Nitrospira synthesizing osmoprotectants in a wastewater microbiome treating high-ammonia brackish landfill leachate.

BACKGROUND: In the late stages of landfill operation, leachate becomes brackish and contains high concentrations of ammonia with limited organic carbon. At leachate treatment facilities, it is typically subjected to nitrification followed by denitrification, with methanol supplied as an external electron donor. This unique environment may harbor novel microorganisms, including nitrifiers. Although a variety of microorganisms are involved in nitrification, their substrate specificity and salinity tolerance remain insufficiently understood. In this study, a genome-centric metagenome analysis was conducted on the microbiome from a leachate treatment facility at a closed landfill. RESULTS: A total of 68 metagenome-assembled genomes (MAGs) were reconstructed, including 64 putative novel species. Among these, two Nitrospira MAGs were recovered: a novel complete ammonia-oxidizing bacterium (comammox), Nitrospira LAS72 (88.72% completeness, 2.10% contamination), and canonical nitrite-oxidizing Nitrospira LAS18 (99.98% completeness, 2.29% contamination). Comparative genomic analysis with 260 publicly available Nitrospira genomes revealed that LAS18 represents a new sub-lineage within lineage VII of the Nitrospira genus. Two ammonia-oxidizing archaea (AOA), Candidatus Nitrosocosmicus LAS21 and Nitrosarchaeum LAS73, were also identified, while canonical ammonia-oxidizing bacteria were not detected. Given the brackish conditions (1.23% salinity) and the methanol-fed operation of the treatment facility, the genomic potential for osmotic stress adaptation and methanol metabolism was investigated. Comammox Nitrospira LAS72 harbors biosynthetic pathways for several compatible solutes (osmoprotectants), including glycine betaine, proline, trehalose, and L-glutamate. Moreover, comammox Nitrospira LAS72 possesses genetic potential for oxidizing formaldehyde, suggesting that it may exploit these methanol-derived intermediates as energy sources. These features indicate that LAS72 may withstand osmotic fluctuations through the production of various osmoprotectants and thrive under the unique conditions of a methanol-fed environment. CONCLUSIONS: The discovery of novel comammox Nitrospira and canonical Nitrospira forming a new sub-lineage within lineage VII of the Nitrospira genus in an ammonia-rich brackish environment provides the first genomic evidence for evolutionary adaptation among nitrifiers to saline, methanol-fed environments. These findings enhance our understanding of the ecological and evolutionary dynamics shaping nitrifier communities in complex treatment ecosystems. Video Abstract.

Ammonia↗

PubMed2007

[Directed evolution of antibody molecules in phage-displayed combinatorial libraries].

Advances in methods for conformational prediction, structural analysis and site-directed mutagenesis of proteins and peptides have contributed to the understanding of their structure and function. However, with the exception of a few successes, the generation of practical functional molecules solely by rational design remains a difficult challenge. The aim of our study is to investigate molecular design relying on evolutionary processes, called as "directed evolution", to generate a novel class of biofunctional molecules. This evolutionary approach consists of three steps; 1) constructions of protein/peptide libraries based on structural information, 2) expressions of the libraries on phage particles, and 3) selections with investigator-imposed selective pressures. In this work, we study on directed evolution with antibody libraries. We have succeeded in generating highly active catalytic antibodies in phage-displayed antibody (Fab) libraries. To evolve catalytic antibodies toward higher catalytic activity, we have mimicked an enzyme-evolutional process, in which an enzyme has evolved their ability to use binding energies for catalysis by increasing the affinity for the transition state of a reaction and decreasing the affinity for the ground state. Thus, phage-displayed libraries derived from an original catalytic antibody were selected against a newly-devised TSA, which was programmed to optimize the differential affinity for the transition state relative to the ground state, to provide variants with improved reaction rates (k(cat)). The in vitro evolution has great potential for generating novel catalysts as well as for providing opportunities to examine the evolutionary dynamics of enzymes.

Antibodies↗

PubMed2002

Molecular phylogeny of Japanese Leporidae, the Amami rabbit Pentalagus furnessi, the Japanese hare Lepus brachyurus, and the mountain hare Lepus timidus, inferred from mitochondrial DNA sequences.

We determined mitochondrial 12S ribosomal RNA (rRNA) and cytochrome b (cyt b) gene sequences in three leporid species of Japan, the Amami rabbit Pentalagus furnessi from the Ryukyu Islands, the Japanese hare Lepus brachyurus from Honshu, and a Japanese form of the mountain hare Lepus timidus ainu from Hokkaido. We compared the sequences with those of other taxa of leporids available in databases. Phylogenetic trees of the 12S rRNA gene sequences indicated that the lineage of P. furnessi diversified during the generic radiation of the leporids at an ancient time, which was estimated to have been the middle Miocene. Cyt-b gene trees revealed that the lineage of L. brachyurus branched off at an early stage in the speciation of Lepus, probably at the beginning of the Pliocene. The cyt b sequences of L. t. ainu were somewhat distinct from those of continental conspecific populations; this lineage divergence is likely to have occurred during the middle or late Pleistocene. The results show that the three regions of the Japanese archipelago, Ryukyu, Honshu-Shikoku-Kyushu, and Hokkaido, now preserve their own leporid taxa, each with a different extent of genetic endemicity. It is possible that the zoogeographic traits of the Japanese leporids are a consequence of the evolutionary dynamics of leporids in East Asia, in that the radiation centers of leporids are likely to have shifted from tropical, through temperate, to arctic zones.

Animals↗

PubMed2006

Prevalence and evolution of core photosystem II genes in marine cyanobacterial viruses and their hosts.

Cyanophages (cyanobacterial viruses) are important agents of horizontal gene transfer among marine cyanobacteria, the numerically dominant photosynthetic organisms in the oceans. Some cyanophage genomes carry and express host-like photosynthesis genes, presumably to augment the host photosynthetic machinery during infection. To study the prevalence and evolutionary dynamics of this phenomenon, 33 cultured cyanophages of known family and host range and viral DNA from field samples were screened for the presence of two core photosystem reaction center genes, psbA and psbD. Combining this expanded dataset with published data for nine other cyanophages, we found that 88% of the phage genomes contain psbA, and 50% contain both psbA and psbD. The psbA gene was found in all myoviruses and Prochlorococcus podoviruses, but could not be amplified from Prochlorococcus siphoviruses or Synechococcus podoviruses. Nearly all of the phages that encoded both psbA and psbD had broad host ranges. We speculate that the presence or absence of psbA in a phage genome may be determined by the length of the latent period of infection. Whether it also carries psbD may reflect constraints on coupling of viral- and host-encoded PsbA-PsbD in the photosynthetic reaction center across divergent hosts. Phylogenetic clustering patterns of these genes from cultured phages suggest that whole genes have been transferred from host to phage in a discrete number of events over the course of evolution (four for psbA, and two for psbD), followed by horizontal and vertical transfer between cyanophages. Clustering patterns of psbA and psbD from Synechococcus cells were inconsistent with other molecular phylogenetic markers, suggesting genetic exchanges involving Synechococcus lineages. Signatures of intragenic recombination, detected within the cyanophage gene pool as well as between hosts and phages in both directions, support this hypothesis. The analysis of cyanophage psbA and psbD genes from field populations revealed significant sequence diversity, much of which is represented in our cultured isolates. Collectively, these findings show that photosynthesis genes are common in cyanophages and that significant genetic exchanges occur from host to phage, phage to host, and within the phage gene pool. This generates genetic diversity among the phage, which serves as a reservoir for their hosts, and in turn influences photosystem evolution.

Base Composition↗

PubMed2025

A metabolic atlas of the Klebsiella pneumoniae species complex reveals lineage-specific metabolism and capacity for intra-species co-operation.

The Klebsiella pneumoniae species complex inhabits a wide variety of hosts and environments, and is a major cause of antimicrobial resistant infections. Genomics has revealed the population comprises multiple species/sub-species and hundreds of distinct co-circulating sub-lineage (SLs) that are associated with distinct gene complements. A substantial fraction of the pan-genome is predicted to be involved in metabolic functions and hence these data are consistent with metabolic differentiation at the SL level. However, this has so far remained unsubstantiated because in the past it was not possible to explore metabolic variation at scale. Here, we used a combination of comparative genomics and high-throughput genome-scale metabolic modeling to systematically explore metabolic diversity across the K. pneumoniae species complex (n = 7,835 genomes). We simulated growth outcomes for each isolate using carbon, nitrogen, phosphorus, and sulfur sources under aerobic and anaerobic conditions (n = 1,278 conditions per isolate). We showed that the distributions of metabolic genes and growth capabilities are structured in the population, and confirmed that SLs exhibit unique metabolic profiles. In vitro co-culture experiments demonstrated reciprocal commensalistic cross-feeding between SLs, effectively extending the range of conditions supporting individual growth. We propose that these substrate specializations may promote the existence and persistence of co-circulating SLs by reducing nutrient competition and facilitating commensal interactions. Our findings have implications for understanding the eco-evolutionary dynamics of K. pneumoniae and for the design of novel strategies to prevent opportunistic infections caused by this World Health Organization priority antimicrobial resistant pathogen.

Klebsiella pneumoniae↗

PubMed2006

Modeling somatic evolution in tumorigenesis.

Tumorigenesis in humans is thought to be a multistep process where certain mutations confer a selective advantage, allowing lineages derived from the mutated cell to outcompete other cells. Although molecular cell biology has substantially advanced cancer research, our understanding of the evolutionary dynamics that govern tumorigenesis is limited. This paper analyzes the computational implications of cancer progression presented by Hanahan and Weinberg in The Hallmarks of Cancer. We model the complexities of tumor progression as a small set of underlying rules that govern the transformation of normal cells to tumor cells. The rules are implemented in a stochastic multistep model. The model predicts that (i) early-onset cancers proceed through a different sequence of mutation acquisition than late-onset cancers; (ii) tumor heterogeneity varies with acquisition of genetic instability, mutation pathway, and selective pressures during tumorigenesis; (iii) there exists an optimal initial telomere length which lowers cancer incidence and raises time of cancer onset; and (iv) the ability to initiate angiogenesis is an important stage-setting mutation, which is often exploited by other cells. The model offers insight into how the sequence of acquired mutations affects the timing and cellular makeup of the resulting tumor and how the cellular-level population dynamics drive neoplastic evolution.

Cell Transformation, Neoplastic↗

PubMed2006

Large-scale turnover of functional transcription factor binding sites in Drosophila.

The gain and loss of functional transcription factor binding sites has been proposed as a major source of evolutionary change in cis-regulatory DNA and gene expression. We have developed an evolutionary model to study binding-site turnover that uses multiple sequence alignments to assess the evolutionary constraint on individual binding sites, and to map gain and loss events along a phylogenetic tree. We apply this model to study the evolutionary dynamics of binding sites of the Drosophila melanogaster transcription factor Zeste, using genome-wide in vivo (ChIP-chip) binding data to identify functional Zeste binding sites, and the genome sequences of D. melanogaster, D. simulans, D. erecta, and D. yakuba to study their evolution. We estimate that more than 5% of functional Zeste binding sites in D. melanogaster were gained along the D. melanogaster lineage or lost along one of the other lineages. We find that Zeste-bound regions have a reduced rate of binding-site loss and an increased rate of binding-site gain relative to flanking sequences. Finally, we show that binding-site gains and losses are asymmetrically distributed with respect to D. melanogaster, consistent with lineage-specific acquisition and loss of Zeste-responsive regulatory elements.

Animals↗

PubMed2025

Phylogenetic diversity and molecular evolution of Hantaan virus harbored by Apodemus chejuensis on Jeju Island, Republic of Korea, 2022-2023.

BACKGROUND: Hantaan virus (HTNV), hosted by Apodemus spp., is a well-recognized causative agent of hemorrhagic fever with renal syndrome (HFRS) and poses a crucial global public health concern. Based on the current evidence, HTNV carried by A. chejuensis is proposed as the likely etiological agent of HFRS on Jeju Island, Republic of Korea (ROK). METHODOLOGY/PRINCIPAL FINDINGS: In this study, 50 small mammals were collected from five locations in Seogwipo-si and Jeju-si on Jeju Island, ROK, during 2022-2023. Serological and molecular analyses revealed HTNV prevalence rates of 34% (16/47) and 27.7% (13/47), respectively. Using a multiplex polymerase chain reaction-based nanopore sequencing approach, nine complete HTNV genomes were sequenced from the lung tissues of A. chejuensis, representing the first comprehensive genomic characterization of HTNV from Seogwipo-si (Hogeun-dong) and Jeju-si (Sangdae-ri). Phylodynamic analyses suggest evolutionary divergence and phylogeographic diversity, with four unique amino acid substitutions identified in HTNV genomes from Seogwipo-si. CONCLUSION/SIGNIFICANCE: This study provides important insights into the genomic surveillance, genetic diversity, and evolutionary dynamics of orthohantaviruses, which are essential for guiding effective public health strategies to control and prevent future HFRS outbreaks in the ROK.

Animals↗

PubMed2006

Experimental rugged fitness landscape in protein sequence space.

The fitness landscape in sequence space determines the process of biomolecular evolution. To plot the fitness landscape of protein function, we carried out in vitro molecular evolution beginning with a defective fd phage carrying a random polypeptide of 139 amino acids in place of the g3p minor coat protein D2 domain, which is essential for phage infection. After 20 cycles of random substitution at sites 12-130 of the initial random polypeptide and selection for infectivity, the selected phage showed a 1.7x10(4)-fold increase in infectivity, defined as the number of infected cells per ml of phage suspension. Fitness was defined as the logarithm of infectivity, and we analyzed (1) the dependence of stationary fitness on library size, which increased gradually, and (2) the time course of changes in fitness in transitional phases, based on an original theory regarding the evolutionary dynamics in Kauffman's n-k fitness landscape model. In the landscape model, single mutations at single sites among n sites affect the contribution of k other sites to fitness. Based on the results of these analyses, k was estimated to be 18-24. According to the estimated parameters, the landscape was plotted as a smooth surface up to a relative fitness of 0.4 of the global peak, whereas the landscape had a highly rugged surface with many local peaks above this relative fitness value. Based on the landscapes of these two different surfaces, it appears possible for adaptive walks with only random substitutions to climb with relative ease up to the middle region of the fitness landscape from any primordial or random sequence, whereas an enormous range of sequence diversity is required to climb further up the rugged surface above the middle region.

Amino Acid Sequence↗

PubMed2024

Assessing the emergence time of SARS-CoV-2 zoonotic spillover.

Understanding the evolution of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) and its relationship to other coronaviruses in the wild is crucial for preventing future virus outbreaks. While the origin of the SARS-CoV-2 pandemic remains uncertain, mounting evidence suggests the direct involvement of the bat and pangolin coronaviruses in the evolution of the SARS-CoV-2 genome. To unravel the early days of a probable zoonotic spillover event, we analyzed genomic data from various coronavirus strains from both human and wild hosts. Bayesian phylogenetic analysis was performed using multiple datasets, using strict and relaxed clock evolutionary models to estimate the occurrence times of key speciation, gene transfer, and recombination events affecting the evolution of SARS-CoV-2 and its closest relatives. We found strong evidence supporting the presence of temporal structure in datasets containing SARS-CoV-2 variants, enabling us to estimate the time of SARS-CoV-2 zoonotic spillover between August and early October 2019. In contrast, datasets without SARS-CoV-2 variants provided mixed results in terms of temporal structure. However, they allowed us to establish that the presence of a statistically robust clade in the phylogenies of gene S and its receptor-binding (RBD) domain, including two bat (BANAL) and two Guangdong pangolin coronaviruses (CoVs), is due to the horizontal gene transfer of this gene from the bat CoV to the pangolin CoV that occurred in the middle of 2018. Importantly, this clade is closely located to SARS-CoV-2 in both phylogenies. This phylogenetic proximity had been explained by an RBD gene transfer from the Guangdong pangolin CoV to a very recent ancestor of SARS-CoV-2 in some earlier works in the field before the BANAL coronaviruses were discovered. Overall, our study provides valuable insights into the timeline and evolutionary dynamics of the SARS-CoV-2 pandemic.

Animals↗

PubMed2005

Replication of an Egfr-wing shape association in a wild-caught cohort of Drosophila melanogaster.

Linkage disequilibrium mapping has been used extensively in medical and evolutionary genetics to map causal polymorphisms within genes associated with disease status or phenotypic variation for a trait. However, the initial findings of most nonhuman studies have not been replicated in subsequent studies, due in part to false positives, as well as additional factors that can render true positives unreplicable. These factors may be more severe when the initial study is performed using an experimental population of organisms reared under controlled lab conditions. We demonstrate that despite considerable phenotypic differences for wing shape between a lab-reared experimental population and a wild-caught cohort of Drosophila melanogaster, an association between a putative regulatory polymorphism in Egfr and wing shape can be replicated. These results are discussed both within the framework of future association-mapping studies and within the context of the evolutionary dynamics of alleles in populations.

Alleles↗

PubMed2005

Cross-species comparison of Drosophila male accessory gland protein genes.

Drosophila melanogaster males transfer seminal fluid proteins along with sperm during mating. Among these proteins, ACPs (Accessory gland proteins) from the male's accessory gland induce behavioral, physiological, and life span reduction in mated females and mediate sperm storage and utilization. A previous evolutionary EST screen in D. simulans identified partial cDNAs for 57 new candidate ACPs. Here we report the annotation and confirmation of the corresponding Acp genes in D. melanogaster. Of 57 new candidate Acp genes previously reported in D. melanogaster, 34 conform to our more stringent criteria for encoding putative male accessory gland extracellular proteins, thus bringing the total number of ACPs identified to 52 (34 plus 18 previously identified). This comprehensive set of Acp genes allows us to dissect the patterns of evolutionary change in a suite of proteins from a single male-specific reproductive tissue. We used sequence-based analysis to examine codon bias, gene duplications, and levels of divergence (via dN/dS values and ortholog detection) of the 52 D. melanogaster ACPs in D. simulans, D. yakuba, and D. pseudoobscura. We show that 58% of the 52 D. melanogaster Acp genes are detectable in D. pseudoobscura. Sequence comparisons of ACPs shared and not shared between D. melanogaster and D. pseudoobscura show that there are separate classes undergoing distinctly dissimilar evolutionary dynamics.

Animals↗

PubMed2005

Monitoring the mode and tempo of concerted evolution in the Drosophila melanogaster rDNA locus.

Non-LTR retrotransposons R1 and R2 have persisted in rRNA gene loci (rDNA) since the origin of arthropods despite their continued elimination by the recombinational mechanisms of concerted evolution. This study evaluated the short-term evolutionary dynamics of the rDNA locus by measuring the divergence among replicate Drosophila melanogaster lines after 400 generations. The total number of rDNA units on the X chromosome of each line varied from 140 to 310, while the fraction of units inserted with R1 and R2 retrotransposons ranged from 37 to 65%. This level of variation is comparable to that found in natural population surveys. Variation in locus size and retrotransposon load was correlated with large changes in the number of uninserted and R1-inserted units, yet the numbers of R2-inserted units were relatively unchanged. Intergenic spacer (IGS) region length variants were also used to evaluate changes in the rDNA loci. All IGS length variants present in the lines showed significant increases and decreases of copy number. These studies, combined with previous data following specific R1 and R2 insertions in these lines, help to define the type and distribution, both within the locus and within the individual units, of recombinational events that give rise to the concerted evolution of the rDNA locus.

Animals↗

PubMed2004

The effective size of mixed sexually and asexually reproducing populations.

Using the transition matrix of inbreeding and coancestry coefficients, the inbreeding (N(eI)), variance (N(eV)), and asymptotic (N(e lambda)) effective sizes of mixed sexual and asexual populations are formulated in terms of asexuality rate (delta), variance of asexual (C) and sexual (K) reproductive contributions of individuals, correlation between asexual and sexual contributions (rho(ck)), selfing rate (beta), and census population size (N). The trajectory of N(eI) toward N(e lambda) changes crucially depending on delta, N, and beta, whereas that of N(eV) is rather consistent. With increasing asexuality, N(e lambda) either increases or decreases depending on C, K, and rho(ck). The parameter space in which a partially asexual population has a larger N(e lambda) than a fully sexual population is delineated. This structure is destroyed when N(1 - delta) < 1 or delta > 1 - 1/N. With such a high asexuality, tremendously many generations are required for the asymptotic size N(e lambda) to be established, and N(e lambda) is extremely large with any value of C, K, and rho(ck) because the population is dominated eventually by individuals of the same genotype and the allelic diversity within the individuals decays quite slowly. In reality, the asymptotic state would occur only occasionally, and instantaneous rather than asymptotic effective sizes should be practical when predicting evolutionary dynamics of highly asexual populations.

Alleles↗

PubMed2005

Habitat avoidance: overlooking an important aspect of host-specific mating and sympatric speciation?

Understanding speciation requires discerning how reproductive barriers to gene flow evolve between previously interbreeding populations. Models of sympatric speciation for phytophagous insects posit that reproductive isolation can evolve in the absence of geographic isolation as a consequence of an insect shifting and ecologically adapting to a new host plant. One important adaptation contributing to sympatric differentiation is host-specific mating. When organisms mate in preferred habitats, a system of positive assortative mating is established that facilitates sympatric divergence. Models of host fidelity generally assume that host choice is determined by the aggregate effect of alleles imparting positive preferences for different plant species. But negative effect genes for avoiding nonnatal plants may also influence host use. Previous studies have shown that apple and hawthorn-infesting races of Rhagoletis pomonella flies use volatile compounds emitted from the surface of fruit as key chemosensory cues to recognize and distinguish between their host plants. Here, we report results from field trials indicating that in addition to preferring the odor of their natal fruit, apple and hawthorn flies, and their undescribed sister species infesting flowering dogwood (Cornus florida), also avoid the odors of nonnatal fruit. We discuss the implications of nonnatal fruit avoidance for the evolutionary dynamics and genetics of sympatric speciation. Our findings reveal an underappreciated role for habitat avoidance as a potential postmating, as well as prezygotic, barrier to gene flow.

Animals↗

PubMed2001

Evolution of the DUT gene: horizontal transfer between host and pathogen in all three domains of life.

The ubiquity of the dut gene in Eukarya, Eubacteria, and Archaea implies its existence in the last common ancestor of the three domains of life. The dut gene exists as single, tandemly duplicated, and tandemly triplicated copies. The dUTPase is encoded as an auxiliary gene in the genomes of several DNA viruses and two distinct lineages of retroviruses. A comprehensive analysis of dUTPase amino acid sequence relationships explores the evolutionary dynamics of dut genes in viruses and their hosts. The data set was comprised of representative sequences from available Eukaryotes, Archaea, Eubacteria cells and viruses. A multiple alignment of these protein sequences was generated using a hidden Markov model (HMM) approach developed to align divergent data. Phylogenetic analysis revealed that horizontal transfer from hosts to virus genomes has occurred in all three domains of life. The evidence for horizontal transfers is particularly interesting in Eukaryotes as these dut genes have introns, while DNA virus dut genes do not. This implies an intermediary Retroid Agent facilitated the horizontal transfer process, via reverse transcription, between host mRNA and DNA viruses. The horizontal transfer of the dut gene from Eukaryotic, Eubacterial, and Archaeal organisms to both DNA and RNA viruses is the first documented case of host to pathogen transfer that has occurred in all three domains of life.

Amino Acid Sequence↗

PubMed1993

Cytochrome oxidase subunit III from Arbacia lixula: detection of functional constraints by comparison with homologous sequences.

In this paper we report the comparison of the sequences of the cytochrome oxidase subunit III from three different sea urchin species. Both nucleotide and amino acid sequences have been analyzed. The nucleotide sequence analysis reveals that the sea urchin sequences obey some rules already found in mammals. The base substitution analysis carried out on the sequences of the three species pairs, shows that the evolutionary dynamics of the first and the second codon positions are so slow that do not allow a quantitative measurement of their genetic distances, thus demonstrating that also in these species the COIII gene is strongly conserved during evolution. Changes occurring at the third codon positions indicate that the three species evolved from a common ancestor under different directional mutational pressure. The multi-alignment of the sea urchin proteins indicates the existence of the amino acid sequence motif N R T that represents a possible glycosylation site. Another glycosylation site has been detected in the mammalian cytochrome oxidase subunit III, in a position slightly different. Such an analysis revealed, for the first time, a new functional aspect of this sequence.

Amino Acid Sequence↗

PubMed2026

Mitogenome assembly and phylogenetic relationships of Phalaris arundinacea.

INTRODUCTION: As a perennial herb of Poaceae, Phalaris arundinacea plays key roles in grazing, production, and soil and water conservation because of its well-developed rhizomes and seed dispersal. We assembled and annotated the first mitogenome of P. arundinacea to support evolutionary and taxonomic research. METHODS: We assembled and annotated the first complete mitochondrial genome of P. arundinacea by integrating Illumina short reads with Nanopore long reads via a hybrid assembly strategy. The genome architecture was comprehensively characterized, encompassing codon usage bias, repetitive sequence organization, and inter-organellar genetic exchange with the chloroplast genome. RESULTS AND DISCUSSION: Assembly of the P. arundinacea mitogenome revealed two circular structures with a combined length of 526,717 bp. The genome comprised a set of 37 protein-coding genes (PCGs), 27 tRNAs, and 8 rRNAs, with the rRNA genes exhibiting full assembly (100% coverage). The mitochondrial genome contained 154 forward and 164 palindromic repeats, along with 25 tandem repeats and 124 simple sequence repeats (SSRs). Notably, 102 SSRs were distributed on contig1, predominantly in tetrameric form. Furthermore, 376 RNA editing sites were predicted. A total of 104 fragments were integrated into the mitochondrial genome from the chloroplast, amounting to 55,866 bp of transferred sequence. Finally, phylogenetic analysis of 28 plant mitogenomes placed P. arundinacea closest to species within the genus Poa (P. chaixii and P. pratensis). Comparative analysis of non-synonymous-to-synonymous substitution rate (Ka/Ks) ratios across divergent species revealed that the mitochondrial genome of P. arundinacea underwent stabilizing evolutionary dynamics, characterized by predominant purifying selection with several lineage-specific variations in selective pressure. Our findings support the close phylogenetic relationship between P. arundinacea and species of the genus Poa and provide a reference mitochondrial genome resource for future comparative studies within Phalaris that incorporate broader taxon sampling. These results support deeper phylogenetic investigations of P. arundinacea and facilitate future work on its germplasm characterization and applied use.

Phalaris arundinacea↗