Widespread sex differences in gene expression and splicing in the adult human brain

There is strong evidence to show that men and women differ in terms of neurodevelopment, neurochemistry and susceptibility to neurodegenerative and neuropsychiatric disease. The molecular basis of these differences remains unclear. Progress in this field has been hampered by the lack of genome-wide information on sex differences in gene expression and in particular splicing in the human brain. Here we address this issue by using post-mortem adult human brain and spinal cord samples originating from 137 neuropathologically confirmed control individuals to study whole-genome gene expression and splicing in 12 CNS regions. We show that sex differences in gene expression and splicing are widespread in adult human brain, being detectable in all major brain regions and involving 2.5% of all expressed genes. We give examples of genes where sex-biased expression is both disease-relevant and likely to have functional consequences, and provide evidence suggesting that sex biases in expression may reflect sex-biased gene regulatory structures

Of Humans and Gerbils— Independent Diversification of Neuroligin-4 Into X- and Y-Specific Genes in Primates and Rodents

The neural cell adhesion protein neuroligin-4 has puzzled neuroscientists and geneticist alike for almost two decades. Its clinical association with autism spectrum disorders (ASD) is well established, however, its diversification into sex chromosome-specific copies, NLGN4X and NLGN4Y, remains uncharted territory. Just recently, the presence of substantial neuroligin-4 sequence differences between humans and laboratory mice, in which Nlgn4 is a pseudoautosomal gene, could be explained as a consequence of dramatic changes affecting the pseudoautosomal region on both sex chromosomes in a subset of rodents, the clade eumuroida. In this study, we describe the presence of sex chromosome-specific copies of neuroligin-4 genes in the Mongolian gerbil (Meriones unguiculatus) marking the first encounter of its kind in rodents. Gerbils are members of the family Muridae and are closely related to mice and rats. Our results have been incorporated into an extended evolutionary analysis covering primates, rodents, lagomorphs, treeshrews and culogos comprising together the mammalian superorder euarchontoglires. We gathered evidence that substantial changes in neuroligin-4 genes have also occurred outside eumuroida in other rodent species as well as in lagomorphs. These changes feature, e.g., a general reduction of its gene size, an increase in its average GC-content as well as in the third position (GC3) of synonymous codons, and the accumulation of repetitive sequences in line with previous observations. We further show conclusively that the diversification of neuroligin-4 in sex chromosome-specific copies has happened multiple times independently during mammal evolution proving that Y-chromosomal NLGN4Y genes do not originate from a single common NLGN4Y ancestor

A paradigmatic autistic phenotype associated with loss of PCDH11Y and NLGN4Y genes

Background: Most studies relative to Y chromosome abnormalities are focused on the sexual developmental disorders. Recently, a few studies suggest that some genes located on Y chromosome may be related to different neurodevelopment disorders. Case presentation: We report a child with sexual developmental disorder associated with a peculiar phenotype characterized by severe language impairment and autistic behaviour associated with a mosaicism [45,X(11)/46,XY(89)] and a partial deletion of the short and long arm of Y chromosome (del Yp11.31q11.23) that also involves the loss of both PCDH11Y and NLGN4Y genes. To our knowledge no study has ever reported the occurrence of the lack of both PCDH11Y and NLGN4Y located in the Y chromosome in the same patient. Conclusions: We hypothesized a functional complementary role of PCDH11Y and NLGN4Y within formation/maturation of the cerebral cortex. The impairment of early language development may be mainly related to the lack of PCDH11Y that underlies the early language network development and the later appearance of the autistic behaviour may be mainly related to deficit of inhibitory glicinergic neurotransmission NLGN4Y-linked

Novel human sex-typing strategies based on the autism candidate gene NLGN4X and its male-specific gametologue NLGN4Y

