RNA sequencing of transcriptomes in human brain regions: protein-coding and non-coding RNAs, isoforms and alleles.

BACKGROUND: We used RNA sequencing to analyze transcript profiles of ten autopsy brain regions from ten subjects. RNA sequencing techniques were designed to detect both coding and non-coding RNA, splice isoform composition, and allelic expression. Brain regions were selected from five subjects with a documented history of smoking and five non-smokers. Paired-end RNA sequencing was performed on SOLiD instruments to a depth of >40 million reads, using linearly amplified, ribosomally depleted RNA. Sequencing libraries were prepared with both poly-dT and random hexamer primers to detect all RNA classes, including long non-coding (lncRNA), intronic and intergenic transcripts, and transcripts lacking poly-A tails, providing additional data not previously available. The study was designed to generate a database of the complete transcriptomes in brain region for gene network analyses and discovery of regulatory variants. RESULTS: Of 20,318 protein coding and 18,080 lncRNA genes annotated from GENCODE and lncipedia, 12 thousand protein coding and 2 thousand lncRNA transcripts were detectable at a conservative threshold. Of the aligned reads, 52 % were exonic, 34 % intronic and 14 % intergenic. A majority of protein coding genes (65 %) was expressed in all regions, whereas ncRNAs displayed a more restricted distribution. Profiles of RNA isoforms varied across brain regions and subjects at multiple gene loci, with neurexin 3 (NRXN3) a prominent example. Allelic RNA ratios deviating from unity were identified in > 400 genes, detectable in both protein-coding and non-coding genes, indicating the presence of cis-acting regulatory variants. Mathematical modeling was used to identify RNAs stably expressed in all brain regions (serving as potential markers for normalizing expression levels), linked to basic cellular functions. An initial analysis of differential expression analysis between smokers and nonsmokers implicated a number of genes, several previously associated with nicotine exposure. CONCLUSIONS: RNA sequencing identifies distinct and consistent differences in gene expression between brain regions, with non-coding RNA displaying greater diversity between brain regions than mRNAs. Numerous RNAs exhibit robust allele selective expression, proving a means for discovery of cis-acting regulatory factors with potential clinical relevance.

Delineation of candidate genes responsible for structural brain abnormalities in patients with terminal deletions of chromosome 6q27.

Patients with terminal deletions of chromosome 6q present with structural brain abnormalities including agenesis of corpus callosum, hydrocephalus, periventricular nodular heterotopia, and cerebellar malformations. The 6q27 region harbors genes that are important for the normal development of brain and delineation of a critical deletion region for structural brain abnormalities may lead to a better genotype-phenotype correlation. We conducted a detailed clinical and molecular characterization of seven unrelated patients with deletions involving chromosome 6q27. All patients had structural brain abnormalities. Using array comparative genomic hybridization, we mapped the size, extent, and genomic content of these deletions. The smallest region of overlap spans 1.7 Mb and contains DLL1, THBS2, PHF10, and C6orf70 (ERMARD) that are plausible candidates for the causation of structural brain abnormalities. Our study reiterates the importance of 6q27 region in normal development of brain and helps identify putative genes in causation of structural brain anomalies.

Genetic variants affecting alternative splicing of human cholesteryl ester transfer protein.

Cholesteryl ester transfer protein (CETP) plays an important role in reverse cholesterol transport, with decreased CETP activity increasing HDL levels. Formation of an alternative splice form lacking exon 9 (Δ9-CETP) has been associated with two single nucleotide polymorphisms (SNPs) in high linkage disequilibrium with each other, namely rs9930761 T>C located in intron 8 in a putative splicing branch site and rs5883 C>T in a possible exonic splicing enhancer (ESE) site in exon 9. To assess the relative effect of rs9930761 and rs5883 on splicing, mini-gene constructs spanning CETP exons 8 to 10, carrying all four possible allele combinations, were transfected into HEK293 and HepG2 cells. The minor T allele of rs5883 enhanced splicing significantly in both cell lines whereas the minor C allele of rs9930761 did not. In combination, the two alleles did not yield greater splicing than the rs5883 T allele alone in HepG2 cells. These results indicate that the genetic effect on CETP splicing is largely attributable to rs5883. We also confirm that Δ9-CETP protein is expressed in the liver but fails to circulate in the blood.

Whole transcriptome RNA-Seq allelic expression in human brain.

