Variants In Regulatory Regions Research Articles

Genetic variant of miR-4293 rs12220909 is associated with susceptibility to non-small cell lung cancer in a Chinese Han population.

Non-small cell lung cancer is one of the most common cancers and the leading cause of cancer death worldwide. Genetic variants in regulatory regions of some miRNAs might be involved in non-small cell lung cancer susceptibility and survival. rs12220909 (G/C) genetic polymorphism in miR-4293 has been shown to be associated with decreased risk of esophageal squamous cell carcinoma. However, the influence of rs12220909 genetic variation on non-small cell lung cancer susceptibility has not been reported. In order to evaluate the potential association between miR-4293 rs12220909 and non-small cell lung cancer risk in a Chinese population, we performed a case-control study among 998 non-small cell lung cancer cases and 1471 controls. The data shows that miR-4293 rs12220909 was significantly associated with decreased susceptibility to non-small cell lung cancer (GC vs.GG: OR = 0.681, 95%CI = 0.555–0.835, P = 2.19E-4; GG vs. GC+CC: OR = 0.687, 95%CI = 0.564–0.837, P = 1.95E-4), which indicates that rs12220909 in miR-4293 may play a significant role in the development of non-small cell lung cancer.

Immunoseq: the identification of functionally relevant variants through targeted capture and sequencing of active regulatory regions in human immune cells.

BackgroundThe observation that the genetic variants identified in genome-wide association studies (GWAS) frequently lie in non-coding regions of the genome that contain cis-regulatory elements suggests that altered gene expression underlies the development of many complex traits. In order to efficiently make a comprehensive assessment of the impact of non-coding genetic variation in immune related diseases we emulated the whole-exome sequencing paradigm and developed a custom capture panel for the known DNase I hypersensitive site (DHS) in immune cells – “Immunoseq”.ResultsWe performed Immunoseq in 30 healthy individuals where we had existing transcriptome data from T cells. We identified a large number of novel non-coding variants in these samples. Relying on allele specific expression measurements, we also showed that our selected capture regions are enriched for functional variants that have an impact on differential allelic gene expression. The results from a replication set with 180 samples confirmed our observations.ConclusionsWe show that Immunoseq is a powerful approach to detect novel rare variants in regulatory regions. We also demonstrate that these novel variants have a potential functional role in immune cells.Electronic supplementary materialThe online version of this article (doi:10.1186/s12920-016-0220-7) contains supplementary material, which is available to authorized users.

Abstract 1859: Bioinformatic and experimental evaluation of regulatory variants in breast cancer susceptibility genes

Abstract Many germline DNA sequence variants in gene regulatory elements are associated with an increased risk of cancer. These include large genomic deletions of the promoter of the BRCA1 gene, through to single nucleotide variants in the promoter, UTRs and long-range enhancers of multiple other cancer related genes. Next generation sequence analysis of early onset and familial breast cancer cases is currently identifying an escalating number of variants in non-coding regions of the genome, however for most of these the significance in unknown. This project aims to determine the functional and clinical significance of breast cancer associated variants in regulatory regions of BRCA1 and BRCA2, as part of an international collaborative project arising from the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium. To date over two hundred variants mapping to promoter, UTR and deep intronic regions have been identified in early onset or familial breast cancer patients that have no identifiable coding or splicing changes in BRCA1 or BRCA2. Through a pipeline of bioinformatics analyses using ENCODE datasets, these variants have been prioritized for experimental evaluation. This has lead to the identification of several variants that significantly alter the regulation of BRCA1 or BRCA2 expression, through mechanisms including altered binding of transcription factors and microRNAs. These studies will lay the groundwork for comprehensive statistical analyses, in which bioinformatics and experimental data is combined with clinical and genetic data to establish multifactorial risk prediction models that can ultimately be used in a clinical setting. Citation Format: Jan Sevcik, Leslie Burke, Gaetana Gambino, Brooke L. Brewster, Emma Tudini, Philip J. Whiley, Siranoush Manoukian, ENIGMA Consortium, Thomas van Overeem Hansen, Marta Santamariña Pena, Ana Vega, Maria A. Caligo, Paolo Radice, Paolo Peterlongo, Etienne Rouleau, Amanda B. Spurdle, Melissa A. Brown. Bioinformatic and experimental evaluation of regulatory variants in breast cancer susceptibility genes. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 1859.

