Extreme context specificity in differential allelic expression

Allelic imbalance of HLA-B expression in human lung cells infected with coronavirus and other respiratory viruses

The human leucocyte antigen (HLA) loci have been widely characterized to be associated with viral infectious diseases using either HLA allele frequency-based association or in silico predicted studies. However, there is less experimental evidence to link the HLA alleles with COVID-19 and other respiratory infectious diseases, particularly in the lung cells. To examine the role of HLA alleles in response to coronavirus and other respiratory viral infections in disease-relevant cells, we designed a two-stage study by integrating publicly accessible RNA-seq data sets, and performed allelic expression (AE) analysis on heterozygous HLA genotypes. We discovered an increased AE pattern accompanied with overexpression of HLA-B gene in SARS-CoV-2-infected human lung epithelial cells. Analysis of independent data sets verified the respiratory virus-induced AE of HLA-B gene in lung cells and tissues. The results were further experimentally validated in cultured lung cells infected with SARS-CoV-2. We further uncovered that the antiviral cytokine IFNβ contribute to AE of the HLA-B gene in lung cells. Our analyses provide a new insight into allelic influence on the HLA expression in association with SARS-CoV-2 and other common viral infectious diseases.

Synonymous GATA2 mutations result in selective loss of mutated RNA and are common in patients with GATA2 deficiency

Deficiency of the transcription factor GATA2 is a highly penetrant genetic disorder predisposing to myelodysplastic syndromes (MDS) and immunodeficiency. It has been recognized as the most common cause underlying primary MDS in children. Triggered by the discovery of a recurrent synonymous GATA2 variant, we systematically investigated 911 patients with phenotype of pediatric MDS or cellular deficiencies for the presence of synonymous alterations in GATA2. In total, we identified nine individuals with five heterozygous synonymous mutations: c.351C>G, p.T117T (N = 4); c.649C>T, p.L217L; c.981G>A, p.G327G; c.1023C>T, p.A341A; and c.1416G>A, p.P472P (N = 2). They accounted for 8.2% (9/110) of cases with GATA2 deficiency in our cohort and resulted in selective loss of mutant RNA. While for the hotspot mutation (c.351C>G) a splicing error leading to RNA and protein reduction was identified, severe, likely late stage RNA loss without splicing disruption was found for other mutations. Finally, the synonymous mutations did not alter protein function or stability. In summary, synonymous GATA2 substitutions are a new common cause of GATA2 deficiency. These findings have broad implications for genetic counseling and pathogenic variant discovery in Mendelian disorders.

Allele-specific expression identified rs2509956 as a novel long-distance cis-regulatory SNP for SCGB1A1, an important gene for multiple pulmonary diseases

SCGB1A1 (secretoglobin family 1A member 1) is an important protein for multiple pulmonary diseases, especially asthma, chronic obstructive pulmonary disease, and lung cancer. One single-nucleotide polymorphism (SNP) at 5'-untranslated region of SCGB1A1, rs3741240, has been suggested to be associated with reduced protein expression and further asthma susceptibility. However, it was still unclear whether there were other cis-regulatory elements for SCGB1A1 that might further contribute to pulmonary diseases. Allele-specific expression (ASE) is a novel approach to identify the functional region in human genome. In the present study, we measured ASE on rs3741240 in lung tissues and observed a consistent excess of G allele over A (P < 10^-6), which indicated that this SNP or the one(s) in linkage disequilibrium (LD) could regulate SCGB1A1 expression. By analyzing 1000 Genomes Project data for Chinese, one SNP locating ~10.2 kb away and downstream of SCGB1A1, rs2509956, was identified to be in strong LD with rs3741240. Reporter gene assay confirmed that both SNPs could regulate gene expression in the lung cell. By chromosome conformation capture, it was verified that the region surrounding rs2509956 could interact with SCGB1A1 promoter region and act as an enhancer. Through chromatin immunoprecipitation and overexpression assay, the related transcription factor RELA (RELA proto-oncogene, NF-kB subunit) was recognized to bind the region spanning rs2509956. Our work identified a novel long-distance cis-regulatory SNP for SCGB1A1, which might contribute to multiple pulmonary diseases.

Expression analysis of the osteoarthritis genetic susceptibility mapping to the matrix Gla protein gene MGP

BACKGROUND

Osteoarthritis (OA) is a common disease of older individuals that impacts detrimentally on the quality and the length of life. It is characterised by the painful loss of articular cartilage and is polygenic and multifactorial. Genome-wide association scans have highlighted over 90 osteoarthritis genetic signals, some of which reside within or close to highly plausible candidate genes. An example is an association to polymorphisms within and adjacent to the matrix Gla protein gene MGP. We set out to undertake a functional study of this gene.

METHODS

Nucleic acid was extracted from cartilage, infrapatellar fat pad, synovium, trabecular bone, trapezium and peripheral whole blood from OA patients and also from mesenchymal stem cells (MSCs) subjected to chondrogenesis. Expression of MGP was measured by quantitative PCR (qPCR), RNA-sequencing and allelic expression imbalance (AEI) analysis. Matrix Gla protein was depleted in chondrocytes by knocking down MGP expression using RNA interference (RNAi) and the effect on a range of genes assessed by qPCR.

RESULTS

MGP is expressed in joint tissues, blood and chondrocytes cultured from MSCs. There is a higher expression in diseased versus non-diseased cartilage. Polymorphisms that are associated with OA also correlate with the expression of MGP, with the OA risk-conferring allele showing significantly reduced expression in cartilage, fat pad and synovium but increased expression in blood. Depletion of Matrix Gla protein had a significant effect on the majority of genes tested, with an increased expression of catabolic genes that encode enzymes that degrade cartilage.

CONCLUSIONS

MGP expression is subject to cis-acting regulators that correlate with the OA association signal. These are active in a range of joint tissues but have effects which are particularly strong in cartilage. An opposite effect is observed in blood, highlighting the context-specific nature of the regulation of this gene's expression. Recapitulation of the genetic deficit in cartilage chondrocytes is pro-catabolic.

Quantification of allelic differential expression using a simple Fluorescence primer PCR-RFLP-based method

Allelic differential expression (ADE) is common in diploid organisms, and is often the key reason for specific phenotype variations. Thus, ADE detection is important for identification of major genes and causal mutations. To date, sensitive and simple methods to detect ADE are still lacking. In this study, we have developed an accurate, simple, and sensitive method, named fluorescence primer PCR-RFLP quantitative method (fPCR-RFLP), for ADE analysis. This method involves two rounds of PCR amplification using a pair of primers, one of which is double-labeled with an overhang 6-FAM. The two alleles are then separated by RFLP and quantified by fluorescence density. fPCR-RFLP could precisely distinguish ADE cross a range of 1- to 32-fold differences. Using this method, we verified PLAG1 and KIT, two candidate genes related to growth rate and immune response traits of pigs, to be ADE both at different developmental stages and in different tissues. Our data demonstrates that fPCR-RFLP is an accurate and sensitive method for detecting ADE on both DNA and RNA level. Therefore, this powerful tool provides a way to analyze mutations that cause ADE.

