Whole blood transcriptional profiling in ankylosing spondylitis identifies novel candidate genes that might contribute to the inflammatory and tissue-destructive disease aspects

Introduction: A number of genetic-association studies have identified genes contributing to ankylosing spondylitis (AS) susceptibility but such approaches provide little information as to the gene activity changes occurring during the disease process. Transcriptional profiling generates a 'snapshot' of the sampled cells' activity and thus can provide insights into the molecular processes driving the disease process. We undertook a whole-genome microarray approach to identify candidate genes associated with AS and validated these gene-expression changes in a larger sample cohort. Methods: A total of 18 active AS patients, classified according to the New York criteria, and 18 gender- and age-matched controls were profiled using Illumina HT-12 whole-genome expression BeadChips which carry cDNAs for 48,000 genes and transcripts. Class comparison analysis identified a number of differentially expressed candidate genes. These candidate genes were then validated in a larger cohort using qPCR-based TaqMan low density arrays (TLDAs). Results: A total of 239 probes corresponding to 221 genes were identified as being significantly different between patients and controls with a P-value <0.0005 (80% confidence level of false discovery rate). Forty-seven genes were then selected for validation studies, using the TLDAs. Thirteen of these genes were validated in the second patient cohort with 12 downregulated 1.3- to 2-fold and only 1 upregulated (1.6-fold). Among a number of identified genes with well-documented inflammatory roles we also validated genes that might be of great interest to the understanding of AS progression such as SPOCK2 (osteonectin) and EP300, which modulate cartilage and bone metabolism. Conclusions: We have validated a gene expression signature for AS from whole blood and identified strong candidate genes that may play roles in both the inflammatory and joint destruction aspects of the disease.

Genetics of ankylosing spondylitis

Purpose of Review Over the past 3 years, several new genes and gene deserts have been identified that are associated with ankylosing spondylitis (AS). The purpose of this review is to discuss the major findings of these studies, and the answers they provide and questions they raise about the pathogenesis of this common condition. Recent Findings: Five genes/genetic regions have now definitively been associated with AS [the major histocompatibility complex (MHC), IL23R, ERAP1, 2p15 and 21q22]. Strong evidence to support association with the disease has been demonstrated for the genes IL1R2, ANTXR2, TNFSF15, TNFR1 and a region on chromosome 16q including the gene TRADD. There is an overrepresentation of genes involved in Th17 lymphocyte differentiation/activation among genes associated with AS and the related diseases inflammatory bowel disease and psoriasis, pointing strongly to this pathway as playing a major causative role in the disease. Increasing information about differential association of HLA-B27 subtypes with disease suggests a hierarchy of strength of association of those alleles with AS, providing a useful test as to the validity of different potential mechanisms of association of HLA-B27 with AS. The mechanism underlying the association of the gene deserts, 2p15 and 21q22, suggests the involvement of noncoding RNA in AS etiopathogenesis. Summary: The increasing list of genes identified as being definitely involved in AS provides a useful platform for hypothesis-driven research in the field, providing a potential alternative route to determining the underlying mechanisms involved in the disease to research focusing on HLA-B27 alone.

Statistical analysis of the associations of hereditary factors with the risk of Type 1 diabetes and Diabetic Nephropathy based on the familial data collected through ascertaiment from population-based registers

In genetic epidemiology, population-based disease registries are commonly used to collect genotype or other risk factor information concerning affected subjects and their relatives. This work presents two new approaches for the statistical inference of ascertained data: a conditional and full likelihood approaches for the disease with variable age at onset phenotype using familial data obtained from population-based registry of incident cases. The aim is to obtain statistically reliable estimates of the general population parameters. The statistical analysis of familial data with variable age at onset becomes more complicated when some of the study subjects are non-susceptible, that is to say these subjects never get the disease. A statistical model for a variable age at onset with long-term survivors is proposed for studies of familial aggregation, using latent variable approach, as well as for prospective studies of genetic association studies with candidate genes. In addition, we explore the possibility of a genetic explanation of the observed increase in the incidence of Type 1 diabetes (T1D) in Finland in recent decades and the hypothesis of non-Mendelian transmission of T1D associated genes. Both classical and Bayesian statistical inference were used in the modelling and estimation. Despite the fact that this work contains five studies with different statistical models, they all concern data obtained from nationwide registries of T1D and genetics of T1D. In the analyses of T1D data, non-Mendelian transmission of T1D susceptibility alleles was not observed. In addition, non-Mendelian transmission of T1D susceptibility genes did not make a plausible explanation for the increase in T1D incidence in Finland. Instead, the Human Leucocyte Antigen associations with T1D were confirmed in the population-based analysis, which combines T1D registry information, reference sample of healthy subjects and birth cohort information of the Finnish population. Finally, a substantial familial variation in the susceptibility of T1D nephropathy was observed. The presented studies show the benefits of sophisticated statistical modelling to explore risk factors for complex diseases.

