Effects of immunomodulatory drugs on depressive symptoms: A mega-analysis of randomized, placebo-controlled clinical trials in inflammatory disorders.

Activation of the innate immune system is commonly associated with depression. Immunomodulatory drugs may have efficacy for depressive symptoms that are co-morbidly associated with inflammatory disorders. We report a large-scale re-analysis by standardized procedures (mega-analysis) of patient-level data combined from 18 randomized clinical trials conducted by Janssen or GlaxoSmithKline for one of nine disorders (N = 10,743 participants). Core depressive symptoms (low mood, anhedonia) were measured by the Short Form Survey (SF-36) or the Hospital Anxiety and Depression Scale (HADS), and participants were stratified into high (N = 1921) versus low-depressive strata based on baseline ratings. Placebo-controlled change from baseline after 4-16 weeks of treatment was estimated by the standardized mean difference (SMD) over all trials and for each subgroup of trials targeting one of 7 mechanisms (IL-6, TNF-α, IL-12/23, CD20, COX2, BLγS, p38/MAPK14). Patients in the high depressive stratum showed modest but significant effects on core depressive symptoms (SMD = 0.29, 95% CI [0.12-0.45]) and related SF-36 measures of mental health and vitality. Anti-IL-6 antibodies (SMD = 0.8, 95% CI [0.20-1.41]) and an anti-IL-12/23 antibody (SMD = 0.48, 95% CI [0.26-0.70]) had larger effects on depressive symptoms than other drug classes. Adjustments for physical health outcome marginally attenuated the average treatment effect on depressive symptoms (SMD = 0.20, 95% CI: 0.06-0.35), but more strongly attenuated effects on mental health and vitality. Effects of anti-IL-12/23 remained significant and anti-IL-6 antibodies became a trend after controlling for physical response to treatment. Novel immune-therapeutics can produce antidepressant effects in depressed patients with primary inflammatory disorders that are not entirely explained by treatment-related changes in physical health.

Group-based variant calling leveraging next-generation supercomputing for large-scale whole-genome sequencing studies.

MOTIVATION: Next-generation sequencing (NGS) technologies have become much more efficient, allowing whole human genomes to be sequenced faster and cheaper than ever before. However, processing the raw sequence reads associated with NGS technologies requires care and sophistication in order to draw compelling inferences about phenotypic consequences of variation in human genomes. It has been shown that different approaches to variant calling from NGS data can lead to different conclusions. Ensuring appropriate accuracy and quality in variant calling can come at a computational cost. RESULTS: We describe our experience implementing and evaluating a group-based approach to calling variants on large numbers of whole human genomes. We explore the influence of many factors that may impact the accuracy and efficiency of group-based variant calling, including group size, the biogeographical backgrounds of the individuals who have been sequenced, and the computing environment used. We make efficient use of the Gordon supercomputer cluster at the San Diego Supercomputer Center by incorporating job-packing and parallelization considerations into our workflow while calling variants on 437 whole human genomes generated as part of large association study. CONCLUSIONS: We ultimately find that our workflow resulted in high-quality variant calls in a computationally efficient manner. We argue that studies like ours should motivate further investigations combining hardware-oriented advances in computing systems with algorithmic developments to tackle emerging 'big data' problems in biomedical research brought on by the expansion of NGS technologies.

Comprehensive analysis of treatment response phenotypes in rheumatoid arthritis for pharmacogenetic studies.

