Role of PATJ in stroke prognosis by modulating endothelial to mesenchymal transition through the Hippo/Notch/PI3K axis.

Through GWAS studies we identified PATJ associated with functional outcome after ischemic stroke (IS). The aim of this study was to determine PATJ role in brain endothelial cells (ECs) in the context of stroke outcome. PATJ expression analyses in patient's blood revealed that: (i) the risk allele of rs76221407 induces higher expression of PATJ, (ii) PATJ is downregulated 24 h after IS, and (iii) its expression is significantly lower in those patients with functional independence, measured at 3 months with the modified Rankin scale ((mRS) ≤2), compared to those patients with marked disability (mRS = 4-5). In mice brains, PATJ was also downregulated in the injured hemisphere at 48 h after ischemia. Oxygen-glucose deprivation and hypoxia-dependent of Hypoxia Inducible Factor-1α also caused PATJ depletion in ECs. To study the effects of PATJ downregulation, we generated PATJ-knockdown human microvascular ECs. Their transcriptomic profile evidenced a complex cell reprogramming involving Notch, TGF-ß, PI3K/Akt, and Hippo signaling that translates in morphological and functional changes compatible with endothelial to mesenchymal transition (EndMT). PATJ depletion caused loss of cell-cell adhesion, upregulation of metalloproteases, actin cytoskeleton remodeling, cytoplasmic accumulation of the signal transducer C-terminal transmembrane Mucin 1 (MUC1-C) and downregulation of Notch and Hippo signaling. The EndMT phenotype of PATJ-depleted cells was associated with the nuclear recruitment of MUC1-C, YAP/TAZ, ß-catenin, and ZEB1. Our results suggest that PATJ downregulation 24 h after IS promotes EndMT, an initial step prior to secondary activation of a pro-angiogenic program. This effect is associated with functional independence suggesting that activation of EndMT shortly after stroke onset is beneficial for stroke recovery.

DNA Methylation and Ischemic Stroke Risk: An Epigenome-Wide Association Study.

BACKGROUND: Ischemic stroke (IS) risk heritability is partly explained by genetics. Other heritable factors, such as epigenetics, could explain an unknown proportion of the IS risk. The objective of this study is to evaluate DNA methylation association with IS using epigenome-wide association studies (EWAS). METHODS: We performed a two-stage EWAS comprising 1,156 subjects. Differentially methylated positions (DMPs) and differentially methylated regions (DMRs) were assessed using the Infinium 450K and EPIC BeadChip in the discovery cohort (252 IS and 43 controls). Significant DMPs were replicated in an independent cohort (618 IS and 243 controls). Stroke subtype associations were also evaluated. Differentially methylated cell-type (DMCT) was analyzed in the replicated CpG sites using EpiDISH. We additionally performed pathway enrichment analysis and causality analysis with Mendelian randomization for the replicated CpG sites. RESULTS: A total of 957 CpG sites were epigenome-wide-significant (p ≤ 10-7) in the discovery cohort, being CpG sites in the top signals (logFC = 0.058, p = 2.35 × 10-22; logFC = 0.035, p = 3.22 × 10-22, respectively). ZFHX3 and MAP3K1 were among the most significant DMRs. In addition, 697 CpG sites were replicated considering Bonferroni-corrected p-values (p < 5.22 × 10-5). All the replicated DMPs were associated with risk of cardioembolic, atherothrombotic, and undetermined stroke. The DMCT analysis demonstrated that the significant associations were driven by natural killer cells. The pathway enrichment analysis showed overrepresentation of genes belonging to certain pathways including oxidative stress. ZFHX3 and MAP3K1 methylation was causally associated with specific stroke-subtype risk. CONCLUSION: Specific DNA methylation pattern is causally associated with IS risk. These results could be useful for specifically predicting stroke occurrence and could potentially be evaluated as therapeutic targets.

Multi-ancestry GWAS reveals excitotoxicity associated with outcome after ischaemic stroke.

During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6 h of stroke onset and NIHSS at 24 h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24 h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke.

Evolutionary diversification of the canonical Wnt signaling effector TCF/LEF in chordates.

