Chromatin modifier HUSH co-operates with RNA decay factor NEXT to restrict transposable element expression.

Transposable elements (TEs) are widespread genetic parasites known to be kept under tight transcriptional control. Here, we describe a functional connection between the mouse-orthologous "nuclear exosome targeting" (NEXT) and "human silencing hub" (HUSH) complexes, involved in nuclear RNA decay and the epigenetic silencing of TEs, respectively. Knocking out the NEXT component ZCCHC8 in embryonic stem cells results in elevated TE RNA levels. We identify a physical interaction between ZCCHC8 and the MPP8 protein of HUSH and establish that HUSH recruits NEXT to chromatin at MPP8-bound TE loci. However, while NEXT and HUSH both dampen TE RNA expression, their activities predominantly affect shorter non-polyadenylated and full-length polyadenylated transcripts, respectively. Indeed, our data suggest that the repressive action of HUSH promotes a condition favoring NEXT RNA decay activity. In this way, transcriptional and post-transcriptional machineries synergize to suppress the genotoxic potential of TE RNAs.

Determinants of enhancer and promoter activities of regulatory elements.

The proper activities of enhancers and gene promoters are essential for coordinated transcription within a cell. Although diverse methodologies have been developed to identify enhancers and promoters, most have tacitly assumed that these elements are distinct. However, studies have unexpectedly shown that regulatory elements may have both enhancer and promoter functions. Here we review these results, focusing on the factors that determine the promoter and/or enhancer activity of regulatory elements. We discuss emerging models that define regulatory elements by accessible DNA and their non-mutually-exclusive abilities to drive transcription initiation (promoter activity) and/or to enhance transcription at other such regions (enhancer activity).

JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profiles.

JASPAR (https://jaspar.elixir.no/) is a widely-used open-access database presenting manually curated high-quality and non-redundant DNA-binding profiles for transcription factors (TFs) across taxa. In this 10th release and 20th-anniversary update, the CORE collection has expanded with 329 new profiles. We updated three existing profiles and provided orthogonal support for 72 profiles from the previous release's UNVALIDATED collection. Altogether, the JASPAR 2024 update provides a 20% increase in CORE profiles from the previous release. A trimming algorithm enhanced profiles by removing low information content flanking base pairs, which were likely uninformative (within the capacity of the PFM models) for TFBS predictions and modelling TF-DNA interactions. This release includes enhanced metadata, featuring a refined classification for plant TFs' structural DNA-binding domains. The new JASPAR collections prompt updates to the genomic tracks of predicted TF binding sites (TFBSs) in 8 organisms, with human and mouse tracks available as native tracks in the UCSC Genome browser. All data are available through the JASPAR web interface and programmatically through its API and the updated Bioconductor and pyJASPAR packages. Finally, a new TFBS extraction tool enables users to retrieve predicted JASPAR TFBSs intersecting their genomic regions of interest.

A data-driven genome annotation approach for cassava.

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave.

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

Identification of a Nuclear Exosome Decay Pathway for Processed Transcripts.

The RNA exosome is fundamental for the degradation of RNA in eukaryotic nuclei. Substrate targeting is facilitated by its co-factor Mtr4p/hMTR4, which links to RNA-binding protein adaptors. One example is the trimeric human nuclear exosome targeting (NEXT) complex, which is composed of hMTR4, the Zn-finger protein ZCCHC8, and the RNA-binding factor RBM7. NEXT primarily targets early and unprocessed transcripts, which demands a rationale for how the nuclear exosome recognizes processed RNAs. Here, we describe the poly(A) tail exosome targeting (PAXT) connection, which comprises the ZFC3H1 Zn-knuckle protein as a central link between hMTR4 and the nuclear poly(A)-binding protein PABPN1. Individual depletion of ZFC3H1 and PABPN1 results in the accumulation of common transcripts that are generally both longer and more extensively polyadenylated than NEXT substrates. Importantly, ZFC3H1/PABPN1 and ZCCHC8/RBM7 contact hMTR4 in a mutually exclusive manner, revealing that the exosome targets nuclear transcripts of different maturation status by substituting its hMTR4-associating adaptors.

JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net/) is an open-access database containing manually curated, non-redundant transcription factor (TF) binding profiles for TFs across six taxonomic groups. In this 9th release, we expanded the CORE collection with 341 new profiles (148 for plants, 101 for vertebrates, 85 for urochordates, and 7 for insects), which corresponds to a 19% expansion over the previous release. We added 298 new profiles to the Unvalidated collection when no orthogonal evidence was found in the literature. All the profiles were clustered to provide familial binding profiles for each taxonomic group. Moreover, we revised the structural classification of DNA binding domains to consider plant-specific TFs. This release introduces word clouds to represent the scientific knowledge associated with each TF. We updated the genome tracks of TFBSs predicted with JASPAR profiles in eight organisms; the human and mouse TFBS predictions can be visualized as native tracks in the UCSC Genome Browser. Finally, we provide a new tool to perform JASPAR TFBS enrichment analysis in user-provided genomic regions. All the data is accessible through the JASPAR website, its associated RESTful API, the R/Bioconductor data package, and a new Python package, pyJASPAR, that facilitates serverless access to the data.

