Regulatory polymorphisms underlying complex disease traits

Ali ME, Abuduhier FM, Ullah MF. Clinical Implications of Krüpple-like Transcription Factor KLF-14 and Certain Micro-RNA (miR-27a, miR-196a2, miR-423) Gene Variations as a Risk Factor in the Genetic Predisposition to PCOS

Polycystic ovary syndrome (PCOS) is a disorder with a symptomatic manifestation of an array of metabolic and endocrine impairments. PCOS has a relatively high prevalence rate among young women of reproductive age and is a risk factor for some severe metabolic diseases such as T2DM, insulin insensitivity, and obesity, while the most dominant endocrine malfunction is an excess of testosterone showing hyperandrogenism and hirsutism. MicroRNAs have been implicated as mediators of metabolic diseases including obesity and insulin resistance, as these can regulate multiple cellular pathways such as insulin signaling and adipogenesis. Genome-wide association studies during the last few years have also linked the Krüpple-like family of transcription factors such as KLF14, which contribute in mechanisms of mammalian gene regulation, with certain altered metabolic traits and risk of atherosclerosis and type-2 DM. This study has characterized the biochemical and endocrine parameters in PCOS patients with a comprehensive serum profiling in comparison to healthy controls and further examined the influence of allelic variations for miRNAs 27a (rs895819 A > G), 196a2 (rs11614913 C > T), 423 (rs6505162C > A), and transcription factor KLF14 (rs972283 A > G) gene polymorphism on the risk and susceptibility to PCOS. The experimental protocol included amplification refractory mutation-specific (ARMS)-PCR to detect and determine the presence of these polymorphic variants in the study subjects. The results in this case−control study showed that most of the serum biomarkers, both biochemical and endocrine, that were analyzed in the study demonstrated statistically significant alterations in PCOS patients, including lipids (LDL, HDL, cholesterol), T2DM markers (fasting glucose, free insulin, HOMA-IR), and hormones (FSH, LH, testosterone, and progesterone). The distribution of Krüppel-like factor 14 rs972283 G > A, miR-27a rs895819 A > G, and miR-196a-2 rs11614913 C > T genotypes analyzed within PCOS patients and healthy controls in the considered population was significant (p < 0.05), except for miR-423 rs6505162 C > A genotypes (p > 0.05). The study found that in the codominant model, KLF14-AA was strongly associated with greater PCOS susceptibility (OR 2.35, 95% CI = 1.128 to 4.893, p < 0.022), miR-27a-GA was linked to an enhanced PCOS susceptibility (OR 2.06, 95% CI = 1.165 to 3.650, p < 0.012), and miR-196a-CT was associated with higher PCOS susceptibility (OR 2.06, 95% CI = 1.191 to 3.58, p < 0.009). Moreover, allele A of KLF-14 and allele T of miR-196a2 were strongly associated with PCOS susceptibility in the considered population.

Whole-genome sequencing reveals that variants in the Interleukin 18 Receptor Accessory Protein 3'UTR protect against ALS

The noncoding genome is substantially larger than the protein-coding genome but has been largely unexplored by genetic association studies. Here, we performed region-based rare variant association analysis of >25,000 variants in untranslated regions of 6,139 amyotrophic lateral sclerosis (ALS) whole genomes and the whole genomes of 70,403 non-ALS controls. We identified interleukin-18 receptor accessory protein (IL18RAP) 3' untranslated region (3'UTR) variants as significantly enriched in non-ALS genomes and associated with a fivefold reduced risk of developing ALS, and this was replicated in an independent cohort. These variants in the IL18RAP 3'UTR reduce mRNA stability and the binding of double-stranded RNA (dsRNA)-binding proteins. Finally, the variants of the IL18RAP 3'UTR confer a survival advantage for motor neurons because they dampen neurotoxicity of human induced pluripotent stem cell (iPSC)-derived microglia bearing an ALS-associated expansion in C9orf72, and this depends on NF-κB signaling. This study reveals genetic variants that protect against ALS by reducing neuroinflammation and emphasizes the importance of noncoding genetic association studies.

Preliminary Evidence That the Short Allele of 5-HTTLPR Moderates the Association of Psychiatric Symptom Severity on Suicide Attempt: The Example in Obsessive-Compulsive Disorder

BACKGROUND

The severity of symptoms represents an important source of distress in patients with a psychiatric disease. However, the extent to which this endogenous stress factor interacts with genetic vulnerability factors for predicting suicide risks remains unclear.

METHODS

We evaluated whether the severity of symptoms interacts with a genetic vulnerability factor (the serotonin transporter gene-linked promoter region variation) in predicting the frequency of lifetime suicide attempts in patients with a psychiatric disease. Symptom severity and 5-HTTLPR polymorphism were collected from a sample of 95 patients with obsessive-compulsive disorder (OCD). Lifetime suicide attempt was the primary outcome, and antecedent of multiple suicide attempts was the secondary outcome.

RESULTS

The gene-by-symptoms interaction was associated with an excess risk of suicide attempts (OR = 4.39, 95CI[1.44, 13.38], p < 0.009) and of multiple suicide attempts (OR = 4.18, 95CI[1.04, 16.77], p = 0.043). Symptom severity (moderate, severe, or extreme) was associated with an approximately five-fold increase in the odds of a lifetime suicide attempt in patients carrying one or two copies of the short allele of 5-HTTLPR. No such relationship was found for patients carrying the long allele.

CONCLUSION

This study provides preliminary evidence for the gene-by-stress interaction on suicide attempt when stress is operationalized as symptom severity. Progress in suicide research may come from efforts to investigate the gene-by-symptoms interaction hypothesis in a variety of diseases.

Assigning Co-Regulated Human Genes and Regulatory Gene Clusters

Elucidating the role of genetic variation in the regulation of gene expression is key to understanding the pathobiology of complex diseases which, in consequence, is crucial in devising targeted treatment options. Expression quantitative trait locus (eQTL) analysis correlates a genetic variant with the strength of gene expression, thus defining thousands of regulated genes in a multitude of human cell types and tissues. Some eQTL may not act independently of each other but instead may be regulated in a coordinated fashion by seemingly independent genetic variants. To address this issue, we combined the approaches of eQTL analysis and colocalization studies. Gene expression was determined in datasets comprising 49 tissues from the Genotype-Tissue Expression (GTEx) project. From about 33,000 regulated genes, over 14,000 were found to be co-regulated in pairs and were assembled across all tissues to almost 15,000 unique clusters containing up to nine regulated genes affected by the same eQTL signal. The distance of co-regulated eGenes was, on average, 112 kilobase pairs. Of 713 genes known to express clinical symptoms upon haploinsufficiency, 231 (32.4%) are part of at least one of the identified clusters. This calls for caution should treatment approaches aim at an upregulation of a haploinsufficient gene. In conclusion, we present an unbiased approach to identifying co-regulated genes in and across multiple tissues. Knowledge of such common effects is crucial to appreciate implications on biological pathways involved, specifically when a treatment option targets a co-regulated disease gene.

miR-155 T/A (rs767649) and miR-146a A/G (rs57095329) single nucleotide polymorphisms as risk factors for chronic hepatitis B virus infection among Egyptian patients

