Recombinant Genomes Research Articles

Evolution of homologous recombination rates across bacteria

Bacteria are nonsexual organisms but are capable of exchanging DNA at diverse degrees through homologous recombination. Intriguingly, the rates of recombination vary immensely across lineages where some species have been described as purely clonal and others as "quasi-sexual." However, estimating recombination rates has proven a difficult endeavor and estimates often vary substantially across studies. It is unclear whether these variations reflect natural variations across populations or are due to differences in methodologies. Consequently, the impact of recombination on bacterial evolution has not been extensively evaluated and the evolution of recombination rate-as a trait-remains to be accurately described. Here, we developed an approach based on Approximate Bayesian Computation that integrates multiple signals of recombination to estimate recombination rates. We inferred the rate of recombination of 162 bacterial species and one archaeon and tested the robustness of our approach. Our results confirm that recombination rates vary drastically across bacteria; however, we found that recombination rate-as a trait-is conserved in several lineages but evolves rapidly in others. Although some traits are thought to be associated with recombination rate (e.g., GC-content), we found no clear association between genomic or phenotypic traits and recombination rate. Overall, our results provide an overview of recombination rate, its evolution, and its impact on bacterial evolution.

Genomic Diversity and Recombination Analysis of the Spike Protein Gene from Selected Human Coronaviruses.

Human coronaviruses (HCoVs) are seriously associated with respiratory diseases in humans and animals. The first human pathogenic SARS-CoV emerged in 2002-2003. The second was MERS-CoV, reported from Jeddah, the Kingdom of Saudi Arabia, in 2012, and the third one was SARS-CoV-2, identified from Wuhan City, China, in late December 2019. The HCoV-Spike (S) gene has the highest mutation/insertion/deletion rate and has been the most utilized target for vaccine/antiviral development. In this manuscript, we discuss the genetic diversity, phylogenetic relationships, and recombination patterns of selected HCoVs with emphasis on the S protein gene of MERS-CoV and SARS-CoV-2 to elucidate the possible emergence of new variants/strains of coronavirus in the near future. The findings showed that MERS-CoV and SARS-CoV-2 have significant sequence identity with the selected HCoVs. The phylogenetic tree analysis formed a separate cluster for each HCoV. The recombination pattern analysis showed that the HCoV-NL63-Japan was a probable recombinant. The HCoV-NL63-USA was identified as a major parent while the HCoV-NL63-Netherland was identified as a minor parent. The recombination breakpoints start in the viral genome at the 142 nucleotide position and end at the 1082 nucleotide position with a 99% CI and Bonferroni-corrected p-value of 0.05. The findings of this study provide insightful information about HCoV-S gene diversity, recombination, and evolutionary patterns. Based on these data, it can be concluded that the possible emergence of new strains/variants of HCoV is imminent.

Data-driven recombination detection in viral genomes.

Recombination is a key molecular mechanism for the evolution and adaptation of viruses. The first recombinant SARS-CoV-2 genomes were recognized in 2021; as of today, more than ninety SARS-CoV-2 lineages are designated as recombinant. In the wake of the COVID-19 pandemic, several methods for detecting recombination in SARS-CoV-2 have been proposed; however, none could faithfully confirm manual analyses by experts in the field. We hereby present RecombinHunt, an original data-driven method for the identification of recombinant genomes, capable of recognizing recombinant SARS-CoV-2 genomes (or lineages) with one or two breakpoints with high accuracy and within reduced turn-around times. ReconbinHunt shows high specificity and sensitivity, compares favorably with other state-of-the-art methods, and faithfully confirms manual analyses by experts. RecombinHunt identifies recombinant viral genomes from the recent monkeypox epidemic in high concordance with manually curated analyses by experts, suggesting that our approach is robust and can be applied to any epidemic/pandemic virus.

Microsatellite break-induced replication generates highly mutagenized extrachromosomal circular DNAs.

