Understanding the influence of gene expression on the molecular mechanisms underpinning human phenotypic diversity is fundamental to being able to predict health outcomes and treat disease. We have carried out whole transcriptome expression analysis on a series of eight normal human postmortem eyes by RNA sequencing. Here we present data showing that ∼80% of the transcriptome is expressed in the posterior layers of the eye and that there is significant differential expression not only between the layers of the posterior part of the eye but also between locations of a tissue layer. These differences in expression also extend to alternative splicing and splicing factors. Differentially expressed genes are enriched for genes associated with psychiatric, immune and cardiovascular disorders. Enrichment categories for gene ontology included ion transport, synaptic transmission and visual and sensory perception. Lastly, allele-specific expression was found to be significant forCFH,C3 andCFB, which are known risk genes for age-related macular degeneration. These expression differences should be useful in determining the underlying biology of associations with common diseases of the human retina, retinal pigment epithelium and choroid and in guiding the analysis of the genomic regions involved in the control of normal gene expression.

Central conducting lymphatic anomaly (CCLA) is one of the complex lymphatic anomalies characterized by dilated lymphatic channels, lymphatic channel dysmotility and distal obstruction affecting lymphatic drainage. We performed whole exome sequencing (WES) of DNA from a four-generation pedigree and examined the consequences of the variant by transfection of mammalian cells and morpholino and rescue studies in zebrafish. WES revealed a heterozygous mutation in EPHB4 (RefSeq NM_004444.4; c.2334 + 1G>C) and RNA-Seq demonstrated that the EPHB4 mutation destroys the normal donor site, which leads to the use of a cryptic splice donor that results in retention of the intervening 12-bp intron sequence. Transient co-expression of the wild-type and mutant EPHB4 proteins showed reduced phosphorylation of tyrosine, consistent with a loss-of-function effect. Zebrafish ephb4a morpholino resulted in vessel misbranching and deformities in the lymphatic vessel development, indicative of possible differentiation defects in lymphatic vessels, mimicking the lymphatic presentations of the patients. Immunoblot analysis using zebrafish lysates demonstrated over-activation of mTORC1 as a consequence of reduced EPHB4 signaling. Strikingly, drugs that inhibit mTOR signaling or RAS-MAPK signaling effectively rescued the misbranching phenotype in a comparable manner. Moreover, knock-in of EPHB4 mutation in HEK293T cells also induced mTORC1 activity. Our data demonstrate the pathogenicity of the identified EPHB4 mutation as a novel cause of CCLA and suggesting that ERK inhibitors may have therapeutic benefits in such patients with complex lymphatic anomalies.

Recent advances in RNA sequencing (RNA-Seq) have enabled the discovery of novel transcriptomic variations that are not possible with traditional microarray-based methods. Tissue and cell specific transcriptome changes during pathophysiological stress in disease cases versus controls and in response to therapies are of particular interest to investigators studying cardiometabolic diseases. Thus, knowledge on the relationships between sequencing depth and detection of transcriptomic variation is needed for designing RNA-Seq experiments and for interpreting results of analyses. Using deeply sequenced Illumina HiSeq 2000 101 bp paired-end RNA-Seq data derived from adipose of a healthy individual before and after systemic administration of endotoxin (LPS), we investigated the sequencing depths needed for studies of gene expression and alternative splicing (AS). In order to detect expressed genes and AS events, we found that ∼100 to 150 million (M) filtered reads were needed. However, the requirement on sequencing depth for the detection of LPS modulated differential expression (DE) and differential alternative splicing (DAS) was much higher. To detect 80% of events, ∼300 M filtered reads were needed for DE analysis whereas at least 400 M filtered reads were necessary for detecting DAS. Although the majority of expressed genes and AS events can be detected with modest sequencing depths (∼100 M filtered reads), the estimated gene expression levels and exon/intron inclusion levels were less accurate. We report the first study that evaluates the relationship between RNA-Seq depth and the ability to detect DE and DAS in human adipose. Our results suggest that a much higher sequencing depth is needed to reliably identify DAS events than for DE genes.

Objective—: Inappropriate transcriptional activation of innate immunity is a pathological feature of several cardiometabolic disorders, but little is known about inflammatory modulation of long intergenic noncoding RNAs (lincRNAs) in disease-relevant human tissues. Approach and Results—: We applied deep RNA sequencing (>500 million filtered reads per sample) to blood and adipose during low-dose experimental endotoxemia (lipopolysaccharide) in a healthy human, with targeted replication in separate individuals undergoing endotoxemia (n=6), to identify inflammatory lincRNAs. A subset of these lincRNAs was examined for expression in adipocytes and monocytes, modulation in adipose of obese humans, and overlap with genome-wide association study signals for inflammatory and cardiometabolic traits. Of a stringent set of 4284 lincRNAs, ≈11% to 22% were expressed with 201 and 56 lincRNAs modulated by lipopolysaccharide in blood or adipose, respectively. Tissue-specific expression of a subset of 6 lipopolysaccharide-lincRNAs was replicated with lipopolysaccharide modulation confirmed for all 3 expressed in blood and 2 of 4 expressed in adipose. The broader generalizability of findings in blood of subject A was confirmed by RNA sequencing in 7 additional subjects. We confirmed adipocytes and monocytes as potential cell-sources of selective lipopolysaccharide-regulated lincRNAs, and 2 of these, linc-DMRT2 ( P =0.002) and linc-TP53I13 ( P =0.01), were suppressed in adipose of obese humans. Finally, we provide examples of lipopolysaccharide-modulated lincRNAs that overlap single nucleotide polymorphisms that are associated with cardiometabolic traits. Conclusions—: Our findings provide novel insights into tissue-level, inflammatory transcriptome regulation in cardiometabolic diseases. These are complementary to more usual approaches limited to interrogation of DNA variations.