Identification of a novel, methylation-dependent, RUNX2 regulatory region associated with osteoarthritis risk

Osteoarthritis (OA) is a common, multifactorial and polygenic skeletal disease that, in its severest form, requires joint replacement surgery to restore mobility and to relieve chronic pain. Using tissues from the articulating joints of 260 patients with OA and a range of in vitro experiments, including CRISPR-Cas9, we have characterized an intergenic regulatory element. Here, genotype at an OA risk locus correlates with differential DNA methylation, with altered gene expression of both a transcriptional regulator (RUNX2), and a chromatin remodelling protein (SUPT3H). RUNX2 is a strong candidate for OA susceptibility, with its encoded protein being essential for skeletogenesis and healthy joint function. The OA risk locus includes single nucleotide polymorphisms (SNPs) located within and flanking the differentially methylated region (DMR). The OA association SNP, rs10948172, demonstrates particularly strong correlation with methylation, and two intergenic SNPs falling within the DMR (rs62435998 and rs62435999) demonstrate genetic and epigenetic effects on the regulatory activity of this region. We therefore posit that the OA signal mediates its effect by modulating the methylation of the regulatory element, which then impacts on gene expression, with RUNX2 being the principal target. Our study highlights the interplay between DNA methylation, OA genetic risk and the downstream regulation of genes critical to normal joint function.

Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage

OBJECTIVE

Multiple single-nucleotide polymorphisms (SNPs) conferring susceptibility to osteoarthritis (OA) mark imbalanced expression of positional genes in articular cartilage, reflected by unequally expressed alleles among heterozygotes (allelic imbalance [AI]). We undertook this study to explore the articular cartilage transcriptome from OA patients for AI events to identify putative disease-driving genetic variation.

METHODS

AI was assessed in 42 preserved and 5 lesioned OA cartilage samples (from the Research Arthritis and Articular Cartilage study) for which RNA sequencing data were available. The count fraction of the alternative alleles among the alternative and reference alleles together (φ) was determined for heterozygous individuals. A meta-analysis was performed to generate a meta-φ and P value for each SNP with a false discovery rate (FDR) correction for multiple comparisons. To further validate AI events, we explored them as a function of multiple additional OA features.

RESULTS

We observed a total of 2,070 SNPs that consistently marked AI of 1,031 unique genes in articular cartilage. Of these genes, 49 were found to be significantly differentially expressed (fold change <0.5 or >2, FDR <0.05) between preserved and paired lesioned cartilage, and 18 had previously been reported to confer susceptibility to OA and/or related phenotypes. Moreover, we identified notable highly significant AI SNPs in the CRLF1, WWP2, and RPS3 genes that were related to multiple OA features.

CONCLUSION

We present a framework and resulting data set for researchers in the OA research field to probe for disease-relevant genetic variation that affects gene expression in pivotal disease-affected tissue. This likely includes putative novel compelling OA risk genes such as CRLF1, WWP2, and RPS3.

Functional Characterization of the Osteoarthritis Genetic Risk Residing at ALDH1A2 Identifies rs12915901 as a Key Target Variant

OBJECTIVE

To identify the functional single-nucleotide polymorphisms (SNPs) and mechanisms conferring increased risk of hand osteoarthritis (OA) at the ALDH1A2 locus, which is a retinoic acid regulatory gene.

METHODS

Tissue samples from 247 patients with knee, hip, or hand OA who had undergone joint surgery were included. RNA-sequencing analysis was used to investigate differential expression of ALDH1A2 and other retinoic acid signaling pathway genes in cartilage. Expression of ALDH1A2 in joint tissues obtained from multiple sites was quantified using quantitative reverse transcription-polymerase chain reaction. Allelic expression imbalance (AEI) was measured by pyrosequencing. The consequences of ALDH1A2 depletion by RNA interference were assessed in primary human chondrocytes. In silico and in vitro analyses were used to pinpoint which, among 62 highly correlated SNPs, could account for the association at the locus.

RESULTS

ALDH1A2 expression was observed across multiple joint tissue samples, including osteochondral tissue from the hand. The expression of ALDH1A2 and of several retinoic acid signaling genes was different in diseased cartilage compared to non-diseased cartilage, with ALDH1A2 showing lower levels in OA cartilage. Experimental depletion of ALDH1A2 resulted in changes in the expression levels of a number of chondrogenic markers, including SOX9. In addition, reduced expression of the OA risk-conferring allele was witnessed in a number of joint tissues, with the strongest effect in cartilage. The intronic SNP rs12915901 recapitulated the AEI observed in patient tissues, while the Ets transcription factors were identified as potential mediators of this effect.

CONCLUSION

The ALDH1A2 locus seems to increase the risk of hand OA through decreased expression of ALDH1A2 in joint tissues, with the effect dependent on rs12915901. These findings indicate a mechanism that may now be targeted to modulate OA risk.

Variation in the Oxytocin Receptor Gene Predicts Brain Region-Specific Expression and Social Attachment

BACKGROUND

Oxytocin (OXT) modulates several aspects of social behavior. Intranasal OXT is a leading candidate for treating social deficits in patients with autism spectrum disorder, and common genetic variants in the human OXTR gene are associated with emotion recognition, relationship quality, and autism spectrum disorder. Animal models have revealed that individual differences in Oxtr expression in the brain drive social behavior variation. Our understanding of how genetic variation contributes to brain OXTR expression is very limited.

METHODS

We investigated Oxtr expression in monogamous prairie voles, which have a well-characterized OXT system. We quantified brain region-specific levels of Oxtr messenger RNA and oxytocin receptor protein with established neuroanatomic methods. We used pyrosequencing to investigate allelic imbalance of Oxtr mRNA, a molecular signature of polymorphic genetic regulatory elements. We performed next-generation sequencing to discover variants in and near the Oxtr gene. We investigated social attachment using the partner preference test.

RESULTS

Our allelic imbalance data demonstrate that genetic variants contribute to individual differences in Oxtr expression, but only in particular brain regions, including the nucleus accumbens, where oxytocin receptor signaling facilitates social attachment. Next-generation sequencing identified one polymorphism in the Oxtr intron, near a putative cis-regulatory element, explaining 74% of the variance in striatal Oxtr expression specifically. Males homozygous for the high expressing allele display enhanced social attachment.

CONCLUSIONS

Taken together, these findings provide convincing evidence for robust genetic influence on Oxtr expression and provide novel insights into how noncoding polymorphisms in OXTR might influence individual differences in human social cognition and behavior.

Correlation of the osteoarthritis susceptibility variants that map to chromosome 20q13 with an expression quantitative trait locus operating on NCOA3 and with functional variation at the polymorphism rs116855380

Objective

To functionally characterize the osteoarthritis (OA) susceptibility variants that map to a region of high linkage disequilibrium (LD) on chromosome 20q13 marked by the single‐nucleotide polymorphism (SNP) rs6094710 and encompassing NCOA3 and SULF2.

