Comprehensive Classification of Retinal Bipolar Neurons by Single-Cell Transcriptomics

Patterns of gene expression can be used to characterize and classify neuronal types. It is challenging, however, to generate taxonomies that fulfill the essential criteria of being comprehensive, harmonizing with conventional classification schemes, and lacking superfluous subdivisions of genuine types. To address these challenges, we used massively parallel single-cell RNA profiling and optimized computational methods on a heterogeneous class of neurons, mouse retinal bipolar cells (BCs). From a population of ∼25,000 BCs, we derived a molecular classification that identified 15 types, including all types observed previously and two novel types, one of which has a non-canonical morphology and position. We validated the classification scheme and identified dozens of novel markers using methods that match molecular expression to cell morphology. This work provides a systematic methodology for achieving comprehensive molecular classification of neurons, identifies novel neuronal types, and uncovers transcriptional differences that distinguish types within a class

βIII-tubulin: a novel mediator of chemoresistance and metastases in pancreatic cancer

Pancreatic cancer is a leading cause of cancer-related deaths in Western societies. This poor prognosis is due to chemotherapeutic drug resistance and metastatic spread. Evidence suggests that microtubule proteins namely, β-tubulins are dysregulated in tumor cells and are involved in regulating chemosensitivity. However, the role of β-tubulins in pancreatic cancer are unknown. We measured the expression of different β-tubulin isotypes in pancreatic adenocarcinoma tissue and pancreatic cancer cells. Next, we used RNAi to silence βIII-tubulin expression in pancreatic cancer cells, and measured cell growth in the absence and presence of chemotherapeutic drugs. Finally, we assessed the role of βIII-tubulin in regulating tumor growth and metastases using an orthotopic pancreatic cancer mouse model. We found that βIII-tubulin is highly expressed in pancreatic adenocarcinoma tissue and pancreatic cancer cells. Further, we demonstrated that silencing βIII-tubulin expression reduced pancreatic cancer cell growth and tumorigenic potential in the absence and presence of chemotherapeutic drugs. Finally, we demonstrated that suppression of βIII-tubulin reduced tumor growth and metastases in vivo. Our novel data demonstrate that βIII-tubulin is a key player in promoting pancreatic cancer growth and survival, and silencing its expression may be a potential therapeutic strategy to increase the long-term survival of pancreatic cancer patients

Specific β-Tubulin Isotypes Can Functionally Enhance or Diminish Epothilone B Sensitivity in Non-Small Cell Lung Cancer Cells

Epothilones are a new class of microtubule stabilizing agents with promising preclinical and clinical activity. Their cellular target is β-tubulin and factors influencing intrinsic sensitivity to epothilones are not well understood. In this study, the functional significance of specific β-tubulin isotypes in intrinsic sensitivity to epothilone B was investigated using siRNA gene knockdown against βII-, βIII- or βIVb-tubulins in two independent non-small cell lung cancer (NSCLC) cell lines, NCI-H460 and Calu-6. Drug-treated clonogenic assays showed that sensitivity to epothilone B was not altered following knockdown of βII-tubulin in both NSCLC cell lines. In contrast, knockdown of βIII-tubulin significantly increased sensitivity to epothilone B. Interestingly, βIVb-tubulin knockdowns were significantly less sensitive to epothilone B, compared to mock- and control siRNA cells. Cell cycle analysis of βIII-tubulin knockdown cells showed a higher percentage of cell death with epothilone B concentrations as low as 0.5 nM. In contrast, βIVb-tubulin knockdown cells displayed a decrease in epothilone B-induced G2-M cell cycle accumulation compared to control siRNA cells. Importantly, βIII-tubulin knockdowns displayed a significant dose-dependent increase in the percentage of apoptotic cells upon treatment with epothilone B, as detected using caspase 3/7 activity and Annexin-V staining. Higher concentrations of epothilone B were required to induce apoptosis in the βIVb-tubulin knockdowns compared to control siRNA, highlighting a potential mechanism underlying decreased sensitivity to this agent. This study demonstrates that specific β-tubulin isotypes can influence sensitivity to epothilone B and may influence differential sensitivity to this promising new agent

Phenotypic screen for oxygen consumption rate identifies an anti-cancer naphthoquinone that induces mitochondrial oxidative stress.

A hallmark of cancer cells is their ability to reprogram nutrient metabolism. Thus, disruption to this phenotype is a potential avenue for anti-cancer therapy. Herein we used a phenotypic chemical library screening approach to identify molecules that disrupted nutrient metabolism (by increasing cellular oxygen consumption rate) and were toxic to cancer cells. From this screen we discovered a 1,4-Naphthoquinone (referred to as BH10) that is toxic to a broad range of cancer cell types. BH10 has improved cancer-selective toxicity compared to doxorubicin, 17-AAG, vitamin K3, and other known anti-cancer quinones. BH10 increases glucose oxidation via both mitochondrial and pentose phosphate pathways, decreases glycolysis, lowers GSH:GSSG and NAPDH/NAPD+ ratios exclusively in cancer cells, and induces necrosis. BH10 targets mitochondrial redox defence as evidenced by increased mitochondrial peroxiredoxin 3 oxidation and decreased mitochondrial aconitase activity, without changes in markers of cytosolic or nuclear damage. Over-expression of mitochondria-targeted catalase protects cells from BH10-mediated toxicity, while the thioredoxin reductase inhibitor auranofin synergistically enhances BH10-induced peroxiredoxin 3 oxidation and cytotoxicity. Overall, BH10 represents a 1,4-Naphthoquinone with an improved cancer-selective cytotoxicity profile via its mitochondrial specificity

Mapping Copy Number Variation by Population Scale Genome Sequencing

Genomic structural variants (SVs) are abundant in humans, differing from other forms of variation in extent, origin and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (that is, copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analysing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.Organismic and Evolutionary Biolog

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes

The Characterization of Twenty Sequenced Human Genomes

We present the analysis of twenty human genomes to evaluate the prospects for identifying rare functional variants that contribute to a phenotype of interest. We sequenced at high coverage ten “case” genomes from individuals with severe hemophilia A and ten “control” genomes. We summarize the number of genetic variants emerging from a study of this magnitude, and provide a proof of concept for the identification of rare and highly-penetrant functional variants by confirming that the cause of hemophilia A is easily recognizable in this data set. We also show that the number of novel single nucleotide variants (SNVs) discovered per genome seems to stabilize at about 144,000 new variants per genome, after the first 15 individuals have been sequenced. Finally, we find that, on average, each genome carries 165 homozygous protein-truncating or stop loss variants in genes representing a diverse set of pathways

Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude

Six new loci associated with body mass index highlight a neuronal influence on body weight regulation

Common variants at only two loci, FTO and MC4R, have been reproducibly associated with body mass index (BMI) in humans. To identify additional loci, we conducted meta-analysis of 15 genome-wide association studies for BMI (n > 32,000) and followed up top signals in 14 additional cohorts (n > 59,000). We strongly confirm FTO and MC4R and identify six additional loci (P < 5 × 10−8): TMEM18, KCTD15, GNPDA2, SH2B1, MTCH2 and NEGR1 (where a 45-kb deletion polymorphism is a candidate causal variant). Several of the likely causal genes are highly expressed or known to act in the central nervous system (CNS), emphasizing, as in rare monogenic forms of obesity, the role of the CNS in predisposition to obesity