BACKGROUND: An individual patient's response to a particular drug is influenced by multiple factors, which may include genetic predisposition. Pharmacogenetic studies attempt to discover and estimate the contributions of genetic variants to the variability in response to a drug treatment. The task of identifying the genetic contribution is often complicated by response phenotypes that are based on imprecise or subjective clinical observations. Because the success of a pharmacogenetic study depends on the analysis of a heritable phenotype, it is important to identify phenotypes with a significant heritable component to ensure reliable and reproducible results in subsequent genetic association studies. METHODS: We retrospectively analyzed data collected from 436 rheumatoid arthritis patients treated with golimumab during the phase III GO-FURTHER study. We investigated the reliability of several potential response outcomes after golimumab treatment. Using whole-genome sequencing of the clinical trial cohort, we estimated the heritability of each potential outcome measure. We further performed a longitudinal analysis of the clinical data to estimate variability of outcome measures over time and the degree to which each response metric could be confounded by placebo response. RESULTS: We determined that the high degree of within-patient variation over time makes a single follow-up visit insufficient to assess an individual patient's response to golimumab treatment. We found that different potential response outcomes had varying degrees of heritability and that averaging across multiple follow-up visits yielded higher heritability estimates than single follow-up estimates. Importantly, we found that the change in swollen and tender joint counts were the most heritable outcome metrics we tested; however, we showed that they are also more likely to be confounded by a placebo response than objective phenotypes like the change in C-reactive protein levels. CONCLUSIONS: Our rigorous approach to finding robust and heritable response phenotypes could be beneficial to all pharmacogenetic studies and may lead to more reliable and reproducible results. TRIAL REGISTRATION: Clinicaltrials.gov NCT00973479 . Registered 4 September 2009.

Integrative genomic deconvolution of rheumatoid arthritis GWAS loci into gene and cell type associations.

BACKGROUND: Although genome-wide association studies (GWAS) have identified over 100 genetic loci associated with rheumatoid arthritis (RA), our ability to translate these results into disease understanding and novel therapeutics is limited. Most RA GWAS loci reside outside of protein-coding regions and likely affect distal transcriptional enhancers. Furthermore, GWAS do not identify the cell types where the associated causal gene functions. Thus, mapping the transcriptional regulatory roles of GWAS hits and the relevant cell types will lead to better understanding of RA pathogenesis. RESULTS: We combine the whole-genome sequences and blood transcription profiles of 377 RA patients and identify over 6000 unique genes with expression quantitative trait loci (eQTLs). We demonstrate the quality of the identified eQTLs through comparison to non-RA individuals. We integrate the eQTLs with immune cell epigenome maps, RA GWAS risk loci, and adjustment for linkage disequilibrium to propose target genes of immune cell enhancers that overlap RA risk loci. We examine 20 immune cell epigenomes and perform a focused analysis on primary monocytes, B cells, and T cells. CONCLUSIONS: We highlight cell-specific gene associations with relevance to RA pathogenesis including the identification of FCGR2B in B cells as possessing both intragenic and enhancer regulatory GWAS hits. We show that our RA patient cohort derived eQTL network is more informative for studying RA than that from a healthy cohort. While not experimentally validated here, the reported eQTLs and cell type-specific RA risk associations can prioritize future experiments with the goal of elucidating the regulatory mechanisms behind genetic risk associations.

Modular analysis of peripheral blood gene expression in rheumatoid arthritis captures reproducible gene expression changes in tumor necrosis factor responders.

OBJECTIVE: To establish whether the analysis of whole-blood gene expression is useful in predicting or monitoring response to anti-tumor necrosis factor (anti-TNF) therapy in patients with rheumatoid arthritis (RA). METHODS: Whole-blood RNA (using a PAXgene system to stabilize whole-blood RNA in the collection tube) was obtained at baseline and at 14 weeks from 3 independent cohorts, consisting of a combined total of 240 RA patients who were beginning therapy with anti-TNF. We used an approach to gene expression analysis that is based on modular patterns of gene expression, or modules. RESULTS: Good and moderate responders according to the European League Against Rheumatism criteria exhibited highly significant and consistent changes in multiple gene expression modules after 14 weeks of therapy, as demonstrated by hypergeometric analysis. Strikingly, nonresponders exhibited very little change in any modules, despite exposure to TNF blockade. These patterns of change were highly consistent across all 3 cohorts, indicating that immunologic changes after TNF treatment are specific to the combination of both drug exposure and responder status. In contrast, modular patterns of gene expression did not exhibit consistent differences between responders and nonresponders at baseline in the 3 study cohorts. CONCLUSION: These data provide evidence that using gene expression modules related to inflammatory disease may provide a valuable method for objective monitoring of the response of RA patients who are treated with TNF inhibitors.

Sex stratification of an inflammatory bowel disease genome search shows male-specific linkage to the HLA region of chromosome 6.