Wnt signaling is essential during animal development and regeneration, but also plays an important role in diseases such as cancer and diabetes. The canonical Wnt signaling pathway is one of the most conserved signaling cascades in the animal kingdom, with the T-cell factor/lymphoid enhancer factor (TCF/LEF) proteins being the major mediators of Wnt/ß-catenin-regulated gene expression. In comparison with invertebrates, vertebrates possess a high diversity of TCF/LEF family genes, implicating this as a possible key change to Wnt signaling at the evolutionary origin of vertebrates. However, the precise nature of this diversification is only poorly understood. The aim of this study is to clarify orthology, paralogy, and isoform relationships within the TCF/LEF gene family within chordates via in silico comparative study of TCF/LEF gene structure, molecular phylogeny, and gene synteny. Our results support the notion that the four TCF/LEF paralog subfamilies in jawed vertebrates (gnathostomes) evolved via the two rounds of whole-genome duplications that occurred during early vertebrate evolution. Importantly, gene structure comparisons and synteny analysis of jawless vertebrate (cyclostome) TCFs suggest that a TCF7L2-like form of gene structure is a close proxy for the ancestral vertebrate structure. In conclusion, we propose a detailed evolutionary path based on a new pre-whole-genome duplication vertebrate TCF gene model. This ancestor gene model highlights the chordate and vertebrate innovations of TCF/LEF gene structure, providing the foundation for understanding the role of Wnt/ß-catenin signaling in vertebrate evolution.

Genome-wide transcriptome study in skin biopsies reveals an association of E2F4 with cadasil and cognitive impairment.

CADASIL is a small vessel disease caused by mutations in NOTCH3 that lead to an odd number of cysteines in the EGF-like repeat domain, causing protein misfolding and aggregation. The main symptoms are migraine, psychiatric disturbances, recurrent strokes and dementia, being executive function characteristically impaired. The molecular pathways altered by this receptor aggregation need to be studied further. A genome-wide transcriptome study (four cases paired with three healthy siblings) was carried out, in addition to a qRT-PCR for validation purposes (ten new cases and eight new controls). To study the expression profile by cell type of the significant mRNAs found, we performed an in situ hybridization (ISH) (nine cases and eight controls) and a research in the Single-nuclei Brain RNA-seq expression browser (SNBREB). Pathway analysis enrichment was carried out with Gene Ontology and Reactome. Neuropsychological tests were performed in five of the qRT-PCR cases. The two most significant differentially expressed mRNAs (BANP, p-value = 7.23 × 10-4 and PDCD6IP, p-value = 8.36 × 10-4) were selected for the validation study by qRT-PCR. Additionally, we selected two more mRNAs (CAMK2G, p-value = 4.52 × 10-3 and E2F4, p-value = 4.77 × 10-3) due to their association with ischemic neuronal death. E2F4 showed differential expression in the genome-wide transcriptome study and in the qRT-PCR (p = 1.23 × 10-3), and it was upregulated in CADASIL cases. Furthermore, higher E2F4 expression was associated with worse executive function (p = 2.04 × 10-2) and attention and information processing speed (IPS) (p = 8.73 × 10-2). In situ hibridization showed E2F4 expression in endothelial and vascular smooth vessel cells. In silico studies indicated that E2F4 is also expressed in brain endothelial cells. Among the most significant pathways analyzed, there was an enrichment of vascular development, cell adhesion and vesicular machinery terms and autophagy process. E2F4 is more highly expressed in the skin biopsy of CADASIL patients compared to controls, and its expression is present in endothelial cells and VSMCs. Further studies are needed to understand whether E2F4 could be useful as a biomarker, to monitor the disease or be used as a therapeutic target.

Single nucleotide variations in ZBTB46 are associated with post-thrombolytic parenchymal haematoma.

Haemorrhagic transformation is a complication of recombinant tissue-plasminogen activator treatment. The most severe form, parenchymal haematoma, can result in neurological deterioration, disability, and death. Our objective was to identify single nucleotide variations associated with a risk of parenchymal haematoma following thrombolytic therapy in patients with acute ischaemic stroke. A fixed-effect genome-wide meta-analysis was performed combining two-stage genome-wide association studies (n = 1904). The discovery stage (three cohorts) comprised 1324 ischaemic stroke individuals, 5.4% of whom had a parenchymal haematoma. Genetic variants yielding a P-value < 0.05 1 × 10-5 were analysed in the validation stage (six cohorts), formed by 580 ischaemic stroke patients with 12.1% haemorrhagic events. All participants received recombinant tissue-plasminogen activator; cases were parenchymal haematoma type 1 or 2 as defined by the European Cooperative Acute Stroke Study (ECASS) criteria. Genome-wide significant findings (P < 5 × 10-8) were characterized by in silico functional annotation, gene expression, and DNA regulatory elements. We analysed 7 989 272 single nucleotide polymorphisms and identified a genome-wide association locus on chromosome 20 in the discovery cohort; functional annotation indicated that the ZBTB46 gene was driving the association for chromosome 20. The top single nucleotide polymorphism was rs76484331 in the ZBTB46 gene [P = 2.49 × 10-8; odds ratio (OR): 11.21; 95% confidence interval (CI): 4.82-26.55]. In the replication cohort (n = 580), the rs76484331 polymorphism was associated with parenchymal haematoma (P = 0.01), and the overall association after meta-analysis increased (P = 1.61 × 10-8; OR: 5.84; 95% CI: 3.16-10.76). ZBTB46 codes the zinc finger and BTB domain-containing protein 46 that acts as a transcription factor. In silico studies indicated that ZBTB46 is expressed in brain tissue by neurons and endothelial cells. Moreover, rs76484331 interacts with the promoter sites located at 20q13. In conclusion, we identified single nucleotide variants in the ZBTB46 gene associated with a higher risk of parenchymal haematoma following recombinant tissue-plasminogen activator treatment.