Extrachromosomal circular DNAs (eccDNAs) are produced from all regions of the eucaryotic genome. We used inverse PCR of non-B microsatellites capable of forming hairpin, triplex, quadruplex and AT-rich structures integrated at a common ectopic chromosomal site to show that these non-B DNAs generate highly mutagenized eccDNAs by replication-dependent mechanisms. Mutagenesis occurs within the non-B DNAs and extends several kilobases bidirectionally into flanking and nonallelic DNA. Each non-B DNA exhibits a different pattern of mutagenesis, while sister clones containing the same non-B DNA also display distinct patterns of recombination, microhomology-mediated template switching and base substitutions. Mutations include mismatches, short duplications, long nontemplated insertions, large deletions and template switches to sister chromatids and nonallelic chromosomes. Drug-induced replication stress or the depletion of DNA repair factors Rad51, the COPS2 signalosome subunit or POLη change the pattern of template switching and alter the eccDNA mutagenic profiles. We propose an asynchronous capture model based on break-induced replication from microsatellite-induced DNA double strand breaks to account for the generation and circularization of mutagenized eccDNAs and the appearance of genomic homologous recombination deficiency (HRD) scars. These results may help to explain the appearance of tumor eccDNAS and their roles in neoantigen production, oncogenesis and resistance to chemotherapy.

Integration and the risk of liver cancer-Is there a real risk?

Adeno-associated virus (AAV)-based gene therapies are in clinical development for haemophilia and other genetic diseases. Since the recombinant AAV genome primarily remains episomal, it provides the opportunity for long-term expression in tissues that are not proliferating and reduces the safety concerns compared with integrating viral vectors. However, AAV integration events are detected at a low frequency. Preclinical studies in mouse models have reported hepatocellular carcinoma (HCC) after systemic AAV administration in some settings, though this has not been reported in large animal models. The risk of HCC or other cancers after AAV gene therapy in clinical studies thus remains theoretical. Potential risk factors for HCC after gene therapy are beginning to be elucidated through animal studies, but their relevance to human studies remains unknown. Studies to investigate the factors that may influence the risk of oncogenesis as well as detailed investigation of cases of cancer in AAV gene therapy patients will be important to define the potential risk of AAV genotoxicity.

Near-full-length genome analysis of two novel HIV second recombinant forms in Hebei, China.

The number of individuals infected with HIV-1 among men who have sex with men (MSM) has risen rapidly in recent years in China, and the subtypes CRF01_AE, CRF07_BC, and B, as well as many novel unique recombinant forms (URFs) are prevalent among them. Co-circulation of strains among MSM populations allows the generation of circulating recombinant forms (CRFs) and URFs. In this study, we identified two new URFs from two HIV-1-positive subjects who were infected through homosexual contact in Hebei, China. Analysis of near-full-length genome sequences, using phylogenetic and recombination analysis showed that the two URFs originated from CRF01_AE, CRF07_BC, and B, and CRF01_AE segments in the backbone of the URFs were derived from cluster 4 of CRF01_AE. The CRF07_BC segments of two URFs were clustered with 07BC_N in a phylogenetic tree. The identification of novel URFs with complex genomic structures shows that it is necessary to strengthen surveillance of HIV-1 variants in MSM populations in this region.

Assembly and characterization of the complete mitochondrial genome of Descurainia sophia: a medical plant producing glucosinolate

As a commonly used traditional Chinese medicine, Descurainia sophia has various reported pharmacological activities such as prevention of chronic diseases, cardiovascular protection, neuroprotection, anticancer, antioxidant, and has high nutritional values, but its mitochondrial genome (mitogenome) has not been reported. In this study, the D. sophia mitogenome was sequenced, assembled and annotated using next generation sequencing technologies on the Illumina HiSeq and PacBio Sequel platform. The results showed that the D. sophia mitogenome was a circular DNA molecule with the length of 265,457 bp and the GC content of 44.78%, encoding unique 56 genes, including 35 protein coding genes, 18 tRNA genes and 3 rRNA genes. We analyzed codon preference, repeated sequences, DNA sequence transformation, phylogenetic relationship, RNA editing and synteny in the D. sophia mitogenome. Based on the DNA sequences of mitogenomes, the phylogenetic and synteny analysis indicated that D. sophia was more closely related to Boechera stricta, and the mitogenome of D. sophia underwent frequent genome recombination. In this study, we obtained the complete mitogenome of D. sophia, which can provide valuable data support for future studies on genetic characteristics, phylogenetic relationships, species identification and development of new molecular markers in D. sophia.

Identification and evolution of ICE-PmuST394: a novel integrative conjugative element in Pasteurella multocida ST394.

