Identifying breast cancer risk loci by global differential allele-specific expression (DASE) analysis in mammary epithelial transcriptome

Identification of candidate causal variants and target genes at 41 breast cancer risk loci through differential allelic expression analysis

Understanding breast cancer genetic risk relies on identifying causal variants and candidate target genes in risk loci identified by genome-wide association studies (GWAS), which remains challenging. Since most loci fall in active gene regulatory regions, we developed a novel approach facilitated by pinpointing the variants with greater regulatory potential in the disease's tissue of origin. Through genome-wide differential allelic expression (DAE) analysis, using microarray data from 64 normal breast tissue samples, we mapped the variants associated with DAE (daeQTLs). Then, we intersected these with GWAS data to reveal candidate risk regulatory variants and analysed their cis-acting regulatory potential. Finally, we validated our approach by extensive functional analysis of the 5q14.1 breast cancer risk locus. We observed widespread gene expression regulation by cis-acting variants in breast tissue, with 65% of coding and noncoding expressed genes displaying DAE (daeGenes). We identified over 54 K daeQTLs for 6761 (26%) daeGenes, including 385 daeGenes harbouring variants previously associated with BC risk. We found 1431 daeQTLs mapped to 93 different loci in strong linkage disequilibrium with risk-associated variants (risk-daeQTLs), suggesting a link between risk-causing variants and cis-regulation. There were 122 risk-daeQTL with stronger cis-acting potential in active regulatory regions with protein binding evidence. These variants mapped to 41 risk loci, of which 29 had no previous report of target genes and were candidates for regulating the expression levels of 65 genes. As validation, we identified and functionally characterised five candidate causal variants at the 5q14.1 risk locus targeting the ATG10 and ATP6AP1L genes, likely acting via modulation of alternative transcription and transcription factor binding. Our study demonstrates the power of DAE analysis and daeQTL mapping to identify causal regulatory variants and target genes at breast cancer risk loci, including those with complex regulatory landscapes. It additionally provides a genome-wide resource of variants associated with DAE for future functional studies.

Caveolin-1, a Determinant of the Fate of MCF-7 Breast Cancer Cells

Background

The scaffolding protein, caveolin-1 (Cav-1), participates in multiple cellular functions including promotion of sodium excretion from the kidney. Loss of expression of Cav-1 is associated with tumorigenesis of various types of cancer. We have shown the potential link between hypertension and breast cancer via abnormal function of the G protein-coupled receptor kinase type 4 (GRK4).

Objective

The current studies tested the hypothesis that Cav-1 acts as a tumor-suppressive factor in breast cancer cells and enhances the sensitivity to the inhibitory effect of the type 1 dopaminergic receptor (D1R).

Methods

Michigan Cancer Foundation (MCF) MCF-7 cells stably expressing a Cav-1/mCherry fusion protein or mCherry alone were used as models to examine the effect of Cav-1 on cell growth, apoptosis, and senescence. Cell proliferation was determined by cell counting, cell cycle analysis (flow cytometry), and BrdU incorporation. Apoptosis was determined using the Cell Death Detection ELISA kit from Roche Diagnosis. Senescence was determined using the senescence associated beta galactosidase (SA-β-gal) assay. Reactive oxygen species (ROS) was measured using 2',7'-dichlorodihydrofluorescein diacetate. Western blot analysis was used to measure activation of signaling pathway molecules. All statistical analyses were conducted with Microsoft Excel.

Results

Overexpression of Cav-1 in MCF-7 cells reduced cellular growth rate. Both inhibition of proliferation and induction of cell death are contributing factors. Multiple signaling pathways were activated in Cav-1-expressing MCF-7 cells. Activation of Akt was prominent. In MCF-7-expressing Cav-1 (MCF-7 Cav-1) cells, the levels of phosphorylated Akt at S473 and T308 were increased 28- and 8.7-fold, respectively. Instead of protecting cells from apoptosis, extremely high levels of activated Akt resulted in increased levels of ROS which led to apoptosis and senescence. The tumor-suppressive effect plus downregulation of GRK4 makes Cav-1-expressing MCF-7 cells significantly more sensitive to the inhibitory effect of the D1R agonist, SKF38393.

