Search

Motile cilia of multiciliated epithelial cells have important roles in animal development and cell homeostasis. Although several studies have identified and reported proteins localized in this complex organelle and the related immotile primary cilia from various cell types, it is still challenging to isolate high quantities of ciliary proteins for proteomic analysis. In this study, African clawed frog (Xenopus laevis) embryos, which have many multiciliated cells in the epidermis, were treated with a simple ionic buffer to identify 1009 proteins conserved across vertebrates; these proteins were putatively localized in motile cilia. Using two ciliary proteome databases, we confirmed that previously validated cilia-associated proteins are highly enriched in our ciliary proteome. Proteins localized at the transition zone and Ellis-van Creveld zone, which are distinct regions at the base of cilia, near the junction with the apical cell surface, were isolated using our method. Among the newly identified ciliary proteins, we report that KRT17 may have an unrecognized function in motile cilia. Hence, the method developed in this study would be useful for understanding the ciliary proteome.

Exposure to particulate matter (PM) in ambient air is known to increase the risk of cardiovascular disorders and mortality. The cytotoxicity of PM is mainly due to the abnormal increase of reactive oxygen species (ROS), which damage cellular components such as DNA, RNA, and proteins. The correlation between PM exposure and human disorders, including mortality, is based on long-term exposure. In this study we have investigated acute responses of mucus-secreting goblet cells upon exposure to PM derived from a heavy diesel engine. To this end, we employed the mucociliary epithelium of amphibian embryos and human Calu-3 cells to examine PM mucotoxicity. Our data suggest that acute exposure to PM significantly impairs mucus secretion and results in the accumulation of mucus vesicles in the cytoplasm of goblet cells. RNA-seq analysis revealed that acute responses to PM exposure significantly altered gene expression patterns; however, known regulators of mucus production and the secretory pathway were not significantly altered. Interestingly, pretreatment with α-tocopherol nearly recovered the hyposecretion of mucus from both amphibian and human goblet cells. We believe this study demonstrates the mucotoxicity of PM and the protective function of α-tocopherol on mucotoxicity caused by acute PM exposure from heavy diesel engines.

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. Electronic supplementary material The online version of this article (10.1186/s13041-019-0452-5) contains supplementary material, which is available to authorized users.

Figure S1. Prenatal FLX treatment does not induce other behavioral deficits. Figure S2. The prenatal FLX treatment does not change the morphology or spine density of L5 PrL neurons. Figure S3. Passive membrane properties of FS interneurons in the PrL of SAL- and FLX- treated mice. Figure S4. Expression of the mRNAs encoding serotonergic receptors and transporters and effects of the 5-HT treatment on IPSC frequency and amplitude. Figure S5. Effect of 5-HT1AR antagonists on FS interneurons in the PrL of SAL- and FLX-treated mice. Figure S6. Effects of the 5-HT treatment on pyramidal neurons in the PrL of SAL- and FLX-treated mice. Table S1. Summary of studies investigating the effects of perinatal serotonin reuptake inhibitors (SSRI) on adult male mice. Table S2. Intrinsic properties of fast-spiking interneurons of SAL and FLX treated mice before and after 5HT-treatment. Table S3. Statistical analysis conducted for each behavioral test. (PDF 1403 kb)

In design optimization, as the number of input variables increases, the convergence rate of optimization tends to decrease, and the number of function calls and design change costs tend to increase. Neighborhood component feature selection (NCFS) was adopted to select significant input variables. However, the parameter determination process of the NCFS incurs a high computational cost and weakens robustness. Therefore, this study proposes a normalized NCFS (nNCFS) by normalizing scales between mean loss and regularization terms via the initial dataset Additionally, in the case of a multi-response system, complex decision-making processes that involve the allocation of weights for multiple responses are required. It is possible to allocate weights by using conventional methods such as the analytic hierarchy process and entropy methods. However, the analytic hierarchy process method is highly influenced by the designer’s subjectivity, and the entropy method is unable to consider a design optimization problem. Accordingly, the feasible-improved weight allocation (FIWA) method is now proposed by considering a design optimization problem objectively. Comparing the NCFS with the nNCFS through mathematical examples, we found that the nNCFS significantly improved the computational cost and robustness. Moreover, the FIWA method selected significant input variables that yielded feasible and improved designs. Then, the nNCFS and the FIWA methods were applied to the design of the body-in-white of a vehicle. The significance of input variables was analyzed using the nNCFS, and feasible and improved designs were obtained on the basis of the significant input variables selected using the FIWA method.