The emergence of macrolide and tetracycline resistance within Pasteurella multocida isolated from feedlot cattle and the dominance of ST394 in Australia was reported recently. To establish the genetic context of the resistance genes in P. multocida 17BRD-035, the ST394 reference genome, and conduct a molecular risk assessment of their ability to disperse laterally. A bioinformatic analysis of the P. multocida 17BRD-035 genome was conducted to determine if integrative conjugative elements (ICEs) carrying resistance genes, which hamper antibiotic treatment options locally, are in circulation in Australian feedlots. A novel element, ICE-PmuST394, was characterized in P. multocida 17BRD-035. It was also identified in three other isolates (two ST394s and a ST125) in Australia and is likely present in a genome representing P. multocida ST79 from the USA. ICE-PmuST394 houses a resistance module carrying two variants of the blaROB gene, blaROB-1 and blaROB-13, and the macrolide esterase gene, estT. The resistance gene combination on ICE-PmuST394 confers resistance to ampicillin and tilmicosin, but not to tulathromycin and tildipirosin. Our analysis suggests that ICE-PmuST394 is circulating both by clonal expansion and horizontal transfer but is currently restricted to a single feedlot in Australia. ICE-PmuST394 carries a limited number of unusual antimicrobial resistance genes but has hotspots that facilitate genomic recombination. The element is therefore amenable to hosting more resistance genes, and therefore its presence (or dispersal) should be regularly monitored. The element has a unique molecular marker, which could be exploited for genomic surveillance purposes locally and globally.

Small RNAs from mitochondrial genome recombination sites are incorporated into T. gondii mitoribosomes.

The mitochondrial genomes of apicomplexans comprise merely three protein-coding genes, alongside a set of thirty to forty genes encoding small RNAs (sRNAs), many of which exhibit homologies to rRNA from E. coli. The expression status and integration of these short RNAs into ribosomes remains unclear and direct evidence for active ribosomes within apicomplexan mitochondria is still lacking. In this study, we conducted small RNA sequencing on the apicomplexan Toxoplasma gondii to investigate the occurrence and function of mitochondrial sRNAs. To enhance the analysis of sRNA sequencing outcomes, we also re-sequenced the T. gondii mitochondrial genome using an improved organelle enrichment protocol and Nanopore sequencing. It has been established previously that the T. gondii genome comprises 21 sequence blocks that undergo recombination among themselves but that their order is not entirely random. The enhanced coverage of the mitochondrial genome allowed us to characterize block combinations at increased resolution. Employing this refined genome for sRNA mapping, we find that many small RNAs originated from the junction sites between protein-coding blocks and rRNA sequence blocks. Surprisingly, such block border sRNAs were incorporated into polysomes together with canonical rRNA fragments and mRNAs. In conclusion, apicomplexan ribosomes are active within polysomes and are indeed assembled through the integration of sRNAs, including previously undetected sRNAs with merged mRNA-rRNA sequences. Our findings lead to the hypothesis that T. gondii's block-based genome organization enables the dual utilization of mitochondrial sequences as both messenger RNAs and ribosomal RNAs, potentially establishing a link between the regulation of rRNA and mRNA expression.

The development of a rapid, high-throughput neutralization assay using a SARS-CoV-2 reporter

Many methods have been developed to measure the neutralizing capacity of antibodies to SARS-CoV-2. However, these methods are low throughput and can be difficult to quickly modify in response to emerging variants. Therefore, an experimental system for rapid and easy measurement of the neutralizing capacity of antibodies against various variants is needed. In this study, we developed an experimental system that can efficiently measure the neutralizing capacity of sera by using a GFP-carrying recombinant SARS-CoV-2 with spike proteins of multiple variants (B.1.1, BA.5, or XBB.1.5). For all 3 recombinant chimeric genomes generated, neutralizing antibody titers determined by measuring GFP fluorescence intensity correlated significantly with those calculated from viral RNA levels measured by RT-qPCR in the supernatant of infected cells. Furthermore, neutralizing antibody titers determined by visually assessing GFP fluorescence using microscopy were also significantly correlated with those determined by RT-qPCR. By using this high-throughput method, it is now possible to quickly and easily determine the neutralizing capacity of antibodies against SARS-CoV-2 variants.

Genomic surveillance and evolutionary dynamics of type 2 porcine reproductive and respiratory syndrome virus in China spanning the African swine fever outbreak.

