Human basonuclin 2 up-regulates a cascade set of interferon-stimulated genes with anti-cancerous properties in a lung cancer model

BACKGROUND

Human basonuclin 2 (BNC2) acts as a tumor suppressor in multiple cancers in an as yet unidentified manner. The role and expression of the BNC2 gene in lung cancer has not yet been investigated.

METHODS

BNC2 expression was studied in the A549 and BEAS-2B cell lines, as well as in lung cancer tissue. Illumina array analysis and a viability assay were used to study the effects of transient transfection of BNC2 in A549 cells. Ingenuity pathway analysis and g:Profiler were applied to identify affected pathways and networks. RT-qPCR was used to validate the array results.

RESULTS

We showed the reduced mRNA expression of BNC2 in non-small cell lung cancer tissue and lung cancer cell line A549 compared to non-cancerous lung tissue and BEAS-2B cells, respectively. Further array analysis demonstrated that the transfection of BNC2 into A549 cells resulted in the increased expression of 139 genes and the down-regulation of 13 genes. Pathway analysis revealed that half of the up-regulated genes were from the interferon/signal transducer and activator of transcription signaling pathways. The differential expression of selected sets of genes, including interferon-stimulated and tumor suppressor genes of the XAF1 and OAS families, was confirmed by RT-qPCR. In addition, we showed that the over-expression of BNC2 inhibited the proliferation of A549 cells.

CONCLUSION

Our data suggest that human BNC2 is an activator of a subset of IFN-regulated genes and might thereby act as a tumor suppressor.

The transcriptional landscape of age in human peripheral blood

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the ‘transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Ageing increases the risk of many diseases. Here the authors compare blood cell transcriptomes of over 14,000 individuals and identify a set of about 1,500 genes that are differently expressed with age, shedding light on transcriptional programs linked to the ageing process and age-associated diseases.

Extensive alterations of the whole-blood transcriptome are associated with body mass index: results of an mRNA profiling study involving two large population-based cohorts

Background

Obesity, defined as pathologically increased body mass index (BMI), is strongly related to an increased risk for numerous common cardiovascular and metabolic diseases. It is particularly associated with insulin resistance, hyperglycemia, and systemic oxidative stress and represents the most important risk factor for type 2 diabetes (T2D). However, the pathophysiological mechanisms underlying these associations are still not completely understood. Therefore, in order to identify potentially disease-relevant BMI-associated gene expression signatures, a transcriptome-wide association study (TWAS) on BMI was performed.

Methods

Whole-blood mRNA levels determined by array-based transcriptional profiling were correlated with BMI in two large independent population-based cohort studies (KORA F4 and SHIP-TREND) comprising a total of 1977 individuals.

Results

Extensive alterations of the whole-blood transcriptome were associated with BMI: More than 3500 transcripts exhibited significant positive or negative BMI-correlation. Three major whole-blood gene expression signatures associated with increased BMI were identified. The three signatures suggested: i) a ratio shift from mature erythrocytes towards reticulocytes, ii) decreased expression of several genes essentially involved in the transmission and amplification of the insulin signal, and iii) reduced expression of several key genes involved in the defence against reactive oxygen species (ROS).

Conclusions

Whereas the first signature confirms published results, the other two provide possible mechanistic explanations for well-known epidemiological findings under conditions of increased BMI, namely attenuated insulin signaling and increased oxidative stress. The putatively causative BMI-dependent down-regulation of the expression of numerous genes on the mRNA level represents a novel finding. BMI-associated negative transcriptional regulation of insulin signaling and oxidative stress management provide new insights into the pathogenesis of metabolic syndrome and T2D.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-015-0141-x) contains supplementary material, which is available to authorized users.

A meta-analysis of gene expression signatures of blood pressure and hypertension

Genome-wide association studies (GWAS) have uncovered numerous genetic variants (SNPs) that are associated with blood pressure (BP). Genetic variants may lead to BP changes by acting on intermediate molecular phenotypes such as coded protein sequence or gene expression, which in turn affect BP variability. Therefore, characterizing genes whose expression is associated with BP may reveal cellular processes involved in BP regulation and uncover how transcripts mediate genetic and environmental effects on BP variability. A meta-analysis of results from six studies of global gene expression profiles of BP and hypertension in whole blood was performed in 7017 individuals who were not receiving antihypertensive drug treatment. We identified 34 genes that were differentially expressed in relation to BP (Bonferroni-corrected p<0.05). Among these genes, FOS and PTGS2 have been previously reported to be involved in BP-related processes; the others are novel. The top BP signature genes in aggregate explain 5%–9% of inter-individual variance in BP. Of note, rs3184504 in SH2B3, which was also reported in GWAS to be associated with BP, was found to be a trans regulator of the expression of 6 of the transcripts we found to be associated with BP (FOS, MYADM, PP1R15A, TAGAP, S100A10, and FGBP2). Gene set enrichment analysis suggested that the BP-related global gene expression changes include genes involved in inflammatory response and apoptosis pathways. Our study provides new insights into molecular mechanisms underlying BP regulation, and suggests novel transcriptomic markers for the treatment and prevention of hypertension.