Evaluation of polymorphisms in predicted target sites for micro RNAs differentially expressed in endometriosis

Previous microarray analyses identified 22 microRNAs (miRNAs) differentially expressed in paired ectopic and eutopic endometrium of women with and without endometriosis. To investigate further the role of these miRNAs in women with endometriosis, we conducted an association study aiming to explore the relationship between endometriosis risk and single-nucleotide polymorphisms (SNPs) in miRNA target sites for these differentially expressed miRNAs. A panel of 102 SNPs in the predicted miRNA binding sites were evaluated for an endometriosis association study and an ingenuity pathway analysis was performed. Fourteen rare variants were identified in this study. We found SNP rs14647 in the Wolf-Hirschhorn syndrome candidate gene1 (WHSC1) 3'UTR (untranslated region) was associated with endometriosis-related infertility presenting an odds ratio of 12.2 (95% confidence interval = 2.4-60.7, P = 9.03 x 10(-5)). SNP haplotype AGG in the solute carrier family 22, member 23 (SLC22A23) 3'UTR was associated with endometriosis-related infertility and more severe disease. With the individual genotyping data, ingenuity pathways analysis identified the tumour necrosis factor and cyclin-dependant kinase inhibitor as major factors in the molecular pathways. Significant associations between WHSC1 alleles and endometriosis-related infertility and SLC22A23 haplotypes and the disease severe stage were identified. These findings may help focus future research on subphenotypes of this disease. Replication studies in independent large sample sets to confirm and characterize the involvement of the gene variation in the pathogenesis of endometriosis are needed.

Variation in BMPR1B, TGFRB1 and BMPR2 and control of dizygotic twinning

Genes in the TGF9 signaling pathway play important roles in the regulation of ovarian follicle growth and ovulation rate. Mutations in three genes in this pathway, growth differentiation factor 9 (GDF9), bone morphogenetic protein 15 (BMP15) and the bone morphogenetic protein receptor B 1 (BMPRB1), influence dizygotic (DZ) twinning rates in sheep. To date, only variants in GDF9 and BMP15, but not their receptors transforming growth factor ss receptor 1 (TGFBR1), bone morphogenetic protein receptor 2 (BMPR2) and BMPR1B, have been investigated with respect to their roles in human DZ twinning. We screened for rare and novel variants in TGFBR1, BMPR2 and BMPR1B in mothers of dizygotic twins (MODZT) from twin-dense families, and assessed association between genotyped and imputed variants and DZ twinning in another large sample of MODZT. Three novel variants were found: a deep intronic variant in BMPR2, and one intronic and one non-synonymous exonic variant in BMPRB1 which would result in the replacement of glutamine by glutamic acid at amino acid position 294 (p.Gln294Glu). None of these variants were predicted to have major impacts on gene function. However, the p.Gln294Glu variant changes the same amino acid as a sheep BMPR1B functional variant and may have functional consequences. Six BMPR1B variants were marginally associated with DZ twinning in the larger case-control sample, but these were no longer significant once multiple testing was taken into account. Our results suggest that variation in the TGF9 signaling pathway type II receptors has limited effects on DZ twinning rates in humans.

A rare P2X7 variant Arg307Gln with absent pore formation function protects against neuroinflammation in multiple sclerosis

