Analysis of intron sequence features associated with transcriptional regulation in human genes

Genome-wide sequencing identified extrachromosomal circular DNA as a transcription factor-binding motif of the senescence genes that govern replicative senescence in human mesenchymal stem cells

Introduction

Mesenchymal stem cells (MSCs) have long been postulated as an important source cell in regenerative medicine. During subculture expansion, mesenchymal stem cell (MSC) senescence diminishes their multi-differentiation capabilities, leading to a loss of therapeutic potential. Up to date, the extrachromosomal circular DNAs (eccDNAs) have been demonstrated to be involved in senescence but the roles of eccDNAs during MSC.

Methods

Here we explored eccDNA profiles in human bone marrow MSCs (BM-MSCs). EccDNA and mRNA was purified and sequenced, followed by quantification and functional annotation. Moreover, we mapped our datasets with the downloading enhancer and transcription factor-regulated genes to explore the potential role of eccDNAs.

Results

Sequentially, gene annotation analysis revealed that the majority of eccDNA were mapped in the intron regions with limited BM-MSC enhancer overlaps. We discovered that these eccDNA motifs in senescent BMSCs acted as motifs for binding transcription factors (TFs) of senescence-related genes.

Discussion

These findings are highly significant for identifying biomarkers of senescence and therapeutic targets in mesenchymal stem cells (MSCs) for future clinical applications. The potential of eccDNA as a stable therapeutic target for senescence-related disorders warrants further investigation, particularly exploring chemically synthesized eccDNAs as transcription factor regulatory elements to reverse cellular senescence.

Identification of the Efficient Enhancer Elements in FVIII-Padua for Gene Therapy Study of Hemophilia A

Hemophilia A (HA) is a common X-linked recessive hereditary bleeding disorder. Coagulation factor VIII (FVIII) is insufficient in patients with HA due to the mutations in the F8 gene. The restoration of plasma levels of FVIII via both recombinant B-domain-deleted FVIII (BDD-FVIII) and B-domain-deleted F8 (BDDF8) transgenes was proven to be helpful. FVIII-Padua is a 23.4 kb tandem repeat mutation in the F8 associated with a high F8 gene expression and thrombogenesis. Here we screened a core enhancer element in FVIII-Padua for improving the F8 expression. In detail, we identified a 400 bp efficient enhancer element, C400, in FVIII-Padua for the first time. The core enhancer C400 extensively improved the transcription of BDDF8 driven by human elongation factor-1 alpha in HepG2, HeLa, HEK-293T and induced pluripotent stem cells (iPSCs) with different genetic backgrounds, as well as iPSCs-derived endothelial progenitor cells (iEPCs) and iPSCs-derived mesenchymal stem cells (iMSCs). The expression of FVIII protein was increased by C400, especially in iEPCs. Our research provides a novel molecular target to enhance expression of FVIII protein, which has scientific value and application prospects in both viral and nonviral HA gene therapy strategies.

The role of salinity on genome-wide DNA methylation dynamics in European sea bass gills

Epigenetic modifications, like DNA methylation, generate phenotypic diversity in fish and ultimately lead to adaptive evolutionary processes. Euryhaline marine species that migrate between salinity-contrasted habitats have received little attention regarding the role of salinity on whole-genome DNA methylation. Investigation of salinity-induced DNA methylation in fish will help to better understand the potential role of this process in salinity acclimation. Using whole-genome bisulfite sequencing, we compared DNA methylation patterns in European sea bass (Dicentrarchus labrax) juveniles in seawater and after freshwater transfer. We targeted the gill as a crucial organ involved in plastic responses to environmental changes. To investigate the function of DNA methylation in gills, we performed RNAseq and assessed DNA methylome-transcriptome correlations. We showed a negative correlation between gene expression levels and DNA methylation levels in promoters, first introns and first exons. A significant effect of salinity on DNA methylation dynamics with an overall DNA hypomethylation in freshwater-transferred fish compared to seawater controls was demonstrated. This suggests a role of DNA methylation changes in salinity acclimation. Genes involved in key functions as metabolism, ion transport and transepithelial permeability (junctional complexes) were differentially methylated and expressed between salinity conditions. Expression of genes involved in mitochondrial metabolism (tricarboxylic acid cycle) was increased, whereas the expression of DNA methyltransferases 3a was repressed. This study reveals novel links between DNA methylation, mainly in promoters and first exons/introns, and gene expression patterns following salinity change.

Analysis of BRCA1 and BRCA2 alternative splicing in predisposition to ovarian cancer

BACKGROUND

The mRNA splicing is regulated on multiple levels, resulting in the proper distribution of genes' transcripts in each cell and maintaining cell homeostasis. At the same time, the expression of alternative transcripts can change in response to underlying genetic variants, often missed during routine diagnostics.

AIM

The main aim of this study was to define the frequency of aberrant splicing in BRCA1 and BRCA2 genes in blood RNA extracted from ovarian cancer patients who were previously found negative for the presence of pathogenic alterations in the 25 most commonly analysed ovarian cancer genes, including BRCA1 and BRCA2.

MATERIAL AND METHODS

Frequency and spectrum of splicing alterations in BRCA1 and BRCA2 genes were analysed in blood RNA from 101 ovarian cancer patients and healthy controls (80 healthy women) using PCR followed by gel electrophoresis and Sanger sequencing. The expression of splicing events was examined using RT-qPCR.

RESULTS

