Transposon Insertion Research Articles

Identification of two mutant JASON-RELATED genes associated with unreduced pollen production in potato

Key messageMultiple QTLs control unreduced pollen production in potato. Two major-effect QTLs co-locate with mutant alleles of genes with homology to AtJAS, a known regulator of meiotic spindle orientation.In diploid potato the production of unreduced gametes with a diploid (2n) rather than a haploid (n) number of chromosomes has been widely reported. Besides their evolutionary important role in sexual polyploidisation, unreduced gametes also have a practical value for potato breeding as a bridge between diploid and tetraploid germplasm. Although early articles argued for a monogenic recessive inheritance, the genetic basis of unreduced pollen production in potato has remained elusive. Here, three diploid full-sib populations were genotyped with an amplicon sequencing approach and phenotyped for unreduced pollen production across two growing seasons. We identified two minor-effect and three major-effect QTLs regulating this trait. The two QTLs with the largest effect displayed a recessive inheritance and an additive interaction. Both QTLs co-localised with genes encoding for putative AtJAS homologs, a key regulator of meiosis II spindle orientation in Arabidopsis thaliana. The function of these candidate genes is consistent with the cytological phenotype of mis-oriented metaphase II plates observed in the parental clones. The alleles associated with elevated levels of unreduced pollen showed deleterious mutation events: an exonic transposon insert causing a premature stop, and an amino acid change within a highly conserved domain. Taken together, our findings shed light on the natural variation underlying unreduced pollen production in potato and will facilitate interploidy breeding by enabling marker-assisted selection for this trait.

CRISPR RNA-Guided Transposases Facilitate Dispensable Gene Study in Phage.

Phages provide a potential therapy for multi-drug-resistant (MDR) bacteria. However, a significant portion of viral genes often remains unknown, posing potential dangers. The identification of non-essential genes helps dissect and simplify phage genomes, but current methods have various limitations. In this study, we present an in vivo two-plasmid transposon insertion system to assess the importance of phage genes, which is based on the V. cholerae transposon Tn6677, encoding a nuclease-deficient type I-F CRISPR-Cas system. We first validated the system in Pseudomonas aeruginosa PAO1 and its phage S1. We then used the selection marker AcrVA1 to protect transposon-inserted phages from CRISPR-Cas12a and enriched the transposon-inserted phages. For a pool of selected 10 open-reading frames (2 known functional protein genes and 8 hypothetical protein genes) of phage S1, we identified 5 (2 known functional protein genes and 3 hypothetical protein genes) as indispensable genes and the remaining 5 (all hypothetical protein genes) as dispensable genes. This approach offers a convenient, site-specific method that does not depend on homologous arms and double-strand breaks (DSBs), holding promise for future applications across a broader range of phages and facilitating the identification of the importance of phage genes and the insertion of genetic cargos.

Genome-wide analysis of acid tolerance genes of Enterococcus faecalis with RNA-seq and Tn-seq

Enterococcus faecalis, a formidable nosocomial and community-acquired opportunistic pathogen, can persist a wide range of extreme environments, including low pH and nutrient deficiency. Clarifying the survival mechanism of E. faecalis in low-pH conditions is the key to combating the infectious diseases caused by E. faecalis. In this study, we combined transcriptome profiling (RNA-seq) and transposon insertion sequencing (TIS) to comprehensively understand the genes that confer these features on E. faecalis. The metadata showed that genes whose products are involved in cation transportation and amino acid biosynthesis were predominantly differentially expressed under acid conditions. The products of genes such as opp1C and copY reduced the hydrion concentration in the cell, whereas those of gldA2, gnd2, ubiD, and ubiD2 mainly participated in amino metabolism, increasing matters to neutralize excess acid. These, together with the folE and hexB genes, which are involved in mismatch repair, form a network of E. faecalis genes necessary for its survival under acid conditions.

ProTInSeq: transposon insertion tracking by ultra-deep DNA sequencing to identify translated large and small ORFs

Identifying open reading frames (ORFs) being translated is not a trivial task. ProTInSeq is a technique designed to characterize proteomes by sequencing transposon insertions engineered to express a selection marker when they occur in-frame within a protein-coding gene. In the bacterium Mycoplasma pneumoniae, ProTInSeq identifies 83% of its annotated proteins, along with 5 proteins and 153 small ORF-encoded proteins (SEPs; ≤100 aa) that were not previously annotated. Moreover, ProTInSeq can be utilized for detecting translational noise, as well as for relative quantification and transmembrane topology estimation of fitness and non-essential proteins. By integrating various identification approaches, the number of initially annotated SEPs in this bacterium increases from 27 to 329, with a quarter of them predicted to possess antimicrobial potential. Herein, we describe a methodology complementary to Ribo-Seq and mass spectroscopy that can identify SEPs while providing other insights in a proteome with a flexible and cost-effective DNA ultra-deep sequencing approach.

