Search

Abstract Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Abstract Multiciliated cells (MCCs) are epithelial cells that control body fluid flow and contribute to the clearance of pathogenic microbes and other particles from the airways, egg transport in oviducts, and circulation of cerebrospinal fluid in the central nervous system. Although MCCs have shared functions to control fluid flow via coordinated motility of multiple ciliary structures, they are found in multiple mammalian tissues originating from distinct germ layers and differentiate via distinct developmental pathways. To understand the similarities and differences of MCCs in multiple tissues, we investigated single-cell transcriptome data of nasal epithelial cells, bronchial tubes, fallopian tubes, and ependymal cells in the subventricular zone from humans and mice by cross-species data integration. Expression of cilia-associated genes was indistinguishable between these MCCs, although cell populations had unique properties by the species and tissue, demonstrating that they share the same final differentiation status for ciliary functions. We further analyzed the final differentiation step of MCCs from their distinctive progenitors and confirmed their convergent gene set expression for ciliogenesis at the final step. These results may provide new insight into understanding ciliogenesis during the developmental process.

The cilium is a microtubule-based organelle that plays a pivotal role in embryonic development and maintenance of physiological functions in the human body. In addition to their function as sensors that transduce diverse extracellular signals, including growth factors, fluid flow, and physical forces, cilia are intricately involved in cell cycle regulation and preservation of DNA integrity, as their formation and resorption dynamics are tightly linked to cell cycle progression. Recently, several studies have linked defects in specific ciliary proteins to the DNA damage response. However, it remains unclear whether and how primary cilia contribute to cancer development. Mebendazole (MBZ) is an anthelmintic drug with anticancer properties in some cancer cells. MBZ is continuously being tested for clinical studies, but the precise mechanism of its anticancer activities remains unknown. Here, using Xenopus laevis embryos as a model system, we discovered that MBZ significantly hinders cilia formation and induces DNA damage. Remarkably, primary cilium-bearing cancer cells exhibited heightened vulnerability to combined treatment with MBZ and conventional anticancer drugs. Our findings shed light on the specific influence of MBZ on cilia, rather than cytosolic microtubules, in triggering DNA damage, elucidating a previously unidentified mechanism underlying potential MBZ-mediated cancer therapy.

Perfluorooctanesulfonate (PFOS) is a ubiquitous environmental pollutant associated with increasing health concerns and environmental hazards. Toxicological analyses of PFOS exposure are hampered by large interspecies variations and limited studies on the mechanistic details of PFOS-induced toxicity. We investigated the effects of PFOS exposure on Xenopus laevis embryos based on the reported developmental effects in zebrafish. X. laevis was selected to further our understanding of interspecies variation in response to PFOS, and we built upon previous studies by including transcriptomics and an assessment of ciliogenic effects. Midblastula-stage X. laevis embryos were exposed to PFOS using the frog embryo teratogenesis assay Xenopus (FETAX). Results showed teratogenic effects of PFOS in a time- and dose-dependent manner. The morphological abnormalities of skeleton deformities, a small head, and a miscoiled gut were associated with changes in gene expression evidenced by whole-mount in situ hybridization and transcriptomics. The transcriptomic profile of PFOS-exposed embryos indicated the perturbation in the expression of genes associated with cell death, and downregulation in adenosine triphosphate (ATP) biosynthesis. Moreover, we observed the effects of PFOS exposure on cilia development as a reduction in the number of multiciliated cells and changes in the directionality and velocity of the cilia-driven flow. Collectively, these data broaden the molecular understanding of PFOS-induced developmental effects, whereby ciliary dysfunction and disrupted ATP synthesis are implicated as the probable modes of action of embryotoxicity. Furthermore, our findings present a new challenge to understand the links between PFOS-induced developmental toxicity and vital biological processes.

Patient-derived xenograft (PDX) models, which can retain the characteristics of original tumors in an in vivo-mimicking environment, have been developed to identify better treatment options. However, although original tumors and xenograft tissues mostly share oncogenic mutations and global gene expression patterns, their detailed mutation profiles occasionally do not overlap, indicating that selection occurs in the xenograft environment. To understand this mutational alteration in xenografts, we established 13 PDX models derived from 11 brain tumor patients and confirmed their histopathological similarity. Surprisingly, only a limited number of somatic mutations were shared between the original tumor and xenograft tissue. By analyzing deleteriously mutated genes in tumors and xenografts, we found that previously reported brain tumor-related genes were enriched in PDX samples, demonstrating that xenografts are a valuable platform for studying brain tumors. Furthermore, mutated genes involved in cilium movement, microtubule depolymerization, and histone methylation were enriched in PDX samples compared with the original tumors. Even with the limitations of the heterogeneity of clinical lesions with a heterotropic model, our study demonstrates that PDX models can provide more information in genetic analysis using samples with high heterogeneity, such as brain tumors.

