From the cauldron of conflict: Endogenous gene regulation by piRNA and other modes of adaptation enabled by selfish transposable elements.

Transposable elements (TEs) provide a prime example of genetic conflict because they can proliferate in genomes and populations even if they harm the host. However, numerous studies have shown that TEs, though typically harmful, can also provide fuel for adaptation. This is because they code functional sequences that can be useful for the host in which they reside. In this review, I summarize the "how" and "why" of adaptation enabled by the genetic conflict between TEs and hosts. In addition, focusing on mechanisms of TE control by small piwi-interacting RNAs (piRNAs), I highlight an indirect form of adaptation enabled by conflict. In this case, mechanisms of host defense that regulate TEs have been redeployed for endogenous gene regulation. I propose that the genetic conflict released by meiosis in early eukaryotes may have been important because, among other reasons, it spurred evolutionary innovation on multiple interwoven trajectories - on the part of hosts and also embedded genetic parasites. This form of evolution may function as a complexity generating engine that was a critical player in eukaryotic evolution.

Editorial: Genome Editing Goes Beyond CRISPR with the Emergence of 'Bridge' RNA Editing.

Therapeutic human gene editing technologies continue to advance, with the endonuclease, clustered regularly interspaced short palindromic repeats (CRISPR) being one of the most rapidly developing technologies. Recently, in 2024, a method of RNA editing called 'bridge editing' has been described in bacteria, which is more powerful and has broader applications than CRISPR to reshape the genome. The term 'bridge editing' is used because the method physically links, or bridges, two sections of DNA and can alter large sections of a genome. 'Bridge editing' relies on insertion sequence (IS) elements, the simplest autonomous transposable elements in prokaryotic genomes. This method provides a unified mechanism for the three fundamental types of DNA rearrangement required for genome design: inversion, insertion, and excision. The 'bridge' recombination system could expand the range and diversity of nucleic acid-guided therapeutic systems beyond RNA interference and CRISPR. This editorial aims to introduce new developments in 'bridge' RNA editing that have the increased potential to reshape the genome.

Local Genomic Surveillance of Invasive Streptococcus pyogenes in Eastern North Carolina (ENC) in 2022-2023.

The recent increase in Group A Streptococcus (GAS) incidences in several countries across Europe and some areas of the Unites States (U.S.) has raised concerns. To understand GAS diversity and prevalence, we conducted a local genomic surveillance in Eastern North Carolina (ENC) in 2022-2023 with 95 isolates and compared its results to those of the existing national genomic surveillance in the U.S. in 2015-2021 with 13,064 isolates. We observed their epidemiological changes before and during the COVID-19 pandemic and detected a unique sub-lineage in ENC among the most common invasive GAS strain, ST28/emm1. We further discovered a multiple-copy insertion sequence, ISLgar5, in ST399/emm77 and its single-copy variants in some other GAS strains. We discovered ISLgar5 was linked to a Tn5801-like tetM-carrying integrative and conjugative element, and its copy number was associated with an ermT-carrying pRW35-like plasmid. The dynamic insertions of ISLgar5 may play a vital role in genome fitness and adaptation, driving GAS evolution relevant to antimicrobial resistance and potentially GAS virulence.

Distinct horizontal transfer mechanisms for type I and type V CRISPR-associated transposons.

CASTs use both CRISPR-associated proteins and Tn7-family transposons for RNA-guided vertical and horizontal transmission. CASTs encode minimal CRISPR arrays but can't acquire new spacers. Here, we report that CASTs can co-opt defense-associated CRISPR arrays for horizontal transmission. A bioinformatic analysis shows that CASTs co-occur with defense-associated CRISPR systems, with the highest prevalence for type I-B and type V CAST sub-types. Using an E. coli quantitative transposition assay and in vitro reconstitution, we show that CASTs can use CRISPR RNAs from these defense systems. A high-resolution structure of the type I-F CAST-Cascade in complex with a type III-B CRISPR RNA reveals that Cas6 recognizes direct repeats via sequence-independent π - π interactions. In addition to using heterologous CRISPR arrays, type V CASTs can also transpose via an unguided mechanism, even when the S15 co-factor is over-expressed. Over-expressing S15 and the trans-activating CRISPR RNA or a single guide RNA reduces, but does not abrogate, off-target integration for type V CASTs. Our findings suggest that some CASTs may exploit defense-associated CRISPR arrays and that this fact must be considered when porting CASTs to heterologous bacterial hosts. More broadly, this work will guide further efforts to engineer the activity and specificity of CASTs for gene editing applications.