Background: Since the early days of PCR techniques, sex identification, “sex-typing,” of genomic DNA samples has been a fundamental part of human forensic analysis but also in animal genetics aiming at strategic livestock breeding. Most analyses are employing the AMELX/AMELY gene loci on the X and Y chromosomes present in most mammals. We hypothesize that sex-typing in humans is also possible based on the genes NLGN4X and NLGN4Y, which represent X and Y chromosome-specific copies of a common ancestral neuroligin-4 orthologue. Methods: Genomic DNA was isolated from human blood and buccal cell samples (total n = 111) and submitted to two different strategies: (a) a traditional two-primer PCR approach detecting an insertion/deletion (indel) polymorphism immediately upstream of the translational start on exon 1 and (b) detection of a single nucleotide polymorphism, SNP, on the translational stop carrying exon 7. The SNP detection was based on a quantitative PCR approach (rhAMP genotyping) employing DNA/RNA hybrid oligonucleotides that were blocked and which could only be activated upon perfect annealing to the target DNA sequence. Results: All indel PCR-tested human DNA samples showed two bands for males representing X- and Y-specific copies of NLGN4 and a single band for female samples, i.e., homozygosity of NLGN4X and absence of NLGN4Y, in accordance with the self-reported sex of the donors. These results were in perfect agreement with the results of the rhAMP-based SNP-detection method: all males were consequently positive for both alleles, representing either SNP variant, and females were interpreted as homozygous regarding the SNP variant found in NLGN4X. Both methods have shown reliable and consistent results that enabled us to infer the sex of donor DNA samples across different ethnicities. Conclusions: These results indicate that the detection of human NLGN4X/Y is a suitable alternative to previously reported methods employing gene loci such as AMELX/Y. Furthermore, this is the first report applying successfully the rhAMP-genotyping strategy as a means for SNP-based sex-typing, which consequently will be applicable to other gene loci or different species as well

A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments

Background: High-throughput RNA sequencing (RNA-seq) offers unprecedented power to capture the real dynamics of gene expression. Experimental designs with extensive biological replication present a unique opportunity to exploit this feature and distinguish expression profiles with higher resolution. RNA-seq data analysis methods so far have been mostly applied to data sets with few replicates and their default settings try to provide the best performance under this constraint. These methods are based on two well-known count data distributions: the Poisson and the negative binomial. The way to properly calibrate them with large RNA-seq data sets is not trivial for the non-expert bioinformatics user. Results: Here we show that expression profiles produced by extensively-replicated RNA-seq experiments lead to a rich diversity of count data distributions beyond the Poisson and the negative binomial, such as Poisson-Inverse Gaussian or Pólya-Aeppli, which can be captured by a more general family of count data distributions called the Poisson-Tweedie. The flexibility of the Poisson-Tweedie family enables a direct fitting of emerging features of large expression profiles, such as heavy-tails or zero-inflation, without the need to alter a single configuration parameter. We provide a software package for R called tweeDEseq implementing a new test for differential expression based on the Poisson-Tweedie family. Using simulations on synthetic and real RNA-seq data we show that tweeDEseq yields P-values that are equally or more accurate than competing methods under different configuration parameters. By surveying the tiny fraction of sex-specific gene expression changes in human lymphoblastoid cell lines, we also show that tweeDEseq accurately detects differentially expressed genes in a real large RNA-seq data set with improved performance and reproducibility over the previously compared methodologies. Finally, we compared the results with those obtained from microarrays in order to check for reproducibility. Conclusions: RNA-seq data with many replicates leads to a handful of count data distributions which can be accurately estimated with the statistical model illustrated in this paper. This method provides a better fit to the underlying biological variability; this may be critical when comparing groups of RNA-seq samples with markedly different count data distributions. The tweeDEseq package forms part of the Bioconductor project and it is available for download at http://www.bioconductor.or

Evidence of novel finescale structural variation at autism spectrum disorder candidate loci

Background: Autism spectrum disorders (ASD) represent a group of neurodevelopmental disorders characterized by a core set of social-communicative and behavioral impairments. Gamma-aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the brain, acting primarily via the GABA receptors (GABR). Multiple lines of evidence, including altered GABA and GABA receptor expression in autistic patients, indicate that the GABAergic system may be involved in the etiology of autism. Methods: As copy number variations (CNVs), particularly rare and de novo CNVs, have now been implicated in ASD risk, we examined the GABA receptors and genes in related pathways for structural variation that may be associated with autism. We further extended our candidate gene set to include 19 genes and regions that had either been directly implicated in the autism literature or were directly related (via function or ancestry) to these primary candidates. For the high resolution CNV screen we employed custom-designed 244 k comparative genomic hybridization (CGH) arrays. Collectively, our probes spanned a total of 11 Mb of GABA-related and additional candidate regions with a density of approximately one probe every 200 nucleotides, allowing a theoretical resolution for detection of CNVs of approximately 1 kb or greater on average. One hundred and sixty-eight autism cases and 149 control individuals were screened for structural variants. Prioritized CNV events were confirmed using quantitative PCR, and confirmed loci were evaluated on an additional set of 170 cases and 170 control individuals that were not included in the original discovery set. Loci that remained interesting were subsequently screened via quantitative PCR on an additional set of 755 cases and 1,809 unaffected family members. Results: Results include rare deletions in autistic individuals at JAKMIP1, NRXN1, Neuroligin4Y, OXTR, and ABAT. Common insertion/deletion polymorphisms were detected at several loci, including GABBR2 and NRXN3. Overall, statistically significant enrichment in affected vs. unaffected individuals was observed for NRXN1 deletions. Conclusions: These results provide additional support for the role of rare structural variation in ASD

