Functional intronic ERCC1 polymorphism from regulomeDB can predict survival in lung cancer after surgery

The association of ARRB1 polymorphisms with response to antidepressant treatment in depressed patients

Introduction: β-arrestin 1, a protein encoded by ARRB1 involved in receptor signaling, is a potential biomarker for the response to antidepressant drug (ATD) treatment in depression. We examined ARRB1 genetic variants for their association with response following ATD treatment in METADAP, a cohort of 6-month ATD-treated depressed patients. Methods: Patients (n = 388) were assessed at baseline (M0) and after 1 (M1), 3 (M3), and 6 months (M6) of treatment for Hamilton Depression Rating Scale (HDRS) changes, response, and remission. Whole-gene ARRB1 variants identified from high-throughput sequencing were separated by a minor allele frequency (MAF)≥5%. Frequent variants (i.e., MAF≥5%) annotated by RegulomeDB as likely affecting transcription factor binding were analyzed using mixed-effects models. Rare variants (i.e., MAF<5%) were analyzed using a variant set analysis. Results: The variant set analysis of rare variants was significant in explaining HDRS score changes (T = 878.9; p = 0.0033) and remission (T = -1974.1; p = 0.034). Rare variant counts were significant in explaining response (p = 0.016), remission (p = 0.022), and HDRS scores at M1 (p = 0.0021) and M3 (p=<0.001). rs553664 and rs536852 were significantly associated with the HDRS score (rs553664: p = 0.0055 | rs536852: p = 0.046) and remission (rs553664: p = 0.026 | rs536852: p = 0.012) through their interactions with time. At M6, significantly higher HDRS scores were observed in rs553664 AA homozygotes (13.98 ± 1.06) compared to AG heterozygotes (10.59 ± 0.86; p = 0.014) and in rs536852 GG homozygotes (14.88 ± 1.10) compared to AG heterozygotes (11.26 ± 0.95; p = 0.0061). Significantly lower remitter rates were observed in rs536852 GG homozygotes (8%, n = 56) compared to AG heterozygotes (42%, n = 105) at M6 (p = 0.0018). Conclusion: Our results suggest ARRB1 variants may influence the response to ATD treatment in depressed patients. Further analysis of functional ARRB1 variants and rare variant burden in other populations would help corroborate our exploratory analysis. β-arrestin 1 and genetic variants of ARRB1 may be useful clinical biomarkers for clinical improvement following ATD treatment in depressed individuals. Clinical Trial Registration: clinicaltrials.gov; identifier NCT00526383.

Beyond the Big Five: Investigating Myostatin Structure, Polymorphism and Expression in Camelus dromedarius

Myostatin, a negative regulator of skeletal muscle mass in animals, has been shown to play a role in determining muscular hypertrophy in several livestock species, and a high degree of polymorphism has been previously reported for this gene in humans and cattle. In this study, we provide a characterization of the myostatin gene in the dromedary (Camelus dromedarius) at the genomic, transcript and protein level. The gene was found to share high structural and sequence similarity with other mammals, notably Old World camelids. 3D modeling highlighted several non-conservative SNP variants compared to the bovine, as well as putative functional variants involved in the stability of the myostatin dimer. NGS data for nine dromedaries from various countries revealed 66 novel SNPs, all of them falling either upstream or downstream the coding region. The analysis also confirmed the presence of three previously described SNPs in intron 1, predicted here to alter both splicing and transcription factor binding sites (TFBS), thus possibly impacting myostatin processing and/or regulation. Several putative TFBS were identified in the myostatin upstream region, some of them belonging to the myogenic regulatory factor family. Patterns of SNP distribution across countries, as suggested by Bayesian clustering of the nine dromedaries using the 69 SNPs, pointed to weak geographic differentiation, in line with known recurrent gene flow at ancient trading centers along caravan routes. Myostatin expression was investigated in a set of 8 skeletal muscles, both at transcript and protein level, via Digital Droplet PCR and Western Blotting, respectively. No significant differences were observed at the transcript level, while, at the protein level, the only significant differences concerned the promyostatin dimer (75 kDa), in four pair-wise comparisons, all involving the tensor fasciae latae muscle. Beside the mentioned band at 75 kDa, additional bands were observed at around 40 and 25 kDa, corresponding to the promyostatin monomer and the active C-terminal myostatin dimer, respectively. Their weaker intensity suggests that the unprocessed myostatin dimers could act as important reservoirs of slowly available myostatin forms. Under this assumption, the sequential cleavage steps may contribute additional layers of control within an already complex regulatory framework.

Innovative strategies for annotating the "relationSNP" between variants and molecular phenotypes