We did not identify any novel, potentially pathogenic splicing alterations. Nevertheless, we detected six naturally occurring transcripts, named BRCA1ΔE9-10, BRCA1ΔE11, BRCA1ΔE11q, and BRCA2ΔE3, BRCA2ΔE12 and BRCA2ΔE17-18 of which three (BRCA1ΔE11q, BRCA1ΔE11 and BRCA2ΔE3) were significantly higher expressed in the ovarian cancer cohort than in healthy controls (p ≤ 0.0001).

CONCLUSIONS

This observation indicates that the upregulation of selected naturally occurring transcripts can be stimulated by non-genetic mechanisms and be a potential systemic response to disease progression and/or treatment. However, this hypothesis requires further examination.

The position of the longest intron is related to biological functions in some human genes

The evidence that introns can influence different levels of transfer of genetic information between DNA and the final product is increasing. Longer first introns were found to be a general property of eukaryotic gene structure and shown to contain a higher fraction of conserved sequence and different functional elements. Our work brings more precise information about the position of the longest introns in human protein-coding genes and possible connection with biological function and gene expression. According to our results, the position of the longest intron can be localized to the first third of introns in 64%, the second third in 19%, and the third in 17%, with notable peaks at the middle and last introns of approximately 5% and 6%, respectively. The median lengths of the longest introns decrease with increasing distance from the start of the gene from approximately 15,000 to 5,000 bp. We have shown that the position of the longest intron is in some cases linked to the biological function of the given gene. For example, DNA repair genes have the longest intron more often in the second or third. In the distribution of gene expression according to the position of the longest intron, tissue-specific profiles can be traced with the highest expression usually at the absolute positions of intron 1 and 2. In this work, we present arguments supporting the hypothesis that the position of the longest intron in a gene is another biological factor modulating the transmission of genetic information. The position of the longest intron is related to biological functions in some human genes.

Dnmt1/Tet2-mediated changes in Cmip methylation regulate the development of nonalcoholic fatty liver disease by controlling the Gbp2-Pparγ-CD36 axis

Dynamic alteration of DNA methylation leads to various human diseases, including nonalcoholic fatty liver disease (NAFLD). Although C-Maf-inducing protein (Cmip) has been reported to be associated with NAFLD, its exact underlying mechanism remains unclear. Here, we aimed to elucidate this mechanism in NAFLD in vitro and in vivo. We first identified alterations in the methylation status of the Cmip intron 1 region in mouse liver tissues with high-fat high-sucrose diet-induced NAFLD. Knockdown of DNA methyltransferase (Dnmt) 1 significantly increased Cmip expression. Chromatin immunoprecipitation assays of AML12 cells treated with oleic and palmitic acid (OPA) revealed that Dnmt1 was dissociated and that methylation of H3K27me3 was significantly decreased in the Cmip intron 1 region. Conversely, the knockdown of Tet methylcytosine dioxygenase 2 (Tet2) decreased Cmip expression. Following OPA treatment, the CCCTC-binding factor (Ctcf) was recruited, and H3K4me3 was significantly hypermethylated. Intravenous Cmip siRNA injection ameliorated NAFLD pathogenic features in ob/ob mice. Additionally, Pparγ and Cd36 expression levels were dramatically decreased in the livers of ob/ob mice administered siCmip, and RNA sequencing revealed that Gbp2 was involved. Gbp2 knockdown also induced a decrease in Pparγ and Cd36 expression, resulting in the abrogation of fatty acid uptake into cells. Our data demonstrate that Cmip and Gbp2 expression levels are enhanced in human liver tissues bearing NAFLD features. We also show that Dnmt1-Trt2/Ctcf-mediated reversible modulation of Cmip methylation regulates the Gbp2-Pparγ-Cd36 signaling pathway, indicating the potential of Cmip as a novel therapeutic target for NAFLD.

Genetic variations in AURORA cell cycle kinases are associated with glioblastoma multiforme

Glioblastoma multiforme (GBM) is the most frequent type of primary astrocytomas. We examined the association between single nucleotide polymorphisms (SNPs) in Aurora kinase A (AURKA), Aurora kinase B (AURKB), Aurora kinase C (AURKC) and Polo-like kinase 1 (PLK1) mitotic checkpoint genes and GBM risk by qPCR genotyping. In silico analysis was performed to evaluate effects of polymorphic biological sequences on protein binding motifs. Chi-square and Fisher statistics revealed a significant difference in genotypes frequencies between GBM patients and controls for AURKB rs2289590 variant (p = 0.038). Association with decreased GBM risk was demonstrated for AURKB rs2289590 AC genotype (OR = 0.54; 95% CI = 0.33-0.88; p = 0.015). Furthermore, AURKC rs11084490 CG genotype was associated with lower GBM risk (OR = 0.57; 95% CI = 0.34-0.95; p = 0.031). Bioinformatic analysis of rs2289590 polymorphic region identified additional binding site for the Yin-Yang 1 (YY1) transcription factor in the presence of C allele. Our results indicated that rs2289590 in AURKB and rs11084490 in AURKC were associated with a reduced GBM risk. The present study was performed on a less numerous but ethnically homogeneous population. Hence, future investigations in larger and multiethnic groups are needed to strengthen these results.

Partial F8 gene duplication (factor VIII Padua) associated with high factor VIII levels and familial thrombophilia

