Matthew S. Forrest

Population genomics of human gene expression

Genetic variation influences gene expression, and this variation in gene expression can be efficiently mapped to specific genomic regions and variants. Here we have used gene expression profiling of Epstein-Barr virus–transformed lymphoblastoid cell lines of all 270 individuals genotyped in the HapMap Consortium to elucidate the detailed features of genetic variation underlying gene expression variation. We find that gene expression is heritable and that differentiation between populations is in agreement with earlier small-scale studies. A detailed association analysis of over 2.2 million common SNPs per population (5% frequency in HapMap) with gene expression identified at least 1,348 genes with association signals in cis and at least 180 in trans. Replication in at least one independent population was achieved for 37% of cis signals and 15% of trans signals, respectively. Our results strongly support an abundance of cis-regulatory variation in the human genome. Detection of trans effects is limited but suggests that regulatory variation may be the key primary effect contributing to phenotypic variation in humans. We also explore several methodologies that improve the current state of analysis of gene expression variation.

Discovery of Novel Biomarkers by Microarray Analysis of Peripheral Blood Mononuclear Cell Gene Expression in Benzene-Exposed Workers

Benzene is an industrial chemical and component of gasoline that is an established cause of leukemia. To better understand the risk benzene poses, we examined the effect of benzene exposure on peripheral blood mononuclear cell (PBMC) gene expression in a population of shoe-factory workers with well-characterized occupational exposures using microarrays and real-time polymerase chain reaction (PCR). PBMC RNA was stabilized in the field and analyzed using a comprehensive human array, the U133A/B Affymetrix GeneChip set. A matched analysis of six exposed–control pairs was performed. A combination of robust multiarray analysis and ordering of genes using paired t-statistics, along with bootstrapping to control for a 5% familywise error rate, was used to identify differentially expressed genes in a global analysis. This resulted in a set of 29 known genes being identified that were highly likely to be differentially expressed. We also repeated these analyses on a smaller subset of 508 cytokine probe sets and found that the expression of 19 known cytokine genes was significantly different between the exposed and the control subjects. Six genes were selected for confirmation by real-time PCR, and of these, CXCL16, ZNF331, JUN, and PF4 were the most significantly affected by benzene exposure, a finding that was confirmed in a larger data set from 28 subjects. The altered expression was not caused by changes in the makeup of the PBMC fraction. Thus, microarray analysis along with real-time PCR confirmation reveals that altered expressions of CXCL16, ZNF331, JUN, and PF4 are potential biomarkers of benzene exposure.

Polymorphisms in the estrogen receptor 1 and vitamin C and matrix metalloproteinase gene families are associated with susceptibility to lymphoma.

Background Non-Hodgkin lymphoma (NHL) is the fifth most common cancer in the U.S. and few causes have been identified. Genetic association studies may help identify environmental risk factors and enhance our understanding of disease mechanisms. Methodology/Principal Findings 768 coding and haplotype tagging SNPs in 146 genes were examined using Illumina GoldenGate technology in a large population-based case-control study of NHL in the San Francisco Bay Area (1,292 cases 1,375 controls are included here). Statistical analyses were restricted to HIV- participants of white non-Hispanic origin. Genes involved in steroidogenesis, immune function, cell signaling, sunlight exposure, xenobiotic metabolism/oxidative stress, energy balance, and uptake and metabolism of cholesterol, folate and vitamin C were investigated. Sixteen SNPs in eight pathways and nine haplotypes were associated with NHL after correction for multiple testing at the adjusted q<0.10 level. Eight SNPs were tested in an independent case-control study of lymphoma in Germany (494 NHL cases and 494 matched controls). Novel associations with common variants in estrogen receptor 1 (ESR1) and in the vitamin C receptor and matrix metalloproteinase gene families were observed. Four ESR1 SNPs were associated with follicular lymphoma (FL) in the U.S. study, with rs3020314 remaining associated with reduced risk of FL after multiple testing adjustments [odds ratio (OR) = 0.42, 95% confidence interval (CI) = 0.23–0.77) and replication in the German study (OR = 0.24, 95% CI = 0.06–0.94). Several SNPs and haplotypes in the matrix metalloproteinase-3 (MMP3) and MMP9 genes and in the vitamin C receptor genes, solute carrier family 23 member 1 (SLC23A1) and SLC23A2, showed associations with NHL risk. Conclusions/Significance Our findings suggest a role for estrogen, vitamin C and matrix metalloproteinases in the pathogenesis of NHL that will require further validation.

A human-driven decline in global burned area