Abstract PR-011: The role of tumor microenvironment in metastasis and relapses of medulloblastoma using immunocompetent mouse models

Abstract Medulloblastoma (MB) is an embryonal malignant tumor that arises in the cerebellum of infants, children, and adolescents. Despite highly aggressive therapies, including surgery, radiation, and chemotherapy, metastasis remains the primary cause of death in these patients. The tumor microenvironment (TME) has an important role in metastasis and relapse in several cancer types, but it’s role in high-risk MBs (metastatic and relapse) has been understudied. Here, we investigated the process of relapse in an immunocompetent Shh MB mouse model based on somatic mutagenesis induced by Sleeping Beauty (SB)-mediated transposition on a genetically engineered strain carrying a loss-of-function Patched 1 gene (Ptch1) mutation. Lack of one functional copy of the Ptch1 gene predisposes this strain to develop MB. Somatic mutations induced by the transposon increase penetrance of tumor formation in the triple mutant mice carrying the SB transposon system and the Ptch1 mutation. We used this model to study the mechanisms of chemoresistance to cisplatin and cyclophosphamide, that constitute the backbone of the high-dose chemotherapy used to treat patients under the age of 3, who cannot receive craniospinal irradiation. By treating the SB mouse model with surgery followed by high-dose fractionated chemotherapy, we developed a model of MB chemoresistance which is relevant to the infant SHH-MB patient population. While chemotherapy improved mouse survival, only a small number of animals were cured, mirroring the outcome seen in infant patients with Shh MB. Deep sequencing and mapping of transposon-gDNA junctions, along with gene-centric common insertion site analysis (gCIS), revealed that Trp53 transposon insertions are found exclusively at relapse. In the metastatic compartment we observed convergence of gCISes on genes coding for epigenetic regulators, among them a prominent enrichment was found for PRC2 complex genes. RNA-seq analysis on the transcriptome of primary and recurrent tumors indicated that a majority of the differentially expressed genes are associated with inflammatory responses, IL6, and TNF-α/NF-κB signaling. Immunohistochemistry to visualize the expression levels of inflammatory markers: IL-6, CCL2, and F4-80 and spatial targeted multiplexed assays (Geomx) confirmed the increase of immune cells at progression. ChIP-qPCR on cytokine loci revealed binding of PRC2 and deposition of PRC2 mediated repressive histone mark H3K27me3 on the promoter of CCL2 highlighting LOF mutations in PRC2 components as potential epigenetic mechanisms for immune recruitment at progression. Inhibition of immune activation through pharmacological strategies might be a possible therapeutic option for infant SHH-MB. Citation Format: Niusha Khazaei, Ana Castillo Orozco, Wajih Jawhar, Geoffroy Danieau, Kiran Narta, Michael Taylor, Sorana Morrissy, Ana Stuecklin, Livia Garzia. The role of tumor microenvironment in metastasis and relapses of medulloblastoma using immunocompetent mouse models [abstract]. In: Proceedings of the AACR Special Conference on Brain Cancer; 2023 Oct 19-22; Minneapolis, Minnesota. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_1):Abstract nr PR-011.

A CRISPR interference system for engineering biological nitrogen fixation.