High coagulation factor VIII (FVIII) levels comprise a common risk factor for venous thromboembolism (VTE), but the underlying genetic determinants are largely unknown. We investigated the molecular bases of high FVIII levels in 2 Italian families with severe thrombophilia. The proband of the first family had a history of recurrent VTE before age 50 years, with extremely and persistently elevated FVIII antigen and activity levels (>400%) as the only thrombophilic defects. Genetic analysis revealed a 23.4-kb tandem duplication of the proximal portion of the F8 gene (promoter, exon 1, and a large part of intron 1), which cosegregated with high FVIII levels in the family and was absent in 103 normal controls. Targeted screening of 50 unrelated VTE patients with FVIII levels ≥250% identified a second thrombophilic family with the same F8 rearrangement on the same genetic background, suggesting a founder effect. Carriers of the duplication from both families showed a twofold or greater upregulation of F8 messenger RNA, consistent with the presence of open chromatin signatures and enhancer elements within the duplicated region. Testing of these sequences in a luciferase reporter assay pinpointed a 927-bp region of F8 intron 1 associated with >45-fold increased reporter activity in endothelial cells, potentially mediating the F8 transcriptional enhancement observed in carriers of the duplication. In summary, we report the first thrombophilic defect in the F8 gene (designated FVIII Padua) associated with markedly elevated FVIII levels and severe thrombophilia in 2 Italian families.

Intronic regions of the human coagulation factor VIII gene harboring transcription factor binding sites with a strong bias towards the short-interspersed elements

Increasing data show that intronic derived regulatory elements, such as transcription factor binding sites (TFBs), play key roles in gene regulation, and malfunction. Accordingly, characterizing the sequence context of the intronic regions of the human coagulation factor VIII (hFVIII) gene can be important. In this study, the intronic regions of the hFVIII gene were scrutinized based on in-silico methods. The results disclosed that these regions harbor a rich array of functional elements such as repetitive elements (REs), splicing sites, and transcription factor binding sites (TFBs). Among these elements, TFBs and REs showed a significant distribution and correlation to each other. This survey indicated that 31% of TFBs are localized in the intronic regions of the gene. Moreover, TFBs indicate a strong bias in the regions far from splice sites of introns with mapping to different REs. Accordingly, TFBs showed highly bias toward Short Interspersed Elements (SINEs), which in turn they covering about 12% of the total of REs. However, the distribution pattern of TFBs-REs showed different bias in the intronic regions, spatially into the Introns 13 and 25. The rich array of SINE-TFBs and CR1-TFBs were situated within 5′UTR of the gene that may be an important driving force for regulatory innovation of the hFVIII gene. Taken together, these data may lead to revealing intronic regions with the capacity to renewing gene regulatory networks of the hFVIII gene. On the other hand, these correlations might provide the novel idea for a new hypothesis of molecular evolution of the FVIII gene, and treatment of Hemophilia A which should be considered in future studies.

Comparative sequence analysis across Brassicaceae, regulatory diversity in KCS5 and KCS6 homologs from Arabidopsis thaliana and Brassica juncea, and intronic fragment as a negative transcriptional regulator

Intra- and epicuticular-waxes primarily comprising of very long chain aliphatic lipid (VLCFA), terpenoids and secondary metabolites such as sterol and flavonoids played a major role in successful colonization of terrestrial ecosystem by aquatic plants and are thus considered as a key evolutionary innovation. The key rate limiting step of Fatty Acid (FA) biosynthesis of condensation/elongation are catalyzed by the enzyme, β-ketoacyl coenzyme A synthase (KCS), part of FAE (Fatty Acid Elongase) complex. KCS6 has been shown to be responsible for elongation using C22 fatty acid as substrate and is considered essential for synthesis of VLCFA for cuticular waxes. Earlier studies have established KCS5 as a close paralog of KCS6 in Arabidopsis thaliana, albeit with non-redundant function. We subsequently established segmental duplication responsible for origin of KCS6-KCS5 paralogy which is exclusive to Brassicaceae. In the present study, we aim to understand impact of duplication on regulatory diversification and evolution, through sequence and functional analysis of cis-regulatory element of KCS5 and KCS6. High level of sequence variation leading to conservation of only the proximal end of the promoter corresponding to the core promoter was observed among Brassicaceae members; such high diversity was also revealed when sliding window analysis revealed only two to three phylogenetic footprints. Profiling of transcription factor binding sites (TFBS) across Brassicaceae shows presence of light, hormone and stress responsive motifs; a few motifs involved in tissue specific expression (Skn-1; endosperm) were also detected. Functional characterization using transcriptional fusion constructs revealed regulatory diversification when promoter activity of homologs from A. thaliana and Brassica juncea were compared. When subjected to 5-Azacytidine, altered promoter activity was observed, implying role of DNA methylation in transcriptional regulation. Finally, investigation of the role of an 87 bp fragment from first intron that is retained in a splice variant, revealed it to be a transcriptional repressor. This is a first report on comparative sequence and functional analysis of transcriptional regulation of KCS5 and KCS6; further studies are required before manipulation of cuticular waxes as a strategy for mitigating stress.

Friend or foe, the role of EGR-1 in cancer

Early growth response-1 (EGR-1), also termed NEFI-A and Krox-24, as a multi-domain protein is implicated in several vital physiological processes, including development, metabolism, cell growth and proliferation. Previous studies have implied that EGR-1 was producing in response to the tissue injury, immune response and fibrosis. Meanwhile, emerging studies stressed the pronounced correlation of EGR-1 and human cancers. Nevertheless, the intricate mechanisms of cancer-reduce EGR-1 alteration still poorly characterized. In the review, we evaluated the effects of EGR-1 in tumor cell proliferation, apoptosis, migration, invasion and tumor microenvironment, and then, we dwell on the intricate signaling pathways that EGR-1 involved in. The aberrantly expressed of EGR-1 in cancers are expected to provide a new cancer therapy strategy or a new marker for assessing treatment efficacy.