Early Neurological Change After Ischemic Stroke Is Associated With 90-Day Outcome.

BACKGROUND AND PURPOSE: Large-scale observational studies of acute ischemic stroke (AIS) promise to reveal mechanisms underlying cerebral ischemia. However, meaningful quantitative phenotypes attainable in large patient populations are needed. We characterize a dynamic metric of AIS instability, defined by change in National Institutes of Health Stroke Scale score (NIHSS) from baseline to 24 hours baseline to 24 hours (NIHSSbaseline - NIHSS24hours = ΔNIHSS6-24h), to examine its relevance to AIS mechanisms and long-term outcomes. METHODS: Patients with NIHSS prospectively recorded within 6 hours after onset and then 24 hours later were enrolled in the GENISIS study (Genetics of Early Neurological Instability After Ischemic Stroke). Stepwise linear regression determined variables that independently influenced ΔNIHSS6-24h. In a subcohort of tPA (alteplase)-treated patients with large vessel occlusion, the influence of early sustained recanalization and hemorrhagic transformation on ΔNIHSS6-24h was examined. Finally, the association of ΔNIHSS6-24h with 90-day favorable outcomes (modified Rankin Scale score 0-2) was assessed. Independent analysis was performed using data from the 2 NINDS-tPA stroke trials (National Institute of Neurological Disorders and Stroke rt-PA). RESULTS: For 2555 patients with AIS, median baseline NIHSS was 9 (interquartile range, 4-16), and median ΔNIHSS6-24h was 2 (interquartile range, 0-5). In a multivariable model, baseline NIHSS, tPA-treatment, age, glucose, site, and systolic blood pressure independently predicted ΔNIHSS6-24h (R2=0.15). In the large vessel occlusion subcohort, early sustained recanalization and hemorrhagic transformation increased the explained variance (R2=0.27), but much of the variance remained unexplained. ΔNIHSS6-24h had a significant and independent association with 90-day favorable outcome. For the subjects in the 2 NINDS-tPA trials, ΔNIHSS3-24h was similarly associated with 90-day outcomes. CONCLUSIONS: The dynamic phenotype, ΔNIHSS6-24h, captures both explained and unexplained mechanisms involved in AIS and is significantly and independently associated with long-term outcomes. Thus, ΔNIHSS6-24h promises to be an easily obtainable and meaningful quantitative phenotype for large-scale genomic studies of AIS.

Multi-ancestry genetic study in 5,876 patients identifies an association between excitotoxic genes and early outcomes after acute ischemic stroke.