Genetic variants in microRNAs (miRNAs) can alter the miRNAs expression and/or function, accordingly, affecting the related biological pathways and disease risk. Dysregulation of miR-155 and miR-146a expression levels has been well-described in viral hepatitis B (HBV). In the current study, we aimed to assess rs767649 T/A and rs57095329 A/G polymorphisms in miR-155, and miR-146a genes, respectively, as risk factors for Chronic HBV (CHBV) in the Egyptian population. Also, we aimed to do in silico analysis to investigate the molecules that primarily target these miRNAs. One hundred patients diagnosed as CHBV and one hundred age and sex-matched controls with evidence of past HBV infection were genotyped for miR-155 (rs767649) and miR-146a (rs57095329) using real-time polymerase chain reaction. The rs767649 AT and AA genotypes in CHBV patients confer four folds and ten folds risk respectively, as compared to control subjects [(AOR = 4.245 (95%CI 2.009–8.970), p<0.0001) and AOR = 10.583 (95%CI 4.012–27.919), p<0.0001, respectively)]. The rs767649 A allele was associated with an increased risk of developing CHBV (AOR = 2.777 (95%CI 1.847–4.175), p<0.0001). There was a significant difference in the frequency of rs57095329 AG and GG genotypes in CHBV patients compared to controls. AG and GG genotypes showed an increase in the risk of developing CHBV by about three and six folds respectively [AOR = 2.610 (95%CI 1.362–5.000), p = 0.004] and [AOR = 5.604 (95%CI 2.157–14.563), p<0.0001].We concluded that rs57095329 and rs767649 SNPs can act as potential risk factors for the development of CHBV in the Egyptian population.

Single Nucleotide Polymorphism rs17173608 in the Chemerin Encoding Gene: Is It a Predictor of Insulin Resistance and Severity of Coronary Artery Disease in Non-Obese Type 2 Diabetes?

(1) Background: Chemerin, or the RARRES2 (Retinoic Acid Receptor Responder 2) gene, is found to be associated with an increased incidence of insulin resistance, endothelial dysfunction, type 2 diabetes (T2D), and coronary artery disease (CAD). This study investigates associations of RARRES2rs17173608 with insulin resistance and the severity of CAD in non-obese T2D patients in relation to the clinical and genetic factors. (2) Methods: A total of 300 patients with T2D and CAD were recruited in this study. The associations of insulin resistance and the severity of CAD with RARRES2rs17173608 and clinical factors were assessed. The genotyping procedures were performed using the TaqMan method. The significant associations (p ≤ 0.05) from preliminary tests were employed to carry out the secondary analysis. (3) Results: RARRES2rs17173608 (TT, TG, and GG polymorphisms in the preliminary analysis; TG and GG polymorphisms in a secondary analysis) was associated with insulin resistance and the severity of CAD in both the preliminary and secondary analysis (all p-values were < 0.05). Additionally, in the secondary analysis, FPG and ACEI were also associated with insulin resistance and the severity of CAD (all p-values were < 0.05). (4) Conclusion: From the preliminary findings, rs17173608 is a significant predictor of insulin resistance and the severity of CAD.

Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method

The highly polymorphic human major histocompatibility complex (MHC) also known as the human leukocyte antigen (HLA) encodes class I and II genes that are the cornerstone of the adaptive immune system. Their unique diversity (>25,000 alleles) might affect the outcome of any transplant, infection, and susceptibility to autoimmune diseases. The recent rapid development of new next-generation sequencing (NGS) methods provides the opportunity to study the influence/correlation of this high level of HLA diversity on allele expression levels in health and disease. Here, we describe the NGS capture RNA-Seq method that we developed for genotyping all 12 classical HLA loci (HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRA, HLA-DRB1, HLA-DRB3, HLA-DRB4, and HLA-DRB5) and assessing their allelic imbalance by quantifying their allele RNA levels. This is a target enrichment method where total RNA is converted to a sequencing-ready complementary DNA (cDNA) library and hybridized to a complex pool of RNA-specific HLA biotinylated oligonucleotide capture probes, prior to NGS. This method was applied to 161 peripheral blood mononuclear cells and 48 umbilical cord blood cells of healthy donors. The differential allelic expression of 10 HLA loci (except for HLA-DRA and HLA-DPA1) showed strong significant differences (P < 2.1 × 10-15). The results were corroborated by independent methods. This newly developed NGS method could be applied to a wide range of biological and medical questions including graft rejections and HLA-related diseases.

The role of ESCO2, SALL4 and TBX5 genes in the susceptibility to thalidomide teratogenesis

Thalidomide is widely used for several diseases; however, it causes malformations in embryos exposed during pregnancy. The complete understanding of the mechanisms by which thalidomide affects the embryo development has not yet been obtained. The phenotypic similarity makes TE a phenocopy of syndromes caused by mutations in ESCO2, SALL4 and TBX5 genes. Recently, SALL4 and TBX5 were demonstrated to be thalidomide targets. To understand if these genes act in the TE development, we sequenced them in 27 individuals with TE; we verified how thalidomide affect them in human pluripotent stem cells (hPSCs) through a differential gene expression (DGE) analysis from GSE63935; and we evaluated how these genes are functionally related through an interaction network analysis. We identified 8 variants in ESCO2, 15 in SALL4 and 15 in TBX5. We compared allelic frequencies with data from ExAC, 1000 Genomes and ABraOM databases; eight variants were significantly different (p < 0.05). Eleven variants in SALL4 and TBX5 were previously associated with cardiac diseases or malformations; however, in TE sample there was no association. Variant effect prediction tools showed 97% of the variants with potential to influence in these genes regulation. DGE analysis showed a significant reduction of ESCO2 in hPSCs after thalidomide exposure.

A panel of miRNAs is involved in the effect of sertraline on panic disorder, as implicated by a microarray-based analysis

BACKGROUND

MiRNAs are considered to be significant contributors to the pathogenesis of psychiatric diseases, but little is known about the potential roles of miRNAs in the treatment effect of panic disorder (PD). Therefore, we aimed to identify the miRNAs association with PD over the course of sertraline treatment.

METHODS

Sixty-seven patients were collected for a 6-week period of sertraline treatment, and evaluated using HAMD-17, HAMA-14 and PDSS both at the baseline and 6 weeks later. Blood samples were collected before and after treatment, respectively. Ten pairs of samples were analyzed using miRNA array, and the differentially expressed miRNAs were further validated using RT-PCR in the whole sample.

RESULTS

miR-451a, miR-144-5p, miR-25-3p and miR-660-5p were found to be significantly up-regulated, while miR-1 and miR-148-5p significantly down-regulated after sertraline treatment. The change of miR-25-3p before and after treatment (△miR-25-3p) was positively related to both the changes of PDSS3 scores (△PDSS3) (p = 0.017, 31.5% contribution) and △ PDSS7 (p = 0.016, 32.3% contribution). The △miR-660-5p was positively related to both the △HAMA5 (p = 0.03, 26% contribution) and △PDSS7 (p = 0.032). The △miR-148-5p was positively related to the △PDSS4 (p = 0.046, 21.5% contribution), but negatively related to the △HAMA13 (p = 0.005, 41.9% contribution). The △miR-144-5p was negatively related to the △HAMA9 (p = 0.032, 25.3% contribution).

CONCLUSIONS

These findings might provide some evidences to the involvement of miRNA in the effect of anti-anxiety agents, which contributed to the better understanding the disease and developing new therapeutic genetic targets.