Breakpoint junction features of seven DMD deletion mutations

Duchenne muscular dystrophy is an inherited muscle wasting disease with severe symptoms and onset in early childhood. Duchenne muscular dystrophy is caused by loss-of-function mutations, most commonly deletions, within the DMD gene. Characterizing the junction points of large genomic deletions facilitates a more detailed model of the origins of these mutations and allows for a greater understanding of phenotypic variations associated with particular genotypes, potentially providing insights into the deletion mechanism. Here, we report sequencing of breakpoint junctions for seven patients with intragenic, whole-exon DMD deletions. Of the seven junction sequences identified, we found one instance of a “clean” break, three instances of microhomology (2–5 bp) at the junction site, and three complex rearrangements involving local sequences. Bioinformatics analysis of the upstream and downstream breakpoint regions revealed a possible role of short inverted repeats in the initiation of some of these deletion events.

Researchers in Australia have identified new examples of the genomic factors and mechanisms that lead to deletions linked with Duchenne muscular dystrophy (DMD). DMD is an inherited neuromuscular disease which causes progressive deterioration of muscles and, in some cases, intellectual impairment. Using samples from seven DMD patients, Niall Keegan of Murdoch University in Perth and colleagues sequenced the DNA left behind around the deletions in the DMD gene which cause the disease. They found one clean break, three sections with short repeated sequences, and three with more complex rearrangements. The diversity of these findings led them to suggest that the deletions resulted from a diversity of genomic factors and repair mechanisms. Future work could incorporate these findings into a model to predict where deletions will occur, expanding our understanding of DMD and its causes.

Matters of life and death: How estrogen and estrogen receptor binding to the immunoglobulin heavy chain locus may influence outcomes of infection, allergy, and autoimmune disease