BACKGROUND: Measuring allelic RNA expression ratios is a powerful approach for detecting cis-acting regulatory variants, RNA editing, loss of heterozygosity in cancer, copy number variation, and allele-specific epigenetic gene silencing. Whole transcriptome RNA sequencing (RNA-Seq) has emerged as a genome-wide tool for identifying allelic expression imbalance (AEI), but numerous factors bias allelic RNA ratio measurements. Here, we compare RNA-Seq allelic ratios measured in nine different human brain regions with a highly sensitive and accurate SNaPshot measure of allelic RNA ratios, identifying factors affecting reliable allelic ratio measurement. Accounting for these factors, we subsequently surveyed the variability of RNA editing across brain regions and across individuals. RESULTS: We find that RNA-Seq allelic ratios from standard alignment methods correlate poorly with SNaPshot, but applying alternative alignment strategies and correcting for observed biases significantly improves correlations. Deploying these methods on a transcriptome-wide basis in nine brain regions from a single individual, we identified genes with AEI across all regions (SLC1A3, NHP2L1) and many others with region-specific AEI. In dorsolateral prefrontal cortex (DLPFC) tissues from 14 individuals, we found evidence for frequent regulatory variants affecting RNA expression in tens to hundreds of genes, depending on stringency for assigning AEI. Further, we find that the extent and variability of RNA editing is similar across brain regions and across individuals. CONCLUSIONS: These results identify critical factors affecting allelic ratios measured by RNA-Seq and provide a foundation for using this technology to screen allelic RNA expression on a transcriptome-wide basis. Using this technology as a screening tool reveals tens to hundreds of genes harboring frequent functional variants affecting RNA expression in the human brain. With respect to RNA editing, the similarities within and between individuals leads us to conclude that this post-transcriptional process is under heavy regulatory influence to maintain an optimal degree of editing for normal biological function.

Cholesteryl Ester Transfer Protein (CETP) polymorphisms affect mRNA splicing, HDL levels, and sex-dependent cardiovascular risk.

Polymorphisms in and around the Cholesteryl Ester Transfer Protein (CETP) gene have been associated with HDL levels, risk for coronary artery disease (CAD), and response to therapy. The mechanism of action of these polymorphisms has yet to be defined. We used mRNA allelic expression and splice isoform measurements in human liver tissues to identify the genetic variants affecting CETP levels. Allelic CETP mRNA expression ratios in 56 human livers were strongly associated with several variants 2.5-7 kb upstream of the transcription start site (e.g., rs247616 p = 6.4 × 10(-5), allele frequency 33%). In addition, a common alternatively spliced CETP isoform lacking exon 9 (Δ9), has been shown to prevent CETP secretion in a dominant-negative manner. The Δ 9 expression ranged from 10 to 48% of total CETP mRNA in 94 livers. Increased formation of this isoform was exclusively associated with an exon 9 polymorphism rs5883-C>T (p = 6.8 × 10(-10)) and intron 8 polymorphism rs9930761-T>C (5.6 × 10(-8)) (in high linkage disequilibrium with allele frequencies 6-7%). rs9930761 changes a key splicing branch point nucleotide in intron 8, while rs5883 alters an exonic splicing enhancer sequence in exon 9.The effect of these polymorphisms was evaluated in two clinical studies. In the Whitehall II study of 4745 subjects, both rs247616 and rs5883T/rs9930761C were independently associated with increased HDL-C levels in males with similar effect size (rs247616 p = 9.6 × 10(-28) and rs5883 p = 8.6 × 10(-10), adjusted for rs247616). In an independent multiethnic US cohort of hypertensive subjects with CAD (INVEST-GENE), rs5883T/rs9930761C alone were significantly associated with increased incidence of MI, stroke, and all-cause mortality in males (rs5883: OR 2.36 (CI 1.29-4.30), p = 0.005, n = 866). These variants did not reach significance in females in either study. Similar to earlier results linking low CETP activity with poor outcomes in males, our results suggest genetic, sex-dependent CETP splicing effects on cardiovascular risk by a mechanism independent of circulating HDL-C levels.

Distinct cellular pools of perilipin 5 point to roles in lipid trafficking.

