SNP rs1049430 in the 3'-UTR of SH3GL2 regulates its expression: Clinical and prognostic implications in head and neck squamous cell carcinoma

Sex-Biased CD3ζ 3'-UTR SNP Increased Incidence Risk in Aplastic Anemia

Purpose

Aplastic anemia (AA) is a bone marrow failure syndrome with an unclear pathogenesis. Abnormal T cell immunity is one of the mechanisms involved in AA, and CD3ζ is an important signaling molecule for T cell activation. Single-nucleotide polymorphisms (SNPs) in CD3ζ 3'-untranslated region (3'-UTR) were associated with some immune-related disease occurrence and affect CD3ζ protein level. In this study, our aim was to analyze whether CD3ζ 3'-UTR SNPs were associated with AA susceptibility and had influence on CD3ζ protein level and provide new research data for exploring the pathogenesis of aplastic anemia.

Patients and Methods

We screened the genotypes of SNPs in 101 healthy individuals and 91 AA patients by PCR-RFLP and sequencing. In addition, the effect of specific CD3ζ 3'-UTR SNPs was analyzed by flow cytometry and dual luciferase assay.

Results

Four SNPs of CD3ζ 3'-UTR, 1184 C >G (rs3738212), 1292 delG (rs3831958), 1403 G >C (rs1052230) and 1410 A >T (rs1052231) were identified from Chinese healthy individuals and AA patients in which rs3738212 was not previously reported. Increased risk of AA was observed in female AA who with heterozygous genotype of linkage disequilibrium SNP (rs3831958, rs1052230 and rs1052231). Different genotypes of rs3738212 have sex-biases feature in AA, higher 1184 CC frequency in male AA and higher 1184 CG frequency in female AA. Furthermore, rs3738212 could upregulate CD3ζ protein level.

Conclusion

This study first identified sex-specific CD3ζ 3'-UTR SNPs that were associated with risk of AA. Our data also demonstrated that rs3738212 could upregulate CD3ζ protein level.

Impact of hyperglycemia on the expression of GLUT1 during oral carcinogenesis in rats

BACKGROUND

On the background of the epidemiological link between diabetes and oral cancer, the present study aimed to analyze the potential involvement of selected glucose transporters (GLUT1/GLUT3/GLUT4), if any, in such putative association.

METHODS AND RESULTS

Oral carcinogenesis was induced by 4-nitroquinoline N-oxide in 10 non-diabetic and 10 diabetic rats; 8 non-diabetic and 7 diabetic rats served as controls. Expressions of selected GLUTs at mRNA and protein levels were analyzed in oral tissue (normal/lesion) by quantitative real-time PCR and immunohistochemistry respectively. Premalignant lesions (hyperplasia/dysplasia/carcinoma-in-situ) appeared on tongues of carcinogen-treated animals. Significant increase of GLUT1mRNA level was seen from normal to lesion tongues, along increasing lesion grades (from hyperplasia/mild dysplasia to moderate/severe dysplasia) and in lesions induced under hyperglycemic condition than that induced under normoglycemic one; a similar trend was found in transcript variant-1 of GLUT1, but not in variant-2. GLUT3 and GLUT4 mRNA levels were comparable among lesions irrespective of grades and glycemic status. Concordant to mRNA level, overall expression of GLUT1 protein was higher in tongue lesions in presence of hyperglycemia than in absence of such condition; non-lesion portions of tongues exposed to carcinogen showed a similar trend. Moreover in carcinogen-treated groups, non-lesion and lesion portions of tongues under hyperglycemic condition showed predominantly membranous expression for GLUT1 which was again significantly higher than equivalent portions of tongue under normoglycemic condition.

CONCLUSION

Hyperglycemia seemed to favor GLUT1 over-expression and membrane localization of the protein during oral carcinogenesis. GLUT1 transcript variant-1 appeared to be more important than variant-2 in disease pathogenesis.