Sex hormones are best known for their influences on reproduction, but they also have profound influences on the immune response. Examples of sex-specific differences include: (i) the relatively poor control of influenza virus infections in males compared to females, (ii) allergic asthma, an IgE-associated hypersensitivity reaction that is exacerbated in adolescent females compared to males, and (iii) systemic lupus erythematosus, a life-threatening autoimmune disease with a 9:1 female:male bias. Here we consider how estrogen and estrogen receptor α (ERα) may influence the immune response by modifying class switch recombination (CSR) and immunoglobulin expression patterns. We focus on ERα binding to enhancers (Eμ and the 3' regulatory region) and switch sites (Sµ and Sε) in the immunoglobulin heavy chain locus. Our preliminary data from ChIP-seq analyses of purified, activated B cells show estrogen-mediated changes in the positioning of ERα binding within and near Sµ and Sε. In the presence of estrogen, ERα is bound not only to estrogen response elements (ERE), but also to adenosine-cytidine (AC)-repeats and poly adenosine (poly A) sequences, in some cases within constant region gene introns. We propose that by binding these sites, estrogen and ERα directly participate in the DNA loop formation required for CSR. We further suggest that estrogen regulates immunoglobulin expression patterns and can thereby influence life-and-death outcomes of infection, hypersensitivity, and autoimmune disease.

Characterization and risk association of polymorphisms in Aurora kinases A, B and C with genetic susceptibility to gastric cancer development

BACKGROUND

Single nucleotide polymorphisms (SNPs) in genes encoding mitotic kinases could influence development and progression of gastric cancer (GC).

METHODS

Case-control study of nine SNPs in mitotic genes was conducted using qPCR. The study included 116 GC patients and 203 controls. In silico analysis was performed to evaluate the effects of polymorphisms on transcription factors binding sites.

RESULTS

The AURKA rs1047972 genotypes (CT vs. CC: OR, 1.96; 95% CI, 1.05-3.65; p = 0.033; CC + TT vs. CT: OR, 1.94; 95% CI, 1.04-3.60; p = 0.036) and rs911160 (CC vs. GG: OR, 5.56; 95% CI, 1.24-24.81; p = 0.025; GG + CG vs. CC: OR, 5.26; 95% CI, 1.19-23.22; p = 0.028), were associated with increased GC risk, whereas certain rs8173 genotypes (CG vs. CC: OR, 0.60; 95% CI, 0.36-0.99; p = 0.049; GG vs. CC: OR, 0.38; 95% CI, 0.18-0.79; p = 0.010; CC + CG vs. GG: OR, 0.49; 95% CI, 0.25-0.98; p = 0.043) were protective. Association with increased GC risk was demonstrated for AURKB rs2241909 (GG + AG vs. AA: OR, 1.61; 95% CI, 1.01-2.56; p = 0.041) and rs2289590 (AC vs. AA: OR, 2.41; 95% CI, 1.47-3.98; p = 0.001; CC vs. AA: OR, 6.77; 95% CI, 2.24-20.47; p = 0.001; AA+AC vs. CC: OR, 4.23; 95% CI, 1.44-12.40; p = 0.009). Furthermore, AURKC rs11084490 (GG + CG vs. CC: OR, 1.71; 95% CI, 1.04-2.81; p = 0.033) was associated with increased GC risk. A combined analysis of five SNPs, associated with an increased GC risk, detected polymorphism profiles where all the combinations contribute to the higher GC risk, with an OR increased 1.51-fold for the rs1047972(CT)/rs11084490(CG + GG) to 2.29-fold for the rs1047972(CT)/rs911160(CC) combinations. In silico analysis for rs911160 and rs2289590 demonstrated that different transcription factors preferentially bind to polymorphic sites, indicating that AURKA and AURKB could be regulated differently depending on the presence of particular allele.

CONCLUSIONS

Our results revealed that AURKA (rs1047972 and rs911160), AURKB (rs2241909 and rs2289590) and AURKC (rs11084490) are associated with a higher risk of GC susceptibility. Our findings also showed that the combined effect of these SNPs may influence GC risk, thus indicating the significance of assessing multiple polymorphisms, jointly. The study was conducted on a less numerous but ethnically homogeneous Bosnian population, therefore further investigations in larger and multiethnic groups and the assessment of functional impact of the results are needed to strengthen the findings.

Genetic association of complement component 2 variants with chronic hepatitis B in a Korean population

BACKGROUND & AIMS

Numerous single nucleotide polymorphisms associated with an increased risk of liver diseases, chronic hepatitis B and chronic hepatitis B-related hepatocellular carcinoma have been identified. In this study, we scrutinized the genetic effects of C2 variants, which were conflicting in previous results, on the risk of chronic hepatitis B in a Korean population.

METHODS

We genotyped 22 common C2 genetic variants of 977 chronic hepatitis B cases including 302 chronic hepatitis B-related hepatocellular carcinoma cases and 785 population controls. Statistical analysis was performed to examine the effects of genotype on the risk of chronic hepatitis B and chronic hepatitis B-related hepatocellular carcinoma.

RESULTS

Logistic regression analyses showed that six C2 single nucleotide polymorphisms had significant associations with the risk of chronic hepatitis B and chronic hepatitis B-related hepatocellular carcinoma among the Korean subjects. Stepwise analysis revealed that causal markers (rs9267665 and rs10947223) were identified among the C2 variants (stepwise P = 3.32 × 10^-9 and 2.04 × 10^-5 respectively). In further conditional analysis with previous chronic hepatitis B-associated loci, these two single nucleotide polymorphisms were independently associated with the risk of chronic hepatitis B. In addition, we investigated the ability of genetic risk scores combining 12 multi-chronic hepatitis B loci to predict the risk of chronic hepatitis B. Individuals with higher genetic risk scores showed increased risk for chronic hepatitis B.

CONCLUSIONS

Our results suggested that the C2 gene might be a susceptibility locus for chronic hepatitis B in Korean populations. The cumulative genetic effects may contribute to future etiological explanations for chronic hepatitis B.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

The intracellular pyrimidine 5'-nucleotidase NT5C3A is a negative epigenetic factor in interferon and cytokine signaling

The enzyme pyrimidine 5'-nucleotidase (NT5C3A), which mediates nucleotide catabolism, was previously thought to be restricted to blood cells. We showed that expression of the gene encoding NT5C3A was induced by type I interferons (IFNs) in multiple cell types and that NT5C3A suppressed cytokine production through inhibition of the nuclear factor κB (NF-κB) pathway. NT5C3A expression required both an intronic IFN-stimulated response element and the IFN-stimulated transcription factor IRF1. Overexpression of NT5C3A, but not of its catalytic mutants, suppressed IL-8 production by HEK293 cells. Whereas knockdown of NT5C3A enhanced tumor necrosis factor (TNF)-stimulated IL-8 production, it reduced the IFN-mediated suppression of Il8 expression. Overexpression of NT5C3A increased the abundance of NAD⁺ and the activation of the sirtuins SIRT1 and SIRT6, which are NAD⁺-dependent deacetylases. NT5C3A-stimulated sirtuin activity resulted in deacetylation of histone H3 and the NF-κB subunit RelA (also known as p65), both of which were associated with the proximal region of the Il8 promoter, thus repressing the transcription of Il8 Together, these data identify an anti-inflammatory pathway that depends on the catalytic activity of NT5C3A and functions as a negative feedback regulator of inflammatory cytokine signaling.

An evaluation of the evolution of the gene structure of dystroglycan

Background

Dystroglycan (DG) is an adhesion receptor complex composed of two non-covalently associated subunits, transcribed from a single gene. The extracellular α-DG is highly and heterogeneously glycosylated and binds with high affinity to laminins, and the transmembrane β-DG binds intracellular dystrophin. Multiple cellular functions have been proposed for DG, notwithstanding that its role in skeletal muscle appears central as demonstrated by both primary and secondary severe muscular dystrophic phenotypes collectively known as dystroglycanopathies. We recently analysed the molecular phylogeny of the DG core protein and identified the α/β interface, transmembrane and cytoplasmic domains of β-DG as the most conserved region. It was also identified that the IG2_MAT_NU region has been independently duplicated in multiple lineages.

Results

To understand the evolution of dystroglycan in more depth, we investigated dystroglycan gene structure in 35 species representative of the phyla in which dystroglycan has been identified (i.e., all metazoan phyla except Ctenophora). The gene structure of three exons and two introns is remarkably conserved. However, additional lineage-specific introns were identified, which interrupt the coding sequence at distinct points, were identified in multiple metazoan groups, most prominently in ecdysozoans.

Conclusions

A coding DNA sequence (CDS) intron that interrupts the encoding of the IG1 domain is universally conserved and this intron is longer in gnathostomes (jawed vertebrates) than in other metazoans. Lineage-specific gain of additional introns has occurred notably in ecdysozoans, where multiple introns interrupt the large 3′ exon. More limited intron gain has also occurred in placozoa, cnidarians, urochordates and the DG paralogues of lamprey and teleost fish.

Electronic supplementary material

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

Customized Array Comparative Genomic Hybridization Analysis of 25 Phosphatase-encoding Genes in Colorectal Cancer Tissues

BACKGROUND/AIM

Molecular mechanisms of alterations in protein tyrosine phosphatases (PTPs) genes in cancer have been previously described and include chromosomal aberrations, gene mutations, and epigenetic silencing. However, little is known about small intragenic gains and losses that may lead to either changes in expression or enzyme activity and even loss of protein function.

MATERIALS AND METHODS

The aim of this study was to investigate 25 phosphatase genes using customized array comparative genomic hybridization in 16 sporadic colorectal cancer tissues.

RESULTS

The analysis revealed two unique small alterations: of 2 kb in PTPN14 intron 1 and of 1 kb in PTPRJ intron 1. We also found gains and losses of whole PTPs gene sequences covered by large chromosome aberrations.

CONCLUSION

In our preliminary studies using high-resolution custom microarray we confirmed that PTPs are frequently subjected to whole-gene rearrangements in colorectal cancer, and we revealed that non-polymorphic intragenic changes are rare.

HPV16 integration probably contributes to cervical oncogenesis through interrupting tumor suppressor genes and inducing chromosome instability

Background

The integration of human papilloma virus (HPV) into host genome is one of the critical steps that lead to the progression of precancerous lesion into cancer. However, the mechanisms and consequences of such integration events are poorly understood. This study aims to explore those questions by studying high risk HPV16 integration in women with cervical intraepithelial neoplasia (CIN) and cervical squamous cell carcinoma (SCC).

Methods

Specifically, HPV integration status of 13 HPV16-infected patients were investigated by ligation-mediated PCR (DIPS-PCR) followed by DNA sequencing.

Results

In total, 8 HPV16 integration sites were identified inside or around genes associated with cancer development. In particular, the well-studied tumor suppressor genes SCAI was found to be integrated by HPV16, which would likely disrupt its expression and therefore facilitate the migration of tumor. On top of that, we observed several cases of chromosome translocation events coincide with HPV integration, which suggests the existence of chromosome instability. Additionally, short overlapping sequences were observed between viral derived and host derived fragments in viral-cellular junctions, indicating that integration was mediated by micro homology-mediated DNA repair pathway.

Conclusions

Overall, our study suggests a model in which HPV16 might contribute to oncogenesis not only by disrupting tumor suppressor genes, but also by inducing chromosome instability.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-016-0454-4) contains supplementary material, which is available to authorized users.

