Sequence Repeat Research Articles

"Off-pore" nucleoporins relocalize heterochromatic breaks through phase separation.

Phase separation forms membraneless compartments in the nuclei, including by establishing heterochromatin "domains" and repair foci. Pericentromeric heterochromatin mostly comprises repeated sequences prone to aberrant recombination, and "safe" homologous recombination (HR) repair of these sequences requires the movement of repair sites to the nuclear periphery before Rad51 recruitment and strand invasion. How this mobilization initiates is unknown, and the contribution of phase separation to these dynamics is unclear. Here, we show that Nup98 nucleoporin is recruited to heterochromatic repair sites before relocalization through Sec13 or Nup88 nucleoporins, and downstream from the Smc5/6 complex and SUMOylation. Remarkably, the phase separation properties of Nup98 are required and sufficient to mobilize repair sites and exclude Rad51, thus preventing aberrant recombination while promoting HR repair. Disrupting this pathway results in heterochromatin repair defects and widespread chromosome rearrangements, revealing a novel "off-pore" role for nucleoporins and phase separation in nuclear dynamics and genome integrity in a multicellular eukaryote.

Genome Characterization and Phylogenetic Analysis of Scale Drop Disease Virus Isolated from Asian Seabass (Lates calcarifer).

Scale drop disease virus (SDDV), a double-stranded DNA virus in the family Iridoviridae, has been reported widely in southeast Asian countries as a causative agent of scale drop syndrome (SDS) in Asian seabass. SDS has resulted in high mortality and significant economic losses to the aquaculture industry. This study demonstrated the use of metagenomic methods to investigate bacterial and viral communities present in infected fish tissues and recover a complete genome of the causative agent named SDDV TH7_2019. Characterization of the TH7_2019 genome revealed a genome size of 131 kb with 134 putative ORFs encoding viral proteins potentially associated with host apoptosis manipulation. A comparative genome analysis showed a high degree of amino acid identity across SDDV strains, with variations in number of repeat sequences and mutations within core genes. Phylogenetic analyses indicate a close relationship among SDDV genomes. This research enhances our understanding of the genetic diversity and evolutionary relationship of SDDV, contributing valuable insights for further development of effective control strategies of SDDV.

GWarrange: a pre- and post- genome-wide association studies pipeline for detecting phenotype-associated genome rearrangement events.

The use of k-mers to capture genetic variation in bacterial genome-wide association studies (bGWAS) has demonstrated its effectiveness in overcoming the plasticity of bacterial genomes by providing a comprehensive array of genetic variants in a genome set that is not confined to a single reference genome. However, little attempt has been made to interpret k-mers in the context of genome rearrangements, partly due to challenges in the exhaustive and high-throughput identification of genome structure and individual rearrangement events. Here, we present GWarrange, a pre- and post-bGWAS processing methodology that leverages the unique properties of k-mers to facilitate bGWAS for genome rearrangements. Repeat sequences are common instigators of genome rearrangements through intragenomic homologous recombination, and they are commonly found at rearrangement boundaries. Using whole-genome sequences, repeat sequences are replaced by short placeholder sequences, allowing the regions flanking repeats to be incorporated into relatively short k-mers. Then, locations of flanking regions in significant k-mers are mapped back to complete genome sequences to visualise genome rearrangements. Four case studies based on two bacterial species (Bordetella pertussis and Enterococcus faecium) and a simulated genome set are presented to demonstrate the ability to identify phenotype-associated rearrangements. GWarrange is available at https://github.com/DorothyTamYiLing/GWarrange.

GenRepAI: Utilizing Artificial Intelligence to Identify Repeats in Genomic Suffix Trees

Background: The human genome is densely populated with repetitive DNA sequences that play crucial roles in genomic functions and structures but are also implicated in over 40 human diseases. The computational challenge of identifying and characterizing these repeats is significant due to the complexity and size of the genome, which are overwhelming traditional algorithms. Methods: To address these challenges, we propose GenRepAI, a deep learning framework to navigate and analyze genomic suffix trees. GenRepAI employs supervised machine learning classifiers trained on labeled datasets of repeat annotations and unsupervised anomaly detection to identify novel repeat sequences. The models are trained using convolutional neural networks (CNNs), long short-term memory networks (LSTMs), and vision transformers to classify and annotate repeats within the human genome. Results: GenRepAI is designed to comprehensively profile repeats that underlie various neurological diseases, allowing researchers to identify pathogenic expansions. The framework will integrate into existing genomic analysis pipelines, with the capability to screen patient genomes and highlight potential causal variants for further validation. Conclusion: GenRepAI is set to become a foundational tool in genomics, leveraging artificial intelligence to enhance the characterization of repetitive sequences. It promises significant advancements in the molecular diagnosis of repeat expansion disorders and contributes to a deeper understanding of genomic structure and function, with broad applications in personalized medicine.

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.