Conclusion

Caveolin-1 acts as a tumor-suppressing factor via extreme activation of Akt and down regulation of survival factors such as GRK4, survivin, and cyclin D1.

Differential Allele-Specific Expression Revealed Functional Variants and Candidate Genes Related to Meat Quality Traits in B. indicus Muscle

Traditional transcriptomics approaches have been used to identify candidate genes affecting economically important livestock traits. Regulatory variants affecting these traits, however, remain under covered. Genomic regions showing allele-specific expression (ASE) are under the effect of cis-regulatory variants, being useful for improving the accuracy of genomic selection models. Taking advantage of the better of these two methods, we investigated single nucleotide polymorphisms (SNPs) in regions showing differential ASE (DASE SNPs) between contrasting groups for beef quality traits. For these analyses, we used RNA sequencing data, imputed genotypes and genomic estimated breeding values of muscle-related traits from 190 Nelore (Bos indicus) steers. We selected 40 contrasting unrelated samples for the analysis (N = 20 animals per contrasting group) and used a beta-binomial model to identify ASE SNPs in only one group (i.e., DASE SNPs). We found 1479 DASE SNPs (FDR ≤ 0.05) associated with 55 beef-quality traits. Most DASE genes were involved with tenderness and muscle homeostasis, presenting a co-expression module enriched for the protein ubiquitination process. The results overlapped with epigenetics and phenotype-associated data, suggesting that DASE SNPs are potentially linked to cis-regulatory variants affecting simultaneously the transcription and phenotype through chromatin state modulation.

Ubiquitin-specific protease TRE17/USP6 promotes tumor cell invasion through the regulation of glycoprotein CD147 intracellular trafficking

Disordered expression and distribution of plasma membrane proteins at the cell surface leads to diverse malignant phenotypes in tumors, including cell invasion. The ubiquitin-specific protease TRE17/USP6, an oncogene identified in Ewing sarcoma, is highly expressed in several cancers and locally aggressive tumor-like lesions. We have previously demonstrated that TRE17 regulates the trafficking of plasma membrane proteins that enter cells via clathrin-independent endocytosis (CIE); TRE17 prevents CIE cargo proteins from being targeted to lysosomes for degradation by deubiquitylating them. However, functional insights into the effects of TRE17-mediated CIE cargo trafficking on cell invasion remain unknown. Here, we show that increased expression of TRE17 enhances invasiveness of the human sarcoma cell line HT-1080 by elevating the cell surface levels of the membrane glycoprotein CD147, which plays a central role in tumor progression. We demonstrate overexpression of TRE17 decreases ubiquitylated CD147, which is accompanied by suppression of CD147 transport to lysosomes, resulting in the stabilization and increase of cell surface-localized CD147. On the other hand, we show knockdown of TRE17 decreases cell surface CD147, which is coupled with reduced production of matrix metalloproteinases (MMPs), the enzymes responsible for extracellular matrix degradation. Furthermore, we demonstrate that inhibition of CD147 by a specific inhibitor alleviated the TRE17-promoted tumor cell invasion. We therefore propose a model for the pathogenesis of TRE17-driven tumors in which TRE17 increases CD147 at the cell surface by preventing its lysosomal degradation, which in turn enhances MMP synthesis and matrix degradation, thereby promoting tumor cell invasion.

GRK4, A Potential Link between Hypertension and Breast Cancer

Hypertension and breast cancer are two common diseases occurring in women. Clinical studies have shown increased breast cancer incidence in hypertensive women. Several lines of evidence demonstrate that G protein-coupled Receptor Kinase 4 (GRK4) could be a common risk factor for hypertension and breast cancer. This article reviews our current understanding of molecular mechanisms of GRK4 in hypertension and breast cancer.

Accelerated aging in normal breast tissue of women with breast cancer

Background

DNA methylation alterations have similar patterns in normal aging tissue and in cancer. In this study, we investigated breast tissue-specific age-related DNA methylation alterations and used those methylation sites to identify individuals with outlier phenotypes. Outlier phenotype is identified by unsupervised anomaly detection algorithms and is defined by individuals who have normal tissue age-dependent DNA methylation levels that vary dramatically from the population mean.

