Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes

Genome-wide association studies have successfully discovered thousands of common variants associated with human diseases and traits, but the landscape of rare variations in human disease has not been explored at scale. Exome-sequencing studies of population biobanks provide an opportunity to systematically evaluate the impact of rare coding variations across a wide range of phenotypes to discover genes and allelic series relevant to human health and disease. Here, we present results from systematic association analyses of 4,529 phenotypes using single-variant and gene tests of 394,841 individuals in the UK Biobank with exome-sequence data. We find that the discovery of genetic associations is tightly linked to frequency and is correlated with metrics of deleteriousness and natural selection. We highlight biological findings elucidated by these data and release the dataset as a public resource alongside the Genebass browser for rapidly exploring rare-variant association results.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

SAIGEgds-an efficient statistical tool for large-scale PheWAS with mixed models

SUMMARY

Phenome-wide association studies (PheWASs) are known to be a powerful tool in discovery and replication of genetic association studies. To reduce the computational burden of PheWAS in the large cohorts, such as the UK Biobank, the SAIGE method has been proposed to control for case-control imbalance and sample relatedness in a tractable manner. However, SAIGE is still computationally intensive when deployed in analyzing the associations of thousands of ICD10-coded phenotypes with whole-genome imputed genotype data. Here, we present a new high-performance statistical R package (SAIGEgds) for large-scale PheWAS using generalized linear mixed models. The package implements the SAIGE method in optimized C++ codes, taking advantage of sparse genotype dosages and integrating the efficient genomic data structure file format. Benchmarks using the UK Biobank White British genotype data (N ≈ 430 K) with coronary heart disease and simulated cases show that the implementation in SAIGEgds is 5-6 times faster than the SAIGE R package. When used in conjunction with high-performance computing clusters, SAIGEgds provides an efficient analysis pipeline for biobank-scale PheWAS.

AVAILABILITY AND IMPLEMENTATION

https://bioconductor.org/packages/SAIGEgds; vignettes included.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A catalog of associations between rare coding variants and COVID-19 outcomes

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) causes coronavirus disease-19 (COVID-19), a respiratory illness that can result in hospitalization or death. We investigated associations between rare genetic variants and seven COVID-19 outcomes in 543,213 individuals, including 8,248 with COVID-19. After accounting for multiple testing, we did not identify any clear associations with rare variants either exome-wide or when specifically focusing on (i) 14 interferon pathway genes in which rare deleterious variants have been reported in severe COVID-19 patients; (ii) 167 genes located in COVID-19 GWAS risk loci; or (iii) 32 additional genes of immunologic relevance and/or therapeutic potential. Our analyses indicate there are no significant associations with rare protein-coding variants with detectable effect sizes at our current sample sizes. Analyses will be updated as additional data become available, with results publicly browsable at https://rgc-covid19.regeneron.com.

Checkpoint Inhibitor Colitis Masked by Multi-drug Effect

Casestudy: The utilization of checkpoint inhibitors such as programmed cell death protein 1 (PD-1/PD-L1) inhibitors (nivolumab) and cytotoxic T-lymphocyte antigen 4 inhibitors (ipilimumab) for treatment of certain malignancies has steadily gained popularity. Medication related colitis is uncommon, with a reported incidence of 1–9% depending on the checkpoint inhibitor used, and the histologic features have been characterized in recent literature. Because of the immunomodulating effect of these drugs, infectious colitis is in the differential diagnosis of enteritis. Multi-drug therapy in many of these patients further complicates identification of the culprit drug. We present the case of a 63-year-old male with metastatic renal cell carcinoma being treated with both nivolumab and ipilimumab who presented with acute on chronic non-bloody diarrhea. His clinical course was complicated by hypotension, acidosis and coagulopathy. The clinical differential for his colitis was cytomegalovirus infection versus a checkpoint inhibitor colitis. Colonoscopy revealed continuous circumferential loss of vascularity and diffuse erythema throughout the colon. Histology showed acute colitis with prominent apoptosis, cryptitis, crypt abscesses, and rare ringed mitotic figures, but without architectural distortion. Occasional smooth purple crystals consistent with pill material were present in the mucosa, but without significant tissue reaction. No pathogenic organisms were identified, and a cytomegalovirus immunostain was negative. These histologic findings in concert with the clinical history are consistent with checkpoint inhibitor colitis and multi-drug effect. A review of the patient’s chart showed cholestyramine was added to the patient’s regimen during hospitalization, which was consistent with the morphologic appearance of the crystals. Given the acute complications of checkpoint inhibitor induced enterocolitis and potential for increased morbidity (rare cases of bowel perforation and subsequent resection), accurate diagnosis is imperative. Management of checkpoint inhibitor associated colitis ranges from initiation of immunosuppression to checkpoint inhibitor cessation. When these findings are masked by multi-drug effect, accurate diagnosis can be difficult.

Are drug targets with genetic support twice as likely to be approved? Revised estimates of the impact of genetic support for drug mechanisms on the probability of drug approval

Despite strong vetting for disease activity, only 10% of candidate new molecular entities in early stage clinical trials are eventually approved. Analyzing historical pipeline data, Nelson et al. 2015 (Nat. Genet.) concluded pipeline drug targets with human genetic evidence of disease association are twice as likely to lead to approved drugs. Taking advantage of recent clinical development advances and rapid growth in GWAS datasets, we extend the original work using updated data, test whether genetic evidence predicts future successes and introduce statistical models adjusting for target and indication-level properties. Our work confirms drugs with genetically supported targets were more likely to be successful in Phases II and III. When causal genes are clear (Mendelian traits and GWAS associations linked to coding variants), we find the use of human genetic evidence increases approval by greater than two-fold, and, for Mendelian associations, the positive association holds prospectively. Our findings suggest investments into genomics and genetics are likely to be beneficial to companies deploying this strategy.

The growth of human genetics resources has the potential to help us develop better drugs. By looking at whether and how historical drug approvals could have been predicted from our current knowledge of human genetics, we can validate this approach and assess which types of genetic evidence are most likely to be useful in guiding drug discovery. Validation is important because we are often uncertain about the biological mechanisms behind genetic variants linked to disease. Most associated variants do not occur within protein-coding regions of the genome, and it is difficult to tell which of many nearby genes is contributing to disease risk. In this paper, we confirm previous correlations between genetic evidence and historical drug approvals. We find genetic evidence from severe genetic disorders and from genetic variants that alter protein sequence is more strongly associated with historical approvals. We offer statistical approaches for prioritizing new drug candidates based on whether their mechanisms are supported by human genetic evidence.

Mineralocorticoid receptor antagonism reverses diabetes-related coronary vasodilator dysfunction: A unique vascular transcriptomic signature