Complete Chloroplast Genome of Megacarpaea megalocarpa and Comparative Analysis with Related Species from Brassicaceae.

Megacarpaea megalocarpa, a perennial herbaceous species belonging to the Brassicaceae family, has potential medicinal value. We isolated and characterized the chloroplast (cp) genome of M. megalocarpa and compared it with closely related species. The chloroplast genome displayed a typical quadripartite structure, spanning 154,877 bp, with an overall guanine-cytosine (GC) content of 36.20%. Additionally, this genome contained 129 genes, 105 simple sequence repeats (SSRs), and 48 long repeat sequences. Significantly, the ycf1 gene exhibited a high degree of polymorphism at the small single copy (SSC) region and the inverted repeat a (IRa) boundary. Despite this polymorphism, relative synonymous codon usage (RSCU) values were found to be similar across species, and no large segment rearrangements or inversions were detected. The large single copy (LSC) and SSC regions showed higher sequence variations and nucleotide polymorphisms compared to the IR region. Thirteen distinct hotspot regions were identified as potential molecular markers. Our selection pressure analysis revealed that the protein-coding gene rpl20 is subjected to different selection pressures in various species. Phylogenetic analysis positioned M. megalocarpa within the expanded lineage II of the Brassicaceae family. The estimated divergence time suggests that M. megalocarpa diverged approximately 4.97 million years ago. In summary, this study provides crucial baseline information for the molecular identification, phylogenetic relationships, conservation efforts, and utilization of wild resources in Megacarpaea.

Characterization of a MHYT domain-coupled transcriptional regulator that responds to carbon monoxide.

The MHYT domain, identified over two decades ago for its potential to detect diatomic gases like CO, O2or NO, has awaited experimental validation as a protein sensory domain. Here, we characterize the MHYT domain-containing transcriptional regulator CoxC, which governs the expression of the cox genes responsible for aerobic CO oxidation in the carboxidotrophic bacterium Afipia carboxidovorans OM5. The C-terminal LytTR-type DNA-binding domain of CoxC binds to an operator region consisting of three direct repeats sequences overlapping the -35 box at the target PcoxB promoter, which is consistent with the role of CoxC as a specific transcriptional repressor of the cox genes. Notably, the N-terminal transmembrane MHYT domain endows CoxC with the ability to sense CO as an effector molecule, as demonstrated by the relief of CoxC-mediated repression and binding to the PcoxB promoter upon CO exposure. Furthermore, copper serves as the essential divalent cation for the interaction of CO with CoxC, thereby confirming previous hypothesis regarding the role of copper in the gas-sensing mechanism of MHYT domains. CoxC represents the prototype of a novel subfamily of single-component LytTR transcriptional regulators, characterized by the fusion of a DNA-binding domain with a membrane-bound MHYT sensor domain.

Phylogenomic analysis and dynamic evolution of chloroplast genomes of Clematis nannophylla

Clematis nannophylla is a perennial shrub of Clematis with ecological, ornamental, and medicinal value, distributed in the arid and semi-arid areas of northwest China. This study successfully determined the chloroplast (cp) genome of C. nannophylla, reconstructing a phylogenetic tree of Clematis. This cp genome is 159,801 bp in length and has a typical tetrad structure, including a large single-copy, a small single-copy, and a pair of reverse repeats (IRa and IRb). It contains 133 unique genes, including 89 protein-coding, 36 tRNA, and 8 rRNA genes. Additionally, 66 simple repeat sequences, 50 dispersed repeats, and 24 tandem repeats were found; many of the dispersed and tandem repeats were between 20–30 bp and 10–20 bp, respectively, and the abundant repeats were located in the large single copy region. The cp genome was relatively conserved, especially in the IR region, where no inversion or rearrangement was observed, further revealing that the coding regions were more conserved than the noncoding regions. Phylogenetic analysis showed that C. nannophylla is more closely related to C. fruticosa and C. songorica. Our analysis provides reference data for molecular marker development, phylogenetic analysis, population studies, and cp genome processes to better utilise C. nannophylla.

Repetitive Sequences, Codon Usage Bias and Phylogenetic Analysis of the Plastome of Miliusa glochidioides.