Methods

Nucleic acids were extracted from the cartilage of OA patients. Overall and allelic expression of NCOA3 and SULF2 were measured by quantitative reverse transcription–polymerase chain reaction and pyrosequencing, respectively. The functional effect of SNPs within the 20q13 locus was assessed in vitro using luciferase reporter constructs and electrophoretic mobility shift assays (EMSAs). The in vivo effect of nuclear receptor coactivator 3 (NCOA3) protein depletion on primary human OA articular cartilage chondrocytes was assessed using RNA interference.

Results

Expression of NCOA3 correlated with the genotype at rs6094710 (P = 0.006), and the gene demonstrated allelic expression imbalance (AEI) in individuals heterozygous for the SNP (mean AEI 1.21; P < 0.0001). In both instances, expression of the OA‐associated allele was reduced. In addition, there was reduced enhancer activity of the OA‐associated allele of rs116855380, a SNP in perfect LD with rs6094710 in luciferase assays (P < 0.001). EMSAs demonstrated a protein complex binding with reduced affinity to this allele. Depletion of NCOA3 led to significant changes (all P < 0.05) in the expression of genes involved in cartilage homeostasis.

Conclusion

NCOA3 is subject to a cis‐acting expression quantitative trait locus in articular cartilage, which correlates with the OA association signal and with the OA‐associated allele of the functional SNP rs116855380, a SNP that is located only 10.3 kb upstream of NCOA3. These findings elucidate the effect of the association of the 20q13 region on OA cartilage and provide compelling evidence of a potentially causal candidate SNP.

Differential allelic expression of SOS1 and hyperexpression of the activating SOS1 c.755C variant in a Noonan syndrome family

Noonan syndrome (NS) is a genetic condition characterized by congenital heart defects, short stature and characteristic facial features. We here present the case of a girl with moderate learning disabilities, delayed language development, craniofacial features and skin anomalies reminiscent of NS. After a mutation screening of the known NS genes PTPN11, SOS1, RAF1, KRAS, GRB2, BRAF and SHOC2 we found the heterozygous c.755T>C variant in SOS1 causing the p.I252T amino-acid substitution, which was considered possibly pathogenetic by bioinformatic predictions. The same variant was present in the proband's mother, displaying some NS features, and maternal grandfather showing no NS traits, but also by a healthy subject in 1000 genomes project database without phenotype informations. The functional analysis revealed that SOS1 c.755C activated the RAS-ERK intracellular pathway, whereas no effects on RAC-JNK cascade have been detected. After a comparison between the sequence of SOS1 cDNA from peripheral blood and SOS1 genomic DNA, we showed for the first time a differential allelic expression of the SOS1 gene in healthy individuals, thus occurring as a physiologic condition. Interestingly, we found that the mutated allele C was 50% more expressed than the wild-type allele T in all familial carriers. The comparable amount of SOS1 mRNA between mutated individuals and the controls indicates that the variant does not affect SOS1 expression. The present study provides a first evidence of allelic imbalance of SOS1 and pinpoints this condition as a possible mechanism underlying a different penetrance of some SOS1-mutated alleles in unrelated carriers.

Characterization of the cartilage DNA methylome in knee and hip osteoarthritis

Objective

The aim of this study was to characterize the genome-wide DNA methylation profile of chondrocytes from knee and hip cartilage obtained from patients with osteoarthritis (OA) and hip cartilage obtained from patients with femoral neck fracture, providing the first comparison of DNA methylation between OA and non-OA hip cartilage, and between OA hip and OA knee cartilage.

Methods

The study was performed using the Illumina Infinium HumanMethylation450 BeadChip array, which allows the annotation of ∼480,000 CpG sites. Genome-wide methylation was assessed in chondrocyte DNA extracted from 23 hip OA patients, 73 knee OA patients, and 21 healthy hip control patients with femoral neck fracture.

Results

Analysis revealed that chondrocytes from the hip cartilage of OA patients and healthy controls have unique methylation profiles, with 5,322 differentially methylated loci (DMLs) identified between the 2 groups. In addition, a comparison between hip and knee OA chondrocytes revealed 5,547 DMLs between the 2 groups, including DMLs in several genes known to be involved in the pathogenesis of OA. Hip OA samples were found to cluster into 2 groups. A total of 15,239 DMLs were identified between the 2 clusters, with an enrichment of genes involved in inflammation and immunity. Similarly, we confirmed a previous report of knee OA samples that also clustered into 2 groups.

Conclusion

We demonstrated that global DNA methylation using a high-density array can be a powerful tool in the characterization of OA at the molecular level. Identification of pathways enriched in DMLs between OA and OA-free cartilage highlight potential etiologic mechanisms that are involved in the initiation and/or progression of the disease and that could be therapeutically targeted.

Allelic expression analysis of the osteoarthritis susceptibility locus that maps to chromosome 3p21 reveals cis-acting eQTLs at GNL3 and SPCS1

Background

An osteoarthritis (OA) susceptibility locus has been mapped to chromosome 3p21, to a region of high linkage disequilibrium encompassing twelve genes. Six of these genes are expressed in joint tissues and we therefore assessed whether any of the six were subject to cis-acting regulatory polymorphisms active in these tissues and which could therefore account for the association signal.

Methods

We measured allelic expression using pyrosequencing assays that can distinguish mRNA output from each allele of a transcript single nucleotide polymorphism. We assessed RNA extracted from the cartilage and other joint tissues of OA patients who had undergone elective joint replacement surgery. A two-tailed Mann–Whitney exact test was used to test the significance of any allelic differences.

Results

GNL3 and SPCS1 demonstrated significant allelic expression imbalance (AEI) in OA cartilage (GNL3, mean AEI = 1.04, p = 0.0002; SPCS1, mean AEI = 1.07, p < 0.0001). Similar results were observed in other tissues. Expression of the OA-associated allele was lower than that of the non-associated allele for both genes.

Conclusions

cis-acting regulatory polymorphisms acting on GNL3 and SPCS1 contribute to the OA association signal at chromosome 3p21, and these genes therefore merit further investigation.

Identification and analysis of a SMAD3 cis-acting eQTL operating in primary osteoarthritis and in the aneurysms and osteoarthritis syndrome

OBJECTIVE

The TGF-β pathway plays a central role in joint development with polymorphism in TGF-β pathway genes implicated in osteoarthritis susceptibility. One association is to rs12901499, within intron 1 of SMAD3. Since rs12901499 is not in linkage disequilibrium with a non-synonymous polymorphism, it is likely the association is operating by influencing expression of SMAD3.

DESIGN

Using tissues from the joints of primary osteoarthritis patients who had undergone joint replacement we measured the overall expression of SMAD3 by quantitative real-time PCR. We also measured allelic expression of SMAD3 using these tissues and vascular smooth muscle cells from patients with aneurysms and osteoarthritis syndrome, a rare condition featuring early-onset osteoarthritis. We tested the functional effect of SNPs in vitro using luciferase assays and assessed association with osteoarthritis using a large osteoarthritis case-control dataset.