Coronary microvascular dysfunction predicts and may be a proximate cause of cardiac dysfunction and mortality in diabetes; however, few effective treatments exist for these conditions. We recently demonstrated that mineralocorticoid receptor (MR) antagonism reversed cardiovascular dysfunction in early-stage obesity/insulin resistance. The mechanisms underlying this benefit of MR antagonism and its relevance in the setting of long-term obesity complications like diabetes; however, remain unclear. Thus, the present study evaluated the impact of MR antagonism on diabetes-related coronary dysfunction and defines the MR-dependent vascular transcriptome in the Otsuka Long-Evans Tokushima Fatty (OLETF) rat recapitulating later stages of human diabetes. OLETF rats were treated with spironolactone (Sp) and compared to untreated OLETF and lean Long-Evans Tokushima Otsuka rats. Sp treatment attenuated diabetes-associated adipose and cardiac inflammation/fibrosis and improved coronary endothelium-dependent vasodilation but did not alter enhanced coronary vasoconstriction, blood pressure, or metabolic parameters in OLETF rats. Further mechanistic studies using RNA deep sequencing of OLETF rat aortas revealed 157 differentially expressed genes following Sp including upregulation of genes involved in the molecular regulation of nitric oxide bioavailability (Hsp90ab1, Ahsa1, Ahsa2) as well as novel changes in α_1D adrenergic receptors (Adra1d), cyclooxygenase-2 (Ptgs2), and modulatory factors of these pathways (Ackr3, Acsl4). Further, Ingenuity Pathway Analysis predicted inhibition of upstream inflammatory regulators by Sp and inhibition of 'migration of endothelial cells', 'differentiation of smooth muscle', and 'angiogenesis' biological functions by Sp in diabetes. Thus, this study is the first to define the MR-dependent vascular transcriptome underlying treatment of diabetes-related coronary microvascular dysfunction by Sp.

Epigenetic analysis of the IFNλ3 gene identifies a novel marker for response to therapy in HCV-infected subjects

Chronic hepatitis C virus (HCV) infection is characterized by high interindividual variability in response to pegylated interferon and ribavirin. A genetic polymorphism on chromosome 19 (rs12979860) upstream of interferon-λ3 (IFNλ3) is associated with a twofold change in sustained virologic response rate after 48 weeks of treatment with pegylated interferon/ribavirin in HCV genotype 1 (GT1) treatment-naïve patients. We conducted epigenetic analysis on the IFNλ3 promoter to investigate whether DNA methylation is associated with response to HCV therapy. DNA samples from HCV GT1-infected subjects receiving an interferon-free paritaprevir-containing combination regimen (N=540) and from HCV-uninfected, healthy controls (N=124) were analysed for IFNλ3 methylation levels. Methylation was strongly associated with rs12979860 allele status whether adjusting for HCV status (r=65.0%, 95% CI: [60.2%, 69.5%]), or not (r=64.4%), both with P<2.2×10-16 . In HCV GT1-infected subjects, C/C genotypes had significantly lower methylation levels relative to C/T or T/T genotypes (P<1×10-14 ), with each T allele resulting in a nine-unit increase in mean methylation level. Methylation levels did not correlate with response in subjects treated for 12 or 24 weeks. However, non-C/C subjects with higher methylation levels were more likely to relapse when treatment duration was 8 weeks. This analysis suggests that methylation status of the IFNλ3 promoter region may be a useful parameter that identifies patients more likely to relapse following HCV therapy; however, continuing therapy for a sufficient duration can overcome this difference. These findings may provide mechanistic insight into the role of IFNλ3 genetic variants in HCV treatment response.

Seasonal Population Dynamics of Megacopta cribraria (Hemiptera: Plataspidae) in Kudzu and Soybean, and Implication for Insecticidal Management in Soybean

Megacopta cribraria (F.), an invasive species introduced from Asia in 2009, is now prolific in the southeastern United States. Megacopta cribraria develops primarily on kudzu and soybean completing two generations. It is not well understood how this economic pest is affected by changes in geographic distribution in the United States or how population levels have changed since its establishment. The effect of insecticide application timing on field populations of M. cribraria is not well documented. These studies seek to understand how population dynamics of M. cribraria vary with geographic regions in Georgia. Effect of application timing on populations throughout the growing season was also examined. Weekly from 2012 to 2013, all life stages were enumerated from kudzu and soybean environments at several locations throughout Georgia from sweeps samples and flight intercept captures. Coordinates were recorded for locations, and classified as belonging to the Piedmont or Coastal Plain region of the state. Single spray trials were conducted from 2011-2014, and applications were made to soybean at intervals throughout the season. From 2012 to 2015, two kudzu patches near Griffin, GA, were monitored to detect population changes. Differences in population dynamics from locations around the state were found, but no clear effect of latitude, longitude, or region was observed. Insecticide applications applied in July suppressed nymph populations significantly better than treatments made earlier or later. Megacopta cribraria populations declined in 2014 and 2015 compared with 2012 and 2013. These studies provide the critical information for M. cribraria management in soybean in the southeastern United States.

Transcriptomic effects of metformin in skeletal muscle arteries of obese insulin-resistant rats

We examined the effects of metformin, a commonly used antidiabetic drug, on gene expression in multiple arteries. Specifically, transcriptional profiles of feed arteries and second branch order arterioles in the soleus, gastrocnemius, and diaphragm muscles as well as aortic endothelial scrapes were examined from obese insulin-resistant Otsuka Long-Evans Tokushima Fatty rats treated with ( n = 9) or without ( n = 10) metformin from 20 to 32 weeks of age. Metformin-treated rats exhibited a reduction in body weight, adiposity, and HbA1c ( P < 0.05). The greatest number of differentially expressed genes (FDR < 15%) between those treated with and without metformin was found in the red gastrocnemius 2a arterioles (93 genes), followed by the diaphragm 2a arterioles (62 genes), and soleus 2a arterioles (15 genes). We also found that two genes were differentially expressed in aortic endothelial cells (LETMD1 and HMGCS2, both downregulated), one gene in the gastrocnemius feed artery (BLNK, downregulated), and no genes in the soleus and diaphragm feed arteries and white gastrocnemius 2a arterioles. No single gene was altered by metformin across all vessels examined. This study provides evidence that metformin treatment produces distinct gene expression effects throughout the arterial tree in a rat model of obesity and insulin resistance. Genes whose expression was modulated with metformin do not appear to have a clear connection with its known mechanisms of action. These findings support the notion that vascular gene regulation in response to oral pharmacological therapy, such as metformin, is vessel specific. Impact statement This study provides evidence that metformin treatment produces artery-specific gene expression effects. The genes whose expression was modulated with metformin do not appear to have a clear connection with its known mechanisms of action.

Inferring a role for methylation of intergenic DNA in the regulation of genes aberrantly expressed in precursor B-cell acute lymphoblastic leukemia

A complete understanding of the mechanisms involved in the development of pre-B ALL is lacking. In this study, we integrated DNA methylation data and gene expression data to elucidate the impact of aberrant intergenic DNA methylation on gene expression in pre-B ALL. We found a subset of differentially methylated intergenic loci that were associated with altered gene expression in pre-B ALL patients. Notably, 84% of these regions were also bound by transcription factors (TF) known to play roles in differentiation and B-cell development in a lymphoblastoid cell line. Further, an overall downregulation of eRNA transcripts was observed in pre-B ALL patients and these transcripts were associated with the downregulation of putative target genes involved in B-cell migration, proliferation, and apoptosis. The identification of novel putative regulatory regions highlights the significance of intergenic DNA sequences and may contribute to the identification of new therapeutic targets for the treatment of pre-B ALL.

Integrated methylome and transcriptome analysis reveals novel regulatory elements in pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common cancer diagnosed in children under the age of 15. In addition to genetic aberrations, epigenetic modifications such as DNA methylation are altered in cancer and impact gene expression. To identify epigenetic alterations in ALL, genome-wide methylation profiles were generated using the methylated CpG island recovery assay followed by next-generation sequencing. More than 25,000 differentially methylated regions (DMR) were observed in ALL patients with ∼90% present within intronic or intergenic regions. To determine the regulatory potential of the DMR, whole-transcriptome analysis was performed and integrated with methylation data. Aberrant promoter methylation was associated with the altered expression of genes involved in transcriptional regulation, apoptosis, and proliferation. Novel enhancer-like sequences were identified within intronic and intergenic DMR. Aberrant methylation in these regions was associated with the altered expression of neighboring genes involved in cell cycle processes, lymphocyte activation and apoptosis. These genes include potential epi-driver genes, such as SYNE1, PTPRS, PAWR, HDAC9, RGCC, MCOLN2, LYN, TRAF3, FLT1, and MELK, which may provide a selective advantage to leukemic cells. In addition, the differential expression of epigenetic modifier genes, pseudogenes, and non-coding RNAs was also observed accentuating the role of erroneous epigenetic gene regulation in ALL.