Detection of cooperatively bound transcription factor pairs using ChIP-seq peak intensities and expectation maximization

Transcription factors (TFs) often work cooperatively, where the binding of one TF to DNA enhances the binding affinity of a second TF to a nearby location. Such cooperative binding is important for activating gene expression from promoters and enhancers in both prokaryotic and eukaryotic cells. Existing methods to detect cooperative binding of a TF pair rely on analyzing the sequence that is bound. We propose a method that uses, instead, only ChIP-seq peak intensities and an expectation maximization (CPI-EM) algorithm. We validate our method using ChIP-seq data from cells where one of a pair of TFs under consideration has been genetically knocked out. Our algorithm relies on our observation that cooperative TF-TF binding is correlated with weak binding of one of the TFs, which we demonstrate in a variety of cell types, including E. coli, S. cerevisiae and M. musculus cells. We show that this method performs significantly better than a predictor based only on the ChIP-seq peak distance of the TFs under consideration. This suggests that peak intensities contain information that can help detect the cooperative binding of a TF pair. CPI-EM also outperforms an existing sequence-based algorithm in detecting cooperative binding. The CPI-EM algorithm is available at https://github.com/vishakad/cpi-em.

Genetic variations in LTA gene and PDCD1 gene and intrauterine infection of hepatitis B virus: a case-control study in China

Intrauterine infection with hepatitis B virus (HBV) has been suggested to accounting for most cases of chronic HBV infection, which cannot be blocked by combined immunoprophylaxis. The fact that the genetic background might impact the susceptibility to intrauterine infection of HBV has been identified by recent researches. A case-control study included sixty-nine HBsAg-positive mother-newborn pairs with intrauterine infection as cases compared to 138 mother-newborn pairs without intrauterine infection as controls. We studied the correlations between HBV intrauterine transmission and 15 maternal SNPs in eight genes (LTA, LTBR, TNFSF14, PDCD1, APOBEC3B, CD274, CD40 and CD40LG). There was a substantially significantly decreased risk of intrauterine infection of HBV in mothers with the rs2227981 TT genotype in PDCD1 gene compared to those with the rs2227981 GG genotype (OR 0.11, 95% CI 0.01-0.95, P = 0.045). Under recessive model (OR 0.51, 95% CI 0.26-1, P = 0.050) and additive model (OR 0.50, 95% CI 0.28-0.88, P = 0.017), we also found a marginally significantly decreased risk of intrauterine infection of HBV. Furthermore, under additive model, maternal genotype for rs2239704 in LTA gene was marginally significantly related to an increased risk of intrauterine HBV infection (OR 1.62, 95% CI 1-6.66, P = 0.055). However, there were no statistically significant associations among the remaining 13 SNPs and the risk of intrauterine infection of HBV. The examination implied that hereditary variants of PDCD1 and LTA genes were associated with intrauterine infection of HBV.

Genome-Wide Association Studies: Progress in Identifying Genetic Biomarkers in Common, Complex Diseases

Novel, comprehensive approaches for biomarker discovery and validation are urgently needed. One particular area of methodologic need is for discovery of novel genetic biomarkers in complex diseases and traits. Here, we review recent successes in the use of genome wide association (GWA) approaches to identify genetic biomarkers in common human diseases and traits. Such studies are yielding initial insights into the allelic architecture of complex traits. In general, it appears that complex diseases are associated with many common polymorphisms, implying profound genetic heterogeneity between affected individuals.

Alterations in the expression of a neurodevelopmental gene exert long-lasting effects on cognitive-emotional phenotypes and functional brain networks: translational evidence from the stress-resilient Ahi1 knockout mouse

Many psychiatric disorders are highly heritable and may represent the clinical outcome of early aberrations in the formation of neural networks. The placement of brain connectivity as an 'intermediate phenotype' renders it an attractive target for exploring its interaction with genomics and behavior. Given the complexity of genetic make up and phenotypic heterogeneity in humans, translational studies are indicated. Recently, we demonstrated that a mouse model with heterozygous knockout of the key neurodevelopmental gene Ahi1 displays a consistent stress-resilient phenotype. Extending these data, the current research describes our multi-faceted effort to link early variations in Ahi1 expression with long-term consequences for functional brain networks and cognitive-emotional phenotypes. By combining behavioral paradigms with graph-based analysis of whole-brain functional networks, and then cross-validating the data with robust neuroinformatic data sets, our research suggests that physiological variation in gene expression during neurodevelopment is eventually translated into a continuum of global network metrics that serve as intermediate phenotypes. Within this framework, we suggest that organization of functional brain networks may result, in part, from an adaptive trade-off between efficiency and resilience, ultimately culminating in a phenotypic diversity that encompasses dimensions such as emotional regulation and cognitive function.

NFATC1 genotypes affect acute rejection and long-term graft function in cyclosporine-treated renal transplant recipients

AIM

To investigate the effects of SNPs in the cyclophilin A/calcineurin/nuclear factor of activated T-cells (NFATs) pathway genes (PPIA, PPP3CB, PPP3R1, NFATC1 and NFATC2) on cyclosporine (CsA) efficacy in renal transplant recipients.

MATERIALS & METHODS

Seventy-six tag SNPs were detected in 155 CsA-treated renal recipients with at least a 5-year follow-up. The associations of SNPs with acute rejection, nephrotoxicity, pneumonia and estimated glomerular filtration rate post transplant were explored.

RESULTS

NFATC1 rs3894049 GC was a risk factor for acute rejection compared with CC carriers (p = 0.0005). NFATC1 rs2280055 TT carriers had a more stable estimated glomerular filtration rate level than CC (p = 0.0004).

CONCLUSION

Detecting NFATC1 polymorphisms could help predict CsA efficacy in renal transplant patients.

Genetic variation in leptin and leptin receptor genes is a risk factor for idiopathic recurrent spontaneous abortion

AIM

To determine whether maternal leptin (LEP) and leptin receptor (LEPR) gene polymorphisms are associated with idiopathic recurrent spontaneous abortion (IRSA).

METHODS

This case-control association study conducted from 2010 to 2012 at the Department of Gynecology and Obstetrics, University Hospital Center Osijek and Clinical Institute of Medical Genetics Ljubljana included 178 women with a history of three or more IRSAs before the 22nd week of gestation and 145 women with at least two live births and no history of pathologic pregnancies during reproductive period. Polymorphisms of maternal LEP (rs7799039, rs2122627, rs11761556, rs10244329) and LEPR (rs1137101, rs7516341, rs1186403, rs12062820) were assessed by allele specific real-time polymerase chain reaction. Genotype distribution, allele frequencies, and frequency of haplotypes at LEP and LEPR genetic loci were determined.

RESULTS

We observed more frequent genotype for rs7516341 (nominal P=0.034, odds ratio [OR] 0.61, 95% confidence interval [CI] 0.38-0.97) and rs1137101 (nominal P=0.048, OR 1.66, 95% CI 1.00-2.80) in the LEPR gene in patients than in controls, but these results did not remain significant after correction for multiple testing according to Bonferroni (adjusted P value threshold was set at 0.05). We did not observe differential distribution of genotype frequencies in the LEP gene between cases and controls. In patients with IRSA, GTCC haplotype in the LEPR gene locus was significantly less frequent than in controls (PP=0.00865, OR 0.45), contrary to ACTC haplotype (PP=0.0087, OR 1.98).