Methods

We generated whole-genome DNA methylation profiles (GSE160233) on purified epithelial cells and used publicly available Infinium HumanMethylation 450K array datasets (TCGA, GSE88883, GSE69914, GSE101961, and GSE74214) for discovery and validation.

Results

We found that hypermethylation in normal breast tissue is the best predictor of hypermethylation in cancer. Using unsupervised anomaly detection approaches, we found that about 10% of the individuals (39/427) were outliers for DNA methylation from 6 DNA methylation datasets. We also found that there were significantly more outlier samples in normal-adjacent to cancer (24/139, 17.3%) than in normal samples (15/228, 5.2%). Additionally, we found significant differences between the predicted ages based on DNA methylation and the chronological ages among outliers and not-outliers. Additionally, we found that accelerated outliers (older predicted age) were more frequent in normal-adjacent to cancer (14/17, 82%) compared to normal samples from individuals without cancer (3/17, 18%). Furthermore, in matched samples, we found that the epigenome of the outliers in the pre-malignant tissue was as severely altered as in cancer.

Conclusions

A subset of patients with breast cancer has severely altered epigenomes which are characterized by accelerated aging in their normal-appearing tissue. In the future, these DNA methylation sites should be studied further such as in cell-free DNA to determine their potential use as biomarkers for early detection of malignant transformation and preventive intervention in breast cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01434-7.

The Hypertension Related Gene G-Protein Coupled Receptor Kinase 4 Contributes to Breast Cancer Proliferation

PURPOSE

Clinical studies have shown that breast cancer risk is increased in hypertensive women. The underlying molecular mechanism remains undetermined. The current study tests our hypothesis that G protein coupled receptor kinase 4 (GRK4) is a molecule that links hypertension and breast cancer. GRK4 plays an important role in regulation of renal sodium excretion. Sustained activation of GRK4 as in the circumstances of single nucleotide polymorphism (SNPs) causes hypertension. Expression of GRK4 in the kidney is regulated by cMyc, an oncogene that is amplified in breast cancer.

METHODS

Western analysis was used to evaluate GRK4 protein expression in seven breast cancer cell lines. GRK4 gene single nucleotide polymorphism in breast cancer cell lines and in breast cancer cDNA arrays were determined using TaqMan Genotyping qPRC. The function of GRK4 was evaluated in MCF-7 cells with cMyc knock-down and GRK4 re-expression and in MDA-MB-468 cells expressing inducible GRK4 shRNA lentivirus constructs. Nuclei counting and 5-Bromo-2'-deoxy-uridine (BrdU) labeling were used to evaluate cell growth and proliferation.

RESULTS

Genotyping of GRK4 SNPs in breast cancer cDNA arrays (n = 94) revealed that the frequency of carrying two hypertension related SNPs A142 V or R65 L is threefold higher in breast cancer patients than in healthy people (P = 7.53E-11). GRK4 protein is expressed in seven breast cancer cell lines but not the benign mammary epithelial cell line, MCF-10A. Three hypertension related SNPs in the GRK4 gene were identified in the breast cancer cell lines. Except for BT20, all other breast cancer lines have 1-3 GRK4 SNPs of which A142 V occurs in all 6 lines. MDA-MB-468 cells contain homozygous A142 V and R65 L SNPs. Knocking down cMyc in MCF-7 cells significantly reduced the growth rate, which was rescued by re-expression of GRK4. We then generated three stable GRK4 knock-down MDA-MB-468 lines using inducible lentiviral shRNA vectors. Doxycycline (Dox) induced GRK4 silencing significantly reduced GRK4 mRNA and protein levels, growth rates, and proliferation. As a marker of cell proliferation, the percentage of BrdU-labeled cells decreased from 45 ± 3% in the cells without Dox to 32 ± 5% in the cells treated with 0.1 µg/mL Dox.

CONCLUSIONS

GRK4 acts as an independent proliferation promotor in breast cancer. Our results suggest that targeted inhibition of GRK4 could be a new therapy for both hypertension and breast cancer.

NLRP1 variant M1184V decreases inflammasome activation in the context of DPP9 inhibition and asthma severity