RESULTS

We observed that genotype at rs12901499 did not correlate with overall SMAD3 expression or allelic expression. However, genotype at a 3'UTR SNP, rs8031440, did correlate with SMAD3 expression in cartilage (P = 0.005) which was supported by allelic expression data showing that the G allele correlated with decreased SMAD3 expression in joint tissues and vascular smooth muscle cells. This G allele was underrepresented in osteoarthritis cases vs controls (P = 0.027, odds ratio = 0.921). rs8031440 is in perfect linkage disequilibrium with five other SMAD3 3'UTR SNPs and our luciferase analysis identified rs3743342 and rs12595334 as being functional.

CONCLUSION

SMAD3 is subject to cis-acting regulatory polymorphism in the tissues of relevance to both primary osteoarthritis and the aneurysms-osteoarthritis syndrome.

Evidence of allelic imbalance in the schizophrenia susceptibility gene ZNF804A in human dorsolateral prefrontal cortex

The rs1344706, an intronic SNP within the zinc-finger protein 804A gene (ZNF804A), was identified as one of the most compelling risk SNPs for schizophrenia (SZ) and bipolar disorder (BD). It is however not clear by which molecular mechanisms ZNF804A increases disease risk. We evaluated the role of ZNF804A in SZ and BD by genotyping the originally associated rs1344706 SNP and an exonic SNP (rs12476147) located in exon four of ZNF804A in a sample of 422 SZ, 382 BD, and 507 controls from the isolated population of the Costa Rica Central Valley. We also investigated the rs1344706 SNP for allelic specific expression (ASE) imbalance in the dorsolateral prefrontal cortex (DLPFC) of 46 heterozygous postmortem brains. While no significant association between rs1344706 and SZ or BD was observed in the Costa Rica sample, we observed an increased risk of SZ for the minor allele (A) of the exonic rs12476147 SNP (p=0.026). Our ASE assay detected a significant over-expression of the rs12476147 A allele in DLPFC of rs1344706 heterozygous subjects. Interestingly, cDNA allele ratios were significantly different according to the intronic rs1344706 genotypes (p-value=0.03), with the rs1344706 A allele associated with increased ZNF804A rs12476147 A allele expression (average 1.06, p-value=0.02, for heterozygous subjects vs. genomic DNA). In conclusion, we have demonstrated a significant association of rs12476147 with SZ, and using a powerful within-subject design, an allelic expression imbalance of ZNF804A exonic SNP rs12476147 in the DLPFC. Although this data does not preclude the possibility of other functional variants in ZNF804A, it provides evidence that the rs1344706 SZ risk allele is the cis-regulatory variant directly responsible for this allelic expression imbalance in adult cortex.

Allelic expression analysis of the osteoarthritis susceptibility gene COL11A1 in human joint tissues

Background

The single nucleotide polymorphism (SNP) rs2615977 is associated with osteoarthritis (OA) and is located in intron 31 of COL11A1, a strong candidate gene for this degenerative musculoskeletal disease. Furthermore, the common non-synonymous COL11A1 SNP rs1676486 is associated with another degenerative musculoskeletal disease, lumbar disc herniation (LDH). rs1676486 is a C-T transition mediating its affect on LDH susceptibility by modulating COL11A1 expression. The risk T-allele of rs1676486 leads to reduced expression of the COL11A1 transcript, a phenomenon known as allelic expression imbalance (AEI). We were keen therefore to assess whether the effect that rs1676486 has on COL11A1 expression in LDH is also observed in OA and whether the rs2615977 association to OA also marked AEI.

Methods

Using RNA from OA cartilage, we assessed whether either SNP correlated with COL11A1 AEI by 1) measuring COL11A1 expression and stratifying the data by genotype at each SNP; and 2) quantifying the mRNA transcribed from each allele of the two SNPs. We also assessed whether rs1676486 was associated with OA susceptibility using a case–control cohort of over 18,000 individuals.

Results

We observed significant AEI at rs1676486 (p < 0.0001) with the T-allele correlating with reduced COL11A1 expression. This corresponded with observations in LDH but the SNP was not associated with OA. We did not observe AEI at rs2615977.

Conclusions

COL11A1 is subject to AEI in OA cartilage. AEI at rs1676486 is a risk factor for LDH, but not for OA. These two diseases therefore share a common functional phenotype, namely AEI of COL11A1, but this appears to be a disease risk only in LDH. Other functional effects on COL11A1 presumably account for the OA susceptibility that maps to this gene.

Multiple isoforms and differential allelic expression of CHRNA5 in lung tissue and lung adenocarcinoma

CHRNA5 gene expression variation may play a role in individual susceptibility to lung cancer. Analysis of CHRNA5 transcripts expressed in normal lung tissue detected the full-length transcript (isoform-1) and four splicing transcripts (isoform-2 to isoform-5), derived from the recognition of other splice sites in exon 5. Isoforms-2, -3 and -4 were found by protein modeling to form a completely folded, potentially functional extracellular domain and were observed at the protein level, whereas isoform-5 lacked a consistent part of the distorted β sandwich and was not seen at the protein level. Only isoform-1 appeared to encode a complete, functional subunit able to fulfill the ion channel function. We previously reported that CHRNA5 expression is associated with genetic polymorphisms at this locus and that three haplotypes in its promoter region show functional regulation in vitro. Analysis of differential allelic expression (DAE) of three single nucleotide polymorphisms (rs503464, rs55853698 and rs55781567) tagging the expression haplotypes of the CHRNA5 promoter indicated statistically significant DAE at rs55853698 and rs55781567, in both normal lung and lung adenocarcinoma. Overall, our findings provide evidence for the presence of multiple CHRNA5 messenger RNA (mRNA) isoforms that may modulate the multimeric nicotine receptor and cis-regulatory variations in the CHRNA5 locus that act in vivo in the control of CHRNA5 mRNA expression, in normal lung tissue and in lung adenocarcinoma.

Association of a BMP5 microsatellite with knee osteoarthritis: case-control study

Introduction

We aimed to explore the involvement of a multiallelic functional polymorphism in knee osteoarthritis (OA) susceptibility as a prototype of possible genetic factors escaping GWAS detection.

Methods

OA patients and controls from three European populations (Greece, Spain and the UK) adding up to 1003 patients (716 women, 287 men) that had undergone total knee joint replacement (TKR) due to severe primary OA and 1543 controls (758 women, 785 men) lacking clinical signs or symptoms of OA were genotyped for the D6S1276 microsatellite in intron 1 of BMP5. Genotype and mutiallelic trend tests were used to compare cases and controls.

Results

Significant association was found between the microsatellite and knee OA in women (P from 3.1 x10^-4 to 4.1 x10^-4 depending on the test), but not in men. Three of the alleles showed significant differences between patients and controls, one of them of increased risk and two of protection. The gender association and the allele direction of change were very concordant with those previously reported for hip OA.