Inherited defects of piRNA biogenesis cause transposon de-repression, impaired spermatogenesis, and human male infertility.

piRNAs are crucial for transposon silencing, germ cell maturation, and fertility in male mice. Here, we report on the genetic landscape of piRNA dysfunction in humans and present 39 infertile men carrying biallelic variants in 14 different piRNA pathway genes, including PIWIL1, GTSF1, GPAT2, MAEL, TDRD1, and DDX4. In some affected men, the testicular phenotypes differ from those of the respective knockout mice and range from complete germ cell loss to the production of a few morphologically abnormal sperm. A reduced number of pachytene piRNAs was detected in the testicular tissue of variant carriers, demonstrating impaired piRNA biogenesis. Furthermore, LINE1 expression in spermatogonia links impaired piRNA biogenesis to transposon de-silencing and serves to classify variants as functionally relevant. These results establish the disrupted piRNA pathway as a major cause of human spermatogenic failure and provide insights into transposon silencing in human male germ cells.

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes.

Long-read sequencing, exemplified by PacBio, revolutionizes genomics, overcoming challenges like repetitive sequences. However, the high DNA requirement ( > 1 µg) is prohibitive for small organisms. We develop a low-input (100 ng), low-cost, and amplification-free library-generation method for PacBio sequencing (LILAP) using Tn5-based tagmentation and DNA circularization within one tube. We test LILAP with two Drosophila melanogaster individuals, and generate near-complete genomes, surpassing preexisting single-fly genomes. By analyzing variations in these two genomes, we characterize mutational processes: complex transpositions (transposon insertions together with extra duplications and/or deletions) prefer regions characterized by non-B DNA structures, and gene conversion of transposons occurs on both DNA and RNA levels. Concurrently, we generate two complete assemblies for the endosymbiotic bacterium Wolbachia in these flies and similarly detect transposon conversion. Thus, LILAP promises a broad PacBio sequencing adoption for not only mutational studies of flies and their symbionts but also explorations of other small organisms or precious samples.

A deoxyviolacein-based transposon insertion vector for pigmented tracer studies.

Pigments provide a simple means to rapidly visually ascertain the quantities or presence of specific microbes in a complex community. The selection of pigment-producing colonies that are simple to differentiate from common colony phenotypes provides a high degree of certainty for the identity of pigment-tagged strains. Successful employment of pigment production is dependent on various intrinsic factors related to proper levels of gene expression and pigment production that are not always easy to predict and vary within each microbe. We have constructed a simple transposon system that incorporates the genes for the production of deoxyviolacein, a pigment produced from intracellular reserves of the amino acid tryptophan, to randomly insert these genes throughout the genome. This tool allows the user to select from many thousands of potential sites throughout a bacterial genome for an ideal location to generate the desired amount of pigment. We have applied this system to a small selection of endophytes and other model bacteria to differentiate these strains from complex communities and confirm their presence after several weeks in natural environments. We provide two examples of applications using the pigments to trace strains following introduction into plant tissues or to produce a reporter strain for extracellular nitrogen compound sensing. We recognize that this tool could have far broader utility in other applications and microbes, and describe the methodology for use by the greater scientific community.

DNA methylation enables recurrent endogenization of giant viruses in an animal relative.

5-Methylcytosine (5mC) is a widespread silencing mechanism that controls genomic parasites. In eukaryotes, 5mC has gained complex roles in gene regulation beyond parasite control, yet 5mC has also been lost in many lineages. The causes for 5mC retention and its genomic consequences are still poorly understood. Here, we show that the protist closely related to animals Amoebidium appalachense features both transposon and gene body methylation, a pattern reminiscent of invertebrates and plants. Unexpectedly, hypermethylated genomic regions in Amoebidium derive from viral insertions, including hundreds of endogenized giant viruses, contributing 14% of the proteome. Using a combination of inhibitors and genomic assays, we demonstrate that 5mC silences these giant virus insertions. Moreover, alternative Amoebidium isolates show polymorphic giant virus insertions, highlighting a dynamic process of infection, endogenization, and purging. Our results indicate that 5mC is critical for the controlled coexistence of newly acquired viral DNA into eukaryotic genomes, making Amoebidium a unique model to understand the hybrid origins of eukaryotic DNA.