Burn less, baby, burn less Humans have, and always have had, a major impact on wildfire activity, which is expected to increase in our warming world. Andela et al. use satellite data to show that, unexpectedly, global burned area declined by ∼25% over the past 18 years, despite the influence of climate. The decrease has been largest in savannas and grasslands because of agricultural expansion and intensification. The decline of burned area has consequences for predictions of future changes to the atmosphere, vegetation, and the terrestrial carbon sink. Science, this issue p. 1356 Global burned area has declined by ~25% over the past 18 years. Fire is an essential Earth system process that alters ecosystem and atmospheric composition. Here we assessed long-term fire trends using multiple satellite data sets. We found that global burned area declined by 24.3 ± 8.8% over the past 18 years. The estimated decrease in burned area remained robust after adjusting for precipitation variability and was largest in savannas. Agricultural expansion and intensification were primary drivers of declining fire activity. Fewer and smaller fires reduced aerosol concentrations, modified vegetation structure, and increased the magnitude of the terrestrial carbon sink. Fire models were unable to reproduce the pattern and magnitude of observed declines, suggesting that they may overestimate fire emissions in future projections. Using economic and demographic variables, we developed a conceptual model for predicting fire in human-dominated landscapes.

Rapid Detection of Mycobacterium tuberculosis by Recombinase Polymerase Amplification

Improved access to effective tests for diagnosing tuberculosis (TB) has been designated a public health priority by the World Health Organisation. In high burden TB countries nucleic acid based TB tests have been restricted to centralised laboratories and specialised research settings. Requirements such as a constant electrical supply, air conditioning and skilled, computer literate operators prevent implementation of such tests in many settings. Isothermal DNA amplification technologies permit the use of simpler, less energy intensive detection platforms more suited to low resource settings that allow the accurate diagnosis of a disease within a short timeframe. Recombinase Polymerase Amplification (RPA) is a rapid, low temperature isothermal DNA amplification reaction. We report here RPA-based detection of Mycobacterium tuberculosis complex (MTC) DNA in <20 minutes at 39°C. Assays for two MTC specific targets were investigated, IS6110 and IS1081. When testing purified MTC genomic DNA, limits of detection of 6.25 fg (IS6110) and 20 fg (IS1081)were consistently achieved. When testing a convenience sample of pulmonary specimens from suspected TB patients, RPA demonstrated superior accuracy to indirect fluorescence microscopy. Compared to culture, sensitivities for the IS1081 RPA and microscopy were 91.4% (95%CI: 85, 97.9) and 86.1% (95%CI: 78.1, 94.1) respectively (n = 71). Specificities were 100% and 88.6% (95% CI: 80.8, 96.1) respectively. For the IS6110 RPA and microscopy sensitivities of 87.5% (95%CI: 81.7, 93.2) and 70.8% (95%CI: 62.9, 78.7) were obtained (n = 90). Specificities were 95.4 (95% CI: 92.3,98.1) and 88% (95% CI: 83.6, 92.4) respectively. The superior specificity of RPA for detecting tuberculosis was due to the reduced ability of fluorescence microscopy to distinguish Mtb complex from other acid fast bacteria. The rapid nature of the RPA assay and its low energy requirement compared to other amplification technologies suggest RPA-based TB assays could be of use for integration into a point-of-care test for use in resource constrained settings.

Polymorphisms in innate immunity genes and risk of non-Hodgkin lymphoma.

Genetic variation in innate immunity may alter host‐pathogen defence mechanisms and promote aberrant immune responses and non‐Hodgkin lymphoma (NHL). To test this hypothesis, we investigated polymorphisms in innate immune genes in a pooled analysis of two population‐based case‐control studies of NHL from the San Francisco Bay Area (308 cases, 684 controls) and UK (596 cases, 758 controls). The caspase recruitment domain‐containing protein 1007fs homozygote variant was positively associated with NHL risk (odds ratios (OR) = 3·1, 95% confidence intervals (CI) 1·1–8·8), whereas the toll‐like receptor 4 1063A>G variant allele was inversely associated with diffuse large cell lymphoma (OR = 0·67, 95% CI 0·45–0·99). These results suggest that variation in innate immune genes may alter NHL susceptibility.

Polymorphisms in Cytokine and Cellular Adhesion Molecule Genes and Susceptibility to Hematotoxicity among Workers Exposed to Benzene