Conclusions

We have found association of knee OA in women with the D6S1276 functional microsatellite that modifies in cis the expression of BMP5 making this a sounder OA genetic factor and extending its involvement to other joints. This result also shows the interest of analysing other multiallelic polymorphisms.

Evidence of cis-acting regulatory variation in PTPN22 in patients with rheumatoid arthritis

OBJECTIVES

To assess whether there are cis-regulatory polymorphisms that regulate protein tyrosine phosphatase, non-receptor type 22 (PTPN22) expression in rheumatoid arthritis (RA).

METHODS

RNA was extracted from positively selected CD56+, CD8+, and CD4+ mononuclear cells and the 'residual' cells from 12 RA patients heterozygous for the PTPN22 C1858T single nucleotide polymorphism (SNP) (rs2476601). Relative allelic expression was measured by single base extension (SBE) assay.

RESULTS

There was relative differential allelic expression (DAE ≥ 20%) in eight patients (p < 10(-5)); seven patients demonstrated DAE in more than one cell type; four patients had statistically significant differences between these cell populations (p(corrected) < 0.05).

CONCLUSIONS

We have demonstrated significant differences in expression of PTPN22 alleles in RA patients, indicating the probable existence of cis-acting regulatory elements.

Allelic expression analysis of the osteoarthritis susceptibility locus that maps to MICAL3

Background

A genome-wide association scan with subsequent replication study that involved over 67,000 individuals of European ancestry has produced evidence of association of single nucleotide polymorphism rs2277831 to primary osteoarthritis (OA) with a P-value of 2.9 × 10^-5. rs2277831, an A/G transition, is located in an intron of MICAL3. This gene is located on chromosome 22q11.21 and the association signal encompasses two additional genes, BCL2L13 and BID. It is becoming increasingly apparent that many common complex traits are mediated by cis-acting regulatory polymorphisms that influence, in a tissue-specific manner, gene expression or transcript stability.

Methods

We used total and allelic expression analysis to assess whether the OA association to rs2277831 is mediated by an influence on MICAL3, BCL2L13 or BID expression. Using RNA extracted from joint tissues of 60 patients who had undergone elective joint replacement surgery, we assessed whether rs2277831 correlated with allelic expression of either of the three genes by: 1) measuring the expression of each gene by quantitative PCR and then stratifying the data by genotype at rs2277831 and 2) accurately discriminating and quantifying the mRNA synthesised from the alleles of OA patients using allelic-quantitative PCR.

Results

We found no evidence for a correlation between gene expression and genotype at rs2277831, with P-values of 0.09 for BCL2L13, 0.07 for BID and 0.33 for MICAL3. In the allelic expression analysis we observed several examples of significant (p < 0.05) allelic imbalances, with an allelic expression ratio of 2.82 observed in BCL2L13 (P = 0.004), 2.09 at BID (P = 0.001) and the most extreme case being at MICAL3, with an allelic expression ratio of 5.47 (P = 0.001). However, there was no correlation observed between the pattern of allelic expression and the genotype at rs2277831.

Conclusions

In the tissues that we have studied, our data do not support our hypothesis that the association between rs2277831 and OA is due to the effect this SNP has on MICAL3, BCL2L13 or BID gene expression. Instead, our data point towards other functional effects accounting for the OA associated signal.

Comparing methods for mapping cis acting polymorphisms using allelic expression ratios

Genome wide association studies frequently reveal associations between disease susceptibility and polymorphisms outside coding regions. Such associations cannot always be explained by linkage disequilibrium with changes affecting the transcription products. This has stimulated the interest in characterising sequence variation influencing gene expression levels, in particular in changes acting in cis. Differences in transcription between the two alleles at an autosomal locus can be used to test the association between candidate polymorphisms and the modulation of gene expression in cis. This type of approach requires at least one transcribed polymorphism and one candidate polymorphism. In the past five years, different methods have been proposed to analyse such data. Here we use simulations and real data sets to compare the power of some of these methods. The results show that when it is not possible to determine the phase between the transcribed and potentially cis acting allele there is some advantage in using methods that estimate phased genotype and effect on expression simultaneously. However when the phase can be determined, simple regression models seem preferable because of their simplicity and flexibility. The simulations and the analysis of experimental data suggest that in the majority of situations, methods that assume a lognormal distribution of the allelic expression ratios are both robust to deviations from this assumption and more powerful than alternatives that do not make these assumptions.

Analysis of allelic differential expression in the human genome using allele-specific serial analysis of gene expression tags

Recent reports have demonstrated that a significant proportion of human genes display allelic differential expression (ADE). ADE is associated with phenotypic variability and may contribute to complex genetic diseases. Here, we present a computational analysis of ADE using allele-specific serial analysis of gene expression (SAGE) tags representing 1295 human genes. We identified 472 genes for which unequal representation (>3-fold) of allele-specific SAGE tags was observed in at least one SAGE library, suggesting the occurrence of ADE. For 235 out of these 472 genes, the difference in the expression level between both allele-specific SAGE tags was statistically significant (p < 0.05). Eleven candidate genes were then subjected to experimental validation and ADE was confirmed for 8 out of these 11 genes. Our results suggest that at least 25% of the human genes display ADE and that allele-specific SAGE tags can be efficiently used for the identification of such genes.

Detecting differential allelic expression using high-resolution melting curve analysis: application to the breast cancer susceptibility gene CHEK2

BACKGROUND

The gene CHEK2 encodes a checkpoint kinase playing a key role in the DNA damage pathway. Though CHEK2 has been identified as an intermediate breast cancer susceptibility gene, only a small proportion of high-risk families have been explained by genetic variants located in its coding region. Alteration in gene expression regulation provides a potential mechanism for generating disease susceptibility. The detection of differential allelic expression (DAE) represents a sensitive assay to direct the search for a functional sequence variant within the transcriptional regulatory elements of a candidate gene. We aimed to assess whether CHEK2 was subject to DAE in lymphoblastoid cell lines (LCLs) from high-risk breast cancer patients for whom no mutation in BRCA1 or BRCA2 had been identified.

METHODS

We implemented an assay based on high-resolution melting (HRM) curve analysis and developed an analysis tool for DAE assessment.

RESULTS

We observed allelic expression imbalance in 4 of the 41 LCLs examined. All four were carriers of the truncating mutation 1100delC. We confirmed previous findings that this mutation induces non-sense mediated mRNA decay. In our series, we ruled out the possibility of a functional sequence variant located in the promoter region or in a regulatory element of CHEK2 that would lead to DAE in the transcriptional regulatory milieu of freely proliferating LCLs.

CONCLUSIONS

Our results support that HRM is a sensitive and accurate method for DAE assessment. This approach would be of great interest for high-throughput mutation screening projects aiming to identify genes carrying functional regulatory polymorphisms.

No association between germline allele-specific expression of TGFBR1 and colorectal cancer risk in Caucasian and Ashkenazi populations

