Genomic Elements Research Articles

EPCO-19. MULTIOMIC PROFILING OF EPIGENETIC HETEROGENEITY IN GLIOBLASTOMA AT SINGLE-CELL RESOLUTION USING MULTIOMIC SINGLE-CELL CUT&TAG

Abstract INTRODUCTION Recent evidence suggests epigenetic heterogeneity may be the primary driver of cell phenotype and treatment resistance in glioblastoma (GBM). Epigenetic heterogeneity is defined by alterations of chromatin histone modifications with activating marks such as H3K27 acetylation and repressive marks such as H3K27 trimethylation. Combinations of these histone marks govern the activity of functional genomic elements such as enhancers or promoters. Single-cell (“sc-”) technologies have only recently been adapted to the study of epigenetic cell states and significant technical challenges have limited their widespread adoption, especially in complex tissues. Here, we report the development of an optimized method for performing tissue processing and simultaneous single-cell epigenetic and transcriptomic analyses in flash frozen GBM samples. METHODS Building on our prior work, we developed a new protocol, multiomic single-cell cleavage-under-targets and tagmentation, or “mscCUT&TAG”. We optimized nuclei isolation, cleanup and permeabilization from flash frozen GBM and reaction conditions for the scCUT&TAG portion of the assay. Critically, these changes maintained high-quality, intact nuclei in suspension while limiting well-described technical issues such as nuclei clumping and sample loss. Furthermore, we incorporated the ability to multiplex both patient samples and histone marks. RESULTS We applied 10X scRNA-seq, 10X Multiome (RNA+ATAC-seq), and mscCUT&TAG for five histone marks to the same GBM sample, enabling an integrated evaluation of the tumor’s transcriptome, chromatin accessibility and epigenetic regulation at single-cell resolution. mscCUT&TAG enabled the first assessment of single-cell genome segmentation in GBM. We integrated across modalities to identify specific tumor cell populations, cell-type specific coordinated changes in epigenetics and transcriptome, and define functional regulatory elements. CONCLUSIONS We developed a novel multiomic protocol, mscCUT&TAG, that allows for a comprehensive analysis of epigenetic heterogeneity in GBM. We anticipate this protocol will enhance our understanding of fundamental biology and potential therapeutic avenues in GBM.

From the genome's perspective: Bearing somatic retrotransposition to leverage the regulatory potential of L1 RNAs

AbstractTransposable elements (TEs) are mobile genomic elements constituting a big fraction of eukaryotic genomes. They ignite an evolutionary arms race with host genomes, which in turn evolve strategies to restrict their activity. Despite being tightly repressed, TEs display precisely regulated expression patterns during specific stages of mammalian development, suggesting potential benefits for the host. Among TEs, the long interspersed nuclear element (LINE‐1 or L1) has been found to be active in neurons. This activity prompted extensive research into its possible role in cognition. So far, no specific cause‐effect relationship between L1 retrotransposition and brain functions has been conclusively identified. Nevertheless, accumulating evidence suggests that interactions between L1 RNAs and RNA/DNA binding proteins encode specific messages that cells utilize to activate or repress entire transcriptional programs. We summarize recent findings highlighting the activity of L1 RNAs at the non‐coding level during early embryonic and brain development. We propose a hypothesis suggesting a mutualistic relationship between L1 mRNAs and the host cell. In this scenario, cells tolerate a certain rate of retrotransposition to leverage the regulatory effects of L1s as non‐coding RNAs on potentiating their mitotic potential. In turn, L1s benefit from the cell's proliferative state to increase their chance to mobilize.

Genomic landscape of adult testicular germ cell tumours in the 100,000 Genomes Project

