A conserved salt bridge in the G loop of multiple protein kinases is important for catalysis and for in vivo Lyn function.

The glycine-rich G loop controls ATP binding and phosphate transfer in protein kinases. Here we show that the functions of Src family and Abl protein tyrosine kinases require an electrostatic interaction between oppositely charged amino acids within their G loops that is conserved in multiple other phylogenetically distinct protein kinases, from plants to humans. By limiting G loop flexibility, it controls ATP binding, catalysis, and inhibition by ATP-competitive compounds such as Imatinib. In WeeB mice, mutational disruption of the interaction results in expression of a Lyn protein with reduced catalytic activity, and in perturbed B cell receptor signaling. Like Lyn(-/-) mice, WeeB mice show profound defects in B cell development and function and succumb to autoimmune glomerulonephritis. This demonstrates the physiological importance of the conserved G loop salt bridge and at the same time distinguishes the in vivo requirement for the Lyn kinase activity from other potential functions of the protein.

Fitness costs of rifampicin resistance in Mycobacterium tuberculosis are amplified under conditions of nutrient starvation and compensated by mutation in the ß' subunit of RNA polymerase.

Rifampicin resistance, a defining attribute of multidrug-resistant tuberculosis, is conferred by mutations in the ß subunit of RNA polymerase. Sequencing of rifampicin-resistant (RIF-R) clinical isolates of Mycobacterium tuberculosis revealed, in addition to RIF-R mutations, enrichment of potential compensatory mutations around the double-psi ß-barrel domain of the ß' subunit comprising the catalytic site and the exit tunnel for newly synthesized RNA. Sequential introduction of the resistance allele followed by the compensatory allele in isogenic Mycobacterium smegmatis showed that these mutations respectively caused and compensated a starvation enhanced growth defect by altering RNA polymerase activity. While specific combinations of resistance and compensatory alleles converged in divergent lineages, other combinations recurred among related isolates suggesting transmission of compensated RIF-R strains. These findings suggest nutrient poor growth conditions impose larger selective pressure on RIF-R organisms that results in the selection of compensatory mutations in a domain involved in catalysis and starvation control of RNA polymerase transcription.

Identification of synthetic lethality of PRKDC in MYC-dependent human cancers by pooled shRNA screening.

BACKGROUND: MYC family members are among the most frequently deregulated oncogenes in human cancers, yet direct therapeutic targeting of MYC in cancer has been challenging thus far. Synthetic lethality provides an opportunity for therapeutic intervention of MYC-driven cancers. METHODS: A pooled kinase shRNA library screen was performed and next-generation deep sequencing efforts identified that PRKDC was synthetically lethal in cells overexpressing MYC. Genes and proteins of interest were knocked down or inhibited using RNAi technology and small molecule inhibitors, respectively. Quantitative RT-PCR using TaqMan probes examined mRNA expression levels and cell viability was assessed using CellTiter-Glo (Promega). Western blotting was performed to monitor different protein levels in the presence or absence of RNAi or compound treatment. Statistical significance of differences among data sets were determined using unpaired t test (Mann-Whitney test) or ANOVA. RESULTS: Inhibition of PRKDC using RNAi (RNA interference) or small molecular inhibitors preferentially killed MYC-overexpressing human lung fibroblasts. Moreover, inducible PRKDC knockdown decreased cell viability selectively in high MYC-expressing human small cell lung cancer cell lines. At the molecular level, we found that inhibition of PRKDC downregulated MYC mRNA and protein expression in multiple cancer cell lines. In addition, we confirmed that overexpression of MYC family proteins induced DNA double-strand breaks; our results also revealed that PRKDC inhibition in these cells led to an increase in DNA damage levels. CONCLUSIONS: Our data suggest that the synthetic lethality between PRKDC and MYC may in part be due to PRKDC dependent modulation of MYC expression, as well as MYC-induced DNA damage where PRKDC plays a key role in DNA damage repair.