Background

NLRP1 is an innate immune sensor that can form cytoplasmic inflammasome complexes. Polymorphisms in NLRP1 are linked to asthma; however, there is currently no functional or mechanistic explanation for this.

Objective

We sought to clarify the role of NLRP1 in asthma pathogenesis.

Methods

Results from the GALA II cohort study were used to identify a link between NLRP1 and asthma in Mexican Americans. In vitro and in vivo models for NLRP1 activation were applied to investigate the role of this inflammasome in asthma at the molecular level.

Results

We document the association of an NLRP1 haplotype with asthma for which the single nucleotide polymorphism rs11651270 (M1184V) individually is the most significant. Surprisingly, M1184V increases NLRP1 activation in the context of N-terminal destabilization, but decreases NLRP1 activation on dipeptidyl peptidase 9 inhibition. In vitro studies demonstrate that M1184V increases binding to dipeptidyl peptidase 9, which can account for its inhibitory role in this context. In addition, in vivo data from a mouse model of airway inflammation reveal a protective role for NLRP1 inflammasome activation reducing eosinophilia in this setting.

Conclusions

Linking our in vitro and in vivo results, we found that the NLRP1 variant M1184V reduces inflammasome activation in the context of dipeptidyl peptidase 9 inhibition and could thereby increase asthma severity. Our studies may have implications for the treatment of asthma in patients carrying this variant of NLRP1.

Tumor suppressor genes and allele-specific expression: mechanisms and significance

Recent findings indicate that allele-specific expression (ASE) at specific cancer driver gene loci may be of importance in onset/progression of the disease. Of particular interest are loss-of-function (LOF) of tumor suppressor gene (TSGs) alleles. While LOF tumor suppressor mutations are typically considered to be recessive, if these mutant alleles can be significantly differentially expressed relative to wild-type alleles in heterozygotes, the clinical consequences could be significant. LOF TSG alleles are shown to be segregating at high frequencies in world-wide populations of normal/healthy individuals. Matched sets of normal and tumor tissues isolated from 233 cancer patients representing four diverse tumor types demonstrate functionally important changes in patterns of ASE in individuals heterozygous for LOF TSG alleles associated with cancer onset/progression. While a variety of molecular mechanisms were identified as potentially contributing to changes in ASE patterns in cancer, changes in DNA copy number and allele-specific alternative splicing possibly mediated by antisense RNA emerged as predominant factors. In conclusion, LOF TSGs are segregating in human populations at significant frequencies indicating that many otherwise healthy individuals are at elevated risk of developing cancer. Changes in ASE between normal and cancer tissues indicates that LOF TSG alleles may contribute to cancer onset/progression even when heterozygous with wild-type functional alleles.

A Case/Control Study: AGBL1 Polymorphism Related to Lung Cancer Risk in Chinese Nonsmoking Females

ATP/GTP binding protein like 1 (AGBL1) plays a role in controlling the length of polyglutamate side chains. Polymorphism rs4513061 in AGBL1 is suspected to influence the risk of lung cancer. A case/control study was performed involving 556 cases and 563 controls from a hospital participating in donation. The relationship between rs4513061 and the risk of lung cancer and the interaction between rs4513061 and environmental exposure were determined by the chi-square tests, logistic regression analysis, and crossover analysis. The survival analysis was conducted by Cox proportional hazard regression. The results showed that rs4513061 polymorphism is associated with the risk of lung cancer. The stratified analysis suggested the protective effect of rs4513061 to different histological types of lung cancer, including lung adenocarcinoma (AA vs. GG: odds ratio [OR] = 0.505, 95% confidence interval [CI] = 0.337-0.756, p < 0.001), squamous cell lung cancer (AG vs. GG: OR = 0.488, 95% CI = 0.269-0.883, p = 0.018), and small-cell lung cancer (AA vs. GG: OR = 0.421, 95% CI = 0.216-0.819, p = 0.011). Nevertheless, there was no significant interaction between rs4513061 and cooking oil fume. Significant impact was not observed between the rs4513061 polymorphism and survival time of lung cancer. Our study indicated that rs4513061 in AGBL1 decreases the risk of lung cancer in nonsmoking females from northeast China.