Placental changes caused by food restriction during early pregnancy in mice are reversible

In a previous study, 50% calorie restriction in mice from d1.5 to 11.5 of pregnancy resulted in reduced placental weights and areas,relative sparing of labyrinth zone area compared to junctional zone area, and dramatic changes in global gene expression profiles.However, little lasting effect was seen on adult offspring of these pregnancies, with a slight reduction in adiposity in males and some changes in liver gene expression in both sexes. The goals of the present study were to determine whether the placental changes induced by caloric restriction in early pregnancy had permanent, irreversible effects on the placenta, and whether the changes in liver gene expression in adult offspring were present before birth. There were no differences in placental weights or areas, or the areas of individual placental zones near term in mice that had previously been food restricted. Global gene expression profiles at d18.5 were indistinguishable in placentas from control and previously food-restricted mothers. In fetuses from restricted dams at d18.5, liver expression of Gck, a key regulator of glycogen synthesis, was reduced, whereas its expression was increased in livers from adult offspring of restricted dams. Ppara expression was also reduced in fetal livers from restricted dams at d18.5, but not in adult offspring livers. We conclude that alterations in the placenta caused by nutrient restriction in early pregnancy are reversible, and that alterations in gene expression in livers of adult offspring are not a result of changes initiated during pregnancy and maintained through adulthood.

Serum miR-122 may serve as a biomarker for response to direct acting antivirals: effect of paritaprevir/R with dasabuvir or ombitasvir on miR-122 in HCV-infected subjects

Circulating microRNAs (miRNA) have been intensely investigated as biomarkers in disease and therapy. Several studies have identified miR-122 as an important regulator of HCV replication. The effect of new therapies that directly target the HCV replication life cycle on circulating microRNA levels has not been elucidated. We performed expression profiling of circulating miRNA in serum in subjects treated with HCV direct-acting antiviral agents (DAAs). Serum miRNA levels were evaluated from two studies in HCV GT1-infected treatment-naïve subjects and prior nonresponders to pegylated interferon (pegIFN) and ribavirin (RBV) who received paritaprevir/ritonavir + dasabuvir + RBV for 12 weeks, and in treatment-naïve genotype (GT)1-3-infected subjects who received paritaprevir/ritonavir + ombitasvir ± RBV for 12 weeks. Over 100 different miRNA species were detected in serum. Of these, levels of miR-122 showed the most consistent change in response to treatment across all HCV genotypes. In all subjects, miR-122 showed an average four-fold reduction between baseline and week 2, and remained below baseline through post-treatment week 12 in subjects who achieved sustained virological response. In contrast, in subjects who did not achieve SVR, miR-122 levels began to return to baseline levels after the second week of treatment. The change in miR-122 levels was similar across genotypes, and was comparable with or without RBV. This is the first report comparing expression levels of circulating miRNA in HCV GT1-3 subjects treated with IFN-free combinations of DAAs. The results suggest that serum levels of miR-122 are reduced following treatment in subjects who achieve SVR, and correlate with HCV RNA levels across genotypes.

Improved Murine Blastocyst Quality and Development in a Single Culture Medium Compared to Sequential Culture Media

OBJECTIVE

Validate single versus sequential culture media for murine embryo development.

DESIGN

Prospective laboratory experiment.

SETTING

Assisted Reproduction Laboratory.

ANIMALS

Murine embryos.

INTERVENTIONS

Thawed murine zygotes cultured for 3 or 5 days (d3 or d5) in single or sequential embryo culture media developed for human in vitro fertilization.

MAIN OUTCOME MEASURES

On d3, zygotes developing to the 8 cell (8C) stage or greater were quantified using 4',6-diamidino-2-phenylindole (DAPI), and quality was assessed by morphological analysis. On d5, the number of embryos reaching the blastocyst stage was counted. DAPI was used to quantify total nuclei and inner cell mass nuclei. Localization of ubiquitin C-terminal hydrolase L1 (UCHL1) and ubiquitin C-terminal hydrolase L3 (UCHL3) was reference points for evaluating cell quality.

RESULTS

Comparing outcomes in single versus to sequential media, the odds of embryos developing to the 8C stage on d3 were 2.34 time greater (P = .06). On d5, more embryos reached the blastocyst stage (P = <.0001), hatched, and had significantly more trophoblast cells (P = .005) contributing to the increased total cell number. Also at d5, localization of distinct cytoplasmic UCHL1 and nuclear UCHL3 was found in high-quality hatching blastocysts. Localization of UCHL1 and UCHL3 was diffuse and inappropriately dispersed throughout the cytoplasm in low-quality nonhatching blastocysts.

CONCLUSIONS

Single medium yields greater cell numbers, an increased growth rate, and more hatching of murine embryos. Cytoplasmic UCHL1 and nuclear UHCL3 localization patterns were indicative of embryo quality. Our conclusions are limited to murine embryos but one might speculate that single medium may also be more beneficial for human embryo culture. Human embryo studies are needed.

Soil bioretention protects juvenile salmon and their prey from the toxic impacts of urban stormwater runoff

Green stormwater infrastructure (GSI), or low impact development, encompasses a diverse and expanding portfolio of strategies to reduce the impacts of stormwater runoff on natural systems. Benchmarks for GSI success are usually framed in terms of hydrology and water chemistry, with reduced flow and loadings of toxic chemical contaminants as primary metrics. Despite the central goal of protecting aquatic species abundance and diversity, the effectiveness of GSI treatments in maintaining diverse assemblages of sensitive aquatic taxa has not been widely evaluated. In the present study we characterized the baseline toxicity of untreated urban runoff from a highway in Seattle, WA, across six storm events. For all storms, first flush runoff was toxic to the daphniid Ceriodaphnia dubia, causing up to 100% mortality or impairing reproduction among survivors. We then evaluated whether soil media used in bioretention, a conventional GSI method, could reduce or eliminate toxicity to juvenile coho salmon (Oncorhynchus kisutch) as well as their macroinvertebrate prey, including cultured C. dubia and wild-collected mayfly nymphs (Baetis spp.). Untreated highway runoff was generally lethal to salmon and invertebrates, and this acute mortality was eliminated when the runoff was filtered through soil media in bioretention columns. Soil treatment also protected against sublethal reproductive toxicity in C. dubia. Thus, a relatively inexpensive GSI technology can be highly effective at reversing the acutely lethal and sublethal effects of urban runoff on multiple aquatic species.

Genome-wide DNA methylation analysis in precursor B-cells

DNA methylation is responsible for regulating gene expression and cellular differentiation and for maintaining genomic stability during normal human development. Furthermore, it plays a significant role in the regulation of hematopoiesis. In order to elucidate the influence of DNA methylation during B-cell development, genome-wide DNA methylation status of pro-B, pre-BI, pre-BII, and naïve-B-cells isolated from human umbilical cord blood was determined using the methylated CpG island recovery assay followed by next generation sequencing. On average, 182–200 million sequences were generated for each precursor B-cell subset in 10 biological replicates. An overall decrease in methylation was observed during the transition from pro-B to pre-BI, whereas no differential methylation was observed in the pre-BI to pre-BII transition or in the pre-BII to naïve B-cell transition. Most of the methylated regions were located within intergenic and intronic regions not present in a CpG island context. Putative novel enhancers were identified in these regions that were differentially methylated between pro-B and pre-BI cells. The genome-wide methylation profiles are publically available and may be used to gain a better understanding of the involvement of atypical DNA methylation in the pathogenesis of malignancies associated with precursor B-cells.