Benzene is a recognized hematotoxin and leukemogen but its mechanism of action and the role of genetic susceptibility are still unclear. Cytokines, chemokines, and cellular adhesion molecules are soluble proteins that play an important regulatory role in hematopoiesis. We therefore hypothesized that variation in these genes could influence benzene-induced hematotoxicity. We analyzed common, well-studied single-nucleotide polymorphisms (SNPs) in 20 candidate genes drawn from these pathways in a study of 250 workers exposed to benzene and 140 unexposed controls in China. After accounting for multiple comparisons, SNPs in five genes were associated with a statistically significant decrease in total WBC counts among exposed workers [IL-1A (-889C>T), IL-4 (-1098T>G), IL-10 (-819T>C), IL-12A (8685G>A), and VCAM1 (-1591T>C)], and one SNP [CSF3 (Ex4-165C>T)] was associated with an increase in WBC counts. The adhesion molecule VCAM1 variant was particularly noteworthy as it was associated with a decrease in B cells, natural killer cells, CD4+ T cells, and monocytes. Further, VCAM1 (-1591T>C) and CSF3 (Ex4-165C>T) were associated, respectively, with decreased (P = 0.041) and increased (P = 0.076) CFU-GEMM progenitor cell colony formation in 29 benzene-exposed workers. This is the first report to provide evidence that SNPs in genes that regulate hematopoiesis influence benzene-induced hematotoxicity.

Changes in the peripheral blood transcriptome associated with occupational benzene exposure identified by cross-comparison on two microarray platforms

Benzene is an established cause of leukemia, and possibly lymphoma, in humans, but the underlying molecular pathways remain largely undetermined. We used two microarray platforms to identify global gene expression changes associated with well-characterized occupational benzene exposure in the peripheral blood mononuclear cells (PBMC) of a population of shoe-factory workers. Differential expression of 2692 genes (Affymetrix) and 1828 genes (Illumina) was found and the concordance was 50% (based on an average fold-change > or =1.3 from the two platforms), with similar expression ratios among the concordant genes. Four genes (CXCL16, ZNF331, JUN and PF4), which we previously identified by microarray and confirmed by real-time PCR, were among the top 100 genes identified by both platforms in the current study. Gene ontology analysis showed overrepresentation of genes involved in apoptosis among the concordant genes while pathway analysis identified pathways related to lipid metabolism. The two-platform approach allows for robust changes in the PBMC transcriptome of benzene-exposed individuals to be identified.

Modifier Effects between Regulatory and Protein-Coding Variation

Genome-wide associations have shown a lot of promise in dissecting the genetics of complex traits in humans with single variants, yet a large fraction of the genetic effects is still unaccounted for. Analyzing genetic interactions between variants (epistasis) is one of the potential ways forward. We investigated the abundance and functional impact of a specific type of epistasis, namely the interaction between regulatory and protein-coding variants. Using genotype and gene expression data from the 210 unrelated individuals of the original four HapMap populations, we have explored the combined effects of regulatory and protein-coding single nucleotide polymorphisms (SNPs). We predict that about 18% (1,502 out of 8,233 nsSNPs) of protein-coding variants are differentially expressed among individuals and demonstrate that regulatory variants can modify the functional effect of a coding variant in cis. Furthermore, we show that such interactions in cis can affect the expression of downstream targets of the gene containing the protein-coding SNP. In this way, a cis interaction between regulatory and protein-coding variants has a trans impact on gene expression. Given the abundance of both types of variants in human populations, we propose that joint consideration of regulatory and protein-coding variants may reveal additional genetic effects underlying complex traits and disease and may shed light on causes of differential penetrance of known disease variants.

Recombinations in Staphylococcal Cassette Chromosome mec Elements Compromise the Molecular Detection of Methicillin Resistance in Staphylococcus aureus

Clinical laboratories are increasingly using molecular tests for methicillin-resistant Staphylococcus aureus (MRSA) screening. However, primers have to be targeted to a variable chromosomal region, the staphylococcal cassette chromosome mec (SCCmec). We initially screened 726 MRSA isolates from a single UK hospital trust by recombinase polymerase amplification (RPA), a novel, isothermal alternative to PCR. Undetected isolates were further characterised using multilocus sequence, spa typing and whole genome sequencing. 96% of our tested phenotypically MRSA isolates contained one of the six orfX-SCCmec junctions our RPA test and commercially available molecular tests target. However 30 isolates could not be detected. Sequencing of 24 of these isolates demonstrated recombinations within the SCCmec element with novel insertions that interfered with the RPA, preventing identification as MRSA. This result suggests that clinical laboratories cannot rely solely upon molecular assays to reliably detect all methicillin-resistance. The presence of significant recombinations in the SCCmec element, where the majority of assays target their primers, suggests that there will continue to be isolates that escape identification. We caution that dependence on amplification-based molecular assays will continue to result in failure to diagnose a small proportion (∼4%) of MRSA isolates, unless the true level of SCCmec natural diversity is determined by whole genome sequencing of a large collection of MRSA isolates.