Inflammatory bowel disease (IBD) is a multifactorial disorder, with both genetic and environmental factors contributing to the two clinical phenotypes of Crohn's disease (CD) and ulcerative colitis (UC). The underlying genetic model is thought to involve multiple genes with complex interactions between disease loci, and the NOD2 gene on chromosome 16 has recently been identified as a CD susceptibility locus. Several genome-wide linkage studies have identified candidate regions, but there has been little replication across studies. Here we investigate the role of sex-specific loci in susceptibility to IBD. Linkage data from our previously reported genome search and follow-up study were stratified by the sex of the affected sib pair. Non-parametric linkage analysis was performed using Genehunter Plus. Simulation studies were used to assess the significance of differences in LOD scores between male and female families for each chromosome. Several regions of sex-specific linkage were identified, including existing and novel candidate loci. The major histocompatibility region on chromosome 6p, referred to as IBD3, showed evidence of male-specific linkage with a maximum LOD score of 5.9 in both CD and UC male-affected families. Regions on chromosomes 11, 14 and 18 showed strong evidence of linkage in male-affected families but not in female-affected families. No evidence of sex-specific linkage was found in the IBD1 or IBD2 candidate regions of chromosomes 16 and 12. The existence of sex-specific linkage is further evidence of the complex mechanisms involved in IBD and will facilitate future studies to identify susceptibility genes.

Ethnic differences in cardiac potassium channel variants: implications for genetic susceptibility to sudden cardiac death and genetic testing for congenital long QT syndrome.

OBJECTIVE: To determine the spectrum, frequency, and ethnic-specificity of channel variants in the potassium channel genes implicated in congenital long QT syndrome (LQTS) among healthy subjects. SUBJECTS AND METHODS: Genomic DNA from 744 apparently healthy individuals-305 black, 187 white, 134 Asian, and 118 Hispanic--was subject to a comprehensive mutational analysis of the 4 LQTS-causing potassium channel genes: KCNQ1 (LQT1), KCNH2 (LQT2), KCNE1 (LQT5), and KCNE2 (LQT6). RESULTS: Overall, 49 distinct amino acid-altering variants (36 novel) were identified: KCNQ1 (n = 16), KCNH2 (n = 25),KCNE1 (n = 5), and KCNE2 (n = 3). More than half of these variants (26/49) were found exclusively in black subjects. The known K897T-HERG and the G38S-min K common polymorphisms were identified in all 4 ethnic groups. Excluding these 2 common polymorphisms, 25% of black subjects had at least 1 nonsynonymous potassium channel variant compared with 14% of white subjects (P < .01). CONCLUSIONS: To our knowledge, this study represents the first comprehensive determination of the frequency and spectrum of cardiac channel variants found among healthy subjects from 4 major ethnic groups. Defining the population burden of genetic variants in these critical cardiac ion channels is crucial for proper interpretation of genetic test results of individuals at risk for LQTS. This compendium provides a resource for epidemiological and functional investigation of variant effects on the repolarization properties of cardiac tissues, including susceptibility to lethal cardiac arrhythmias.

Spectrum and prevalence of cardiac sodium channel variants among black, white, Asian, and Hispanic individuals: implications for arrhythmogenic susceptibility and Brugada/long QT syndrome genetic testing.

OBJECTIVES: The purpose of this study was to determine the prevalence and spectrum of nonsynonymous polymorphisms (amino acid variants) in the cardiac sodium channel among healthy subjects. BACKGROUND: Pathogenic mutations in the cardiac sodium channel gene, SCN5A, cause approximately 15 to 20% of Brugada syndrome (BrS1), 5 to 10% of long QT syndrome (LQT3), and 2 to 5% of sudden infant death syndrome. METHODS: Using single-stranded conformation polymorphism, denaturing high-performance liquid chromatography, and/or direct DNA sequencing, mutational analysis of the protein-encoding exons of SCN5A was performed on 829 unrelated, anonymous healthy subjects: 319 black, 295 white, 112 Asian, and 103 Hispanic. RESULTS: In addition to the four known common polymorphisms (R34C, H558R, S1103Y, and R1193Q), four relatively ethnic-specific polymorphisms were identified: R481W, S524Y, P1090L, and V1951L. Overall, 39 distinct missense variants (28 novel) were elucidated. Nineteen variants (49%) were found only in the black cohort. Only seven variants (18%) localized to transmembrane-spanning domains. Four variants (F1293S, R1512W, and V1951L cited previously as BrS1-causing mutations and S1787N previously published as a possible LQT3-causing mutation) were identified in this healthy cohort. CONCLUSIONS: This study provides the first comprehensive determination of the prevalence and spectrum of cardiac sodium channel variants in healthy subjects from four distinct ethnic groups. This compendium of SCN5A variants is critical for proper interpretation of SCN5A genetic testing and provides an essential hit list of targets for future functional studies to determine whether or not any of these variants mediate genetic susceptibility for arrhythmias in the setting of either drugs or disease.

