MAP3K6 Mutations in a Neurovascular Disease Causing Stroke, Cognitive Impairment, and Tremor.

OBJECTIVE: To describe a possible novel genetic mechanism for cerebral small vessel disease (cSVD) and stroke. METHODS: We studied a Swedish kindred with ischemic stroke and intracerebral hemorrhage, tremor, dysautonomia, and mild cognitive decline. Members were examined clinically, radiologically, and by histopathology. Genetic workup included whole-exome sequencing (WES) and whole-genome sequencing (WGS) and intrafamilial cosegregation analyses. RESULTS: Fifteen family members were examined clinically. Twelve affected individuals had white matter hyperintensities and 1 or more of (1) stroke episodes, (2) clinically silent lacunar ischemic lesions, and (3) cognitive dysfunction. All affected individuals had tremor and/or atactic gait disturbance. Mild symmetric basal ganglia calcifications were seen in 3 affected members. Postmortem examination of 1 affected member showed pathologic alterations in both small and large arteries the brain. Skin biopsies of 3 affected members showed extracellular amorphous deposits within the subepidermal zone, which may represent degenerated arterioles. WES or WGS did not reveal any potentially disease-causing variants in known genes for cSVDs or idiopathic basal ganglia calcification, but identified 1 heterozygous variant, NM_004672.4 MAP3K6 c.322G>A p.(Asp108Asn), that cosegregated with the disease in this large family. MAP3K6 has known functions in angiogenesis and affects vascular endothelial growth factor expression, which may be implicated in cerebrovascular disease. CONCLUSIONS: Our data strongly suggest the MAP3K6 variant to be causative for this novel disease phenotype, but the absence of functional data and the present lack of additional families with this disease and MAP3K6 mutations still limit the formal evidence for the variant's pathogenicity.

Identification of candidate genetic variants and altered protein expression in neural stem and mature neural cells support altered microtubule function to be an essential component in bipolar disorder.

Identification of causative genetic variants leading to the development of bipolar disorder (BD) could result in genetic tests that would facilitate diagnosis. A better understanding of affected genes and pathways is also necessary for targeting of genes that may improve treatment strategies. To date several susceptibility genes have been reported from genome-wide association studies (GWAS), but little is known about specific variants that affect disease development. Here, we performed quantitative proteomics and whole-genome sequencing (WGS). Quantitative proteomics revealed NLRP2 as the most significantly up-regulated protein in neural stem cells and mature neural cells obtained from BD-patient cell samples. These results are in concordance with our previously published transcriptome analysis. Furthermore, the levels of FEZ2 and CADM2 proteins were also significantly differentially expressed in BD compared to control derived cells. The levels of FEZ2 were significantly downregulated in neural stem cells (NSC) while CADM2 was significantly up-regulated in mature neuronal cell culture. Promising novel candidate mutations were identified in the ANK3, NEK3, NEK7, TUBB, ANKRD1, and BRD2 genes. A literature search of candidate variants and deregulated proteins revealed that there are several connections to microtubule function for the molecules putatively involved. Microtubule function in neurons is critical for axon structure and axonal transport. A functional dynamic microtubule is also needed for an advocate response to cellular and environmental stress. If microtubule dynamics is compromised by mutations, it could be followed by deregulated expression forming a possible explanation for the inherited vulnerability to stressful life events that have been proposed to trigger mood episodes in BD patients.

Whole-Exome Sequencing in 22 Young Ischemic Stroke Patients With Familial Clustering of Stroke.