Genetic Variants in the Transcriptional Regulatory Region of the ALOX5AP gene and Susceptibility to Ischemic Stroke in Chinese Populations

No coding sequence variants of the ALOX5AP gene that lead to amino acid substitutions have been identified. A two-stage study design was used to explore the relationship between variants in the transcriptional regulatory region of ALOX5AP gene and ischemic stroke (IS) risk in Chinese populations. IS was determined using CT and/or MRI. First, 18 SNPs, located in the upstream promoter region of ALOX5AP gene, were genotyped in 200 IS patients and 200 controls. And one potential associated SNP (rs17222919) was identified (P = 0.005,OR = 0.623, 95% CI: 0.448~0.866). Next, another independent case-control cohort comprising 810 IS patients and 825 matched controls was recruited to investigate the role of rs17222919, rs9579646 polymorphisms and their haplotypes in IS risk. The G allele frequency of rs17222919 in the IS group was significantly lower than that in control group (P = 0.007, OR = 0.792, 95% CI: 0.669~0.937). T-A and G-A haplotypes were associated with IS (P = 0.001,OR = 1.282, 95% CI:1.100~1.495; P = 0.0001, OR = 0.712, 95% CI: 0.598~0.848; respectively). Our study providesevidence that rs17222919 is a potential genetic protective factor against IS. Furthermore, the T-A haplotype is a risk factor and the G-A haplotype is a protective factor against IS in Chinese population.

Construction of an integrative regulatory element and variation map of the murine Tst locus

BACKGROUND

Given the abundance of new genomic projects and gene annotations, researchers trying to pinpoint causal genetic variants are faced with a challenging task of how to efficiently integrate all current genomic information. The objective of the study was to develop an approach to integrate various genomic annotations for a recently positionally-cloned Tst gene (Thiosulfate Sulfur Transferase, synonym Rhodanese) responsible for the Fob3b2 QTL effect on leanness and improved metabolic parameters. The second aim was to identify and prioritize Tst genetic variants that may be causal for the phenotypic effects.

RESULTS

A bioinformatics approach was developed to integrate existing knowledge of regulatory elements of the Tst gene. The entire Tst locus along with flanking segments was sequenced between our unique polygenic mouse Fat and Lean strains that were generated by divergent selection on adiposity for over 60 generations. The bioinformatics-generated regulatory element map of the Tst locus was then combined with genetic variants between the Fat and Lean mice and with comparative analyses of polymorphisms across 17 mouse strains in order to prioritise likely causal polymorphisms. Two candidate regulatory variants were identified, one overlapping an evolutionary constrained Tst intronic element and the other residing in the seed region of a predicted 3'UTR miRNA binding site.

CONCLUSIONS