Mammalian target of rapamycin (mTOR) signaling is a core pathway in cellular metabolism, and control of the mTOR pathway by rapamycin shows potential for the treatment of metabolic diseases. In this study, we employed a new proximity biotin-labeling method using promiscuous biotin ligase (pBirA) to identify unknown elements in the rapamycin-induced interactome on the FK506-rapamycin binding (FRB) domain in living cells. FKBP25 showed the strongest biotin labeling by FRB–pBirA in the presence of rapamycin. Immunoprecipitation and immunofluorescence experiments confirmed that endogenous FKBP25 has a rapamycin-induced physical interaction with the FRB domain. Furthermore, the crystal structure of the ternary complex of FRB–rapamycin–FKBP25 was determined at 1.67-Å resolution. In this crystal structure we found that the conformational changes of FRB generate a hole where there is a methionine-rich space, and covalent metalloid coordination was observed at C2085 of FRB located at the bottom of the hole. Our results imply that FKBP25 might have a unique physiological role related to metallomics in mTOR signaling., A proximity biotin-labeling method revealed that FKBP25 is a strong rapamycin-induced interaction partner of FK506-rapamycin binding (FRB). This interaction was supported by the crystal structure of the ternary complex of FKBP25−rapamycin−FRB.

Developing and mature chondrocytes constantly interact with and remodel the surrounding extracellular matrix (ECM). Recent research indicates that integrin-ECM interaction is differentially regulated during cartilage formation (chondrogenesis). Integrin signaling is also a key source of the catabolic reactions responsible for joint destruction in both rheumatoid arthritis and osteoarthritis. However, we do not understand how chondrocytes dynamically regulate integrin signaling in such an ECM-rich environment. Here, we found that developing chondrocytes express integrin-β-like 1 (Itgbl1) at specific stages, inhibiting integrin signaling and promoting chondrogenesis. Unlike cytosolic integrin inhibitors, ITGBL1 is secreted and physically interacts with integrins to down-regulate activity. We observed that Itgbl1 expression was strongly reduced in the damaged articular cartilage of patients with osteoarthritis (OA). Ectopic expression of Itgbl1 protected joint cartilage against OA development in the destabilization of the medial meniscus-induced OA mouse model. Our results reveal ITGBL1 signaling as an underlying mechanism of protection against destructive cartilage disorders and suggest the potential therapeutic utility of targeting ITGBL1 to modulate integrin signaling in human disease.

Peroxiredoxin5 (Prdx5), a thioredoxin peroxidase, is an antioxidant enzyme that is widely studied for its antioxidant properties and protective roles in neurological and cardiovascular disorders. This study is aimed at investigating the functional significance of Prdx5 in mitochondria and at analyzing its roles in ciliogenesis during the process of vertebrate development.We found that several Prdx genes were strongly expressed in multiciliated cells in developing Xenopus embryos, and their peroxidatic functions were crucial for normal cilia development. Depletion of Prdx5 increased levels of cellular reactive oxygen species (ROS), consequently leading to mitochondrial dysfunction and abnormal cilia formation. Proteomic and transcriptomic approaches revealed that excessive ROS accumulation on Prdx5 depletion subsequently reduced the expression level of pyruvate kinase (PK), a key metabolic enzyme in energy production. We further confirmed that the promotor activity of PK was significantly reduced on Prdx5 depletion and that the reduction in PK expression and its promoter activity led to ciliary defects observed in Prdx5-depleted cells.Our data revealed the novel relationship between ROS and Prdx5 and the consequent effects of this interaction on vertebrate ciliogenesis. The normal process of ciliogenesis is interrupted by the Prdx5 depletion, resulting in excessive ROS levels and suggesting cilia as vulnerable targets of ROS.Prdx5 plays protective roles in mitochondria and is critical for normal cilia development by regulating the levels of ROS. The loss of Prdx5 is associated with excessive production of ROS, resulting in mitochondrial dysfunction and aberrant ciliogenesis.