Explorative analysis of the gene expression profile during liver regeneration of mouse: a microarray-based study

The liver is an amazing organ due to its powerful regenerative capacity. Although many studies on liver regeneration have been documented, the detailed mechanisms remain unclear. Two-third partial hepatectomy (PH) in rodents plays a crucial role in the study of liver regeneration. In this study, the time series data of gene expression during liver regeneration in mouse were analyzed using the gene set numbered GSE6998 in GEO. A variety of bioinformatics methods, including masigPro, Weighted Gene Co-expression Network Analysis (WGCNA), spatial analysis of functional enrichment (SAFE) and ingenuity canonical pathway analysis (IPA) were used to identify and compare the significantly changed pathways, potential upstream regulators and key genes during liver regeneration. Our study showed that liver regeneration in the mouse is a coordinated process, which cell-cycle-related progress are at the centre of the interaction network involved in liver regeneration. Several candidate upstream regulators including PPARA, NFE2L2, MAD1 and CNR1 and some key genes such as Cdk1, Plk1, Cdc20, Aurka, Racgap1, Cenpa, Rrm1, Rrm2 were identified. In conclusion, these findings could contribute to revealing the molecular mechanism of liver regeneration after PH, which could provide new ideas and treatment methods for regenerative medicine, oncological drug development and oncological treatment.

The Protein Encoded by the CCDC170 Breast Cancer Gene Functions to Organize the Golgi-Microtubule Network

Genome-Wide Association Studies (GWAS) and subsequent fine-mapping studies (>50) have implicated single nucleotide polymorphisms (SNPs) located at the CCDC170/C6ORF97-ESR1 locus (6q25.1) as being associated with the risk of breast cancer. Surprisingly, our analysis using genome-wide differential allele-specific expression (DASE), an indicator for breast cancer susceptibility, suggested that the genetic alterations of CCDC170, but not ESR1, account for GWAS-associated breast cancer risk at this locus. Breast cancer-associated CCDC170 nonsense mutations and rearrangements have also been detected, with the latter being specifically implicated in driving breast cancer. Here we report that the wild type CCDC170 protein localizes to the region of the Golgi apparatus and binds Golgi-associated microtubules (MTs), and that breast cancer-linked truncations of CCDC170 result in loss of Golgi localization. Overexpression of wild type CCDC170 triggers Golgi reorganization, and enhances Golgi-associated MT stabilization and acetyltransferase ATAT1-dependent α-tubulin acetylation. Golgi-derived MTs regulate cellular polarity and motility, and we provide evidence that dysregulation of CCDC170 affects polarized cell migration. Taken together, our findings demonstrate that CCDC170 plays an essential role in Golgi-associated MT organization and stabilization, and implicate a mechanism for how perturbations in the CCDC170 gene may contribute to the hallmark changes in cell polarity and motility seen in breast cancer.

LincIN, a novel NF90-binding long non-coding RNA, is overexpressed in advanced breast tumors and involved in metastasis

BACKGROUND

Recent genome-wide profiling by sequencing and distinctive chromatin signatures has identified thousands of long non-coding RNA (lncRNA) species (>200 nt). LncRNAs have emerged as important regulators of gene expression, involving in both developmental and pathological processes. While altered expression of lncRNAs has been observed in breast cancer development, their roles in breast cancer progression and metastasis are still poorly understood.

METHODS

To identify novel breast cancer-associated lncRNA candidates, we employed a high-density SNP array-based approach to uncover intergenic lncRNA genes that are aberrantly expressed in breast cancer. We first evaluated the potential value as a breast cancer prognostic biomarker for one breast cancer-associated lncRNA, LincIN, using a breast cancer cohort retrieved from The Cancer Genome Atlas (TCGA) Data Portal. Then we characterized the role of LincIN in breast cancer progression and metastasis by in vitro invasion assay and a mouse tail vein injection metastasis model. To study the action of LincIN, we identified LincIN-interacting protein partner(s) by RNA pull-down experiments followed with protein identification by mass spectrometry.

RESULTS