Gene expression in human brain implicates sexually dimorphic pathways in autism spectrum disorders.

Autism spectrum disorder (ASD) is more prevalent in males, and the mechanisms behind this sex-differential risk are not fully understood. Two competing, but not mutually exclusive, hypotheses are that ASD risk genes are sex-differentially regulated, or alternatively, that they interact with characteristic sexually dimorphic pathways. Here we characterized sexually dimorphic gene expression in multiple data sets from neurotypical adult and prenatal human neocortical tissue, and evaluated ASD risk genes for evidence of sex-biased expression. We find no evidence for systematic sex-differential expression of ASD risk genes. Instead, we observe that genes expressed at higher levels in males are significantly enriched for genes upregulated in post-mortem autistic brain, including astrocyte and microglia markers. This suggests that it is not sex-differential regulation of ASD risk genes, but rather naturally occurring sexually dimorphic processes, potentially including neuron-glial interactions, that modulate the impact of risk variants and contribute to the sex-skewed prevalence of ASD

Detection of 1p19q Deletion by Real-Time Comparative Quantitative PCR

1p/19q (1p and/or 19q) deletions are prognostic factors in oligodendroglial tumors (OT) and predict better survival after both chemotherapy and radiotherapy. While studying 1p/19q status as a potential variable within multivariate prognosis models for OT, we have frequently encountered unknown 1p/19q status within our glioma sample database due to lack of paired blood samples for loss of heterozygosity (LOH) assay and/or failure to perform fluorescence in situ hybridization (FISH). We realized that a 1p and 19q deletion assay that could be reliably performed solely on tumor DNA samples would allow us to fill in these molecular biology data “holes”. We built recombinant DNA with fragments of the selected “marker” genes in 1p (E2F2, NOTCH2), and 19q (PLAUR) and “reference” genes (ERC2, SPOCK1, and SPAG16 ) and used it as quantification standard in real-time PCR to gain absolute ratios of marker/reference gene copy numbers in tumor DNA samples, thus called comparative quantitative PCR (CQ-PCR). Using CQ-PCR, we identified 1p and/ or 19q deletions in majority of pure low-grade oligodenroglioma (OG) tumors (17/21, 81%), a large portion of anaplastic oligodendroglioma (AO) tumors (6/15, 47%), but rarely found in mixed oligoastrcytomas (OA) tumors (1/8, 13%). These data are consistent with results of LOH and FISH assays generally reported for these tumor types. In addition, 15 out 18 samples showed concordant results between FISH and CQ-PCR. We conclude that CQ-PCR is a potential means to gain 1p/19q deletion information, which prognostic and predictive values of CQ-PCR-derived 1p/19q status will be determined in a future study

Experiences with array-based sequence capture; toward clinical applications

Although sequencing of a human genome gradually becomes an option, zooming in on the region of interest remains attractive and cost saving. We performed array-based sequence capture using 385K Roche NimbleGen, Inc. arrays to zoom in on the protein-coding and immediate intron-flanking sequences of 112 genes, potentially involved in mental retardation and congenital malformation. Captured material was sequenced using Illumina technology. A data analysis pipeline was built that detects sequence variants, positions them in relation to the gene, checks for presence in databases (eg, db single-nucleotide polymorphism (SNP)) and predicts the potential consequences at the level of RNA splicing and protein translation. In the samples analyzed, all known variants were reliably detected, including pathogenic variants from control cases and SNPs derived from array experiments. Although overall coverage varied considerably, it was reproducible per region and facilitated the detection of large deletions and duplications (copy number variations), including a partial deletion in the B3GALTL gene from a patient sample. For ultimate diagnostic application, overall results need to be improved. Future arrays should contain probes from both DNA strands, and to obtain a more even coverage, one could add fewer probes from densely and more probes from sparsely covered regions