During the first hours after stroke onset neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between NIH stroke scale (NIHSS) within six hours of stroke onset and NIHSS at 24h (ΔNIHSS). A total of 5,876 individuals from seven countries (Spain, Finland, Poland, United States, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of ΔNIHSS variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture than that of stroke risk. Seven loci (2p25.1, 2q31.2, 2q33.3, 4q34.3, 5q33.2, 6q26 and 7p21.1) were genome-wide significant and explained 2.1% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each loci. eQTL mapping and SMR indicate that ADAM23 (log Bayes Factor (LBF)=6.34) was driving the association for 2q33.3. Gene based analyses suggested that GRIA1 (LBF=5.26), which is predominantly expressed in brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated PARK2 (LBF=5.30) and ABCB5 (LBF=5.70) for the 6q26 and 7p21.1 loci. Human brain single nuclei RNA-seq indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23 , a pre-synaptic protein, and GRIA1 , a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provides the first evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischemic stroke. RESEARCH INTO CONTEXT: Evidence before this study: No previous genome-wide association studies have investigated the genetic architecture of early outcomes after ischemic stroke.Added Value of this study: This is the first study that investigated genetic influences on early outcomes after ischemic stroke using a genome-wide approach, revealing seven genome-wide significant loci. A unique aspect of this genetic study is the inclusion of all of the major ethnicities by recruiting from participants throughout the world. Most genetic studies to date have been limited to populations of European ancestry.Implications of all available evidence: The findings provide the first evidence that genes implicating excitotoxicity contribute to human acute ischemic stroke, and demonstrates proof of principle that GWAS of acute ischemic stroke patients can reveal mechanisms involved in ischemic brain injury.

DNA methylation of MMPs and TIMPs in atherothrombosis process in carotid plaques and blood tissues.

BACKGROUND AND PURPOSE: Polymorphisms and serum levels of Matrix Metalloproteinases (MMP) and Tissue Inhibitor of Metalloproteinases (TIMP) have been studied with regard to atheromatous plaques and ischemic stroke, while no studies of DNA methylation (DNAm) patterns of MMP or TIMP have been performed to that end. Here, we evaluate DNAm levels of the MMP and TIMP gene families in human carotid plaques and blood samples of atherothrombotic stroke patients. METHODS: We profiled the DNAm status of stable and ulcerated atherosclerotic plaques obtained as pair sets from three patients who underwent carotid endarterectomy surgery. We selected 415 CpG sites, mapping into MMPs and TIMPs genes for further study. Secondly, the statistically associated CpG sites were analyzed in blood samples from two separate atherothrombotic stroke cohorts (total sample size = 307), ischemic stroke-cohort 1 (ISC-1): 37 atherothrombotic patients and 6 controls, ischemic stroke-cohort 2 (ISC-2): 80 atherothrombotic patients and 184 controls. DNAm levels from plaque tissue and blood samples were evaluated using a high-density microarray Infinium, HumanMethylation450 BeadChip and Infinium MethylationEPIC BeadChip. RESULTS: Three CpG sites were statistically significantly associated with unstable plaque portions; cg02969624, q-value = 0.035 (TIMP2), and cg04316754, q-value = 0.037 (MMP24) were hypermethylated, while cg24211657 q-value = 0.035 (TIMP2) was hypomethylated. Association of cg04316754 (MMP24) methylation levels with atherothrombotic risk was also observed in blood tissue: ISC-1 p-values = 0.03, ISC-2 p-value = 1.9 × 10-04. CONCLUSIONS: The results suggest different DNAm status of MMP24 between stable and unstable atherothrombotic carotid plaques, and between atherothrombotic stroke and controls in blood samples.

Genome-Wide Association Study of VKORC1 and CYP2C9 on acenocoumarol dose, stroke recurrence and intracranial haemorrhage in Spain.

Acenocoumarol is an oral anticoagulant with significant interindividual dose variations. Variants in CYP2C9 and VKORC1 have been associated with acenocoumarol maintenance dose. We analysed whether any of the 49 polymorphisms in CYP2C9 and VKORC1 previously associated with acenocoumarol maintenance dose in a Genome-Wide Association study (GWAs) in Dutch population are associated with stroke recurrence, intracranial haemorrhage (ICH) and acenocoumarol maintenance dose in a Spanish population. We performed a GWAs using Human Core Exome-chip (Illumina) in 78 patients stroke patients treated with acenocoumarol for secondary prevention enrolled as part of the prospective investigator-initiated study (IIS) SEDMAN Study. Patients were followed-up a median of 12.8 months. Three and eight patients had recurrent stroke and ICH events, respectively. We found 14 of the 49 published variants associated with acenocoumarol maintenance dose (p < 0.05). Six polymorphisms were associated with stroke recurrence and four variants with ICH (p < 0.05). In conclusion, variants in VKORC1 and CYP2C9 are associated with acenocoumarol maintenance dose, stroke recurrence and ICH in a Spanish cohort. These results highlight the relevance of studying pharmacogenetics associated with efficacy and safety of anticoagulant drugs and justify studies with larger sample size and different ethnic populations.

Clinical Variables and Genetic Risk Factors Associated with the Acute Outcome of Ischemic Stroke: A Systematic Review.

Stroke is a complex disease and one of the main causes of morbidity and mortality among the adult population. A huge variety of factors is known to influence patient outcome, including demographic variables, comorbidities or genetics. In this review, we expound what is known about the influence of clinical variables and related genetic risk factors on ischemic stroke outcome, focusing on acute and subacute outcome (within 24 to 48 hours after stroke and until day 10, respectively), as they are the first indicators of stroke damage. We searched the PubMed data base for articles that investigated the interaction between clinical variables or genetic factors and acute or subacute stroke outcome. A total of 61 studies were finally included in this review. Regarding the data collected, the variables consistently associated with acute stroke outcome are: glucose levels, blood pressure, presence of atrial fibrillation, prior statin treatment, stroke severity, type of acute treatment performed, severe neurological complications, leukocyte levels, and genetic risk factors. Further research and international efforts are required in this field, which should include genome-wide association studies.

Genome-Wide Association Study of White Blood Cell Counts in Patients With Ischemic Stroke.

Background and Purpose- Immune cells play a key role in the first 24h poststroke (acute phase), being associated with stroke outcome. We aimed to find genetic risk factors associated with leukocyte counts during the acute phase of stroke. Methods- Ischemic stroke patients with leukocyte counts data during the first 24h were included. Genome-wide association study and gene expression studies were performed. Results- Our genome-wide association study, which included 2064 (Discovery) and 407 (Replication) patients, revealed a new locus (14q24.3) associated with leukocyte counts. After Joint analysis (n=2471) 5 more polymorphisms reached genome-wide significance (P<5×10-8). The 14q24.3 locus was associated with acute stroke outcome (rs112809786, P=0.036) and with ACOT1 and PTGR2 gene expression. Previous polymorphisms associated with leukocyte counts in general-population did not show any significance in our study. Conclusions- We have found the first locus associated with leukocyte counts in ischemic stroke, also associated with acute outcome. Genetic analysis of acute endophenotypes could be useful to find the genetic factors associated with stroke outcome. Our findings suggested a different modulation of immune cells in stroke compared with healthy conditions.

Validation of a clinical-genetics score to predict hemorrhagic transformations after rtPA.

OBJECTIVE: To validate the Genot-PA score, a clinical-genetic logistic regression score that stratifies the thrombolytic therapy safety, in a new cohort of patients with stroke. METHODS: We enrolled 1,482 recombinant tissue plasminogen activator (rtPA)-treated patients with stroke in Spain and Finland from 2003 to 2016. Cohorts were analyzed on the basis of ethnicity and therapy: Spanish patients treated with IV rtPA within 4.5 hours of onset (cohort A and B) or rtPA in combination with mechanical thrombectomy within 6 hours of onset (cohort C) and Finnish participants treated with IV rtPA within 4.5 hours of onset (cohort D). The Genot-PA score was calculated, and hemorrhagic transformation (HT) and parenchymal hematoma (PH) risks were determined for each score stratum. RESULTS: Genot-PA score was tested in 1,324 (cohort A, n = 726; B, n = 334; C, n = 54; and D, n = 210) patients who had enough information to complete the score. Of these, 213 (16.1%) participants developed HT and 85 (6.4%) developed PH. In cohorts A, B, and D, HT occurrence was predicted by the score (p = 2.02 × 10-6, p = 0.023, p = 0.033); PH prediction was associated in cohorts A through C (p = 0.012, p = 0.034, p = 5.32 × 10-4). Increased frequency of PH events from the lowest to the highest risk group was found (cohort A 4%-15.7%, cohort B 1.5%-18.2%, cohort C 0%-100%). The best odds ratio for PH prediction in the highest-risk group was obtained in cohort A (odds ratio 5.16, 95% confidence interval 1.46-18.08, p = 0.009). CONCLUSION: The Genot-PA score predicts HT in patients with stroke treated with IV rtPA. Moreover, in an exploratory study, the score was associated with PH risk in mechanical thrombectomy-treated patients.

PATJ Low Frequency Variants Are Associated With Worse Ischemic Stroke Functional Outcome.

RATIONALE: Ischemic stroke is among the leading causes of adult disability. Part of the variability in functional outcome after stroke has been attributed to genetic factors but no locus has been consistently associated with stroke outcome. OBJECTIVE: Our aim was to identify genetic loci influencing the recovery process using accurate phenotyping to produce the largest GWAS (genome-wide association study) in ischemic stroke recovery to date. METHODS AND RESULTS: A 12-cohort, 2-phase (discovery-replication and joint) meta-analysis of GWAS included anterior-territory and previously independent ischemic stroke cases. Functional outcome was recorded using 3-month modified Rankin Scale. Analyses were adjusted for confounders such as discharge National Institutes of Health Stroke Scale. A gene-based burden test was performed. The discovery phase (n=1225) was followed by open (n=2482) and stringent joint-analyses (n=1791). Those cohorts with modified Rankin Scale recorded at time points other than 3-month or incomplete data on previous functional status were excluded in the stringent analyses. Novel variants in PATJ (Pals1-associated tight junction) gene were associated with worse functional outcome at 3-month after stroke. The top variant was rs76221407 (G allele, ß=0.40, P=1.70×10-9). CONCLUSIONS: Our results identify a set of common variants in PATJ gene associated with 3-month functional outcome at genome-wide significance level. Future studies should examine the role of PATJ in stroke recovery and consider stringent phenotyping to enrich the information captured to unveil additional stroke outcome loci.

Diatom bloom-derived biotoxins cause aberrant development and gene expression in the appendicularian chordate Oikopleura dioica.

Investigating environmental hazards than could affect appendicularians is of prime ecological interest because they are among the most abundant components of the mesozooplankton. This work shows that embryo development of the appendicularian Oikopleura dioica is compromised by diatom bloom-derived biotoxins, even at concentrations in the same range as those measured after blooms. Developmental gene expression analysis of biotoxin-treated embryos uncovers an aberrant golf ball-like phenotype affecting morphogenesis, midline convergence, and tail elongation. Biotoxins induce a rapid upregulation of defensome genes, and considerable delay and silencing of zygotic transcription of developmental genes. Upon a possible future intensification of blooms associated with ocean warming and acidification, our work puts an alert on the potential impact that an increase of biotoxins may have on marine food webs, and points to defensome genes as molecular biosensors that marine ecologists could use to monitor the genetic stress of natural populations exposed to microalgal blooms.

Coelimination and Survival in Gene Network Evolution: Dismantling the RA-Signaling in a Chordate.

The bloom of genomics is revealing gene loss as a pervasive evolutionary force generating genetic diversity that shapes the evolution of species. Outside bacteria and yeast, however, the understanding of the process of gene loss remains elusive, especially in the evolution of animal species. Here, using the dismantling of the retinoic acid metabolic gene network (RA-MGN) in the chordate Oikopleura dioica as a case study, we combine approaches of comparative genomics, phylogenetics, biochemistry, and developmental biology to investigate the mutational robustness associated to biased patterns of gene loss. We demonstrate the absence of alternative pathways for RA-synthesis in O. dioica, which suggests that gene losses of RA-MGN were not compensated by mutational robustness, but occurred in a scenario of regressive evolution. In addition, the lack of drastic phenotypic changes associated to the loss of RA-signaling provides an example of the inverse paradox of Evo-Devo. This work illustrates how the identification of patterns of gene coelimination-in our case five losses (Rdh10, Rdh16, Bco1, Aldh1a, and Cyp26)-is a useful strategy to recognize gene network modules associated to distinct functions. Our work also illustrates how the identification of survival genes helps to recognize neofunctionalization events and ancestral functions. Thus, the survival and extensive duplication of Cco and RdhE2 in O. dioica correlated with the acquisition of complex compartmentalization of expression domains in the digestive system and a process of enzymatic neofunctionalization of the Cco, while the surviving Aldh8 could be related to its ancestral housekeeping role against toxic aldehydes.

Oikopleura dioica culturing made easy: a low-cost facility for an emerging animal model in EvoDevo.

The genome sequencing and the development of RNAi knockdown technologies in the urochordate Oikopleura dioica are making this organism an attractive emergent model in the field of EvoDevo. To succeed as a new animal model, however, an organism needs to be easily and affordably cultured in the laboratory. Nowadays, there are only two facilities in the world capable to indefinitely maintain Oikopleura dioica, one in the SARS institute (Bergen, Norway) and the other in the Osaka University (Japan). Here, we describe the setup of a new facility in the University of Barcelona (Spain) in which we have modified previously published husbandry protocols to optimize the weekly production of thousands of embryos and hundreds of mature animals using the minimum amount of space, human resources, and technical equipment. This optimization includes novel protocols of cryopreservation and solid cultures for long-term maintenance of microalgal stocks-Chaetoceros calcitrans, Isochrysis sp., Rhinomonas reticulate, and Synechococcus sp.-needed for Oikopleura dioica feeding. Our culture system maintains partially inbred lines healthy with similar characteristics to wild animals, and it is easily expandable to satisfy on demand the needs of any laboratory that may wish to use Oikopleura dioica as a model organism.