Somatic Mutations in Renal Cyst Epithelium in Autosomal Dominant Polycystic Kidney Disease.

BACKGROUND: Autosomal dominant polycystic kidney disease (ADPKD) is a ciliopathy caused by mutations in PKD1 and PKD2 that is characterized by renal tubular epithelial cell proliferation and progressive CKD. Although the molecular mechanisms involved in cystogenesis are not established, concurrent inactivating constitutional and somatic mutations in ADPKD genes in cyst epithelium have been proposed as a cellular recessive mechanism. METHODS: We characterized, by whole-exome sequencing (WES) and long-range PCR techniques, the somatic mutations in PKD1 and PKD2 genes in renal epithelial cells from 83 kidney cysts obtained from nine patients with ADPKD, for whom a constitutional mutation in PKD1 or PKD2 was identified. RESULTS: Complete sequencing data by long-range PCR and WES was available for 63 and 65 cysts, respectively. Private somatic mutations of PKD1 or PKD2 were identified in all patients and in 90% of the cysts analyzed; 90% of these mutations were truncating, splice site, or in-frame variations predicted to be pathogenic mutations. No trans-heterozygous mutations of PKD1 or PKD2 genes were identified. Copy number changes of PKD1 ranging from 151 bp to 28 kb were observed in 12% of the cysts. WES also identified significant mutations in 53 non-PKD1/2 genes, including other ciliopathy genes and cancer-related genes. CONCLUSIONS: These findings support a cellular recessive mechanism for cyst formation in ADPKD caused primarily by inactivating constitutional and somatic mutations of PKD1 or PKD2 in kidney cyst epithelium. The potential interactions of these genes with other ciliopathy- and cancer-related genes to influence ADPKD severity merits further evaluation.

Drug repurposing to target Ebola virus replication and virulence using structural systems pharmacology.

BACKGROUND: The recent outbreak of Ebola has been cited as the largest in history. Despite this global health crisis, few drugs are available to efficiently treat Ebola infections. Drug repurposing provides a potentially efficient solution to accelerating the development of therapeutic approaches in response to Ebola outbreak. To identify such candidates, we use an integrated structural systems pharmacology pipeline which combines proteome-scale ligand binding site comparison, protein-ligand docking, and Molecular Dynamics (MD) simulation. RESULTS: One thousand seven hundred and sixty-six FDA-approved drugs and 259 experimental drugs were screened to identify those with the potential to inhibit the replication and virulence of Ebola, and to determine the binding modes with their respective targets. Initial screening has identified a number of promising hits. Notably, Indinavir; an HIV protease inhibitor, may be effective in reducing the virulence of Ebola. Additionally, an antifungal (Sinefungin) and several anti-viral drugs (e.g. Maraviroc, Abacavir, Telbivudine, and Cidofovir) may inhibit Ebola RNA-directed RNA polymerase through targeting the MTase domain. CONCLUSIONS: Identification of safe drug candidates is a crucial first step toward the determination of timely and effective therapeutic approaches to address and mitigate the impact of the Ebola global crisis and future outbreaks of pathogenic diseases. Further in vitro and in vivo testing to evaluate the anti-Ebola activity of these drugs is warranted.

Identification, modeling, and characterization studies of Tetrahymena thermophila myosin FERM domains suggests a conserved core fold but functional differences.

Myosins (MYO) define a superfamily of motor proteins which facilitate movement along cytoskeletal actin filaments in an ATP-dependent manner. To date, over 30 classes of myosin have been defined that vary in their roles and distribution across different taxa. The multidomain tail of myosin is responsible for the observed functional differences in different myosin classes facilitating differential binding to different cargos. One domain found in this region, the FERM domain, is found in several diverse proteins and is involved in many biological functions ranging from cell adhesion and actin-driven cytoskeleton assembly to cell signaling. Recently, new classes of unconventional myosin have been identified in Tetrahymena thermophila. In this study, we have identified, modeled, and characterized eight FERM domains from the unconventional T. thermophila myosins as their complete functional MyTH4-FERM cassettes. Our results reveal notable sequence, structural, and electrostatic differences between T. thermophila and other characterized FERM domains. Specifically, T. thermophila FERM domains contain helical inserts or extensions, which contribute to significant differences in surface electrostatic profiles of T. thermophila myosin FERMs when compared to the conventional FERM domains. Analyses of the modeled domains reveal differences in key functional residues as well as phosphoinositide-binding signatures and affinities. The work presented here broadens the scope of our understanding of myosin classes and their inherent functions, and provides a platform for experimentalists to design rational experimental studies to test the functional roles for T. thermophila myosins.