Testicular germ cell tumours (TGCT), which comprise seminoma and non-seminoma subtypes, are the most common cancers in young men. In this study, we present a comprehensive whole genome sequencing analysis of adult TGCTs. Leveraging samples from participants recruited via the UK National Health Service and data from the Genomics England 100,000 Genomes Project, our results provide an extended description of genomic elements underlying TGCT pathogenesis. This catalogue offers a comprehensive, high-resolution map of copy number alterations, structural variation, and key global genome features, including mutational signatures and analysis of extrachromosomal DNA amplification. This study establishes correlations between genomic alterations and histological diversification, revealing divergent evolutionary trajectories among TGCT subtypes. By reconstructing the chronological order of driver events, we identify a subgroup of adult TGCTs undergoing relatively late whole genome duplication. Additionally, we present evidence that human leukocyte antigen loss is a more prevalent mechanism of immune disruption in seminomas. Collectively, our findings provide valuable insights into the developmental and immune modulatory processes implicated in TGCT pathogenesis and progression.

Endogenous nege-like viral elements in arthropod genomes reveal virus-host coevolution and ancient history of two plant virus families.

Although negevirus is phylogenetically related to plant virus, the evolutionary history of negevirus-host and its relationship with plant virus remain largely unknown. In this study, we used endogenous nege-like viral elements (ENVEs) as the molecular fossil records to investigate the history of nege-like viral infection in arthropod hosts and the evolution of two related plant virus families (Virgaviridae and Kitaviridae). Our results showed the infection of nege-like viruses for over 60 My during the arthropod evolution. ENVEs highly homologous to viral sequences in Virgaviridae and Kitaviridae were present in a wide range of arthropod genomes but were absent in plant genomes, indicating that plant viruses in these two families possibly evolved from the nege-like virus ancestor through cross-species horizontal virus transmission. Our findings provide a new perspective on the virus-host coevolution and the origins of plant viruses.

Deep mining reveals the diversity of endogenous viral elements in vertebrate genomes

Integration of viruses into host genomes can give rise to endogenous viral elements (EVEs), which provide insights into viral diversity, host range and evolution. A systematic search for EVEs is becoming computationally challenging given the available genomic data. We used a cloud-computing approach to perform a comprehensive search for EVEs in the kingdoms Shotokuvirae and Orthornavirae across vertebrates. We identified 2,040 EVEs in 295 vertebrate genomes and provide evidence for EVEs belonging to the families Chuviridae, Paramyxoviridae, Nairoviridae and Benyviridae. We also find an EVE from the Hepacivirus genus of flaviviruses with orthology across murine rodents. In addition, our analyses revealed that reptarenaviruses and filoviruses probably acquired their glycoprotein ectodomains three times independently from retroviral elements. Taken together, these findings encourage the addition of 4 virus families and the Hepacivirus genus to the growing virus fossil record of vertebrates, providing key insights into their natural history and evolution.

Human Satellite 3 DNA encodes megabase-scale transcription factor binding platforms.

Eukaryotic genomes are frequently littered with large arrays of tandem repeats, called satellite DNA, which underlie the constitutive heterochromatin often found around centromeric regions. While some satellite DNA types have well-established roles in centromere biology, other abundant satellite DNAs have poorly characterized functions. For example, Human Satellite 3 (HSat3), which makes up roughly 2% of the human genome, forms enormous arrays up to tens of megabases, but these arrays play no known roles in centromere function and were almost fully excluded from genome assemblies until recently. As a result, these massive genomic regions have remained relatively understudied, and the potential functional roles of HSat3 have remained largely unknown. To address this, we performed a systematic screen for novel HSat3 binding factors. Our work revealed HSat3 arrays contain high densities of transcription factor (TF) motifs that are bound by factors related to multiple, highly conserved signaling pathways. Unexpectedly, the most enriched TFs in HSat3 belong to the Hippo pathway transcription effector family TEAD. We found that TEAD recruits the co-activator YAP to HSat3 regions in a cell-state specific manner. Using RNA polymerase-I reporter assays, targeted repression of HSat3, inducible degradation of YAP, and super-resolution microscopy, we show that HSat3 arrays can localize YAP/TEAD inside the nucleolus, where YAP regulates RNA Polymerase-I activity. Beyond revealing a direct relationship between the Hippo pathway and ribosomal DNA regulation, this work demonstrates that satellite DNA can encode multiple transcription factor binding motifs, defining a new role for these enormous genomic elements.

Characterization of non-standard viral genomes during arenavirus infections identifies prominent S RNA intergenic region deletions.

