Adolescent exposure to Δ9-tetrahydrocannabinol alters the transcriptional trajectory and dendritic architecture of prefrontal pyramidal neurons.

Neuronal circuits within the prefrontal cortex (PFC) mediate higher cognitive functions and emotional regulation that are disrupted in psychiatric disorders. The PFC undergoes significant maturation during adolescence, a period when cannabis use in humans has been linked to subsequent vulnerability to psychiatric disorders such as addiction and schizophrenia. Here, we investigated in a rat model the effects of adolescent exposure to Δ9-tetrahydrocannabinol (THC), a psychoactive component of cannabis, on the morphological architecture and transcriptional profile of layer III pyramidal neurons-using cell type- and layer-specific high-resolution microscopy, laser capture microdissection and next-generation RNA-sequencing. The results confirmed known normal expansions in basal dendritic arborization and dendritic spine pruning during the transition from late adolescence to early adulthood that were accompanied by differential expression of gene networks associated with neurodevelopment in control animals. In contrast, THC exposure disrupted the normal developmental process by inducing premature pruning of dendritic spines and allostatic atrophy of dendritic arborization in early adulthood. Surprisingly, there was minimal overlap of the developmental transcriptomes between THC- and vehicle-exposed rats. THC altered functional gene networks related to cell morphogenesis, dendritic development, and cytoskeleton organization. Marked developmental network disturbances were evident for epigenetic regulators with enhanced co-expression of chromatin- and dendrite-related genes in THC-treated animals. Dysregulated PFC co-expression networks common to both the THC-treated animals and patients with schizophrenia were enriched for cytoskeletal and neurite development. Overall, adolescent THC exposure altered the morphological and transcriptional trajectory of PFC pyramidal neurons, which could enhance vulnerability to psychiatric disorders.

Identification of rare de novo epigenetic variations in congenital disorders.

Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND-CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND-CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND-CAs, and as such likely have diagnostic relevance.

Mutation screening of MIR146A/B and BRCA1/2 3'-UTRs in the GENESIS study.

Although a wide number of breast cancer susceptibility alleles associated with various levels of risk have been identified to date, about 50% of the heritability is still missing. Although the major BRCA1 and BRCA2 genes are being extensively screened for truncating and missense variants in breast and/or ovarian cancer families, potential regulatory variants affecting their expression remain largely unexplored. In an attempt to identify such variants, we focused our attention on gene regulation mediated by microRNAs (miRs). We screened two genes, MIR146A and MIR146B, producing miR-146a and miR-146b-5p, respectively, that regulate BRCA1, and the 3'- untranslated regions (3'-UTRs) of BRCA1 and BRCA2 in the GENESIS French national case/control study (BRCA1- and BRCA2-negative breast cancer cases with at least one sister with breast cancer and matched controls). We identified one rare variant in MIR146A, four in MIR146B, five in BRCA1 3'-UTR and one in BRCA2 3'-UTR in 716 index cases and 619 controls. Among these 11 rare variants, 7 were identified each in 1 index case. None of the three relevant MIR146A/MIR146B variants affected the pre-miR sequences. The potential causality of the four relevant BRCA1/BRCA2 3'-UTRs variants was evaluated with luciferase reporter assays and co-segregation studies, as well as with bioinformatics analyses to predict miRs-binding sites, RNA secondary structures and RNA accessibility. This is the first study to report the screening of miR genes and of BRCA2 3'-UTR in a large series of familial breast cancer cases. None of the variant identified in this study gave convincing evidence of potential pathogenicity.

Occurrence of a non deleterious gene conversion event in the BRCA1 gene.

The duplication in the primate lineage of a portion of the breast and ovarian cancer susceptibility gene BRCA1 has created a BRCA1 pseudogene 45 kb away. Non-allelic homologous recombination (NAHR) between BRCA1 and BRCA1P1 has generated recurrent deleterious germ-line 37-kb deletions encompassing the first two exons of BRCA1, accounting for several breast and ovarian cancer families in various populations. In principle, NAHR intermediates resolution could also lead through a non-crossover configuration to interlocus gene conversion (IGC), but none had been described as yet. Here, we report for the first time an IGC event identified in a breast and ovarian cancer family involving exactly the same segment as that involved in the 37-kb deletions. Close examination of the consequences of this IGC event showed that it does not impact BRCA1 expression. Detailed analysis of the regions of homology between BRCA1 and its pseudogene revealed the specificity of the segment where recombination systematically occurs.