The gap junction complex functions as a transport channel across the membrane. Among gap junction subunits, gap junction protein α1 (GJA1) is the most commonly expressed subunit. A recent study showed that GJA1 is necessary for the maintenance of motile cilia; however, the molecular mechanism and function of GJA1 in ciliogenesis remain unknown. Here, we examined the functions of GJA1 during ciliogenesis in human retinal pigment epithelium-1 and Xenopus laevis embryonic multiciliated-cells. GJA1 localizes to the motile ciliary axonemes or pericentriolar regions beneath the primary cilium. GJA1 depletion caused malformation of both the primary cilium and motile cilia. Further study revealed that GJA1 depletion affected several ciliary proteins such as BBS4, CP110, and Rab11 in the pericentriolar region and basal body. Interestingly, CP110 removal from the mother centriole was significantly reduced by GJA1 depletion. Importantly, Rab11, a key regulator during ciliogenesis, was immunoprecipitated with GJA1 and GJA1 knockdown caused the mislocalization of Rab11. These findings suggest that GJA1 regulates ciliogenesis by interacting with the Rab11-Rab8 ciliary trafficking pathway.

Thyroid hormone receptor-interacting protein 13 (TRIP13) participates in various regulatory steps related to the cell cycle, such as the mitotic spindle assembly checkpoint and meiotic recombination, possibly by interacting with members of the HORMA domain protein family. Recently, it was reported that TRIP13 could regulate the choice of the DNA repair pathway, i.e., homologous recombination (HR) or nonhomologous end-joining (NHEJ). However, TRIP13 is recruited to DNA damage sites within a few seconds after damage and may therefore have another function in DNA repair other than regulation of the pathway choice. Furthermore, the depletion of TRIP13 inhibited both HR and NHEJ, suggesting that TRIP13 plays other roles besides regulation of choice between HR and NHEJ. To explore the unidentified functions of TRIP13 in the DNA damage response, we investigated its genome-wide interaction partners in the context of DNA damage using quantitative proteomics with proximity labeling. We identified MRE11 as a novel interacting partner of TRIP13. TRIP13 controlled the recruitment of MDC1 to DNA damage sites by regulating the interaction between MDC1 and the MRN complex. Consistently, TRIP13 was involved in ATM signaling amplification. Our study provides new insight into the function of TRIP13 in immediate-early DNA damage sensing and ATM signaling activation.

Recent studies have suggested the existence of a blood microbiome in the healthy host. However, changes in the blood microbiome upon bloodstream infection are not known. Here, we analyzed the dynamics of the blood microbiome in a porcine model of polymicrobial bacteremia induced by fecal peritonitis. Surprisingly, we detected bacterial populations in the bloodstream even before the infection, and these populations were maintained over time. The native blood microbiome was notably taxonomically different from the fecal microbiome that was used to induce peritonitis, reflecting microbial tropism for the blood. Although the population composition after the infection was similar to that of the native blood microbiome, new bacterial strains entered the bloodstream upon peritonitis induction as clinical symptoms relevant to sepsis developed. This indicates that the bacteria detected in the blood before peritonitis induction were derived from the blood rather than a contamination. Comparison of the functional pathways enriched in the blood and fecal microbiomes revealed that communication and stress management pathways are essential for the survival of the blood microbiome.

Abstract Selective serotonin reuptake inhibitors (SSRIs) are commonly prescribed antidepressant drugs in pregnant women. Infants born following prenatal exposure to SSRIs have a higher risk for behavioral abnormalities, however, the underlying mechanisms remains unknown. Therefore, we examined the effects of prenatal fluoxetine, the most commonly prescribed SSRI, in mice. Intriguingly, chronic in utero fluoxetine treatment impaired working memory and social novelty recognition in adult males. In the medial prefrontal cortex (mPFC), a key region regulating these behaviors, we found augmented spontaneous inhibitory synaptic transmission onto the layer 5 pyramidal neurons. Fast-spiking interneurons in mPFC exhibited enhanced intrinsic excitability and serotonin-induced excitability due to upregulated serotonin (5-HT) 2A receptor (5-HT2AR) signaling. More importantly, the behavioral deficits in prenatal fluoxetine treated mice were reversed by the application of a 5-HT2AR antagonist. Taken together, our findings suggest that alterations in inhibitory neuronal modulation are responsible for the behavioral alterations following prenatal exposure to SSRIs.