Copy number variation introduced by a massive mobile element facilitates global thermal adaptation in a fungal wheat pathogen.

Copy number variation (CNV) can drive rapid evolution in changing environments. In microbial pathogens, such adaptation is a key factor underpinning epidemics and colonization of new niches. However, the genomic determinants of such adaptation remain poorly understood. Here, we systematically investigate CNVs in a large genome sequencing dataset spanning a worldwide collection of 1104 genomes from the major wheat pathogen Zymoseptoria tritici. We found overall strong purifying selection acting on most CNVs. Genomic defense mechanisms likely accelerated gene loss over episodes of continental colonization. Local adaptation along climatic gradients was likely facilitated by CNVs affecting secondary metabolite production and gene loss in general. One of the strongest loci for climatic adaptation is a highly conserved gene of the NAD-dependent Sirtuin family. The Sirtuin CNV locus localizes to an ~68-kb Starship mobile element unique to the species carrying genes highly expressed during plant infection. The element has likely lost the ability to transpose, demonstrating how the ongoing domestication of cargo-carrying selfish elements can contribute to selectable variation within populations. Our work highlights how standing variation in gene copy numbers at the global scale can be a major factor driving climatic and metabolic adaptation in microbial species.

Rapid evolution of promoters from germline-specifically expressed genes including transposon silencing factors.

BACKGROUND: The piRNA pathway in animal gonads functions as an 'RNA-based immune system', serving to silence transposable elements and prevent inheritance of novel invaders. In Drosophila, this pathway relies on three gonad-specific Argonaute proteins (Argonaute-3, Aubergine and Piwi) that associate with 23-28 nucleotide piRNAs, directing the silencing of transposon-derived transcripts. Transposons constitute a primary driver of genome evolution, yet the evolution of piRNA pathway factors has not received in-depth exploration. Specifically, channel nuclear pore proteins, which impact piRNA processing, exhibit regions of rapid evolution in their promoters. Consequently, the question arises whether such a mode of evolution is a general feature of transposon silencing pathways. RESULTS: By employing genomic analysis of coding and promoter regions within genes that function in transposon silencing in Drosophila, we demonstrate that the promoters of germ cell-specific piRNA factors are undergoing rapid evolution. Our findings indicate that rapid promoter evolution is a common trait among piRNA factors engaged in germline silencing across insect species, potentially contributing to gene expression divergence in closely related taxa. Furthermore, we observe that the promoters of genes exclusively expressed in germ cells generally exhibit rapid evolution, with some divergence in gene expression. CONCLUSION: Our results suggest that increased germline promoter evolution, in partnership with other factors, could contribute to transposon silencing and evolution of species through differential expression of genes driven by invading transposons.

High-Throughput Mapping of Chromosomal Conformations in E. coli Under Physiological Conditions Using Massively Multiplexed Mu Transposition.

Many approaches for measuring three-dimensional chromosomal conformations rely upon formaldehyde crosslinking followed by subsequent proximity ligation, a family of methods exemplified by 3C, Hi-C, etc. Here we provide an alternative crosslinking-free procedure for high-throughput identification of long-range contacts in the chromosomes of enterobacteria, making use of contact-dependent transposition of phage Mu to identify distant loci in close contact. The procedure described here will suffice to provide a comprehensive map of transposition frequencies between tens of thousands of loci in a bacterial genome, with the resolution limited by the diversity of the insertion site library used and the sequencing depth applied.

Comprehensive analysis of the Xya riparia genome uncovers the dominance of DNA transposons, LTR/Gypsy elements, and their evolutionary dynamics.