Multiple sclerosis (MS) is a chronic relapsing-remitting inflammatory disease of the central nervous system characterized by oligodendrocyte damage, demyelination and neuronal death. Genetic association studies have shown a 2-fold or greater prevalence of the HLA-DRB1*1501 allele in the MS population compared with normal Caucasians. In discovery cohorts of Australasian patients with MS (total 2941 patients and 3008 controls), we examined the associations of 12 functional polymorphisms of P2X7, a microglial/macrophage receptor with proinflammatory effects when activated by extracellular adenosine triphosphate (ATP). In discovery cohorts, rs28360457, coding for Arg307Gln was associated with MS and combined analysis showed a 2-fold lower minor allele frequency compared with controls (1.11% for MS and 2.15% for controls, P = 0.0000071). Replication analysis of four independent European MS case–control cohorts (total 2140 cases and 2634 controls) confirmed this association [odds ratio (OR) = 0.69, P = 0.026]. A meta-analysis of all Australasian and European cohorts indicated that Arg307Gln confers a 1.8-fold protective effect on MS risk (OR = 0.57, P = 0.0000024). Fresh human monocytes heterozygous for Arg307Gln have >85% loss of ‘pore’ function of the P2X7 receptor measured by ATP-induced ethidium uptake. Analysis shows Arg307Gln always occurred with 270His suggesting a single 307Gln–270His haplotype that confers dominant negative effects on P2X7 function and protection against MS. Modeling based on the homologous zP2X4 receptor showed Arg307 is located in a region rich in basic residues located only 12 Å from the ligand binding site. Our data show the protective effect against MS of a rare genetic variant of P2RX7 with heterozygotes showing near absent proinflammatory ‘pore’ function.

Heritability of the network architecture of intrinsic brain functional connectivity

The brain's functional network exhibits many features facilitating functional specialization, integration, and robustness to attack. Using graph theory to characterize brain networks, studies demonstrate their small-world, modular, and "rich-club" properties, with deviations reported in many common neuropathological conditions. Here we estimate the heritability of five widely used graph theoretical metrics (mean clustering coefficient (γ), modularity (Q), rich-club coefficient (ϕnorm), global efficiency (λ), small-worldness (σ)) over a range of connection densities (k=5-25%) in a large cohort of twins (N=592, 84 MZ and 89 DZ twin pairs, 246 single twins, age 23±2.5). We also considered the effects of global signal regression (GSR). We found that the graph metrics were moderately influenced by genetic factors h2 (γ=47-59%, Q=38-59%, ϕnorm=0-29%, λ=52-64%, σ=51-59%) at lower connection densities (≤15%), and when global signal regression was implemented, heritability estimates decreased substantially h2 (γ=0-26%, Q=0-28%, ϕnorm=0%, λ=23-30%, σ=0-27%). Distinct network features were phenotypically correlated (|r|=0.15-0.81), and γ, Q, and λ were found to be influenced by overlapping genetic factors. Our findings suggest that these metrics may be potential endophenotypes for psychiatric disease and suitable for genetic association studies, but that genetic effects must be interpreted with respect to methodological choices.

Novel associations of UDP-glucuronosyltransferase 2B gene variants with prostate cancer risk in a multiethnic study.

BACKGROUND: We have previously shown that a functional polymorphism of the UGT2B15 gene (rs1902023) was associated with increased risk of prostate cancer (PC). Novel functional polymorphisms of the UGT2B17 and UGT2B15 genes have been recently characterized by in vitro assays but have not been evaluated in epidemiologic studies. METHODS: Fifteen functional SNPs of the UGT2B17 and UGT2B15 genes, including cis-acting UGT2B gene SNPs, were genotyped in African American and Caucasian men (233 PC cases and 342 controls). Regression models were used to analyze the association between SNPs and PC risk. RESULTS: After adjusting for race, age and BMI, we found that six UGT2B15 SNPs (rs4148269, rs3100, rs9994887, rs13112099, rs7686914 and rs7696472) were associated with an increased risk of PC in log-additive models (p < 0.05). A SNP cis-acting on UGT2B17 and UGT2B15 expression (rs17147338) was also associated with increased risk of prostate cancer (OR = 1.65, 95% CI = 1.00-2.70); while a stronger association among men with high Gleason sum was observed for SNPs rs4148269 and rs3100. CONCLUSIONS: Although small sample size limits inference, we report novel associations between UGT2B15 and UGT2B17 variants and PC risk. These associations with PC risk in men with high Gleason sum, more frequently found in African American men, support the relevance of genetic differences in the androgen metabolism pathway, which could explain, in part, the high incidence of PC among African American men. Larger studies are required.