Abstract Background Lysine-specific histone demethylase 5C (KDM5C) belongs to the jumonji family of demethylases and is specific for the di- and tri-demethylation of lysine 4 residues on histone 3 (H3K4 me2/3). KDM5C is expressed in the brain and skeletal muscles of humans and is associated with various biologically significant processes. KDM5C is known to be associated with X-linked mental retardation and is also involved in the development of cancer. However, the developmental significance of KDM5C has not been explored yet. In the present study, we investigated the physiological roles of KDM5C during Xenopus laevis embryonic development. Results Loss-of-function analysis using kdm5c antisense morpholino oligonucleotides indicated that kdm5c knockdown led to small-sized heads, reduced cartilage size, and malformed eyes (i.e., small-sized and deformed eyes). Molecular analyses of KDM5C functional roles using whole-mount in situ hybridization, β-galactosidase staining, and reverse transcription-polymerase chain reaction revealed that loss of kdm5c resulted in reduced expression levels of neural crest specifiers and genes involved in eye development. Furthermore, transcriptome analysis indicated the significance of KDM5C in morphogenesis and organogenesis. Conclusion Our findings indicated that KDM5C is associated with embryonic development and provided additional information regarding the complex and dynamic gene network that regulates neural crest formation and eye development. This study emphasizes the functional significance of KDM5C in Xenopus embryogenesis; however, further analysis is needed to explore the interactions of KDM5C with specific developmental genes.

While comparatively few amphibian species have been described on the North East Asian mainland in the last decades, several species have been the subject of taxonomical debates in relation to the Yellow sea. Here, we sampled Dryophytes sp. treefrogs from the Republic of Korea, the Democratic People's Republic of Korea and the People's Republic of China to clarify the status of this clade around the Yellow sea and determine the impact of sea level change on treefrogs' phylogenetic relationships. Based on genetics, call properties, adult morphology, tadpole morphology and niche modelling, we determined the segregated status species of D. suweonensis and D. immaculatus. We then proceeded to describe a new treefrog species, D. flaviventris sp. nov., from the central lowlands of the Republic of Korea. The new species is geographically segregated from D. suweonensis by the Chilgap mountain range and known to occur only in the area of Buyeo, Nonsan and Iksan in the Republic of Korea. While the Yellow sea is the principal element to the current isolation of the three clades, the paleorivers of the Yellow sea basin are likely to have been the major factor for the divergences within this clade. We recommend conducting rapid conservation assessments as these species are present on very narrow and declining ranges.

Abstract The lysine-specific histone demethylase 1A (KDM1A) was the first demethylase to challenge the concept of the irreversible nature of methylation marks. KDM1A, containing a flavin adenine dinucleotide (FAD)-dependent amine oxidase domain, demethylates histone 3 lysine 4 and histone 3 lysine 9 (H3K4me1/2 and H3K9me1/2). It has emerged as an epigenetic developmental regulator and was shown to be involved in carcinogenesis. The functional diversity of KDM1A originates from its complex structure and interactions with transcription factors, promoters, enhancers, oncoproteins, and tumor-associated genes (tumor suppressors and activators). In this review, we discuss the microenvironment of KDM1A in cancer progression that enables this protein to activate or repress target gene expression, thus making it an important epigenetic modifier that regulates the growth and differentiation potential of cells. A detailed analysis of the mechanisms underlying the interactions between KDM1A and the associated complexes will help to improve our understanding of epigenetic regulation, which may enable the discovery of more effective anticancer drugs.