Genetic variants in regulatory regions of microRNAs are associated with lung cancer risk.

Genetic variants in regulatory regions of some miRNAs might be associated with lung cancer risk and survival. We performed a case-control study including 1341 non-small cell lung cancer (NSCLC) cases and 1982 controls to evaluate the associations of 7 potentially functional polymorphisms in several differently expressed miRNAs with NSCLC risk. Each SNP was also tested for the association with overall survival of 1001 NSCLC patients. We identified that rs9660710 in miR-200b/200a/429 cluster and rs763354 in miR-30a were significantly associated with NSCLC risk [odds ratio (OR) = 1.17, 95% confidence interval (CI) = 1.06–1.30, P = 0.002; OR = 0.88, 95% CI = 0.80–0.98, P = 0.017; respectively]. However, no significant association between variants and NSCLC death risk was observed in survival analysis. Functional annotation showed that both rs9660710 and rs763354 were located in regulatory elements in lung cancer cells. Compared to normal tissues, miR-200a-3p, miR-200a-5p, miR-200b-3p, miR-200b-5p and miR-429 were significantly increased in The Cancer Genome Atlas (TCGA) Lung Adenocarcinoma (LUAD) tumors, whereas miR-30a-3p and miR-30a-5p were significantly decreased in tumors (all P < 0.05). Furthermore, we observed that rs9660710 is an expression quantitative trait locus (eQTL) or methylation eQTL for miR-429 expression in TCGA normal tissues. Our results indicated that rs9660710 in miR-200b/200a/429 cluster and rs763354 in miR-30a might modify the susceptibility to NSCLC.

Functional characterization of annexin A5 gene promoter allelic variants

BackgroundM1 and M2 haplotypes are defined by 4 consecutive allelic variants in regulatory regions of the annexin A5 gene and have been found to reduce promoter activity. To date, no research has been carried out to investigate differential and individual impact each of the allelic variants has on promoter activity. In the current study, we functionally characterized the M1 and M2 haplotype allelic variants (c.-467G>A, c.-448A>C, c.-422T>C, c.-373G>A). We also characterized two other allelic variants located in the same regulatory region (c.-628C>T, c.-302T>G). Materials and methodsTheir impact on the ANXA5 promoter activity was examined using a luciferase reporter assay in BeWo cells. Electrophoretic mobility shift assay with probes centered around each polymorphism was used to examine the binding ability of the allelic variants to nuclear proteins from BeWo cells. ResultsOnly the c.-467G>A and c.-628C>T allelic variants influenced the activity of the ANXA5 promoter, as measured by luciferase activity. Differential specific interactions with nuclear proteins were obtained for all allelic variants, except for the c.-302T>G, indicating that these polymorphisms could have an impact on the ANXA5 expression. ConclusionsWe have functionally characterized allelic variants in the ANXA5 promoter, both alone and in combinations, and the results suggest that combinations of several individual variants contribute to modulate the ANXA5 transcriptional activity, most likely through binding of nuclear factors. These results provide new knowledge and insight into the mechanisms underlying the regulation of annexin A5 levels in healthy controls.

Sequence analysis of candidate genes in two Roma families with severe tooth agenesis

Selective tooth agenesis is the most common congenital disorder affecting the formation of dentition in humans. Both its forms (hypodontia and more severe oligodontia) can be found either in isolated form and they can be associated with systemic condition (syndromic tooth agenesis). In addition to previously known genes (PAX9, MSX1 and AXIN2) mutations in EDA, EDARADD and WNT10 gene were recently found to be involved in isolated forms of tooth agenesis. The objective of this study was to characterize the phenotype of affected members in two large families of Roma origin segregating severe isolated tooth agenesis with very variable phenotype and to perform mutation analysis of seven genes with aim to find causal mutation. 26 family members were clinically examined and coding regions of seven genes (MSX1, PAX9, AXIN2, EDA, EDAR, EDARADD and WNT10A) were sequenced. With exclusion of third molars, average number of missing teeth was 8.2 ? 4.9 in family 1 and 7.1 ? 2.3 in family 2. The most frequently missing teeth were maxillary lateral incisors and first premolars and mandibular central incisors. Sequencing revealed four potentially damaging variants (g.Ala40Gly in MSX1, g.Ala240Pro in PAX9, g.Pro50Ser in AXIN2 and g.Met9Ile in EDARADD); however, none of them was present in all affected family members. Variable phenotype in both families examined in this study is in favour of heterogeneous genetic cause of tooth agenesis in these families: possible interaction of several defected genes, sequence variants in regulatory regions and additional environmental factors is assumed.