This study developed a map of regulatory elements for the Tst locus in mice and identified candidate genetic variants with increased causal likelihood. This map provides a basis for experimental validation and functional analyses of this novel candidate leanness and antidiabetic gene. Our methodological approach is of general utility for analyzing regulation of loci that have limited annotations and experimental evidence and for identifying candidate causal regulatory genetic variants in post-GWAS or post-QTL- cloning studies.

Phorbol ester-mediated re-expression of endogenous LAT adapter in J.CaM2 cells: a model for dissecting drivers and blockers of LAT transcription

Linker for activation of T cells (LAT) is a raft-associated, transmembrane adapter protein critical for T-cell development and function. LAT expression is transiently upregulated upon T-cell receptor (TCR) engagement, but molecular mechanisms conveying TCR signaling to enhanced LAT transcription are not fully understood. Here we found that a Jurkat subline J.CaM2, initially characterized as LAT deficient, conditionally re-expressed LAT upon the treatment with a protein kinase C activator, phorbol 12-myristate 13-acetate (PMA). We took advantage of the above observation for studying cis-elements and trans-acting factors contributing to the activation-induced expression of LAT. We identified a LAT gene region spanning nucleotide position −14 to +357 relative to the ATG start codon as containing novel cis-regulatory elements that were able to promote PMA-induced reporter transcription in the absence of the core LAT promoter. Interestingly, a point mutation in LAT intron 1, identified in J.CaM2 cells, downmodulated LAT promoter activity by 50%. Mithramycin A, a selective Sp1 DNA-binding inhibitor, abolished LAT expression upon PMA treatment as did calcium ionophore ionomycin (Iono) and valproic acid (VPA), widely used as an anti-epileptic drug. Our data introduce J.CaM2 cells as a model for dissecting drivers and blockers of activation induced expression of LAT.

Histone demethylase KDM4B regulates otic vesicle invagination via epigenetic control of Dlx3 expression

In vertebrate embryos, the histone demethylase KDM4B affects otic placode proliferation, intercellular adhesion, and invagination by directly regulating Dlx3 expression.

In vertebrates, the inner ear arises from the otic placode, a thickened swathe of ectoderm that invaginates to form the otic vesicle. We report that histone demethylase KDM4B is dynamically expressed during early stages of chick inner ear formation. A loss of KDM4B results in defective invagination and striking morphological changes in the otic epithelium, characterized by abnormal localization of adhesion and cytoskeletal molecules and reduced expression of several inner ear markers, including Dlx3. In vivo chromatin immunoprecipitation reveals direct and dynamic occupancy of KDM4B and its target, H3K9me3, at regulatory regions of the Dlx3 locus. Accordingly, coelectroporations of DLX3 or KDM4B encoding constructs, but not a catalytically dead mutant of KDM4B, rescue the ear invagination phenotype caused by KDM4B knockdown. Moreover, a loss of DLX3 phenocopies a loss of KDM4B. Collectively, our findings suggest that KDM4B play a critical role during inner ear invagination via modulating histone methylation of the direct target Dlx3.

Downregulation of GLS2 in glioblastoma cells is related to DNA hypermethylation but not to the p53 status

Human phosphate-activated glutaminase (GA) is encoded by two genes: GLS and GLS2. Glioblastomas (GB) usually lack GLS2 transcripts, and their reintroduction inhibits GB growth. The GLS2 gene in peripheral tumors may be i) methylation- controlled and ii) a target of tumor suppressor p53 often mutated in gliomas. Here we assessed the relation of GLS2 downregulation in GB to its methylation and TP53 status. DNA demethylation with 5-aza-2'-deoxycytidine restored GLS2 mRNA and protein content in human GB cell lines with both mutated (T98G) and wild-type (U87MG) p53 and reduced the methylation of CpG1 (promoter region island), and CpG2 (first intron island) in both cell lines. In cell lines and clinical GB samples alike, methylated CpG islands were detected both in the GLS2 promoter (as reported earlier) and in the first intron of this gene. CpG methylation of either island was absent in GLS2-expressing non-tumoros brain tissues. Screening for mutation in the exons 5-8 of TP53 revealed a point mutation in only one out of seven GB examined. In conclusion, aberrant methylation of CpG islands, appear to contribute to silencing of GLS2 in GB by a mechanism bypassing TP53 mutations. © 2015 Wiley Periodicals, Inc.

Polymorphisms in the AOX2 gene are associated with the rooting ability of olive cuttings

Different rooting ability candidate genes were tested on an olive cross progeny. Our results demonstrated that only the AOX2 gene was strongly induced. OeAOX2 was fully characterised and correlated to phenotypical traits. The formation of adventitious roots is a key step in the vegetative propagation of trees crop species, and this ability is under strict genetic control. While numerous studies have been carried out to identify genes controlling adventitious root formation, only a few loci have been characterised. In this work, candidate genes that were putatively involved in rooting ability were identified in olive (Olea europaea L.) by similarity with orthologs identified in other plant species. The mRNA levels of these genes were analysed by real-time PCR during root induction in high- (HR) and low-rooting (LR) individuals. Interestingly, alternative oxidase 2 (AOX2), which was previously reported to be a functional marker for rooting in olive cuttings, showed a strong induction in HR individuals. From the OeAOX2 full-length gene, alleles and effective polymorphisms were distinguished and analysed in the cross progeny, which were segregated based on rooting. The results revealed a possible correlation between two single nucleotide polymorphisms of OeAOX2 gene and rooting ability.

Comparative analyses of the cholinergic locus of ChAT and VAChT and its expression in the silkworm Bombyx mori