Glutathione peroxidase 1 (Gpx1) and peroxiredoxin 2 (Prdx2) belong to the thiol peroxidase family of antioxidants, and have been studied for their antioxidant functions and roles in cancers. However, the physiological significance of Gpx1 and Prdx2 during vertebrate embryogenesis are lacking. Currently, we investigated the functional roles of Gpx1 and Prdx2 during vertebrate embryogenesis using Xenopus laevis as a vertebrate model. Our investigations revealed the zygotic nature of gpx1 having its localization in the eye region of developing embryos, whereas prdx2 exhibited a maternal nature and were localized in embryonic ventral blood islands. Furthermore, the gpx1-morphants exhibited malformed eyes with incompletely detached lenses. However, the depletion of prdx2 has not established its involvement with embryogenesis. A molecular analysis of gpx1-depleted embryos revealed the perturbed expression of a cryba1-lens-specific marker and also exhibited reactive oxygen species (ROS) accumulation in the eye regions of gpx1-morphants. Additionally, transcriptomics analysis of gpx1-knockout embryos demonstrated the involvement of Wnt, cadherin, and integrin signaling pathways in the development of malformed eyes. Conclusively, our findings indicate the association of gpx1 with a complex network of embryonic developmental pathways and ROS responses, but detailed investigation is a prerequisite in order to pinpoint the mechanistic details of these interactions.

Abstract Peroxiredoxin1 (Prdx1) is an antioxidant enzyme belonging to the peroxiredoxin family of proteins. Prdx1 catalyzes the reduction of H2O2 and alkyl hydroperoxide and plays an important role in different biological processes. Prdx1 also participates in various age-related diseases and cancers. In this study, we investigated the role of Prdx1 in pronephros development during embryogenesis. Prdx1 knockdown markedly inhibited proximal tubule formation in the pronephros and significantly increased the cellular levels of reactive oxygen species (ROS), which impaired primary cilia formation. Additionally, treatment with ROS (H2O2) severely disrupted proximal tubule formation, whereas Prdx1 overexpression reversed the ROS-mediated inhibition in proximal tubule formation. Epistatic analysis revealed that Prdx1 has a crucial role in retinoic acid and Wnt signaling pathways during pronephrogenesis. In conclusion, Prdx1 facilitates proximal tubule formation during pronephrogenesis by regulating ROS levels.

Determining how size is controlled is a fundamental question in biology that is poorly understood at the organismal, cellular, and subcellular levels. The Xenopus species, X. laevis and X. tropicalis differ in size at all three of these levels. Despite these differences, fertilization of X. laevis eggs with X. tropicalis sperm gives rise to viable hybrid animals that are intermediate in size. We observed that although hybrid and X. laevis embryogenesis initiates from the same sized zygote and proceeds synchronously through development, hybrid animals were smaller by the tailbud stage, and a change in the ratio of nuclear size to cell size was observed shortly after zygotic genome activation (ZGA), suggesting that differential gene expression contributes to size differences. Transcriptome analysis at the onset of ZGA identified twelve transcription factors paternally expressed in hybrids. A screen of these X. tropicalis factors by expression in X. laevis embryos revealed that Hes7 and Ventx2 significantly reduced X. laevis body length size by the tailbud stage, although nuclear to cell size scaling relationships were not affected as in the hybrid. Together, these results suggest that transcriptional regulation contributes to biological size control in Xenopus.

Glutathione peroxidase 3 (GPx3) belongs to the glutathione peroxidase family of selenoproteins and is a key antioxidant enzyme in multicellular organisms against oxidative damage. Downregulation of GPx3 affects tumor progression and metastasis and is associated with liver and heart disease. However, the physiological significance of GPx3 in vertebrate embryonic development remains poorly understood. The current study aimed to investigate the functional roles of gpx3 during embryogenesis. To this end, we determined gpx3’s spatiotemporal expression using Xenopus laevis as a model organism. Using reverse transcription polymerase chain reaction (RT-PCR), we demonstrated the zygotic nature of this gene. Interestingly, the expression of gpx3 enhanced during the tailbud stage of development, and whole mount in situ hybridization (WISH) analysis revealed gpx3 localization in prospective tail region of developing embryo. gpx3 knockdown using antisense morpholino oligonucleotides (MOs) resulted in short post-anal tails, and these malformed tails were significantly rescued by glutathione peroxidase mimic ebselen. The gene expression analysis indicated that gpx3 knockdown significantly altered the expression of genes associated with Wnt, Notch, and bone morphogenetic protein (BMP) signaling pathways involved in tailbud development. Moreover, RNA sequencing identified that gpx3 plays a role in regulation of cell death in the developing embryo. Terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) and phospho-histone 3 (PH3) staining confirmed the association of gpx3 knockdown with increased cell death and decreased cell proliferation in tail region of developing embryos, establishing the involvement of gpx3 in tailbud development by regulating the cell death. Furthermore, these findings are inter-related with increased reactive oxygen species (ROS) levels in gpx3 knockdown embryos, as measured by using a redox-sensitive fluorescent probe HyPer. Taken together, our results suggest that gpx3 plays a critical role in posterior embryonic development by regulating cell death and proliferation during vertebrate embryogenesis.

