Mangalathu S. Rajeevan

Assuring the quality of next-generation sequencing in clinical laboratory practice

Amy S Gargis, Centers for Disease Control and Prevention Lisa Kalman, Centers for Disease Control and Prevention Meredith W Berry, SeqWright Inc David P Bick, Medical College of Wisconsin David P Dimmock, Medical College of Wisconsin Tina Hambuch, Illumina Clinical Services Fei Lu, SeqWright Inc Elaine Lyon, University of Utah Karl V Voelkerding, University of Utah Barbara Zehnbauer, Emory University

Utility of the blood for gene expression profiling and biomarker discovery in chronic fatigue syndrome

Chronic fatigue syndrome (CFS) is a debilitating illness lacking consistent anatomic lesions and eluding conventional laboratory diagnosis. Demonstration of the utility of the blood for gene expression profiling and biomarker discovery would have implications into the pathophysiology of CFS. The objective of this study was to determine if gene expression profiles of peripheral blood mononuclear cells (PMBCs) could distinguish between subjects with CFS and healthy controls. Total RNA from PBMCs of five CFS cases and seventeen controls was labeled and hybridized to 1764 genes on filter arrays. Gene intensity values were analyzed by various classification algorithms and nonparametric statistical methods. The classification algorithms grouped the majority of the CFS cases together, and distinguished them from the healthy controls. Eight genes were differentially expressed in both an age-matched case-control analysis and when comparing all CFS cases to all controls. Several of the diffrentially expressed genes are associated with immunologic functions (e.g., CMRF35 antigen, IL-8, HD protein) and implicate immune dysfunction in the pathophysiology of CFS. These results successfully demonstrate the utility of the blood for gene expression profiling to distinguish subjects with CFS from healthy controls and for identifying genes that could serve as CFS biomarkers.

Glucocorticoid receptor polymorphisms and haplotypes associated with chronic fatigue syndrome

Chronic fatigue syndrome (CFS) is a significant public health problem of unknown etiology, the pathophysiology has not been elucidated, and there are no characteristic physical signs or laboratory abnormalities. Some studies have indicated an association of CFS with deregulation of immune functions and hypothalamic–pituitary–adrenal (HPA) axis activity. In this study, we examined the association of sequence variations in the glucocorticoid receptor gene (NR3C1) with CFS because NR3C1 is a major effector of the HPA axis. There were 137 study participants (40 with CFS, 55 with insufficient symptoms or fatigue, termed as ISF, and 42 non‐fatigued controls) who were clinically evaluated and identified from the general population of Wichita, KS. Nine single nucleotide polymorphisms (SNPs) in NR3C1 were tested for association of polymorphisms and haplotypes with CFS. We observed an association of multiple SNPs with chronic fatigue compared to non‐fatigued (NF) subjects (P < 0.05) and found similar associations with quantitative assessments of functional impairment (by the SF‐36), with fatigue (by the Multidimensional Fatigue Inventory) and with symptoms (assessed by the Centers for Disease Control Symptom Inventory). Subjects homozygous for the major allele of all associated SNPs were at increased risk for CFS with odds ratios ranging from 2.61 (CI 1.05–6.45) to 3.00 (CI 1.12–8.05). Five SNPs, covering a region of approximately 80 kb, demonstrated high linkage disequilibrium (LD) in CFS, but LD gradually declined in ISF to NF subjects. Furthermore, haplotype analysis of the region in LD identified two associated haplotypes with opposite alleles: one protective and the other conferring risk of CFS. These results demonstrate NR3C1 as a potential mediator of chronic fatigue, and implicate variations in the 5′ region of NR3C1 as a possible mechanism through which the alterations in HPA axis regulation and behavioural characteristics of CFS may manifest.

Polymorphisms in genes regulating the HPA axis associated with empirically delineated classes of unexplained chronic fatigue