Arenaviruses, a family of negative-sense RNA viruses spread by rodents, are a leading cause of severe hemorrhagic fever in humans. Due to a paucity of antivirals and vaccines for arenaviruses, there is a need to identify new mechanisms for interfering with arenavirus replication. In several negative-sense RNA viruses, natural viral interference results from the production of non-standard viral genomes (nsVGs) that activate the innate immune system and/or compete for essential viral products. Although it is well established that arenaviruses produce strong interfering activities, it is unknown if they produce interfering nsVGs. Here, we show that arenaviruses produce deletions within the intergenic region of their small (S) RNA genome, and these deletions inhibit viral glycoprotein production during minigenome replication. S RNA deletions are more abundant when arenaviruses are grown in high-interfering conditions and are associated with reduced viral replication. Overall, we found that arenaviruses produce internal deletions within the S RNA intergenic region that are capable of decreasing glycoprotein production. These natural arenavirus interfering molecules provide a new target for the generation of therapeutics against arenaviruses.IMPORTANCEArenaviruses are hemorrhagic fever-causing pathogens that infect millions of people a year. There are currently no approved antivirals that target arenaviruses, and understanding natural mechanisms that inhibit arenavirus replication is crucial for the development of effective therapeutics. Here, we identified multiple deletions within arenavirus genomes that remove major replicative elements of the viral genomes. We show that deletions that remove the intergenic region of the viral genome can prevent viral protein production. These deletions were found in all arenaviruses tested in this study representing a mechanism that could be harnessed for the development of antivirals that broadly target the arenavirus family.

Integrating Biomarker Testing into EHR Systems: Implications for the Laboratory Based on Findings from a Multistakeholder Summit

Abstract Introduction/Objective Guideline-driven tumor biomarker testing remains standard of care for precision oncology due to the increasing number of treatment decisions that can be made based on tumor-specific alterations. However, most electronic health record (EHR) systems lack bidirectional interfaces with labs that enable computerized order entry and tumor biomarker testing results. Without these integrations, cancer centers may experience delays in test ordering and resulting, manual data entry errors, and operational inefficiencies that impact patient care and survival outcomes. Methods/Case Report To explore how cancer centers, hospitals, oncology practices, and labs integrate EHRs to streamline cancer biomarker testing processes, the Association of Cancer Care Centers (ACCC), in partnership with LUNGevity, convened a multistakeholder summit on Oct. 4, 2023. The summit brought together clinicians, pathologists, genetic counselors, and administrators to discuss practical solutions to issues surrounding interoperability, clinical workflows, and access to clinical data. Results (if a Case Study enter NA) Key findings from the summit include: -EHR-integrated orders and results save time, reduce errors, and streamline clinical workflow. These electronic interfaces enable clear communication between multidisciplinary teams system-wide and ensure that required patient information is readily available at the time of therapy planning. -Standardized workflows that define test ordering, tracking, and resulting protocols for in-house and send-out testing offer important downstream benefits. Staff should clearly understand how specimens are handled for in-house and send-out testing. -Integrations should be planned and built based on the EHR/laboratory information systems (LIS) in place can incorporate solutions to integrate with reference labs using HL7 standards and/or add-on modules to enable discrete genomic data elements. Conclusion Clinicians, pathologists, and administrators should prioritize EHR/LIS integrations to enable more effective cancer biomarker test ordering and resulting. These pre/post-analytic processes have significant implications for patient care. Findings from this summit can enable cancer programs and labs to evaluate current processes, work collaboratively to build integrations, and streamline biomarker testing workflows to enhance the quality and efficiency of delivering precision oncology.

DeOri 10.0: An Updated Database of Experimentally Identified Eukaryotic Replication Origins.