Whole-genome sequencing and microarray analysis of ex vivo Plasmodium vivax reveal selective pressure on putative drug resistance genes.

Plasmodium vivax causes 25-40% of malaria cases worldwide, yet research on this human malaria parasite has been neglected. Nevertheless, the recent publication of the P. vivax reference genome now allows genomics and systems biology approaches to be applied to this pathogen. We show here that whole-genome analysis of the parasite can be achieved directly from ex vivo-isolated parasites, without the need for in vitro propagation. A single isolate of P. vivax obtained from a febrile patient with clinical malaria from Peru was subjected to whole-genome sequencing (30× coverage). This analysis revealed over 18,261 single-nucleotide polymorphisms (SNPs), 6,257 of which were further validated using a tiling microarray. Within core chromosomal genes we find that one SNP per every 985 bases of coding sequence distinguishes this recent Peruvian isolate, designated IQ07, from the reference Salvador I strain obtained in 1972. This full-genome sequence of an uncultured P. vivax isolate shows that the same regions with low numbers of aligned sequencing reads are also highly variable by genomic microarray analysis. Finally, we show that the genes containing the largest ratio of nonsynonymous-to-synonymous SNPs include two AP2 transcription factors and the P. vivax multidrug resistance-associated protein (PvMRP1), an ABC transporter shown to be associated with quinoline and antifolate tolerance in Plasmodium falciparum. This analysis provides a data set for comparative analysis with important potential for identifying markers for global parasite diversity and drug resistance mapping studies.

Screening the mammalian extracellular proteome for regulators of embryonic human stem cell pluripotency.

Approximately 3,500 mammalian genes are predicted to be secreted or single-pass transmembrane proteins. The function of the majority of these genes is still unknown, and a number of the encoded proteins might find use as new therapeutic agents themselves or as targets for small molecule or antibody drug development. To analyze the physiological activities of the extracellular proteome, we developed a large-scale, high-throughput protein expression, purification, and screening platform. For this study, the complete human extracellular proteome was analyzed and prioritized based on genome-wide disease association studies to select 529 initial target genes. These genes were cloned into three expression vectors as native sequences and as N-terminal and C-terminal Fc fusions to create an initial collection of 806 purified secreted proteins. To determine its utility, this library was screened in an OCT4-based cellular assay to identify regulators of human embryonic stem-cell self-renewal. We found that the pigment epithelium-derived factor can promote long-term pluripotent growth of human embryonic stem cells without bFGF or TGFbeta/Activin/Nodal ligand supplementation. Our results further indicate that activation of the pigment epithelium-derived factor receptor-Erk1/2 signaling pathway by the pigment epithelium-derived factor is sufficient to maintain the self-renewal of pluripotent human embryonic stem cells. These experiments illustrate the potential for discovering novel biological functions by directly screening protein diversity in cell-based phenotypic or reporter assays.

Approaches to long-read sequencing in a clinical setting to improve diagnostic rate.

Over the past decade, advances in genetic testing, particularly the advent of next-generation sequencing, have led to a paradigm shift in the diagnosis of molecular diseases and disorders. Despite our present collective ability to interrogate more than 90% of the human genome, portions of the genome have eluded us, resulting in stagnation of diagnostic yield with existing methodologies. Here we show how application of a new technology, long-read sequencing, has the potential to improve molecular diagnostic rates. Whole genome sequencing by long reads was able to cover 98% of next-generation sequencing dead zones, which are areas of the genome that are not interpretable by conventional industry-standard short-read sequencing. Through the ability of long-read sequencing to unambiguously call variants in these regions, we discovered an immunodeficiency due to a variant in IKBKG in a subject who had previously received a negative genome sequencing result. Additionally, we demonstrate the ability of long-read sequencing to detect small variants on par with short-read sequencing, its superior performance in identifying structural variants, and thirdly, its capacity to determine genomic methylation defects in native DNA. Though the latter technical abilities have been demonstrated, we demonstrate the clinical application of this technology to successfully identify multiple types of variants using a single test.