Background:

Germline allele-specific expression (ASE) of the TGFBR1 gene has been reported as a strong risk factor for colorectal cancer (CRC) with an odds ratio close to 9. Considering the potential implications of the finding, we undertook the task of validating the initial results in this study.

Methods:

Allele-specific expression was measured using the highly quantitative and robust technique of pyrosequencing. Individuals from two different populations were studied, one Caucasian-dominated and the other of Ashkenazi Jewish descent, with different sources of non-tumoral genetic material in each.

Results:

Our results showed no statistically significant differences in the degree of ASE between CRC patients and controls, considering ASE as either a quantitative or a binary trait. Using defined cutoff values to categorise ASE, 1.0% of blood lymphocytes from informative Israeli cases (total n=96) were ASE positive (median 1.00; range 0.76–1.31) and 2.2% of informative matched controls (total n=90) were ASE positive (median 1.00; range 0.76–1.87). Likewise, normal mucosae from Spanish patients (median 1.03; range: 0.68–1.43; n=75) did not show significant differences in the degree of ASE when compared with the Israeli patients or controls.

Conclusions:

Taken together, these results suggest that ASE of TGFBR1 does not confer an increased risk of CRC.

Functionality of promoter microsatellites of arginine vasopressin receptor 1A (AVPR1A): implications for autism

BACKGROUND

Arginine vasopressin (AVP) has been hypothesized to play a role in aetiology of autism based on a demonstrated involvement in the regulation of social behaviours. The arginine vasopressin receptor 1A gene (AVPR1A) is widely expressed in the brain and is considered to be a key receptor for regulation of social behaviour. Moreover, genetic variation at AVPR1A has been reported to be associated with autism. Evidence from non-human mammals implicates variation in the 5'-flanking region of AVPR1A in variable gene expression and social behaviour.

METHODS

We examined four tagging single nucleotide polymorphisms (SNPs) (rs3803107, rs1042615, rs3741865, rs11174815) and three microsatellites (RS3, RS1 and AVR) at the AVPR1A gene for association in an autism cohort from Ireland. Two 5'-flanking region polymorphisms in the human AVPR1A, RS3 and RS1, were also tested for their effect on relative promoter activity.

RESULTS

The short alleles of RS1 and the SNP rs11174815 show weak association with autism in the Irish population (P = 0.036 and P = 0.008, respectively). Both RS1 and RS3 showed differences in relative promoter activity by length. Shorter repeat alleles of RS1 and RS3 decreased relative promoter activity in the human neuroblastoma cell line SH-SY5Y.

CONCLUSIONS

These aligning results can be interpreted as a functional route for this association, namely that shorter alleles of RS1 lead to decreased AVPR1A transcription, which may proffer increased susceptibility to the autism phenotype.

Effects of cis-regulatory variation differ across regions of the adult human brain

Cis-regulatory variation is considered to be an important determinant of human phenotypic variability, including susceptibility to complex disease. Recent studies have shown that the effects of cis-regulatory polymorphism on gene expression can differ widely between tissues. In the present study, we tested whether the effects of cis-regulatory variation can also differ between regions of the adult human brain. We used relative allelic expression to measure cis-effects on the RNA expression of five candidate genes for neuropsychiatric illness (ZNF804A, NOS1, RGS4, AKT1 and TCF4) across multiple discrete brain regions within individual subjects. For all five genes, we observed significant differences in allelic expression between brain regions in several individual subjects, suggesting regional differences in the effects of cis-regulatory polymorphism to be a common phenomenon. As well as highlighting an important caveat for studies of regulatory polymorphism in the brain, our findings indicate that it is possible to delineate brain areas in which cis-regulatory variants are active. This may provide important insights into the fundamental biology of neuropsychiatric phenotypes with which such variants are associated.

BMP5 activates multiple signaling pathways and promotes chondrogenic differentiation in the ATDC5 growth plate model

The bone morphogenetic protein 5 (BMP5) participates in skeletal development but its direct effects on the function of growth plate chondrocytes during chondrogenesis have not been explored. We have investigated the signaling pathways activated by BMP5 and its effect on chondrogenic differentiation in the ATDC5 growth plate chondrocyte model. BMP5 transiently activated p38 mitogen-activated protein kinase (MAPK) and extracellular signal-regulated kinase signaling after 10 days of differentiation; sustained Smad and p38 MAPK signaling were seen after 15 days differentiation. All three pathways were activated by BMP5 in human adult articular chondrocytes. BMP5 alone and in combination with the chondrogenic enhancer, insulin, induced proteoglycan synthesis, aggrecan core protein 1 expression, and alkaline phosphatase activity. Upregulation of hypertrophic markers parathyroid receptor 1 and collagen type X alpha 1 occurred in BMP5-treated ATDC5 cultures. BMP5 is clearly chondrogenic and exhibits stage-specific regulation of multiple signaling pathways in this growth plate model. In particular, BMP5 accelerates expression of hypertrophy markers which is of relevance in both development and diseases such as osteoarthritis.

Allelic imbalance (AI) identifies novel tissue-specific cis-regulatory variation for human UGT2B15

Allelic imbalance (AI) is a powerful tool to identify cis-regulatory variation for gene expression. UGT2B15 is an important enzyme involved in the metabolism of multiple endobiotics and xenobiotics. In this study, we measured the relative expression of two alleles at this gene by using SNP rs1902023:G>T. An excess of the G over the T allele was consistently observed in liver (P<0.001), but not in breast (P=0.06) samples, suggesting that SNPs in strong linkage disequilibrium with G253T regulate UGT2B15 expression in liver. Seven such SNPs were identified by resequencing the promoter and exon 1, which define two distinct haplotypes. Reporter gene assays confirmed that one haplotype displayed approximately 20% higher promoter activity compared to the other major haplotype in liver HepG2 (P<0.001), but not in breast MCF-7 (P=0.540) cells. Reporter gene assays with additional constructs pointed to rs34010522:G>T and rs35513228:C>T as the cis-regulatory variants; both SNPs were also evaluated in LNCaP and Caco-2 cells. By ChIP, we showed that the transcription factor Nrf2 binds to the region spanning rs34010522:G>T in all four cell lines. Our results provide a good example for how AI can be used to identify cis-regulatory variation and gain insights into the tissue specific regulation of gene expression.

Association of a functional microsatellite within intron 1 of the BMP5 gene with susceptibility to osteoarthritis

Background

In a previous study carried out by our group, the genotyping of 36 microsatellite markers from within a narrow interval of chromosome 6p12.3-q13 generated evidence for linkage and for association to female hip osteoarthritis (OA), with the most compelling association found for a marker within intron 1 of the bone morphogenetic protein 5 gene (BMP5). In this study, we aimed to further categorize the association of variants within intron 1 of BMP5 with OA through an expanded genetic association study of the intron and subsequent functional analysis of associated polymorphisms.

Methods