DNA replication is a complex and crucial biological process in eukaryotes. To facilitate the study of eukaryotic replication events, we present a database of eukaryotic DNA replication origins (DeOri), which collects scattered data and integrates extensive sequencing data on eukaryotic DNA replication origins. With continuous updates of DeOri, the number of datasets in the new release increased from 10 to 151 and the number of sequences increased from 16,145 to 9,742,396. Besides nucleotide sequences and bed files, corresponding annotation files, such as coding sequences (CDS), mRNA, and other biological elements within replication origins, are also provided. The experimental techniques used for each dataset, as well as other statistical data, are also presented on web page. Differences in experimental methods, cell lines, and sequencing technologies have resulted in distinct replication origins, making it challenging to differentiate between cell-specific and non-specific replication. We combined multiple replication origins at the species level, scored them, and screened them. The screened regions were considered as species-conservative origins. They are integrated and presented as reference replication origins (rORIs), including Homo sapiens, Gallus gallus, Mus musculus, Drosophila melanogaster, and Caenorhabditis elegans. Additionally, we analyzed the distribution of relevant genomic elements associated with replication origins at the genome level, such as CpG island (CGI), transcription start site (TSS), and G-quadruplex (G4). These analysis results allow users to select the desired data based on it. DeOri is available at http://tubic.tju.edu.cn/deori/.

Template switching during DNA replication is a prevalent source of adaptive gene amplification.

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.

Machine learning classification of archaea and bacteria identifies novel predictive genomic features

BackgroundArchaea and Bacteria are distinct domains of life that are adapted to a variety of ecological niches. Several genome-based methods have been developed for their accurate classification, yet many aspects of the specific genomic features that determine these differences are not fully understood. In this study, we used publicly available whole-genome sequences from bacteria (N=2546\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=2546$$\\end{document}) and archaea (N=109\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$N=109$$\\end{document}). From these, a set of genomic features (nucleotide frequencies and proportions, coding sequences (CDS), non-coding, ribosomal and transfer RNA genes (ncRNA, rRNA, tRNA), Chargaff’s, topological entropy and Shannon’s entropy scores) was extracted and used as input data to develop machine learning models for the classification of archaea and bacteria.ResultsThe classification accuracy ranged from 0.993 (Random Forest) to 0.998 (Neural Networks). Over the four models, only 11 examples were misclassified, especially those belonging to the minority class (Archaea). From variable importance, tRNA topological and Shannon’s entropy, nucleotide frequencies in tRNA, rRNA and ncRNA, CDS, tRNA and rRNA Chargaff’s scores have emerged as the top discriminating factors. In particular, tRNA entropy (both topological and Shannon’s) was the most important genomic feature for classification, pointing at the complex interactions between the genetic code, tRNAs and the translational machinery.ConclusionstRNA, rRNA and ncRNA genes emerged as the key genomic elements that underpin the classification of archaea and bacteria. In particular, higher nucleotide diversity was found in tRNA from bacteria compared to archaea. The analysis of the few classification errors reflects the complex phylogenetic relationships between bacteria, archaea and eukaryotes.

ScNanoSeq-CUT&Tag: a single-cell long-read CUT&Tag sequencing method for efficient chromatin modification profiling within individual cells.

Chromatin modifications are fundamental epigenetic marks that determine genome functions, but it remains challenging to profile those of repetitive elements and complex genomic regions. Here, we develop scNanoSeq-CUT&Tag, a streamlined method, by adapting modified cleavage under targets and tagmentation (CUT&Tag) to the nanopore sequencing platform for genome-wide chromatin modification profiling within individual cells. We show that scNanoSeq-CUT&Tag can accurately profile histone marks and transcription factor occupancy patterns at single-cell resolution as well as distinguish different cell types. scNanoSeq-CUT&Tag efficiently maps the allele-specific chromatin modifications and allows analysis of their neighboring region co-occupancy patterns within individual cells. Moreover, scNanoSeq-CUT&Tag can accurately detect chromatin modifications for individual copies of repetitive elements in both human and mouse genomes. Overall, we prove that scNanoSeq-CUT&Tag is a valuable single-cell tool for efficiently profiling histone marks and transcription factor occupancies, especially for previously poorly studied complex genomic regions and blacklist genomic regions.

A global genome dataset for Salmonella Gallinarum recovered between 1920 and 2024