Gene set enrichment in eQTL data identifies novel annotations and pathway regulators.

Genome-wide gene expression profiling has been extensively used to generate biological hypotheses based on differential expression. Recently, many studies have used microarrays to measure gene expression levels across genetic mapping populations. These gene expression phenotypes have been used for genome-wide association analyses, an analysis referred to as expression QTL (eQTL) mapping. Here, eQTL analysis was performed in adipose tissue from 28 inbred strains of mice. We focused our analysis on "trans-eQTL bands", defined as instances in which the expression patterns of many genes were all associated to a common genetic locus. Genes comprising trans-eQTL bands were screened for enrichments in functional gene sets representing known biological pathways, and genes located at associated trans-eQTL band loci were considered candidate transcriptional modulators. We demonstrate that these patterns were enriched for previously characterized relationships between known upstream transcriptional regulators and their downstream target genes. Moreover, we used this strategy to identify both novel regulators and novel members of known pathways. Finally, based on a putative regulatory relationship identified in our analysis, we identified and validated a previously uncharacterized role for cyclin H in the regulation of oxidative phosphorylation. We believe that the specific molecular hypotheses generated in this study will reveal many additional pathway members and regulators, and that the analysis approaches described herein will be broadly applicable to other eQTL data sets.

A coactivator trap identifies NONO (p54nrb) as a component of the cAMP-signaling pathway.

Signal transduction pathways often use a transcriptional component to mediate adaptive cellular responses. Coactivator proteins function prominently in these pathways as the conduit to the basic transcriptional machinery. Here we present a high-throughput cell-based screening strategy, termed the "coactivator trap," to study the functional interactions of coactivators with transcription factors. We applied this strategy to the cAMP signaling pathway, which utilizes two families of coactivators, the cAMP response element binding protein (CREB) binding protein (CBP)/p300 family and the recently identified transducers of regulated CREB activity family (TORCs1-3). In addition to identifying numerous known interactions of these coactivators, this analysis identified NONO (p54(nrb)) as a TORC-interacting protein. RNA interference experiments demonstrate that NONO is necessary for cAMP-dependent activation of CREB target genes in vivo. Furthermore, TORC2 and NONO complex on cAMP-responsive promoters, and NONO acts as a bridge between the CREB/TORC complex and RNA polymerase II. These data demonstrate the utility of the coactivator trap by identification of a component of cAMP-mediated transcription.

Complex Loci in human and mouse genomes.

Mammalian genomes harbor a larger than expected number of complex loci, in which multiple genes are coupled by shared transcribed regions in antisense orientation and/or by bidirectional core promoters. To determine the incidence, functional significance, and evolutionary context of mammalian complex loci, we identified and characterized 5,248 cis-antisense pairs, 1,638 bidirectional promoters, and 1,153 chains of multiple cis-antisense and/or bidirectionally promoted pairs from 36,606 mouse transcriptional units (TUs), along with 6,141 cis-antisense pairs, 2,113 bidirectional promoters, and 1,480 chains from 42,887 human TUs. In both human and mouse, 25% of TUs resided in cis-antisense pairs, only 17% of which were conserved between the two organisms, indicating frequent species specificity of antisense gene arrangements. A sampling approach indicated that over 40% of all TUs might actually be in cis-antisense pairs, and that only a minority of these arrangements are likely to be conserved between human and mouse. Bidirectional promoters were characterized by variable transcriptional start sites and an identifiable midpoint at which overall sequence composition changed strand and the direction of transcriptional initiation switched. In microarray data covering a wide range of mouse tissues, genes in cis-antisense and bidirectionally promoted arrangement showed a higher probability of being coordinately expressed than random pairs of genes. In a case study on homeotic loci, we observed extensive transcription of nonconserved sequences on the noncoding strand, implying that the presence rather than the sequence of these transcripts is of functional importance. Complex loci are ubiquitous, host numerous nonconserved gene structures and lineage-specific exonification events, and may have a cis-regulatory impact on the member genes.