Additional file 2: Fig. S2. Transcriptome analysis revealed that KDM5C is essential for Xenopus embryonic development. We performed RNA-seq and analyzed groups of genes that are essential for several biological processes. The downregulation of specific genes by kdm5c knockdown indicated that KDM5C plays significant roles in organ development and structure morphogenesis.

This is the open peer reviewers comments and recommendations regarding the submitted GigaScience article and/or dataset.

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.

Summary There are six members of the tubulin superfamily in eukaryotes [1]. Alpha- and beta-tubulin form a heterodimer that polymerizes to form microtubules, and gamma-tubulin nucleates microtubules as a component of the gamma-tubulin ring complex. Alpha-, beta-, and gamma-tubulin are conserved in all eukaryotes. In contrast, delta- and epsilon-tubulin are conserved in many, but not all, eukaryotes and are associated with centrioles, although their molecular function is unclear [2–7]. Zeta-tubulin is the sixth and final member of the tubulin superfamily and is largely uncharacterized. We find that zeta-, epsilon-, and delta-tubulin form an evolutionarily co-conserved module, the ZED module, that has been lost at several junctions in eukaryotic evolution and that zeta- and delta-tubulin are evolutionarily interchangeable. Humans lack zeta-tubulin but have delta-tubulin. In Xenopus multiciliated cells, zeta-tubulin is a component of the basal foot, a centriolar appendage that connects centrioles to the apical cytoskeleton, and co-localizes there with epsilon-tubulin. Depletion of zeta-tubulin results in disorganization of centriole distribution and polarity in multiciliated cells. In contrast with multiciliated cells, zeta-tubulin in cycling cells does not localize to centrioles and is associated with the TRiC/CCT cytoplasmic chaperone complex. We conclude that zeta-tubulin facilitates interactions between the centrioles and the apical cytoskeleton as a component of the basal foot in differentiated cells and propose that the ZED tubulins are important for centriole functionalization and orientation of centrioles with respect to cellular polarity axes.

Xenopus is an important model organism for the study of genome duplication in vertebrates. With the full genome sequence of diploid Xenopus tropicalis available, and that of allotetraploid X. laevis close to being finished, we will be able to expand our understanding of how duplicated genes have evolved. One of the key features in the study of the functional consequence of gene duplication is how their expression patterns vary across different conditions, and RNA-seq seems to have enough resolution to discriminate the expression of highly similar duplicated genes. However, most of the current RNA-seq analysis methods were not designed to study samples with duplicate genes such as in X. laevis. Here, various computational methods to quantify gene expression in RNA-seq data were evaluated, using 2 independent X. laevis egg RNA-seq datasets and 2 reference databases for duplicated genes. The fact that RNA-seq can measure expression levels of similar duplicated genes was confirmed, but long paired-end reads are more informative than short single-end reads to discriminate duplicated genes. Also, it was found that bowtie, one of the most popular mappers in RNA-seq analysis, reports significantly smaller numbers of unique hits according to a mapping quality score compared to other mappers tested (BWA, GSNAP, STAR). Calculated from unique hits based on a mapping quality score, both expression levels and the expression ratio of duplicated genes can be estimated consistently among biological replicates, demonstrating that this method can successfully discriminate the expression of each copy of a duplicated gene pair. This comprehensive evaluation will be a useful guideline for studying gene expression of organisms with genome duplication using RNA-seq in the future.

Genome duplication creates redundancy in proteins and their interaction networks, and subsequent smaller-scale gene duplication can further amplify genetic redundancy. Mutations then lead to the loss, maintenance or functional divergence of duplicated genes. Genome duplication occurred many times in African clawed frogs (genus Xenopus), and almost all extant species in this group evolved from a polyploid ancestor. To better understand the nature of selective constraints in a polyploid genome, we examined molecular polymorphism and divergence of duplicates and single-copy genes in 2 tetraploid African clawed frog species, Xenopus laevis and X. victorianus. We found that molecular polymorphism in the coding regions of putative duplicated genes was higher than in singletons, but not significantly so. Our findings also suggest that transcriptome evolution in polyploids is influenced by variation in the genome-wide mutation rate, and do not reject the hypothesis that gene dosage balance is also important.

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.