Altered Gene Expression Associated with microRNA Binding Site Polymorphisms.

Allele-specific gene expression associated with genetic variation in regulatory regions can play an important role in the development of complex traits. We hypothesized that polymorphisms in microRNA (miRNA) response elements (MRE-SNPs) that either disrupt a miRNA binding site or create a new miRNA binding site can affect the allele-specific expression of target genes. By integrating public expression quantitative trait locus (eQTL) data, miRNA binding site predictions, small RNA sequencing, and Argonaute crosslinking immunoprecipitation (AGO-CLIP) datasets, we identified genetic variants that can affect gene expression by modulating miRNA binding efficiency. We also identified MRE-SNPs located in regions associated with complex traits, indicating possible causative mechanisms associated with these loci. The results of this study expand the current understanding of gene expression regulation and help to interpret the mechanisms underlying eQTL effects.

Lactoferrin gene promoter variants and their association with clinical and subclinical mastitis in indigenous and crossbred cattle.

Lactoferrin (Lf) gene promoter was screened for the presence of single nucleotide polymphism in indigenous and crossbred cattle from North India and to evaluate its association with Mastitis. Study revealed the presence of genetic variation in regulatory region of bovine Lactoferrin gene using PCR-RFLP technique. Three genotypes namely GG, GH and HH were identified. A single nucleotide change, from guanine to adenine at 25th position was found to be significantly associated (p<0.05) with clinical mastitis in indigenous Sahiwal and crossbred Karan Fries cattle maintained at organised herd of National Dairy Research Institute, Karnal. A non-significant association was observed between subclinical mastitis, somatic cell score (SCS), and GG genotype in Karan Fries cattle, however, a lower SCS was observed in animals having GG genotype. Overall a lower incidence of clinical mastitis was recorded in those animals having GG genotype of Lf in Sahiwal and Karan Fries (KF) cattle. The SNP identified in the promoter region may effect expression lactoferrin protein, which may lead to different levels of antibacterial and anti-inflammatory activity of Lf gene. Results from this study indicated the probable role played by Lactoferrin promoter to serve as candidate gene for mastitis susceptibility among indigenous and crossbred milch cattle.

Inducible nitric oxide synthase (iNOS) regulatory region variation in non-human primates

Inducible nitric oxide synthase (iNOS) is an enzyme that plays a key role in intracellular immune response against respiratory infections. Since various species of nonhuman primates exhibit different levels of susceptibility to infectious respiratory diseases, and since variation in regulatory regions of genes is thought to play a key role in expression levels of genes, two candidate regulatory regions of iNOS were mapped, sequenced, and compared across five species of nonhuman primates: African green monkeys (Chlorocebus sabaeus), pig-tailed macaques (Macaca nemestrina), cynomolgus macaques (Macaca fascicularis), Indian rhesus macaques (Macaca mulatta), and Chinese rhesus macaques (M. mulatta). In addition, we conducted an in silico analysis of the transcription factor binding sites associated with genetic variation in these two candidate regulatory regions across species. We found that only one of the two candidate regions showed strong evidence of involvement in iNOS regulation. Specifically, we found evidence of 13 conserved binding site candidates linked to iNOS regulation: AP-1, C/EBPB, CREB, GATA-1, GATA-3, NF-AT, NF-AT5, NF-κB, KLF4, Oct-1, PEA3, SMAD3, and TCF11. Additionally, we found evidence of interspecies variation in binding sites for several regulatory elements linked to iNOS (GATA-3, GATA-4, KLF6, SRF, STAT-1, STAT-3, OLF-1 and HIF-1) across species, especially in African green monkeys relative to other species. Given the key role of iNOS in respiratory immune response, the findings of this study might help guide the direction of future studies aimed to uncover the molecular mechanisms underlying the increased susceptibility of African green monkeys to several viral and bacterial respiratory infections.

