OBJECTIVE: How hexanucleotide (GGGGCC) repeat expansions in C9ORF72 cause amyotrophic lateral sclerosis (ALS) remains poorly understood. Both gain- and loss-of-function mechanisms have been proposed. Evidence supporting these mechanisms in vivo is, however, incomplete. Here we determined the effect of C9orf72 loss-of-function in mice. METHODS: We generated and analyzed a conditional C9orf72 knockout mouse model. C9orf72(fl/fl) mice were crossed with Nestin-Cre mice to selectively remove C9orf72 from neurons and glial cells. Immunohistochemistry was performed to study motor neurons and neuromuscular integrity, as well as several pathological hallmarks of ALS, such as gliosis and TDP-43 mislocalization. In addition, motor function and survival were assessed. RESULTS: Neural-specific ablation of C9orf72 in conditional C9orf72 knockout mice resulted in significantly reduced body weight but did not induce motor neuron degeneration, defects in motor function, or altered survival. INTERPRETATION: Our data suggest that C9orf72 loss-of-function, by itself, is insufficient to cause motor neuron disease. These results may have important implications for the development of therapeutic strategies for C9orf72-associated ALS.
OBJECTIVE: How hexanucleotide (GGGGCC) repeat expansions in C9ORF72 cause amyotrophic lateral sclerosis (ALS) remains poorly understood. Both gain- and loss-of-function mechanisms have been proposed. Evidence supporting these mechanisms in vivo is, however, incomplete. Here we determined the effect of C9orf72 loss-of-function in mice. METHODS: We generated and analyzed a conditional C9orf72 knockout mouse model. C9orf72(fl/fl) mice were crossed with Nestin-Cre mice to selectively remove C9orf72 from neurons and glial cells. Immunohistochemistry was performed to study motor neurons and neuromuscular integrity, as well as several pathological hallmarks of ALS, such as gliosis and TDP-43 mislocalization. In addition, motor function and survival were assessed. RESULTS: Neural-specific ablation of C9orf72 in conditional C9orf72 knockout mice resulted in significantly reduced body weight but did not induce motor neuron degeneration, defects in motor function, or altered survival. INTERPRETATION: Our data suggest that C9orf72 loss-of-function, by itself, is insufficient to cause motor neuron disease. These results may have important implications for the development of therapeutic strategies for C9orf72-associated ALS.
A hexanucleotide repeat expansion (GGGGCC) in a noncoding region of C9ORF72 is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD; C9ALS/FTD).1, 2, 3 Three possible pathogenic mechanisms have been linked to C9ORF72 repeat expansions: sequestration of RNA‐binding proteins, toxicity mediated by dipeptides formed as a result of repeat‐associated non‐ATG translation, and haploinsufficiency.1, 2, 4, 5, 6 Experimental evidence supporting each of these mechanisms is accumulating.7 However, how these mechanisms cause ALS/FTD and their relevance in vivo remain poorly understood.The function of C9ORF72 is not known, but it has been suggested to play a role in protein trafficking.8, 9 Several observations identify C9ORF72haploinsufficiency as a hallmark of C9ALS/FTD. First, different studies report decreased C9ORF72 mRNA expression in brain tissue, lymphoblast cells, and induced pluripotent stem cell–derived neurons of patients.1, 10, 11 Second, one study shows decreased C9ORF72 protein expression in frontal cortex of C9ALS/FTDpatients.12 Third, knockdown of C9orf72 in model organisms such as Caenorhabditis elegans and zebrafish embryos causes motor deficits.10, 13 However, knockdown of C9orf72 caused by a single intracerebroventricular injection of antisense oligonucleotides (ASOs) in mice does not affect motor function or anxiety.14To test the haploinsufficiency model and to determine whether lack of C9orf72 expression leads to motor neuron degeneration or abnormal motor function, we generated and analyzed a conditional C9orf72 knockout mouse model.
Materials and Methods
Mouse Husbandry, Breeding, and Genotyping
All animal use and care were in accordance with local institution guidelines. Mice were kept on a 12‐hour light/dark cycle with food and water available ad libitum. B6;SJL‐Tg(ACTFLPe)9205Dym/J mice and B6.Cg‐Tg(Nes‐cre)1Kln/J mice were obtained from Jackson Laboratory (Bar Harbor, ME; 003800; 003771) and C57Bl/6J mice from Charles River Laboratories (Wilmington, MA). To generate C9orf72mice, a targeting construct was designed to insert an Frt‐flanked neomycin cassette and 1 loxP site upstream of exon 4 and 1 loxP site downstream of exon 5 of C9orf72. Deletion of exons 4 and 5 targets all reported C9orf72 isoforms in the mouse. This construct was electroporated into C57Bl/6 embryonic stem cells. Correctly targeted stem cells, as determined by polymerase chain reaction (PCR) and Southern blot analysis, were injected into blastocysts, and chimeric mice were bred with C57Bl/6J mice. The resulting C9orf72mice were then bred with mice expressing Flp recombinase in their germline to remove the Frt‐flanked neomycin cassette, generating C9orf72 offspring. Female C9orf72 or C9orf72mice were crossed with male Nestin‐Cre mice to generate neural‐specific C9orf72 conditional knockout mice. Mice were genotyped using primers to detect the Cre gene (forward = 5′‐GCGGTCTGGCAGTAAAAACTATC‐3′, reverse = 5‐′GTGAAACAGCATTGCTGTCACTT‐3′) and the genomic region containing the loxP sequences (forward = 5′‐CCACGGAGGGATGTTCTTTA‐3′, reverse = 5′‐GAAACCAGACCCAAACACAGA‐3′).
Antibody Generation
The anti‐C9orf72rabbit polyclonal antibody C9‐2034 was generated from an N‐terminal thioredoxin fusion of a stretch of 58 amino acids present in all humanC9orf72 isoforms comprising MEDQGQSIIPMLTGEVIPVMELLSSMKSHSVPEEIDIADTVLNDDDIGDSCHEGFLLK. This was generated from a C9orf72 short isoform expression construct using the primers 5′‐CCCGAATTCGAGAGAATGGAAGATCAGGGT‐3′ and 5′‐GAAGCGGCCGCATCTGCTTCATCCAGCTTTTATGA‐3′. The PCR product was digested and cloned into the EcoRI/NotI sites of pET32a (Clontech Laboratories, Mountain View, CA). Thioredoxin (Thx)‐C9orf72‐tail expression was induced in transformed BL21(DE3) Escherichia coli (Agilent Technologies, Santa Clara, CA) using 1mM Isopropyl β‐D‐1‐thiogalactopyranoside (IPTG) for 3 hours at 37°C. Thx‐C9orf72‐Short was purified using HisPur resin according to the manufacturer's instructions (Thermo Fisher Scientific, Waltham, MA) following solubilization in sonication buffer (20mM Tris pH 8.0, 100mM NaCl) using a Vibra‐Cell Ultrasonic Processor (Sonics & Materials, Newtown, CT). The anti‐C9orf72rabbit polyclonal antibody C9‐2074 was generated by reimmunization of the previously described N‐terminal thioredoxin fusion of full‐length C9orf72 short isoform.12 The respective purified fusion proteins were used as antigens for custom rabbit polyclonal antibody generation (Covalab, Villeurbanne, France). Custom antisera were immunoaffinity‐purified against their antigen following preabsorption against Thx and glial fibrillary acidic protein (Thx‐GFAP expressed in BL21[DE3] transformed with GFAP pET32a). Both antibodies detected a 50 to 55kDa protein in mouse and human brain lysates, as well as myc‐tagged C9orf72 from transfected HEK293T cells.The anti‐C9orf72rabbit polyclonal antibodies were validated for Western blotting using siRNA‐mediated knockdown of endogenous C9orf72 in HEK293T cells using 27‐mer oligonucleotide duplexes as previously described.12 Additional validation in SH‐SY5Y neuroblastoma cells was performed following knockdown of endogenous C9orf72 using Dicer substrate RNAi interference using custom dicer substrate siRNA (DsiRNA; Integrated DNA Technologies, Coralville, IA). Custom DsiRNA oligonucleotide duplex sequences are C9ORF72‐DsiRNA1F 5‐GGAAAGAAUAUGGAUGCAUAAGGAAA‐3′/C9ORF72‐DsiRNA1R 5′‐UUUCCUUAUGCAUCCAUAUUCUUCCUU‐3′ and C9ORF72‐DsiRNA2F 5′‐GUACUCAAUGAUGAUGAUAUUGGTG‐3′/C9ORF72‐DsiRNA2R 5′‐CACCAAUAUCAUCAUCAUUGAGUACUG‐3′. The inventoried FLuc‐S1 exogenous reporter gene DsiRNA duplex was used as a nontargeting control. DsiRNA was transfected using the same conditions as 27mer oligonucleotide duplexes. Whole cell and tissue lysate preparation and Western blotting were performed as described previously.12 Western blotting using the anti‐C9orf72rabbit polyclonal antibodies detected no C9orf72 expression in knockout mice (Fig 1) or reduced expression following knockdown of C9orf72 (in HEK293 or SHSY5Y cells; data not shown).
Figure 1
Expression of C9orf72 and generation of C9orf72 conditional knockout mice. (A–G) In situ hybridization shows expression of C9orf72 in coronal sections of the adult brain and spinal cord. (F, G) Double staining combining in situ hybridization for C9orf72 and immunohistochemistry for SMI132, a marker for motor neurons (indicated by arrowheads). (E, M) Sections incubated with sense probes do not display specific staining. (H–O) At embryonic and early postnatal stages, C9orf72 is expressed in spinal cord motor neurons (arrowheads in H–J), the retina, olfactory system, and other brain regions, in addition to non‐neuronal tissues such as kidney. (P) Schematic representation of the C9orf72 targeting vector, targeted allele, floxed allele following Flp recombination, and deleted allele following Cre recombination. Black arrows indicate transcription start site, gray arrows indicate location of loxP primers (P1, P2), and numbers indicate exons. Please note that exons 4 and 5 are shared by all 3 reported mouse C9orf72 isoforms. (Q) Genotyping using loxP primers on DNA isolated from ear and brain tissue of Nestin‐Cre and Nestin‐Cre mice shows brain‐specific recombination. (R) Quantitative real‐time polymerase chain reaction shows tissue‐specific reduction of C9orf72 levels in Nestin‐Cre and/or Nestin‐Cre in brain (0.493, p = 0.0004 and 0.044, p = 0.00007) and spinal cord (0.684, p = 0.083 and 0.145, p = 0.002), but not in muscle (0.789, p = 0.594 and 0.585, p = 0.169) or liver (1.94, p = 0.219 and 1.74, p = 0.473). Results are shown with Tbp as a housekeeping gene. Data are presented as mean ± standard deviation, n = 3 for each genotype. (S) In situ hybridization for C9orf72 shows loss of C9orf72 expression in spinal cord of Nestin‐Cre mice. **p < 0.01, ***p < 0.001. CA1 = cornu ammonis 1; Cx = cortex; DG = dentate gyrus; DRG = dorsal root ganglia; F = Frt site; GL = granular cell layer; GM = gray matter; L = loxP site; Mol = molecular layer; N = neomycin cassette; OE = olfactory epithelium; P1 = forward loxP primer; P2 = reverse loxP primer; Ret = retina; RN = red nucleus; VNO = vomeronasal organ; WM = white matter; WT = wild type. Scale bars: A, 200μm; H, 100μm; K, 100μm, S, 150μm.
Expression of C9orf72 and generation of C9orf72 conditional knockout mice. (A–G) In situ hybridization shows expression of C9orf72 in coronal sections of the adult brain and spinal cord. (F, G) Double staining combining in situ hybridization for C9orf72 and immunohistochemistry for SMI132, a marker for motor neurons (indicated by arrowheads). (E, M) Sections incubated with sense probes do not display specific staining. (H–O) At embryonic and early postnatal stages, C9orf72 is expressed in spinal cord motor neurons (arrowheads in H–J), the retina, olfactory system, and other brain regions, in addition to non‐neuronal tissues such as kidney. (P) Schematic representation of the C9orf72 targeting vector, targeted allele, floxed allele following Flp recombination, and deleted allele following Cre recombination. Black arrows indicate transcription start site, gray arrows indicate location of loxP primers (P1, P2), and numbers indicate exons. Please note that exons 4 and 5 are shared by all 3 reported mouseC9orf72 isoforms. (Q) Genotyping using loxP primers on DNA isolated from ear and brain tissue of Nestin‐Cre and Nestin‐Cre mice shows brain‐specific recombination. (R) Quantitative real‐time polymerase chain reaction shows tissue‐specific reduction of C9orf72 levels in Nestin‐Cre and/or Nestin‐Cre in brain (0.493, p = 0.0004 and 0.044, p = 0.00007) and spinal cord (0.684, p = 0.083 and 0.145, p = 0.002), but not in muscle (0.789, p = 0.594 and 0.585, p = 0.169) or liver (1.94, p = 0.219 and 1.74, p = 0.473). Results are shown with Tbp as a housekeeping gene. Data are presented as mean ± standard deviation, n = 3 for each genotype. (S) In situ hybridization for C9orf72 shows loss of C9orf72 expression in spinal cord of Nestin‐Cre mice. **p < 0.01, ***p < 0.001. CA1 = cornu ammonis 1; Cx = cortex; DG = dentate gyrus; DRG = dorsal root ganglia; F = Frt site; GL = granular cell layer; GM = gray matter; L = loxP site; Mol = molecular layer; N = neomycin cassette; OE = olfactory epithelium; P1 = forward loxP primer; P2 = reverse loxP primer; Ret = retina; RN = red nucleus; VNO = vomeronasal organ; WM = white matter; WT = wild type. Scale bars: A, 200μm; H, 100μm; K, 100μm, S, 150μm.
In Situ Hybridization
Fresh frozen mouse embryos, or spinal cord and brain tissue from adult mice were cut at 16 to 20µm on a cryostat. Nonradioactive in situ hybridization was performed as described previously.15 A digoxigenin‐labeled C9orf72 probe was transcribed from mouse cDNA (a 736bp fragment corresponding to nucleotides 298–1,034 of 3110043O21Rik) using a DIG RNA labeling SP6/T7 kit (Roche, Basel, Switzerland). Sense probes were used to confirm specificity. Images were taken on a Zeiss (Oberkochen, Germany) AxioSkop2 microscope.
Quantitative Real‐Time PCR
Fresh frozen tissue was homogenized in TRIzol reagent and total RNA was isolated. After DNAse treatment, cDNA was synthesized using superscript II reverse transcriptase (Life Technologies, Carlsbad, CA) and oligo‐dT primers. Quantitative real‐time PCR was performed using SYBR green reagent (Roche) on a QuantStudio 6 Flex real‐time PCR system (Applied Biosystems, Foster City, CA) using primers against 3110043O21Rik (Q1 [exon 2]: forward = 5′‐GACCTGGAATGCAGTGAGAC‐3′, reverse = 5′‐TCCCAGTAAGCAAAGGTAGC; Q2 [exons 3–6]: forward = 5′‐GGCACAGAGAGGATGGAAGA‐3′, reverse = 5′‐CTGCCAACTACAACGGAACA‐3′; Q3 [exons 6–8]: forward = 5′‐TGGCTGTTCCGTTGTAGTTG‐3′, reverse = 5′‐AACTGCCTGTTGCATCCTTT‐3′; Q4 [exons 7–8]: FW = 5′‐GTACGAATCGGGACTCTTTG‐3′, reverse = 5′‐GCTTGACAGTGTTGACATCC‐3′; Q5 [exons 8–9]: forward = 5′‐CAACGCAGATACATGAGGTC‐3′, reverse = 5′‐GGAAGGCTTTCACTAGAGTGTC‐3′; Q6 [exons 10–11]: forward = 5′‐CCCTTTAAGTCTCTTCGGAACC‐3′, reverse = 5′‐AGACAGACACAGCAAACCAC‐3′) and mRNA levels were normalized to Gapdh (forward = 5′‐GAGACGGCCGCATCTTCTTGT‐3′, reverse = 5‐′ CACACCGACCTTCACCATTTT‐3′) and Tbp (QT00198443; Qiagen, Valencia, CA) housekeeping genes. Samples from 3 mice per genotype were amplified in triplicate, and relative quantification using comparative ΔΔCT values was statically analyzed using Student t test.
Western Blotting
Western blotting was performed as described previously.16 In brief, fresh frozen tissue was homogenized in radioimmunoprecipitation assay buffer (50mM Tris pH 7.6, 150mM NaCl, 0.5mM ethylenediaminetetraacetic acid, 1% sodium deoxycholate, 0.1% sodium dodecyl sulfate, 1% Triton‐X‐100, and protease inhibitor cocktail), centrifuged twice at 10,000 × g, and the supernatant was used as total protein. Western blotting was performed using 50µg of total protein per lane as measured by bicinchoninic acid assay (Pierce Biotechnology, Rockford, IL), and blots were probed with the following antibodies: anti‐C9orf72 (2034; 1:250), anti‐C9orf72 (2074; 1:250), and anti–β‐actin (Sigma, St Louis, MO; 1:2,000). After overnight incubation, blots were probed with horseradish peroxidase–conjugated secondary antibodies and developed with SuperSignal West Femto chemiluminescent substrate (Pierce Biotechnology).
Motor Function Testing
Mice were weighed and tested for their hindlimb extension reflex weekly. Every 3 weeks, mice performed 5 trials on an accelerating rotarod (Ugo Basile, Varese, Italy) at 4rpm with an acceleration of 0.2rpm/s for 180 seconds and the latency to fall was recorded. Mice that rotated passively were scored as fallen. Grip strength of front paws and of front and hind paws together was measured every 3 months using a grip strength meter (Bioseb, Vitrolles, France). Strength was measured 5 times for each mouse both for front and for front and hind paws. All behavioral tests were performed on ≥11 mice per genotype for each time point. Statistical analysis was carried out using R software (http://www.r‐project.org). (Non‐)linear mixed models were used to account for repeated measurements and other sources of correlation within the longitudinal data and to provide a correct model in case of missing data. For statistical analysis of body weight, a nonlinear mixed effects model was used with gender as a covariate. Interactions between time and behavioral tests were analyzed to study possible changes during follow‐up. For rotarod and grip strength statistical analysis, a nonlinear mixed effects model was used with gender and body weight as covariates. Survival analysis was carried out using a Cox regression analysis with mean body weight as covariate. Cox proportional hazard assumptions were tested.
Immunohistochemistry
Immunohistochemistry was performed as described previously.17 Tissues were removed and cryoprotected in 30% sucrose after transcardial perfusion with 4% paraformaldehyde. Spinal cord and brain were cut at 20µm and mounted on SuperFrost Plus slides (Thermo Fisher Scientific). For immunofluorescent stainings, sections were blocked and incubated overnight at 4°C with the following antibodies: anti‐GFAP (Dako, Carpinteria, CA), anti‐Iba1 (Wako Chemicals, Richmond, VA), anti–choline acetyltransferase (ChAT; Millipore, Billerica, MA), and anti–TDP‐43 (Proteintech, Chicago, IL). Alexa Fluor–conjugated secondary antibodies were incubated for 2 hours at room temperature (RT) and 4′,6‐diamidino‐2‐phenylindole was used to stain nuclei. For diaminobenzidine (DAB) stainings, sections were incubated overnight at 4°C with anti‐ubiquitin antibodies (Millipore) and for 1 h at RT with goat anti‐mousebiotin. Slides were incubated with ABC reagent (Vector Laboratories, Burlingame, CA) for 1 h and developed with DAB reagent (Sigma). For motor neuron counting, at least 10 spinal cord sections (L1–L6) per mouse using 3 mice per genotype were analyzed and only ChAT‐positive neurons with a clear nucleus were counted and measured. For GFAP and Iba1 quantification, at least 10 spinal cord sections (L1–L6) from 3 mice per genotype were analyzed for positive staining area using ImageJ (NIH, Bethesda, MD). For neuromuscular junction staining, gastrocnemius muscle was cut longitudinally at 40µm and stained with anti–α‐bungarotoxin‐488 (Life Technologies), anti‐sv2 (Developmental Studies Hybridoma Bank [DSHB], Iowa City, IA), and anti‐neurofilament (DSHB) antibodies. At least 150 neuromuscular junctions were counted per mouse. Images were taken on an Olympus (Tokyo, Japan) FluoView FV1000 confocal microscope or a Zeiss AxioScope microscope followed by image analysis using ImageJ. Statistical analysis was performed using Student t test (SPSS Statistics; IBM, Armonk, NY).
Results
Genetic Ablation of C9orf72 in Mice
Previous work analyzing C9orf72‐LacZ knockin mice reported C9orf72 expression in neurons known to degenerate in ALS and FTD, such as cortical and spinal cord motor neurons.18 To confirm and extend these findings, we performed in situ hybridization for the mouseC9orf72 orthologue, 3110043O21Rik (here referred to as C9orf72), at embryonic, postnatal, and adult stages. As reported previously, neuronal C9orf72 expression was detected in several regions of the adult brain and spinal cord (see Fig 1; Table1).18 In contrast to previous work, however, prominent expression of C9orf72 was also observed in embryonic and early postnatal neurons. These included retinal ganglion cells, sensory neurons in the olfactory epithelium and in dorsal root ganglia, and spinal motor neurons. In addition, expression in non‐neuronal tissues such as kidney and tooth was detected. These results indicate that C9orf72 may also function during embryonic and postnatal development.
Table 1
Overview of C9orf72Expression in the Developing Mouse Embryo and Adult Central Nervous System
Location
E13.5
E15.5
E17.5
P1
P15
Adult
Brain
Cortex
+
+
+
+
++
+
Hippocampal formation
+
+
++
++
++
++
Dentate gyrus
N/D
N/D
++
++
++
++
Cornu ammonis
N/D
N/D
++
++
++
++
Thalamus
N/D
N/D
N/D
+
++
+
Midbrain/pons
N/D
+
+
+
++
++
Ventral tegmental area
N/D
N/D
N/D
−
−
−
Pontine nuclei
N/D
N/D
N/D
+
++
++
Red nucleus
N/D
N/D
N/D
+
++
++
Facial motor nucleus
N/D
N/D
N/D
+
++
++
Cerebellum
+
+
+
+
+
+
Molecular layer
N/D
N/D
N/D
+
+
+
Granule cell layer
N/D
N/D
N/D
+
+
+
Spinal cord
Gray matter
++
++
++
++
++
++
White matter
−
−
−
−
−
−
Dorsal root ganglia
++
++
++
++
++
++
Other
Eye
++
++
++
++
N/D
N/D
Olfactory epithelium
++
++
++
++
N/D
N/D
Tooth
++
++
++
++
N/D
N/D
Kidney
+
+
+
+
N/D
N/D
N/D = not determined, − = no expression, + = weak expression, ++ = moderate/strong expression.
Overview of C9orf72Expression in the Developing Mouse Embryo and Adult Central Nervous SystemN/D = not determined, − = no expression, + = weak expression, ++ = moderate/strong expression.To circumvent problems such as embryonic lethality, we generated conditional C9orf72 knockout mice using the Cre‐loxP system and bred the resulting C9orf72mice to Nestin‐Cre mice (see Fig 1). In Nestin‐Cre mice, Cre recombination occurs in neuronal and glial precursors from E10.5 onward.19 The specific ablation of C9orf72 mRNA and protein in Nestin‐Cre mice was confirmed using different approaches. First, DNA genotyping on ear clips or brain tissue showed excision of exons 4 and 5 only in brain tissue. Second, quantitative real‐time PCR detected a reduction in C9orf72 mRNA levels in brain and spinal cord tissue but not in other tissues of Nestin‐Cre mice. Third, in situ hybridization detected C9orf72 in the spinal cord of Nestin‐Cre but not in Nestin‐Cre mice.Three C9orf72 isoforms have been reported in mice. These isoforms all contain exons 4 and 5 and are targeted in the C9orf72 conditional knockout mouse (Fig 2). To confirm that the expression of all isoforms is affected by the excision of exons 4 and 5, we performed quantitative real‐time PCR to detect expression of different exonic regions present in 1 or multiple isoforms. As expected, a reduction of C9orf72 mRNA levels was found for all exonic regions analyzed. In addition, we generated a new anti‐C9orf72 antibody. This antibody was raised against a 58 amino acid stretch shared by all isoforms (C9‐2034). The other antibody used was reported previously and was raised using a large part of the C9orf72 protein. The region used is present in 2 of the 3 isoforms (C9‐2074).12 Western blotting using these antibodies showed a predicted ∼55kDa band, corresponding to isoform 1, that was lost in brain tissue of Nestin‐Cre mice. The other 2 C9orf72 isoforms (isoforms 2 and 3) were not detected in wild‐type mouse brain tissue (data not shown). Together, these results confirm the neural‐specific and successful ablation of C9orf72 in Nestin‐Cre mice.
Figure 2
Cre‐mediated excision of C9orf72 affects all C9orf72 isoforms. (A) Schematic overview of the 3 reported mouse C9orf72 isoforms. Q1 to Q6 indicate regions amplified by quantitative real‐time polymerase chain reaction (qPCR). The in situ probes used are indicated. Dotted lines indicate antigen regions used to generate the 2 anti‐C9orf72 antibodies (C9‐2034 and C9‐2074). (B) Amino acid sequence of the antigen used to generate the novel anti‐C9orf72 antibody C9‐2034, and its comparison to human C9ORF72 isoforms and mouse C9orf72 isoforms. (C) qPCR using primers that target different exonic regions shared by 1 or multiple mouse isoforms shows a significant reduction of brain C9orf72 mRNA levels in Nestin‐Cre and/or Nestin‐Cre mice. Results are shown using Tbp as a housekeeping gene. Data are presented as mean ± standard deviation, n = 2 for each genotype. (D) Western blot analysis of whole brain lysates using 2 different anti‐C9orf72 antibodies shows a ∼55kDa band in control mice that is reduced or absent in Nestin‐Cre or Nestin‐Cre mice, respectively. Note that only isoform 1 is detected at the protein level in mouse brain. β‐Actin was used as a loading control. ***p < 0.001. L = loxP sites.
Cre‐mediated excision of C9orf72 affects all C9orf72 isoforms. (A) Schematic overview of the 3 reported mouseC9orf72 isoforms. Q1 to Q6 indicate regions amplified by quantitative real‐time polymerase chain reaction (qPCR). The in situ probes used are indicated. Dotted lines indicate antigen regions used to generate the 2 anti‐C9orf72 antibodies (C9‐2034 and C9‐2074). (B) Amino acid sequence of the antigen used to generate the novel anti‐C9orf72 antibody C9‐2034, and its comparison to humanC9ORF72 isoforms and mouseC9orf72 isoforms. (C) qPCR using primers that target different exonic regions shared by 1 or multiple mouse isoforms shows a significant reduction of brain C9orf72 mRNA levels in Nestin‐Cre and/or Nestin‐Cre mice. Results are shown using Tbp as a housekeeping gene. Data are presented as mean ± standard deviation, n = 2 for each genotype. (D) Western blot analysis of whole brain lysates using 2 different anti‐C9orf72 antibodies shows a ∼55kDa band in control mice that is reduced or absent in Nestin‐Cre or Nestin‐Cre mice, respectively. Note that only isoform 1 is detected at the protein level in mouse brain. β‐Actin was used as a loading control. ***p < 0.001. L = loxP sites.
Mice Lacking C9orf72 in Neurons and Glial Cells Do Not Display Motor Neuron Degeneration or Defects in Motor Function
ALS is characterized by progressive motor neuron degeneration and a decline in motor performance. To examine a possible role for C9orf72haploinsufficiency in ALS pathogenesis, we therefore assessed motor function, motor neuron number, and neuromuscular integrity in 12‐month‐old (data not shown) and 18‐month‐old Nestin‐Cre mice and littermate controls. Nestin‐Cre mice were born at the expected Mendelian ratio and were indistinguishable from Nestin‐Cre control littermates in viability, appearance, fertility, and the gross anatomical and histological appearance of the brain. Immunohistochemistry for ChAT, to visualize spinal motor neurons, did not reveal changes in neuron number or size (Fig 3A, B). Furthermore, innervation of the gastrocnemius muscle by motor axons, as assessed by overlap between the pre‐ and postsynaptic compartment of the neuromuscular junction, was similar in Nestin‐Cre mice and controls (see Fig 3C, D). Astro‐ and microgliosis are hallmarks of ALS. However, no differences in immunostaining for GFAP or Iba1, to label astrocytes or microglia, respectively, were detected between spinal cord or brain tissue from Nestin‐Cre mice and littermate controls (Figs 4A–D, 5A–F″, 6A–F″). Similarly, other hallmarks of ALS pathogenesis, such as TDP‐43 mislocalization or changes in ubiquitin staining, were absent in Nestin‐Cre mice (see Figs 4E–H, 7). Similar results were obtained in 12‐month‐old Nestin‐Cre mice, that is, absence of defects in motor neuron number, size, and connectivity and no signs of gliosis or other pathological hallmarks of ALS (data not shown). Thus, C9orf72deficiency in mice fails to cause several of the pathological hallmarks reported in ALSpatients or mouse models.
Figure 3
Neural‐specific ablation of C9orf72 does not cause changes in motor neuron number, size, or neuromuscular junction (NMJ) connectivity in 18‐month‐old Nestin‐Cre mice. (A, B) Immunofluorescent staining of adult spinal cord using anti–choline acetyltransferase (ChAT) antibodies. Graphs show quantification of motor neuron number and size. Dotted lines indicate border between gray matter (GM) and white matter (WM). (C, D) Representative images of NMJs in gastrocnemius muscle visualized by immunohistochemistry for neurofilament (NF) and synaptic vesicle protein 2 (SV2) in combination with Alexa Fluor 488–conjugated bungarotoxin (BTX) binding. Graph shows overlap between presynaptic (gray arrowhead) and postsynaptic staining (white arrowhead). Motor neurons and their NMJs are intact in Nestin‐Cre mice. Data in graphs represent mean ± standard deviation, n = 3 per genotype. Scale bars: A, 100μm; C, 30μm. DAPI = 4′,6‐diamidino‐2‐phenylindole.
Figure 4
Neural‐specific ablation of C9orf72 does not cause signs of neuroinflammation or other pathological hallmarks of amyotrophic lateral sclerosis in 18‐month‐old Nestin‐Cre mice. (A–D) Immunofluorescent staining of adult spinal cord using anti–glial fibrillary acidic protein (GFAP) or anti‐Iba1 antibodies. Graphs show quantification of GFAP‐ or Iba1‐positive areas. No signs of gliosis are observed in Nestin‐Cre mice. (E, F) Immunofluorescent staining of adult spinal cord using anti–TDP‐43 antibodies. Nuclear TDP‐43 staining is detected in Nestin‐Cre mice and control littermates. (G, H) Representative images showing immunostaining for ubiquitin in adult spinal cord. Ubiquitin expression is not affected by loss of C9orf72. Data in graphs represent mean ± standard deviation, n = 3 per genotype. Scale bars: A, C, G, 100μm; E, 30μm. DAPI = 4′,6‐diamidino‐2‐phenylindole; GM = gray matter; WM = white matter.
Figure 5
Immunostaining for glial fibrillary acidic protein (GFAP) in the brain of 18‐month‐old mice. Immunofluorescent staining is shown for GFAP in frontal cortex, motor cortex, and hippocampus of Nestin‐Cre mice (A–C″) and Nestin‐Cre mice (D–F″). Boxes in A–C and D–F are shown at higher magnification in A′–C″ and D′–F″, respectively. No changes in the distribution, number, or appearance of GFAP‐positive astrocytes were detected in Nestin‐Cre mice. Scale bars: A, D, 200μm; A′, D′, 100μm. DAPI = 4′,6‐diamidino‐2‐phenylindole.
Figure 7
Immunostaining for TDP‐43 and ubiquitin in several brain regions of 18‐month‐old mice. (A–C, G–I) Immunofluorescent staining of adult spinal cord using anti–TDP‐43 antibodies is shown. Nuclear TDP‐43 staining is detected in Nestin‐Cre mice and control littermates. (D–F, J–L) Representative images show immunostaining for ubiquitin in adult spinal cord. Ubiquitin expression is not affected by loss of C9orf72. DAPI = 4′,6‐diamidino‐2‐phenylindole; GCL = granule cell layer; H = hilus. Scale bars = 60μm.
Neural‐specific ablation of C9orf72 does not cause changes in motor neuron number, size, or neuromuscular junction (NMJ) connectivity in 18‐month‐old Nestin‐Cre mice. (A, B) Immunofluorescent staining of adult spinal cord using anti–choline acetyltransferase (ChAT) antibodies. Graphs show quantification of motor neuron number and size. Dotted lines indicate border between gray matter (GM) and white matter (WM). (C, D) Representative images of NMJs in gastrocnemius muscle visualized by immunohistochemistry for neurofilament (NF) and synaptic vesicle protein 2 (SV2) in combination with Alexa Fluor 488–conjugated bungarotoxin (BTX) binding. Graph shows overlap between presynaptic (gray arrowhead) and postsynaptic staining (white arrowhead). Motor neurons and their NMJs are intact in Nestin‐Cre mice. Data in graphs represent mean ± standard deviation, n = 3 per genotype. Scale bars: A, 100μm; C, 30μm. DAPI = 4′,6‐diamidino‐2‐phenylindole.Neural‐specific ablation of C9orf72 does not cause signs of neuroinflammation or other pathological hallmarks of amyotrophic lateral sclerosis in 18‐month‐old Nestin‐Cre mice. (A–D) Immunofluorescent staining of adult spinal cord using anti–glial fibrillary acidic protein (GFAP) or anti‐Iba1 antibodies. Graphs show quantification of GFAP‐ or Iba1‐positive areas. No signs of gliosis are observed in Nestin‐Cre mice. (E, F) Immunofluorescent staining of adult spinal cord using anti–TDP‐43 antibodies. Nuclear TDP‐43 staining is detected in Nestin‐Cre mice and control littermates. (G, H) Representative images showing immunostaining for ubiquitin in adult spinal cord. Ubiquitin expression is not affected by loss of C9orf72. Data in graphs represent mean ± standard deviation, n = 3 per genotype. Scale bars: A, C, G, 100μm; E, 30μm. DAPI = 4′,6‐diamidino‐2‐phenylindole; GM = gray matter; WM = white matter.Immunostaining for glial fibrillary acidic protein (GFAP) in the brain of 18‐month‐old mice. Immunofluorescent staining is shown for GFAP in frontal cortex, motor cortex, and hippocampus of Nestin‐Cre mice (A–C″) and Nestin‐Cre mice (D–F″). Boxes in A–C and D–F are shown at higher magnification in A′–C″ and D′–F″, respectively. No changes in the distribution, number, or appearance of GFAP‐positive astrocytes were detected in Nestin‐Cre mice. Scale bars: A, D, 200μm; A′, D′, 100μm. DAPI = 4′,6‐diamidino‐2‐phenylindole.Immunostaining for Iba1 in the brain of 18‐month‐old mice. Immunofluorescent staining is shown for Iba1 in frontal cortex, motor cortex, and hippocampus of Nestin‐Cre mice (A–C″) and Nestin‐Cre mice (D–F″). Boxes in A–C and D–F are shown at higher magnification in A′–C″ and D′–F″, respectively. No changes in the distribution, number, or appearance of Iba1‐positive microglia are detected in Nestin‐Cre mice. Scale bars: A,D, 200μm; A′, D′, 100μm. DAPI = 4′,6‐diamidino‐2‐phenylindole.Immunostaining for TDP‐43 and ubiquitin in several brain regions of 18‐month‐old mice. (A–C, G–I) Immunofluorescent staining of adult spinal cord using anti–TDP‐43 antibodies is shown. Nuclear TDP‐43 staining is detected in Nestin‐Cre mice and control littermates. (D–F, J–L) Representative images show immunostaining for ubiquitin in adult spinal cord. Ubiquitin expression is not affected by loss of C9orf72. DAPI = 4′,6‐diamidino‐2‐phenylindole; GCL = granule cell layer; H = hilus. Scale bars = 60μm.A reduction in body weight was detected in all Nestin‐Cre–positive mice, as has been previously reported for other studies using this Cre driver.20 However, even when considering this effect, Nestin‐Cre mice displayed a small but significant decrease in body weight, as compared to Nestin‐Cre control mice (−6.02%, p = 0.019). This difference was detected at the first time point of measurements and did not significantly change in time (Fig 8). To assess motor performance, the hindlimb extension reflex was tested (data not shown) and mice were assessed using an accelerating rotarod and a grip strength meter for a period of 18 months after birth. No differences in motor performance or grip strength were found between Nestin‐Cre mice and various control mice. A similar result was obtained in a smaller cohort of 2‐year‐old mice (n ≥ 4; data not shown). Finally, survival was not significantly altered in Nestin‐Cre mice, the oldest mice living >24 months. Thus, Nestin‐Cre mice do not display overt defects in motor function or decreased survival.
Figure 8
Loss of C9orf72 in mice does not affect motor function and survival. (A) A small decrease in body weight is detected in Nestin‐Cre mice as compared to Nestin‐Cre controls (−6.02%, p = 0.019, linear mixed‐effect model for repeated measurements [LME]; n ≥ 11 for each genotype at each time point). (B–D) Rotarod performance and grip strength is similar in Nestin‐Cre, Nestin‐Cre, and Nestin‐Cre mice (rotarod, p = 0.084 and p = 0.899; grip strength front, p = 0.346 and p = 0.295; front and hind, p = 0.731 and p = 0.305, LME). Data are shown as mean ± standard deviation, n ≥ 11 for each genotype at each time point. (E) Survival rates as shown by a Kaplan–Meier curve. Survival of Nestin‐Cre mice is similar to control mice.
Loss of C9orf72 in mice does not affect motor function and survival. (A) A small decrease in body weight is detected in Nestin‐Cre mice as compared to Nestin‐Cre controls (−6.02%, p = 0.019, linear mixed‐effect model for repeated measurements [LME]; n ≥ 11 for each genotype at each time point). (B–D) Rotarod performance and grip strength is similar in Nestin‐Cre, Nestin‐Cre, and Nestin‐Cre mice (rotarod, p = 0.084 and p = 0.899; grip strength front, p = 0.346 and p = 0.295; front and hind, p = 0.731 and p = 0.305, LME). Data are shown as mean ± standard deviation, n ≥ 11 for each genotype at each time point. (E) Survival rates as shown by a Kaplan–Meier curve. Survival of Nestin‐Cre mice is similar to control mice.
Discussion
C9orf72haploinsufficiency may contribute to the pathogenic mechanism in C9ALS/FTDpatients. To test this hypothesis, we generated a conditional C9orf72 knockout mouse model. Two previous studies have reported the generation of C9orf72 knockout mice, but neither study provided a histological or behavioral characterization.18, 21 In one of these studies, using C9orf72‐LacZ knockin reporter mice, C9orf72 expression was not detected at embryonic and early postnatal stages in the mouse.18 In contrast, we show neuronal C9orf72 expression at several sites in the embryonic nervous system, as well as in non‐neuronal tissues. This observation is in line with a recent study showing expression of C9orf72 both inside and outside the nervous system of zebrafish embryos.10 The cause of this apparent discrepancy is unknown but may include the inability of heterozygous LacZ reporter mice to detect low levels of gene expression.22 Together, our data suggest that C9orf72 plays a role during embryonic and postnatal development.The expression of C9orf72 at embryonic stages and in several tissues outside the nervous system prompted us to selectively ablate C9orf72 in neurons and glial cells by crossing C9orf72mice with Nestin‐Cre mice. The resulting Nestin‐Cre mice did not show overt defects in motor function or pathological hallmarks of ALS, such as reduced motor neuron number, gliosis, TDP‐43 mislocalization, or enhanced ubiquitination. A small but significant decrease in body weight was found in Nestin‐Cre mice. Our results are in apparent contrast with recent work in model organisms such as zebrafish and C. elegans. Knockdown or knockout, respectively, of C9orf72 orthologues in zebrafish and C. elegans causes motor deficits. These include age‐dependent motility defects and γ‐aminobutyric acidergic motor neuron degeneration in C. elegans, and motor axon and behavior deficits in zebrafish.10, 13 The discrepancy with our results may be caused by the low homology of C9orf72 orthologues in humans and in zebrafish and C. elegans (76% and 23% compared to humanC9ORF72, respectively) as compared to the high homology between human and mouseC9orf72 (98%). Another possible explanation could be redundancy of the C9orf72 gene in mammalian species. In line with our results, a single treatment of wild‐type mice with ASOs directed against C9orf72 did not lead to motor abnormalities.14Together, our data suggest that loss of C9orf72 on its own is insufficient to cause ALS. Our results, however, do not rule out the possibility that C9ORF72 loss‐of‐function modulates the disease process and influences disease onset, severity, and duration in C9ALS/FTDpatients. The disease mechanisms caused by repeat expansions in C9ORF72 are likely to be complex and may ultimately result from interplay between C9ORF72 loss‐ and gain‐of‐function mechanisms.23, 24 In support of a gain‐of‐function mechanism, recent work shows that AAV‐mediated expression of 66 GGGGCC repeats in the adult mouse brain induces key features of ALS/FTD, such as TDP‐43 pathology and behavioral deficits.6 A different study, using transgenic mice with inducible expression of hexanucleotide (GGGGCC) repeats, found ubiquitinated inclusions, but these inclusions were not detected in neuronal populations relevant for ALS/FTD and no motor phenotype was found.25 Our finding that loss of C9orf72 function does not lead to motor neuron disease in mice has consequences for therapeutic strategies, because it indicates that ASOs directed at the repeat expansion that also lower C9ORF72 expression are unlikely to have negative side effects due to reduced C9ORF72 expression.
Authorship
M.K., A.M.B., H.‐J.W., M.L.T., C.A.C.Z., R.V.d.S., and R.D.S. performed experiments and data analysis. A.J.W. and D.J.B. provided unpublished reagents. M.K., A.M.B., H.‐J.W., J.H.V., L.H.v.d.B., and R.J.P. designed and coordinated the study. M.K. and R.J.P. prepared the manuscript with input from all coauthors. L.H.v.d.B. and R.J.P. are joint senior authors.
Authors: Clotilde Lagier-Tourenne; Michael Baughn; Frank Rigo; Shuying Sun; Patrick Liu; Hai-Ri Li; Jie Jiang; Andrew T Watt; Seung Chun; Melanie Katz; Jinsong Qiu; Ying Sun; Shuo-Chien Ling; Qiang Zhu; Magdalini Polymenidou; Kevin Drenner; Jonathan W Artates; Melissa McAlonis-Downes; Sebastian Markmiller; Kasey R Hutt; Donald P Pizzo; Janet Cady; Matthew B Harms; Robert H Baloh; Scott R Vandenberg; Gene W Yeo; Xiang-Dong Fu; C Frank Bennett; Don W Cleveland; John Ravits Journal: Proc Natl Acad Sci U S A Date: 2013-10-29 Impact factor: 11.205
Authors: Peter E A Ash; Kevin F Bieniek; Tania F Gendron; Thomas Caulfield; Wen-Lang Lin; Mariely Dejesus-Hernandez; Marka M van Blitterswijk; Karen Jansen-West; Joseph W Paul; Rosa Rademakers; Kevin B Boylan; Dennis W Dickson; Leonard Petrucelli Journal: Neuron Date: 2013-02-12 Impact factor: 17.173
Authors: Ewout J N Groen; Katsumi Fumoto; Anna M Blokhuis; Jooyeon Engelen-Lee; Yeping Zhou; Dianne M A van den Heuvel; Max Koppers; Femke van Diggelen; Jessica van Heest; Jeroen A A Demmers; Janine Kirby; Pamela J Shaw; Eleonora Aronica; Wim G M Spliet; Jan H Veldink; Leonard H van den Berg; R Jeroen Pasterkamp Journal: Hum Mol Genet Date: 2013-05-15 Impact factor: 6.150
Authors: Kohji Mori; Shih-Ming Weng; Thomas Arzberger; Stephanie May; Kristin Rentzsch; Elisabeth Kremmer; Bettina Schmid; Hans A Kretzschmar; Marc Cruts; Christine Van Broeckhoven; Christian Haass; Dieter Edbauer Journal: Science Date: 2013-02-07 Impact factor: 47.728
Authors: Christopher J Donnelly; Ping-Wu Zhang; Jacqueline T Pham; Aaron R Haeusler; Aaron R Heusler; Nipun A Mistry; Svetlana Vidensky; Elizabeth L Daley; Erin M Poth; Benjamin Hoover; Daniel M Fines; Nicholas Maragakis; Pentti J Tienari; Leonard Petrucelli; Bryan J Traynor; Jiou Wang; Frank Rigo; C Frank Bennett; Seth Blackshaw; Rita Sattler; Jeffrey D Rothstein Journal: Neuron Date: 2013-10-16 Impact factor: 17.173
Authors: Alan E Renton; Elisa Majounie; Adrian Waite; Javier Simón-Sánchez; Sara Rollinson; J Raphael Gibbs; Jennifer C Schymick; Hannu Laaksovirta; John C van Swieten; Liisa Myllykangas; Hannu Kalimo; Anders Paetau; Yevgeniya Abramzon; Anne M Remes; Alice Kaganovich; Sonja W Scholz; Jamie Duckworth; Jinhui Ding; Daniel W Harmer; Dena G Hernandez; Janel O Johnson; Kin Mok; Mina Ryten; Danyah Trabzuni; Rita J Guerreiro; Richard W Orrell; James Neal; Alex Murray; Justin Pearson; Iris E Jansen; David Sondervan; Harro Seelaar; Derek Blake; Kate Young; Nicola Halliwell; Janis Bennion Callister; Greg Toulson; Anna Richardson; Alex Gerhard; Julie Snowden; David Mann; David Neary; Michael A Nalls; Terhi Peuralinna; Lilja Jansson; Veli-Matti Isoviita; Anna-Lotta Kaivorinne; Maarit Hölttä-Vuori; Elina Ikonen; Raimo Sulkava; Michael Benatar; Joanne Wuu; Adriano Chiò; Gabriella Restagno; Giuseppe Borghero; Mario Sabatelli; David Heckerman; Ekaterina Rogaeva; Lorne Zinman; Jeffrey D Rothstein; Michael Sendtner; Carsten Drepper; Evan E Eichler; Can Alkan; Ziedulla Abdullaev; Svetlana D Pack; Amalia Dutra; Evgenia Pak; John Hardy; Andrew Singleton; Nigel M Williams; Peter Heutink; Stuart Pickering-Brown; Huw R Morris; Pentti J Tienari; Bryan J Traynor Journal: Neuron Date: 2011-09-21 Impact factor: 17.173
Authors: Elisa Majounie; Alan E Renton; Kin Mok; Elise G P Dopper; Adrian Waite; Sara Rollinson; Adriano Chiò; Gabriella Restagno; Nayia Nicolaou; Javier Simon-Sanchez; John C van Swieten; Yevgeniya Abramzon; Janel O Johnson; Michael Sendtner; Roger Pamphlett; Richard W Orrell; Simon Mead; Katie C Sidle; Henry Houlden; Jonathan D Rohrer; Karen E Morrison; Hardev Pall; Kevin Talbot; Olaf Ansorge; Dena G Hernandez; Sampath Arepalli; Mario Sabatelli; Gabriele Mora; Massimo Corbo; Fabio Giannini; Andrea Calvo; Elisabet Englund; Giuseppe Borghero; Gian Luca Floris; Anne M Remes; Hannu Laaksovirta; Leo McCluskey; John Q Trojanowski; Vivianna M Van Deerlin; Gerard D Schellenberg; Michael A Nalls; Vivian E Drory; Chin-Song Lu; Tu-Hsueh Yeh; Hiroyuki Ishiura; Yuji Takahashi; Shoji Tsuji; Isabelle Le Ber; Alexis Brice; Carsten Drepper; Nigel Williams; Janine Kirby; Pamela Shaw; John Hardy; Pentti J Tienari; Peter Heutink; Huw R Morris; Stuart Pickering-Brown; Bryan J Traynor Journal: Lancet Neurol Date: 2012-03-09 Impact factor: 44.182
Authors: Timothy P Levine; Rachel D Daniels; Alberto T Gatta; Louise H Wong; Matthew J Hayes Journal: Bioinformatics Date: 2013-01-16 Impact factor: 6.937
Authors: Yingxiao Shi; Shaoyu Lin; Kim A Staats; Yichen Li; Wen-Hsuan Chang; Shu-Ting Hung; Eric Hendricks; Gabriel R Linares; Yaoming Wang; Esther Y Son; Xinmei Wen; Kassandra Kisler; Brent Wilkinson; Louise Menendez; Tohru Sugawara; Phillip Woolwine; Mickey Huang; Michael J Cowan; Brandon Ge; Nicole Koutsodendris; Kaitlin P Sandor; Jacob Komberg; Vamshidhar R Vangoor; Ketharini Senthilkumar; Valerie Hennes; Carina Seah; Amy R Nelson; Tze-Yuan Cheng; Shih-Jong J Lee; Paul R August; Jason A Chen; Nicholas Wisniewski; Victor Hanson-Smith; T Grant Belgard; Alice Zhang; Marcelo Coba; Chris Grunseich; Michael E Ward; Leonard H van den Berg; R Jeroen Pasterkamp; Davide Trotti; Berislav V Zlokovic; Justin K Ichida Journal: Nat Med Date: 2018-02-05 Impact factor: 53.440
Authors: Vincent Picher-Martel; Paul N Valdmanis; Peter V Gould; Jean-Pierre Julien; Nicolas Dupré Journal: Acta Neuropathol Commun Date: 2016-07-11 Impact factor: 7.801