The inner mitochondrial membrane (IMM) proteome plays a central role in maintaining mitochondrial physiology and cellular metabolism. Various important biochemical reactions such as oxidative phosphorylation, metabolite production, and mitochondrial biogenesis are conducted by the IMM proteome, and mitochondria-targeted therapeutics have been developed for IMM proteins, which is deeply related for various human metabolic diseases including cancer and neurodegenerative diseases. However, the membrane topology of the IMM proteome remains largely unclear because of the lack of methods to evaluate it in live cells in a high-throughput manner. In this article, we reveal the in vivo topological direction of 135 IMM proteins, using an in situ-generated radical probe with genetically targeted peroxidase (APEX). Owing to the short lifetime of phenoxyl radicals generated in situ by submitochondrial targeted APEX and the impermeability of the IMM to small molecules, the solvent-exposed tyrosine residues of both the matrix and intermembrane space (IMS) sides of IMM proteins were exclusively labeled with the radical probe in live cells by Matrix-APEX and IMS-APEX, respectively and identified by mass spectrometry. From this analysis, we confirmed 58 IMM protein topologies and we could determine the topological direction of 77 IMM proteins whose topology at the IMM has not been fully characterized. We also found several IMM proteins (e.g., LETM1 and OXA1) whose topological information should be revised on the basis of our results. Overall, our identification of structural information on the mitochondrial inner-membrane proteome can provide valuable insights for the architecture and connectome of the IMM proteome in live cells.

This is the open peer reviewers comments and recommendations regarding the submitted GigaScience article and/or dataset.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

During vertebrate retinogenesis, the precise balance between retinoblast proliferation and differentiation is spatially and temporally regulated through a number of intrinsic factors and extrinsic signaling pathways. Moreover, there are complex gene regulatory network interactions between these intrinsic factors and extrinsic pathways, which ultimately function to determine when retinoblasts exit the cell cycle and terminally differentiate. We recently uncovered a cell non-autonomous role for the intrinsic HLH factor, Id2a, in regulating retinoblast proliferation and differentiation, with Id2a-deficient retinae containing an abundance of proliferative retinoblasts and an absence of terminally differentiated retinal neurons and glia. Here, we report that Id2a function is necessary and sufficient to limit Notch pathway activity during retinogenesis. Id2a-deficient retinae possess elevated levels of Notch pathway component gene expression, while retinae overexpressing id2a possess reduced expression of Notch pathway component genes. Attenuation of Notch signaling activity by DAPT or by morpholino knockdown of Notch1a is sufficient to rescue both the proliferative and differentiation defects in Id2a-deficient retinae. In addition to regulating Notch pathway activity, through a novel RNA-Seq and differential gene expression analysis of Id2a-deficient retinae, we identify a number of additional intrinsic and extrinsic regulatory pathway components whose expression is regulated by Id2a. These data highlight the integral role played by Id2a in the gene regulatory network governing the transition from retinoblast proliferation to terminal differentiation during vertebrate retinogenesis.

Shotgun proteomics using mass spectrometry is a powerful method for protein identification but suffers limited sensitivity in complex samples. Integrating peptide identifications from multiple database search engines is a promising strategy to increase the number of peptide identifications and reduce the volume of unassigned tandem mass spectra. Existing methods pool statistical significance scores such as p-values or posterior probabilities of peptide-spectrum matches (PSMs) from multiple search engines after high scoring peptides have been assigned to spectra, but these methods lack reliable control of identification error rates as data are integrated from different search engines. We developed a statistically coherent method for integrative analysis, termed MSblender. MSblender converts raw search scores from search engines into a probability score for all possible PSMs and properly accounts for the correlation between search scores. The method reliably estimates false discovery rates and identifies more PSMs than any single search engine at the same false discovery rate. Increased identifications increment spectral counts for all detected proteins and allow quantification of proteins that would not have been quantified by individual search engines. We also demonstrate that enhanced quantification contributes to improve sensitivity in differential expression analyses.

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.