Exercise training causes differential changes in gene expression in diaphragm arteries and 2A arterioles of obese rats

We employed next-generation, transcriptome-wide RNA sequencing (RNA-Seq) technology to assess the effects of two different exercise training protocols on transcriptional profiles in diaphragm second-order arterioles (D2a) and in the diaphragm feed artery (DFA) from Otsuka Long Evans Tokushima Fatty (OLETF) rats. Arterioles were isolated from the diaphragm of OLETF rats that underwent an endurance exercise training program (EX; n = 13), interval sprint training program (SPRINT; n = 14), or remained sedentary (Sed; n = 12). Our hypothesis was that exercise training would have similar effects on gene expression in the diaphragm and soleus muscle arterioles because diaphragm blood flow increases during exercise to a similar extent as in soleus. Results reveal that several canonical pathways that were significantly altered by exercise in limb skeletal muscles were not among the pathways significantly changed in the diaphragm arterioles including actin cytoskeleton signaling, role of NFAT in regulation of immune response, protein kinase A signaling, and protein ubiquitination pathway. EX training altered the expression of a smaller number of genes than did SPRINT in the DFA but induced a larger number of genes with altered expression in the D2a than did SPRINT. In fact, FDR differential expression analysis (FDR, 10%) indicated that only two genes exhibited altered expression in D2a of SPRINT rats. Very few of the genes that exhibited altered expression in the DFA or D2a were also altered in limb muscle arterioles. Finally, results indicate that the 2a arterioles of soleus muscle (S2a) from endurance-trained animals and the DFA of SPRINT animals exhibited the largest number of genes with altered expression.

Exercise-induced differential changes in gene expression among arterioles of skeletal muscles of obese rats

Using next-generation, transcriptome-wide RNA sequencing (RNA-Seq) technology we assessed the effects of exercise training on transcriptional profiles in skeletal muscle arterioles isolated from the soleus and gastrocnemius muscles of Otsuka Long Evans Tokushima Fatty (OLETF) rats that underwent an endurance exercise training program (EX; n = 13), interval sprint training program (SPRINT; n = 14), or remained sedentary (Sed; n = 12). We hypothesized that the greatest effects of exercise would be in the gastrocnemius arterioles. Results show that EX caused the largest number of changes in gene expression in the soleus and white gastrocnemius 2a arterioles with little to no changes in the feed arteries. In contrast, SPRINT caused substantial changes in gene expression in the feed arteries. IPA canonical pathway analysis revealed 18 pathways with significant changes in gene expression when analyzed across vessels and revealed that EX induces increased expression of the following genes in all arterioles examined: Shc1, desert hedgehog protein (Dhh), adenylate cyclase 4 (Adcy4), G protein binding protein, alpha (Gnat1), and Bcl2l1 and decreased expression of ubiquitin D (Ubd) and cAMP response element modulator (Crem). EX increased expression of endothelin converting enzyme (Ece1), Hsp90b, Fkbp5, and Cdcl4b in four of five arterioles. SPRINT had effects on expression of Crem, Dhh, Bcl2l1, and Ubd that were similar to EX. SPRINT also increased expression of Nfkbia, Hspa5, Tubb 2a and Tubb 2b, and Fkbp5 in all five arterioles and increased expression of Gnat1 in all but the soleus second-order arterioles. Many contractile and/or structural protein genes were increased by SPRINT in the gastrocnemius feed artery, but the same genes exhibited decreased expression in red gastrocnemius arterioles. We conclude that training-induced changes in arteriolar gene expression patterns differ by muscle fiber type composition and along the arteriolar tree.

Transcriptomic-Wide Discovery of Direct and Indirect HuR RNA Targets in Activated CD4+ T Cells

Due to poor correlation between steady state mRNA levels and protein product, purely transcriptomic profiling methods may miss genes posttranscriptionally regulated by RNA binding proteins (RBPs) and microRNAs (miRNAs). RNA immunoprecipitation (RIP) methods developed to identify in vivo targets of RBPs have greatly elucidated those mRNAs which may be regulated via transcript stability and translation. The RBP HuR (ELAVL1) and family members are major stabilizers of mRNA. Many labs have identified HuR mRNA targets; however, many of these analyses have been performed in cell lines and oftentimes are not independent biological replicates. Little is known about how HuR target mRNAs behave in conditional knock-out models. In the present work, we performed HuR RIP-Seq and RNA-Seq to investigate HuR direct and indirect targets using a novel conditional knock-out model of HuR genetic ablation during CD4+ T activation and Th2 differentiation. Using independent biological replicates, we generated a high coverage RIP-Seq data set (>160 million reads) that was analyzed using bioinformatics methods specifically designed to find direct mRNA targets in RIP-Seq data. Simultaneously, another set of independent biological replicates were sequenced by RNA-Seq (>425 million reads) to identify indirect HuR targets. These direct and indirect targets were combined to determine canonical pathways in CD4+ T cell activation and differentiation for which HuR plays an important role. We show that HuR may regulate genes in multiple canonical pathways involved in T cell activation especially the CD28 family signaling pathway. These data provide insights into potential HuR-regulated genes during T cell activation and immune mechanisms.

Integrating chemohormonal therapy and surgery in known or suspected lymph node metastatic prostate cancer

BACKGROUND

Prostate cancer persisting in the primary site after systemic therapy may contribute to emergence of resistance and progression. We previously demonstrated molecular characteristics of lethal cancer in the prostatectomy specimens of patients presenting with lymph node metastasis after chemohormonal treatment. Here we report the post-treatment outcomes of these patients and assess whether a link exists between surgery and treatment-free/cancer-free survival.

METHODS

Patients with either clinically detected lymph node metastasis or primaries at high risk for nodal dissemination were treated with androgen ablation and docetaxel. Those responding with PSA concentration <1 ng ml(-1) were recommended surgery 1 year from enrollment. ADT was withheld postoperatively. The rate of survival without biochemical progression 1 year after surgery was measured to screen for efficacy.

RESULTS

Forty patients were enrolled and 39 were evaluable. Three patients (7.7%) declined surgery. Of the remaining 36, 4 patients experienced disease progression during treatment and 4 more did not reach PSA <1. Twenty-six patients (67%) completed surgery, and 13 (33%) were also progression-free 1 year postoperatively (8 with undetectable PSA). With a median follow-up of 61 months, time to treatment failure was 27 months in the patients undergoing surgery. The most frequent patterns of first disease recurrence were biochemical (10 patients) and systemic (5).

CONCLUSIONS

Half of the patients undergoing surgery were off treatment and progression-free 1 year following completion of all therapy. These results suggest that integration of surgery is feasible and may be superior to systemic therapy alone for selected prostate cancer patients presenting with nodal metastasis.

Effects of vendor and genetic background on the composition of the fecal microbiota of inbred mice