Chronic fatigue syndrome (CFS) is characterized by persistent or relapsing fatigue that is not alleviated by rest, causes substantial reduction in activities and is accompanied by a variety of symptoms. Its unknown etiology may reflect that CFS is heterogeneous. Latent class analyses of symptoms and physiological systems were used to delineate subgroups within a population-based sample of fatigued and nonfatigued subjects [1] . This study examined whether genetic differences underlie the individual subgroups of the latent class solution. Polymorphisms in 11 candidate genes related to both hypothalamic-pituitary-adrenal (HPA) axis function and mood-related neurotransmitter systems were evaluated by comparing each of the five ill classes (Class 1, n = 33; Class 3, n = 22; Class 4, n = 22; Class 5, n = 17; Class 6, n = 11) of fatigued subjects with subjects defined as well (Class 2, n = 35). Of the five classes of subjects with unexplained fatigue, three classes were distinguished by gene polymorphsims involved in either HPA axis function or neurotransmitter systems, including proopiomelanocortin (POMC), nuclear receptor subfamily 3, group C, member 1 (NR3C1), monoamine oxidase A (MAOA), monoamine oxidase B (MAOB), and tryptophan hydroxylase 2 (TPH2). These data support the hypothesis that medically unexplained chronic fatigue is heterogeneous and presents preliminary evidence of the genetic mechanisms underlying some of the putative conditions.

Epidemiologic and viral factors associated with cervical neoplasia in HPV-16-positive women

While infection with high‐risk HPV is the most important risk factor for cervical cancer, HPV alone is insufficient. Our purpose was to identify viral and epidemiologic factors associated with cervical disease in HPV‐16 DNA‐positive women referred to colposcopy. We used a standardized interview to collect epidemiologic data from consenting women. Total nucleic acids from exfoliated cervical cells were used for all viral assays (HPV detection and typing using L1 consensus PCR with line probe hybridization, variant classification by sequencing, viral load and transcript copy determination by quantitative PCR and transcript pattern by nested RT‐PCR). Cervical disease was based on colposcopic biopsy. Logistic regression was used to calculate ORs with 95% CIs. There were 115 HPV‐16 positive women among 839 enrollees. By univariate analyses, age >25 years (OR = 3.05, 95% CI 1.20–7.76), smoking (OR = 3.0, 95% CI 1.19–7.56), high viral load (OR = 5.27, 95% CI 2.05–13.60), detection of both E6 and E6*I transcripts (OR = 10.0, 95% CI 2.1–47.58) and high transcript copies (OR = 5.56, 95% CI 2.05–13.60) were significant risk factors for CIN III with reference to No CIN/CIN I. Less than a third of the women (31.5%) had prototype HPV‐16 detected, and variants showed no association with disease, viral load or transcription. Viral DNA and transcript copies were highly correlated, and the ratio of transcript copies to DNA copies was not changed with disease status. While viral load, transcript copies and transcript pattern were statistically associated with CIN III, none of these measures effectively discriminated between HPV‐16 women with disease requiring treatment and those who could be followed. Cellular proliferation and differentiation pathways affected by HPV should be investigated as biomarkers for cervical cancer screening.

Reproducibility of Alternative Probe Synthesis Approaches for Gene Expression Profiling with Arrays

Before gene expression profiling with microarray technology can be transferred to the diagnostic setting, we must have alternative approaches for synthesizing probe from limited RNA samples, and we must understand the limits of reproducibility in interpreting gene expression results. The current gold standard of probes for use with both microarrays and high-density filter arrays are synthesized from 1 microg of purified poly(A)+ RNA. We evaluated two approaches for synthesizing cDNA probes from total RNA with subsequent hybridization to high-density filter arrays: 1) reverse transcription (RT) of 5 microg total RNA and 2) RT-polymerase chain reaction (RT-PCR) of 1 microg total RNA, using the SMART system. The reproducibility of these two approaches was compared to the current gold standard. All three methods were highly reproducible. Triplicate experiments resulted in the following concordance correlation coefficients to evaluate reproducibility: 0.88 for the gold standard, 0.86 for cDNA probe synthesized by RT from total RNA, and 0.96 for the SMART cDNA probe synthesized from total RNA. We also compared the expression profile of 588 genes for the total RNA methods to that obtained with the gold standard. Of 150 positive genes detected by the gold standard, 97 (65%) were detected by cDNA probe synthesized by RT of total RNA, and 122 (81%) were detected by the SMART cDNA probe. We conclude that SMART cDNA probe produces highly reproducible results and yields gene expression profiles that represent the majority of transcripts detected with the gold standard.