Transposable elements (TEs) are DNA sequences that can move or replicate within a genome, and their study has become increasingly important in understanding genome evolution and function. The Tridactylidae family, including Xya riparia (pygmy mole cricket), harbors a variety of transposable elements (TEs) that have been insufficiently investigated. Further research is required to fully understand their diversity and evolutionary characteristics. Hence, we conducted a comprehensive repeatome analysis of X. riparia species using the chromosome-level assembled genome. The study aimed to comprehensively analyze the abundance, distribution, and age of transposable elements (TEs) in the genome. The results indicated that the genome was 1.67 Gb, with 731.63 Mb of repetitive sequences, comprising 27% of Class II (443.25 Mb), 16% of Class I (268.45 Mb), and 1% of unknown TEs (19.92 Mb). The study found that DNA transposons dominate the genome, accounting for approximately 60% of the total repeat size, with retrotransposons and unknown elements accounting for 37% and 3% of the genome, respectively. The members of the Gypsy superfamily were the most abundant amongst retrotransposons, accounting for 63% of them. The transposable superfamilies (LTR/Gypsy, DNA/nMITE, DNA/hAT, and DNA/Helitron) collectively constituted almost 70% of the total repeat size of all six chromosomes. The study further unveiled a significant linear correlation (Pearson correlation: r = 0.99, p-value = 0.00003) between the size of the chromosomes and the repetitive sequences. The average age of DNA transposon and retrotransposon insertions ranges from 25 My (million years) to 5 My. The satellitome analysis discovered 13 satellite DNA families that comprise about 0.15% of the entire genome. In addition, the transcriptional analysis of TEs found that DNA transposons were more transcriptionally active than retrotransposons. Overall, the study suggests that the genome of X. riparia is complex, characterized by a substantial portion of repetitive elements. These findings not only enhance our understanding of TE evolution within the Tridactylidae family but also provide a foundation for future investigations into the genomic intricacies of related species.

Structural Differences between the Genomes of Deinococcus radiodurans Strains from Different Laboratories.

The bacterium Deinococcus radiodurans is known to efficiently and accurately reassemble its genome after hundreds of DNA double-strand breaks (DSBs). Only at very large amounts of radiation-induced DSBs is this accuracy affected in the wild-type D. radiodurans, causing rearrangements in its genome structure. However, changes in its genome structure may also be possible during the propagation and storage of cell cultures. We investigate this possibility by listing structural differences between three completely sequenced genomes of D. radiodurans strains with a recent common ancestor-the type strain stored and sequenced in two different laboratories (of the ATCC 13939 lineage) and the first sequenced strain historically used as the reference (ATCC BAA-816). We detected a number of structural differences and found the most likely mechanisms behind them: (i) transposition/copy number change in mobile interspersed repeats-insertion sequences and small non-coding repeats, (ii) variable number of monomers within tandem repeats, (iii) deletions between long direct DNA repeats, and (iv) deletions between short (4-10 bp) direct DNA repeats. The most surprising finding was the deletions between short repeats because it indicates the utilization of a less accurate DSB repair mechanism in conditions in which a more accurate one should be both available and preferred. The detected structural differences, as well as SNPs and short indels, while being important footprints of deinococcal DNA metabolism and repair, are also a valuable resource for researchers using these D. radiodurans strains.

Transposable elements-mediated recruitment of KDM1A epigenetically silences HNF4A expression to promote hepatocellular carcinoma.

Transposable elements (TEs) contribute to gene expression regulation by acting as cis-regulatory elements that attract transcription factors and epigenetic regulators. This research aims to explore the functional and clinical implications of transposable element-related molecular events in hepatocellular carcinoma, focusing on the mechanism through which liver-specific accessible TEs (liver-TEs) regulate adjacent gene expression. Our findings reveal that the expression of HNF4A is inversely regulated by proximate liver-TEs, which facilitates liver cancer cell proliferation. Mechanistically, liver-TEs are predominantly occupied by the histone demethylase, KDM1A. KDM1A negatively influences the methylation of histone H3 Lys4 (H3K4) of liver-TEs, resulting in the epigenetic silencing of HNF4A expression. The suppression of HNF4A mediated by KDM1A promotes liver cancer cell proliferation. In conclusion, this study uncovers a liver-TE/KDM1A/HNF4A regulatory axis that promotes liver cancer growth and highlights KDM1A as a promising therapeutic target. Our findings provide insight into the transposable element-related molecular mechanisms underlying liver cancer progression.

HiTE: a fast and accurate dynamic boundary adjustment approach for full-length transposable element detection and annotation.

Recent advancements in genome assembly have greatly improved the prospects for comprehensive annotation of Transposable Elements (TEs). However, existing methods for TE annotation using genome assemblies suffer from limited accuracy and robustness, requiring extensive manual editing. In addition, the currently available gold-standard TE databases are not comprehensive, even for extensively studied species, highlighting the critical need for an automated TE detection method to supplement existing repositories. In this study, we introduce HiTE, a fast and accurate dynamic boundary adjustment approach designed to detect full-length TEs. The experimental results demonstrate that HiTE outperforms RepeatModeler2, the state-of-the-art tool, across various species. Furthermore, HiTE has identified numerous novel transposons with well-defined structures containing protein-coding domains, some of which are directly inserted within crucial genes, leading to direct alterations in gene expression. A Nextflow version of HiTE is also available, with enhanced parallelism, reproducibility, and portability.