Phylogenomic identification of regulatory sequences in bacteria: an analysis of statistical power and an application to Borrelia burgdorferi sensu lato.

UNLABELLED: Phylogenomic footprinting is an approach for ab initio identification of genome-wide regulatory elements in bacterial species based on sequence conservation. The statistical power of the phylogenomic approach depends on the degree of sequence conservation, the length of regulatory elements, and the level of phylogenetic divergence among genomes. Building on an earlier model, we propose a binomial model that uses synonymous tree lengths as neutral expectations for determining the statistical significance of conserved intergenic spacer (IGS) sequences. Simulations show that the binomial model is robust to variations in the value of evolutionary parameters, including base frequencies and the transition-to-transversion ratio. We used the model to search for regulatory sequences in the Lyme disease species group (Borrelia burgdorferi sensu lato) using 23 genomes. The model indicates that the currently available set of Borrelia genomes would not yield regulatory sequences shorter than five bases, suggesting that genome sequences of additional B. burgdorferi sensu lato species are needed. Nevertheless, we show that previously known regulatory elements are indeed strongly conserved in sequence or structure across these Borrelia species. Further, we predict with sufficient confidence two new RpoS binding sites, 39 promoters, 19 transcription terminators, 28 noncoding RNAs, and four sets of coregulated genes. These putative cis- and trans-regulatory elements suggest novel, Borrelia-specific mechanisms regulating the transition between the tick and host environments, a key adaptation and virulence mechanism of B. burgdorferi. Alignments of IGS sequences are available on BorreliaBase.org, an online database of orthologous open reading frame (ORF) and IGS sequences in Borrelia. IMPORTANCE: While bacterial genomes contain mostly protein-coding genes, they also house DNA sequences regulating the expression of these genes. Gene regulatory sequences tend to be conserved during evolution. By sequencing and comparing related genomes, one can therefore identify regulatory sequences in bacteria based on sequence conservation. Here, we describe a statistical framework by which one may determine how many genomes need to be sequenced and at what level of evolutionary relatedness in order to achieve a high level of statistical significance. We applied the framework to Borrelia burgdorferi, the Lyme disease agent, and identified a large number of candidate regulatory sequences, many of which are known to be involved in regulating the phase transition between the tick vector and mammalian hosts.

Proteome-wide analysis of mutant p53 targets in breast cancer identifies new levels of gain-of-function that influence PARP, PCNA, and MCM4.

The gain-of-function mutant p53 (mtp53) transcriptome has been studied, but, to date, no detailed analysis of the mtp53-associated proteome has been described. We coupled cell fractionation with stable isotope labeling with amino acids in cell culture (SILAC) and inducible knockdown of endogenous mtp53 to determine the mtp53-driven proteome. Our fractionation data highlight the underappreciated biology that missense mtp53 proteins R273H, R280K, and L194F are tightly associated with chromatin. Using SILAC coupled to tandem MS, we identified that R273H mtp53 expression in MDA-MB-468 breast cancer cells up- and down-regulated multiple proteins and metabolic pathways. Here we provide the data set obtained from sequencing 73,154 peptide pairs that then corresponded to 3,010 proteins detected under reciprocal labeling conditions. Importantly, the high impact regulated targets included the previously identified transcriptionally regulated mevalonate pathway proteins but also identified two new levels of mtp53 protein regulation for nontranscriptional targets. Interestingly, mtp53 depletion profoundly influenced poly(ADP ribose) polymerase 1 (PARP1) localization, with increased cytoplasmic and decreased chromatin-associated protein. An enzymatic PARP shift occurred with high mtp53 expression, resulting in increased poly-ADP-ribosylated proteins in the nucleus. Mtp53 increased the level of proliferating cell nuclear antigen (PCNA) and minichromosome maintenance 4 (MCM4) proteins without changing the amount of pcna and mcm4 transcripts. Pathway enrichment analysis ranked the DNA replication pathway above the cholesterol biosynthesis pathway as a R273H mtp53 activated proteomic target. Knowledge of the proteome diversity driven by mtp53 suggests that DNA replication and repair pathways are major targets of mtp53 and highlights consideration of combination chemotherapeutic strategies targeting cholesterol biosynthesis and PARP inhibition.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes.