Annonaceae is the largest family in Magnoliales, exhibiting the greatest diversity among and within genera. In this study, we conducted an analysis of repetitive sequences and codon usage bias in the previously acquired plastome of Miliusa glochidioides. Using a concatenated dataset of shared genes, we constructed the phylogenetic relationships among 27 Annonaceae species. The results showed that the size of the plastomes in the Annonaceae ranged from 159 to 202kb, with the size of the inverted repeat region ranging from 40 to 65kb. Within the plastome of M. glochidioides, we identified 42 SSRs, 36 tandem repeats, and 9 dispersed repeats. These SSRs consist of three nucleotide types and eight motif types, with a preference for A/T bases, primarily located in the large single-copy regions and intergenic spacers. Tandem and dispersed repeat sequences were predominantly detected in the IR region. Through codon usage bias analysis, we identified 30 high-frequency codons and 11 optimal codons. The plastome of M. glochidioides demonstrated relatively weak codon usage bias, favoring codons with A/T endings, primarily influenced by natural selection. Phylogenetic analysis revealed that all four subfamilies formed monophyletic groups, with Cananga odorata (Ambavioideae) and Anaxagorea javanica (Anaxagoreoideae) successively nested outside Annonoideae + Malmeoideae. These findings improve our understanding of the plastome of M. glochidioides and provide additional insights for studying plastome evolution in Annonaceae.

Potential to use of viral reverse transcriptase inhibitors in oncology

In preparing the review, articles on the functioning of the reverse transcriptase enzyme of endogenous repeat sequences LINE1, the mechanisms of action and antitumor activity of viral reverse transcriptase inhibitors. Articles available in the biomedical literature information databases SciVerse Scopus, PubMed, Web of Science, Russian Science CitationIndex (RSCI) were analyzed. The review used information from 140 publications, of which 95 and 39 were published, respectively, over the last ten and three years, 2 articles present the results of clinical studies, and 45 articles refer to results demonstrating the anticancer properties of the studied compounds in various models in vitro and in vivo. Aim. Based on data on the functional properties of the reverse transcriptase enzyme of endogenous repeat sequences LINE1 (long interspersed nuclear elements 1), analyze the potential use of viral reverse transcriptase inhibitors in oncology, presenting their classification and main mechanisms of action. About 98 % of the human genome consists of repetitive sequences, most of which are represented by mobil genetic elements, the activation of which leads to increased genome instability. These include long (LINE) and short (SINE) interspersed nuclear element repeated DNA sequences interspersed nuclear elements, respectively, which occupy about 45 % of the human genome. Increased expression levels of these sequences in the genome have been identified in many forms of malignant neoplasms. Their transposition occurs due to the expression of LINE1-encoded reverse transcriptase, whichis homologous to viral reverse transcriptase. To date, reverse transcriptase inhibitors of viruses of nucleoside and non-nucleoside structure have been developed and are successfully used in the clinic. These drugs demonstrate an inhibitory effect on both LINE1 reverse transcriptase and telomerase, which provides the tumor cell with the ability to overcome replicative senescence. Due to these properties, these compounds are expected to exhibit both their own antitumor activity and increase the sensitivity of tumor cells to the therapy of malignant neoplasms, which is experimentally confirmed in models of malignant tumors in vitro and in vivo. Use of reverse transcriptase inhibitors in combination therapy seems advisable both to prevent further genome rearrangements caused by LINE1 and to suppress the survival of tumor cells by inhibiting telomerase activity.

Molecular evolution of the mammalian kinetochore complex.

Mammalian centromeres are satellite-rich chromatin domains that serve as sites for kinetochore complex assembly. Centromeres are highly variable in sequence and satellite organization across species, but the processes that govern the co-evolutionary dynamics between rapidly evolving centromeres and their associated kinetochore proteins remain poorly understood. Here, we pursue a course of phylogenetic analyses to investigate the molecular evolution of the complete kinetochore complex across primate and rodent species with divergent centromere repeat sequences and features. We show that many protein components of the core centromere associated network (CCAN) harbor signals of adaptive evolution, consistent with their intimate association with centromere satellite DNA and roles in the stability and recruitment of additional kinetochore proteins. Surprisingly, CCAN and outer kinetochore proteins exhibit comparable rates of adaptive divergence, suggesting that changes in centromere DNA can ripple across the kinetochore to drive adaptive protein evolution within distant domains of the complex. Our work further identifies kinetochore proteins subject to lineage-specific adaptive evolution, including rapidly evolving proteins in species with centromere satellites characterized by higher-order repeat structure and lacking CENP-B boxes. Thus, features of centromeric chromatin beyond the linear DNA sequence may drive selection on kinetochore proteins. Overall, our work spotlights adaptively evolving proteins with diverse centromere-associated functions, including centromere chromatin structure, kinetochore protein assembly, kinetochore-microtubule association, cohesion maintenance, and DNA damage response pathways. These adaptively evolving kinetochore protein candidates present compelling opportunities for future functional investigations exploring how their concerted changes with centromere DNA ensure the maintenance of genome stability.