Comparison of Bordetella species identification among differing rt-PCR assays in the United States.

In the United States, the general laboratory method for diagnosing pertussis, caused by Bordetella pertussis, is real-time PCR (rt-PCR) targeting insertion sequence 481 (IS481). Other Bordetella species (parapertussis, holmesii, and bronchiseptica) can also cause a pertussis-like syndrome, and some commercial laboratory assays include the insertion sequence 1001 (pIS1001) that can detect B. parapertussis/B. bronchiseptica (BppBb). Because IS481 exists in B. pertussis and B. holmesii, current commercial assays cannot differentiate these two species. We used a multiplex rt-PCR assay containing species-specific targets to Bordetella to evaluate clinical specimens detected as B. pertussis/B. holmesii (BpBh) or BppBb by commercial laboratories. A sample of 3,984 clinical specimens positive for IS481 or pIS1001 from two commercial laboratories during 2012-2019 were re-tested at CDC. Agreement of Bordetella species between the CDC and commercial laboratory assays, and the proportion of commercial laboratory specimens that were non-B. pertussis by CDC's assay was assessed. Overall agreement in Bordetella species detection and identification between the CDC and commercial lab assays was 85%. Agreement for identifying B. pertussis was 87% for 3,663 BpBh specimens and 98% for identifying B. parapertussis in 310 BppBb specimens. CDC's assay detected B. holmesii in 55/3,984 (1.4%) specimens. Most discrepant results (410/490, 82%) were BpBh specimens interpreted as indeterminate B. pertussis at CDC. We found a small portion of B. holmesii in a sample of IS481-positive clinical specimens originally identified by commercial laboratory rt-PCR assays, suggesting that commercial PCR assays are a reliable diagnostic tool for correctly identifying Bordetella species in most patients with suspected pertussis. IMPORTANCE: When testing specimens collected from patients with suspected pertussis, large-scale commercial laboratories in the United States employ an IS481-based assay that cannot differentiate between Bordetella pertussis and Bordetella holmseii. The level of B. holmesii causing pertussis-like illness in the United States is not well-understood given that only B. pertussis is nationally notifiable. After re-testing with a multiplex, species-specific rt-PCR assay, our data show low levels of B. holmesii identified in a sample of IS481-positive clinical specimens originally identified by commercial laboratory rt-PCR assays. These results reinforce the validity of large-scale commercial rt-PCR testing as a reliable diagnostic tool for pertussis in the United States.

Association of a rapidly selected 4.3kb transposon-containing structural variation with a P450-based resistance to pyrethroids in the African malaria vector Anopheles funestus.

Deciphering the evolutionary forces controlling insecticide resistance in malaria vectors remains a prerequisite to designing molecular tools to detect and assess resistance impact on control tools. Here, we demonstrate that a 4.3kb transposon-containing structural variation is associated with pyrethroid resistance in central/eastern African populations of the malaria vector Anopheles funestus. In this study, we analysed Pooled template sequencing data and direct sequencing to identify an insertion of 4.3kb containing a putative retro-transposon in the intergenic region of two P450s CYP6P5-CYP6P9b in mosquitoes of the malaria vector Anopheles funestus from Uganda. We then designed a PCR assay to track its spread temporally and regionally and decipher its role in insecticide resistance. The insertion originates in or near Uganda in East Africa, where it is fixed and has spread to high frequencies in the Central African nation of Cameroon but is still at low frequency in West Africa and absent in Southern Africa. A marked and rapid selection was observed with the 4.3kb-SV frequency increasing from 3% in 2014 to 98% in 2021 in Cameroon. A strong association was established between this SV and pyrethroid resistance in field populations and is reducing pyrethroid-only nets' efficacy. Genetic crosses and qRT-PCR revealed that this SV enhances the expression of CYP6P9a/b but not CYP6P5. Within this structural variant (SV), we identified putative binding sites for transcription factors associated with the regulation of detoxification genes. An inverse correlation was observed between the 4.3kb SV and malaria parasite infection, indicating that mosquitoes lacking the 4.3kb SV were more frequently infected compared to those possessing it. Our findings highlight the underexplored role and rapid spread of SVs in the evolution of insecticide resistance and provide additional tools for molecular surveillance of insecticide resistance.

