Rearrangement Events Research Articles

New mitochondrial gene order arrangements and evolutionary implications in the class Octocorallia.

The complete mitochondrial genomes of octocorals typically range from 18.5 kb to 20.5 kb in length and include 14 protein-coding genes (PCGs), two ribosomal RNA genes and one tRNA. To date, seven different gene orders (A-G) have been described, yet comprehensive investigations of the actual number of arrangements, as well as comparative analyses and evolutionary reconstructions of mitochondrial genome evolution within the whole class Octocorallia, have been often overlooked. Here, we considered the complete mitochondrial genomes available for octocorals and explored their structure and gene order variability. Our results updated the actual number of mitochondrial gene order arrangements so far known for octocorals from 7 to 14 and allowed us to explore and preliminarily discuss the role of some of the structural and functional factors in the mitogenomes. We performed comparative mitogenomic analyses on the existing and novel octocoral gene orders, considering different mitogenomic structural features such as genome size, GC percentage, AT and GC skewness. The mitochondrial gene order history mapped on a recently published nuclear loci phylogeny showed that the most common rearrangement events in octocorals are inversions, inverted transpositions and transpositions. Furthermore, gene order rearrangement events were restricted only to some regions of the tree. Overall, different rearrangement events arose independently and from the ancestral and most common gene order, instead of being derived from other rearranged orders. Finally, our data demonstrate how the study of mitochondrial gene orders can be used to explore the evolution of octocorals and in some cases can be used to assess the phylogenetic placement of certain taxa.

Applying rearrangement distances to enable plasmid epidemiology with pling.