Diagnostics and personalized medicine: observations from the World Companion Diagnostics Summit.

On 1-2 December 2010 worldwide leaders in the field of co-diagnostics gathered in Boston, MA, USA, to discuss the state of predictive medicine and the complexities of biomarker implementation and partnering. Flanking the primary conference were three workshops where participants could increase their understanding of companion diagnostic commercialization, best practices from well-known case studies and regulatory guidance and practice. The primary conference reviewed here consisted of 18 lectures and one round table discussion. Topics for the summit included biomarker discovery and validation, partnering between pharmaceutical development and diagnostics companies, commercialization models and evolving regulatory perspectives and guidance. While it is clear that there are many difficulties remaining in making personalized medicine a broad-based reality for all therapeutic areas, it is equally clear that much progress has already been made and that the conference participants are optimistic and dedicated to fulfilling the personalized medicine promise of bringing the right drug to the right patient at the right time.

Stratification by CARD15 variant genotype in a genome-wide search for inflammatory bowel disease susceptibility loci.

Previously we have conducted a genome-wide search for inflammatory bowel disease susceptibility loci in a large European cohort. Results from this study demonstrated suggestive evidence of linkage to loci at chromosomes 1q, 6p, and 10p and replicated linkages on chromosomes 12 and 16. Recently, NOD2/CARD15 on chromosome 16q12 has been found to be strongly associated with Crohn's disease. In order to determine if there are other loci in the genome that interact with the three associated functional variants in CARD15 (R702W, G908R, 1007fs), we have stratified our large inflammatory bowel disease genome scan cohort by dividing pedigrees into two groups stratified by CARD15 variant genotype. The two pedigree groups were analysed using non-parametric allele sharing methods. The group of pedigrees that contained one of the three CARD15 variants had two suggestive linkage results occurring in 6p (lod = 3.06 at D6S197, IBD phenotype) and 10p (lod=2.29 at D10S197, CD phenotype). In addition, at 16q12 where CARD15 is located, the original genome scan had a peak lod score of 2.18 at D16S415 (CD phenotype). The stratified pedigree cohort containing one of three CARD15 variants had a peak lod score of 0.90 at D16S415 (CD phenotype), accounting for approximately less than half of the genetic evidence for linkage at this locus. This result is in agreement with the existence of a substantial number of private variants at the NOD2/CARD15 locus. Interaction with NOD2/CARD15 needs to be considered in future gene identification efforts on chromosomes 6 and 10.

Sudden infant death syndrome: association with a promoter polymorphism of the serotonin transporter gene.

Serotonergic receptor binding in the arcuate nucleus, n. raphé obscurus, and other medullary regions is decreased in sudden infant death syndrome (SIDS) cases. Further, a variable tandem repeat sequence polymorphism in the promoter region of the serotonin transporter protein (5-HTT) gene has recently been associated with risk of SIDS in a Japanese cohort. This polymorphism differentially regulates 5-HTT expression, with the long allele (L), the SIDS-associated allele, being a more effective promoter than the short allele (S). We therefore investigated the 5-HTT promoter polymorphism in a cohort of 87 SIDS cases (43 African American and 44 Caucasian) and gender/ethnicity-matched controls. Significant positive associations were found between SIDS and the 5-HTT genotype distribution (P = 0.022), specifically with the L/L genotype (P = 0.048), and between SIDS and the 5-HTT L allele (P = 0.005). There was also a significant negative association between SIDS and the S/S genotype (P = 0.011). The comparisons were repeated in the African American and Caucasian subgroups. The data patterns were consistent in the subgroups, i.e., the L/L genotype and L allele were increased in the cases, but not all subgroup comparisons were statistically significant. These results indicate a relationship between SIDS and the L allele of the 5-HTT gene in African Americans and Caucasians, and if confirmed, will provide an important tool for identifying at-risk individuals and estimating the risk of recurrence.