CONCLUSIONS

We showed that genetic variability in the LEPR gene was associated with IRSA, warranting confirmation on a greater number of patients and pathogenesis investigation.

Identifying changes in punitive transcriptional factor binding sites from regulatory single nucleotide polymorphisms that are significantly associated with disease or sickness

AIM

To identify punitive transcriptional factor binding sites (TFBS) from regulatory single nucleotide polymorphisms (rSNPs) that are significantly associated with disease.

METHODS

The genome-wide association studies have provided us with nearly 6500 disease or trait-predisposing SNPs where 93% are located within non-coding regions such as gene regulatory or intergenic areas of the genome. In the regulatory region of a gene, a SNP can change the DNA sequence of a transcriptional factor (TF) motif and in turn may affect the process of gene regulation. SNP changes that affect gene expression and impact gene regulatory sequences such as promoters, enhancers, and silencers are known as rSNPs. Computational tools can be used to identify unique punitive TFBS created by rSNPs that are associated with disease or sickness. Computational analysis was used to identify punitive TFBS generated by the alleles of these rSNPs.

RESULTS

rSNPs within nine genes that have been significantly associated with disease or sickness were used to illustrate the tremendous diversity of punitive unique TFBS that can be generated by their alleles. The genes studied are the adrenergic, beta, receptor kinase 1, the v-akt murine thymoma viral oncogene homolog 3, the activating transcription factor 3, the type 2 demodkinase gene, the endothetal Per-Arnt-Sim domain protein 1, the lysosomal acid lipase A, the signal Transducer and Activator of Transcription 4, the thromboxane A2 receptor and the vascular endothelial growth factor A. From this sampling of SNPs among the nine genes, there are 73 potential unique TFBS generated by the common alleles compared to 124 generated by the minor alleles indicating the tremendous diversity of potential TFs that are capable of regulating these genes.

CONCLUSION

From the diversity of unique punitive binding sites for TFs, it was found that some TFs play a role in the disease or sickness being studied.

Possible role of intronic polymorphisms in the PHACTR1 gene on the development of cardiovascular disease

Cardiovascular disease (CVD) is a complex multifactorial and polygenetic disease in which the interaction of numerous genes, genetic variants, and environmental factors plays a major role in its development. In an attempt to demonstrate the association between certain genetic variants and CVD, researchers have run large genomic wild association studies (GWAS) in recent decades. These studies have correlated several genomic variants with the presence of CVD. Recently, certain polymorphisms in the phosphatase and actin regulator 1 (PHACTR1) gene have been shown to be associated with CVD (i.e., coronary artery disease, coronary artery calcification, early onset myocardial infarction, cervical artery dissection and hypertension) in different ethnic groups. It is important to state that all of the described PHACTR1 genetic variants associated with CVD are located in non-translating gene regions known as introns. Thus, the purpose of this article is to hypothesize the effect of certain intronic polymorphisms in the PHACTR1 gene on pathological processes in the cardiovascular system. In addition, we present compelling evidence that supports this hypothesis as well as a methodology that could be used to assess the allelic effect using in vitro and in vivo models, which will ultimately demonstrate the pathophysiological contribution of PHACTR1 intronic polymorphisms to the development of CVD.

The genetic architecture of autism spectrum disorders (ASDs) and the potential importance of common regulatory genetic variants

Currently, there is great interest in identifying genetic variants that contribute to the risk of developing autism spectrum disorders (ASDs), due in part to recent increases in the frequency of diagnosis of these disorders worldwide. While there is nearly universal agreement that ASDs are complex diseases, with multiple genetic and environmental contributing factors, there is less agreement concerning the relative importance of common vs rare genetic variants in ASD liability. Recent observations that rare mutations and copy number variants (CNVs) are frequently associated with ASDs, combined with reduced fecundity of individuals with these disorders, has led to the hypothesis that ASDs are caused primarily by de novo or rare genetic mutations. Based on this model, large-scale whole-genome DNA sequencing has been proposed as the most appropriate method for discovering ASD liability genes. While this approach will undoubtedly identify many novel candidate genes and produce important new insights concerning the genetic causes of these disorders, a full accounting of the genetics of ASDs will be incomplete absent an understanding of the contributions of common regulatory variants, which are likely to influence ASD liability by modifying the effects of rare variants or, by assuming unfavorable combinations, directly produce these disorders. Because it is not yet possible to identify regulatory genetic variants by examination of DNA sequences alone, their identification will require experimentation. In this essay, I discuss these issues and describe the advantages of measurements of allelic expression imbalance (AEI) of mRNA expression for identifying cis-acting regulatory variants that contribute to ASDs.

Correlation between the level of microRNA expression in peripheral blood mononuclear cells and symptomatology in patients with generalized anxiety disorder

This study investigated the correlation between the level of microRNA expression in peripheral blood mononuclear cells (PBMCs) and symptomatology in patients with generalized anxiety disorder (GAD). MicroRNA array was performed in peripheral blood mononuclear cells (PBMCs) obtained from GAD patients with gender, age, ethnicity-matched healthy controls. Then real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to verify the top 7 miRNAs with the highest fold-change values in 76 GAD patients and 39 healthy controls. It demonstrated that 5 miRNAs showed significantly differences in expression levels (P<0.01). These 5 GAD-associated miRNAs were finally selected into our study to analyze the association between the plasma level of miRNAs expression and symptomatology scores in Hamilton Anxiety Scale (HAMA). Results showed that the level of miR-4505 and miR-663 was negatively correlated with the total HAMA scores in GAD patients (r=0.2228, r=0.264 P<0.05). MiR-663 was selected into the regression equation of HAMA total scores and psychic anxiety symptomatology scores, and it could explain 5.3% of the HAMA total scores and 15.3% of the anxiety symptomatology scores. This study analyzed preliminarily possible circulating miRNAs expression changes in GAD patients, and the expression level of miR-663 highly correlated with psychic anxiety symptoms, further molecular mechanism of which needs to be explored.

The handedness-associated PCSK6 locus spans an intronic promoter regulating novel transcripts

We recently reported the association of the PCSK6 gene with handedness through a quantitative genome-wide association study (GWAS; P < 0.5 × 10⁻⁸) for a relative hand skill measure in individuals with dyslexia. PCSK6 activates Nodal, a morphogen involved in regulating left–right body axis determination. Therefore, the GWAS data suggest that the biology underlying the patterning of structural asymmetries may also contribute to behavioural laterality, e.g. handedness. The association is further supported by an independent study reporting a variable number tandem repeat (VNTR) within the same PCSK6 locus to be associated with degree of handedness in a general population cohort. Here, we have conducted a functional analysis of the PCSK6 locus combining further genetic analysis, in silico predictions and molecular assays. We have shown that the previous GWAS signal was not tagging a VNTR effect, suggesting that the two markers have independent effects. We demonstrated experimentally that one of the top GWAS-associated markers, rs11855145, directly alters the binding site for a nuclear factor. Furthermore, we have shown that the predicted regulatory region adjacent to rs11855415 acts as a bidirectional promoter controlling the expression of novel RNA transcripts. These include both an antisense long non-coding RNA (lncRNA) and a short PCSK6 isoform predicted to be coding. This is the first molecular characterization of a handedness-associated locus that supports the role of common variants in non-coding sequences in influencing complex phenotypes through gene expression regulation.