Estimation of the RNU2 macrosatellite mutation rate by BRCA1 mutation tracing.

Large tandem repeat sequences have been poorly investigated as severe technical limitations and their frequent absence from the genome reference hinder their analysis. Extensive allelotyping of this class of variation has not been possible until now and their mutational dynamics are still poorly known. In order to estimate the mutation rate of a macrosatellite, we analysed in detail the RNU2 locus, which displays at least 50 different alleles containing 5-82 copies of a 6.1 kb repeat unit. Mining data from the 1000 Genomes Project allowed us to precisely estimate copy numbers of the RNU2 repeat unit using read depth of coverage. This further revealed significantly different mean values in various recent modern human populations, favoring a scenario of fast evolution of this locus. Its proximity to a disease gene with numerous founder mutations, BRCA1, within the same linkage disequilibrium block, offered the unique opportunity to trace RNU2 arrays over a large timescale. Analysis of the transmission of RNU2 arrays associated with one 'private' mutation in an extended kindred and four founder mutations in multiple kindreds gave an estimation by maximum likelihood of 5 × 10(-3) mutations per generation, which is close to that of microsatellites.

Direct visualization of the highly polymorphic RNU2 locus in proximity to the BRCA1 gene.

Although the breast cancer susceptibility gene BRCA1 is one of the most extensively characterized genetic loci, much less is known about its upstream variable number tandem repeat element, the RNU2 locus. RNU2 encodes the U2 small nuclear RNA, an essential splicing element, but this locus is missing from the human genome assembly due to the inherent difficulty in the assembly of repetitive sequences. To fill the gap between RNU2 and BRCA1, we have reconstructed the physical map of this region by re-examining genomic clone sequences of public databases, which allowed us to precisely localize the RNU2 array 124 kb telomeric to BRCA1. We measured by performing FISH analyses on combed DNA for the first time the exact number of repeats carried by each of the two alleles in 41 individuals and found a range of 6-82 copies and a level of heterozygosity of 98%. The precise localisation of the RNU2 locus in the genome reference assembly and the implementation of a new technical tool to study it will make the detailed exploration of this locus possible. This recently neglected macrosatellite could be valuable for evaluating the potential role of structural variations in disease due to its location next to a major cancer susceptibility gene.

Calibration of multiple in silico tools for predicting pathogenicity of mismatch repair gene missense substitutions.

Classification of rare missense substitutions observed during genetic testing for patient management is a considerable problem in clinical genetics. The Bayesian integrated evaluation of unclassified variants is a solution originally developed for BRCA1/2. Here, we take a step toward an analogous system for the mismatch repair (MMR) genes (MLH1, MSH2, MSH6, and PMS2) that confer colon cancer susceptibility in Lynch syndrome by calibrating in silico tools to estimate prior probabilities of pathogenicity for MMR gene missense substitutions. A qualitative five-class classification system was developed and applied to 143 MMR missense variants. This identified 74 missense substitutions suitable for calibration. These substitutions were scored using six different in silico tools (Align-Grantham Variation Grantham Deviation, multivariate analysis of protein polymorphisms [MAPP], MutPred, PolyPhen-2.1, Sorting Intolerant From Tolerant, and Xvar), using curated MMR multiple sequence alignments where possible. The output from each tool was calibrated by regression against the classifications of the 74 missense substitutions; these calibrated outputs are interpretable as prior probabilities of pathogenicity. MAPP was the most accurate tool and MAPP + PolyPhen-2.1 provided the best-combined model (R(2)  = 0.62 and area under receiver operating characteristic = 0.93). The MAPP + PolyPhen-2.1 output is sufficiently predictive to feed as a continuous variable into the quantitative Bayesian integrated evaluation for clinical classification of MMR gene missense substitutions.