The commensal gut microbiota has been implicated as a determinant in several human diseases and conditions. There is mounting evidence that the gut microbiota of laboratory mice (Mus musculus) similarly modulates the phenotype of mouse models used to study human disease and development. While differing model phenotypes have been reported using mice purchased from different vendors, the composition and uniformity of the fecal microbiota in mice of various genetic backgrounds from different vendors is unclear. Using culture-independent methods and robust statistical analysis, we demonstrate significant differences in the richness and diversity of fecal microbial populations in mice purchased from two large commercial vendors. Moreover, the abundance of many operational taxonomic units, often identified to the species level, as well as several higher taxa, differed in vendor- and strain-dependent manners. Such differences were evident in the fecal microbiota of weanling mice and persisted throughout the study, to twenty-four weeks of age. These data provide the first in-depth analysis of the developmental trajectory of the fecal microbiota in mice from different vendors, and a starting point from which researchers may be able to refine animal models affected by differences in the gut microbiota and thus possibly reduce the number of animals required to perform studies with sufficient statistical power.

Identification of a human papillomavirus-associated oncogenic miRNA panel in human oropharyngeal squamous cell carcinoma validated by bioinformatics analysis of the Cancer Genome Atlas

High-risk human papillomavirus (HPV) is a causative agent for an increasing subset of oropharyngeal squamous cell carcinomas (OPSCCs), and current evidence supports these tumors as having identifiable risk factors and improved response to therapy. However, the biochemical and molecular alterations underlying the pathobiology of HPV-associated OPSCC (designated HPV(+) OPSCC) remain unclear. Herein, we profile miRNA expression patterns in HPV(+) OPSCC to provide a more detailed understanding of pathologic molecular events and to identify biomarkers that may have applicability for early diagnosis, improved staging, and prognostic stratification. Differentially expressed miRNAs were identified in RNA isolated from an initial clinical cohort of HPV(+/-) OPSCC tumors by quantitative PCR-based miRNA profiling. This oncogenic miRNA panel was validated using miRNA sequencing and clinical data from The Cancer Genome Atlas and miRNA in situ hybridization. The HPV-associated oncogenic miRNA panel has potential utility in diagnosis and disease stratification and in mechanistic elucidation of molecular factors that contribute to OPSCC development, progression, and differential response to therapy.

Zebrafish and clean water technology: assessing soil bioretention as a protective treatment for toxic urban runoff

Urban stormwater contains a complex mixture of contaminants that can be acutely toxic to aquatic biota. Green stormwater infrastructure (GSI) is a set of evolving technologies intended to reduce impacts on natural systems by slowing and filtering runoff. The extent to which GSI methods work as intended is usually assessed in terms of water quantity (hydrology) and quality (chemistry). Biological indicators of GSI effectiveness have received less attention, despite an overarching goal of protecting the health of aquatic species. Here we use the zebrafish (Danio rerio) experimental model to evaluate bioinfiltration as a relatively inexpensive technology for treating runoff from an urban highway with dense motor vehicle traffic. Zebrafish embryos exposed to untreated runoff (48-96h; six storm events) displayed an array of developmental abnormalities, including delayed hatching, reduced growth, pericardial edema, microphthalmia (small eyes), and reduced swim bladder inflation. Three of the six storms were acutely lethal, and sublethal toxicity was evident across all storms, even when stormwater was diluted by as much as 95% in clean water. As anticipated from exposure to cardiotoxic polycyclic aromatic hydrocarbons (PAHs), untreated runoff also caused heart failure, as indicated by circulatory stasis, pericardial edema, and looping defects. Bioretention treatment dramatically improved stormwater quality and reversed nearly all forms of developmental toxicity. The zebrafish model therefore provides a versatile experimental platform for rapidly assessing GSI effectiveness.

Novel microRNA correlations in the severely injured

PURPOSE

Severe injury initiates an inflammatory response that can perpetuate immunological dysfunction, uncontrolled inflammation, and subsequent multisystem organ failure. MicroRNAs (miRNAs) have recently been identified as regulators of this inflammatory response. Our study sought to identify the differential expression of unique miRNAs and their correlations with genes of the Toll-like receptor (TLR) pathways, and clinical parameters in the severely injured.

METHODS

Fourteen trauma patients requiring transfusion were prospectively enrolled in this institutional review board-approved study. Inclusion criteria consisted of adult patients deemed clinically to be in hemorrhagic shock necessitating transfusion in the acute phase of their injury care. Peripheral blood samples were obtained after admission to the surgical intensive care unit. Expression of circulating mature miRNA from each patient, as well as from 10 healthy, age-matched controls, was determined and compared using the HiSeq 2500 sequencing system and the R software system. Gene expression of TLR signaling pathways for each patient was examined using custom gene expression polymerase chain reaction arrays. Statistical analyses were performed using general linear models and empirical Bayes methods to determine differential expression and Spearman's nonparametric correlation analysis.

RESULTS

Subjects were 21-77 years old (mean, 42), 80% male, Injury Severity Score 11-43 (mean, 26), with 11 blunt and 3 penetrating injuries. Three were intubated and 5 received blood products before arrival. Base deficit upon hospital admission was 3 to 20 (mean, 9). All patients required blood transfusion secondary to blood loss sustained during injury. Survival to discharge was 93%. Controls were 27-64 years old (mean, 40) and 60% male. Sequencing analysis revealed 69 differentially expressed miRNAs (P < .05) in the severely injured. Within the differentially expressed miRNAs, there were 12 direct and 6 indirect correlations with multiple genes involved in the TLR3 and TLR4 signaling pathways. The relationships between these same miRNAs and clinical parameters were also analyzed. We discovered 4 direct correlations with base deficit and HCO3, and 7 indirect correlations involving total fresh frozen plasma transfused, base deficit, HCO3, and PaCO2 levels.

CONCLUSION

Differential expression and correlations between miRNAs, genes of the TLR pathways, and clinical parameters are unique findings in the severely injured and may lead to a greater understanding of the regulation of sterile inflammation after severe injury.

Identification of genes whose expression is altered by obesity throughout the arterial tree

We used next-generation RNA sequencing (RNA-Seq) technology on the whole transcriptome to identify genes whose expression is consistently affected by obesity across multiple arteries. Specifically, we examined transcriptional profiles of the iliac artery as well as the feed artery, first, second, and third branch order arterioles in the soleus, gastrocnemius, and diaphragm muscles from obese Otsuka Long-Evans Tokushima Fatty (OLETF) and lean Long-Evans Tokushima Otsuka (LETO) rats. Within the gastrocnemius and soleus muscles, the number of genes differentially expressed with obesity tended to increase with increasing branch order arteriole number (i.e., decreasing size of the artery). This trend was opposite in the diaphragm. We found a total of 15 genes that were consistently upregulated with obesity (MIS18A, CTRB1, FAM151B, FOLR2, PXMP4, OAS1B, SREBF2, KLRA17, SLC25A44, SNX10, SLFN3, MEF2BNB, IRF7, RAD23A, LGALS3BP) and five genes that were consistently downregulated with obesity (C2, GOLGA7, RIN3, PCP4, CYP2E1). A small fraction (∼9%) of the genes affected by obesity was modulated across all arteries examined. In conclusion, the present study identifies a select number of genes (i.e., 20 genes) whose expression is consistently altered throughout the arterial network in response to obesity and provides further insight into the heterogeneous vascular effects of obesity. Although there is no known direct function of the majority of 20 genes related to vascular health, the obesity-associated upregulation of SREBF2, LGALS3BP, IRF7, and FOLR2 across all arteries is suggestive of an unfavorable vascular phenotypic alteration with obesity. These data may serve as an important resource for identifying novel therapeutic targets against obesity-related vascular complications.

The role of estrogen signaling in a mouse model of inflammatory bowel disease: a Helicobacter hepaticus model