scReQTL: an approach to correlate SNVs to gene expression from individual scRNA-seq datasets

BACKGROUND

Recently, pioneering expression quantitative trait loci (eQTL) studies on single cell RNA sequencing (scRNA-seq) data have revealed new and cell-specific regulatory single nucleotide variants (SNVs). Here, we present an alternative QTL-related approach applicable to transcribed SNV loci from scRNA-seq data: scReQTL. ScReQTL uses Variant Allele Fraction (VAFRNA) at expressed biallelic loci, and corelates it to gene expression from the corresponding cell.

RESULTS

Our approach employs the advantage that, when estimated from multiple cells, VAFRNA can be used to assess effects of SNVs in a single sample or individual. In this setting scReQTL operates in the context of identical genotypes, where it is likely to capture RNA-mediated genetic interactions with cell-specific and transient effects. Applying scReQTL on scRNA-seq data generated on the 10 × Genomics Chromium platform using 26,640 mesenchymal cells derived from adipose tissue obtained from three healthy female donors, we identified 1272 unique scReQTLs. ScReQTLs common between individuals or cell types were consistent in terms of the directionality of the relationship and the effect size. Comparative assessment with eQTLs from bulk sequencing data showed that scReQTL analysis identifies a distinct set of SNV-gene correlations, that are substantially enriched in known gene-gene interactions and significant genome-wide association studies (GWAS) loci.

CONCLUSION

ScReQTL is relevant to the rapidly growing source of scRNA-seq data and can be applied to outline SNVs potentially contributing to cell type-specific and/or dynamic genetic interactions from an individual scRNA-seq dataset.

AVAILABILITY

https://github.com/HorvathLab/NGS/tree/master/scReQTL.

Estimating the Allele-Specific Expression of SNVs From 10× Genomics Single-Cell RNA-Sequencing Data

With the recent advances in single-cell RNA-sequencing (scRNA-seq) technologies, the estimation of allele expression from single cells is becoming increasingly reliable. Allele expression is both quantitative and dynamic and is an essential component of the genomic interactome. Here, we systematically estimate the allele expression from heterozygous single nucleotide variant (SNV) loci using scRNA-seq data generated on the 10×Genomics Chromium platform. We analyzed 26,640 human adipose-derived mesenchymal stem cells (from three healthy donors), sequenced to an average of 150K sequencing reads per cell (more than 4 billion scRNA-seq reads in total). High-quality SNV calls assessed in our study contained approximately 15% exonic and >50% intronic loci. To analyze the allele expression, we estimated the expressed variant allele fraction (VAF_RNA) from SNV-aware alignments and analyzed its variance and distribution (mono- and bi-allelic) at different minimum sequencing read thresholds. Our analysis shows that when assessing positions covered by a minimum of three unique sequencing reads, over 50% of the heterozygous SNVs show bi-allelic expression, while at a threshold of 10 reads, nearly 90% of the SNVs are bi-allelic. In addition, our analysis demonstrates the feasibility of scVAF_RNA estimation from current scRNA-seq datasets and shows that the 3′-based library generation protocol of 10×Genomics scRNA-seq data can be informative in SNV-based studies, including analyses of transcriptional kinetics.

Association of rs8444 polymorphism in the LASS2 3'-UTR and bladder cancer risk in Chinese population