Regulatory SNPs Alter the Gene Expression of Diabetic Retinopathy Associated Secretary Factors

OBJECTIVES

Diabetic retinopathy (DR) is a common microvascular complication in both type I and type II diabetes. Several previous reports indicated the serum centration of some secretary factors were highly associated with DR. Therefore, we hypothesis regulatory SNPs (rSNPs) genotype in secretary factors may alter these gene expression and lead to DR.

METHODS

At first, pyrosequencing were applying to screen the SNPs which present allele frequency different in DR and DNR. Then individual genotyping was processed by Taqman assays in Taiwanese DR and DNR patients. To evaluate the effect of SNP allele on transcriptional activity, we measured promoter activity using luciferase reporter constructs.

RESULTS

We found the frequencies of the CC, CG, and GG genotype of the rs2010963 polymorphism were 15.09%, 47.14%, and 37.74% in DR and 12.90%, 19.35%, and 67.74% in DNR, respectively (p = 0.0205). The prevalence of DR was higher (p = 0.00793) in patients with the CC or CG genotype (62.26% and 32.26% for DR and DNR, respectively) compared with the patients with the GG genotype. To evaluate the effect of rs2010963-C allele on transcriptional activity, we measured promoter activity using luciferase reporter constructs. The rs2010963-C reporter showed 1.6 to 2-fold higher luciferase activity than rs2010963-G in 3 cell lines.

CONCLUSION

Our data proposed rs2010963-C altered the expression level of VEGFA in different tissues. We suggested small increase but long term exposure to VEGFA may lead to DR finally.

Genome-wide analysis of the genetic regulation of gene expression in human neutrophils

Neutrophils are an abundant immune cell type involved in both antimicrobial defence and autoimmunity. The regulation of their gene expression, however, is still largely unknown. Here we report an eQTL study on isolated neutrophils from 114 healthy individuals of Chinese ethnicity, identifying 21,210 eQTLs on 832 unique genes. Unsupervised clustering analysis of these eQTLs confirms their role in inflammatory responses and immunological diseases but also indicates strong involvement in dermatological pathologies. One of the strongest eQTL identified (rs2058660) is also the tagSNP of a linkage block reported to affect leprosy and Crohn's disease in opposite directions. In a functional study, we can link the C allele with low expression of the β-chain of IL18-receptor (IL18RAP). In neutrophils, this results in a reduced responsiveness to IL-18, detected both on the RNA and protein level. Thus, the polymorphic regulation of human neutrophils can impact beneficial as well as pathological inflammatory responses.

Neutrophils are abundant immune cells important for antimicrobial defence and in autoimmunity. Here, by mapping expression quantitative trait loci (eQTL) in neutrophils of Chinese ethnicity from Singapore, Andiappan et al. provide a resource for understanding immune-related trait associated genetic variants.

RNA-Seq alignment to individualized genomes improves transcript abundance estimates in multiparent populations

Massively parallel RNA sequencing (RNA-seq) has yielded a wealth of new insights into transcriptional regulation. A first step in the analysis of RNA-seq data is the alignment of short sequence reads to a common reference genome or transcriptome. Genetic variants that distinguish individual genomes from the reference sequence can cause reads to be misaligned, resulting in biased estimates of transcript abundance. Fine-tuning of read alignment algorithms does not correct this problem. We have developed Seqnature software to construct individualized diploid genomes and transcriptomes for multiparent populations and have implemented a complete analysis pipeline that incorporates other existing software tools. We demonstrate in simulated and real data sets that alignment to individualized transcriptomes increases read mapping accuracy, improves estimation of transcript abundance, and enables the direct estimation of allele-specific expression. Moreover, when applied to expression QTL mapping we find that our individualized alignment strategy corrects false-positive linkage signals and unmasks hidden associations. We recommend the use of individualized diploid genomes over reference sequence alignment for all applications of high-throughput sequencing technology in genetically diverse populations.

PTK2 rs7460 and rs7843014 polymorphisms and exceptional longevity: a functional replication study

Focal adhesion is critical for cell survival. The focal adhesion kinase (FAK, or PTK2) is an important component of the human interactome and thus is a potential longevity-related protein. Here we studied the association between two PTK2 gene single-nucleotide polymorphisms (SNPs) (rs7843014, rs7460) and exceptional longevity (EL). In addition to gaining insight into their functionality by determining luciferase gene reporter activity, we studied the genotype/allele frequency of these two SNPs among three different cohorts: (1) Spanish centenarians (n=175, 100-111 years, 144 women) and healthy controls (n=355, 20-50 years, 284 women); (2) Italian centenarians (n=79, 100-104 years, 40 women) and controls (n=316, 29-50 years, 156 women); and (3) Japanese centenarians (n=742, 100-116 years, 623 women) and healthy controls (n=499, 23-59 years, 356 women). Both SNPs had functional significance, with the A allele up-regulating luciferase activity compared to the other allele (rs7460 T allele and rs7843014 C allele, respectively). The A allele of both SNPs was negatively associated with EL in the Spanish cohort (rs7460, odds ratio [OR] adjusted by sex=0.40, 95% confidence intervals [CI] 0.3, 0.6, p<0.001); rs7843014, OR=0.37, 95% CI 0.3, 0.5, p<0.001). The OR of being a centenarian if having the rs7460-TT genotype was 6.68 (95% CI 4.1, 10.8, p<0.001). The rs7843014 CC genotype was also positively associated with EL (OR=7.58, 95% CI 4.6, 12.3, p<0.001]. No association was, however, found for the Italian or Japanese cohorts. Thus, two genotypes of the FAK gene, rs7460 TT and rs7843014 CC, are possibly associated with lower gene expression and might favor the likelihood of reaching EL in the Spanish population. Further research is needed to unveil the mechanisms by which FAK expression could perhaps influence the rate of aging.

Exceptional longevity and muscle and fitness related genotypes: a functional in vitro analysis and case-control association replication study with SNPs THRH rs7832552, IL6 rs1800795, and ACSL1 rs6552828

There are several gene variants that are candidates to influence functional capacity in long-lived individuals. As such, their potential association with exceptional longevity (EL, i.e., reaching 100+ years) deserves analysis. Among them are rs7832552 in the thyrotropin-releasing hormone receptor (TRHR) gene, rs1800795 in the interleukin-6 (IL6) gene and rs6552828 in the coenzyme A synthetase long-chain 1 (ACSL1) gene. To gain insight into their functionality (which is yet unknown), here we determined for the first time luciferase gene reporter activity at the muscle tissue level in rs7832552 and rs6552828. We then compared allele/genotype frequencies of the 3 abovementioned variants among centenarians [n = 138, age range 100–111 years (114 women)] and healthy controls [n = 334, 20–50 years (141 women)] of the same ethnic and geographic origin (Spain). We also studied healthy centenarians [n = 79, 100–104 years (40 women)] and controls [n = 316, 27–81 years (156 women)] from Italy, and centenarians [n = 742, 100–116 years (623 women)] and healthy controls [n = 499, 23–59 years (356 women)] from Japan. The THRH rs7832552 T-allele and ACSL1 rs6552828 A-allele up-regulated luciferase activity compared to the C and G-allele, respectively (P = 0.001). Yet we found no significant association of EL with rs7832552, rs1800795 or rs6552828 in any of the 3 cohorts. Further research is needed with larger cohorts of centenarians of different origin as well as with younger old people.