Finding our way through phenotypes.

Despite a large and multifaceted effort to understand the vast landscape of phenotypic data, their current form inhibits productive data analysis. The lack of a community-wide, consensus-based, human- and machine-interpretable language for describing phenotypes and their genomic and environmental contexts is perhaps the most pressing scientific bottleneck to integration across many key fields in biology, including genomics, systems biology, development, medicine, evolution, ecology, and systematics. Here we survey the current phenomics landscape, including data resources and handling, and the progress that has been made to accurately capture relevant data descriptions for phenotypes. We present an example of the kind of integration across domains that computable phenotypes would enable, and we call upon the broader biology community, publishers, and relevant funding agencies to support efforts to surmount today's data barriers and facilitate analytical reproducibility.

Polymorphisms in the kinesin-like factor 1 B gene and risk of epithelial ovarian cancer in Eastern Chinese women.

The kinesin-like factor 1 B (KIF1B) gene plays an important role in the process of apoptosis and the transformation and progression of malignant cells. Genetic variations in KIF1B may contribute to risk of epithelial ovarian cancer (EOC). In this study of 1,324 EOC patients and 1,386 cancer-free female controls, we investigated associations between two potentially functional single nucleotide polymorphisms in KIF1B and EOC risk by the conditional logistic regression analysis. General linear regression model was used to evaluate the correlation between the number of variant alleles and KIF1B mRNA expression levels. We found that the rs17401966 variant AG/GG genotypes were significantly associated with a decreased risk of EOC (adjusted odds ratio (OR) = 0.81, 95 % confidence interval (CI) = 0.68-0.97), compared with the AA genotype, but no associations were observed for rs1002076. Women who carried both rs17401966 AG/GG and rs1002076 AG/AA genotypes of KIF1B had a 0.82-fold decreased risk (adjusted 95 % CI = 0.69-0.97), compared with others. Additionally, there was no evidence of possible interactions between about-mentioned co-variants. Further genotype-phenotype correlation analysis indicated that the number of rs17401966 variant G allele was significantly associated with KIF1B mRNA expression levels (P for GLM = 0.003 and 0.001 in all and Chinese subjects, respectively), with GG carriers having the lowest level of KIF1B mRNA expression. Taken together, the rs17401966 polymorphism likely regulates KIF1B mRNA expression and thus may be associated with EOC risk in Eastern Chinese women. Larger, independent studies are warranted to validate our findings.

Apoptotic variants as predictors of risk of oropharyngeal cancer recurrence after definitive radiotherapy.

Single nucleotide polymorphisms (SNPs) in the promoter region of FAS and FASLG may alter their transcriptional activity. Thus, we determined the associations between four FAS and FASLG promoter variants (FAS1377G>A, rs2234767; 670A>G, rs1800682; FASLG844T>C, rs763110 and 124A>G, rs5030772) and the risk of recurrence of squamous cell carcinoma of the oropharynx (SCCOP). We evaluated the associations between FAS and FASLG genetic variants and the risk of recurrence in a cohort of 1,008 patients. The log-rank test and multivariate Cox models were used to evaluate the associations. Compared with patients with common homozygous genotypes of FAS670 and FASLG844 polymorphisms, patients with variant genotypes had lower disease-free survival rates (log-rank p < 0.0001 and p < 0.0001, respectively) and an approximately threefold higher risk of SCCOP recurrence (HR, 3.2;95% CI, 2.2-4.6; and HR, 3.1; 95% CI, 2.2-4.4, respectively) after multivariate adjustment. Furthermore, among patients with HPV16-positive tumors, those with variant genotypes of these two polymorphisms had lower disease-free survival rates (log-rank, p < 0.0001 and p < 0.0001, respectively) and a higher recurrence risk than did patients with common homozygous genotypes (HR, 12.9; 95% CI, 3.8-43.6; and HR, 8.1; 95% CI, 3.6-18.6, respectively), whereas no significant associations were found for FAS1377 and FASLG124 polymorphisms. Our findings suggest that FAS670 and FASLG844 polymorphisms modulate the risk of recurrence of SCCOP, particularly in patients with HPV16-positive tumors. Larger studies are needed to validate these results.