Backgrounds and Purpose- Although new methods for genetic analyses are rapidly evolving, there are currently knowledge gaps in how to detect Mendelian forms of stroke. Methods- We performed whole-exome sequencing in 22 probands, under 56 years at their first ischemic stroke episode, from multi-incident stroke families. With the use of a comprehensive stroke-gene panel, we searched for variants in stroke-related genes. The probands' clinical stroke subtype was related to clinical characteristics previously associated with pathogenic variants in these genes. Relatives were genotyped in 7 families to evaluate stroke-gene variants of unknown significance. In 2 larger families with embolic stroke of unknown source, whole-exome sequencing was performed in additional members to examine the possibility of identifying new stroke genes. Results- Six of 22 probands carried pathogenic or possibly pathogenic variants in genes reported to be associated with their stroke subtype. A known pathogenic variant in NOTCH3 and a possibly pathogenic variant in ACAD9 gene were identified. A novel JAK2:c.3188G>A (p.Arg1063His) mutation was seen in a proband with embolic stroke of undetermined source and prothrombotic status. However, penetrance in the family was incomplete. COL4A2:c.3368A>G (p.Glu1123Gly) was detected in 2 probands but did not cosegregate with the disease in their families. Whole-exome sequencing in multiple members of 2 pedigrees with embolic stroke of undetermined source revealed possibly pathogenic variants in genes not previously associated with stroke, GPR142:c.148C>G (p.Leu50Val), and PTPRN2:c.2416A>G (p.Ile806Val); LRRC1 c.808A>G (p.Ile270Val), SLC7A10c.1294dupG (p.Val432fs), IKBKB: c.1070C>T (p.Ala357Val), and OXGR1 c.392G>A (p.Arg131His), respectively. Conclusions- Screening with whole-exome sequencing using a comprehensive stroke-gene panel may identify rare monogenic forms of stroke, but careful evaluation of clinical characteristics and potential pathogenicity of novel variants remain important. In our study, the majority of individuals with familial aggregation of stroke lacked any identified genetic causes.

A spontaneous mitonuclear epistasis converging on Rieske Fe-S protein exacerbates complex III deficiency in mice.

We previously observed an unexpected fivefold (35 vs. 200 days) difference in the survival of respiratory chain complex III (CIII) deficient Bcs1lp.S78G mice between two congenic backgrounds. Here, we identify a spontaneous homoplasmic mtDNA variant (m.G14904A, mt-Cybp.D254N), affecting the CIII subunit cytochrome b (MT-CYB), in the background with short survival. We utilize maternal inheritance of mtDNA to confirm this as the causative variant and show that it further decreases the low CIII activity in Bcs1lp.S78G tissues to below survival threshold by 35 days of age. Molecular dynamics simulations predict D254N to restrict the flexibility of MT-CYB ef loop, potentially affecting RISP dynamics. In Rhodobacter cytochrome bc1 complex the equivalent substitution causes a kinetics defect with longer occupancy of RISP head domain towards the quinol oxidation site. These findings represent a unique case of spontaneous mitonuclear epistasis and highlight the role of mtDNA variation as modifier of mitochondrial disease phenotypes.

Choroid plexus transcriptome and ultrastructure analysis reveals a TLR2-specific chemotaxis signature and cytoskeleton remodeling in leukocyte trafficking.

Perinatal infection and inflammation are major risk factors for injury in the developing brain, however, underlying mechanisms are not fully understood. Leukocyte migration to the cerebrospinal fluid (CSF) and brain is a hallmark of many pathologies of the central nervous system including those in neonates. We previously reported that systemic activation of Toll-like receptor (TLR) 2, a major receptor for gram-positive bacteria, by agonist Pam3CSK4 (P3C) resulted in dramatic neutrophil and monocyte infiltration to the CSF and periventricular brain of neonatal mice, an effect that was absent by the TLR4 agonist, LPS. Here we first report that choroid plexus is a route of TLR2-mediated leukocyte infiltration to the CSF by performing flow cytometry and transmission electron microscopy (TEM) of the choroid plexus. Next, we exploited the striking discrepancy between P3C and LPS effects on cell migration to determine the pathways regulating leukocyte trafficking through the choroid plexus. We performed RNA sequencing on the choroid plexus after administration of P3C and LPS to postnatal day 8 mice. A cluster gene analysis revealed a TLR2-specific signature of chemotaxis represented by 80-fold increased expression of the gene Ccl3 and 1000-fold increased expression of the gene Cxcl2. Ingenuity pathway analysis (IPA) revealed TLR2-specific molecular signaling related to cytoskeleton organization (e.g. actin signaling) as well as inositol phospholipids biosynthesis and degradation. This included upregulation of genes such as Rac2 and Micall2. In support of IPA results, ultrastructural analysis by TEM revealed clefting and perforations in the basement membrane of the choroid plexus epithelial cells in P3C-treated mice. In summary, we show that the choroid plexus is a route of TLR2-mediated transmigration of neutrophils and monocytes to the developing brain, and reveal previously unrecognized mechanisms that includes a specific chemotaxis profile as well as pathways regulating cytoskeleton and basement membrane remodeling.