Idiopathic congenital central hypoventilation syndrome: analysis of genes pertinent to early autonomic nervous system embryologic development and identification of mutations in PHOX2b.

Idiopathic congenital central hypoventilation syndrome (CCHS) has been linked to autonomic nervous system dysregulation and/or dysfunction (ANSD) since it was first described in 1970. A genetic basis of CCHS has been proposed because of the reports of four families with two affected children, because of mother-child transmission, and because of a recent report of a polyalanine expansion mutation in PHOX2b in a subset of CCHS subjects. We, therefore, studied genes pertinent to early embryologic development of the ANS including mammalian achaete-scute homolog-1 (MASH1), bone morphogenic protein-2 (BMP2), engrailed-1 (EN1), TLX3, endothelin converting enzyme-1 (ECE1), endothelin-1 (EDN1), PHOX2a, and PHOX2b in 67 probands with CCHS, and gender- and ethnicity-matched controls. No disease-defining mutations were identified in MASH1, BMP2, EN1, TLX3, ECE1, EDN1, or PHOX2a. The 65/67 CCHS probands (97%) were found to be heterozygous for the exon 3 polyalanine expansion mutation identified previously in PHOX2b. Further, there was an association between repeat mutation length and severity of the CCHS/ANSD phenotype. Of the two probands who did not carry the expansion mutation, one had a nonsense mutation in exon 3 which truncated the protein and the other had no mutation in PHOX2b but had a previously reported EDN3 frameshift point mutation. The polyalanine expansion mutation was not found in any of 67 matched controls. Of 54 available families (including 97 unaffected parents), whose child carried the PHOX2b mutation, 4 parents demonstrated mosaicism for an expansion mutation identical to that seen in the CCHS cases, suggesting that not all mutations in affected probands with unaffected parents are de novo. We also studied four women with CCHS who were heterozygous for the PHOX2b mutation, each with one child. Three of the four children were also affected and had the same mutation, demonstrating autosomal dominant inheritance of the mutation. Assay of the PHOX2b polyalanine repeat mutation represents a highly sensitive and specific technique for confirming the diagnosis of CCHS. Identification of the CCHS mutation will lead to clarification of the phenotype, allow for prenatal diagnosis for parents of CCHS probands and adults with CCHS in future pregnancies, and potentially direct intervention strategies for the treatment of CCHS.

Sudden infant death syndrome: case-control frequency differences at genes pertinent to early autonomic nervous system embryologic development.

We have previously identified polymorphisms in the serotonin transporter gene promoter region and in intron 2 that were more common among sudden infant death syndrome (SIDS) cases compared with control subjects. To elucidate further the genetic profile that might increase an infant's vulnerability to SIDS, we focused on the recognized relationship between autonomic nervous system (ANS) dysregulation and SIDS. We therefore studied genes pertinent to early embryologic development of the ANS, including MASH1, BMP2, PHOX2a, PHOX2b, RET, ECE1, EDN1, TLX3, and EN1 in 92 probands with SIDS and 92 gender- and ethnicity-matched control subjects. Eleven protein-changing rare mutations were identified in 14 of 92 SIDS cases among the PHOX2a, RET, ECE1, TLX3, and EN1 genes. Only 1 of these mutations (TLX3) was identified in 2 of 92 control subjects. Black infants accounted for 10 of these mutations in SIDS cases and 2 control subjects. Four protein-changing common polymorphisms were identified in BMP2, RET, ECE1, and EDN1, but the allele frequency did not differ between SIDS cases and control subjects. However, among SIDS cases, the allele frequency for the BMP2 common polymorphism demonstrated ethnic differences; among control subjects, the allele frequency for the BMP2 and the ECE1 common polymorphisms also demonstrated ethnic differences. These data represent further refinement of the genetic profile that might place an infant at risk for SIDS.