The airway epithelium in human plays a central role as the first line of defense against environmental contaminants. Most respiratory diseases such as chronic obstructive pulmonary disease (COPD), asthma, and respiratory infections, disturb normal muco-ciliary functions by stimulating the hypersecretion of mucus. Several muco-active agents have been used to treat hypersecretion symptoms in patients. Current muco-active reagents control mucus secretion by modulating either airway inflammation, cholinergic parasympathetic nerve activities or by reducing the viscosity by cleaving crosslinking in mucin and digesting DNAs in mucus. However, none of the current medication regulates mucus secretion by directly targeting airway goblet cells. The major hurdle for screening potential muco-active agents that directly affect the goblet cells, is the unavailability of in vivo model systems suitable for high-throughput screening. In this study, we developed a high-throughput in vivo model system for identifying muco-active reagents using Xenopus laevis embryos. We tested mucus secretion under various conditions and developed a screening strategy to identify potential muco-regulators. Using this novel screening technique, we identified narasin as a potential muco-regulator. Narasin treatment of developing Xenopus embryos significantly reduced mucus secretion. Furthermore, the human lung epithelial cell line, Calu-3, responded similarly to narasin treatment, validating our technique for discovering muco-active reagents.

Recently, using the frog Xenopus laevis as a model system, we showed that the transcription factor Rfx2 coordinates many genes involved in ciliogenesis and cell movement in multiciliated cells (Chung et al., 2014). To our knowledge, it was the first paper to utilize the genomic resources, including genome sequences and interim gene annotations, from the ongoing X. laevis genome project. For researchers who are interested in the application of genomics and systems biology approaches in Xenopus studies, here we provide additional details about our dataset (NCBI GEO accession number GSE50593) and describe how we analyzed RNA-seq and ChIP-seq data to identify direct targets of Rfx2.

ABSTRACT Antimicrobial-resistant bacteria pose a serious threat in the clinic. This is particularly true for opportunistic pathogens that possess high intrinsic resistance. Though many studies have focused on understanding the acquisition of bacterial resistance upon exposure to antimicrobials, the mechanisms controlling intrinsic resistance are not well understood. In this study, we subjected the model opportunistic superbug Pseudomonas aeruginosa to 14 antimicrobials under highly controlled conditions and assessed its response using expression- and fitness-based genomic approaches. Our results reveal that gene expression changes and mutant fitness in response to sub-MIC antimicrobials do not correlate on a genomewide scale, indicating that gene expression is not a good predictor of fitness determinants. In general, fewer fitness determinants were identified for antiseptics and disinfectants than for antibiotics. Analysis of gene expression and fitness data together allowed the prediction of antagonistic interactions between antimicrobials and insight into the molecular mechanisms controlling these interactions. IMPORTANCE Infections involving multidrug-resistant pathogens are difficult to treat because the therapeutic options are limited. These infections impose a significant financial burden on infected patients and on health care systems. Despite years of antimicrobial resistance research, we lack a comprehensive understanding of the intrinsic mechanisms controlling antimicrobial resistance. This work uses two fine-scale genomic approaches to identify genetic loci important for antimicrobial resistance of the opportunistic pathogen Pseudomonas aeruginosa. Our results reveal that antibiotics have more resistance determinants than antiseptics/disinfectants and that gene expression upon exposure to antimicrobials is not a good predictor of these resistance determinants. In addition, we show that when used together, genomewide gene expression and fitness profiling can provide mechanistic insights into multidrug resistance mechanisms.

The mechanisms linking systems-level programs of gene expression to discrete cell biological processes in vivo remain poorly understood. In this study, we have defined such a program for multi-ciliated epithelial cells (MCCs), a cell type critical for proper development and homeostasis of the airway, brain and reproductive tracts. Starting from genomic analysis of the cilia-associated transcription factor Rfx2, we used bioinformatics and in vivo cell biological approaches to gain insights into the molecular basis of cilia assembly and function. Moreover, we discovered a previously un-recognized role for an Rfx factor in cell movement, finding that Rfx2 cell-autonomously controls apical surface expansion in nascent MCCs. Thus, Rfx2 coordinates multiple, distinct gene expression programs in MCCs, regulating genes that control cell movement, ciliogenesis, and cilia function. As such, the work serves as a paradigm for understanding genomic control of cell biological processes that span from early cell morphogenetic events to terminally differentiated cellular functions.