Porcine reproductive and respiratory syndrome virus (PRRSV) poses a serious threat to the pig industry in China. Our previous study demonstrated that PRRSV persists with local circulations and overseas imports in China and has formed a relatively stable epidemic pattern. However, the sudden African swine fever (ASF) outbreak in 2018 caused serious damage to China's pig industry structure, which resulted in about 40 per cent of pigs being slaughtered. The pig yields recovered by the end of 2019. Thus, whether the ASF outbreak reframed PRRSV evolution with changes in pig populations and further posed new threats to the pig industry becomes a matter of concern. For this purpose, we conducted genomic surveillance and recombination, NSP2 polymorphism, population dynamics, and geographical spread analysis of PRRSV-2, which is dominant in China. The results showed that the prevalence of ASF had no significant effects on genetic diversities like lineage composition, recombination patterns, and NSP2 insertion and deletion patterns but was likely to lead to changes in PRRSV-2 recombination frequency. As for circulation of the two major sub-lineages of Lineage 1, there was no apparent transmission of NADC30-like among provinces, while NADC34-like had obvious signs of inter-provincial transmission and foreign importation during the ASF epidemic. In addition, two suspected vaccine recombinant epidemic strains suggest a slight safety issue of vaccine use. Herein, the interference of ASF to the PRRSV-2 evolutionary pattern was evaluated and vaccine safety was analyzed, in order to monitor the potential threat of PRRSV-2 to China's pig industry in the post-epidemic era of ASF.

Genomic mutations and homologous recombination deficiency (HRD) score in Japanese patients with metastatic hormone-sensitive prostate cancer (mHSPC): The impact of intraductal carcinoma of the prostate (IDC-P).

Genomic mutations and homologous recombination deficiency (HRD) score in Japanese patients with metastatic hormone-sensitive prostate cancer (mHSPC): The impact of intraductal carcinoma of the prostate (IDC-P).

Real-world effectiveness of PARP inhibitors (PARPi) in metastatic castration-resistant prostate cancer (mCRPC) by genomic homologous recombination repair (HRR) alterations and homologous recombination deficiency signature (HRDsig).

186 Background: Multiple randomized controlled trials have demonstrated efficacy of poly ADP-ribose polymerase (PARP) inhibitors in patients with metastatic castration resistant prostate cancer (mCRPC) harboring deleterious BRCA1 or BRCA2 alterations (BRCAalt). Alterations in other DNA HRR gene pathways are common and the efficacy of PARPi in these populations is less clear. We sought to evaluate the real-world effectiveness of PARPi in commonly defined HRR groups and explore clinical validity of HRDsig to additionally predict outcomes on PARPi. Methods: This study used the nationwide (US-based) de-identified Flatiron Health-Foundation Medicine Prostate Cancer clinico-genomic database (FH-FMI CGDB), originating from approximately 280 US cancer clinics (~800 sites of care). Patients with mCRPC and tumor genomic profiling (tissue or liquid) who underwent single agent PARPi treatment were included, grouped by biomarkers with deleterious alterations detected: BRCA, ATM, other HRR ( ATR, ATRX, BAP1, BARD1, BRIP1, CHEK1/2, CDK12, FANCA, FANCL, MRE11, RAD51B/C/D, RAD54L, PALB2), or no HRR. Separately, patients with mCRPC and tumor tissue profiling were grouped by HRDsig status (+ vs -) using a machine learning- based algorithm built using a diverse set of copy number features. Kaplan Meier estimates and multivariable Cox PH models assessed time to next therapy (TTNT), time to treatment discontinuation (TTD), and real-world overall survival (rwOS). Results: 313 patients with mCRPC met inclusion criteria, of which 162 had tissue specimens for genomics evaluations and 151 had liquid profiling. In the tissue cohort, BRCAalt had more favorable TTNT (median 10.0 months, 95% CI 8.0 – 14.1) relative to ATM (median 3.8 months, 95% CI 2.7 – 6.9), other HRR (median 4.9 months, 95% CI 3.6 – 7.2) and no HRR (median 4.9, 95% CI 3.6 – 7.2 months). Other proxies for drug effectiveness (TTD, PSA kinetics) and rwOS had consistent results. More granular mutational analyses and liquid profiling data will be presented at the conference. Among evaluable patients without BRCAalt on tissue biopsy receiving PARPi, 5 of 75 (7%) were HRDsig(+) and had directionally more favorable (vs. HRDsig(-)) TTNT (HR: 0.39, 95%CI 0.14 – 1.12), TTD (0.46, 95%CI 0.16 – 1.31) and rwOS (HR: 0.27, 95%CI 0.07 – 1.14). Conclusions: We observe no significant outcomes difference between non- BRCAalt groups (defined as ATM, other HRR, and no HRR) in this cohort with respect to proxies of drug effectiveness. These results are largely consistent with biomarker-defined subgroup analyses from completed clinical trials. However, HRDsig may be able to identify a small non-BRCAalt subgroup with enhanced benefit. Associations of BRCAalt(-)/HRDsig(+) and PARPi performance deserve further attention in additional cohorts.