Functional Genomics of Serotonin Receptor 2A (HTR2A): Interaction of Polymorphism, Methylation, Expression and Disease Association

Serotonergic neurotransmission plays a key role in the pathophysiology of neuropsychiatric illnesses. The functional significance of a promoter polymorphism, −1438G/A (rs6311), in one of the major genes of this system (serotonin receptor 2A, HTR2A) remains poorly understood in the context of epigenetic factors, transcription factors and endocrine influences. We used functional and structural equation modeling (SEM) approaches to assess the contributions of the polymorphism (rs6311), DNA methylation and clinical variables to HTR2A expression in chronic fatigue syndrome (CFS) subjects from a population-based study. HTR2A was up-regulated in CFS through allele-specific expression modulated by transcription factors at critical sites in its promoter: an E47 binding site at position −1,438, (created by the A-allele of rs6311 polymorphism), a glucocorticoid receptor (GR) binding site encompassing a CpG at position −1,420, and Sp1 binding at CpG methylation site −1,224. Methylation at −1,420 was strongly correlated with methylation at −1,439, a CpG site that is dependent upon the G-allele of rs6311 at position −1,438. SEM revealed a strong negative interaction between E47 and GR binding (in conjunction with cortisol level) on HTR2A expression. This study suggests that the promoter polymorphism (rs6311) can affect both transcription factor binding and promoter methylation, and this along with an individual’s stress response can impact the rate of HTR2A transcription in a genotype and methylation-dependent manner. This study can serve as an example for deciphering the molecular determinants of transcriptional regulation of major genes of medical importance by integrating functional genomics and SEM approaches. Confirmation in an independent study population is required.

Genetic evaluation of the serotonergic system in chronic fatigue syndrome.

Chronic fatigue syndrome (CFS) is a debilitating disorder of unknown etiology with no known lesions, diagnostic markers or therapeutic intervention. The pathophysiology of CFS remains elusive, although abnormalities in the central nervous system (CNS) have been implicated, particularly hyperactivity of the serotonergic (5-hydroxytryptamine; 5-HT) system and hypoactivity of the hypothalamic-pituitary-adrenal (HPA) axis. Since alterations in 5-HT signaling can lead to physiologic and behavioral changes, a genetic evaluation of the 5-HT system was undertaken to identify serotonergic markers associated with CFS and potential mechanisms for CNS abnormality. A total of 77 polymorphisms in genes related to serotonin synthesis (TPH2), signaling (HTR1A, HTR1E, HTR2A, HTR2B, HTR2C, HTR3A, HTR3B, HTR4, HTR5A, HTR6, and HTR7), transport (SLC6A4), and catabolism (MAOA) were examined in 137 clinically evaluated subjects (40 CFS, 55 with insufficient fatigue, and 42 non-fatigued, NF, controls) derived from a population-based CFS surveillance study in Wichita, Kansas. Of the polymorphisms examined, three markers (-1438G/A, C102T, and rs1923884) all located in the 5-HT receptor subtype HTR2A were associated with CFS when compared to NF controls. Additionally, consistent associations were observed between HTR2A variants and quantitative measures of disability and fatigue in all subjects. The most compelling of these associations was with the A allele of -1438G/A (rs6311) which is suggested to have increased promoter activity in functional studies. Further, in silico analysis revealed that the -1438 A allele creates a consensus binding site for Th1/E47, a transcription factor implicated in the development of the nervous system. Electrophoretic mobility shift assay supports allele-specific binding of E47 to the A allele but not the G allele at this locus. These data indicate that sequence variation in HTR2A, potentially resulting in its enhanced activity, may be involved in the pathophysiology of CFS.