Motivation: High-throughput protein identification experiments based on tandem mass spectrometry (MS/MS) often suffer from low sensitivity and low-confidence protein identifications. In a typical shotgun proteomics experiment, it is assumed that all proteins are equally likely to be present. However, there is often other evidence to suggest that a protein is present and confidence in individual protein identification can be updated accordingly. Results: We develop a method that analyzes MS/MS experiments in the larger context of the biological processes active in a cell. Our method, MSNet, improves protein identification in shotgun proteomics experiments by considering information on functional associations from a gene functional network. MSNet substantially increases the number of proteins identified in the sample at a given error rate. We identify 8–29% more proteins than the original MS experiment when applied to yeast grown in different experimental conditions analyzed on different MS/MS instruments, and 37% more proteins in a human sample. We validate up to 94% of our identifications in yeast by presence in ground-truth reference sets. Availability and Implementation: Software and datasets are available at http://aug.csres.utexas.edu/msnet Contact: miranker@cs.utexas.edu, marcotte@icmb.utexas.edu Supplementary information: Supplementary data are available at Bioinformatics online.

For a long time, it has been assumed that the only role of sperm at fertilization is to introduce the male genome into the egg. Recently, ideas have emerged that the epigenetic state of the sperm nucleus could influence transcription in the embryo. However, conflicting reports have challenged the existence of epigenetic marks on sperm genes, and there are no functional tests supporting the role of sperm epigenetic marking on embryonic gene expression. Here, we show that sperm is epigenetically programmed to regulate embryonic gene expression. By comparing the development of sperm- and spermatid-derived frog embryos, we show that the programming of sperm for successful development relates to its ability to regulate transcription of a set of developmentally important genes. During spermatid maturation into sperm, these genes lose H3K4me2/3 and retain H3K27me3 marks. Experimental removal of these epigenetic marks at fertilization de-regulates gene expression in the resulting embryos in a paternal chromatin-dependent manner. This demonstrates that epigenetic instructions delivered by the sperm at fertilization are required for correct regulation of gene expression in the future embryos. The epigenetic mechanisms of developmental programming revealed here are likely to relate to the mechanisms involved in transgenerational transmission of acquired traits. Understanding how parental experience can influence development of the progeny has broad potential for improving human health.

Genome duplication creates redundancy in proteins and their interaction networks, and subsequent smaller-scale gene duplication can further amplify genetic redundancy. Mutations then lead to the loss, maintenance or functional divergence of duplicated genes. Genome duplication occurred many times in African clawed frogs (genus Xenopus), and almost all extant species in this group evolved from a polyploid ancestor. To better understand the nature of selective constraints in a polyploid genome, we examined molecular polymorphism and divergence of duplicates and single-copy genes in 2 tetraploid African clawed frog species, Xenopus laevis and X. victorianus. We found that molecular polymorphism in the coding regions of putative duplicated genes was higher than in singletons, but not significantly so. Our findings also suggest that transcriptome evolution in polyploids is influenced by variation in the genome-wide mutation rate, and do not reject the hypothesis that gene dosage balance is also important.

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. DOI: http://dx.doi.org/10.7554/eLife.01439.001, eLife digest Cells that have hundreds of tiny hair-like structures called cilia on their surface have important roles in our airways and also in the brain and reproductive system. By beating in a coordinated manner, the cilia cause fluid to flow in a particular direction. The development of these multiciliated cells is a complex process in which genes are expressed as proteins, with this gene expression being regulated by other proteins called transcription factors. In invertebrates the development of the cilia is controlled by transcription factors from the RFX family, which also appear to be important for development of cilia in vertebrates. However, the details of this process—in particular, the identities of the genes that are involved and how their functions are related—are not well understood in vertebrates. Chung et al. have sought to remedy this by analyzing the network of genes whose expression is controlled by the transcription factor Rfx2 in vertebrates. The results showed that the genes controlled by Rfx2 were involved in all aspects of cilia, including several genes that are known to be mutated in diseases caused by abnormal cilia. Chung et al. also identified genes that were not previously thought to be relevant to cilia. As multiciliated cells are developing, but before they can generate cilia, they must first migrate from the bottom of the epithelium, the layer of tissue in which they function, to the top of this layer. Chung et al. found that Rfx2 was also involved in this process. The approach taken by Chung et al.—which involved a combination of RNA sequence analysis, examination of Rfx2 binding sites on chromosomes, computational predictions of protein interactions and in vivo cellular imaging—could be used to perform similar systems-level analyses of other developmental and biological processes. DOI: http://dx.doi.org/10.7554/eLife.01439.002