The aim of the present study was to explore the correlations between single nucleotide polymorphisms in LASS2 gene 3'-untranslated regions and bladder cancer risk in Chinese population. We first performed PCR and sequence for LASS2-3'-UTR in 105 bladder cancer patients and 100 control subjects. Next, multivariate logistic regression analysis was used to determine the relationship between single nucleotide polymorphisms frequency and susceptibility of bladder cancer, and clinical features in 105 cases. In addition, survival curves and Cox Regression analysis were used to investigate the effect of single nucleotide polymorphisms on clinical outcome in 58 cases. Finally, quantitative reverse-transcription PCR and immunohistochemical were performed to explore the influence of single nucleotide polymorphisms on LASS2 expression. We found that a single nucleotide polymorphism (rs8444 C>T) located in the 3'-UTR of LASS2 was significantly associated with the risk of bladder cancer. We also showed the frequency of rs8444 T genotype was higher in bladder cancer group and correlated with the risk of clinical prognosis. Yet, there were no significant correlations between T/C allele frequencies and the distributions of rs8444 genotype and tumor-node-metastasis stage, histological grade and distant metastasis in bladder cancer. Furthermore, we demonstrated that rs8444 C>T could affect LASS2 expression by single nucleotide polymorphism-related mRNA stability. Our results showed that LASS2-3'-UTR rs8444 C>T polymorphism was significantly associated with the individual risk and the poor overall survival of bladder cancer, suggesting that rs8444 TT genotype maybe act as an independent risk factor of susceptibility and clinical prognosis for bladder cancer in Chinese population.

Comparative transcriptome analysis reveals genes related to the yolk ratio of duck eggs

Yolk ratio is an important production index in the salted duck egg industry. Yolk constituents are deposited during development of follicles. However, the molecular mechanism responsible for different yolk ratios in laying ducks remains elusive. In this study, Shaoxing ducks laying eggs with different yolk ratios were chosen for an analysis of liver and ovary transcriptome information. Twelve libraries were constructed and generated an average of 58.5 million clean reads per library, of which 69% of clean reads from liver and 65% of clean reads from ovary were mapped to a reference genome. Between cross-phenotype groups, a total of 250 and 230 differently expressed genes (DEGs) were identified in liver and ovary respectively, of which 101 and 50 DEGs respectively were characterized. Several DEGs were detected, among which HMGCS1, HMGCR, FDFT1, (DHCR7), (STARD4), CYP46A1 and LPIN3 are involved in cholesterol metabolism-related pathways; KIAA0319, STARD4, AP1S3, SH3GL2 and CAV2 are involved in vesicular transport in the liver; and ELOVL2 and PSD2 are involved in fatty acid elongation and endocytosis in the ovary. High yolk-ratio ducks had higher activity for cholesterol synthesis and molecular trafficking. The identification of candidate genes greatly advances the understanding of the genetic basis of the formation of different yolk ratios.

Genome-Wide Association Study Identifies a Genetic Prediction Model for Postoperative Survival in Patients with Hepatocellular Carcinoma

Background

As an important aspect of tumor heterogeneity, genetic variation may influence susceptibility and prognosis in different types of cancer. By exploring the prognostic value of genetic variation, this study aimed to establish a model for predicting postoperative survival and assessing the impact of variation on clinical outcomes in patients with hepatocellular carcinoma (HCC).

Material/Methods

A genome-wide association study of 367 patients with HCC was conducted to identify single nucleotide polymorphisms (SNPs) associated with prognosis. Identified predictors were further evaluated in 758 patients. Two prognostic models were established using Cox proportional hazards regression and Nomogram strategy, and validated in another 316 patients. The effect of the SNP rs2431 was analyzed in detail.

Results

A prognostic model including 5 SNPs (rs10893585, rs2431, rs34675408, rs6078460, and rs6766361) was established and exhibited high predictive accuracy for HCC prognosis. The panel combined with tumor node metastasis (TNM) stage resulted in a significantly higher c-index (0.723) than the individual c-index values. Stratified by the Nomogram prediction model, the median overall survival for the low-risk and high-risk groups were 100.1 versus 30.8 months (P<0.001) in the training set and 82.2 versus 22.5 months (P<0.001) in the validation set. A closer examination of rs2431 revealed that it may regulate the expression of FNDC3B by disrupting a microRNA-binding site.

Conclusions

This study established prediction models based on genetic factors alone or in combination with TNM stage for postoperative survival in patients with HCC, and identified FNDC3B as a potential therapeutic target for combating HCC metastasis.