SNP-SIG 2013: the state of the art of genomic variant interpretation

Abstract Contact: yanab@rci.rutgers.edu or emidio@uab.edu .

A comparison of whole genome sequencing with exome sequencing for family-based association studies.

As the cost of DNA sequencing decreases, association studies based on whole genome sequencing are now becoming feasible. It is still unclear, however, how much more we could gain from whole genome sequencing compared to exome sequencing, which has been widely used to study a variety of diseases. In this project, we performed a comparison between whole genome sequencing and exome sequencing for family-based association analysis using data from Genetic Analysis Workshop 18. Whole genome sequencing was able to identify several significant hits within intergenic regions. However, the increased cost of multiple testing counteracted the benefits and resulted in a higher false discovery rate. Our results suggest that exome sequencing is a cost-effective way to identify disease-related variants. With the decreasing sequencing cost and accumulating knowledge of the human genome, whole genome sequencing has the potential to identify important variants in regulatory regions typically inaccessible for exome sequencing.

A Rb1 promoter variant with reduced activity contributes to osteosarcoma susceptibility in irradiated mice.

BackgroundSyndromic forms of osteosarcoma (OS) account for less than 10% of all recorded cases of this malignancy. An individual OS predisposition is also possible by the inheritance of low penetrance alleles of tumor susceptibility genes, usually without evidence of a syndromic condition. Genetic variants involved in such a non-syndromic form of tumor predisposition are difficult to identify, given the low incidence of osteosarcoma cases and the genetic heterogeneity of patients. We recently mapped a major OS susceptibility QTL to mouse chromosome 14 by comparing alpha-radiation induced osteosarcoma in mouse strains which differ in their tumor susceptibility.MethodsTumor-specific allelic losses in murine osteosacoma were mapped along chromosome 14 using microsatellite markers and SNP allelotyping. Candidate gene search in the mapped interval was refined using PosMed data mining and mRNA expression analysis in normal osteoblasts. A strain-specific promoter variant in Rb1 was tested for its influence on mRNA expression using reporter assay.ResultsA common Rb1 allele derived from the BALB/cHeNhg strain was identified as the major determinant of radiation-induced OS risk at this locus. Increased OS-risk is linked with a hexanucleotide deletion in the promoter region which is predicted to change WT1 and SP1 transcription factor-binding sites. Both in-vitro reporter and in-vivo expression assays confirmed an approx. 1.5 fold reduced gene expression by this promoter variant. Concordantly, the 50% reduction in Rb1 expression in mice bearing a conditional hemizygous Rb1 deletion causes a significant rise of OS incidence following alpha-irradiation.ConclusionThis is the first experimental demonstration of a functional and genetic link between reduced Rb1 expression from a common promoter variant and increased tumor risk after radiation exposure. We propose that a reduced Rb1 expression by common variants in regulatory regions can modify the risk for a malignant transformation of bone cells after radiation exposure.

Extensive variation in chromatin states across humans.

The majority of disease-associated variants lie outside protein-coding regions, suggesting a link between variation in regulatory regions and disease predisposition. We studied differences in chromatin states using five histone modifications, cohesin, and CTCF in lymphoblastoid lines from 19 individuals of diverse ancestry. We found extensive signal variation in regulatory regions, which often switch between active and repressed states across individuals. Enhancer activity is particularly diverse among individuals, whereas gene expression remains relatively stable. Chromatin variability shows genetic inheritance in trios, correlates with genetic variation and population divergence, and is associated with disruptions of transcription factor binding motifs. Overall, our results provide insights into chromatin variation among humans.

Sex-specific associations of variants in regulatory regions of NADPH oxidase-2 (CYBB) and glutathione peroxidase 4 (GPX4) genes with kidney disease in type 1 diabetes