Linear DNA-driven recombination in mammalian mitochondria.

Mitochondrial DNA (mtDNA) recombination in animals has remained enigmatic due to its uniparental inheritance and subsequent homoplasmic state, which excludes the biological need for genetic recombination, as well as limits tools to study it. However, molecular recombination is an important genome maintenance mechanism for all organisms, most notably being required for double-strand break repair. To demonstrate the existence of mtDNA recombination, we took advantage of a cell model with two different types of mitochondrial genomes and impaired its ability to degrade broken mtDNA. The resulting excess of linear DNA fragments caused increased formation of cruciform mtDNA, appearance of heterodimeric mtDNA complexes and recombinant mtDNA genomes, detectable by Southern blot and by long range PacBio® HiFi sequencing approach. Besides utilizing different electrophoretic methods, we also directly observed molecular complexes between different mtDNA haplotypes and recombination intermediates using transmission electron microscopy. We propose that the known copy-choice recombination by mitochondrial replisome could be sufficient for the needs of the small genome, thus removing the requirement for a specialized mitochondrial recombinase. The error-proneness of this system is likely to contribute to the formation of pathological mtDNA rearrangements.

Detection of multiple recombinations of avian coronavirus in South Korea by whole-genome analysis

Infectious bronchitis virus (IBV), an avian coronavirus, has caused considerable damage to the poultry industry. In Korea, indigenous KM91-like and newly introduced QX-like lineages belonging to the GI-19 lineage have been prevalent despite constant vaccination. In this study, complete genome sequences of 23 IBV isolates in Korea from 2010 to 2020 were obtained using next-generation sequencing, and their phylogenetic relationship and recombination events were analyzed. Phylogenetic analysis based on the S1 gene showed that all isolates belonged to the GI-19 lineage and were divided into five subgroups (KM91-like, K40/09-like, and QX-like II to IV). Among the 23 isolates, 14 recombinants were found, including frequent recombination between KM91-like and QX-like strains. In addition, it was observed that other lineages, such as GI-1, GI-13, and GI-16, were involved in recombination. Most recombination breakpoints were detected in the ORF1ab gene, particularly nsp3. However, when considering the size of each genome, recombination occurred more frequently in the 3a, E and 5a genes. Taken together, genetic recombination frequently occurred throughout the entire genome between various IBV strains in Korea, including live attenuated vaccine strain. Our study suggests the necessity of further research on the contribution of recombination of genomes outside the spike region to the biological characteristics of IBV.

Two mitochondrial DNA polymorphisms modulate cardiolipin binding and lead to synthetic lethality

Genetic screens have been used extensively to probe interactions between nuclear genes and their impact on phenotypes. Probing interactions between mitochondrial genes and their phenotypic outcome, however, has not been possible due to a lack of tools to map the responsible polymorphisms. Here, using a toolkit we previously established in Drosophila, we isolate over 300 recombinant mitochondrial genomes and map a naturally occurring polymorphism at the cytochrome c oxidase III residue 109 (CoIII109) that fully rescues the lethality and other defects associated with a point mutation in cytochrome c oxidase I (CoIT300I). Through lipidomics profiling, biochemical assays and phenotypic analyses, we show that the CoIII109 polymorphism modulates cardiolipin binding to prevent complex IV instability caused by the CoIT300I mutation. This study demonstrates the feasibility of genetic interaction screens in animal mitochondrial DNA. It unwraps the complex intra-genomic interplays underlying disorders linked to mitochondrial DNA and how they influence disease expression.