The PAT family of lipid storage droplet proteins comprised five members, each of which has become an established regulator of cellular neutral lipid metabolism. Perilipin 5 (also known as lsdp-5, MLDP, PAT-1, and OXPAT), the most recently discovered member of the family, has been shown to localize to two distinct intracellular pools: the lipid storage droplet (LD), and a poorly characterized cytosolic fraction. We have characterized the denser of these intracellular pools and find that a population of perilipin 5 not associated with large LDs resides in complexes with a discrete density (~1.15 g/ml) and size (~575 kDa). Using immunofluorescence, western blotting of isolated sucrose density fractions, native gradient gel electrophoresis, and co-immunoprecipitation, we have shown that these small (~15 nm), perilipin 5-encoated structures do not contain the PAT protein perilipin 2 (ADRP), but do contain perilipin 3 and several other as of yet uncharacterized proteins. The size and density of these particles as well as their susceptibility to degradation by lipases suggest that like larger LDs, they have a neutral lipid rich core. When treated with oleic acid to promote neutral lipid deposition, cells ectopically expressing perilipin 5 experienced a reorganization of LDs in the cell, resulting in fewer, larger droplets at the expense of smaller ones. Collectively, these data demonstrate that a portion of cytosolic perilipin 5 resides in high density lipid droplet complexes that participate in cellular neutral lipid accumulation.

Aflatoxicosis in nine dogs after exposure to contaminated commercial dog food.

The purpose of this study was to characterize light and electron microscopic findings from 9 dogs that had consumed aflatoxin-contaminated commercial dog food from recalled batches. Four dogs died and 5 were euthanized after signs of liver failure. Analysis of feed and liver samples confirmed exposure to aflatoxin. Of the 9 dogs, 8 had classic signs of liver failure, and 1 had signs of liver failure. Enlarged, pale yellow livers were seen macroscopically at necropsy in the dogs with subacute hepatopathy, and cirrhosis was noted in the dog with chronic hepatopathy. Histopathologic findings included hepatic lipidosis, portal fibroplasia, and biliary hyperplasia, which supported a diagnosis of subacute toxic hepatopathy in the 8 symptomatic animals. Marked lobular atrophy, bridging portal fibrosis, and regenerative hepatocellular nodules characterized the dog with chronic hepatopathy. Electron microscopy revealed marked hepatocellular lipid vacuolation and early fibroplasia in the dogs with acute hepatopathy and marked fibrosis and regeneration in the dog with chronic hepatopathy. Analysis of feed for aflatoxin consistently revealed high levels of aflatoxin B1 (range of 223-579 ppb), and hepatic tissue contained elevated levels of aflatoxin B1 metabolite M1 (0.6-4.4 ppb). Although dogs are not commonly affected by aflatoxicosis, they are highly susceptible and can present with classic signs of acute or chronic hepatopathy. Characteristic gross, histologic, and electron microscopic changes help pathologists determine a presumptive toxic insult. Detecting aflatoxins or their metabolites in feed or liver specimens can help confirm the diagnosis of aflatoxicosis.

Characterization of mu, kappa, and delta opioid binding in amphibian whole brain tissue homogenates.

Opioid agonists produce analgesia in mammals through the activation of mu, kappa, or delta opioid receptors. Previous behavioral and binding studies from our laboratory using an amphibian model suggested that mu, kappa, or delta opioid agonists may activate a single type of opioid receptor in the grass frog, Rana pipiens. In the present study, kinetic, saturation, and competitive binding profiles for three opioid radioligands, [(3)H]DAMGO ([D-Ala(2),N-Me-Phe(4),Gly(5)-ol]-enkephalin) (mu-selective), [(3)H]U65953 [(5 alpha, 7 alpha,8 beta)-(+)-N-methyl-N-[7-(1-pyrrolidinyl)-1-oxaspiro[4.5]dec-8-yl]-benzeneacetamide] (kappa-selective), and [(3)H]DPDPE ([D-Pen(2),D-Pen(5)]-enkephalin) (delta-selective) were determined using frog whole brain homogenates. Kinetic analyses and experimentally derived values from saturation experiments gave affinity constants (K(D)) in the low nanomolar range. The density of opioid binding sites (B(max)) was 224.4, 118.6, and 268.9 fmol/mg for mu, kappa, and delta opioid radioligands, respectively. The affinity values did not significantly differ among the three opioid radioligands, but the kappa radioligand bound to significantly fewer sites than did the mu or delta radioligands. K(i) values for unlabeled mu, kappa, and delta competitors, including highly selective opioid antagonists, were consistent with each radioligand selectivity profile. The present data suggest that mu, kappa, and delta opioid radioligands bind to distinct opioid receptors in amphibians that are surprisingly similar to those found in mammalian brain.