Transcript annotation in FANTOM3: mouse gene catalog based on physical cDNAs.

The international FANTOM consortium aims to produce a comprehensive picture of the mammalian transcriptome, based upon an extensive cDNA collection and functional annotation of full-length enriched cDNAs. The previous dataset, FANTOM2, comprised 60,770 full-length enriched cDNAs. Functional annotation revealed that this cDNA dataset contained only about half of the estimated number of mouse protein-coding genes, indicating that a number of cDNAs still remained to be collected and identified. To pursue the complete gene catalog that covers all predicted mouse genes, cloning and sequencing of full-length enriched cDNAs has been continued since FANTOM2. In FANTOM3, 42,031 newly isolated cDNAs were subjected to functional annotation, and the annotation of 4,347 FANTOM2 cDNAs was updated. To accomplish accurate functional annotation, we improved our automated annotation pipeline by introducing new coding sequence prediction programs and developed a Web-based annotation interface for simplifying the annotation procedures to reduce manual annotation errors. Automated coding sequence and function prediction was followed with manual curation and review by expert curators. A total of 102,801 full-length enriched mouse cDNAs were annotated. Out of 102,801 transcripts, 56,722 were functionally annotated as protein coding (including partial or truncated transcripts), providing to our knowledge the greatest current coverage of the mouse proteome by full-length cDNAs. The total number of distinct non-protein-coding transcripts increased to 34,030. The FANTOM3 annotation system, consisting of automated computational prediction, manual curation, and final expert curation, facilitated the comprehensive characterization of the mouse transcriptome, and could be applied to the transcriptomes of other species.

Genomewide association analysis in diverse inbred mice: power and population structure.

The discovery of quantitative trait loci (QTL) in model organisms has relied heavily on the ability to perform controlled breeding to generate genotypic and phenotypic diversity. Recently, we and others have demonstrated the use of an existing set of diverse inbred mice (referred to here as the mouse diversity panel, MDP) as a QTL mapping population. The use of the MDP population has many advantages relative to traditional F(2) mapping populations, including increased phenotypic diversity, a higher recombination frequency, and the ability to collect genotype and phenotype data in community databases. However, these methods are complicated by population structure inherent in the MDP and the lack of an analytical framework to assess statistical power. To address these issues, we measured gene expression levels in hypothalamus across the MDP. We then mapped these phenotypes as quantitative traits with our association algorithm, resulting in a large set of expression QTL (eQTL). We utilized these eQTL, and specifically cis-eQTL, to develop a novel nonparametric method for association analysis in structured populations like the MDP. These eQTL data confirmed that the MDP is a suitable mapping population for QTL discovery and that eQTL results can serve as a gold standard for relative measures of statistical power.

A systematic map of genetic variation in Plasmodium falciparum.

Discovering novel genes involved in immune evasion and drug resistance in the human malaria parasite, Plasmodium falciparum, is of critical importance to global health. Such knowledge may assist in the development of new effective vaccines and in the appropriate use of antimalarial drugs. By performing a full-genome scan of allelic variability in 14 field and laboratory strains of P. falciparum, we comprehensively identified approximately 500 genes evolving at higher than neutral rates. The majority of the most variable genes have paralogs within the P. falciparum genome and may be subject to a different evolutionary clock than those without. The group of 211 variable genes without paralogs contains most known immunogens and a few drug targets, consistent with the idea that the human immune system and drug use is driving parasite evolution. We also reveal gene-amplification events including one surrounding pfmdr1, the P. falciparum multidrug-resistance gene, and a previously uncharacterized amplification centered around the P. falciparum GTP cyclohydrolase gene, the first enzyme in the folate biosynthesis pathway. Although GTP cyclohydrolase is not the known target of any current drugs, downstream members of the pathway are targeted by several widely used antimalarials. We speculate that an amplification of the GTP cyclohydrolase enzyme in the folate biosynthesis pathway may increase flux through this pathway and facilitate parasite resistance to antifolate drugs.