A grand challenge for the next century is in facing a changing climate through bioengineering solutions. Biological nitrogen fixation, the globally consequential, nitrogenase-catalyzed reduction of atmospheric nitrogen to bioavailable ammonia, is a vital area of focus. Nitrogen fixation engineering relies upon extensive understanding of underlying genetics in microbial models, including the broadly utilized gammaproteobacterium, Azotobacter vinelandii (A. vinelandii). Here, we report the first CRISPR interference (CRISPRi) system for targeted gene silencing in A. vinelandii that integrates genomically via site-specific transposon insertion. We demonstrate that CRISPRi can repress transcription of an essential nitrogen fixation gene by ~60%. Further, we show that nitrogenase genes are suitably expressed from the transposon insertion site, indicating that CRISPRi and engineered nitrogen fixation genes can be co-integrated for combinatorial studies of gene expression and engineering. Our established CRISPRi system fills an important gap for engineering microbial nitrogen fixation for desired purposes.IMPORTANCEAll life on Earth requires nitrogen to survive. About 78% of the atmosphere alone is nitrogen, yet humans cannot use it directly. Instead, we obtain the nitrogen we need for our survival through the food we eat. For more than 100 years, a substantial portion of agricultural productivity has relied on industrial methods for nitrogen fertilizer synthesis, which consumes significant amounts of nonrenewable energy resources and exacerbates environmental degradation and human-induced climate change. Promising alternatives to these industrial methods rely on engineering the only biological pathway for generating bioaccessible nitrogen: microbial nitrogen fixation. Bioengineering strategies require an extensive understanding of underlying genetics in nitrogen-fixing microbes, but genetic tools for this critical goal remain lacking. The CRISPRi gene silencing system that we report, developed in the broadly utilized nitrogen-fixing bacterial model, Azotobacter vinelandii, is an important step toward elucidating the complexity of nitrogen fixation genetics and enabling their manipulation.

Regulatory mechanism of heat-active retrotransposons by the SET domain protein SUVH2.

New transposon insertions are deleterious to genome stability. The RNA-directed DNA methylation (RdDM) pathway evolved to regulate transposon activity via DNA methylation. However, current studies have not yet clearly described the transposition regulation. ONSEN is a heat-activated retrotransposon that is activated at 37°C. The plant-specific SUPPRESSOR OF VARIEGATION 3-9 HOMOLOG (SUVH) family proteins function downstream of the RdDM pathway. The SUVH protein families are linked to TE silencing by two pathways, one through DNA methylation and the other through chromatin remodeling. In this study, we analyzed the regulation of ONSEN activity by SUVH2. We observed that ONSEN transcripts were increased; however, there was no transpositional activity in Arabidopsis suvh2 mutant. The suvh2 mutant produced siRNAs from the ONSEN locus under heat stress, suggesting that siRNAs are involved in suppressing transposition. These results provide new insights into the regulatory mechanisms of retrotransposons that involve siRNA in the RdDM pathway.

Machine learning analysis of RB-TnSeq fitness data predicts functional gene modules in Pseudomonas putida KT2440.

This study demonstrates a rapid, automated approach for elucidating functional modules within complex genetic networks. While Pseudomonas putida randomly barcoded transposon insertion sequencing data were used as a proof of concept, this approach is applicable to any organism with existing functional genomics data sets and may serve as a useful tool for many valuable applications, such as guiding metabolic engineering efforts in other microbes or understanding functional relationships between virulence-associated genes in pathogenic microbes. Furthermore, this work demonstrates that comparison of data obtained from independent component analysis of transcriptomics and gene fitness datasets can elucidate regulatory-functional relationships between genes, which may have utility in a variety of applications, such as metabolic modeling, strain engineering, or identification of antimicrobial drug targets.

The maize striate leaves2 (sr2) gene encodes a conserved DUF3732 domain and is homologous to the rice yss1 gene.

Maize striate leaves2 (sr2) is a mutant that causes white stripes on leaves that has been used in mapping studies for decades though the underlying gene has not been identified. The sr2 locus has been previously mapped to small regions of normal chromosome 10 (N10) and a rearranged variant called abnormal chromosome 10 (Ab10). A comparison of assembled genomes carrying N10 and Ab10 revealed only five candidate sr2 genes. Analysis of a stock carrying the sr2 reference allele (sr2-ref) showed that one of the five genes has a transposon insertion that disrupts its protein sequence and has a severe reduction in mRNA. An independent Mutator transposon insertion in the gene (sr2-Mu) failed to complement the sr2-ref mutation, and plants homozygous for sr2-Mu showed white striped leaf margins. The sr2 gene encodes a DUF3732 protein with strong homology to a rice gene with a similar mutant phenotype called young seedling stripe1 (yss1). These and other published data suggest that sr2 may have a function in plastid gene expression.

MBARq: a versatile and user-friendly framework for the analysis of DNA barcodes from transposon insertion libraries, knockout mutants, and isogenic strain populations.