Polymorphic rs9471643 and rs6458238 upregulate PGC transcription and protein expression in overdominant or dominant models

The pepsinogen C (PGC) gene encodes a major differentiation biomarker for gastric mucosa and has two single nucleotide polymorphisms, rs9471643 G>C and rs6458238 G>A, within its 5' upstream region that are involved in gastric carcinogenesis. However, in what genetic models the two polymorphisms modulate disease risk and how they relate to gastric carcinogenesis needs further study. We fitted the most appropriate genetic models to the PGC polymorphisms and validated their robustness; then with knowledge of the genetic model, we investigated the influence of functional variant alleles or genotypes on gene expression in vitro and in vivo. We confirmed that rs9471643 CG genotype was stably associated with reduced gastric cancer risk in complete overdominant model. This favorable CG genotype was also associated with reduced atrophic gastritis risk in subjects carrying rs6458238 AG/AA genotype. The G>C transition at rs9471643 enhanced promoter activity and transcription factor binding ability, and the CG genotype was consistently associated with elevated levels of PGC mRNA, in situ protein and serum protein in complete overdominant model based-analyses. Additionally, rs6458238 AG/AA genotype was associated with reduced atrophic gastritis risk in dominant model. Its favorable A allele was related to higher promoter activity and lower transcription factor binding ability, and the AG/AA genotype showed association with elevated levels of serum PGC protein in dominant model based-analyses. Our results suggest that rs9471643 CG and rs6458238 AG/AA genotypes have important roles in up-regulating PGC expression, which may partially explain why individuals with these favorable genotypes have decreased risks of getting gastric cancer.

PGC-related gene variants and elite endurance athletic status in a Chinese cohort: a functional study

This study aims to examine the association between proliferator-activated receptor γ (PGC)-gene family-related single nucleotide polymorphisms (SNPs) and elite endurance runners' status in a Chinese cohort, and to gain insights into the functionality of a subset of SNPs. Genotype distributions of 133 SNPs in PPARGC1A, PPARGC1B, PPRC1, TFAM, TFB1M, TFB2M, NRF1, GABPA, GABPB1, ERRα, and SIRT1 genes were compared between 235 elite Chinese (Han) endurance runners (127 women) and 504 healthy non-athletic controls (237 women). Luciferase gene reporter activity was determined in 20 SNPs. After adjusting for multiple comparisons (in which threshold P-value was set at 0.00041), no significant differences were found in allele/genotype frequencies between athletes and controls (when both sexes were analyzed either together or separately). The lowest P-value was found in PPARGC1A rs4697425 (P = 0.001 for the comparison of allele frequencies between elite female endurance runners and their gender-matched controls). However, no association (all P > 0.05) was observed for this SNP in a replication cohort from Poland (194 endurance athletes and 190 controls). Using functional genomics tool, the following SNPs were found to have functional significance: PPARGC1A rs6821591, rs12650562, rs12374310, rs4697425, rs13113110, and rs4452416; PPARGC1B rs251466 and rs17110586; and PPRC1 rs17114388 (all P < 0.001). This study found no significant association between PGC-related SNPs and elite endurance athlete status in the Chinese population, despite some SNPs showing potential functional significance and the strong biological rationale to hypothesize that this gene pathway is a candidate to influence endurance exercise capacity.

FNDC5 (irisin) gene and exceptional longevity: a functional replication study with rs16835198 and rs726344 SNPs

Irisin might play an important role in reducing the risk of obesity, insulin resistance, or several related diseases, and high irisin levels may contribute to successful aging. Thus, the irisin precursor (FNDC5) gene is a candidate to influence exceptional longevity (EL), i.e., being a centenarian. It has been recently shown that two single-nucleotide polymorphisms (SNPs) in the FNDC5 gene, rs16835198 and rs726344, are associated with in vivo insulin sensitivity in adults. We determined luciferase gene reporter activity in the two above-mentioned SNPs and studied genotype distributions among centenarians (n = 175, 144 women) and healthy controls (n = 347, 142 women) from Spain. We also studied an Italian [79 healthy centenarians (40 women) and 316 healthy controls (156 women)] and a Japanese cohort [742 centenarians (623 women) and 499 healthy controls (356 women)]. The rs726344 SNP had functional significance, as shown by differences in luciferase activity between the constructs of this SNP (all P ≤ 0.05), with the variant A-allele having higher luciferase activity compared with the G-allele (P = 0.04). For the rs16835198 SNP, the variant T-allele tended to show higher luciferase activity compared with the G-allele (P = 0.07). However, we found no differences between genotype/allele frequencies of the two SNPs in centenarians versus controls in any cohort, and no significant association (using logistic regression adjusted by sex) between the two SNPs and EL. Further research is needed with different cohorts as well as with additional variants in the FNDC5 gene or in other genes involved in irisin signaling.

Association of polymorphism in microRNA 604 with susceptibility to persistent hepatitis B virus infection and development of hepatocellular carcinoma

MicroRNA polymorphisms may be associated with carcinogenesis or immunopathogenesis of infection. We evaluated whether the mircoRNA-604 (miR-604) polymorphism can affect the persistence of hepatitis B virus (HBV) infection, and the development to hepatocellular carcinoma (HCC) in patients with chronic HBV infection. A total of 1,439 subjects, who have either past or present HBV infection, were enrolled and divided into four groups (spontaneous recovery, chronic HBV carrier without cirrhosis, liver cirrhosis and HCC). We genotyped the precursor miR-604 genome region polymorphism. The CC genotype of miR-604 rs2368392 was most frequently observed and T allele frequency was 0.326 in all study subjects. The HBV persistence after infection was higher in those subjects with miR-604 T allele (P=0.05 in a co-dominant and dominant model), which implied that the patients with miR-604 T allele may have a higher risk for HBV chronicity. In contrast, there was a higher rate of the miR-604 T allele in the chronic carrier without HCC patients, compared to those of the HCC patients (P=0.03 in a co-dominant model, P=0.02 in a recessive model). The T allele at miR-604 rs2368392 may be a risk allele for the chronicity of HBV infection, but may be a protective allele for the progression to HCC in chronic HBV carriers.

Effects of genetic variations on microRNA: target interactions

Genetic variations within microRNA (miRNA) binding sites can affect miRNA-mediated gene regulation, which may lead to phenotypes and diseases. We perform a transcriptome-scale analysis of genetic variants and miRNA:target interactions identified by CLASH. This analysis reveals that rare variants tend to reside in CDSs, whereas common variants tend to reside in the 3′ UTRs. miRNA binding sites are more likely to reside within those targets in the transcriptome with lower variant densities, especially target regions in which nucleotides have low mutation frequencies. Furthermore, an overwhelming majority of genetic variants within or near miRNA binding sites can alter not only the potential of miRNA:target hybridization but also the structural accessibility of the binding sites and flanking regions. These suggest an interpretation for certain associations between genetic variants and diseases, i.e. modulation of miRNA-mediated gene regulation by common or rare variants within or near miRNA binding sites, likely through target structure alterations. Our data will be valuable for discovering new associations among miRNAs, genetic variations and human diseases.

Neural mechanisms underlying stress resilience in Ahi1 knockout mice: relevance to neuropsychiatric disorders