Concomitant diagnosis of immune deficiency and Pseudomonas sepsis in a 19 month old with ecthyma gangrenosum by host whole-genome sequencing.

X-linked agammaglobulinemia (XLA, OMIM#300300) is a rare monogenic primary immunodeficiency caused by mutations in the Bruton tyrosine kinase (BTK) gene. XLA is characterized by insufficient immunoglobulin levels and susceptibility to life-threatening bacterial infections. We report on a patient that presented with ecthyma gangrenosum and septicemia. Rapid trio whole-genome sequencing (rWGS) revealed an apparently de novo hemizygous pathogenic variant (c.726dupT; p.Ile243TyrfsTer15) in the BTK gene. Metagenomic analysis of rWGS sequences that did not align to the human genome revealed 770 aligned to the Pseudomonas aeruginosa PAO1 genome. The patient was diagnosed with XLA and pseudomonal sepsis.

The NSIGHT1-randomized controlled trial: rapid whole-genome sequencing for accelerated etiologic diagnosis in critically ill infants.

Genetic disorders are a leading cause of morbidity and mortality in infants in neonatal and pediatric intensive care units (NICU/PICU). While genomic sequencing is useful for genetic disease diagnosis, results are usually reported too late to guide inpatient management. We performed an investigator-initiated, partially blinded, pragmatic, randomized, controlled trial to test the hypothesis that rapid whole-genome sequencing (rWGS) increased the proportion of NICU/PICU infants receiving a genetic diagnosis within 28 days. The participants were families with infants aged <4 months in a regional NICU and PICU, with illnesses of unknown etiology. The intervention was trio rWGS. Enrollment from October 2014 to June 2016, and follow-up until November 2016. Of all, 26 female infants, 37 male infants, and 2 infants of undetermined sex were randomized to receive rWGS plus standard genetic tests (n = 32, cases) or standard genetic tests alone (n = 33, controls). The study was terminated early due to loss of equipoise: 73% (24) controls received genomic sequencing as standard tests, and 15% (five) controls underwent compassionate cross-over to receive rWGS. Nevertheless, intention to treat analysis showed the rate of genetic diagnosis within 28 days of enrollment (the primary end-point) to be higher in cases (31%, 10 of 32) than controls (3%, 1 of 33; difference, 28% [95% CI, 10-46%]; p = 0.003). Among infants enrolled in the first 25 days of life, the rate of neonatal diagnosis was higher in cases (32%, 7 of 22) than controls (0%, 0 of 23; difference, 32% [95% CI, 11-53%];p = 0.004). Median age at diagnosis (25 days [range 14-90] in cases vs. 130 days [range 37-451] in controls) and median time to diagnosis (13 days [range 1-84] in cases, vs. 107 days [range 21-429] in controls) were significantly less in cases than controls (p = 0.04). In conclusion, rWGS increased the proportion of NICU/PICU infants who received timely diagnoses of genetic diseases.

Use of a dense single nucleotide polymorphism map for in silico mapping in the mouse.

Rapid expansion of available data, both phenotypic and genotypic, for multiple strains of mice has enabled the development of new methods to interrogate the mouse genome for functional genetic perturbations. In silico mapping provides an expedient way to associate the natural diversity of phenotypic traits with ancestrally inherited polymorphisms for the purpose of dissecting genetic traits. In mouse, the current single nucleotide polymorphism (SNP) data have lacked the density across the genome and coverage of enough strains to properly achieve this goal. To remedy this, 470,407 allele calls were produced for 10,990 evenly spaced SNP loci across 48 inbred mouse strains. Use of the SNP set with statistical models that considered unique patterns within blocks of three SNPs as an inferred haplotype could successfully map known single gene traits and a cloned quantitative trait gene. Application of this method to high-density lipoprotein and gallstone phenotypes reproduced previously characterized quantitative trait loci (QTL). The inferred haplotype data also facilitates the refinement of QTL regions such that candidate genes can be more easily identified and characterized as shown for adenylate cyclase 7.

Rapid whole-genome sequencing identifies a novel GABRA1 variant associated with West syndrome.

A 9-mo-old infant was admitted with infantile spasms that improved on administration of topiramate and steroids. He also had developmental delay, esotropia, and hypsarrhythmia on interictal electroencephalogram (EEG), and normal brain magnetic resonance imaging (MRI). West syndrome is the triad of infantile spasms, interictal hypsarrhythmia, and mental retardation. Rapid trio whole-genome sequencing (WGS) revealed a novel, likely pathogenic, de novo variant in the gene encoding γ-aminobutyric acid (GABA) type A receptor, α1 polypeptide (GABRA1 c.789G>A, p.Met263Ile) in the proband. GABRA1 mutations have been associated with early infantile epileptic encephalopathy type 19 (EIEE19). We suggest that GABRA1 p.Met263Ile is associated with a distinct West syndrome phenotype.

Mst1 directs Myosin IIa partitioning of low and higher affinity integrins during T cell migration.

Chemokines promote T cell migration by transmitting signals that induce T cell polarization and integrin activation and adhesion. Mst1 kinase is a key signal mediator required for both of these processes; however, its molecular mechanism remains unclear. Here, we present a mouse model in which Mst1 function is disrupted by a hypomorphic mutation. Microscopic analysis of Mst1-deficient CD4 T cells revealed a necessary role for Mst1 in controlling the localization and activity of Myosin IIa, a molecular motor that moves along actin filaments. Using affinity specific LFA-1 antibodies, we identified a requirement for Myosin IIa-dependent contraction in the precise spatial distribution of low and higher affinity LFA-1 on the membrane of migrating T cells. Mst1 deficiency or Myosin inhibition resulted in multipolar cells, difficulties in uropod detachment and mis-localization of low affinity LFA-1. Thus, Mst1 regulates Myosin IIa dynamics to organize high and low affinity LFA-1 to the anterior and posterior membrane during T cell migration.

Evolution of the human ion channel set.

Ion channels are intimately involved in virtually every physiological process of consequence in humans. Their importance is underscored by the identification of numerous "channelopathies", human diseases caused by ion channel mutations. Ion Channels have consequently been viewed as fertile ground for drug discovery and, indeed, they represent one of the largest target classes for current medicines. The future prospects of ion channels as a target class are tied to the functional characterization of the human ion channel set on a genomic scale. The focus of this review is to describe the molecular diversity and conservation of human ion channels. The human genome contains at least 232 genes that encode the pore-forming subunits of plasma membrane ion channels. Comparative genome analysis shows that most human ion channel gene families have their origins in the earliest metazoans but the human genes are largely derived from duplications that took place in the vertebrate lineage. The mouse and human ion channel gene sets are virtually identical, but differ significantly from fish channel sets. Genome comparisons highlight a number of highly conserved channel families that do not yet have specifically defined functional roles in vivo. These channel families are likely to have non-redundant functions in metazoans and represent some of the best new opportunities for channel target prospecting. Furthermore, genome-wide patterns of sequence conservation can now be used to refine strategies for the identification of gene-specific channel probes.

BioGPS: an extensible and customizable portal for querying and organizing gene annotation resources.

Online gene annotation resources are indispensable for analysis of genomics data. However, the landscape of these online resources is highly fragmented, and scientists often visit dozens of these sites for each gene in a candidate gene list. Here, we introduce BioGPS http://biogps.gnf.org, a centralized gene portal for aggregating distributed gene annotation resources. Moreover, BioGPS embraces the principle of community intelligence, enabling any user to easily and directly contribute to the BioGPS platform.