Low Frequency ALK Hotspots Mutations In Neuroblastoma Tumours Detected By Ultra-deep Sequencing: Implications For ALK Inhibitor Treatment.

The ALK tyrosine kinase receptor is oncogenically activated in neuroblastoma. Whereas numerous ALK fusion genes have been reported in different malignancies, in neuroblastoma ALK is mainly activated through point mutations. Three hotspot residues (F1174, F1245, and R1275) account for 85% of mutant ALK seen in neuroblastoma. In a cohort of 105 Swedish neuroblastoma cases of all stages, these hotspot regions were re-sequenced (>5000X). ALK mutations were detected in 16 of 105 patients (range of variant allele fraction: 2.7-60%). Mutations at the F1174 and F1245 hotspot were observed in eleven and three cases respectively. ALK mutations were also detected at the I1171 and L1240 codons in one tumor each. No mutations were detected at R1275. Sanger sequencing could confirm ALK status for all mutated samples with variant allele fraction above 15%. Four of the samples with subclonal ALK mutation fraction below this would have gone undetected relying on Sanger sequencing only. No distinct mutation spectrum in relation to neuroblastoma tumours genomic subtypes could be detected although there was a paucity of ALK mutations among 11q-deleted tumors. As ALK mutations status opens up an excellent opportunity for application of small molecule inhibitors targeting ALK, early and sensitive detection of ALK alterations is clinically important considering its potential role in tumour progression.

Genome-wide multi-omics profiling of the 8p11-p12 amplicon in breast carcinoma.

Genomic instability contributes to the neoplastic phenotype by deregulating key cancer-related genes, which in turn can have a detrimental effect on patient outcome. DNA amplification of the 8p11-p12 genomic region has clinical and biological implications in multiple malignancies, including breast carcinoma where the amplicon has been associated with tumor progression and poor prognosis. However, oncogenes driving increased cancer-related death and recurrent genetic features associated with the 8p11-p12 amplicon remain to be identified. In this study, DNA copy number and transcriptome profiling data for 229 primary invasive breast carcinomas (corresponding to 185 patients) were evaluated in conjunction with clinicopathological features to identify putative oncogenes in 8p11-p12 amplified samples. Illumina paired-end whole transcriptome sequencing and whole-genome SNP genotyping were subsequently performed on 23 samples showing high-level regional 8p11-p12 amplification to characterize recurrent genetic variants (SNPs and indels), expressed gene fusions, gene expression profiles and allelic imbalances. We now show previously undescribed chromothripsis-like patterns spanning the 8p11-p12 genomic region and allele-specific DNA amplification events. In addition, recurrent amplification-specific genetic features were identified, including genetic variants in the HIST1H1E and UQCRHL genes and fusion transcripts containing MALAT1 non-coding RNA, which is known to be a prognostic indicator for breast cancer and stimulated by estrogen. In summary, these findings highlight novel candidate targets for improved treatment of 8p11-p12 amplified breast carcinomas.

Serum neutrophil gelatinase-associated lipocalin (NGAL) concentration is independently associated with mortality in patients with acute coronary syndrome.