Crosstalk between tumor acidosis, p53 and extracellular matrix regulates pancreatic cancer aggressiveness.

Pancreatic ductal adenocarcinoma (PDAC) is an extremely aggressive malignancy with minimal treatment options and a global rise in prevalence. PDAC is characterized by frequent driver mutations including KRAS and TP53 (p53), and a dense, acidic tumor microenvironment (TME). The relation between genotype and TME in PDAC development is unknown. Strikingly, when wild type (WT) Panc02 PDAC cells were adapted to growth in an acidic TME and returned to normal pH to mimic invasive cells escaping acidic regions, they displayed a strong increase of aggressive traits such as increased growth in 3-dimensional (3D) culture, adhesion-independent colony formation and invasive outgrowth. This pattern of acidosis-induced aggressiveness was observed in 3D spheroid culture as well as upon organotypic growth in matrigel, collagen-I and combination thereof, mimicking early and later stages of PDAC development. Acid-adaptation-induced gain of cancerous traits was further increased by p53 knockout (KO), but only in specific extracellular matrix (ECM) compositions. Akt- and Transforming growth factor-ß (TGFß) signaling, as well as expression of the Na+ /H+ exchanger NHE1, were increased by acid adaptation. Whereas Akt inhibition decreased spheroid growth regardless of treatment and genotype, stimulation with TGFßI increased growth of WT control spheroids, and inhibition of TGFß signaling tended to limit growth under acidic conditions only. Our results indicate that a complex crosstalk between tumor acidosis, ECM composition and genotype contributes to PDAC development. The findings may guide future strategies for acidosis-targeted therapies.

Chronic acidosis rewires cancer cell metabolism through PPARα signaling.

The mechanisms linking tumor microenvironment acidosis to disease progression are not understood. Here, we used mammary, pancreatic, and colon cancer cells to show that adaptation to growth at an extracellular pH (pHe ) mimicking acidic tumor niches is associated with upregulated net acid extrusion capacity and elevated intracellular pH at physiological pHe , but not at acidic pHe . Using metabolic profiling, shotgun lipidomics, imaging and biochemical analyses, we show that the acid adaptation-induced phenotype is characterized by a shift toward oxidative metabolism, increased lipid droplet-, triacylglycerol-, peroxisome content and mitochondrial hyperfusion. Peroxisome proliferator-activated receptor-α (PPARA, PPARα) expression and activity are upregulated, at least in part by increased fatty acid uptake. PPARα upregulates genes driving increased mitochondrial and peroxisomal mass and ß-oxidation capacity, including mitochondrial lipid import proteins CPT1A, CPT2 and SLC25A20, electron transport chain components, peroxisomal proteins PEX11A and ACOX1, and thioredoxin-interacting protein (TXNIP), a negative regulator of glycolysis. This endows acid-adapted cancer cells with increased capacity for utilizing fatty acids for metabolic needs, while limiting glycolysis. As a consequence, the acid-adapted cells exhibit increased sensitivity to PPARα inhibition. We conclude that PPARα is a key upstream regulator of metabolic changes favoring cancer cell survival in acidic tumor niches.

Characterization of Arabidopsis thaliana Promoter Bidirectionality and Antisense RNAs by Inactivation of Nuclear RNA Decay Pathways.