Comparative molecular characterization of typical and exceptional responders in glioblastoma

Glioblastoma (GBM) is the most common and the deadliest type of primary brain tumor, with a median survival time of only 15 months despite aggressive treatment. Although most patients have an extremely poor prognosis, a relatively small number of patients survive far beyond the median survival time. Investigation of these exceptional responders has sparked a great deal of interest and is becoming an important focus in the field of cancer research. To investigate the molecular differences between typical and exceptional responders in GBM, comparative analyses of somatic mutations, copy number, methylation, and gene expression datasets from The Cancer Genome Atlas were performed, and the results of these analyses were integrated via gene ontology and pathway analyses to assess the functional significance of the differential aberrations. Less severe copy number loss of CDKN2A, lower expression of CXCL8, and FLG mutations are all associated with an exceptional response. Typical responders are characterized by upregulation of NF-κB signaling and of pro-inflammatory cytokines, while exceptional responders are characterized by upregulation of Alzheimer's and Parkinson's disease pathways as well as of genes involved in synaptic transmission. The upregulated pathways and processes in typical responders are consistently associated with more aggressive tumor phenotypes, while those in the exceptional responders suggest a retained ability in tumor cells to undergo cell death in response to treatment. With the upcoming launch of the National Cancer Institute's Exceptional Responders Initiative, similar studies with much larger sample sizes will likely become possible, hopefully providing even more insight into the molecular differences between typical and exceptional responders.

Structural Characterization and Association of Ovine Dickkopf-1 Gene with Wool Production and Quality Traits in Chinese Merino

Dickkopf-1 (DKK1) is an inhibitor of canonical Wnt signaling pathway and regulates hair follicle morphogenesis and cycling. To investigate the potential involvement of DKK1 in wool production and quality traits, we characterized the genomic structure of ovine DKK1, performed polymorphism detection and association analysis of ovine DKK1 with wool production and quality traits in Chinese Merino. Our results showed that ovine DKK1 consists of four exons and three introns, which encodes a protein of 262 amino acids. The coding sequence of ovine DKK1 and its deduced amino acid sequence were highly conserved in mammals. Eleven single nucleotide polymorphisms (SNPs) were identified within the ovine DKK1 genomic region. Gene-wide association analysis showed that SNP5 was significantly associated with mean fiber diameter (MFD) in the B (selected for long wool fiber and high-quality wool), PW (selected for high reproductive capacity, high clean wool yield and high-quality wool) and U (selected for long wool fiber with good uniformity, high wool yield and lower fiber diameter) strains (p < 4.55 × 10-3 = 0.05/11). Single Nucleotide Polymorphisms wide association analysis showed that SNP8 was significantly associated with MFD in A strain and fleece weight in A (selected for large body size), PM (selected for large body size, high reproductive capacity and high meat yield) and SF (selected for mean fiber diameter less than 18 μm and wool fiber length between 5 and 9 cm) strains (p < 0.05), SNP9 was significantly associated with curvature in B and U strains (p < 0.05) and SNP10 was significantly associated with coefficient of variation of fiber diameter in A, PW and PM strains and standard deviation of fiber diameter in A and PM strains (p < 0.05). The haplotypes derived from these 11 identified SNPs were significantly associated with MFD (p < 0.05). In conclusion, our results suggest that DKK1 may be a major gene controlling wool production and quality traits, also the identified SNPs (SNPs5, 8, 9 and 10) might be used as potential molecular markers for improving sheep wool production and quality in sheep breeding.

Double hit of NEMO gene in preeclampsia