The focus of blood pressure (BP) GWAS has been the identification of common DNA sequence variants associated with the phenotype; this approach provides only one dimension of molecular information about BP. While it is a critical dimension, analyzing DNA variation alone is not sufficient for achieving an understanding of the multidimensional complexity of BP physiology. The top loci identified by GWAS explain only about 1 percent of inter-individual BP variability. In this study, we performed a meta-analysis of gene expression profiles in relation to BP and hypertension in 7017 individuals from six studies. We identified 34 differentially expressed genes for BP, and discovered that the top BP signature genes explain 5%–9% of BP variability. We further linked BP gene expression signature genes with BP GWAS results by integrating expression associated SNPs (eSNPs) and discovered that one of the top BP loci from GWAS, rs3184504 in SH2B3, is a trans regulator of expression of 6 of the top 34 BP signature genes. Our study, in conjunction with prior GWAS, provides a deeper understanding of the molecular and genetic basis of BP regulation, and identifies several potential targets and pathways for the treatment and prevention of hypertension and its sequelae.

New genetic loci link adipose and insulin biology to body fat distribution

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, we conducted genome-wide association meta-analyses of waist and hip circumference-related traits in up to 224,459 individuals. We identified 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (WHRadjBMI) and an additional 19 loci newly associated with related waist and hip circumference measures (P<5×10⁻⁸). Twenty of the 49 WHRadjBMI loci showed significant sexual dimorphism, 19 of which displayed a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation, and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

Mapping the genetic architecture of gene regulation in whole blood

BACKGROUND

We aimed to assess whether whole blood expression quantitative trait loci (eQTLs) with effects in cis and trans are robust and can be used to identify regulatory pathways affecting disease susceptibility.

MATERIALS AND METHODS

We performed whole-genome eQTL analyses in 890 participants of the KORA F4 study and in two independent replication samples (SHIP-TREND, N = 976 and EGCUT, N = 842) using linear regression models and Bonferroni correction.

RESULTS

In the KORA F4 study, 4,116 cis-eQTLs (defined as SNP-probe pairs where the SNP is located within a 500 kb window around the transcription unit) and 94 trans-eQTLs reached genome-wide significance and overall 91% (92% of cis-, 84% of trans-eQTLs) were confirmed in at least one of the two replication studies. Different study designs including distinct laboratory reagents (PAXgene™ vs. Tempus™ tubes) did not affect reproducibility (separate overall replication overlap: 78% and 82%). Immune response pathways were enriched in cis- and trans-eQTLs and significant cis-eQTLs were partly coexistent in other tissues (cross-tissue similarity 40-70%). Furthermore, four chromosomal regions displayed simultaneous impact on multiple gene expression levels in trans, and 746 eQTL-SNPs have been previously reported to have clinical relevance. We demonstrated cross-associations between eQTL-SNPs, gene expression levels in trans, and clinical phenotypes as well as a link between eQTLs and human metabolic traits via modification of gene regulation in cis.

CONCLUSIONS

Our data suggest that whole blood is a robust tissue for eQTL analysis and may be used both for biomarker studies and to enhance our understanding of molecular mechanisms underlying gene-disease associations.

A common 16p11.2 inversion underlies the joint susceptibility to asthma and obesity

The prevalence of asthma and obesity is increasing worldwide, and obesity is a well-documented risk factor for asthma. The mechanisms underlying this association and parallel time trends remain largely unknown but genetic factors may be involved. Here, we report on a common ~0.45 Mb genomic inversion at 16p11.2 that can be accurately genotyped via SNP array data. We show that the inversion allele protects against the joint occurrence of asthma and obesity in five large independent studies (combined sample size of 317 cases and 543 controls drawn from a total of 5,809 samples; combined OR = 0.48, p = 5.5 × 10(-6)). Allele frequencies show remarkable worldwide population stratification, ranging from 10% in East Africa to 49% in Northern Europe, consistent with discordant and extreme genetic drifts or adaptive selections after human migration out of Africa. Inversion alleles strongly correlate with expression levels of neighboring genes, especially TUFM (p = 3.0 × 10(-40)) that encodes a mitochondrial protein regulator of energy balance and inhibitor of type 1 interferon, and other candidates for asthma (IL27) and obesity (APOB48R and SH2B1). Therefore, by affecting gene expression, the ~0.45 Mb 16p11.2 inversion provides a genetic basis for the joint susceptibility to asthma and obesity, with a population attributable risk of 39.7%. Differential mitochondrial function and basal energy balance of inversion alleles might also underlie the potential selection signature that led to their uneven distribution in world populations.