Plasmids are a key vector of antibiotic resistance, but the current bioinformatics toolkit is not well suited to tracking them. The rapid structural changes seen in plasmid genomes present considerable challenges to evolutionary and epidemiological analysis. Typical approaches are either low resolution (replicon typing) or use shared k-mer content to define a genetic distance. However, this distance can both overestimate plasmid relatedness by ignoring rearrangements, and underestimate by over-penalizing gene gain/loss. Therefore a model is needed which captures the key components of how plasmid genomes evolve structurally - through gene/block gain or loss, and rearrangement. A secondary requirement is to prevent promiscuous transposable elements (TEs) leading to over-clustering of unrelated plasmids. We choose the 'Double Cut and Join Indel' (DCJ-Indel) model, in which plasmids are studied at a coarse level, as a sequence of signed integers (representing genes or aligned blocks), and the distance between two plasmids is the minimum number of rearrangement events or indels needed to transform one into the other. We show how this gives much more meaningful distances between plasmids. We introduce a software workflow pling (https://github.com/iqbal-lab-org/pling), which uses the DCJ-Indel model, to calculate distances between plasmids and then cluster them. In our approach, we combine containment distances and DCJ-Indel distances to build a TE-aware plasmid network. We demonstrate superior performance and interpretability to other plasmid clustering tools on the 'Russian Doll' dataset and a hospital transmission dataset.

Characterization of the Complete Mitochondrial Genome of Pleurogenoides japonicus (Digenea, Pleurogenidae): Comparison With the Members of Microphalloidea and Phylogenetic Implications.

Pleurogenoides japonicus (Trematoda: Microphalloidea) is an important parasite in wood frogs with high infection rates and significant ecological, economic, and societal importance. The scarcity of molecular data for these parasites severely limits population genetics and phylogenetic studies. In the present study, for the first time, we determined and described the entire mitochondrial (mt) genome of P. japonicus as the first representative of the family Pleurogenidae. The entire mt genome of P. japonicus was circular, with 15,043 bp (GenBank accession number OR900118), containing 36 genes, comprising 12 protein-coding genes (cox1-3, nad1-6, nad4L, cytb, and atp6), two ribosomal RNA genes, 22 transfer RNA genes, and two non-coding regions. There were 23 intergenic spacers, ranging from 2 to 162 bp, and only one 40 bp overlap between nad4L and nad4 genes in the P. japonicus mt genome. The nucleotide composition of P. japonicus mt genome exhibited a strong AT bias with a 63.75% A + T content, while the AT- and GC-skews were - 0.435 and 0.407, respectively. Comparative analysis demonstrated that the P. japonicus mt genome shared the most common characteristics with Microphalloidea trematodes, and the cox1 gene was the longest and most conserved gene in Microphalloidea trematodes. The gene arrangements of Xiphidiata trematodes were of the same order based on protein-coding genes and rRNA genes, except for tRNA. More than two gene arrangement types exist in Echinostomata and Xiphidiata, and the gene rearrangement events mainly occurred in "trnE-trnG" and "trnG-trnE". Phylogenetic analysis suggested that trematodes of the family Pleurogenidae clustered more with Prosthogonimidae than Eucotylidae. The mt genome data of P. japonicus provide an accurate genetic marker for further studies of Xiphidiata trematodes.

High-quality genome of a novel Thermosynechococcaceae species from Namibia and characterization of its protein expression patterns at elevated temperatures.

Thermophilic cyanobacteria thrive in extreme environments, making their thermoresistant enzymes valuable for industrial applications. Common habitats include hot springs, which act as evolutionary accelerators for speciation due to geographical isolation. The family Thermosynechococcaceae comprises thermophilic cyanobacteria known for their ability to thrive in high-temperature environments. These bacteria are notable for their photosynthetic capabilities, significantly contributing to primary production in extreme habitats. Members of Thermosynechococcaceae exhibit unique adaptations that allow them to perform photosynthesis efficiently at elevated temperatures, making them subjects of interest for studies on microbial ecology, evolution, and potential biotechnological applications. In this study, the genome of a thermophilic cyanobacterium, isolated from a hot spring near Okahandja in Namibia, was sequenced using a PacBio Sequel IIe long-read platform. Cultivations were performed at elevated temperatures of 40, 50, and 55°C, followed by proteome analyses based on the annotated genome. Phylogenetic investigations, informed by the 16S rRNA gene and aligned nucleotide identity (ANI), suggest that the novel cyanobacterium is a member of the family Thermosynechococcaceae. Furthermore, the new species was assigned to a separate branch, potentially representing a novel genus. Whole-genome alignments supported this finding, revealing few conserved regions and multiple genetic rearrangement events. Additionally, 129 proteins were identified as differentially expressed in a temperature-dependent manner. The results of this study broaden our understanding of cyanobacterial adaptation to extreme environments, providing a novel high-quality genome of Thermosynechococcaceae cyanobacterium sp. Okahandja and several promising candidate proteins for expression and characterization studies.

Combining optical genome mapping and RNA-seq for structural variants detection and interpretation in unsolved neurodevelopmental disorders

BackgroundStructural variations (SVs) are key genetic contributors to neurodevelopmental disorders (NDDs). Exome sequencing (ES), the current first-line tool for genetic testing of NDDs, falls short in SVs detection. This diagnostic gap is being actively addressed by new methods such as optical genome mapping (OGM).MethodsThis study evaluated the utility of combining OGM and RNA-seq in the detection and interpretation of SVs in ES-negative NDDs. OGM was performed in 43 patients with NDDs with inconclusive ES results. Candidate SVs were selected based on disease association and pathogenicity evaluation, and further validated or reconstructed by alternative methods, including long-read sequencing for a complex rearrangement event. RNA-Seq was performed on blood samples from patients with candidate SVs to facilitate interpretation of pathogenicity.ResultsOGM detected four candidate SVs, and RNA-seq confirmed the pathogenicity of three SVs in the patient cohort. This combined approach solved three cases—two cases with de novo SVs in genes associated with autosomal dominant NDDs, including a deletion encompassing the promoter and 5′UTR of MBD5 and an intragenic duplication of PAFAH1B1, and a third case possessing an intragenic duplication in trans with a pathogenic single-nucleotide variant of PLA2G6, associated with autosomal recessive NDDs. The expression alteration of the affected genes and the tandem positioning of two intragenic duplications were confirmed by RNA-seq. In the fourth case, OGM detected a complex rearrangement involving chromosomes 2 and 6, much more complex than the de novo t(2:6)(q13;q15) indicated by conventional cytogenetic analysis. Reconstruction showed that 17 segments of 6q15 spanning 9.3 Mb were disarranged and joined 2q11.2, with four breakpoints detected in the 5′ and 3′ non-coding region of the NDD-associated gene SYNCRIP. RNA-seq revealed largely preserved SYNCRIP expression, leaving the pathogenicity of this complex rearrangement event uncertain.ConclusionsSVs in ES-negative NDDs can be identified by OGM, which is particularly useful for SVs in non-coding regions not covered by ES. OGM helps to construct complex SVs and provides information on the location and orientation of duplications, which is crucial for pathogenicity interpretation. The integration of RNA-seq facilitates the interpretation of the functional consequences of SVs at the transcriptional level. These findings demonstrate the utility and feasibility of combining OGM and RNA-seq in ES-negative cases with NDDs.

Characterization of a multiresistance optrA- and lsa(E)-harbouring unconventional circularizable structure in Streptococcus suis.

To identify novel genetic elements facilitating the horizontal transfer of the oxazolidinone/phenicol resistance gene optrA and the pleuromutilin-lincosamide-streptogramin A resistance gene lsa(E) in Streptococcus suis. The complete genomes of S. suis HB18 and two transconjugants were obtained using both the Illumina and Nanopore platforms. MICs were determined by broth microdilution. Inverse PCR was performed to identify circular forms of the novel unconventional circularizable structure (UCS), genomic island (GI) and integrative and conjugative element (ICE). Conjugation experiments assessed the transferability of optrA and lsa(E) genes in S. suis. S. suis HB18 carried a multiresistance gene cluster optrA-lsa(E)-lnu(B)-aphA-aadE-spw. This gene cluster, flanked by intact and truncated erm(B) in the same orientation, resided on a novel ICESsuHB18. Inverse PCR revealed the existence of a novel UCS, named UCS-optrA + lsa(E), which could excise the gene cluster optrA-lsa(E)-lnu(B)-aphA-aadE-spw and one copy of erm(B) from ICESsuHB18. Two transconjugants with different characteristics were obtained. In transconjugant T-JH-GI, UCS-optrA + lsa(E) excised from ICESsuHB18 inserted into the erm(B)-positive GI, designated GISsuHB18, generating the novel GISsuHB18-1. Meanwhile, in T-JH-ICE, genetic rearrangement events occurred in ICESsuHB18 and GISsuHB18, forming the novel ICESsuHB18-1. This is the first report demonstrating the functionally active UCS-optrA + lsa(E) excising from ICESsuHB18 and inserting into the erm(B)-positive GISsuHB18 during the conjugation process. The location of optrA and lsa(E) on a multiresistance UCS enhances its persistence and dissemination.

E-cadherin tunes tissue mechanical behavior before and during morphogenetic tissue flows

Adhesion between epithelial cells enables the remarkable mechanical behavior of epithelial tissues during morphogenesis. However, it remains unclear how cell-cell adhesion influences mechanics in both static and dynamically flowing confluent epithelial tissues. Here, we systematically modulate E-cadherin-mediated adhesion in the Drosophila embryo and study the effects on the mechanical behavior of the germband epithelium before and during dramatic tissue remodeling and flow associated with body axis elongation.Before axis elongation, we find that increasing E-cadherin levels produces tissue comprising more elongated cells and predicted to be more fluid-like, providing reduced resistance to tissue flow. During axis elongation, we find that the dominant effect of E-cadherin is tuning the speed at which cells proceed through rearrangement events. Before and during axis elongation, E-cadherin levels influence patterns of actomyosin-dependent forces, supporting the notion that E-cadherin tunes tissue mechanics in part through effects on actomyosin. Notably, the effects of ∼4-fold changes in E-cadherin levels on overall tissue structure and flow are relatively weak, suggesting that the system is tolerant to changes in absolute E-cadherin levels over this range where an intact tissue is formed. Taken together, these findings reveal dual-and sometimes opposing-roles for E-cadherin-mediated adhesion in controlling tissue structure and dynamics invivo, which result in unexpected relationships between adhesion and flow in confluent tissues.

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.

Characterization and novel gene rearrangements in the mitochondrial genome of Chelonus munakatae (Hymenoptera: Braconidae): Insights into evolutionary dynamics and genomic variability

Chelonus munakatae Matsumura (Hymenoptera: Braconidae) is an important parasitoid of several lepidopteran species. The mitochondrial genome of this species is a crucial source for molecular identification and genetic studies. Here, we used next-generation sequencing to obtain the C. munakatae mitochondrial genome, which was found to be 14,868 bp in length and contained 13 protein-coding genes, 22 tRNAs and two rRNAs. The nucleotide composition exhibits a significant bias towards A + T content, with a value of 86.8 %. A large number of PCG and tRNA gene rearrangement events were detected (I-V-M-C-W-Q and A-N-S1-E-R-F), which represent a novel rearrangement type in Hymenoptera. Furthermore, synteny analysis based on the whole mitochondrial genomes showed that a significant rearrangement of homologous blocks had occurred within this Braconidae species. This study not only delineates the mitochondrial characteristics of C. munakatae but also significantly contributes to the repository of mitochondrial genetic information for the Braconidae family. Moreover, it highlights the occurrence of gene rearrangement within Hymenoptera species, offering valuable insights into the evolutionary dynamics and genomic variability of these insects.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome.

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution.

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

Rearrangement Distance Problems: An updated survey

One of the challenges in the Comparative Genomics field is to infer how close two organisms are based on the similarities and differences between their genetic materials. Recent advances in DNA sequencing have made complete genomes increasingly available. That said, several new algorithms trying to infer the distance between two organisms based on genome rearrangements have been proposed in the literature. However, given the diversity of approaches, the diversity of genome rearrangement events, or even how each work models the genomes and what assumptions are made by each of them, finding the ideal algorithm for each situation or simply knowing the range of applicable approaches can be challenging. In this work, we review these approaches having the algorithmic and combinatorial advances since 2010 as our main focus. This survey aims to organize the recently published papers using a concise notation and to indicate the gaps filled by each of them in the literature. This makes it easier to understand what still needs to be done and what has room for enhancement.

Abstract 4337: Illuminating the complex rearrangement landscape of sarcoma

Abstract Structural rearrangements and in particular, gene fusions, are key mutational processes in bone and soft tissue tumors, used for both disease classification and as biomarkers. However, the mutational processes and rearrangement processes underlying many of these events remain poorly characterized. Recent data indicate that sarcomas show particularly high frequencies of complex rearrangement events, including patterns which do not fit those of known mutational mechanisms. As the largest whole genome sequencing (WGS) cohort of sarcomas to date, the Genomics England (GE) 100,000 Genomes project represents a unique dataset in which to profile these events. Structural variants (SVs) were identified in whole genome sequencing data for 978 samples from the Genomics England 100,000 Genomes Project using an optimized approach based on 5 specialized SV callers. SVs identified by at least 2 callers were filtered to retain only those which passed the quality filter, and to remove variants present in either matched germline samples or >1% of a panel of normal variants. Chromothriptic events, extrachromosomal DNA (ecDNA) and genes enriched for SVs were identified using established algorithms. The prevalence of SVs varies across sarcoma, in keeping with previous reports, with particularly high rates observed in myxofibrosarcoma and liposarcoma. As described previously, liposarcomas showed the highest rates of chromothripsis. We have examined some tumor types at higher granularity than previously available, demonstrating directly that the rates of complex rearrangement events vary across subtypes. For example, chromothriptic events are identified in 100% of well-differentiated liposarcomas, but less than 10% of myxoid liposarcomas. Similarly, the presence of ecDNA varies by tumor type, with the novel observation of particularly high rates in angiosarcoma. Using these data, we have catalogued the rearrangement patterns generating canonical gene fusions in fusion-driven soft tissue tumor types. Despite more recent advances in histological classification, survival for patients with sarcoma has remained largely unchanged for 40 years. Characterization of the mutational processes underlying these rearrangements will shed light on the pathogenesis of these tumors. Citation Format: Sara Waise, Tom Lesluyes, Jonas Demeulemeester, Maxime Tarabichi, Adrienne Flanagan, Nischalan Pillay, Peter Van Loo. Illuminating the complex rearrangement landscape of sarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 4337.

Description and mitochondrial genome sequencing of a new species of inquiline gall wasp, Synergus nanlingensis (Hymenoptera, Cynipidae, Synergini), from China

A new species of inquiline gall wasp, Synergus nanlingensis Wang & Zeng, sp. nov., which was reared from galls on Castanopsis eyrei Tutch (Fagaceae) collected in Guangdong Province, China, is described and illustrated herein along with its mitochondrial genome. The mitogenome of S. nanlingensis is 16,604 base pairs in length and comprises 37 genes, which is typical of mitogenomes. One large control region was detected in the S. nanlingensis mitogenome, which differed from that reported for other Cynipidae species. Similar to other Cynipidae species, S. nanlingensis has the same four common gene rearrangement events; however, it shows some differences, as follows: trnS1 is downstream of Cytb; trnS2 is upstream of nad1; and trnC is downstream of rrnS. Phylogenetic analysis using COI, CytB, and 28S-D2 sequences confirmed that S. nanlingensis is a distinct species belonging to the genus Synergus Hartig.

Comparative Mitogenomics Analysis Revealed Evolutionary Divergence among Neopestalotiopsis Species Complex (Fungi: Xylariales).

The genus Neopestalotiopsis consists of obligate parasites that cause ring spot, scab, and leaf blight diseases in higher plant species. We assembled the three complete mitogenomes for the guava fruit ring spot pathogen, Neopestalotiopsis cubana. The mitogenomes are circular, with sizes of 38,666 bp, 33,846 bp, and 32,593 bp. The comparative analyses with Pestalotiopsis fici showed that N. cubana differs greatly from it in the length of the mitogenomes and the number of introns. Moreover, they showed significant differences in the gene content and tRNAs. The two genera showed little difference in gene skewness and codon preference for core protein-coding genes (PCGs). We compared gene sequencing in the mitogenomes of the order Xylariales and found large-scale gene rearrangement events, such as gene translocations and the duplication of tRNAs. N. cubana shows a unique evolutionary position in the phylum Ascomycota constructed in phylogenetic analyses. We also found a more concentrated distribution of evolutionary pressures on the PCGs of Neopestalotiopsis in the phylum Ascomycota and that they are under little selective pressure compared to other species and are subjected to purifying selection. This study explores the evolutionary dynamics of the mitogenomes of Neopestalotiopsis and provides important support for genetic and taxonomic studies.

Allergenicity and Conformational Diversity of Allergens

Allergens are substances that cause abnormal immune responses and can originate from various sources. IgE-mediated allergies are one of the most common and severe types of allergies, affecting more than 20% of the population in Western countries. Allergens can be subdivided into a limited number of families based on their structure, but this does not necessarily indicate the origin or the route of administration of the allergen, nor is the molecular basis of allergenicity clearly understood. This review examines how understanding the allergenicity of proteins involves their structural characterization and elucidates the study of conformational diversity by nuclear magnetic resonance spectroscopy. This article also discusses allergen cross-reactivity and the mechanisms by which IgE antibodies recognize and bind to allergens based on their conformational and linear epitopes. In addition, we outline how the pH, the proteolytic susceptibility and the endosomal degradation affect the outcome of allergic reactions, and how this is correlated with conformational changes and secondary structure rearrangement events. We want to emphasize the importance of considering structural diversity and dynamics, proteolytic susceptibility and pH-dependent factors to fully comprehend allergenicity.

Partition Based Algorithms for Rearrangement Distances with Flexible Intergenic Regions.

Genome Rearrangement distance problems are used in Computational Biology to estimate the evolutionary distance between genomes. These problems consist of minimizing the number of rearrangement events necessary to transform one genome into another. Two commonly used rearrangement events are reversal and transposition. The first studied problems ignored nucleotides outside genes (called intergenic regions), or assumed that genomes have a single copy of each gene. Recent works made advancements in more general problems considering the number of nucleotides in intergenic regions, and replicated genes. Nevertheless, genomes tend to have wildly different quantities of nucleotides on their intergenic regions, which poses a problem when comparing these regions exactly. To overcome this limitation, our work considers some flexibility when matching intergenic regions that do not have the same number of nucleotides. We propose new problems seeking the minimum number of reversals, or reversals and transpositions, necessary to transform one genome into another, while considering flexible intergenic region information. We show approximations for these problems by exploring their relationship with the Signed Minimum Common Flexible Intergenic String Partition problem. We also present different heuristics for the partition problem, and conduct experimental tests on simulated genomes to assess the performance of our algorithms.

Impact of HBV Integration on Hepatocellular Carcinoma After Long-Term Antiviral Therapy.

Few studies have reported the integrated characteristics of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) after long-term antiviral therapy. This study aimed to investigate the HBV integration features in HBV-HCC patients who had undergone long-term antiviral therapy, evaluate their impact on clinical indicators, and analyze the potential mechanisms involved. We utilized genome-wide association study (GWAS) to analyze liver cancer tissues and detect the presence of HBV integration. Seventeen patients with HBV integration were included in the integration (Int) group, while the remaining five patients were included in the non-integration (N-int) group. Clinical indicators were regularly monitored and compared between the two groups. The characteristics of HBV integration patterns were analyzed, and differences between the groups were explored at the chromosome and genomic levels. After long-term antiviral therapy, although the frequency of HBV integration in HBV-HCC was reduced, residual HBV integration still accelerated the development of HCC. It affected the diagnosis, treatment, and prognosis of patients. HBV integration events led to changes in chromosome structure, which were closely related to HCC. Novel fusion genes were detected at a high frequency and had the potential to be specific detection sites for HBV-HCC. HBV integration events are synergistically involved in the human genome and HBV, which can lead to chromosome structural instability, gene rearrangement events closely related to HCC production, and the formation of new specific fusion genes.