A self-complementary AAV proviral plasmid that reduces cross-packaging and ITR promoter activity in AAV vector preparations

Adeno-associated viral vectors (AAV) are a leading delivery system for gene therapy in animal models and humans. With several FDA-approved AAV gene therapies on the market, issues related to vector manufacturing have become increasingly important. In this study, we focused on potentially toxic DNA contaminants that can arise from AAV proviral plasmids, the raw materials required for manufacturing recombinant AAV in eukaryotic cells. Typical AAV proviral plasmids are circular DNAs containing a therapeutic gene cassette flanked by natural AAV inverted terminal repeat (ITR) sequences, and a plasmid backbone carrying prokaryotic sequences required for plasmid replication and selection in bacteria. While the majority of AAV particles package the intended therapeutic payload, some capsids instead package the bacterial sequences located on the proviral plasmid backbone. Since ITR sequences also have promoter activity, potentially toxic bacterial open reading frames can be produced in vivo, thereby representing a safety risk. In this study, we describe a new AAV proviral plasmid for vector manufacturing that (1) significantly decreases cross-packaged bacterial sequences; (2) increases correctly packaged AAV payloads; and (3) blunts ITR-driven transcription of cross-packaged material to avoid expressing potentially toxic bacterial sequences. This system may help improve the safety of AAV vector products.

The Whole Mitochondrial Genome Sequence of Dendrobium loddigesii Rolfe, an Endangered Orchid Species in China, Reveals a Complex Multi-Chromosome Structure.

Dendrobium loddigesii is a precious traditional Chinese medicine with high medicinal and ornamental value. However, the characterization of its mitochondrial genome is still pending. Here, we assembled the complete mitochondrial genome of D. loddigesii and discovered that its genome possessed a complex multi-chromosome structure. The mitogenome of D. loddigesii consisted of 17 circular subgenomes, ranging in size from 16,323 bp to 56,781 bp. The total length of the mitogenome was 513,356 bp, with a GC content of 43.41%. The mitogenome contained 70 genes, comprising 36 protein-coding genes (PCGs), 31 tRNA genes, and 3 rRNA genes. Furthermore, we detected 403 repeat sequences as well as identified 482 RNA-editing sites and 8154 codons across all PCGs. Following the sequence similarity analysis, 27 fragments exhibiting homology to both the mitogenome and chloroplast genome were discovered, accounting for 9.86% mitogenome of D. loddigesii. Synteny analysis revealed numerous sequence rearrangements in D. loddigesii and the mitogenomes of related species. Phylogenetic analysis strongly supported that D. loddigesii and D. Amplum formed a single clade with 100% bootstrap support. The outcomes will significantly augment the orchid mitochondrial genome database, offering profound insights into Dendrobium's intricate mitochondrial genome architecture.

Chromosome-level genome of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus

Recent conservation efforts to protect rare and endangered aquatic species have intensified. Nevertheless, the ornate spiny lobster (Panulirus ornatus), which is prevalent in the Indo-Pacific waters, has been largely ignored. In the absence of a detailed genomic reference, the conservation and population genetics of this crustacean are poorly understood. Here, We assembled a comprehensive chromosome-level genome for P. ornatus. This genome—among the most detailed for lobsters—spans 2.65 Gb with a contig N50 of 51.05 Mb, and 99.11% of the sequences with incorporated to 73 chromosomes. The ornate spiny lobster genome comprises 65.67% repeat sequences and 22,752 protein-coding genes with 99.20% of the genes functionally annotated. The assembly of the P. ornatus genome provides valuable insights into comparative crustacean genomics and endangered species conservation, and lays the groundwork for future research on the speciation, ecology, and evolution of the ornate spiny lobster.

Differential amplification and contraction of satellite DNAs in the distinct lineages of the beetle Euchroma gigantea