The pathogenesis of inflammatory bowel diseases (IBD), Crohn's disease and ulcerative colitis, is due in part to interactions between the immune system, genetics, the environment, and endogenous microbiota. Gonadal sex hormones (GSH), such as estrogen, are thought to be involved in the development of IBD as variations in disease severity occur during pregnancy, menopause, or oral contraceptives use. In certain strains of mice, infection with Helicobacter hepaticus triggers IBD-like mucosal inflammation that is more severe in female mice than in males, suggesting a role for GSH in this model. To determine the role of estrogen signaling in microbiota-induced intestinal inflammation, estrogen receptor (ER) α and β knock-out (KO) mice, ER agonists, and adoptive transfers were utilized. We demonstrate that, when signaling is limited to ERβ on a non-CD4⁺ cell subset, disease is less severe and this correlates with decreased expression of pro-inflammatory mediators.

Estrogen and androgen receptor activities of hydraulic fracturing chemicals and surface and ground water in a drilling-dense region

The rapid rise in natural gas extraction using hydraulic fracturing increases the potential for contamination of surface and ground water from chemicals used throughout the process. Hundreds of products containing more than 750 chemicals and components are potentially used throughout the extraction process, including more than 100 known or suspected endocrine-disrupting chemicals. We hypothesized that a selected subset of chemicals used in natural gas drilling operations and also surface and ground water samples collected in a drilling-dense region of Garfield County, Colorado, would exhibit estrogen and androgen receptor activities. Water samples were collected, solid-phase extracted, and measured for estrogen and androgen receptor activities using reporter gene assays in human cell lines. Of the 39 unique water samples, 89%, 41%, 12%, and 46% exhibited estrogenic, antiestrogenic, androgenic, and antiandrogenic activities, respectively. Testing of a subset of natural gas drilling chemicals revealed novel antiestrogenic, novel antiandrogenic, and limited estrogenic activities. The Colorado River, the drainage basin for this region, exhibited moderate levels of estrogenic, antiestrogenic, and antiandrogenic activities, suggesting that higher localized activity at sites with known natural gas-related spills surrounding the river might be contributing to the multiple receptor activities observed in this water source. The majority of water samples collected from sites in a drilling-dense region of Colorado exhibited more estrogenic, antiestrogenic, or antiandrogenic activities than reference sites with limited nearby drilling operations. Our data suggest that natural gas drilling operations may result in elevated endocrine-disrupting chemical activity in surface and ground water.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. I. Impact of obesity

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the influence of obesity on global gene expression in skeletal muscle feed arteries. Transcriptional profiles of the gastrocnemius and soleus muscle feed arteries (GFA and SFA, respectively) and aortic endothelial cell-enriched samples from obese Otsuka Long-Evans Tokushima Fatty (OLETF) and lean Long-Evans Tokushima Otsuka (LETO) rats were examined. Obesity produced 282 upregulated and 133 downregulated genes in SFA and 163 upregulated and 77 downregulated genes in GFA [false discovery rate (FDR) < 10%] with an overlap of 93 genes between the arteries. In LETO rats, there were 89 upregulated and 114 downregulated genes in the GFA compared with the SFA. There were 244 upregulated and 275 downregulated genes in OLETF rats (FDR < 10%) in the GFA compared with the SFA, with an overlap of 76 differentially expressed genes common to both LETO and OLETF rats in both the GFA and SFA. A total of 396 transcripts were found to be differentially expressed between LETO and OLETF in aortic endothelial cell-enriched samples. Overall, we found 1) the existence of heterogeneity in the transcriptional profile of the SFA and GFA within healthy LETO rats, 2) that this between-vessel heterogeneity was markedly exacerbated in the hyperphagic, obese OLETF rat, and 3) a greater number of genes whose expression was altered by obesity in the SFA compared with the GFA. Also, results indicate that in OLETF rats the GFA takes on a relatively more proatherogenic phenotype compared with the SFA.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. II. Impact of exercise training in obesity

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the extent to which exercise training alters global gene expression in skeletal muscle feed arteries and aortic endothelial cells of obese Otsuka Long-Evans Tokushima Fatty (OLETF) rats. Transcriptional profiles of the soleus and gastrocnemius muscle feed arteries (SFA and GFA, respectively) and aortic endothelial cell-enriched samples from rats that underwent an endurance exercise training program (EndEx; n = 12) or a interval sprint training program (IST; n = 12) or remained sedentary (Sed; n = 12) were examined. In response to EndEx, there were 39 upregulated (e.g., MANF) and 20 downregulated (e.g., ALOX15) genes in SFA and 1 upregulated (i.e., Wisp2) and 1 downregulated (i.e., Crem) gene in GFA [false discovery rate (FDR) < 10%]. In response to IST, there were 305 upregulated (e.g., MANF, HSPA12B) and 324 downregulated genes in SFA and 101 upregulated and 66 downregulated genes in GFA, with an overlap of 32 genes between arteries. Furthermore, in aortic endothelial cells, there were 183 upregulated (e.g., eNOS, SOD-3) and 141 downregulated (e.g., ATF3, Clec1b, npy, leptin) genes with EndEx and 71 upregulated and 69 downregulated genes with IST, with an overlap of 35 between exercise programs. Expression of only two genes (Tubb2b and Slc9a3r2) was altered (i.e., increased) by exercise in all three arteries. The finding that both EndEx and IST produced greater transcriptional changes in the SFA compared with the GFA is intriguing when considering the fact that treadmill bouts of exercise are associated with greater relative increases in blood flow to the gastrocnemius muscle compared with the soleus muscle.

Proteomic mapping of proteins released during necrosis and apoptosis from cultured neonatal cardiac myocytes

Cardiac injury induces myocyte apoptosis and necrosis, resulting in the secretion and/or release of intracellular proteins. Currently, myocardial injury can be detected by analysis of a limited number of biomarkers in blood or coronary artery perfusate. However, the complete proteomic signature of protein release from necrotic cardiac myocytes is unknown. Therefore, we undertook a proteomic-based study of proteins released from cultured neonatal rat cardiac myocytes in response to H2O2 (necrosis) or staurosporine (apoptosis) to identify novel specific markers of cardiac myocyte cell death. Necrosis and apoptosis resulted in the identification of 147 and 79 proteins, respectively. Necrosis resulted in a relative increase in the amount of many proteins including the classical necrotic markers lactate dehydrogenase (LDH), high-mobility group B1 (HMGB1), myoglobin, enolase, and 14-3-3 proteins. Additionally, we identified several novel markers of necrosis including HSP90, α-actinin, and Trim72, many of which were elevated over control levels earlier than classical markers of necrotic injury. In contrast, the majority of identified proteins remained at low levels during apoptotic cell death, resulting in no candidate markers for apoptosis being identified. Blotting for a selection of these proteins confirmed their release during necrosis but not apoptosis. We were able to confirm the presence of classical necrotic markers in the extracellular milieu of necrotic myocytes. We also were able to identify novel markers of necrotic cell death with relatively early release profiles compared with classical protein markers of necrosis. These results have implications for the discovery of novel biomarkers of necrotic myocyte injury, especially in the context of ischemia-reperfusion injury.

Effect of home telemonitoring on glycemic and blood pressure control in primary care clinic patients with diabetes

OBJECTIVE

Patient self-management support may be augmented by using home-based technologies that generate data points that providers can potentially use to make more timely changes in the patients' care. The purpose of this study was to evaluate the effectiveness of short-term targeted use of remote data transmission on treatment outcomes in patients with diabetes who had either out-of-range hemoglobin A1c (A1c) and/or blood pressure (BP) measurements.

MATERIALS AND METHODS