When a duplicate gene has no apparent loss-of-function phenotype, it is commonly considered that the phenotype has been masked as a result of functional redundancy with the remaining paralog. This is supported by indirect evidence showing that multi-copy genes show loss-of-function phenotypes less often than single-copy genes and by direct tests of phenotype masking using select gene sets. Here we take a systematic genome-wide RNA interference approach to assess phenotype masking in paralog pairs in the Caenorhabditis elegans genome. Remarkably, in contrast to expectations, we find that phenotype masking makes only a minor contribution to the low knockdown phenotype rate for duplicate genes. Instead, we find that non-essential genes are highly over-represented among duplicates, leading to a low observed loss-of-function phenotype rate. We further find that duplicate pairs derived from essential and non-essential genes have contrasting evolutionary dynamics: whereas non-essential genes are both more often successfully duplicated (fixed) and lost, essential genes are less often duplicated but upon successful duplication are maintained over longer periods. We expect the fundamental evolutionary duplication dynamics presented here to be broadly applicable.

The gram-negative opportunistic pathogen Pseudomonas aeruginosa is the primary cause of chronic respiratory infections in individuals with the heritable disease cystic fibrosis (CF). These infections can last for decades, during which time P. aeruginosa has been proposed to acquire beneficial traits via adaptive evolution. Because CF lacks an animal model that can acquire chronic P. aeruginosa infections, identifying genes important for long-term in vivo fitness remains difficult. However, since clonal, chronological samples can be obtained from chronically infected individuals, traits undergoing adaptive evolution can be identified. Recently we identified 24 P. aeruginosa gene expression traits undergoing parallel evolution in vivo in multiple individuals, suggesting they are beneficial to the bacterium. The goal of this study was to determine if these genes impact P. aeruginosa phenotypes important for survival in the CF lung. By using a gain-of-function genetic screen, we found that 4 genes and 2 operons undergoing parallel evolution in vivo promote P. aeruginosa biofilm formation. These genes/operons promote biofilm formation by increasing levels of the non-alginate exopolysaccharide Psl. One of these genes, phaF, enhances Psl production via a post-transcriptional mechanism, while the other 5 genes/operons do not act on either psl transcription or translation. Together, these data demonstrate that P. aeruginosa has evolved at least two pathways to over-produce a non-alginate exopolysaccharide during long-term colonization of the CF lung. More broadly, this approach allowed us to attribute a biological significance to genes with unknown function, demonstrating the power of using evolution as a guide for targeted genetic studies.

ABSTRACT The Gram-negative bacterium Pseudomonas aeruginosa is a common cause of chronic airway infections in individuals with the heritable disease cystic fibrosis (CF). After prolonged colonization of the CF lung, P. aeruginosa becomes highly resistant to host clearance and antibiotic treatment; therefore, understanding how this bacterium evolves during chronic infection is important for identifying beneficial adaptations that could be targeted therapeutically. To identify potential adaptive traits of P. aeruginosa during chronic infection, we carried out global transcriptomic profiling of chronological clonal isolates obtained from 3 individuals with CF. Isolates were collected sequentially over periods ranging from 3 months to 8 years, representing up to 39,000 in vivo generations. We identified 24 genes that were commonly regulated by all 3 P. aeruginosa lineages, including several genes encoding traits previously shown to be important for in vivo growth. Our results reveal that parallel evolution occurs in the CF lung and that at least a proportion of the traits identified are beneficial for P. aeruginosa chronic colonization of the CF lung. IMPORTANCE Deadly diseases like AIDS, malaria, and tuberculosis are the result of long-term chronic infections. Pathogens that cause chronic infections adapt to the host environment, avoiding the immune response and resisting antimicrobial agents. Studies of pathogen adaptation are therefore important for understanding how the efficacy of current therapeutics may change upon prolonged infection. One notorious chronic pathogen is Pseudomonas aeruginosa, a bacterium that causes long-term infections in individuals with the heritable disease cystic fibrosis (CF). We used gene expression profiles to identify 24 genes that commonly changed expression over time in 3 P. aeruginosa lineages, indicating that these changes occur in parallel in the lungs of individuals with CF. Several of these genes have previously been shown to encode traits critical for in vivo-relevant processes, suggesting that they are likely beneficial adaptations important for chronic colonization of the CF lung.