The Abelson helper integration site 1 (AHI1) gene has a pivotal role in brain development. Studies by our group and others have demonstrated association of AHI1 with schizophrenia and autism. To elucidate the mechanism whereby alteration in AHI1 expression may be implicated in the pathogenesis of neuropsychiatric disorders, we studied Ahi1 heterozygous knockout (Ahi1(+/-)) mice. Although their performance was not different from wild-type mice on tests that model classical schizophrenia-related endophenotypes, Ahi1(+/-) mice displayed an anxiolytic-like phenotype across different converging modalities. Using behavioral paradigms that involve exposure to environmental and social stress, significantly decreased anxiety was evident in the open field, elevated plus maze and dark-light box, as well as during social interaction in pairs. Assessment of core temperature and corticosterone secretion revealed a significantly blunted response of the autonomic nervous system and the hypothalamic-pituitary-adrenal axis in Ahi1(+/-) mice exposed to environmental and visceral stress. However, response to centrally acting anxiogenic compounds was intact. On resting-state functional MRI, connectivity of the amygdala with other brain regions involved in processing of anxiogenic stimuli and inhibitory avoidance learning, such as the lateral entorhinal cortex, ventral hippocampus and ventral tegmental area, was significantly reduced in the mutant mice. Taken together, our data link Ahi1 under-expression with a defect in the process of threat detection. Alternatively, the results could be interpreted as representing an anxiety-related endophenotype, possibly granting the Ahi1(+/-) mouse relative resilience to various types of stress. The current knockout model highlights the contribution of translational approaches to understanding the genetic basis of emotional regulation and its associated neurocircuitry, with possible relevance to neuropsychiatric disorders.

The maternal ITPK1 gene polymorphism is associated with neural tube defects in a high-risk Chinese population

Background

Epidemiological surveys and animal studies have revealed that inositol metabolism is associated with NTDs, but the mechanisms are not clear. Inositol 1,3,4-trisphosphate 5/6-kinase (ITPK1) is a pivotal regulatory enzyme in inositol metabolic pathway. The objective was to assess the potential impact of the maternal ITPK1 genotypes on the inositol parameter and on the NTD risk in a NTD high-risk area in China.

Methodology/Results

A case-control study of pregnant women affected with NTDs (n = 200) and controls (n = 320) was carried out. 13 tag SNPs of ITPK1 were selected and genotyped by the Sequenom MassArray system. We found that 4 tag SNPs were statistically significant in spina bifida group (P<0.05). MACH was used to impute the un-genotyped SNPs in ITPK1 locus and showed that 3 meaningful SNPs in the non-coding regions were significant. We also predicted the binding capacity of transcription factors in the positive SNPs using the bioinformatics method and found that only rs3783903 was located in the conserved sequence of activator protein-1 (AP-1). To further study the association between biochemical values and genotypes, maternal plasma inositol hexakisphosphate (IP₆) levels were also assessed using LC-MS. The maternal plasma IP₆ concentrations in the spina bifida subgroup were 7.1% lower than control (136.67 vs. 147.05 ng mL⁻¹, P<0.05), and significantly lower in rs3783903 GG genotype than others (P<0.05). EMSA showed a different allelic binding capacity of AP-1 in rs3783903, which was affected by an A→G exchange. The RT-PCR suggested the ITPK1 expression was decreased significantly in mutant-type of rs3783903 compared with wild-type in the 60 healthy pregnancies (P<0.05).

Conclusions/Significance

These results suggested that the maternal rs3783903 of ITPK1 might be associated with spina bifida, and the allele G of rs3783903 might affect the binding of AP-1 and the decrease of maternal plasma IP₆ concentration in this Chinese population.

VEGFA rSNPs, transcriptional factor binding sites and human disease

Three regulatory SNPs (rSNPs) in the promoter region of the vascular endothelial growth factor-A (VEGFA) gene have been significantly associated with several human diseases or conditions. The rSNP alleles alter the DNA landscape for potential transcriptional factors to attach, resulting in changes in transcriptional factor binding sites (TFBS). These TFBS changes are examined with respect to the human diseases which have been found to be significantly associated with the rSNPs.

Ultra-rare mutation in long-range enhancer predisposes to thyroid carcinoma with high penetrance

Thyroid cancer shows high heritability but causative genes remain largely unknown. According to a common hypothesis the genetic predisposition to thyroid cancer is highly heterogeneous; being in part due to many different rare alleles. Here we used linkage analysis and targeted deep sequencing to detect a novel single-nucleotide mutation in chromosome 4q32 (4q32A>C) in a large pedigree displaying non-medullary thyroid carcinoma (NMTC). This mutation is generally ultra-rare; it was not found in 38 NMTC families, in 2676 sporadic NMTC cases or 2470 controls. The mutation is located in a long-range enhancer element whose ability to bind the transcription factors POU2F and YY1 is significantly impaired, with decreased activity in the presence of the C- allele compared with the wild type A-allele. An enhancer RNA (eRNA) is transcribed in thyroid tissue from this region and is greatly downregulated in NMTC tumors. We suggest that this is an example of an ultra-rare mutation predisposing to thyroid cancer with high penetrance.

Big effects of small RNAs: a review of microRNAs in anxiety

Epigenetic and regulatory elements provide an additional layer of complexity to the heterogeneity of anxiety disorders. MicroRNAs (miRNAs) are a class of small, noncoding RNAs that have recently drawn interest as epigenetic modulators of gene expression in psychiatric disorders. miRNAs elicit their effects by binding to target messenger RNAs (mRNAs) and hindering translation or accelerating degradation. Considering their role in neuronal differentiation and synaptic plasticity, miRNAs have opened up new investigative avenues in the aetiology and treatment of anxiety disorders. In this review, we provide a thorough analysis of miRNAs, their targets and their functions in the central nervous system (CNS), focusing on their role in anxiety disorders. The involvement of miRNAs in CNS functions (such as neurogenesis, neurite outgrowth, synaptogenesis and synaptic and neural plasticity) and their intricate regulatory role under stressful conditions strongly support their importance in the aetiology of anxiety disorders. Furthermore, miRNAs could provide new avenues for the development of therapeutic targets in anxiety disorders.

The -308 G>A SNP of TNFA is a factor predisposing to chronic rhinosinusitis associated with nasal polyposis in aspirin-sensitive Hungarian individuals: conclusions of a genetic study with multiple stratifications

Single nucleotide polymorphisms (SNPs) of the tumour necrosis factor alpha (TNFα) gene (TNFA) have been extensively studied and shown to be associated with an increased risk of the development of various chronic inflammatory diseases. Inflammation has been demonstrated to play a central role in the pathogenesis of chronic rhinosinusitis (CRS), and TNFα is a key pro-inflammatory cytokine with important functions in these processes. In order to determine whether the well-known TNFA -308 G>A SNP has a role in a genetic predisposition to CRS in the Hungarian population, we analyzed our genomic collection containing control and CRS patient samples in a case-control study, and compared the genotype and allele frequencies. There was no significant difference in the observed genotype or allele frequencies between the controls and the total CRS group. However, after careful stratification of the patient group on the basis of the observed clinical symptoms, we found a significantly higher carriage rate of the rare A allele-containing genotypes among the CRS patients with nasal polyposis (NP) who also exhibited sensitivity to aspirin (acetylsalicylic acid, ASA(+)). It is concluded that genetic variants of the TNFA gene may affect the risk of CRS in a clinically well-defined group of CRSNP(+)ASA(+) patients in the Hungarian population. Our results also emphasize that the group of CRS patients is not homogenous in that patients exhibiting different clinical symptoms exist. Their carried genetic predisposing factors, and as a result, the exact molecular events leading to the development of various forms of CRS, may also differ.