BACKGROUND: Circulating neutrophil gelatinase-associated lipocalin (NGAL) concentration increases in cardiovascular disease, but the long-term prognostic value of NGAL concentration has not been evaluated in acute coronary syndrome (ACS). We examined the association between NGAL concentration and prognosis in patients with ACS after non-ST-elevation myocardial infarction (NSTEMI) or STEMI. METHODS AND RESULTS: NGAL concentration was measured in blood from 1121 consecutive ACS patients (30% women, mean age 65 years) on the first morning after admission. After adjustment for 14 variables, NGAL concentration predicted long-term (median 167 months) mortality (hazard ratio [HR] 1.33, 95% confidence interval [CI] 1.10-1.61, P = 0.003) for quartile (q) 4 of NGAL concentration. NGAL concentrations also predicted long-term mortality (HR = 1.63, 95% CI 1.31-2.03, P < 0.001, N = 741) when adjusting for Global Registry of Acute Coronary Events (GRACE) score, left ventricular ejection fraction (LVEF), and pro-B-type natriuretic peptide (proBNP) and C-reactive protein (CRP) concentrations. With these adjustments, NGAL concentration predicted long-term mortality in NSTEMI patients (HR = 2.02, 95% CI 1.50-2.72, P < 0.001) but not in STEMI patients (HR = 1.32, 95% CI 0.95-1.83, P = 0.100). In all patients, the combination of NGAL concentration and GRACE score yielded an HR of 5.56 (95% CI 4.37-7.06, P < 0.001) for q4/q4 for both variables. CONCLUSION: NGAL concentration in ACS is associated with long-term prognosis after adjustment for clinical confounders. Measuring circulating NGAL concentration may help to identify patients-particularly those with NSTEMI-needing closer follow-up after ACS.

Correction: A new GTF2I-BRAF fusion mediating MAPK pathway activation in pilocytic astrocytoma.

[This corrects the article DOI: 10.1371/journal.pone.0175638.].

A new GTF2I-BRAF fusion mediating MAPK pathway activation in pilocytic astrocytoma.

Pilocytic astrocytoma (PA) is the most common pediatric brain tumor. A recurrent feature of PA is deregulation of the mitogen activated protein kinase (MAPK) pathway most often through KIAA1549-BRAF fusion, but also by other BRAF- or RAF1-gene fusions and point mutations (e.g. BRAFV600E). These features may serve as diagnostic and prognostic markers, and also facilitate development of targeted therapy. The aims of this study were to characterize the genetic alterations underlying the development of PA in six tumor cases, and evaluate methods for fusion oncogene detection. Using a combined analysis of RNA sequencing and copy number variation data we identified a new BRAF fusion involving the 5' gene fusion partner GTF2I (7q11.23), not previously described in PA. The new GTF2I-BRAF 19-10 fusion was found in one case, while the other five cases harbored the frequent KIAA1549-BRAF 16-9 fusion gene. Similar to other BRAF fusions, the GTF2I-BRAF fusion retains an intact BRAF kinase domain while the inhibitory N-terminal domain is lost. Functional studies on GTF2I-BRAF showed elevated MAPK pathway activation compared to BRAFWT. Comparing fusion detection methods, we found Fluorescence in situ hybridization with BRAF break apart probe as the most sensitive method for detection of different BRAF rearrangements (GTF2I-BRAF and KIAA1549-BRAF). Our finding of a new BRAF fusion in PA further emphasis the important role of B-Raf in tumorigenesis of these tumor types. Moreover, the consistency and growing list of BRAF/RAF gene fusions suggests these rearrangements to be informative tumor markers in molecular diagnostics, which could guide future treatment strategies.

Utilizing the Dog Genome in the Search for Novel Candidate Genes Involved in Glioma Development-Genome Wide Association Mapping followed by Targeted Massive Parallel Sequencing Identifies a Strongly Associated Locus.

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10-8). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility.

Melanoma patient-derived xenografts accurately model the disease and develop fast enough to guide treatment decisions.