Gene expression alterations in the cerebellum and granule neurons of Cstb(-/-) mouse are associated with early synaptic changes and inflammation

Progressive myoclonus epilepsy of Unverricht-Lundborg type (EPM1) is an autosomal recessively inherited neurodegenerative disease, manifesting with myoclonus, seizures and ataxia, caused by mutations in the cystatin B (CSTB) gene. With the aim of understanding the molecular basis of pathogenetic events in EPM1 we characterized gene expression changes in the cerebella of pre-symptomatic postnatal day 7 (P7) and symptomatic P30 cystatin B -deficient (Cstb(-/-) ) mice, a model for the disease, and in cultured Cstb(-/-) cerebellar granule cells using a pathway-based approach. Differentially expressed genes in P7 cerebella were connected to synaptic function and plasticity, and in cultured cerebellar granule cells, to cell cycle, cytoskeleton, and intracellular transport. In particular, the gene expression data pinpointed alterations in GABAergic pathway. Electrophysiological recordings from Cstb(-/-) cerebellar Purkinje cells revealed a shift of the balance towards decreased inhibition, yet the amount of inhibitory interneurons was not declined in young animals. Instead, we found diminished number of GABAergic terminals and reduced ligand binding to GABAA receptors in Cstb(-/-) cerebellum. These results suggest that alterations in GABAergic signaling could result in reduced inhibition in Cstb(-/-) cerebellum leading to the hyperexcitable phenotype of Cstb(-/-) mice. At P30, the microarray data revealed a marked upregulation of immune and defense response genes, compatible with the previously reported early glial activation that precedes neuronal degeneration. This further implies the role of early-onset neuroinflammation in the pathogenesis of EPM1.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Human disease-associated genetic variation impacts large intergenic non-coding RNA expression

Recently it has become clear that only a small percentage (7%) of disease-associated single nucleotide polymorphisms (SNPs) are located in protein-coding regions, while the remaining 93% are located in gene regulatory regions or in intergenic regions. Thus, the understanding of how genetic variations control the expression of non-coding RNAs (in a tissue-dependent manner) has far-reaching implications. We tested the association of SNPs with expression levels (eQTLs) of large intergenic non-coding RNAs (lincRNAs), using genome-wide gene expression and genotype data from five different tissues. We identified 112 cis-regulated lincRNAs, of which 45% could be replicated in an independent dataset. We observed that 75% of the SNPs affecting lincRNA expression (lincRNA cis-eQTLs) were specific to lincRNA alone and did not affect the expression of neighboring protein-coding genes. We show that this specific genotype-lincRNA expression correlation is tissue-dependent and that many of these lincRNA cis-eQTL SNPs are also associated with complex traits and diseases.

Large intergenic non-coding RNAs (lincRNAs) are the largest class of non-coding RNA molecules in the human genome. Many genome-wide association studies (GWAS) have mapped disease-associated genetic variants (SNPs) to, or in, the vicinity of such lincRNA regions. However, it is not clear how these SNPs can affect the disease. We tested whether SNPs were also associated with the lincRNA expression levels in five different human primary tissues. We observed that there is a strong genotype-lincRNA expression correlation that is tissue-dependent. Many of the observed lincRNA cis-eQTLs are disease- or trait-associated SNPs. Our results suggest that lincRNA-eQTLs represent a novel link between non-coding SNPs and the expression of protein-coding genes, which can be exploited to understand the process of gene-regulation through lincRNAs in more detail.

Polymorphisms in ESR1, ESR2 and HSD17B1 genes are associated with fertility status in endometriosis

OBJECTIVE

To investigate whether polymorphisms in genes involved in biosynthesis and signalling of sex steroids influence susceptibility to endometriosis and to infertility associated with it.

MATERIALS AND METHODS

Patients with endometriosis (n = 150) and fertile controls (n = 199) were genotyped for polymorphisms in oestrogen receptor genes ESR1 (rs2234693 - T/C single nucleotide polymorphism (SNP), dinucleotide (TA)(n) repeat) and ESR2 (dinucleotide (CA)(n) repeat), progesterone receptor gene PGR (rs10895068 - G/A SNP, 306-bp Alu-insertion), 17β-hydroxysteroid dehydrogenase type 1 gene HSD17B1 (rs605059 - A/G SNP), and aromatase gene CYP19A1 (rs10046 - C/T SNP, (TTTA)(n) tetranucleotide repeat, 3-bp TCT insertion/deletion polymorphism).