A single-center randomized controlled clinical trial design compared in-home monitoring (n=55) and usual care (n=53) in patients with type 2 diabetes and hypertension being treated in primary care clinics. Primary outcomes were A1c and systolic BP after a 12-week intervention.

RESULTS

There were no significant differences between the intervention and control groups on either A1c or systolic BP following the intervention.

CONCLUSIONS

The addition of technology alone is unlikely to lead to improvements in outcomes. Practices need to be selective in their use of telemonitoring with patients, limiting it to patients who have motivation or a significant change in care, such as starting insulin. Attention to the need for effective and responsive clinic processes to optimize the use of the additional data is also important when implementing these types of technology.

Vascular transcriptional alterations produced by juvenile obesity in Ossabaw swine

We adopted a transcriptome-wide microarray analysis approach to determine the extent to which vascular gene expression is altered as a result of juvenile obesity and identify obesity-responsive mRNAs. We examined transcriptional profiles in the left anterior descending coronary artery (LAD), perivascular fat adjacent to the LAD, and descending thoracic aorta between obese (n = 5) and lean (n = 6) juvenile Ossabaw pigs (age = 22 wk). Obesity was experimentally induced by feeding the animals a high-fat/high-fructose corn syrup/high-cholesterol diet for 16 wk. We found that expression of 189 vascular cell genes in the LAD and expression of 165 genes in the thoracic aorta were altered with juvenile obesity (false discovery rate ≤ 10%) with an overlap of only 28 genes between both arteries. Notably, a number of genes found to be markedly upregulated in the LAD of obese pigs are implicated in atherosclerosis, including ACP5, LYZ, CXCL14, APOE, PLA2G7, LGALS3, SPP1, ITGB2, CYBB, and P2RY12. Furthermore, pathway analysis revealed the induction of proinflammatory and pro-oxidant pathways with obesity primarily in the LAD. Gene expression in the LAD perivascular fat was minimally altered with juvenile obesity. Together, we provide new evidence that obesity produces artery-specific changes in pretranslational regulation with a clear upregulation of proatherogenic genes in the LAD. Our data may offer potential viable drug targets and mechanistic insights regarding the molecular precursors involved in the origins of overnutrition and obesity-associated vascular disease. In particular, our results suggest that the oxidized LDL/LOX-1/NF-κB signaling axis may be involved in the early initiation of a juvenile obesity-induced proatherogenic coronary artery phenotype.

Large-scale analysis of DNA methylation in chronic lymphocytic leukemia

AIMS

B-cell chronic lymphocytic leukemia (CLL) is a heterogeneous malignancy that clinically ranges from indolent to rapidly progressive. CLL, like other cancers, can be affected by epigenetic alterations.

MATERIALS & METHODS

A microarray discovery-based study was initiated to determine DNA methylation in CLL cases with a range of CD38 expression (1–92%).

RESULTS

Many loci were either methylated or unmethylated across all CD38 levels, but differential methylation was also observed for some genes. Genomic sequencing of DLEU7 confirmed extensive cytosine methylation preferentially in patient samples with low CD38 expression, whereas NRP2, SFRP2 and ADAM12 were more commonly methylated in those with high CD38 expression.

CONCLUSION

This study demonstrates that CLL is affected by CpG island methylation in some genes that segregate with CD38 expression levels, while most others show similar methylation patterns across all levels. The CpG island methylation in certain functional gene groups and pathway-associated genes that are known to be deregulated in CLL provides additional insights into the CLL methylome and epigenetic contribution to cellular dysfunction. It will now be useful to investigate the effectiveness of epigenetic therapeutic reversal of these alterations to develop effective treatments for the disease.

Early depression of Ankrd2 and Csrp3 mRNAs in the polyribosomal and whole tissue fractions in skeletal muscle with decreased voluntary running

The wheel-lock (WL) model for depressed ambulatory activity in rats has shown metabolic maladies ensuing within 53-173 h after WL begins. We sought to determine if WL beginning after 21-23 days of voluntary running in growing female Wistar rats affected the mRNA profile in the polyribosomal fraction from plantaris muscle shortly following WL. In experiment 1, WL occurred at 0200 and muscles were harvested at 0700 daily at 5 h (WL5h, n = 4), 29 h (WL29h, n = 4), or 53 h (WL53h, n = 4) after WL. Affymetrix Rat Gene 1.0 ST Arrays were used to test the initial question as to whether WL affects mRNA occupancy on skeletal muscle polyribosomes. Using a false discovery rate of 15%, no changes in mRNAs in the polyribosomal fraction were observed at WL29h and eight mRNAs (of over 8,200 identified targets) were altered at WL53h compared with WL5h. Interestingly, two of the six downregulated genes included ankyrin repeat domain 2 (Ankrd2) and cysteine-rich protein 3/muscle LIM protein (Csrp3), both of which encode mechanical stretch sensors and RT-PCR verified their WL-induced decline. In experiment 2, whole muscle mRNA and protein levels were analyzed for Ankrd2 and Csrp3 from the muscles of WL5h (4 original samples + 2 new), WL29h (4 original), WL53h (4 original + 2 new), as well as WL173 h (n = 6 new) and animals that never ran (SED, 4-5 new). Relative to WL5h controls, whole tissue Ankrd2 and Csrp3 mRNAs were lower (P < 0.05) at WL53h, WL173h, and SED; Ankrd2 protein tended to decrease at WL53h (P = 0.054) and Csrp3 protein was less in WL173h and SED rats (P < 0.05). In summary, unique early declines in Ankrd2 and Csrp3 mRNAs were identified with removal of voluntary running, which was subsequently followed by declines in Csrp3 protein levels during longer periods of wheel lock.

Impact of exercise training on endothelial transcriptional profiles in healthy swine: a genome-wide microarray analysis

While the salutary effects of exercise training on conduit artery endothelial cells have been reported in animals and humans with cardiovascular risk factors or disease, whether a healthy endothelium is alterable with exercise training is less certain. The purpose of this study was to evaluate the impact of exercise training on transcriptional profiles in normal endothelial cells using a genome-wide microarray analysis. Brachial and internal mammary endothelial gene expression was compared between a group of healthy pigs that exercise trained for 16-20 wk (n = 8) and a group that remained sedentary (n = 8). We found that a total of 130 genes were upregulated and 84 genes downregulated in brachial artery endothelial cells with exercise training (>1.5-fold and false discovery rate <15%). In contrast, a total of 113 genes were upregulated and 31 genes downregulated in internal mammary artery endothelial cells using the same criteria. Although there was an overlap of 66 genes (59 upregulated and 7 downregulated with exercise training) between the brachial and internal mammary arteries, the identified endothelial gene networks and biological processes influenced by exercise training were distinctly different between the brachial and internal mammary arteries. These data indicate that a healthy endothelium is indeed responsive to exercise training and support the concept that the influence of physical activity on endothelial gene expression is not homogenously distributed throughout the vasculature.

Overexpression of the RNA binding protein HuR impairs tumor growth in triple negative breast cancer associated with deficient angiogenesis