The development of novel therapies against melanoma would benefit from individualized tumor models to ensure the rapid and accurate identification of biomarkers of therapy response. Previous studies have suggested that patient-derived xenografts (PDXes) could be useful. However, the utility of PDXes in guiding real-time treatment decisions has only been reported in anecdotal forms. Here tumor biopsies from patients with stage III and IV metastatic malignant melanoma were transplanted into immunocompromised mice to generate PDXes. 23/26 melanoma biopsies generated serially transplantable PDX models, and their histology, mutation status and expression profile resembled their corresponding patient biopsy. The potential treatment for one patient was revealed by an in vitro drug screen and treating PDXes with the MEK inhibitor trametinib. In another patient, the BRAF mutation predicted the response of both the patient and its corresponding PDXes to MAPK-targeted therapy. Importantly, in this unselected group of patients, the time from biopsy for generation of PDXes until death was significantly longer than the time required to reach the treatment phase of the PDXes. Thus, it could be clinically meaningful to use this type of platform for melanoma patients as a pre-selection tool in clinical trials.

A frameshift mutation in golden retriever dogs with progressive retinal atrophy endorses SLC4A3 as a candidate gene for human retinal degenerations.

Progressive retinal atrophy (PRA) in dogs, the canine equivalent of retinitis pigmentosa (RP) in humans, is characterised by vision loss due to degeneration of the photoreceptor cells in the retina, eventually leading to complete blindness. It affects more than 100 dog breeds, and is caused by numerous mutations. RP affects 1 in 4000 people in the Western world and 70% of causal mutations remain unknown. Canine diseases are natural models for the study of human diseases and are becoming increasingly useful for the development of therapies in humans. One variant, prcd-PRA, only accounts for a small proportion of PRA cases in the Golden Retriever (GR) breed. Using genome-wide association with 27 cases and 19 controls we identified a novel PRA locus on CFA37 (p(raw) = 1.94×10(-10), p(genome) = 1.0×10(-5)), where a 644 kb region was homozygous within cases. A frameshift mutation was identified in a solute carrier anion exchanger gene (SLC4A3) located within this region. This variant was present in 56% of PRA cases and 87% of obligate carriers, and displayed a recessive mode of inheritance with full penetrance within those lineages in which it segregated. Allele frequencies are approximately 4% in the UK, 6% in Sweden and 2% in France, but the variant has not been found in GRs from the US. A large proportion of cases (approximately 44%) remain unexplained, indicating that PRA in this breed is genetically heterogeneous and caused by at least three mutations. SLC4A3 is important for retinal function and has not previously been associated with spontaneously occurring retinal degenerations in any other species, including humans.

A novel unstable duplication upstream of HAS2 predisposes to a breed-defining skin phenotype and a periodic fever syndrome in Chinese Shar-Pei dogs.

Hereditary periodic fever syndromes are characterized by recurrent episodes of fever and inflammation with no known pathogenic or autoimmune cause. In humans, several genes have been implicated in this group of diseases, but the majority of cases remain unexplained. A similar periodic fever syndrome is relatively frequent in the Chinese Shar-Pei breed of dogs. In the western world, Shar-Pei have been strongly selected for a distinctive thick and heavily folded skin. In this study, a mutation affecting both these traits was identified. Using genome-wide SNP analysis of Shar-Pei and other breeds, the strongest signal of a breed-specific selective sweep was located on chromosome 13. The same region also harbored the strongest genome-wide association (GWA) signal for susceptibility to the periodic fever syndrome (p(raw) = 2.3 × 10⁻6, p(genome) = 0.01). Dense targeted resequencing revealed two partially overlapping duplications, 14.3 Kb and 16.1 Kb in size, unique to Shar-Pei and upstream of the Hyaluronic Acid Synthase 2 (HAS2) gene. HAS2 encodes the rate-limiting enzyme synthesizing hyaluronan (HA), a major component of the skin. HA is up-regulated and accumulates in the thickened skin of Shar-Pei. A high copy number of the 16.1 Kb duplication was associated with an increased expression of HAS2 as well as the periodic fever syndrome (p < 0.0001). When fragmented, HA can act as a trigger of the innate immune system and stimulate sterile fever and inflammation. The strong selection for the skin phenotype therefore appears to enrich for a pleiotropic mutation predisposing these dogs to a periodic fever syndrome. The identification of HA as a major risk factor for this canine disease raises the potential of this glycosaminoglycan as a risk factor for human periodic fevers and as an important driver of chronic inflammation.