The cholinergic locus, which encodes choline acetyltransferase (ChAT) and vesicular acetylcholine transporter (VAChT), is specifically expressed in cholinergic neurons, maintaining the cholinergic phenotype. The organization of the locus is conserved in Bilateria. Here we examined the structure of cholinergic locus and cDNA coding for ChAT and VAChT in the silkworm, Bombyx mori. The B. mori ChAT (BmChAT) cDNA encodes a deduced polypeptide including a putative choline/carnitine O-acyltransferase domain and a conserved His residue required for catalysis. The B. mori VAChT (BmVAChT) cDNA encodes a polypeptide including a putative major facilitator superfamily domain and 10 putative transmembrane domains. BmChAT and BmVAChT cDNAs share the 5'-region corresponding to the first and second exon of cholinergic locus. Polymerase chain reaction analyses revealed that BmChAT and BmVAChT mRNAs were specifically expressed in the brain and segmental ganglia. The expression of BmChAT was detected 3 days after oviposition. The expression level was almost constant during the larval stage, decreased in the early pupal stage, and increased toward eclosion. The average ratios of BmChAT mRNA to BmVAChT mRNA in brain-subesophageal ganglion complexes were 0.54±0.10 in the larvae and 1.92±0.11 in adults. In addition, we examined promoter activity of the cholinergic locus and localization of cholinergic neurons, using a baculovirus-mediated gene transfer system. The promoter sequence, located 2kb upstream from the start of transcription, was essential for cholinergic neuron-specific gene õexpression. Cholinergic neurons were found in several regions of the brain and segmental ganglia in the larvae and pharate adults.

Distribution bias of the sequence matching between exons and introns in exon joint and EJC binding region in C. elegans

We propose a mechanism that there are matching relations between mRNA sequences and corresponding post-spliced introns, and introns play a significant role in the process of gene expression. In order to reveal the sequence matching features, Smith-Waterman local alignment method is used on C. elegans mRNA sequences to obtain optimal matched segments between exon-exon sequences and their corresponding introns. Distribution characters of matching frequency on exon-exon sequences and sequence characters of optimal matched segments are studied. Results show that distributions of matching frequency on exon-exon junction region have obvious differences, and the exon boundary is revealed. Distributions of the length and matching rate of optimal matched segments are consistent with sequence features of siRNA and miRNA. The optimal matched segments have special sequence characters compared with their host sequences. As for the first introns and long introns, matching frequency values of optimal matched segments with high GC content, rich CG dinucleotides and high λCG values show the minimum distribution in exon junction complex (EJC) binding region. High λCG values in optimal matched segments are main characters in distinguishing EJC binding region. Results indicate that EJC and introns have competitive and cooperative relations in the process of combining on protein coding sequences. Also intron sequences and protein coding sequences do have concerted evolution relations.

Novel FOXF1 deep intronic deletion causes lethal lung developmental disorder, alveolar capillary dysplasia with misalignment of pulmonary veins

Haploinsufficiency of FOXF1 causes an autosomal dominant neonatally lethal lung disorder, alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV). We identified novel 0.8-kb deletion within the 1.4-kb intron of FOXF1 in a deceased newborn diagnosed with ACDMPV. The deletion arose de novo on the maternal copy of the chromosome 16, and did not affect FOXF1 minigene splicing tested in lung fibroblasts. However, FOXF1 transcript level in the ACDMPV peripheral lung tissue was reduced by almost 40%. We found that, in an in vitro reporter assay, the FOXF1 intron exhibited moderate transcriptional enhancer activity, correlating with the presence of binding sites for expression regulators CTCF and CEBPB, whereas its truncated copy, which lost major CTCF and CEBPB-binding sites, inhibited the FOXF1 promoter. Our data further emphasize the importance of testing the non-protein coding regions of the genome currently not covered by diagnostic chromosomal microarray analyses or whole-exome sequencing.

Associations between KCNJ6 (GIRK2) gene polymorphisms and pain-related phenotypes

G-protein coupled inwardly rectifying potassium (GIRK) channels are effectors determining degree of analgesia experienced upon opioid receptor activation by endogenous and exogenous opioids. The impact of GIRK-related genetic variation on human pain responses has received little research attention. We used a tag single nucleotide polymorphism (SNP) approach to comprehensively examine pain-related effects of KCNJ3 (GIRK1) and KCNJ6 (GIRK2) gene variation. Forty-one KCNJ3 and 69 KCNJ6 tag SNPs were selected, capturing the known variability in each gene. The primary sample included 311 white patients undergoing total knee arthroplasty in whom postsurgical oral opioid analgesic medication order data were available. Primary sample findings were then replicated in an independent white sample of 63 healthy pain-free individuals and 75 individuals with chronic low back pain (CLBP) who provided data regarding laboratory acute pain responsiveness (ischemic task) and chronic pain intensity and unpleasantness (CLBP only). Univariate quantitative trait analyses in the primary sample revealed that 8 KCNJ6 SNPs were significantly associated with the medication order phenotype (P < .05); overall effects of the KCNJ6 gene (gene set-based analysis) just failed to reach significance (P = .054). No significant KCNJ3 effects were observed. A continuous GIRK Related Risk Score (GRRS) was derived in the primary sample to summarize each individual's number of KCNJ6 "pain risk" alleles. This GRRS was applied to the replication sample, which revealed significant associations (P < .05) between higher GRRS values and lower acute pain tolerance and higher CLBP intensity and unpleasantness. Results suggest further exploration of the impact of KCNJ6 genetic variation on pain outcomes is warranted.