Interactions between RNA binding proteins (RBPs) and genes are not well understood, especially in regulation of angiogenesis. The RBP HuR binds to the AU-rich (ARE) regions of labile mRNAs, facilitating their translation into protein and has been hypothesized to be a tumor-maintenance gene. Elevated levels of cytoplasmic HuR directly correlate with increased invasiveness and poor prognosis for many cancers, including those of the breast. HuR controls the expression of multiple genes involved in angiogenesis including VEGFα, HIF1α and thrombospondin 1 (TSP1). We investigated the role of HuR in estrogen receptor negative (ER(-)) breast cancer. MDA-MB-231 cells with higher levels of HuR have alterations in cell cycle kinetics and faster growth. Unexpectedly, HuR overexpression significantly interfered with tumor growth in orthotopic mouse models. The putative mechanism seems to be an anti-angiogenetic effect by increasing expression of TSP1 but also surprisingly, downregulating VEGF, a target which HuR normally increases. Our findings reveal that HuR may be regulating a cluster of genes involved in blood vessel formation which controls tumor angiogenesis. An approach of modulating HuR levels may overcome limitations associated with monotherapies targeting tumor vessel formation.

Inflammatory cytokines differentially up-regulate human endometrial haptoglobin production in women with endometriosis

BACKGROUND

Evidence suggests that eutopic endometrium from women with endometriosis (US-E) has intrinsic functional anomalies compared with women without endometriosis (US-C). We hypothesized that differences in endometrial haptoglobin (eHp) mRNA and protein levels exist between eutopic endometrium from US-E and US-C and that inflammatory mediators may be involved.

METHODS

Endometrial stromal cells and tissue explants from US-E (n = 18) and US-C (n = 18) were cultured (24 h/48 h for cells/explants) with interleukin (IL)-1alpha, -1beta, -6, -8 or tumor necrosis factor-alpha (TNF-alpha) at 0-100 ng/ml. eHp protein in media and mRNA levels were quantified by enzyme-linked immunosorbent assay and quantitative PCR.

RESULTS

In eutopic endometrial stromal cells from US-E, IL-1beta, IL-6 and TNF-alpha (10 ng/ml) increased eHp mRNA levels (P = 0.002, P < 0.001 and P < 0.001, respectively) and eHp protein (P = 0.023, 0.031 and 0.006, respectively) versus control. In endometrial tissues from US-E, IL-1beta, IL-6 and TNF-alpha increased eHp mRNA (P < 0.001, P = 0.017 and P < 0.001, respectively) and eHp protein (P < 0.001, P = 0.007 and 0.039, respectively) versus control. IL-1alpha and IL-8 had small or no effects on isolated endometrial cells or tissues. In US-C, IL-1beta, IL-8 and TNF-alpha each reduced eHp mRNA in endometrial stromal cells (all P < 0.001) versus control; IL-1alpha and IL-6 had no effect. eHp mRNA increased in endometrial tissues from US-C in response to IL-1beta (P = 0.008), IL-6 (P = 0.015) and TNF-alpha (P = 0.031) versus control; IL-1alpha or IL-8 had no effect.

CONCLUSIONS

Endometrium from US-E differentially responds to specific inflammatory cytokines by production of eHp. We propose that up-regulation of endometrial eHp by inflammatory mediators disrupts normal endometrial function and may facilitate the pathogenesis of endometriosis.

The RNA binding protein HuR differentially regulates unique subsets of mRNAs in estrogen receptor negative and estrogen receptor positive breast cancer

BACKGROUND

The discordance between steady-state levels of mRNAs and protein has been attributed to posttranscriptional control mechanisms affecting mRNA stability and translation. Traditional methods of genome wide microarray analysis, profiling steady-state levels of mRNA, may miss important mRNA targets owing to significant posttranscriptional gene regulation by RNA binding proteins (RBPs).

METHODS

The ribonomic approach, utilizing RNA immunoprecipitation hybridized to microarray (RIP-Chip), provides global identification of putative endogenous mRNA targets of different RBPs. HuR is an RBP that binds to the AU-rich elements (ARE) of labile mRNAs, such as proto-oncogenes, facilitating their translation into protein. HuR has been shown to play a role in cancer progression and elevated levels of cytoplasmic HuR directly correlate with increased invasiveness and poor prognosis for many cancers, including those of the breast. HuR has been described to control genes in several of the acquired capabilities of cancer and has been hypothesized to be a tumor-maintenance gene, allowing for cancers to proliferate once they are established.

RESULTS

We used HuR RIP-Chip as a comprehensive and systematic method to survey breast cancer target genes in both MCF-7 (estrogen receptor positive, ER+) and MDA-MB-231 (estrogen receptor negative, ER-) breast cancer cell lines. We identified unique subsets of HuR-associated mRNAs found individually or in both cell types. Two novel HuR targets, CD9 and CALM2 mRNAs, were identified and validated by quantitative RT-PCR and biotin pull-down analysis.

CONCLUSION

This is the first report of a side-by-side genome-wide comparison of HuR-associated targets in wild type ER+ and ER- breast cancer. We found distinct, differentially expressed subsets of cancer related genes in ER+ and ER- breast cancer cell lines, and noted that the differential regulation of two cancer-related genes by HuR was contingent upon the cellular environment.

Urinary-type plasminogen activator receptor/alpha 3 beta 1 integrin signaling, altered gene expression, and oral tumor progression

Oral squamous cell carcinoma (OSCC) has 50% 5-year survival rate, highlighting our limited understanding of the molecular events that contribute to disease progression. Microarray analyses of primary oral tumors have identified urinary-type plasminogen activator (uPA) and its receptor (uPAR) as key genes associated with human OSCC progression. The uPAR functions as both a proteinase receptor and an integrin ligand, modifying proteolysis, migration, integrin signaling, and cellular transcription. In the current study, uPAR expression levels were modified in OSCC cells followed by analysis of tumor growth in an in vivo orthotopic xenograft model and by transcriptional profiling. Overexpression of uPAR resulted in more infiltrative and less differentiated tumors, with ill-defined borders, cytologic atypia, and enhanced vascularity. Analysis of serial sections of both murine experimental tumors and microarrayed human OSCC showed a statistically significant association between uPAR and alpha(3) integrin colocalization in areas exhibiting extracellular signal-regulated kinase phosphorylation, suggesting that uPAR/alpha(3) integrin interaction potentiates extracellular signal-regulated kinase signaling in vivo. This is supported by cDNA microarray analysis, which showed differential expression of 148 genes (113 upregulated and 35 downregulated). Validation of gene expression changes in human OSCC using immunohistochemistry and quantitative real-time PCR showed increased growth factors, proteinases/inhibitors, and matrix components in uPAR-overexpressing tumors. Together, these results support a model wherein increased uPAR expression promotes alpha(3)beta(1) integrin association, resulting in increased mitogen-activated protein kinase signaling and transcriptional activation, leading to the formation of more aggressive tongue tumors. This combined approach has efficacy to identify additional biomarkers and/or prognostic indicators associated with aggressive human OSCC.

Porous plug gas injection systems for studies of hydrocarbon dissociation and transport in the DIII-D tokamak

A probe has been designed, constructed, and successfully used to inject methane into the DIII-D lower divertor in a manner imitating natural release by chemical erosion. This porous plug injector (PPI) probe consists of a self-contained gas reservoir with an integrated pressure gauge and a 3 cm diameter porous surface through which gas is injected into the lower divertor of the tokamak. The probe is positioned flush with the divertor target surface by means of the divertor materials evaluation system. Two gas delivery systems were developed: in the first, gas flow is regulated by a remotely controlled microvalve and in the second by a fixed micro-orifice flow restrictor. Because of the large area of the porous surface through which gas is admitted, the injected hydrocarbon molecules see a local carbon surface (>90% carbon) similar to that seen by hydrocarbons being emitted by chemical sputtering from surrounding carbon tiles. The distributed gas source also reduces the disturbance to the local plasma while providing sufficient signal for spectroscopic detection. In situ spectroscopic measurements with the PPI in DIII-D allow the direct calibration of response for measured plasma conditions from a known influx of gas.