Salmonella enterica serovar Gallinarum (S. Gallinarum) is an avian-specific pathogen responsible for fowl typhoid, a severe systemic disease with high mortality in chickens. This disease poses a substantial burden to the poultry industry, particularly in developing countries like China. However, comprehensive genome datasets on S. Gallinarum are lacking. Here, we present the most extensive S. Gallinarum genome dataset, comprising 574 well-collated samples. This dataset consists of 366 genomes sequenced in our laboratory and 208 publicly available genomes, collected from various continents over the past century. Using in silico prediction, we categorized S. Gallinarum into three distinct biovars. Regarding antimicrobial resistance, 238 strains (41.5%) carried antimicrobial resistance genes (ARGs) with a total of 635 records, while 232 strains (40.4%) exhibited multi-drug resistance. Mobile genomic elements (MGEs) serve as critical drivers for ARGs. Our dataset includes 5,636 MGEs records, with most MGEs belonging to prophages and plasmids. This dataset expands our understanding of the genomic characteristics of S. Gallinarum, providing valuable resources for future genomic studies to improve disease management.

Permeable TAD boundaries and their impact on genome-associated functions.

TAD boundaries are genomic elements that separate biological processes in neighboring domains by blocking DNA loops that are formed through Cohesin-mediated loop extrusion. Most TAD boundaries consist of arrays of binding sites for the CTCF protein, whose interaction with the Cohesin complex blocks loop extrusion. TAD boundaries are not fully impermeable though and allow a limited amount of inter-TAD loop formation. Based on the reanalysis of Nano-C data, a multicontact Chromosome Conformation Capture assay, we propose a model whereby clustered CTCF binding sites promote the successive stalling of Cohesin and subsequent dissociation from the chromatin. A fraction of Cohesin nonetheless achieves boundary read-through. Due to a constant rate of Cohesin dissociation elsewhere in the genome, the maximum length of inter-TAD loops is restricted though. We speculate that the DNA-encoded organization of stalling sites regulates TAD boundary permeability and discuss implications for enhancer-promoter loop formation and other genomic processes.

Evolution of cardiac genomic elements in humans and non-human primates.

Evolution of cardiac genomic elements in humans and non-human primates.

Effects of multisuperovulation on the transcription and genomic methylation of oocytes and offspring

BackgroundControlled ovarian stimulation is a common skill of assisted reproductive technologies (ARTs). In the clinic, some females would undergo more than one controlled ovarian stimulation cycle. However, few studies have focused on the influence of multi-superovulation on oocytes and offspring.ResultsHere, we found that multi-superovulation disrupted the transcriptome of oocytes and that the differentially expressed genes (DEGs) were associated mainly with metabolism and fertilization. The disruption of mRNA degradation via poly (A) size and metabolism might be a reason for the reduced oocyte maturation rate induced by repeated superovulation. Multi-superovulation results in hypo-genomic methylation in oocytes. However, there was an increase in the methylation level of CGIs. The DMRs are not randomly distributed in genome elements. Genes with differentially methylated regions (DMRs) in promoters are enriched in metabolic pathways. With increasing of superovulation cycles, the glucose and insulin tolerance of offspring is also disturbed.ConclusionsThese results suggest that multi-superovulation has adverse effects on oocyte quality and offspring health.

Genomic and Socioeconomic Determinants of Racial Disparities in Breast Cancer Survival: Insights from the All of Us Program.

Background: Breast cancer outcomes are worse among Black women in the U.S. compared to White women. While extensive research has focused on risk factors contributing to breast cancer; the role of genomic elements in health disparities between these racial groups remains unclear. This study aims to identify genomic variants and socioeconomic status (SES) determinants influencing racial disparities in breast cancer survival through multiple mediation analyses. Methods: Our investigation is based on the NIH-supported All of Us (AoU) program and analyzes 7452 female participants with malignant tumors of breast, including 5073 with genomic data. A log-rank test reveals significant racial differences in overall survival time between Black and White participants (p-value = 0.04). Multiple mediation analysis examines the effects of 9481 genetic variables across 23 chromosomes in explaining the racial disparity in survival, adjusting for SES variables. Results: 15 gene mutations, in addition to age, general health, and general quality of life, have significant effects (p-values < 0.001) in explaining the observed racial disparity. Mutations in TMEM132B, NARFL, SALL1, PAD12, RIPK1, ASB14, DCX, GNB1L, ARHGAP32, AL135787.1, WBP11, SLC16A12AS1, AP000345.1, IKBKB, and SUPT20H have significantly different distributions between Black and White participants. The disparity is completely explained by the included variables as the direct effect is insignificant (p-value = 0.73). Conclusions: The combined impact of SES determinants and genetic mutations can explain the observed differences in breast cancer survival among Black and White participants. Future studies will explore pathways and design in vivo and in vitro experiments to validate the functions of these genes.