The precise etiology of preeclampsia is unknown. Family studies indicate that both genetic and environmental factors influence its development. One of these factors is NFkB, whose activation depends on NEMO (NFkB essential modulator. This is the first study to investigate the association between the existence of single nucleotide variant of the NEMO gene and the appearance of preeclampsia. A total of 151 women (72 preeclamptic women and 79 controls) and their children were examined. Sanger sequencing was performed to identify variants in the NEMO gene in the preeclamptic mothers. The maternal identified variants were then sought in the studied groups of children, and in the maternal and child controls, using RFLP-PCR. Real-time RT-PCR was performed to assess NEMO gene expression in maternal blood, umbilical cord blood and placentas. The sequencing process indicated the existence of two different variants in the 3'UTR region of the NEMO gene of preeclamptic women (IKBKG:c.*368C>A and IKBKG:c.*402C>T). The simultaneous occurrence of the TT genotype in the mother and the TT genotype in the daughter or a T allele in the son increased the risk of preeclampsia development 2.59 fold. Additionally, we found that the configuration of maternal/fetal genotypes (maternal TT/ daughter TT or maternal TT/son T) of IKBKG:c.*402C/T variant is associated with the level of NEMO gene expression. Our results showed that, the simultaneous occurrence of the maternal TT genotype (IKBKG:c.*402C>T variants) and TT genotype in the daughter or T allele in the son correlates with the level of NEMO gene expression and increases the risk of preeclampsia development. Our observations may offer a new insight into the genetic etiology and pathogenesis of preeclampsia.

Loss of SH3GL2 promotes the migration and invasion behaviours of glioblastoma cells through activating the STAT3/MMP2 signalling

SH3GL2 (Src homology 3 (SH3) domain GRB2‐like 2) is mainly expressed in the central nervous system and regarded as a tumour suppressor in human glioma. However, the molecular mechanism of the SH3GL2 protein involved in malignant behaviours of human glioma has not been elucidated. In this study, we tried to investigate the role of SH3GL2 in glioma cell migration and invasion and explore its underlined molecular mechanism. Firstly, we discovered that the protein level of SH3GL2 was widely decreased in the human glioma patients, especially in high‐grade glioma tissues. Then, we determined the role of SH3GL2 in migration and invasion of glioma cells upon SH3GL2 knocking down and overexpressing. It was showed that knocking down of SH3GL2 promoted the migration and invasion of glioma cells, whereas overexpression of SH3GL2 inhibited them. Further study on molecular mechanism disclosed that silencing of SH3GL2 obviously activated the STAT3 (signal transducer and activator of transcription 3) signalling thereby promoting the expression and secretion of MMP2. On the contrary, overexpression of SH3GL2 had opposite effect. Taken together, the above results suggest that SH3GL2 suppresses migration and invasion behaviours of glioma cells through negatively regulating STAT3/MMP2 signalling and that loss of SH3GL2 may intensify the STAT3/MMP2 signalling thereby contributing to the migration and invasion of glioma cells.

The Influence of PSCA Gene Variation on Its Expression and Gastric Adenocarcinoma Susceptibility in the Northwest Chinese Population

Gastric adenocarcinoma (GAC) imposes a considerable health burden around the world. Gene variation in prostate stem cell antigen gene (PSCA) has been identified to be associated with GAC risk, while the results showed regional variation. To explore the influence of PSCA gene variation on its expression and GAC risk in the Northwest Chinese population, four single nucleotide polymorphisms (SNPs) of PSCA were genotyped in 476 GAC cases and 481 controls using MassARRAY system. Two SNPs of rs2294008 (C>T) and rs2976392 (G>A) were identified to be associated with GAC risk. rs2294008, rs2976392 and rs10216533 made up two statistically significant haplotypes (Hap-CGG and Hap-TAG). Additionally, PSCA expression was analyzed by quantitative real time PCR, immunohistochemistry and tissue microarray. The results showed that PSCA expression was decreased in GAC tissues compared with adjacent normal tissues. For normal tissues, PSCA expression was higher with Hap-TA than that with Hap-CG. For GAC tissues, the differentiation degree of Hap-TA was higher than that of Hap-CG. The expression distribution of PSCA in multiple human organs showed disparity. These results suggest that PSCA gene variation has a potential effect on its expression and GAC risk in the Northwest Chinese population.