Recombinant AAV genome size effect on viral vector production, purification, and thermostability

Adeno-associated virus (AAV) has shown great promise as a viral vector for gene therapy in clinical applications. This work studied the effect of the genome size on AAV production, purification, and thermostability by producing AAV2-GFP using suspension adapted HEK293 cells via triple transfection using AAV plasmids containing the same green fluorescent protein (GFP) transgene with DNA stuffers for variable size AAV genomes consisting of 1.9, 3.4, and 4.9 kb (ITR to ITR). Production was performed at the small and large shake flask scales and the results showed that the 4.9 kb GFP genome had significantly reduced encapsidation compared to other genomes. The large shake flask productions were purified by AEX chromatography and the results suggest that the triple transfection condition significantly impacts the AEX retention time and resolution between the full and empty capsid peaks. Charge detection mass spectrometry (CD-MS) was performed on all AEX full capsid peak samples showing a wide distribution of empty, partial, full length, and co-packaged DNA in the capsids. The AEX purified samples were then analyzed by differential scanning fluorimetry (DSF) and the results suggest that sample formulation may improve the thermostability of AAV genome ejection melting temperature regardless of the packaged genome content.

Unraveling the genetic evolution of SARS-CoV-2 Recombinants using mutational dynamics across the different lineages.

Recombination serves as a common strategy employed by RNA viruses for their genetic evolution. Extensive genomic surveillance during the COVID-19 pandemic has reported SARS-CoV-2 Recombinant strains indicating recombination events during the viral evolution. This study introspects the phenomenon of genome recombination by tracing the footprint of prominent lineages of SARS-CoV-2 at different time points in the context of on-going evolution and emergence of Recombinants. Whole genome sequencing was carried out for 2,516 SARS-CoV-2 (discovery cohort) and 1,126 (validation cohort) using nasopharyngeal samples collected between the time period of March 2020 to August 2022, as part of the genomic surveillance program. The sequences were classified according to the different lineages of SARS-CoV-2 prevailing in India at respective time points. Mutational diversity and abundance evaluation across the 12 lineages identified 58 Recombinant sequences as harboring the least number of mutations (n = 111), with 14 low-frequency unique mutations with major chunk of mutations coming from the BA.2. The spontaneously/dynamically increasing and decreasing trends of mutations highlight the loss of mutations in the Recombinants that were associated with the SARS-CoV-2 replication efficiency, infectivity, and disease severity, rendering them functionally with low infectivity and pathogenicity. Linkage disequilibrium (LD) analysis revealed that mutations comprising the LD blocks of BA.1, BA.2, and Recombinants were found as minor alleles or as low-frequency alleles in the LD blocks from the previous SARS-CoV-2 variant samples, especially Pre-VOC. Moreover, a dissipation in the size of LD blocks as well as LD decay along with a high negative regression coefficient (R squared) value was demonstrated in the Omicron and BA.1 and BA.2 lineages, which corroborated with the breakpoint analysis. Together, the findings help to understand the evolution and emergence of Recombinants after the Omicron lineages, for sustenance and adaptability, to maintain the epidemic spread of SARS-CoV-2 in the host population already high in immunity levels.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events.

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

Genomic epidemiology and evolutionary analysis during XBB.1.16-predominant periods of SARS-CoV-2 omicron variant in Bangkok, Thailand: December 2022–August 2023

The growing occurrence of novel recombinants, such as XBB.1.16, has emerged and become predominant, raising concerns about the impact of genomic recombination on the evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study investigated the molecular epidemiological trends and evolution of the Omicron XBB.1.16 epidemic in Bangkok between December 2022 and August 2023. Partial spike and complete genome sequencing of SARS-CoV-2 samples collected from collaborating hospitals were performed. The analysis of 491 partial spike sequences identified 15 distinct lineages, with XBB.1.16 dominating the lineages beginning in March 2023. Phylogenetic analysis revealed at least four distinct XBB.1.16 lineages, suggesting multiple independent introductions into Bangkok. The estimated emergence of XBB.1.16 occurred approximately in January 2022, with an evolutionary rate of 0.79 × 10–3 substitutions per site per year. Monitoring the genomic epidemiology and evolution of XBB.1.16 is vital for the early detection of new strains or emerging variants, which may guide vaccine design and the inclusion of new vaccine strains.