A Cryptic Prophage Transcription Factor Drives Phenotypic Changes via Host Gene Regulation.

Cryptic prophages (CPs) are elements of bacterial genomes acquired from bacteriophage that infect the host cell and ultimately become stably integrated within the host genome. While some proteins encoded by CPs can modulate host phenotypes, the potential for Transcription Factors (TFs) encoded by CPs to impact host physiology by regulating host genes has not been thoroughly investigated. In this work, we report hundreds of host genes regulated by DicC, a DNA-binding TF encoded in the Qin prophage of Esherichia coli. We identified host-encoded regulatory targets of DicC that could be linked to known phenotypes of its induction. We also demonstrate that a DicC-induced growth defect is largely independent of other Qin prophage genes. Our data suggest a greater role for cryptic prophage TFs in controlling bacterial host gene expression than previously appreciated.

Maf-family bZIP transcription factor NRL interacts with RNA-binding proteins and R-loops in retinal photoreceptors.

RNA-binding proteins (RBPs) perform diverse functions including the regulation of chromatin dynamics and the coupling of transcription with RNA processing. However, our understanding of their actions in mammalian neurons remains limited. Using affinity purification, yeast-two-hybrid and proximity ligation assays, we identified interactions of multiple RBPs with NRL, a Maf-family bZIP transcription factor critical for retinal rod photoreceptor development and function. In addition to splicing, many NRL-interacting RBPs are associated with R-loops, which form during transcription and increase during photoreceptor maturation. Focusing on DHX9 RNA helicase, we demonstrate that its expression is modulated by NRL and that the NRL-DHX9 interaction is positively influenced by R-loops. ssDRIP-Seq analysis reveals both stranded and unstranded R-loops at distinct genomic elements, characterized by active and inactive epigenetic signatures and enriched at neuronal genes. NRL binds to both types of R-loops, suggesting an epigenetically independent function. Our findings suggest additional functions of NRL during transcription and highlight complex interactions among transcription factors, RBPs, and R-loops in regulating photoreceptor gene expression in the mammalian retina.

Evaluation of Klebsiella pneumoniae pathogenicity through holistic gene content analysis

Klebsiella pneumoniae is a Gram-negative bacterium that causes both community- and healthcare-associated infections. Although various virulence factors and highly pathogenic phenotypes have been reported, the pathogenicity of K. pneumoniae is still not fully understood. In this study, we utilized whole-genome sequencing data of 168 clinical K. pneumoniae strains to assess pathogenicity. This work was based on the concept that the genetic composition of individual genomes (referred to as holistic gene content) of the strains may contribute to their pathogenicity. Holistic gene content analysis revealed two distinct groups of K. pneumoniae strains (‘major group’ and ‘minor group’). The minor group included strains with known highly pathogenic clones (ST23, ST375, ST65 and ST86). The minor group had higher rates of capsular genotype K1 and presence of nine specific virulence genes (rmpA, iucA, iutA, irp2, fyuA, ybtS, iroN, allS and clbA) compared to the major group. Pathogenicity was assessed using Galleria mellonella larvae. Infection experiments revealed lower survival rates of larvae infected with strains from the minor group, indicating higher virulence. In addition, the minor group had a higher string test positivity rate than the major group. Holistic gene content analysis predicted possession of virulence genes, string test positivity and pathogenicity as observed in the G. mellonella infection model. Moreover, the findings suggested the presence of as yet unrecognized genomic elements that are either involved in the acquisition of virulence genes or associated with pathogenicity.