Double blaKPC-2 copies quadrupled minimum inhibitory concentration of ceftazidime-avibactam in hospital-derived Klebsiella pneumoniae.

To illustrate the genomic and drug resistance traits of the Klebsiella pneumoniae Kpn_XM9, which harbors a transposon (Tn) As1 and was barely susceptible to ceftazidime-avibactam (CZA). Whole-genome sequencing, gene deletion, antimicrobial susceptibility, and conjugation tests were carried out to illustrate the traits of Kpn_XM9. As confirmed by whole-genome sequencing, the Kpn_XM9 harbored a 5,523,536 bp chromosome and five plasmids with lengths being 128,129, 196,512, 84,812, 43,695, and 5,596 bp, respectively. Plasmid p1_Kpn_XM9 (128,219 bp) contained four resistance genes, blaCTX-M-65, blaTEM-1B, rmtB, and two copies of blaKPC-2. Genes blaKPC-2 were bracketed by ISKpn17 and ISKpn16 within a new composite Tn3-like TnAs1. The two tandem repeats, positioned opposite each other, were spaced 93,447 bp apart in p1_Kpn_XM9. Kpn_XM9 belonged to K64 and sequence type (ST) 11. The Kpn_XM9 was resistant to amikacin, aztreonam, ticarcillin/clavulanic acid, piperacillin/tazobactam, ceftazidime, cefepime, imipenem, meropenem, tobramycin, ciprofloxacin, levofloxacin, doxycycline, minocycline, tigecycline, colistin, and trimethoprim/sulfamethoxazole; it was barely susceptible to CZA with a minimum inhibitory concentration of 8/4 µg/mL, which declined to 2/4 µg/mL after a 18,555 bp nucleotide was knocked out and one copy of blaKPC-2 was sustained on p1_Kpn_XM9. Kpn_XM9 had virulence genes encoding Types 1 and 3 fimbriae, four siderophores, and capsular polysaccharide anchoring protein but no genes upregulating capsular polysaccharide synthesis. The Kpn_XM9 presented a classical phenotype with extreme drug resistance. The emergence of double copies of blaKPC-2 in a single plasmid from the predominant ST11 K. pneumoniae represents a new therapeutic challenge.IMPORTANCEWith the wide use of ceftazidime-avibactam against carbapenem-resistant organisms, its resistance is increasingly documented; among the corresponding resistance mechanisms, mutations of blaKPC-2 or blaKPC-3 into other subtypes are dominant to date. However, more copies of blaKPC-2 may also greatly increase the minimum inhibitory concentration of ceftazidime-avibactam, which could be conferred by transposon As1 and insertion sequence 26 and should be of concern.

Pangenome evaluation of gene essentiality in Streptococcus pyogenes.

Bacterial species often consist of strains with variable gene content, collectively referred to as the pangenome. Variations in the genetic makeup of strains can alter bacterial physiology and fitness. To define biologically relevant genes of a genome, genome-wide transposon mutant libraries have been used to identify genes essential for survival or virulence in a given strain. Such phenotypic studies have been conducted in four different genotypes of the human pathogen Streptococcus pyogenes, yet challenges exist in comparing results across studies conducted in different genetic backgrounds and conditions. To advance genotype to phenotype inferences across different S. pyogenes strains, we built a pangenome database of 249 S. pyogenes reference genomes. We systematically re-analyzed publicly available transposon sequencing datasets from S. pyogenes using a transposon sequencing-specific analysis pipeline, Transit. Across four genetic backgrounds and nine phenotypic conditions, 355 genes were essential for survival, corresponding to ~24% of the core genome. Clusters of Orthologous Genes (COG) categories related to coenzyme and lipid transport and growth functions were overrepresented as essential. Finally, essential operons across S. pyogenes genotypes were defined, with an increased number of essential operons detected under in vivo conditions. This study provides an extendible database to which new studies can be added, and a searchable html-based resource to direct future investigations into S. pyogenes biology.IMPORTANCEStreptococcus pyogenes is a human-adapted pathogen occupying restricted ecological niches. Understanding the essentiality of genes across different strains and experimental conditions is important to direct research questions and efforts to prevent the large burden of disease caused by S. pyogenes. To this end we systematically reanalyzed transposon sequencing studies in S. pyogenes using transposon sequencing-specific methods, integrating them into an extendible meta-analysis framework. This provides a repository of gene essentiality in S. pyogenes which was used to highlight specific genes of interest and for the community to guide future phenotypic studies.