Localization of canine brachycephaly using an across breed mapping approach.

The domestic dog, Canis familiaris, exhibits profound phenotypic diversity and is an ideal model organism for the genetic dissection of simple and complex traits. However, some of the most interesting phenotypes are fixed in particular breeds and are therefore less tractable to genetic analysis using classical segregation-based mapping approaches. We implemented an across breed mapping approach using a moderately dense SNP array, a low number of animals and breeds carefully selected for the phenotypes of interest to identify genetic variants responsible for breed-defining characteristics. Using a modest number of affected (10-30) and control (20-60) samples from multiple breeds, the correct chromosomal assignment was identified in a proof of concept experiment using three previously defined loci; hyperuricosuria, white spotting and chondrodysplasia. Genome-wide association was performed in a similar manner for one of the most striking morphological traits in dogs: brachycephalic head type. Although candidate gene approaches based on comparable phenotypes in mice and humans have been utilized for this trait, the causative gene has remained elusive using this method. Samples from nine affected breeds and thirteen control breeds identified strong genome-wide associations for brachycephalic head type on Cfa 1. Two independent datasets identified the same genomic region. Levels of relative heterozygosity in the associated region indicate that it has been subjected to a selective sweep, consistent with it being a breed defining morphological characteristic. Genotyping additional dogs in the region confirmed the association. To date, the genetic structure of dog breeds has primarily been exploited for genome wide association for segregating traits. These results demonstrate that non-segregating traits under strong selection are equally tractable to genetic analysis using small sample numbers.

Genome-wide association mapping identifies multiple loci for a canine SLE-related disease complex.

The unique canine breed structure makes dogs an excellent model for studying genetic diseases. Within a dog breed, linkage disequilibrium is extensive, enabling genome-wide association (GWA) with only around 15,000 SNPs and fewer individuals than in human studies. Incidences of specific diseases are elevated in different breeds, indicating that a few genetic risk factors might have accumulated through drift or selective breeding. In this study, a GWA study with 81 affected dogs (cases) and 57 controls from the Nova Scotia duck tolling retriever breed identified five loci associated with a canine systemic lupus erythematosus (SLE)-related disease complex that includes both antinuclear antibody (ANA)-positive immune-mediated rheumatic disease (IMRD) and steroid-responsive meningitis-arteritis (SRMA). Fine mapping with twice as many dogs validated these loci. Our results indicate that the homogeneity of strong genetic risk factors within dog breeds allows multigenic disorders to be mapped with fewer than 100 cases and 100 controls, making dogs an excellent model in which to identify pathways involved in human complex diseases.

Glomerular transcriptome changes associated with lipopolysaccharide-induced proteinuria.

BACKGROUND: Global gene expression patterns have recently been characterized in normal glomeruli, but gene expression changes that accompany glomerular disease remain poorly characterized. METHOD: Here, we mapped global glomerular gene expression profile changes occurring in conjunction with lipopolysaccharide (LPS)-induced proteinuria in mice. RESULTS: We observed dramatic transcriptional reprogramming in glomeruli in response to LPS, representing some 20% of all genes and about 45% of the genes that are normally highly expressed in glomeruli. Bioinformatic analysis revealed significant changes in transcripts encoding proteins involved in the regulation of adherence junctions, actin cytoskeleton and survival in podocytes. In the LPS-treated mice, we observed dysregulation of genes expressed in glomerular endothelial and mesangial cells and in podocytes, there was also a significant decrease in podocyte number. Moreover, collagen alpha1, alpha2 (IV) and laminin 10 (laminin alpha 5 beta 1 gamma 1), which are expressed in immature glomeruli, were upregulated in the glomeruli of LPS-treated mice, suggesting remodeling of the glomerular basement membrane and activation of mesangial cells. By superimposing the LPS-induced changes onto GlomNet, a protein-protein interaction network was predicted for podocyte proteins affected by LPS. CONCLUSIONS: The detected changes in glomerular gene expression and their involvement in protein interaction networks provide putative markers for early and transient glomerular injury and proteinuria.