Characterizing how variation at the level of individual nucleotides contributes to traits and diseases has been an area of growing interest since the completion of sequencing the first human genome. Our understanding of how a single nucleotide polymorphism (SNP) leads to a pathogenic phenotype on a genome-wide scale is a fruitful endeavor for anyone interested in developing diagnostic tests, therapeutics, or simply wanting to understand the etiology of a disease or trait. To this end, many datasets and algorithms have been developed as resources/tools to annotate SNPs. One of the most common practices is to annotate coding SNPs that affect the protein sequence. Synonymous variants are often grouped as one type of variant, however there are in fact many tools available to dissect their effects on gene expression. More recently, large consortiums like ENCODE and GTEx have made it possible to annotate non-coding regions. Although annotating variants is a common technique among human geneticists, the constant advances in tools and biology surrounding SNPs requires an updated summary of what is known and the trajectory of the field. This review will discuss the history behind SNP annotation, commonly used tools, and newer strategies for SNP annotation. Additionally, we will comment on the caveats that distinguish approaches from one another, along with gaps in the current state of knowledge, and potential future directions. We do not intend for this to be a comprehensive review for any specific area of SNP annotation, but rather it will be an excellent resource for those unfamiliar with computational tools used to functionally characterize SNPs. In summary, this review will help illustrate how each SNP annotation method impacts the way in which the genetic and molecular etiology of a disease is explored in-silico.

Association of Smoking and XPG, CYP1A1, OGG1, ERCC5, ERCC1, MMP2, and MMP9 Gene Polymorphisms with the early detection and occurrence of Laryngeal Squamous Carcinoma

A total of 200 smoking-related laryngeal carcinoma patients with pathology confirmation from the Eye and ENT Hospital and 190 high-risk smokers were included in a survey. All of the participants had a smoking index greater than 400 (cigarettes/day*year.) We obtained data on clinical and baseline characteristics, and peripheral blood was obtained and subjected to DNA extraction to analyse the correlation between smoking and the occurrence of laryngeal carcinoma. We selected candidate genes and SNP fragments that were found to be closely associated with smoking-related tumours in preliminary studies. The selected candidate genes were XPG, CYP1A1, OGG1, ERCC5, ERCC1, MMP2, and MMP9. We then performed SNP sequencing using Sequenom SNP detection technology. Target genes and single nucleotide polymorphism (SNP) fragments were evaluated to analyse the correlation between genotype or allele and smoking-related laryngeal carcinoma and to identify susceptibility genes related to laryngeal carcinoma. The results included four main findings: (1) The smoking index differed significantly between laryngeal cancer patients and control subjects (P=0.0035). The risk of laryngeal cancer was increased among individuals with a smoking index greater than 600 cigarettes/day*year (P=0.03). (2) The smoking index was significantly correlated with the T, N and clinical stages (P<0.05). (3) The polymorphisms CYP1A1-rs1048943, rs4646421, and rs4646422 and MMP9-rs17577 were significantly associated with laryngeal squamous carcinoma (P<0.05). (4) After stratifying the subjects by smoking degree, the GT genotype of ERCC1-rs2298881 was associated with a significantly greater risk of laryngeal carcinoma among heavy smokers (P=0.04). The results suggest that smoking plays an important role in the occurrence and development of laryngeal squamous carcinoma; CYP1A1 and MMP9 might be susceptibility gene SNPs for smoking-related laryngeal carcinoma, and ERCC1 might play an important role in heavy smokers. The results of this study might help identify an early marker for the detection and prevention of laryngeal carcinoma.

Regulatory variants in cancer-related pathway genes predict survival of patients with surgically resected non-small cell lung cancer

BACKGROUND

We conducted this study to identify genetic variants in cancer-related pathway genes which can predict prognosis of NSCLC patients after surgery, using a comprehensive list of regulatory single nucleotide polymorphisms (SNPs) prioritized by RegulomeDB.

METHOD

A total of 509 potentially functional SNPs in cancer-related pathway genes selected from RegulomeDB were evaluated. These SNPs were analyzed in a discovery set (n=354), and a replication study was performed in an independent set (n=772). The association of the SNPs with overall survival (OS) and disease-free survival (DFS) were analyzed.

RESULTS

In the discovery set, 76 SNPs were significantly associated with OS or DFS. Among the 76 SNPs, the association was consistently observed for 5 SNPs (ERCC1 rs2298881C>A, BRCA2 rs3092989G>A, NELFE rs440454C>T, PPP2R4 rs2541164G>A, and LTBP4 rs3786527G>A) in the validation set. In combined analysis, ERCC1 rs2298881C>A, BRCA2 rs3092989, NELFE rs440454C>T, and PPP2R4 rs2541164G>A were significantly associated with OS and DFS (adjusted HR ·aHR· for OS=1.46, 0.62, 078, and 0.76, respectively; P=0.003, 0.002, 0.007, and 0.003 respectively; and aHR for DFS=1.27, 0.69, 0.86, and 0.82, respectively; P=0.02, 0.002, 0.03, and 0.008, respectively). The LTBP4 rs3786527G>A was significantly associated with better OS (aHR=0.75; P=0.003).