BACKGROUND: The bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes. DESCRIPTION: We have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers. CONCLUSIONS: With a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.

Evolutionary genomics of Borrelia burgdorferi sensu lato: findings, hypotheses, and the rise of hybrids.

Borrelia burgdorferi sensu lato (B. burgdorferi s.l.), the group of bacterial species represented by Lyme disease pathogens, has one of the most complex and variable genomic architectures among prokaryotes. Showing frequent recombination within and limited gene flow among geographic populations, the B. burgdorferi s.l. genomes provide an excellent window into the processes of bacterial evolution at both within- and between-population levels. Comparative analyses of B. burgdorferi s.l. genomes revealed a highly dynamic plasmid composition but a conservative gene repertoire. Gene duplication and loss as well as sequence variations at loci encoding surface-localized lipoproteins (e.g., the PF54 genes) are strongly associated with adaptive differences between species. There are a great many conserved intergenic spacer sequences that are candidates for cis-regulatory elements and non-coding RNAs. Recombination among coexisting strains occurs at a rate approximately three times the mutation rate. The coexistence of a large number of genomic groups within local B. burgdorferi s.l. populations may be driven by immune-mediated diversifying selection targeting major antigen loci as well as by adaptation to multiple host species. Questions remain regarding the ecological causes (e.g., climate change, host movements, or new adaptations) of the ongoing range expansion of B. burgdorferi s.l. and on the genomic variations associated with its ecological and clinical variability. Anticipating an explosive growth of the number of B. burgdorferi s.l. genomes sampled from both within and among species, we propose genome-based methods to test adaptive mechanisms and to identify molecular bases of phenotypic variations. Genome sequencing is also necessary for monitoring a likely increase of genetic admixture of previously isolated species and populations in North America and elsewhere.

Mutational and phylogenetic analyses of the mycobacterial mbt gene cluster.

The mycobactin siderophore system is present in many Mycobacterium species, including M. tuberculosis and other clinically relevant mycobacteria. This siderophore system is believed to be utilized by both pathogenic and nonpathogenic mycobacteria for iron acquisition in both in vivo and ex vivo iron-limiting environments, respectively. Several M. tuberculosis genes located in a so-called mbt gene cluster have been predicted to be required for the biosynthesis of the core scaffold of mycobactin based on sequence analysis. A systematic and controlled mutational analysis probing the hypothesized essential nature of each of these genes for mycobactin production has been lacking. The degree of conservation of mbt gene cluster orthologs remains to be investigated as well. In this study, we sought to conclusively establish whether each of nine mbt genes was required for mycobactin production and to examine the conservation of gene clusters orthologous to the M. tuberculosis mbt gene cluster in other bacteria. We report a systematic mutational analysis of the mbt gene cluster ortholog found in Mycobacterium smegmatis. This mutational analysis demonstrates that eight of the nine mbt genes investigated are essential for mycobactin production. Our genome mining and phylogenetic analyses reveal the presence of orthologous mbt gene clusters in several bacterial species. These gene clusters display significant organizational differences originating from an intricate evolutionary path that might have included horizontal gene transfers. Altogether, the findings reported herein advance our understanding of the genetic requirements for the biosynthesis of an important mycobacterial secondary metabolite with relevance to virulence.