DNA barcoding has become a powerful tool for assessing the fitness of strains in a variety of studies, including random transposon mutagenesis screens, attenuation of site-directed mutants, and population dynamics of isogenic strain pools. However, the statistical analysis, visualization, and contextualization of the data resulting from such experiments can be complex and require bioinformatic skills. Here, we developed mBARq, a user-friendly tool designed to simplify these steps for diverse experimental setups. The tool is seamlessly integrated with an intuitive web app for interactive data exploration via the STRING and KEGG databases to accelerate scientific discovery. The tool is implemented in Python. The source code is freely available (https://github.com/MicrobiologyETHZ/mbarq) and the web app can be accessed at: https://microbiomics.io/tools/mbarq-app.

Determination of the Allelic Composition of the sdw1/denso (HvGA20ox2), uzu1 (HvBRI1) and ari-e (HvDep1) Genes in Spring Barley Accessions from the VIR Collection

The lodging of barley significantly limits its potential yield, leads to the deterioration of grain quality, and complicates mechanized harvesting. More than 30 dwarfness and semi-dwarfness genes and loci are known for barley, and their involvement in breeding can solve the problem of lodging. The most common dwarfing alleles are of the genes sdw1/denso (HvGA20ox2), uzu1 (HvBRI1), and ari-e (HvDep1). The aim of this study was the design of dCAPS markers for the sdw1.c and ari-e.GP alleles and the molecular screening of barley accessions from the VIR collection for identifying these and other dwarfing alleles commonly used in breeding. Two dCAPS markers have been developed to identify the sdw1.c allele of the HvGA20ox2 gene and ari-e.GP of HvDep1. These dCAPS markers and two known from the literature CAPS and dCAPS markers of the alleles sdw1.a/sdw1.e, sdw1.c, sdw1.d, and uzu1.a were used in the molecular screening of 32 height-contrasting barley accessions. This made it possible to identify the accessions with alleles sdw1.a/sdw1.e, sdw1.c, and sdw1.d of the HvGA20ox2 gene, as well as accessions with a combination of sdw1.c and uzu1.a alleles of the genes HvGA20ox2 and HvBRI1. A comparison of the results of genotyping and phenotyping showed that the presence of dwarfing alleles in all genotypes determines high or medium lodging resistance regardless of the influence of weather conditions. Twelve accessions were found to contain the new allele sdw1.ins of the HvGA20ox2 gene, which differs from the known allele sdw1.c by a larger size of PCR products. It is characterized by the Thalos_2 transposon insertion; the subsequent GTTA insertion, common with the sdw1.c allele; and by a single-nucleotide G→A substitution at the 165 position.

Genetic basis of Arabidopsis thaliana responses to infection by naïve and adapted isolates of turnip mosaic virus.

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

Genetic basis of Arabidopsis thaliana responses to infection by naïve and adapted isolates of turnip mosaic virus

Plant viruses account for enormous agricultural losses worldwide, and the most effective way to combat them is to identify genetic material conferring plant resistance to these pathogens. Aiming to identify genetic associations with responses to infection, we screened a large panel of Arabidopsis thaliana natural inbred lines for four disease-related traits caused by infection by A. thaliana-naïve and -adapted isolates of the natural pathogen turnip mosaic virus (TuMV). We detected a strong, replicable association in a 1.5 Mb region on chromosome 2 with a 10-fold increase in relative risk of systemic necrosis. The region contains several plausible causal genes as well as abundant structural variation, including an insertion of a Copia transposon into a Toll/interleukin receptor (TIR-NBS-LRR) coding for a gene involved in defense, that could be either a driver or a consequence of the disease-resistance locus. When inoculated with TuMV, loss-of-function mutant plants of this gene exhibited different symptoms than wild-type plants. The direction and severity of symptom differences depended on the adaptation history of the virus. This increase in symptom severity was specific for infections with the adapted isolate. Necrosis-associated alleles are found worldwide, and their distribution is consistent with a trade-off between resistance during viral outbreaks and a cost of resistance otherwise, leading to negative frequency-dependent selection.

Genome-wide profiling of piggyBac transposon insertion mutants reveals loss of the F1F0 ATPase complex causes fluconazole resistance in Candida glabrata.

Invasive candidiasis caused by non-albicans species has been on the rise, with Candida glabrata emerging as the second most common etiological agent. Candida glabrata possesses an intrinsically lower susceptibility to azoles and an alarming propensity to rapidly develop high-level azole resistance during treatment. In this study, we have developed an efficient piggyBac (PB) transposon-mediated mutagenesis system in C. glabrata to conduct genome-wide genetic screens and applied it to profile genes that contribute to azole resistance. When challenged with the antifungal drug fluconazole, PB insertion into 270 genes led to significant resistance. A large subset of these genes has a role in the mitochondria, including almost all genes encoding the subunits of the F1F0 ATPase complex. We show that deleting ATP3 or ATP22 results in increased azole resistance but does not affect susceptibility to polyenes and echinocandins. The increased azole resistance is due to increased expression of PDR1 that encodes a transcription factor known to promote drug efflux pump expression. Deleting PDR1 in the atp3Δ or atp22Δ mutant resulted in hypersensitivity to fluconazole. Our results shed light on the mechanisms contributing to azole resistance in C. glabrata. This PB transposon-mediated mutagenesis system can significantly facilitate future genome-wide genetic screens.

Enhancing the Efficacy of Chloramphenicol Therapy for Escherichia coli by Targeting the Secondary Resistome.

The increasing prevalence of antimicrobial resistance and the limited availability of new antimicrobial agents have created an urgent need for new approaches to combat these issues. One such approach involves reevaluating the use of old antibiotics to ensure their appropriate usage and maximize their effectiveness, as older antibiotics could help alleviate the burden on newer agents. An example of such an antibiotic is chloramphenicol (CHL), which is rarely used due to its hematological toxicity. In the current study, we employed a previously published transposon mutant library in MG1655/pTF2::blaCTX-M-1, containing over 315,000 unique transposon insertions, to identify the genetic factors that play an important role during growth in the presence of CHL. The list of conditionally essential genes, collectively referred to as the secondary resistome (SR), included 67 genes. To validate our findings, we conducted gene knockout experiments on six genes: arcA, hfq, acrZ, cls, mdfA, and nlpI. Deleting these genes resulted in increased susceptibility to CHL as demonstrated by MIC estimations and growth experiments, suggesting that targeting the products encoded from these genes may reduce the dose of CHL needed for treatment and hence reduce the toxicity associated with CHL treatment. Thus, the gene products are indicated as targets for antibiotic adjuvants to favor the use of CHL in modern medicine.

Fitness factors impacting survival of a subsurface bacterium in contaminated groundwater.

Many factors contribute to the ability of a microbial species to persist when encountering complexly contaminated environments, including time of exposure, the nature and concentration of contaminants, availability of nutritional resources, and possession of a combination of appropriate molecular mechanisms needed for survival. Herein we sought to identify genes that are most important for survival of Gram-negative Enterobacteriaceae in contaminated groundwater environments containing high concentrations of nitrate and metals using the metal-tolerant Oak Ridge Reservation isolate, Pantoea sp. MT58 (MT58). Survival fitness experiments in which a randomly barcoded transposon insertion (RB-TnSeq) library of MT58 was exposed directly to contaminated Oak Ridge Reservation groundwater samples from across a nitrate and mixed metal contamination plume were used to identify genes important for survival with increasing exposure times and concentrations of contaminants, and availability of a carbon source. Genes involved in controlling and using carbon, encoding transcriptional regulators, and related to Gram-negative outer membrane processes were among those found to be important for survival in contaminated Oak Ridge Reservation groundwater. A comparative genomics analysis of 75 Pantoea genus strains allowed us to further separate the survival determinants into core and non-core genes in the Pantoea pangenome, revealing insights into the survival of subsurface microorganisms during contaminant plume intrusion.

Genome-wide screen of genetic determinants that govern Escherichia coli growth and persistence in lake water.

Although enteric bacteria normally reside within the animal intestine, the ability to persist extraintestinally is an essential part of their overall lifestyle, and it might contribute to transmission between hosts. Despite this potential importance, few genetic determinants of extraintestinal growth and survival have been identified, even for the best-studied model, Escherichia coli. In this work, we thus used a genome-wide library of barcoded transposon insertions to systematically identify functional clusters of genes that are crucial for E. coli fitness in lake water. Our results revealed that inactivation of pathways involved in maintaining outer membrane integrity, nucleotide biosynthesis, and chemotaxis negatively affected E. coli growth or survival in this extraintestinal environment. In contrast, inactivation of another group of genes apparently benefited E. coli growth or persistence in filtered lake water, resulting in higher abundance of these mutants. This group included rpoS, which encodes the general stress response sigma factor, as well as genes encoding several other global transcriptional regulators and RNA chaperones, along with several poorly annotated genes. Based on this co-enrichment, we identified these gene products as novel positive regulators of RpoS activity. We further observed that, despite their enhanced growth, E. coli mutants with inactive RpoS had reduced viability in lake water, and they were not enriched in the presence of the autochthonous microbiota. This highlights the duality of the general stress response pathway for E. coli growth outside the host.

Sleeping Beauty transposon mutagenesis in mouse intestinal organoids identifies genes involved in tumor progression and metastasis.

To identify genes important for colorectal cancer (CRC) development and metastasis, we established a new metastatic mouse organoid model using Sleeping Beauty (SB) transposon mutagenesis. Intestinal organoids derived from mice carrying actively mobilizing SB transposons, an activating KrasG12D, and an inactivating ApcΔ716 allele, were transplanted to immunodeficient mice. While 66.7% of mice developed primary tumors, 7.6% also developed metastatic tumors. Analysis of SB insertion sites in tumors identified numerous candidate cancer genes (CCGs) identified previously in intestinal SB screens performed in vivo, in addition to new CCGs, such as Slit2 and Atxn1. Metastatic tumors from the same mouse were clonally related to each other and to primary tumors, as evidenced by the transposon insertion site. To provide functional validation, we knocked out Slit2, Atxn1, and Cdkn2a in mouse tumor organoids and transplanted to mice. Tumor development was promoted when these gene were knocked out, demonstrating that these are potent tumor suppressors. Cdkn2a knockout cells also metastasized to the liver in 100% of the mice, demonstrating that Cdkn2a loss confers metastatic ability. Our organoid model thus provides a new approach that can be used to understand the evolutionary forces driving CRC metastasis and a rich resource to uncover CCGs promoting CRC.

Transposon insertion in pmel17 rewired skin and muscle transcriptomes in Mozambique tilapia

The pmel17 gene plays a crucial role in melanin pigmentation. Our previous studies showed that in Mozambique tilapia, a transposon inserted into the 3′ untranslated region (3′ UTR) of pmel17 resulted in the silencing of pmel17 and led to the loss of melanin pigments (golden mutant phenotype). Additionally, the transposon insertion caused reduced growth performance and increased locomotion. In this study, to investigate the mechanisms underlying these phenotypic changes, we sequenced transcriptomes of the skin and muscle samples collected from wildtype and mutant tilapias. A total of 51 and 141 differentially expressed genes (DEGs) were identified in the skin and muscle transcriptomes, respectively. DEGs in the skin were primarily down-regulated in golden genotypes and associated with neural crest development and melanin pigmentation pathways. Besides these DEGs involved in the classic melanin pigmentation pathway of vertebrates, 14 DEGs were also observed to be related to melanogenesis. In muscle transcriptomes, there was an enrichment of GO terms associated with growth factors and cellular lipid catabolic processes. Specifically, DEGs related to growth factor binding exhibited a down-regulation, while those related to lipid metabolism showed an up-regulation in mutant genotypes. These findings agree with observed phenotypic changes. Furthermore, several DEGs associated with muscle function and mobility were up-regulated. Our study sheds light on how a single mutation in a gene can modulate multiple phenotypes by rewiring gene regulation networks. The research also provides valuable insights into the complex genetic mechanisms underlying the regulation of diverse phenotypic traits by a single gene.

CRISPR Interference-Based Functional Small RNA Genomics.

Small RNAs (sRNAs) are versatile regulators universally present in species across the prokaryotic kingdom, yet their functional characterization remains a major bottleneck. Gene inactivation through random transposon insertion has proven extremely valuable in discovering hidden gene functions. However, this approach is biased toward long genes and usually results in the underrepresentation of sRNA mutants. In contrast, CRISPR interference (CRISPRi) harnesses guide RNAs to recruit cleavage-deficient Cas nucleases to specific DNA loci. The ensuing steric hindrance inhibits RNA polymerase assembly at-or migration along-predefined genes, allowing for targeted knockdown screens without major length bias. In this chapter, we provide a detailed protocol for CRISPRi-based functional screening of bacterial sRNAs. Using the abundant microbiota species Bacteroides thetaiotaomicron as a model, we describe the design and generation of a guide library targeting the full intergenic sRNA repertoire of this organism and its application to identify sRNA knockdown-associated fitness effects. Our protocol is generic and thus suitable for the systematic assessment of sRNA-associated phenotypes in a wide range of bacterial species and experimental conditions. We expect CRISPRi-based functional genomics to boost sRNA research in understudied bacterial taxa, for instance, members of the gut microbiota.