The role of adiposity in cardiometabolic traits:a Mendelian randomization analysis

The association between adiposity and cardiometabolic traits is well known from epidemiological studies. Whilst the causal relationship is clear for some of these traits, for others it is not. We aimed to determine whether adiposity is causally related to various cardiometabolic traits using the Mendelian randomization approach.

The Cyclin D1 (CCND1) A870G polymorphism predicts clinical outcome to lapatinib and capecitabine in HER2-positive metastatic breast cancer

BACKGROUND: Lapatinib plus capecitabine emerged as an efficacious therapy in metastatic breast cancer (mBC). We aimed to identify germline single-nucleotide polymorphisms (SNPs) in genes involved in capecitabine catabolism and human epidermal receptor signaling that were associated with clinical outcome to assist in selecting patients likely to benefit from this combination.

PATIENTS AND METHODS: DNA was extracted from 240 of 399 patients enrolled in EGF100151 clinical trial (NCT00078572; clinicaltrials.gov) and SNPs were successfully evaluated in 234 patients. The associations between SNPs and clinical outcome were analyzed using Fisher's exact test, Kaplan-Meier curves, log-rank tests, likelihood ratio test within logistic or Cox regression model, as appropriate.

RESULTS: There were significant interactions between CCND1 A870G and clinical outcome. Patients carrying the A-allele were more likely to benefit from lapatinib plus capecitabine versus capecitabine when compared with patients harboring G/G (P = 0.022, 0.024 and 0.04, respectively). In patients with the A-allele, the response rate (RR) was significantly higher with lapatinib plus capecitabine (35%) compared with capecitabine (11%; P = 0.001) but not between treatments in patients with G/G (RR = 24% and 32%, respectively; P = 0.85). Time to tumor progression (TTP) was longer in patients with the A-allele treated with lapatinib plus capecitabine compared with capecitabine (median TTP = 7.9 and 3.4 months; P < 0.001), but not in patients with G/G (median TTP = 6.1 and 6.6 months; P = 0.92).

CONCLUSION: Our findings suggest that CCND1A870G may be useful in predicting clinical outcome in HER2-positive mBC patients treated with lapatinib plus capecitabine.

Prognostic role of lemur tyrosine kinase-3 germline polymorphisms in adjuvant gastric cancer in Japan and the United States

Lemur tyrosine kinase-3 (LMTK3) was recently identified as an estrogen receptor (ER)-α modulator related to endocrine therapy resistance, and its polymorphisms rs9989661 (T>C) T/T genotype and rs8108419 (G>A) G/G or A/G genotype predicted improved outcomes in breast cancer. Because different predominant ER distributions link to breast and gastric cancer and little is known of the prognostic role of LMTK3 in gastric cancer, this study was carried out to clarify the prognostic role of these polymorphisms in gastric cancer. One-hundred and sixty-nine Japanese and 137 U.S. patients with localized gastric adenocarcinoma were enrolled. Genomic DNA was extracted from blood or tissue, and all samples were analyzed by PCR-based direct DNA sequencing. Overall, these polymorphisms were not associated with survival in both cohorts. When gender was considered, in multivariate analysis, harboring rs9989661 T/T genotype was associated with disease-free survival [HR, 4.37; 95% confidence interval (CI), 2.08-9.18; P < 0.0001] and overall survival (OS; HR, 3.69; 95% CI, 1.65-8.24; P = 0.0014) in the Japanese males and time to recurrence (HR, 7.29; 95% CI, 1.07-49.80; P = 0.043) in the U.S. females. Meanwhile, harboring rs8108419 G/G genotype was associated with OS in the Japanese females (HR, 3.04; 95% CI, 1.08-8.56; P = 0.035) and the U.S. males (HR, 3.39; 95% CI, 1.31-8.80; P = 0.012). The prognostic role of these polymorphisms may be negative in gastric cancer. These findings suggest that the estrogen pathway may play a prognostic role in patients with gastric cancer but this may be dependent on the regional differences both in physiology and genetic alterations of gastric cancer.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.