In animals, RNA polymerase II initiates transcription bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are degraded by the nuclear exosome. Previous studies based on nascent RNA approaches concluded that Arabidopsis (Arabidopsis thaliana) does not produce PROMPTs. Here, we used steady-state RNA sequencing in mutants defective in nuclear RNA decay including the exosome to reassess the existence of Arabidopsis PROMPTs. While they are rare, we identified ∼100 cases of exosome-sensitive PROMPTs in Arabidopsis. Such PROMPTs are sources of small interfering RNAs in exosome-deficient mutants, perhaps explaining why plants have evolved mechanisms to suppress PROMPTs. In addition, we found ∼200 long, unspliced and exosome-sensitive antisense RNAs that arise from transcription start sites within parts of the genome encoding 3'-untranslated regions on the sense strand. The previously characterized noncoding RNA that regulates expression of the key seed dormancy regulator, DELAY OF GERMINATION1, is a typical representative of this class of RNAs. Transcription factor genes are overrepresented among loci with exosome-sensitive antisense RNAs, suggesting a potential for widespread control of gene expression via this class of noncoding RNAs. Lastly, we assess the use of alternative promoters in Arabidopsis and compare the accuracy of existing TSS annotations.

Lactate receptor GPR81 drives breast cancer growth and invasiveness through regulation of ECM properties and Notch ligand DLL4.

BACKGROUND: The lactate receptor GPR81 contributes to cancer development through unclear mechanisms. Here, we investigate the roles of GPR81 in three-dimensional (3D) and in vivo growth of breast cancer cells and study the molecular mechanisms involved. METHODS: GPR81 was stably knocked down (KD) in MCF-7 human breast cancer cells which were subjected to RNA-seq analysis, 3D growth, in situ- and immunofluorescence analyses, and cell viability- and motility assays, combined with KD of key GPR81-regulated genes. Key findings were additionally studied in other breast cancer cell lines and in mammary epithelial cells. RESULTS: GPR81 was upregulated in multiple human cancer types and further upregulated by extracellular lactate and 3D growth in breast cancer spheroids. GPR81 KD increased spheroid necrosis, reduced invasion and in vivo tumor growth, and altered expression of genes related to GO/KEGG terms extracellular matrix, cell adhesion, and Notch signaling. Single cell in situ analysis of MCF-7 cells revealed that several GPR81-regulated genes were upregulated in the same cell clusters. Notch signaling, particularly the Notch ligand Delta-like-4 (DLL4), was strikingly downregulated upon GPR81 KD, and DLL4 KD elicited spheroid necrosis and inhibited invasion in a manner similar to GPR81 KD. CONCLUSIONS: GPR81 supports breast cancer aggressiveness, and in MCF-7 cells, this occurs at least in part via DLL4. Our findings reveal a new GPR81-driven mechanism in breast cancer and substantiate GPR81 as a promising treatment target.

Organismal benefits of transcription speed control at gene boundaries.

RNA polymerase II (RNAPII) transcription is crucial for gene expression. RNAPII density peaks at gene boundaries, associating these key regions for gene expression control with limited RNAPII movement. The connections between RNAPII transcription speed and gene regulation in multicellular organisms are poorly understood. Here, we directly modulate RNAPII transcription speed by point mutations in the second largest subunit of RNAPII in Arabidopsis thaliana. A RNAPII mutation predicted to decelerate transcription is inviable, while accelerating RNAPII transcription confers phenotypes resembling auto-immunity. Nascent transcription profiling revealed that RNAPII complexes with accelerated transcription clear stalling sites at both gene ends, resulting in read-through transcription. The accelerated transcription mutant NRPB2-Y732F exhibits increased association with 5' splice site (5'SS) intermediates and enhanced splicing efficiency. Our findings highlight potential advantages of RNAPII stalling through local reduction in transcription speed to optimize gene expression for the development of multicellular organisms.

The RNA exosome shapes the expression of key protein-coding genes.

The ribonucleolytic exosome complex is central for nuclear RNA degradation, primarily targeting non-coding RNAs. Still, the nuclear exosome could have protein-coding (pc) gene-specific regulatory activities. By depleting an exosome core component, or components of exosome adaptor complexes, we identify ∼2900 transcription start sites (TSSs) from within pc genes that produce exosome-sensitive transcripts. At least 1000 of these overlap with annotated mRNA TSSs and a considerable portion of their transcripts share the annotated mRNA 3' end. We identify two types of pc-genes, both employing a single, annotated TSS across cells, but the first type primarily produces full-length, exosome-sensitive transcripts, whereas the second primarily produces prematurely terminated transcripts. Genes within the former type often belong to immediate early response transcription factors, while genes within the latter are likely transcribed as a consequence of their proximity to upstream TSSs on the opposite strand. Conversely, when genes have multiple active TSSs, alternative TSSs that produce exosome-sensitive transcripts typically do not contribute substantially to overall gene expression, and most such transcripts are prematurely terminated. Our results display a complex landscape of sense transcription within pc-genes and imply a direct role for nuclear RNA turnover in the regulation of a subset of pc-genes.