We genotyped 18 common polymorphisms including 8 microsatellites and 9 single nucleotide polymorphisms (SNPs) and 1 insertion/deletion (INDEL) from within highly conserved regions between human and mouse within intron 1 of BMP5. These markers were then tested for association to OA by a two-stage approach in which the polymorphisms were initially genotyped in a case-control cohort comprising 361 individuals with associated polymorphisms (P ≤ 0.05) then genotyped in a second case-control cohort comprising 1185 individuals.

Results

Two BMP5 intron 1 polymorphisms demonstrated association in the combined case-control cohort of 1546 individuals (765 cases and 781 controls): microsatellite D6S1276 (P = 0.018) and SNP rs921126 (P = 0.013). Functional analyses in osteoblastic, chondrocytic, and adipocytic cell lines indicated that allelic variants of D6S1276 have significant effects on the transcriptional activity of the BMP5 promoter in vitro.

Conclusion

Variability in gene expression of BMP5 may be an important contributor to OA genetic susceptibility.

Extent of differential allelic expression of candidate breast cancer genes is similar in blood and breast

Introduction

Normal gene expression variation is thought to play a central role in inter-individual variation and susceptibility to disease. Regulatory polymorphisms in cis-acting elements result in the unequal expression of alleles. Differential allelic expression (DAE) in heterozygote individuals could be used to develop a new approach to discover regulatory breast cancer susceptibility loci. As access to large numbers of fresh breast tissue to perform such studies is difficult, a suitable surrogate test tissue must be identified for future studies.

Methods

We measured differential allelic expression of 12 candidate genes possibly related to breast cancer susceptibility (BRCA1, BRCA2, C1qA, CCND3, EMSY, GPX1, GPX4, MLH3, MTHFR, NBS1, TP53 and TRXR2) in breast tissue (n = 40) and fresh blood (n = 170) of healthy individuals and EBV-transformed lymphoblastoid cells (n = 19). Differential allelic expression ratios were determined by Taqman assay. Ratio distributions were compared using t-test and Wilcoxon rank sum test, for mean ratios and variances respectively.

Results

We show that differential allelic expression is common among these 12 candidate genes and is comparable between breast and blood (fresh and transformed lymphoblasts) in a significant proportion of them. We found that eight out of nine genes with DAE in breast and fresh blood were comparable, as were 10 out of 11 genes between breast and transformed lymphoblasts.

Conclusions

Our findings support the use of differential allelic expression in blood as a surrogate for breast tissue in future studies on predisposition to breast cancer.

Functional analysis of the osteoarthritis susceptibility-associated GDF5 regulatory polymorphism

OBJECTIVE

Single-nucleotide polymorphism (SNP) rs143383 (T to C) in the 5'-untranslated region (5'-UTR) of GDF5 has recently been reported to be associated with osteoarthritis (OA) susceptibility, with lower expression of the risk-associated T allele observed in vitro and in vivo. The in vivo studies were performed on cartilage tissue from OA patients. The present study was undertaken to expand the analysis of the effect of this SNP on GDF5 allelic expression to more joint tissue types, to investigate for cis and trans factors that interact with the SNP, and to examine novel cis-acting GDF5 regulatory polymorphisms.

METHODS

Tissue samples were collected from OA patients undergoing joint replacement of the hip or knee. Nucleic acid was extracted, and, using rs143383 and an assay that discriminates and quantifies allelic expression, the relative amount of GDF5 expression from the T and C alleles was measured. Additional common variants in the GDF5 transcript sequence were interrogated as potential regulatory elements using allelic expression and luciferase reporter assays, and electrophoretic mobility shift assays were used to search for trans factors binding to rs143383.

RESULTS

We observed a consistent allelic expression imbalance of GDF5 in all tissues tested, implying that the functional effect mediated by rs143383 on GDF5 expression is joint-wide. We identified a second polymorphism, located in the 3'-UTR of GDF5, that influenced allelic expression of the gene independent of rs143383. Finally, we observed differential binding of deformed epidermal autoregulatory factor 1 (DEAF-1) to the 2 alleles of rs143383.

CONCLUSION

These findings show that the OA susceptibility mediated by polymorphism in GDF5 is not restricted to cartilage, emphasizing the need to consider the disease as involving the whole joint. The existence of an additional cis-acting regulatory polymorphism highlights the complexity of the regulation of expression of this important OA susceptibility locus. DEAF-1 is a trans-acting factor that merits further investigation as a potential tool for modulating GDF5 expression.

Infrequent detection of germline allele-specific expression of TGFBR1 in lymphoblasts and tissues of colon cancer patients

Recently, germline allele-specific expression (ASE) of the gene encoding for transforming growth factor-beta type I receptor (TGFBR1) has been proposed to be a major risk factor for cancer predisposition in the colon. Germline ASE results in a lowered expression of one of the TGFBR1 alleles (>1.5-fold), and was shown to occur in approximately 20% of informative familial and sporadic colorectal cancer (CRC) cases. In the present study, using the highly quantitative pyrosequencing technique, we estimated the frequency of ASE in TGFBR1 in a cohort of affected individuals from familial clusters of advanced colon neoplasias (cancers and adenomas with high-grade dysplasia), and also from a cohort of individuals with sporadic CRCs. Cases were considered positive for the presence of ASE if demonstrating an allelic expression ratio <0.67 or >1.5. Using RNA derived from lymphoblastoid cell lines, we find that of 46 informative Caucasian advanced colon neoplasia cases with a family history, only 2 individuals display a modest ASE, with allelic ratios of 1.65 and 1.73, respectively. Given that ASE of TGFBR1, if present, would likely be more pronounced in the colon compared with other tissues, we additionally determined the allele ratios of TGFBR1 in the RNA derived from normal-appearing colonic mucosa of sporadic CRC cases. We, however, found no evidence of ASE in any of 44 informative sporadic cases analyzed. Taken together, we find that germline ASE of TGFBR1, as assayed in lymphoblastoid and colon epithelial cells of colon cancer patients, is a relatively rare event.

Polymorphisms affecting gene transcription and mRNA processing in pharmacogenetic candidate genes: detection through allelic expression imbalance in human target tissues

BACKGROUND

Genetic variation in mRNA expression plays a critical role in human phenotypic diversity, but it has proven difficult to detect regulatory polymorphisms - mostly single nucleotide polymorphisms (rSNPs). Additionally, variants in the transcribed region, termed here 'structural RNA SNPs' (srSNPs), can affect mRNA processing and turnover. Both rSNPs and srSNPs cause allelic mRNA expression imbalance (AEI) in heterozygous individuals. We have used AEI to discover and characterize regulatory polymorphisms in OPRM1, TPH2, MDR1, DRD2, and VKORC1. The objective of this study was to use AEI to determine the extent of cis-regulatory factors in pharmacogenetic genes.

METHODS

We applied a rapid and accurate AEI methodology for testing 42 genes implicated in cardiovascular and central nervous system diseases, and affecting drug metabolism and transport. Each gene was analyzed in physiologically relevant human autopsy tissues, including brain, heart, liver, intestines, and lymphocytes.