High levels of LincIN expression are frequently observed in tumors compared to adjacent normal tissues, and are strongly associated with aggressive breast cancer. Importantly, analysis of TCGA data further suggest that high expression of LincIN is associated with poor overall survival in patients with breast cancer (P = 0.044 and P = 0.011 after adjustment for age). The functional experiments demonstrate that knockdown of LincIN inhibits tumor cell migration and invasion in vitro, which is supported by the results of transcriptome analysis in the LincIN-knockdown cells. Furthermore, knockdown of LincIN diminishes lung metastasis in a mouse tail vein injection model. We also identified a LincIN-binding protein, NF90, through which overexpression of LincIN may repress p21 protein expression by inhibiting its translation, and upregulation of p21 by LincIN knockdown may be associated with less aggressive metastasis phenotypes.

CONCLUSIONS

Our studies provide clear evidence to support LincIN as a new regulator of tumor progression-metastasis at both transcriptional and translational levels and as a promising prognostic biomarker for breast cancer.

Transcriptional Selectivity of Epigenetic Therapy in Cancer

A central challenge in the development of epigenetic cancer therapy is the ability to direct selectivity in modulating gene expression for disease-selective efficacy. To address this issue, we characterized by RNA-seq, DNA methylation, and ChIP-seq analyses the epigenetic response of a set of colon, breast, and leukemia cancer cell lines to small-molecule inhibitors against DNA methyltransferases (DAC), histone deacetylases (Depsi), histone demethylases (KDM1A inhibitor S2101), and histone methylases (EHMT2 inhibitor UNC0638 and EZH2 inhibitor GSK343). We also characterized the effects of DAC as combined with the other compounds. Averaged over the cancer cell models used, we found that DAC affected 8.6% of the transcriptome and that 95.4% of the genes affected were upregulated. DAC preferentially regulated genes that were silenced in cancer and that were methylated at their promoters. In contrast, Depsi affected the expression of 30.4% of the transcriptome but showed little selectivity for gene upregulation or silenced genes. S2101, UNC0638, and GSK343 affected only 2% of the transcriptome, with UNC0638 and GSK343 preferentially targeting genes marked with H3K9me2 or H3K27me3, respectively. When combined with histone methylase inhibitors, the extent of gene upregulation by DAC was extended while still maintaining selectivity for DNA-methylated genes and silenced genes. However, the genes upregulated by combination treatment exhibited limited overlap, indicating the possibility of targeting distinct sets of genes based on different epigenetic therapy combinations. Overall, our results demonstrated that DNA methyltransferase inhibitors preferentially target cancer-relevant genes and can be combined with inhibitors targeting histone methylation for synergistic effects while still maintaining selectivity. Cancer Res; 77(2); 470-81. ©2016 AACR.

Extensive variation between tissues in allele specific expression in an outbred mammal

BACKGROUND

Allele specific gene expression (ASE), with the paternal allele more expressed than the maternal allele or vice versa, appears to be a common phenomenon in humans and mice. In other species the extent of ASE is unknown, and even in humans and mice there are several outstanding questions. These include; to what extent is ASE tissue specific? how often does the direction of allele expression imbalance reverse between tissues? how often is only one of the two alleles expressed? is there a genome wide bias towards expression of the paternal or maternal allele; and finally do genes that are nearby on a chromosome share the same direction of ASE? Here we use gene expression data (RNASeq) from 18 tissues from a single cow to investigate each of these questions in turn, and then validate some of these findings in two tissues from 20 cows.

RESULTS