CONCLUSION

Our results suggest that five SNPs in the cancer-related pathway genes may be useful for the prediction of the prognosis in patients with surgically resected NSCLC.

MAPK and ERK polymorphisms are associated with PCOS risk in Chinese women

In this case-control study, we analyzed the association between eight RegulomeDB-annotated single nucleotide polymorphisms (SNPs) in the MEK1, MEK2, ERK1 and ERK2 genes and polycystic ovarian syndrome (PCOS). Logistic regression analysis demonstrated that MEK1 rs12050732 (OR = 1.29 [95%CI: 1.06-1.58], P = 0.012), ERK2 rs2266966 (OR = 0.81 [95%CI: 0.67-0.99], P = 0.040) and ERK2 rs5999521 (OR = 0.66 [95%CI: 0.51-0.86], P = 0.002) were associated with PCOS risk without adjusting for age and body mass index. Moreover, PCOS risk increased with allele dosage when these three polymorphisms were combined (P_trend = 0.001). These findings suggest that genetic variants in key MAPK and ERK genes contribute to PCOS risk in Chinese women.

Polymorphism in ERCC1 confers susceptibility of coronary artery disease and severity of coronary artery atherosclerosis in a Chinese Han population

Excision repair cross-complementing 1 (ERCC1) gene encodes ERCC1 protein, which is mainly responsible for the repair of DNA damage in different diseases including coronary artery atherosclerosis by acting as a rate-limiting element in nucleotide excision repair (NER). Using a three-stage case-control study with 3037 coronary artery disease (CAD) patients and 3002 controls, we investigated associations of three single nucleotide polymorphisms (SNPs) with CAD risk and severity of coronary artery atherosclerosis in Chinese Han population. In the discovery set, the variant allele T of rs11615 was significantly associated with higher CAD risk (adjusted OR = 1.27, P = 0.006) and severity of coronary artery atherosclerosis (adjusted OR = 1.54, P = 0.003). These associations were more remarkable in the merged set (adjusted OR = 1.23, P = 8 × 10^-6 for CAD risk; adjusted OR = 1.36, P = 4.3 × 10^-5 for severity of coronary artery atherosclerosis). And the expression level of ERCC1 was significantly higher in CAD cases than controls. Multiplicative interactions among SNP rs11615, alcohol drinking, history of T2DM, and history of hyperlipidemia could increase 5.06-fold risk of CAD (P = 1.59 × 10^-9). No significant association of rs2298881 and rs3212986 with CAD risk was identified. Taken together, SNP rs11615 in ERCC1 gene might confer susceptibility to CAD and severity of coronary atherosclerosis in a Chinese Han population.

A practical guide to filtering and prioritizing genetic variants

Next-generation sequencing (NGS) of whole genomes and exomes is a powerful tool in biomedical research and clinical diagnostics. However, the vast amount of data produced by NGS introduces new challenges and opportunities, many of which require novel computational and theoretical approaches when it comes to identifying the causal variant(s) for a disease of interest. While workflows and associated software to process raw data and produce high-confidence variant calls have significantly improved, filtering tens of thousands of candidates to identify a subset relevant to a specific study is still a complex exercise best left to bioinformaticists. However, as this prioritization procedure requires biological/biomedical reasoning, biologists and clinicians are increasingly motivated to handle the task themselves. Here, we describe a set of guidelines, tools, and online resources that can be used to identify functional variants from whole-genome and whole-exome variant calls and then prioritize these variants with potential associations to phenotypes of interest. Insights gained from a recently published analysis of protein-coding gene variation in >60,000 humans by the Exome Aggregation Consortium (ExAC) are also taken into account.

Genetic variants of SOX9 contribute to susceptibility of gliomas among Chinese population

Gliomas make up about 80% of all malignant brain tumors, and cause serious public health problem. Genetic factors and environmental factors jointly caused the development of gliomas, and understanding of the genetic basis is a key component of preventive oncology. However, most genetic factors underlying carcinogenesis of gliomas remain largely unclear. In current study, we systematically evaluated whether genetic variants of SOX9 gene, a transcription factor that plays a central role in the development and differentiation of tumors, contribute to susceptibility of gliomas among Chinese population using a two-stage, case-control study. Results showed that SOX9 rs1042667 was significant associated with increased gliomas risk after adjusted by age, gender, family history of cancer, smoking status and alcohol status (Allele C vs A: OR=1.25; 95% CI=1.11-1.40; P=1.2×10-4). Compared with the carriers of genotype AA, both those of genotype AC (OR=1.37; 95% CI=1.13-1.66) and CC (OR=1.53; 95% CI=1.22-1.91) had significantly increased gliomas risk. This should be the first genetic association study which aims to evaluated the association between genetic variants of SOX9 and susceptibility of gliomas. Additional functional and association studies with different ethnic groups included are needed to further confirm our results.