RESULTS

Substantial AEI was observed in approximately 55% of the surveyed genes. Focusing on cardiovascular candidate genes in human hearts, AEI analysis revealed frequent cis-acting regulatory factors in ACE and SOD2 mRNA expression, having potential clinical significance. SNP scanning to locate regulatory polymorphisms in a number of genes failed to support several previously proposed promoter SNPs discovered with use of reporter gene assays in heterologous tissues, while srSNPs appear more frequent than expected. Computational analysis of mRNA folding indicates that approximately 90% of srSNPs affect mRNA folding, and hence potentially function.

CONCLUSION

Our results indicate that both rSNPs and srSNPs represent a still largely untapped reservoir of variants that contribute to human phenotypic diversity.

Dynamic variation in allele-specific gene expression of Paraoxonase-1 in murine and human tissues

Differential allelic expression has been shown to be common in mice, humans and maize, and variability in the expression of polymorphic alleles has been associated with human disease. Here, we describe the differential expression pattern of Paraoxonase-1, a gene involved in lipid metabolism and implicated in the formation of atherosclerotic lesions. We measured the expression of the murine Paraoxonase-1 gene (Pon1) in livers at different stages of embryonic development using F1 hybrid crosses and quantified the transcriptional level of both parental alleles. Using human foetal tissues, we analysed the expression of the human orthologue (PON1) and found monoallelic or preferential allelic expression in 6/7 and 4/4 samples from liver and pancreas, respectively. We observed that Pon1 does not show a parent-of-origin preference in its allelic expression, but has dramatic variations in allele-specific expression occurring throughout development. This study has important repercussions in the analysis of haplotypes at disease loci, since it implies that the expression of polymorphic alleles can be unequal and dynamic.

Genome-wide association scan identifies a prostaglandin-endoperoxide synthase 2 variant involved in risk of knee osteoarthritis

Osteoarthritis (OA), the most prevalent form of arthritis in the elderly, is characterized by the degradation of articular cartilage and has a strong genetic component. Our aim was to identify genetic variants involved in risk of knee OA in women. A pooled genome-wide association scan with the Illumina550 Duo array was performed in 255 controls and 387 cases. Twenty-eight variants with p < 1 x 10(-5) were estimated to have probabilities of being false positives <or=0.5 and were genotyped individually in the original samples and in replication cohorts from the UK and the U.S. (599 and 272 cases, 1530 and 258 controls, respectively). The top seven associations were subsequently tested in samples from the Netherlands (306 cases and 584 controls). rs4140564 on chromosome 1 mapping 5' to both the PTGS2 and PLA2G4A genes was associated with risk of knee OA in all the cohorts studied (overall odds ratio OR(mh) = 1.55 95% C.I. 1.30-1.85, p < 6.9 x 10(-7)). Differential allelic expression analysis of PTGS2 with mRNA extracted from the cartilage of joint-replacement surgery OA patients revealed a significant difference in allelic expression (p < 1.0 x 10(-6)). These results suggest the existence of cis-acting regulatory polymorphisms that are in, or near to, PTGS2 and in modest linkage disequilibrium with rs4140564. Our results and previous studies on the role of the cyclooxygenase 2 enzyme encoded by PTGS2 underscore the importance of this signaling pathway in the pathogenesis of knee OA.

Gene expression in the etiology of schizophrenia

Gene expression represents a fundamental interface between genes and environment in the development and ongoing plasticity of the human brain. Individual differences in gene expression are likely to underpin much of human diversity, including psychiatric illness. In the past decade, the development of microarray and proteomic technology has enabled global description of gene expression in schizophrenia. However, it is difficult on the basis of gene expression assays alone to distinguish between those changes that constitute primary etiology and those that reflect secondary pathology, compensatory mechanisms, or confounding influences. In this respect, tests of genetic association with schizophrenia will be instructive because changes in gene expression that result from gene variants that are associated with the disorder are likely to be of primary etiological significance. However, regulatory polymorphism is extremely difficult to recognize on the basis of sequence interrogation alone. Functional assays at the messenger RNA and/or protein level will be essential in elucidating the molecular mechanisms underlying genetic association with schizophrenia and are likely to become increasingly important in the identification of regulatory variants with which to test for association with the disorder and related traits. Once established, etiologically relevant changes in gene expression can be recapitulated in model systems in order to elucidate the molecular and physiological pathways that may ultimately give rise to the condition.

An SNP in the 5′-UTR of GDF5 is associated with osteoarthritis susceptibility in Europeans and with in vivo differences in allelic expression in articular cartilage

A compelling genetic association with osteoarthritis (OA) of a functional SNP (rs143383, T/C) in the 5'-UTR of the GDF5 gene was recently reported in case-control cohorts from Japan and China. GDF5 is a pro-chondrogenic growth factor. The T-allele frequency of the gene was elevated in cases, with an odds ratio (OR) of 1.79, and in vitro functional studies demonstrated that this allele mediated a moderate but significant reduction in the activity of the GDF5 promoter in several cell lines. Our initial objective was to assess whether the SNP was also associated with OA in a broad European population by genotyping the SNP in 2487 cases and 2018 age-matched controls from the UK and Spain. The T-allele was associated with OA (P = 0.03, OR = 1.10) as was carrier status for this allele (P = 0.004, OR = 1.28), demonstrating that the SNP is associated with OA in two diverse ethnic groups, Asians and Europeans. We subsequently assessed the functional effect of the SNP on GDF5 allelic expression using RNA extracted from the cartilage of OA patients who had undergone joint-replacement surgery. The associated T-allele showed up to a 27% reduction in expression relative to the C-allele (P = 0.00007), revealing that the functional effect mediated by SNP rs143383 on GDF5 expression is active in patients who have severe disease up to the point at which they require surgery. A small but persistent imbalance of GDF5 expression throughout life therefore appears to render an individual more susceptible to OA.

A survey of genomic properties for the detection of regulatory polymorphisms

Advances in the computational identification of functional noncoding polymorphisms will aid in cataloging novel determinants of health and identifying genetic variants that explain human evolution. To date, however, the development and evaluation of such techniques has been limited by the availability of known regulatory polymorphisms. We have attempted to address this by assembling, from the literature, a computationally tractable set of regulatory polymorphisms within the ORegAnno database (http://www.oreganno.org). We have further used 104 regulatory single-nucleotide polymorphisms from this set and 951 polymorphisms of unknown function, from 2-kb and 152-bp noncoding upstream regions of genes, to investigate the discriminatory potential of 23 properties related to gene regulation and population genetics. Among the most important properties detected in this region are distance to transcription start site, local repetitive content, sequence conservation, minor and derived allele frequencies, and presence of a CpG island. We further used the entire set of properties to evaluate their collective performance in detecting regulatory polymorphisms. Using a 10-fold cross-validation approach, we were able to achieve a sensitivity and specificity of 0.82 and 0.71, respectively, and we show that this performance is strongly influenced by the distance to the transcription start site.