Satellite DNA (satDNA) consists of tandem repeat sequences that typically evolve rapidly through evolutionary mechanisms, including unequal crossover, transposition events, and others. The evolutionary history of Euchroma gigantea is marked by complex chromosomal evolution between lineages, making this species an interesting model for understanding satDNA evolution at intraspecies level. Therefore, our aim was to comprehend the potential contribution of satDNAs to the greater chromosomal differentiation of evolutionary lineages in E. gigantea by investigating the differential patterns of amplification and contraction of the repeats. To achieve this, we employed de novo identification of satDNA using RepeatExplorer and TAREAN, allowing the satellitome characterization between lineages. A total of 26 satDNA families were identified, ranging from 18 to 1101 nucleotides in length, with most families being shared between individuals/lineages, as predicted by the library hypothesis, except for the satDNA EgiSat21-168 that was absent for Northeast Lineage. The total satDNA content of the individuals was less than 11.2%, and it appeared to increase in two directions following the chromosomal evolution model. Thirteen satDNAs exhibited different patterns of amplification, and nine ones were contracted among individuals. Additionally, most repeats showed a divergence of about 10% for these satDNAs, indicating satellitome differentiation for each lineage/individual. This scenario suggests that the expansion of the satellitome occurred differentially among individuals/lineages of E. gigantea, with the contribution of various DNA turnover mechanisms after geographical isolation, and that they could be involved with karyotype evolution.

Generation of genetically modified Friedreich’s ataxia induced pluripotent stem cell lines and isogenic control lines carrying an inducible neurogenin-2 expression cassette

Friedreich’s ataxia (FRDA) is a rare neurodegenerative disease caused by an expansion of a GAA repeat sequence within the Frataxin (FXN) gene. Prominent regions of neurodegeneration include sensory neurons within the dorsal root ganglia. Here we present a set of genetically modified FRDA induced pluripotent stem cell (iPSC) lines that carry an inducible neurogenin-2 (NGN2) expression cassette. Exogenous expression of NGN2 in iPSC derived neural crest progenitors efficiently generates functionally mature sensory neurons. These cell lines will provide a streamlined source of FRDA iPSC sensory neurons for studying both disease mechanism and screening potential therapeutics.

Construction of oral squamous cell carcinoma organoids invitro 3D-culture for drug screening.

Patient-derived organoids are potent pre-chemotherapy models. Due to limited research on diverse types of oral squamous cell carcinoma (OSCC) and construction efficiency, our goal was to optimize OSCC organoid models from various sites and assess drug responsiveness. We screened and optimized culture media, employing three-dimensional techniques to construct human-derived oral squamous cell carcinoma (OSCC) organoid models invitro. Morphological validation, immunofluorescence analysis, tissue origin verification, and Short Tandem Repeat (STR) sequencing confirmed the consistency between organoids and source tissues. These organoid models were then subjected to varying concentrations of anticancer drugs, with subsequent assessment of cell viability to calculate IC50 values. Twenty-nine surgical specimens yielded an 86.2% success rate in culturing 25 organoids invitro. Morphological consistency confirmed nuclear atypia and positive expression of K5, P40, and E-cadherin, indicating squamous epithelial origin. Cultured complex organoids included α-SMA+ tumour-associated fibroblasts and tumour stem cells expressing CD44 and Ki67. STR sequencing affirmed genomic homogeneity between cultured organoids and source tissues. Drug sensitivity testing revealed diverse responses among organoids, highlighting their value for assessing drug sensitivity. An efficient OSCC organoid culture system for personalized invitro drug sensitivity screening was established, laying the foundation for precise treatment development.

Characterization and the first complete genome sequence of a novel strain of Bergeyella porcorum isolated from pigs in China

BackgroundBergeyella porcorum is a newly identified bacterium that has an ambiguous relationship with pneumonia in pigs. However, few studies have adequately characterized this species.ResultsIn this study, we analyzed the morphological, physiological, and genomic characteristics of the newly identified B. porcorum sp. nov. strain QD2021 isolated from pigs. The complete genome sequence of the B. porcorum QD2021 strain consists of a single circular chromosome (2,271,736 bp, 38.51% G + C content), which encodes 2,578 genes. One plasmid with a size of 70,040 bp was detected. A total of 121 scattered repeat sequences, 319 tandem repeat sequences, 4 genomic islands, 5 prophages, 3 CRISPR sequences, and 51 ncRNAs were predicted. The coding genes of the B. porcorum genome were successfully annotated across eight databases (NR, GO, KEGG, COG, TCDB, Pfam, Swiss-Prot and CAZy) and four pathogenicity-related databases (PHI, CARD, VFDB and ARDB). In addition, a comparative genome analysis was performed to explore the evolutionary relationships of B. porcorum QD2021.ConclusionsTo our knowledge, this is the first study to provide fundamental phenotypic and whole-genome sequences for B. porcorum. Our results extensively expand the current knowledge and could serve as a valuable genomic resource for future research on B. porcorum.