Oxidative stress is involved in the pathophysiology of diabetic nephropathy. The superoxide-generating nicotinamide adenine dinucleotide phosphate-oxidase 2 (NOX2, encoded by the CYBB gene) and the antioxidant enzyme glutathione peroxidase 4 (GPX4) play opposing roles in the balance of cellular redox status. In the present study, we investigated associations of single nucleotide polymorphisms (SNPs) in the regulatory regions of CYBB and GPX4 with kidney disease in patients with type 1 diabetes. Two functional SNPs, rs6610650 (CYBB promoter region, chromosome X) and rs713041 (GPX4 3'untranslated region, chromosome 19), were genotyped in 451 patients with type 1 diabetes from a Brazilian cohort (diabetic nephropathy: 44.6%) and in 945 French/Belgian patients with type 1 diabetes from Genesis and GENEDIAB cohorts (diabetic nephropathy: 62.3%). The minor A-allele of CYBB rs6610650 was associated with lower estimated glomerular filtration rate (eGFR) in Brazilian women, and with the prevalence of established/advanced nephropathy in French/Belgian women (odds ratio 1.75, 95% CI 1.11–2.78, p = 0.016). The minor T-allele of GPX4 rs713041 was inversely associated with the prevalence of established/advanced nephropathy in Brazilian men (odds ratio 0.30, 95% CI 0.13–0.68, p = 0.004), and associated with higher eGFR in French/Belgian men. In conclusion, these heterogeneous results suggest that neither CYBB nor GPX4 are major genetic determinants of diabetic nephropathy, but nevertheless, they could modulate in a gender-specific manner the risk for renal disease in patients with type 1 diabetes.

Color differences among feral pigeons (Columba livia) are not attributable to sequence variation in the coding region of the melanocortin-1 receptor gene (MC1R)

BackgroundGenetic variation at the melanocortin-1 receptor (MC1R) gene is correlated with melanin color variation in many birds. Feral pigeons (Columba livia) show two major melanin-based colorations: a red coloration due to pheomelanic pigment and a black coloration due to eumelanic pigment. Furthermore, within each color type, feral pigeons display continuous variation in the amount of melanin pigment present in the feathers, with individuals varying from pure white to a full dark melanic color. Coloration is highly heritable and it has been suggested that it is under natural or sexual selection, or both. Our objective was to investigate whether MC1R allelic variants are associated with plumage color in feral pigeons.FindingsWe sequenced 888 bp of the coding sequence of MC1R among pigeons varying both in the type, eumelanin or pheomelanin, and the amount of melanin in their feathers. We detected 10 non-synonymous substitutions and 2 synonymous substitution but none of them were associated with a plumage type. It remains possible that non-synonymous substitutions that influence coloration are present in the short MC1R fragment that we did not sequence but this seems unlikely because we analyzed the entire functionally important region of the gene.ConclusionsOur results show that color differences among feral pigeons are probably not attributable to amino acid variation at the MC1R locus. Therefore, variation in regulatory regions of MC1R or variation in other genes may be responsible for the color polymorphism of feral pigeons.

Single nucleotide polymorphisms in the 3′ untranslated and near gene regions of transthyretin may play a role in senile systemic amyloidosis

Senile Systemic Amyloidosis (SSA) is a disease occurring mainly in elderly Caucasian males, estimated to affect approximately 25% of individuals over 80 years of age. SSA usually presents as cardiomyopathy caused by extracellular deposits of amyloid fibrils composed of the plasma protein, transthyretin (TTR). We hypothesized that genetic variation in regulatory regions of the TTR gene may contribute to SSA. We examined SNP variations in a region of ~800 bases spanning the 3′ untranslated and near gene regions of the TTR gene in a group of patients with SSA (n=88). For comparison, a group of controls (n=10) matched for age, gender, and race were analyzed. Of 24 previously reported SNPs in this region, we identified one novel and 3 common SNPs, rs62093482, rs1791228, and rs75032823 in our samples. The odds ratios (OR) calculated under allele analysis for these SNPs were &gt;;1, indicating at least one SNP may be associated with SSA. Genotype analysis revealed OR of 7 and 5 under the recessive model for rs62093482 and rs1791228, respectively, and OR of 3 under the dominant model for rs75032823. Haplotype analysis revealed a haplotype block unique to the SSA group including the three SNPs. The common haplotype identified had an OR of 2.5, indicative of association with SSA, and included alleles with high OR for each SNP.This research is supported by the Gruss Foundation, Boston University Amyloid Fund, the Young Family Amyloid Research Fund, and NIH grant RO1AG031804 (LHC).

The impact of 3'UTR variants on differential expression of candidate cancer susceptibility genes.

Variants in regulatory regions are predicted to play an important role in disease susceptibility of common diseases. Polymorphisms mapping to microRNA (miRNA) binding sites have been shown to disrupt the ability of miRNAs to target genes resulting in differential mRNA and protein expression. Skin tumor susceptibility 5 (Skts5) was identified as a locus conferring susceptibility to chemically-induced skin cancer in NIH/Ola by SPRET/Outbred F1 backcrosses. To determine if polymorphisms between the strains which mapped to putative miRNA binding sites in the 3′ untranslated region (3′UTR) of genes at Skts5 influenced expression, we conducted a systematic evaluation of 3′UTRs of candidate genes across this locus. Nine genes had polymorphisms in their 3′UTRs which fit the linkage data and eight of these contained polymorphisms suspected to interfere with or introduce miRNA binding. 3′UTRs of six genes, Bcap29, Dgkb, Hbp1, Pik3cg, Twistnb, and Tspan13 differentially affected luciferase expression, but did not appear to be differentially regulated by the evaluated miRNAs predicted to bind to only one of the two isoforms. 3′UTRs from four additional genes chosen from the locus that fit less stringent criteria were evaluated. Ifrd1 and Etv1 showed differences and contained polymorphisms predicted to disrupt or create miRNA binding sites but showed no difference in regulation by the miRNAs tested. In summary, multiple 3′UTRs with putative functional variants between susceptible and resistant strains of mice influenced differential expression independent of predicted miRNA binding.

Cis-regulatory sequence variation and association with Mycoplasma load in natural populations of the house finch (Carpodacus mexicanus).

Characterization of the genetic basis of fitness traits in natural populations is important for understanding how organisms adapt to the changing environment and to novel events, such as epizootics. However, candidate fitness-influencing loci, such as regulatory regions, are usually unavailable in nonmodel species. Here, we analyze sequence data from targeted resequencing of the cis-regulatory regions of three candidate genes for disease resistance (CD74, HSP90α, and LCP1) in populations of the house finch (Carpodacus mexicanus) historically exposed (Alabama) and naïve (Arizona) to Mycoplasma gallisepticum. Our study, the first to quantify variation in regulatory regions in wild birds, reveals that the upstream regions of CD74 and HSP90α are GC-rich, with the former exhibiting unusually low sequence variation for this species. We identified two SNPs, located in a GC-rich region immediately upstream of an inferred promoter site in the gene HSP90α, that were significantly associated with Mycoplasma pathogen load in the two populations. The SNPs are closely linked and situated in potential regulatory sequences: one in a binding site for the transcription factor nuclear NFYα and the other in a dinucleotide microsatellite ((GC)6). The genotype associated with pathogen load in the putative NFYα binding site was significantly overrepresented in the Alabama birds. However, we did not see strong effects of selection at this SNP, perhaps because selection has acted on standing genetic variation over an extremely short time in a highly recombining region. Our study is a useful starting point to explore functional relationships between sequence polymorphisms, gene expression, and phenotypic traits, such as pathogen resistance that affect fitness in the wild.

Whole-Genome Sequencing: Ready for Prime Time?

Rapid technological advances in DNA sequencing are facilitating the perception that whole-genome sequencing (WGS) will soon be available as a routine clinical diagnostic—and predictive—test. Current clinical testing is typically diagnostic, focusing on single genes or small gene panels, identifying a small number of variants per patient, and yet requiring nontrivial effort and expertise to deliver a clinical laboratory report. Interpreting WGS results will increase the necessary effort exponentially. How will this work in actual practice? In a recent Science Policy Forum (1), Radoje Drmanac promotes WGS as a research tool that will enhance disease diagnosis, prevention, and management by providing a more “comprehensive” method of genetic testing. WGS is touted as better for assessing disorders that have an apparent genetic etiology but that are not amenable to current methods, including those caused by noncoding variants, interaction between multiple variants across the genome, or variants in regulatory regions. Unfortunately, our current level of understanding precludes a clinical interpretation of most noncoding variants. Beyond …

DiGSNP: a web tool for Disease-Gene-SNP hierarchical prioritization

DiGSNP: a web tool for Disease-Gene-SNP hierarchical prioritization