Multiple independent origins of duplicated mitochondrial control regions indicate an apomorphy in the Thysanoptera (Insecta).

The mitochondrial genome (mitogenome) of thrips is characterized by the presence of control region (CR) duplication. However, the evolution pattern of duplicated CRs in thrips is still unclear. In this study, the multiple independent origins of duplicated CR indicated that the CR duplication was not an ancestral state for Thysanoptera. The macroevolutionary pattern suggested that the earliest CR duplication event occurred in the middle Cretaceous (94.85 Ma) coincided with rearrangement events forming the ancestors of Aeolothripidae, but much later than that forming the ancestors of the suborder Terebrantia. The mitogenome with duplicated CRs showed a higher rate of gene rearrangement. The sequence similarity of the CR copies and divergence time were negatively correlated, indicating age-related deterioration of mitochondrial function. No significant differences were found in the mitochondrial DNA, the P123 and P4FD between the single and multiple-CR charactered mitogenomes, which suggested that the duplicated CRs may not affect the replication process in thrip mitogenome. The mitogenomes with duplicated CRs (mean: 0.0088 subs/s/my) show a significantly increased evolutionary rate than that with a single one (mean: 0.0058 subs/s/my). However, it seems that this higher evolutionary rate did not have adaptive mechanisms in Terebrantia. We speculated that the duplicated CRs may cause a more intense production of energy by mitochondria, and an accelerated mutation and substitution rate is expected in such mitogenomes. Our study provided new insights into the presence of CR duplications and their evolution in the mitogenomes of thrips.

Circulation, viral diversity and genomic rearrangement in mpox virus in the Netherlands during the 2022 outbreak and beyond.

Mpox is an emerging zoonotic disease which has now spread to over 113 countries as of August 2023, with over 89,500 confirmed human cases. The Netherlands had one of the highest incidence rates in Europe during the peak of the outbreak. In this study, we generated 158 near-complete mpox virus (MPXV) genomes (12.4% of nationwide cases) that were collected throughout the Netherlands from the start of the outbreak in May 2022 to August 2023 to track viral evolution and investigate outbreak dynamics. We detected 14 different viral lineages, suggesting multiple introductions followed by rapid initial spread within the country. The estimated evolutionary rate was relatively high compared to previously described in orthopoxvirus literature, with an estimated 11.58 mutations per year. Genomic rearrangement events occurred at a rate of 0.63% and featured a large deletion event. In addition, based on phylogenetics, we identified multiple potential transmission clusters which could be supported by direct source- and contact tracing data. This led to the identification of at least two main transmission locations at the beginning of the outbreak. We conclude that whole genome sequencing of MPXV is essential to enhance our understanding of outbreak dynamics and evolution of a relatively understudied and emerging zoonotic pathogen.