RESULTS

The HSD17B1 A/G SNP A allele increased overall endometriosis risk and the risk of stage I-II disease, while ESR1 longer (TA)(n) repeats only correlated with susceptibility to stage I-II endometriosis. When considering patients' fertility status, HSD17B1 A/G SNP A allele and ESR1 longer (TA)(n) repeats were associated with endometriosis accompanied by infertility, while ESR2 shorter (CA)(n) repeats were linked with endometriosis without infertility. Other polymorphisms were distributed similarly among patients and controls.

CONCLUSIONS

Genetic variants in ESR1, ESR2, and HSD17B1 genes could modify susceptibility to endometriosis and might influence the fertility status in endometriosis patients.

166 THE DIFFERENTIAL TRANSCRIPTOME AND ONTOLOGY PROFILES OF MURAL AND CUMULUS GRANULOSA CELLS IN STIMULATED HUMAN ANTRAL FOLLICLES

Communication between various cell types in the ovary is a prerequisite for successful folliculogenesis and ovulation. In human antral follicles, granulosa cells divide into 2 distinct cell populations: mural (enveloping the antrum, MGC) and cumulus granulosa cells (surrounding the oocyte, CGC). During infertility treatment using in vitro fertilization (IVF), granulosa cells can be retrieved during puncture of stimulated follicles offering an excellent opportunity for analysing their functional properties. The aim of this study was to compare the transcriptomes of MGC and CGC. Twenty infertile women undergoing IVF-ICSI treatment were enrolled. The MGC were obtained from follicular fluid and CGC were acquired after oocyte denudation before micromanipulation. Gene expression of both cell populations was analysed using a genome-wide transcription array. The expression profile of the 2 granulosa cell populations varied significantly: out of 28 869 transcripts, 4480 were differentially expressed (q-value < 10–4); 623 transcripts differed in their expression levels by at least 2-fold. The transcriptome of MGC showed higher expression of genes involved in immune regulation (toll-like receptors, IL18, IL17R). In CGC, pathways participating in intercellular interactions, tissue remodelling and protein degradation were more clearly distinguished (tenascin C, IGFBP5). Among the identified differentially expressed genes, several are involved in follicle development, oocyte maturation, or ovulatory processes. Our findings fit well with previously published data. The results provide a basis for future studies on intra- and intercellular signaling in the preovulatory follicle leading towards identifying methods for improving oocyte health, embryo selection, and ultimately IVF success rate. This work was supported by the Estonian Science Foundation (6498, 6585, and 7479); the Estonian Ministry of Education and Science (SF0182641s04, SF0180142Cs08, SF0180044s09, and PBGMR07903); the Estonian University of Life Sciences (P8001VLVL); the European Union Sixth Framework Programme (FP6) (LSHB-CT-2004-503243); the European Cooperation 22 in Science and Technology (COST-STSM-FA0702-03734), and the European Union through the European Regional Development Fund through the Centre of Excellence in Genomics, Estonian Biocentre, and University of Tartu. The authors acknowledge our colleagues from the Nova Vita Clinic for their assistance in sample collection and all patient volunteers for their participation in the study.

The differential transcriptome and ontology profiles of floating and cumulus granulosa cells in stimulated human antral follicles

Communication between various ovarian cell types is a prerequisite for folliculogenesis and ovulation. In antral follicles granulosa cells divide into two distinct populations of mural and cumulus granulosa cells (CGC), enveloping the antrum and surrounding the oocyte, respectively. Both cell types, with the mural compartment in excess, contribute to the floating granulosa cell (FGC) population in the follicular fluid. The aim of this study was to compare the transcriptomes of FGC and CGC in stimulated antral follicles obtained from 19 women undergoing IVF-ICSI procedure. FGC were obtained from follicular fluid during the follicle puncture procedure and CGC were acquired after oocyte denudation for micromanipulation. Gene expression analysis was conducted using the genome-wide Affymetrix transcriptome array. The expression profile of the two granulosa cell populations varied significantly. Out of 28 869 analysed transcripts 4480 were differentially expressed (q-value < 10(-4)) and 489 showed > or =2-fold difference in the expression level with 222 genes up-regulated in FGC and 267 in CGC. The transcriptome of FGC showed higher expression of genes involved in immune response, hematological system function and organismal injury, although CGC had genes involved in protein degradation and nervous system function up-regulated. Cell-to-cell signalling and interaction pathways were noted in both cell populations. Furthermore, numerous novel transcripts that have not been previously described in follicular physiology were identified. In conclusion, our results provide a solid basis for future studies in follicular biology that will help to identify molecular markers for oocyte and embryo viability in IVF.