Human papillomavirus prevalence in invasive anal cancers in the United States before vaccine introduction.

Objective This study aimed to conduct a representative survey of human papillomavirus (HPV) prevalence and its genotype distribution in invasive anal cancer specimens in the United States. Materials and Methods Population-based archival anal cancer specimens were identified from Florida, Kentucky, Louisiana, and Michigan cancer registries and Surveillance, Epidemiology, and End Results (SEER) tissue repositories in Hawaii, Iowa, and Los Angeles. Sections from 1 representative block per case were used for DNA extraction. All extracts were assayed first by linear array and retested with INNO-LiPA if inadequate or HPV negative. Results Among 146 unique invasive anal cancer cases, 93 (63.7%) were from women, and 53 (36.3%) were from men. Human papillomavirus (any type) was detected in 133 cases (91.1%) and 129 (88.4%) contained at least 1 high risk-type, most (80.1%) as a single genotype. Human papillomavirus type 16 had the highest prevalence (113 cases, 77.4%); HPV types 6, 11, 18, and 33 were also found multiple times. Among HPV-16–positive cases, 37% were identified as prototype variant Ep, and 63% were nonprototypes: 33% Em, 12% E-G131G, 5% Af1, 4% AA/NA-1, 3% E-C109G, 3% E-G131T, 2% As, and 1% Af2. No significant differences in the distributions of HPV (any), high-risk types, or HPV-16/18 were seen between sex, race, or age group. Conclusions The establishment of prevaccine HPV prevalence in the United States is critical to the surveillance of vaccine efficacy. Almost 80% of anal cancers were positive for the vaccine types HPV-16 or HPV-18, and in 70%, these were the only types detected, suggesting that a high proportion might be preventable by current vaccines.

Transcriptional control of complement activation in an exercise model of chronic fatigue syndrome.

Complement activation resulting in significant increases of C4a split product may be a marker of postexertional malaise in individuals with chronic fatigue syndrome (CFS). This study focused on identification of the transcriptional control that may contribute to the increased C4a in CFS subjects after exercise. We used quantitative reverse-transcription polymerase chain reaction to evaluate differential expression of genes in the classical and lectin pathways in peripheral blood mononuclear cells (PBMCs). Calibrated expression values were normalized to the internal reference gene peptidylpropyl isomerase B (PPIB), the external reference gene ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL), or the geometric mean (GM) of the genes ribosomal protein, large, P0 (RPLP0) and phosphoglycerate kinase 1 (PGK1). All nine genes tested, except mannose-binding lectin 2 (MBL2), were expressed in PBMCs. At 1 hour postexercise, C4, mannan-binding lectin serine protease 2 (MASP2) and ficolin 1 (FCN1) transcripts were detected at higher levels (≥2-fold) in at least 50% (4 of 8) of CFS subjects and were detected in 88% (7 of 8) CFS subjects when subjects with overexpression of either C4 or MASP2were combined. Only an increase in the MASP2 transcript was statistically significant (PPIB, P = 0.001; GM, P = 0.047; rbcL, P = 0.045). This result may be due to the significant but transient downregulation of MASP2 in control subjects (PPIB, P = 0.023; rbcL, P = 0.027). By 6 hours postexercise, MASP2 expression was similar in both groups. In conclusion, lectin pathway responded to exercise differentially in CFS than in control subjects. MASP2 down-regulation may act as an antiinflammatory acute-phase response in healthy subjects, whereas its elevated level may account for increased C4a and inflammation-mediated postexertional malaise in CFS subjects.