Population genetics of cis-regulatory sequences that operate during embryonic development in the sea urchin Strongylocentrotus purpuratus

Despite the fact that noncoding sequences comprise a substantial fraction of functional sites within all genomes, the evolutionary mechanisms that operate on genetic variation within regulatory elements remain poorly understood. In this study, we examine the population genetics of the core, upstream cis-regulatory regions of eight genes (AN, CyIIa, CyIIIa, Endo16, FoxB, HE, SM30 a, and SM50) that function during the early development of the purple sea urchin, Strongylocentrotus purpuratus. Quantitative and qualitative measures of segregating variation are not conspicuously different between cis-regulatory and closely linked "proxy neutral" noncoding regions containing no known functional sites. Length and compound mutations are common in noncoding sequences; conventional descriptive statistics ignore such mutations, under-representing true genetic variation by approximately 28% for these loci in this population. Patterns of variation in the cis-regulatory regions of six of the genes examined (CyIIa, CyIIIa, Endo16, FoxB, AN, and HE) are consistent with directional selection. Genetic variation within annotated transcription factor binding sites is comparable to, and frequently greater than, that of surrounding sequences. Comparisons of two paralog pairs (CyIIa/CyIIIa and AN/HE) suggest that distinct evolutionary processes have operated on their cis-regulatory regions following gene duplication. Together, these analyses provide a detailed view of the evolutionary mechanisms operating on noncoding sequences within a natural population, and underscore how little is known about how these processes operate on cis-regulatory sequences.

Interpreting noncoding genetic variation in complex traits and human disease

Association studies provide genome-wide information about the genetic basis of complex disease, but medical research has primarily focused on protein-coding variants, due to the difficulty of interpreting non-coding mutations. This picture has changed with advances in the systematic annotation of functional non-coding elements. Evolutionary conservation, functional genomics, chromatin state, sequence motifs, and molecular quantitative trait loci all provide complementary information about non-coding function. These functional maps can help prioritize variants on risk haplotypes, filter mutations encountered in the clinic, and perform systems-level analyses to reveal processes underlying disease associations. Advances in predictive modeling can enable dataset integration to reveal pathways shared across loci and alleles, and richer regulatory models can guide the search for epistatic interactions. Lastly, new massively parallel reporter experiments can systematically validate regulatory predictions. Ultimately, advances in regulatory and systems genomics can help unleash the value of whole-genome sequencing for personalized genomic risk assessment, diagnosis, and treatment.

Investigation of allelic heterogeneity of the CCK-A receptor gene in paranoid schizophrenia

The cholecystokinin type A receptor (CCKAR) gene has been found to be associated with positive symptoms in patients with schizophrenia but the results reported to date are inconsistent. Considering the involvement of allelic heterogeneity in poor replication of the CCKAR finding, we genotyped five single nucleotide polymorphisms (SNPs) located in the 5' putative regulatory region of the CCKAR gene in a Chinese case-control sample and then applied the 5-SNP haplotype analysis to extract allelic heterogeneity information. The results showed that three individual haplotypes were strongly associated with increased risk of schizophrenia (corrected P = 2.9 × 10(-4), P = 2.5 × 10(-5), and P = 1.4 × 10(-5), respectively) and their combination gave an odds ratio (OR) of 6.12 with 95% CI 3.67-10.21 (P = 6.7 × 10(-15)). The haplotypes were also associated with some clinical symptoms including hallucination, suspiciousness, and hostility. Our work provided further evidence in support of the CCKAR hypothesis of schizophrenia and also suggested that haplotype-based association analysis may be a powerful approach for identification of allelic heterogeneity of a disease-underlying gene, which is very likely to be attributable to poor replication of an initial finding due to the reduction of sample power and the complexity of genetic architectures.

Identification of functional DNA variants in the constitutive promoter region of MDM2

Although mutations in the oncoprotein murine double minute 2 (MDM2) are rare, MDM2 gene overexpression has been observed in several human tumors. Given that even modest changes in MDM2 levels might influence the p53 tumor suppressor signaling pathway, we postulated that sequence variation in the promoter region of MDM2 could lead to disregulated expression and variation in gene dosage. Two promoters have been reported for MDM2; an internal promoter (P2), which is located near the end of intron 1 and is p53-responsive, and an upstream constitutive promoter (P1), which is p53-independent. Both promoter regions contain DNA variants that could influence the expression levels of MDM2, including the well-studied single nucleotide polymorphism (SNP) SNP309, which is located in the promoter P2; i.e., upstream of exon 2. In this report, we screened the promoter P1 for DNA variants and assessed the functional impact of the corresponding SNPs. Using the dbSNP database and genotyping validation in individuals of European descent, we identified three common SNPs (-1494 G > A; indel 40 bp; and -182 C > G). Three major promoter haplotypes were inferred by using these three promoter SNPs together with rs2279744 (SNP309). Following subcloning into a gene reporter system, we found that two of the haplotypes significantly influenced MDM2 promoter activity in a haplotype-specific manner. Site-directed mutagenesis experiments indicated that the 40 bp insertion/deletion variation is causing the observed allelic promoter activity. This study suggests that part of the variability in the MDM2 expression levels could be explained by allelic p53-independent P1 promoter activity.

Association between single nucleotide polymorphisms in ERCC4 and risk of squamous cell carcinoma of the head and neck

BACKGROUND

Excision repair cross-complementation group 4 gene (ERCC4/XPF) plays an important role in nucleotide excision repair and participates in removal of DNA interstrand cross-links and DNA double-strand breaks. Single nucleotide polymorphisms (SNPs) in ERCC4 may impact repair capacity and affect cancer susceptibility.

METHODOLOGY/PRINCIPAL FINDINGS

In this case-control study, we evaluated associations of four selected potentially functional SNPs in ERCC4 with risk of squamous cell carcinoma of the head and neck (SCCHN) in 1,040 non-Hispanic white patients with SCCHN and 1,046 cancer-free matched controls. We found that the variant GG genotype of rs2276466 was significantly associated with a decreased risk of SCCHN (OR = 0.69, 95% CI 0.50-0.96), and that the variant TT genotype of rs3136038 showed a borderline significant decreased risk with SCCHN (OR = 0.76, 95% CI: 0.58-1.01) in the recessive model. Such protective effects were more evident in oropharyngeal cancer (OR = 0.61, 95% CI: 0.40-0.92 for rs2276466; OR = 0.69, 95% CI: 0.48-0.98 for rs3136038). No significant associations were found for the other two SNPs (rs1800067 and rs1799798). In addition, individuals with the rs2276466 GG or with the rs3136038 TT genotypes had higher levels of ERCC4 mRNA expression than those with the corresponding wild-type genotypes in 90 Epstein-Barr virus-transformed lymphoblastoid cell lines derived from Caucasians.

CONCLUSIONS

These results suggest that these two SNPs (rs2276466 and rs3136038) in ERCC4 may be functional and contribute to SCCHN susceptibility. However, our findings need to be replicated in further large epidemiological and functional studies.