International audience; A major goal in proteomics is the comprehensive and accurate description of a proteome. This task includes not only the identification of proteins in a sample, but also the accurate quantification of their abundance. Although mass spectrometry typically provides information on peptide identity and abundance in a sample, it does not directly measure the concentration of the corresponding proteins. Specifically, most mass-spectrometry-based approaches (e.g. shotgun proteomics or selected reaction monitoring) allow one to quantify peptides using chromatographic peak intensities or spectral counting information. Ultimately, based on these measurements, one wants to infer the concentrations of the corresponding proteins. Inferring properties of the proteins based on experimental peptide evidence is often a complex problem because of the ambiguity of peptide assignments and different chemical properties of the peptides that affect the observed concentrations. We present SCAMPI, a novel generic and statistically sound framework for computing protein abundance scores based on quantified peptides. In contrast to most previous approaches, our model explicitly includes information from shared peptides to improve protein quantitation, especially in eukaryotes with many homologous sequences. The model accounts for uncertainty in the input data, leading to statistical prediction intervals for the protein scores. Furthermore, peptides with extreme abundances can be reassessed and classified as either regular data points or actual outliers. We used the proposed model with several datasets and compared its performance to that of other, previously used approaches for protein quantification in bottom-up mass spectrometry.

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.

Gram-negative bacteria produce outer membrane vesicles (OMVs) that package and deliver proteins, small molecules, and DNA to prokaryotic and eukaryotic cells. The molecular details of OMV biogenesis have not been fully elucidated, but peptidoglycan-associated outer membrane proteins that tether the outer membrane to the underlying peptidoglycan have been shown to be critical for OMV formation in multiple Enterobacteriaceae . In this study, we demonstrate that the peptidoglycan-associated outer membrane proteins OprF and OprI, but not OprL, impact production of OMVs by the opportunistic pathogen Pseudomonas aeruginosa . Interestingly, OprF does not appear to be important for tethering the outer membrane to peptidoglycan but instead impacts OMV formation through modulation of the levels of the Pseudomonas quinolone signal (PQS), a quorum signal previously shown by our laboratory to be critical for OMV formation. Thus, the mechanism by which OprF impacts OMV formation is distinct from that for other peptidoglycan-associated outer membrane proteins, including OprI.

Proteins play major roles in most biological processes; as a consequence, protein expression levels are highly regulated. While extensive post-transcriptional, translational and protein degradation control clearly influence protein concentration and functionality, it is often thought that protein abundances are primarily determined by the abundances of the corresponding mRNAs. Hence surprisingly, a recent study showed that abundances of orthologous nematode and fly proteins correlate better than their corresponding mRNA abundances. We tested if this phenomenon is general by collecting and testing matching large-scale protein and mRNA expression data sets from seven different species: two bacteria, yeast, nematode, fly, human, and rice. We find that steady-state abundances of proteins show significantly higher correlation across these diverse phylogenetic taxa than the abundances of their corresponding mRNAs (p=0.0008, paired Wilcoxon). These data support the presence of strong selective pressure to maintain protein abundances during evolution, even when mRNA abundances diverge.

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.

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., Author Summary Duplicate genes occur in all organisms. It has been found that mutations in duplicate genes cause defects much less often than when single copy genes are mutated. It is widely believed that this is due to functional redundancy—that is, the two genes can carry out similar functions so that the non-mutated duplicate gene can cover for or “mask” the phenotype of the mutation in the first duplicate. To determine whether this hypothesis is true, it is necessary to test systematically whether defects indeed occur in the organism when both duplicate genes are inhibited. We have for the first time carried out such an analysis in a multicellular organism, the nematode Caenorhabditis elegans. In contrast to expectations, we observed that when both copies of duplicate genes are inhibited deleterious effects are very rare. We show that this is because duplicate genes are much more often non-essential compared to genes where there is only a single copy. Non-essential genes are also lost from the genome much more often than essential genes. However, when essential genes are duplicated, they remain present in the genome over longer periods. Our results give a framework to explain the evolutionary dynamics of duplications in the genome.