Between 40 and 100 million sequence reads were generated per tissue across three replicate samples for each of the eighteen tissues from the single cow (the discovery dataset). A bovine gene expression atlas was created (the first from RNASeq data), and differentially expressed genes in each tissue were identified. To analyse ASE, we had access to unambiguously phased genotypes for all heterozygous variants in the cow's whole genome sequence, where these variants were homozygous in the whole genome sequence of her sire, and as a result we were able to map reads to parental genomes, to determine SNP and genes showing ASE in each tissue. In total 25,251 heterozygous SNP within 7985 genes were tested for ASE in at least one tissue. ASE was pervasive, 89 % of genes tested had significant ASE in at least one tissue. This large proportion of genes displaying ASE was confirmed in the two tissues in a validation dataset. For individual tissues the proportion of genes showing significant ASE varied from as low as 8-16 % of those tested in thymus to as high as 71-82 % of those tested in lung. There were a number of cases where the direction of allele expression imbalance reversed between tissues. For example the gene SPTY2D1 showed almost complete paternal allele expression in kidney and thymus, and almost complete maternal allele expression in the brain caudal lobe and brain cerebellum. Mono allelic expression (MAE) was common, with 1349 of 4856 genes (28 %) tested with more than one heterozygous SNP showing MAE. Across all tissues, 54.17 % of all genes with ASE favoured the paternal allele. Genes that are closely linked on the chromosome were more likely to show higher expression of the same allele (paternal or maternal) than expected by chance. We identified several long runs of neighbouring genes that showed either paternal or maternal ASE, one example was five adjacent genes (GIMAP8, GIMAP7 copy1, GIMAP4, GIMAP7 copy 2 and GIMAP5) that showed almost exclusive paternal expression in brain caudal lobe.

CONCLUSIONS

Investigating the extent of ASE across 18 bovine tissues in one cow and two tissues in 20 cows demonstrated 1) ASE is pervasive in cattle, 2) the ASE is often MAE but ranges from MAE to slight overexpression of the major allele, 3) the ASE is most often tissue specific and that more than half the time displays divergent allele specific expression patterns across tissues, 4) across all genes there is a slight bias towards expression of the paternal allele and 5) genes expressing the same parental allele are clustered together more than expected by chance, and there are several runs of large numbers of genes expressing the same parental allele.

Alterations in Salivary Proteome following Single Twenty-Minute Session of Yogic Breathing

Yogic breathing (YB) has been suggested to reduce stress and blood pressure and increase cognitive processes. However, alterations after YB at the molecular level are not well established. Twenty healthy volunteers were randomized into two groups (N = 10 per group): YB or attention controls (AC). The YB group performed two YB exercises, each for ten minutes, for a total of twenty minutes in a single session. AC group read a text of their choice for 20 minutes. Saliva was collected at baseline and at 5, 10, 15, and 20 minutes. Using Mass Spectrometry (MS), we initially found that 22 proteins were differentially expressed and then validated deleted in malignant brain tumor-1 (DMBT1) and Ig lambda-2 chain C region (IGLC2) using Western Blotting. DMBT1 was elevated in 7 of YB group by 10-fold and 11-fold at 10 and 15 minutes, respectively, whereas it was undetectable in the time-matched AC group (P < 0.05). There was a significant interaction between groups and time assessed by two-way ANOVA (P < 0.001). IGLC2 also showed a significant increase in YB group as measured by Western Blotting. These data are the first to demonstrate the feasibility of stimulating and detecting salivary protein biomarkers in response to an acute Yoga exercise. This trial is registered with ClincalTrials.gov NCT02108769.

Tumor haplotype assembly algorithms for cancer genomics

The growing availability of inexpensive high-throughput sequence data is enabling researchers to sequence tumor populations within a single individual at high coverage. But, cancer genome sequence evolution and mutational phenomena like driver mutations and gene fusions are difficult to investigate without first reconstructing tumor haplotype sequences. Haplotype assembly of single individual tumor populations is an exceedingly difficult task complicated by tumor haplotype heterogeneity, tumor or normal cell sequence contamination, polyploidy, and complex patterns of variation. While computational and experimental haplotype phasing of diploid genomes has seen much progress in recent years, haplotype assembly in cancer genomes remains uncharted territory. In this work, we describe HapCompass-Tumor a computational modeling and algorithmic framework for haplotype assembly of copy number variable cancer genomes containing haplotypes at different frequencies and complex variation. We extend our polyploid haplotype assembly model and present novel algorithms for (1) complex variations, including copy number changes, as varying numbers of disjoint paths in an associated graph, (2) variable haplotype frequencies and contamination, and (3) computation of tumor haplotypes using simple cycles of the compass graph which constrain the space of haplotype assembly solutions. The model and algorithm are implemented in the software package HapCompass-Tumor which is available for download from http://www.brown.edu/Research/Istrail_Lab/.

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.