JASPAR 2020: update of the open-access database of transcription factor binding profiles.

JASPAR (http://jaspar.genereg.net) is an open-access database of curated, non-redundant transcription factor (TF)-binding profiles stored as position frequency matrices (PFMs) for TFs across multiple species in six taxonomic groups. In this 8th release of JASPAR, the CORE collection has been expanded with 245 new PFMs (169 for vertebrates, 42 for plants, 17 for nematodes, 10 for insects, and 7 for fungi), and 156 PFMs were updated (125 for vertebrates, 28 for plants and 3 for insects). These new profiles represent an 18% expansion compared to the previous release. JASPAR 2020 comes with a novel collection of unvalidated TF-binding profiles for which our curators did not find orthogonal supporting evidence in the literature. This collection has a dedicated web form to engage the community in the curation of unvalidated TF-binding profiles. Moreover, we created a Q&A forum to ease the communication between the user community and JASPAR curators. Finally, we updated the genomic tracks, inference tool, and TF-binding profile similarity clusters. All the data is available through the JASPAR website, its associated RESTful API, and through the JASPAR2020 R/Bioconductor package.

Human nonsense-mediated RNA decay initiates widely by endonucleolysis and targets snoRNA host genes.

Eukaryotic RNAs with premature termination codons (PTCs) are eliminated by nonsense-mediated decay (NMD). While human nonsense RNA degradation can be initiated either by an endonucleolytic cleavage event near the PTC or through decapping, the individual contribution of these activities on endogenous substrates has remained unresolved. Here we used concurrent transcriptome-wide identification of NMD substrates and their 5'-3' decay intermediates to establish that SMG6-catalyzed endonucleolysis widely initiates the degradation of human nonsense RNAs, whereas decapping is used to a lesser extent. We also show that a large proportion of genes hosting snoRNAs in their introns produce considerable amounts of NMD-sensitive splice variants, indicating that these RNAs are merely by-products of a primary snoRNA production process. Additionally, transcripts from genes encoding multiple snoRNAs often yield alternative transcript isoforms that allow for differential expression of individual coencoded snoRNAs. Based on our findings, we hypothesize that snoRNA host genes need to be highly transcribed to accommodate high levels of snoRNA production and that the expression of individual snoRNAs and their cognate spliced RNA can be uncoupled via alternative splicing and NMD.

TrancriptomeReconstructoR: data-driven annotation of complex transcriptomes.

BACKGROUND: The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data. RESULTS: We developed TranscriptomeReconstructoR, an R package which implements a pipeline for automated transcriptome annotation. It relies on integrating features from independent and complementary datasets: (i) full-length RNA-seq for detection of splicing patterns and (ii) high-throughput 5' and 3' tag sequencing data for accurate definition of gene borders. The pipeline can also take a nascent RNA-seq dataset to supplement the called gene model with transient transcripts. We reconstructed de novo the transcriptional landscape of wild type Arabidopsis thaliana seedlings and Saccharomyces cerevisiae cells as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the most commonly used community gene models, TAIR10 and Araport11 for A.thaliana and SacCer3 for S.cerevisiae. In particular, we identify multiple transient transcripts missing from the existing annotations. Our new annotations promise to improve the quality of A.thaliana and S.cerevisiae genome research. CONCLUSIONS: Our proof-of-concept data suggest a cost-efficient strategy for rapid and accurate annotation of complex eukaryotic transcriptomes. We combine the choice of library preparation methods and sequencing platforms with the dedicated computational pipeline implemented in the TranscriptomeReconstructoR package. The pipeline only requires prior knowledge on the reference genomic DNA sequence, but not the transcriptome. The package seamlessly integrates with Bioconductor packages for downstream analysis.

Genomic Stability and Phenotypic Characteristics of Industrially Produced Lacticaseibacillus rhamnosus GG in a Yogurt Matrix.

Lacticaseibacillus rhamnosus GG is a widely marketed probiotic with well-documented probiotic properties. Previously, deletion of the mucus-adhesive spaCBA-srtC1 genes in dairy isolates was reported. In this study, we examined the genome preservation of industrially produced L. rhamnosus GG (DSM 33156) cofermented in yogurts. In total, DNA of 66 samples, including 60 isolates, was sequenced. Population samples and 59 isolates exhibited an intact genome. One isolate exhibited loss of spaCBA-srtC1. In addition, we examined phenotypes related to the probiotic properties of L. rhamnosus GG either from frozen pellets or cofermented in yogurt. L. rhamnosus GG from frozen pellets induced a response in intestinal barrier function in vitro, in contrast to frozen pellets of the starter culture. Yogurt matrix, containing only the starter culture, induced a response, but cofermentation with L. rhamnosus GG induced a higher response. Conversely, only the starter culture stimulated cytokine secretion in dendritic cells, and it was observed that the addition of L. rhamnosus GG to the starter culture reduced the response. We conclude that the L. rhamnosus GG genome is preserved in yogurt and that common in vitro probiotic effects of L. rhamnosus GG are observed when examined in the yogurt matrix. IMPORTANCELacticaseibacillus rhamnosus GG is a well-documented probiotic strain recognized for its high acid and bile tolerance and properties of adhesion to enterocytes and mucus. The strain exhibits SpaCBA pili, which have been demonstrated to play an important role in adhesion and therefore are relevant for persistence in the gastrointestinal tract. Recently we demonstrated that the genome and phenotypes of L. rhamnosus GG are preserved throughout an industrial production pipeline. However, as gene deletions in L. rhamnosus GG were previously reported for isolates from dairy products, a key question on the genomic stability of L. rhamnosus GG in a yogurt matrix remained. The aim of this study was to analyze genome stability and phenotypic characteristics of L. rhamnosus GG in yogurt. We found that the genome of L. rhamnosus GG is well conserved when the organism is cofermented in yogurt. Some phenotypic characteristics are consistent in all product matrixes, while other characteristics are modulated.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response.

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50-75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

IsoformSwitchAnalyzeR: analysis of changes in genome-wide patterns of alternative splicing and its functional consequences.

SUMMARY: Alternative splicing is an important mechanism involved in health and disease. Recent work highlights the importance of investigating genome-wide changes in splicing patterns and the subsequent functional consequences. Current computational methods only support such analysis on a gene-by-gene basis. Therefore, we extended IsoformSwitchAnalyzeR R library to enable analysis of genome-wide changes in specific types of alternative splicing and predicted functional consequences of the resulting isoform switches. As a case study, we analyzed RNA-seq data from The Cancer Genome Atlas and found systematic changes in alternative splicing and the consequences of the associated isoform switches. AVAILABILITY AND IMPLEMENTATION: Windows, Linux and Mac OS: http://bioconductor.org/packages/IsoformSwitchAnalyzeR. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Lactobacillus rhamnosus GG Genomic and Phenotypic Stability in an Industrial Production Process.

Lactobacillus rhamnosus GG is one of the most widely marketed and studied probiotic strains. In L. rhamnosus GG, the spaCBA-srtC1 gene cluster encodes pili, which are important for some of the probiotic properties of the strain. A previous study showed that the DNA sequence of the spaCBA-srtC1 gene cluster was not present in some L. rhamnosus GG variants isolated from liquid dairy products. To examine the stability of the L. rhamnosus GG genome in an industrial production process, we sequenced the genome of samples of L. rhamnosus GG (DSM 33156) collected at specific steps of the industrial production process, including the culture collection stock, intermediate fermentations, and final freeze-dried products. We found that the L. rhamnosus GG genome sequence was unchanged throughout the production process. Consequently, the spaCBA-srtC1 gene locus was intact and fully conserved in all 31 samples examined. In addition, different production batches of L. rhamnosus GG exhibited consistent phenotypes, including the presence of pili in final freeze-dried products, and consistent characteristics in in vitro assays of probiotic properties. Our data show that L. rhamnosus GG is highly stable in this industrial production process.IMPORTANCELactobacillus rhamnosus GG is one of the best-studied probiotic strains. One of the well-characterized features of the strain is the pili encoded by the spaCBA-srtC1 gene cluster. These pili are involved in persistence in the gastrointestinal tract and are important for the probiotic properties of L. rhamnosus GG. Previous studies demonstrated that the L. rhamnosus GG genome can be unstable under certain conditions and can lose the spaCBA-srtC1 gene cluster. Since in vitro studies have shown that the loss of the spaCBA-srtC1 gene cluster decreases certain L. rhamnosus GG probiotic properties, we assessed both the genomic stability and phenotypic properties of L. rhamnosus GG throughout an industrial production process. We found that neither genomic nor phenotypic changes occurred in the samples. Therefore, we demonstrate that L. rhamnosus GG retains the spaCBA-srtC1 cluster and exhibits excellent genomic and phenotypic stability in the specific industrial process examined here.