Insights into structure, codon usage, repeats, and RNA editing of the complete mitochondrial genome of Perilla frutescens (Lamiaceae)

Perilla frutescens (L.) Britton, a member of the Lamiaceae family, stands out as a versatile plant highly valued for its unique aroma and medicinal properties. Additionally, P. frutescens seeds are rich in Îś-linolenic acid, holding substantial economic importance. While the nuclear and chloroplast genomes of P. frutescens have already been documented, the complete mitochondrial genome sequence remains unreported. To this end, the sequencing, annotation, and assembly of the entire Mitochondrial genome of P. frutescens were hereby conducted using a combination of Illumina and PacBio data. The assembled P. frutescens mitochondrial genome spanned 299,551 bp and exhibited a typical circular structure, involving a GC content of 45.23%. Within the genome, a total of 59 unique genes were identified, encompassing 37 protein-coding genes, 20 tRNA genes, and 2 rRNA genes. Additionally, 18 introns were observed in 8 protein-coding genes. Notably, the codons of the P. frutescens mitochondrial genome displayed a notable A/T bias. The analysis also revealed 293 dispersed repeat sequences, 77 simple sequence repeats (SSRs), and 6 tandem repeat sequences. Moreover, RNA editing sites preferentially produced leucine at amino acid editing sites. Furthermore, 70 sequence fragments (12,680 bp) having been transferred from the chloroplast to the mitochondrial genome were identified, accounting for 4.23% of the entire mitochondrial genome. Phylogenetic analysis indicated that among Lamiaceae plants, P. frutescens is most closely related to Salviamiltiorrhiza and Platostoma chinense. Meanwhile, inter-species Ka/Ks results suggested that Ka/Ks <1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$<1$$\\end{document} for 28 PCGs, indicating that these genes were evolving under purifying selection. Overall, this study enriches the mitochondrial genome data for P. frutescens and forges a theoretical foundation for future molecular breeding research.

The chloroplast genomes of two medicinal species (Veronica anagallis-aquatica L. and Veronica undulata Wall.) and its comparative analysis with related Veronica species

Veronica anagallis-aquatica L. and Veronica undulata Wall. are widely used ethnomedicinal plants in China. The two species have different clinical efficacies, while their extremely similar morphology and unclear interspecific relationship make it difficult to accurately identify them, leading to increased instances of mixed usage. This article reports on the complete chloroplast genomes sequence of these two species and their related Veronica species to conduct a comparative genomics analysis and phylogenetic construction. The results showed that the chloroplast (cp) genomes of Veronica exhibited typical circular quadripartite structures, with total lengths of 149,386 to 152,319 base pairs (bp), and GC content of 37.9 to 38.1%, and the number of genes was between 129–134. The total number of simple sequence repeats (SSRs) in V. anagallis-aquatica and V. undulata is 37 and 36, while V. arvensis had the highest total number of 56, predominantly characterized by A/T single bases. The vast majority of long repeat sequence types are forward repeats and palindromic repeats. Selective Ka/Ks values showed that three genes were under positive selection. Sequence differences often occur in the non-coding regions of the large single-copy region (LSC) and small single-copy region (SSC), with the lowest sequence variation in the inverted repeat regions (IR). Seven highly variable regions (trnT-GGU-psbD, rps8-rpl16, trnQ-UUG, trnN-GUU-ndhF, petL, ycf3, and ycf1) were detected, which may be potential molecular markers for identifying V. anagallis-aquatica and V. undulata. The phylogenetic tree indicates that there is a close genetic relationship between the genera Veronica and Neopicrorhiza, and V. anagallis-aquatica and V. undulata are sister groups. The molecular clock analysis results indicate that the divergence time of Veronica may occur at ∼ 9.09 Ma, and the divergence time of these two species occurs at ∼ 0.48 Ma. It is speculated that climate change may be the cause of Veronica species diversity and promote the radiation of the genus. The chloroplast genome data of nine Veronica specie provides important insights into the characteristics and evolution of the chloroplast genome of this genus, as well as the phylogenetic relationships of the genus Veronica.