Spatio-temporal predictions of COVID-19 test positivity in Uppsala County, Sweden: a comparative approach.

Previous spatio-temporal COVID-19 prediction models have focused on the prediction of subsequent number of cases, and have shown varying accuracy and lack of high geographical resolution. We aimed to predict trends in COVID-19 test positivity, an important marker for planning local testing capacity and accessibility. We included a full year of information (June 29, 2020-July 4, 2021) with both direct and indirect indicators of transmission, e.g. mobility data, number of calls to the national healthcare advice line and vaccination coverage from Uppsala County, Sweden, as potential predictors. We developed four models for a 1-week-window, based on gradient boosting (GB), random forest (RF), autoregressive integrated moving average (ARIMA) and integrated nested laplace approximations (INLA). Three of the models (GB, RF and INLA) outperformed the naïve baseline model after data from a full pandemic wave became available and demonstrated moderate accuracy. An ensemble model of these three models slightly improved the average root mean square error to 0.039 compared to 0.040 for GB, RF and INLA, 0.055 for ARIMA and 0.046 for the naïve model. Our findings indicate that the collection of a wide variety of data can contribute to spatio-temporal predictions of COVID-19 test positivity.

Publisher Correction: Transcriptomic alterations during ageing reflect the shift from cancer to degenerative diseases in the elderly.

The original version of this Article contained an error in the spelling of the author Jule Müller, which was incorrectly given as Julia Müller. Additionally, in Fig. 4a, the blue-red colour scale for fold change in ageing/disease regulation included a blue stripe in place of a red stripe at the right-hand end of the scale. These errors have been corrected in both the PDF and HTML versions of the Article.

Programmed Cell Death Ligand 1 Immunohistochemistry: A Concordance Study Between Surgical Specimen, Biopsy, and Tissue Microarray.

BACKGROUND: The immunohistochemical analysis of programmed cell death ligand 1 (PD-L1) expression in tumor tissue of non-small-cell lung cancer patients has now been integrated in the diagnostic workup. Analysis is commonly done on small tissue biopsy samples representing a minimal fraction of the whole tumor. The aim of the study was to evaluate the correlation of PD-L1 expression on biopsy specimens with corresponding resection specimens. MATERIALS AND METHODS: In total, 58 consecutive cases with preoperative biopsy and resected tumor specimens were selected. From each resection specimen 2 tumor cores were compiled into a tissue microarray (TMA). Immunohistochemical staining with the antibody SP263 was performed on biopsy specimens, resection specimens (whole sections), as well as on the TMA. RESULTS: The proportion of PD-L1-positive stainings were comparable between the resection specimens (48% and 19%), the biopsies (43% and 17%), and the TMAs (47% and 14%), using cutoffs of 1% and 50%, respectively (P > .39 all comparisons). When the resection specimens were considered as reference, PD-L1 status differed in 16%/5% for biopsies and in 9%/9% for TMAs (1%/50% cutoff). The sensitivity of the biopsy analysis was 79%/82% and the specificity was 90%/98% at the 1%/50% cutoff. The Cohens κ value for the agreement between biopsy and tumor. was 0.70 at the 1% cutoff and 0.83 at the 50% cutoff. CONCLUSION: The results indicate a moderate concordance between the analysis of biopsy and whole tumor tissue, resulting in misclassification of samples in particular when the lower 1% cutoff was used. Clinicians should be aware of this uncertainty when interpreting PD-L1 reports for treatment decisions.

Transcriptomic alterations during ageing reflect the shift from cancer to degenerative diseases in the elderly.

Disease epidemiology during ageing shows a transition from cancer to degenerative chronic disorders as dominant contributors to mortality in the old. Nevertheless, it has remained unclear to what extent molecular signatures of ageing reflect this phenomenon. Here we report on the identification of a conserved transcriptomic signature of ageing based on gene expression data from four vertebrate species across four tissues. We find that ageing-associated transcriptomic changes follow trajectories similar to the transcriptional alterations observed in degenerative ageing diseases but are in opposite direction to the transcriptomic alterations observed in cancer. We confirm the existence of a similar antagonism on the genomic level, where a majority of shared risk alleles which increase the risk of cancer decrease the risk of chronic degenerative disorders and vice versa. These results reveal a fundamental trade-off between cancer and degenerative ageing diseases that sheds light on the pronounced shift in their epidemiology during ageing.

Evaluation of Clonostachys rosea for Control of Plant-Parasitic Nematodes in Soil and in Roots of Carrot and Wheat.

Biological control is a promising approach to reduce plant diseases caused by nematodes. We tested the effect of the fungus Clonostachys rosea strain IK726 inoculation on nematode community composition in a naturally nematode infested soil in a pot experiment, and the effect of C. rosea on plant health. The numbers of plant-parasitic nematode genera extracted from soil and plant roots decreased by 40 to 73% when C. rosea was applied, while genera of nonparasitic nematodes were not affected. Soil inoculation of C. rosea increased fresh shoot weight and shoot length of wheat plants by 20 and 24%, respectively, while only shoot dry weight increased by 48% in carrots. Light microscopy of in vitro C. rosea-nematode interactions did not reveal evidence of direct parasitism. However, culture filtrates of C. rosea growing in potato dextrose broth, malt extract broth and synthetic nutrient broth exhibited toxicity toward nematodes and immobilized 57, 62, and 100% of the nematodes, respectively, within 48 h. This study demonstrates that C. rosea can control plant-parasitic nematodes and thereby improve plant growth. The most likely mechanism responsible for the antagonism is antibiosis through production of nematicidal compounds, rather than direct parasitism.

Global Transcriptional Response of Human Liver Cells to Ethanol Stress of Different Strength Reveals Hormetic Behavior.

BACKGROUND: The liver is the major site for alcohol metabolism in the body and therefore the primary target organ for ethanol (EtOH)-induced toxicity. In this study, we investigated the in vitro response of human liver cells to different EtOH concentrations in a perfused bioartificial liver device that mimics the complex architecture of the natural organ. METHODS: Primary human liver cells were cultured in the bioartificial liver device and treated for 24 hours with medium containing 150 mM (low), 300 mM (medium), or 600 mM (high) EtOH, while a control culture was kept untreated. Gene expression patterns for each EtOH concentration were monitored using Affymetrix Human Gene 1.0 ST Gene chips. Scaled expression profiles of differentially expressed genes (DEGs) were clustered using Fuzzy c-means algorithm. In addition, functional classification methods, KEGG pathway mapping and also a machine learning approach (Random Forest) were utilized. RESULTS: A number of 966 (150 mM EtOH), 1,334 (300 mM EtOH), or 4,132 (600 mM EtOH) genes were found to be differentially expressed. Dose-response relationships of the identified clusters of co-expressed genes showed a monotonic, threshold, or nonmonotonic (hormetic) behavior. Functional classification of DEGs revealed that low or medium EtOH concentrations operate adaptation processes, while alterations observed for the high EtOH concentration reflect the response to cellular damage. The genes displaying a hormetic response were functionally characterized by overrepresented "cellular ketone metabolism" and "carboxylic acid metabolism." Altered expression of the genes BAHD1 and H3F3B was identified as sufficient to classify the samples according to the applied EtOH doses. CONCLUSIONS: Different pathways of metabolic and epigenetic regulation are affected by EtOH exposition and partly undergo hormetic regulation in the bioartificial liver device. Gene expression changes observed at high EtOH concentrations reflect in some aspects the situation of alcoholic hepatitis in humans.

Conserved genes and pathways in primary human fibroblast strains undergoing replicative and radiation induced senescence.

BACKGROUND: Cellular senescence is induced either internally, for example by replication exhaustion and cell division, or externally, for example by irradiation. In both cases, cellular damages accumulate which, if not successfully repaired, can result in senescence induction. Recently, we determined the transcriptional changes combined with the transition into replicative senescence in primary human fibroblast strains. Here, by γ-irradiation we induced premature cellular senescence in the fibroblast cell strains (HFF and MRC-5) and determined the corresponding transcriptional changes by high-throughput RNA sequencing. RESULTS: Comparing the transcriptomes, we found a high degree of similarity in differential gene expression in replicative as well as in irradiation induced senescence for both cell strains suggesting, in each cell strain, a common cellular response to error accumulation. On the functional pathway level, "Cell cycle" was the only pathway commonly down-regulated in replicative and irradiation-induced senescence in both fibroblast strains, confirming the tight link between DNA repair and cell cycle regulation. However, "DNA repair" and "replication" pathways were down-regulated more strongly in fibroblasts undergoing replicative exhaustion. We also retrieved genes and pathways in each of the cell strains specific for irradiation induced senescence. CONCLUSION: We found the pathways associated with "DNA repair" and "replication" less stringently regulated in irradiation induced compared to replicative senescence. The strong regulation of these pathways in replicative senescence highlights the importance of replication errors for its induction.

Longitudinal RNA-Seq Analysis of Vertebrate Aging Identifies Mitochondrial Complex I as a Small-Molecule-Sensitive Modifier of Lifespan.

Mutations and genetic variability affect gene expression and lifespan, but the impact of variations in gene expression within individuals on their aging-related mortality is poorly understood. We performed a longitudinal study in the short-lived killifish, Nothobranchius furzeri, and correlated quantitative variations in gene expression during early adult life with lifespan. Shorter- and longer-lived individuals differ in their gene expression before the onset of aging-related mortality; differences in gene expression are more pronounced early in life. We identified mitochondrial respiratory chain complex I as a hub in a module of genes whose expression is negatively correlated with lifespan. Accordingly, partial pharmacological inhibition of complex I by the small molecule rotenone reversed aging-related regulation of gene expression and extended lifespan in N. furzeri by 15%. These results support the use of N. furzeri as a vertebrate model for identifying the protein targets, pharmacological modulators, and individual-to-individual variability associated with aging.

Low sulfide levels and a high degree of cystathionine ß-synthase (CBS) activation by S-adenosylmethionine (SAM) in the long-lived naked mole-rat.

Hydrogen sulfide (H2S) is a gaseous signalling molecule involved in many physiological and pathological processes. There is increasing evidence that H2S is implicated in aging and lifespan control in the diet-induced longevity models. However, blood sulfide concentration of naturally long-lived species is not known. Here we measured blood sulfide in the long-lived naked mole-rat and five other mammalian species considerably differing in lifespan and found a negative correlation between blood sulfide and maximum longevity residual. In addition, we show that the naked mole-rat cystathionine ß-synthase (CBS), an enzyme whose activity in the liver significantly contributes to systemic sulfide levels, has lower activity in the liver and is activated to a higher degree by S-adenosylmethionine compared to other species. These results add complexity to the understanding of the role of H2S in aging and call for detailed research on naked mole-rat transsulfuration.

Polymorphisms of cystathionine beta-synthase gene are associated with susceptibility to sepsis.

Sepsis is the systemic inflammatory host response to infection. Cystathionine beta-synthase (CBS)-dependent homocysteine (Hcy) pathway was demonstrated to affect disease severity and mortality in patients with severe sepsis/septic shock. Independent studies identified a single-nucleotide polymorphism (SNP, rs6586282, hg19 chr21:g.44478497C>T) in intron 14 of the CBS-coding gene (CBS) associated with Hcy plasma levels. We aimed to describe the association of this SNP and variants of a splice donor-affecting variable-number tandem repeat (VNTR, NG_008938.1:g.22763_22793[16_22]) 243 bp downstream of rs6586282 with severe human sepsis. We analyzed the VNTR structure and genotyped variants of rs6586282 and a neighboring SNP (rs34758144, hg19 chr21:g.44478582G>A) in two case-control studies including patients with severe sepsis/septic shock from Germany (n=168) and Greece (n=237). In both studies, we consistently observed an association of CBS VNTR alleles with sepsis susceptibility. Risk linearly increased with number of tandem repeats (per allele odds ratio in the adjusted analysis 1.34; 95% confidence interval (CI)=1.17-1.55; P<0.001). Association had also been shown for rs34758144 whose risk allele is in linkage disequilibrium with one long VNTR allele (19 repeat). In contrast, we observed no evidence for an effect on 28-day survival in patients with severe sepsis/septic shock (per allele hazard ratio in the adjusted analysis for VNTR 1.10; 95% CI=0.95-1.28; P=0.20). In a minigene approach, we demonstrated alternative splicing in distinct VNTR alleles, which, however, was independent of the number of tandem units. In conclusion, there is no ordinary conjunction between human CBS and severe sepsis/septic shock, but CBS genotypes are involved in disease susceptibility.

Conserved genes and pathways in primary human fibroblast strains undergoing replicative and radiation induced senescence

BACKGROUND: Cellular senescence is induced either internally, for example by replication exhaustion and cell division, or externally, for example by irradiation. In both cases, cellular damages accumulate which, if not successfully repaired, can result in senescence induction. Recently, we determined the transcriptional changes combined with the transition into replicative senescence in primary human fibroblast strains. Here, by γ-irradiation we induced premature cellular senescence in the fibroblast cell strains (HFF and MRC-5) and determined the corresponding transcriptional changes by high-throughput RNA sequencing. RESULTS: Comparing the transcriptomes, we found a high degree of similarity in differential gene expression in replicative as well as in irradiation induced senescence for both cell strains suggesting, in each cell strain, a common cellular response to error accumulation. On the functional pathway level, "Cell cycle" was the only pathway commonly down-regulated in replicative and irradiation-induced senescence in both fibroblast strains, confirming the tight link between DNA repair and cell cycle regulation. However, "DNA repair" and "replication" pathways were down-regulated more strongly in fibroblasts undergoing replicative exhaustion. We also retrieved genes and pathways in each of the cell strains specific for irradiation induced senescence. CONCLUSION: We found the pathways associated with "DNA repair" and "replication" less stringently regulated in irradiation induced compared to replicative senescence. The strong regulation of these pathways in replicative senescence highlights the importance of replication errors for its induction.

Alternative splicing of SMPD1 in human sepsis.

Acid sphingomyelinase (ASM or sphingomyelin phosphodiesterase, SMPD) activity engages a critical role for regulation of immune response and development of organ failure in critically ill patients. Beside genetic variation in the human gene encoding ASM (SMPD1), alternative splicing of the mRNA is involved in regulation of enzymatic activity. Here we show that the patterns of alternatively spliced SMPD1 transcripts are significantly different in patients with systemic inflammatory response syndrome and severe sepsis/septic shock compared to control subjects allowing discrimination of respective disease entity. The different splicing patterns might contribute to the better understanding of the pathophysiology of human sepsis.

Extension of life span by impaired glucose metabolism in Caenorhabditis elegans is accompanied by structural rearrangements of the transcriptomic network.

Glucose restriction mimicked by feeding the roundworm Caenorhabditis elegans with 2-deoxy-D-glucose (DOG) - a glucose molecule that lacks the ability to undergo glycolysis - has been found to increase the life span of the nematodes considerably. To facilitate understanding of the molecular mechanisms behind this life extension, we analyzed transcriptomes of DOG-treated and untreated roundworms obtained by RNA-seq at different ages. We found that, depending on age, DOG changes the magnitude of the expression values of about 2 to 24 percent of the genes significantly, although our results reveal that the gross changes introduced by DOG are small compared to the age-induced changes. We found that 27 genes are constantly either up- or down-regulated by DOG over the whole life span, among them several members of the cytochrome P450 family. The monotonic change with age of the temporal expression patterns of the genes was investigated, leading to the result that 21 genes reverse their monotonic behaviour under impaired glycolysis. Put simply, the DOG-treatment reduces the gross transcriptional activity but increases the interconnectedness of gene expression. However, a detailed analysis of network parameters discloses that the introduced changes differ remarkably between individual signalling pathways. We found a reorganization of the hubs of the mTOR pathway when standard diet is replaced by DOG feeding. By constructing correlation based difference networks, we identified those signalling pathways that are most vigorously changed by impaired glycolysis. Taken together, we have found a number of genes and pathways that are potentially involved in the DOG-driven extension of life span of C. elegans. Furthermore, our results demonstrate how the network structure of ageing-relevant signalling pathways is reorganised under impaired glycolysis.

Alternative 5' untranslated regions are involved in expression regulation of human heme oxygenase-1.

The single nucleotide polymorphism rs2071746 and a (GT)n microsatellite within the human gene encoding heme oxygenase-1 (HMOX1) are associated with incidence or outcome in a variety of diseases. Most of these associations involve either release of heme or oxidative stress. Both polymorphisms are localized in the promoter region, but previously reported correlations with heme oxygenase-1 expression remain not coherent. This ambiguity suggests a more complex organization of the 5' gene region which we sought to investigate more fully. We evaluated the 5' end of HMOX1 and found a novel first exon 1a placing the two previously reported polymorphisms in intronic or exonic positions within the 5' untranslated region respectively. Expression of exon 1a can be induced in HepG2 hepatoma cells by hemin and is a repressor of heme oxygenase-1 translation as shown by luciferase reporter assays. Moreover, minigene approaches revealed that the quantitative outcome of alternative splicing within the 5' untranslated region is affected by the (GT)n microsatellite. This data supporting an extended HMOX1 gene model and provide further insights into expression regulation of heme oxygenase-1. Alternative splicing within the HMOX1 5' untranslated region contributes to translational regulation and is a mechanistic feature involved in the interplay between genetic variations, heme oxygenase-1 expression and disease outcome.

Constant splice-isoform ratios in human lymphoblastoid cells support the concept of a splico-stat.

Splicing generates mature transcripts from genes in pieces in eukaryotic cells. Overwhelming evidence has accumulated that alternative routes in splicing are possible for most human and mammalian genes, thereby allowing formation of different transcripts from one gene. No function has been assigned to the majority of identified alternative splice forms, and it has been assumed that they compose inert or tolerated waste from aberrant or noisy splicing. Here we demonstrate that five human transcription units (WT1, NOD2, GNAS, RABL2A, RABL2B) have constant splice-isoform ratios in genetically diverse lymphoblastoid cell lines independent of the type of alternative splicing (exon skipping, alternative donor/acceptor, tandem splice sites) and gene expression level. Even splice events that create premature stop codons and potentially trigger nonsense-mediated mRNA decay are found at constant fractions. The analyzed alternative splicing events were qualitatively but not quantitatively conserved in corresponding chimpanzee cell lines. Additionally, subtle splicing at tandem acceptor splice sites (GNAS, RABL2A/B) was highly constrained and strongly depends on the upstream donor sequence content. These results also demonstrate that unusual and unproductive splice variants are produced in a regulated manner.

Recurrent genomic alterations in benign and malignant pheochromocytomas and paragangliomas revealed by whole-genome array comparative genomic hybridization analysis.

Pheochromocytomas and abdominal paragangliomas are adrenal and extra-adrenal catecholamine-producing tumours. They arise due to heritable cancer syndromes, or more frequently occur sporadically due to an unknown genetic cause. The majority of cases are benign, but malignant tumours are observed. Previous comparative genomic hybridization (CGH) and loss of heterozygosity studies have shown frequent deletions of chromosome arms 1p, 3q and 22q in pheochromocytomas. We applied high-resolution whole-genome array CGH on 53 benign and malignant pheochromocytomas and paragangliomas to narrow down candidate regions as well as to identify chromosomal alterations more specific to malignant tumours. Minimal overlapping regions (MORs) were identified on 16 chromosomes, with the most frequent MORs of deletion (> or = 32%) occurring on chromosome arms 1p, 3q, 11p/q, 17p and 22q, while the chromosome arms 1q, 7p, 12q and 19p harboured the most common MORs of gain (> or = 14%). The most frequent MORs (61-75%) in the pheochromocytomas were identified at 1p, and the four regions of common losses encompassed 1p36, 1p32-31, 1p22-21 and 1p13. Tumours that did not show 1p loss generally demonstrated aberrations on chromosome 11. Gain of chromosomal material was significantly more frequent among the malignant cases. Moreover, gain at 19q, trisomy 12 and loss at 11q were positively associated with malignant pheochromocytomas, while 1q gain was commonly observed in the malignant paragangliomas. Our study revealed novel and narrow recurrent chromosomal regions of loss and gain at several autosomes, a prerequisite for identifying candidate tumour suppressor genes and oncogenes involved in the development of adrenal and extra-adrenal catecholamine-producing tumours.

Frequent genetic differences between matched primary and metastatic breast cancer provide an approach to identification of biomarkers for disease progression.

Breast cancer is a major cause of morbidity and mortality in women and its metastatic spread is the principal reason behind the fatal outcome. Metastasis-related research of breast cancer is however underdeveloped when compared with the abundant literature on primary tumors. We applied an unexplored approach comparing at high resolution the genomic profiles of primary tumors and synchronous axillary lymph node metastases from 13 patients with breast cancer. Overall, primary tumors displayed 20% higher number of aberrations than metastases. In all but two patients, we detected in total 157 statistically significant differences between primary lesions and matched metastases. We further observed differences that can be linked to metastatic disease and there was also an overlapping pattern of changes between different patients. Many of the differences described here have been previously linked to poor patient survival, suggesting that this is a viable approach toward finding biomarkers for disease progression and definition of new targets useful for development of anticancer drugs. Frequent genetic differences between primary tumors and metastases in breast cancer also question, at least to some extent, the role of primary tumors as a surrogate subject of study for the systemic disease.

Focal amplifications are associated with high grade and recurrences in stage Ta bladder carcinoma.

Urinary bladder cancer is a heterogeneous disease with tumors ranging from papillary noninvasive (stage Ta) to solid muscle infiltrating tumors (stage T2+). The risk of progression and death for the most frequent diagnosed type, Ta, is low, but the high incidence of recurrences has a significant effect on the patients' quality of life and poses substantial costs for health care systems. Consequently, the purpose of this study was to search for predictive factors of recurrence on the basis of genetic profiling. A clinically well characterized cohort of Ta bladder carcinomas, selected by the presence or absence of recurrences, was evaluated by an integrated analysis of DNA copy number changes and gene expression (clone-based 32K, respectively, U133Plus2.0 arrays). Only a few chromosomal aberrations have previously been defined in superficial bladder cancer. Surprisingly, the profiling of Ta tumors with a high-resolution array showed that DNA copy alterations are relatively common in this tumor type. Furthermore, we observed an overrepresentation of focal amplifications within high-grade and recurrent cases. Known (FGFR3, CCND1, MYC, MDM2) and novel candidate genes were identified within the loci. For example, MYBL2, a nuclear transcription factor involved in cell-cycle progression; YWHAB, an antiapoptotic protein; and SDC4, an important component of focal adhesions represent interesting candidates detected within two amplicons on chromosome 20, for which DNA amplification correlated with transcript up-regulation. The observed overrepresentation of amplicons within high-grade and recurrent cases may be clinically useful for the identification of patients who will benefit from a more aggressive therapy.

Genome-wide high-resolution analysis of DNA copy number alterations in NF1-associated malignant peripheral nerve sheath tumors using 32K BAC array.

Neurofibromatosis Type I (NF1) is an autosomal dominant disorder characterized by the development of both benign and malignant tumors. The lifetime risk for developing a malignant peripheral nerve sheath tumor (MPNST) in NF1 patients is approximately 10% with poor survival rates. To date, the molecular basis of MPNST development remains unclear. Here, we report the first genome-wide and high-resolution analysis of DNA copy number alterations in MPNST using the 32K bacterial artificial chromosome microarray on a series of 24 MPNSTs and three neurofibroma samples. In the benign neurofibromas, apart from loss of one copy of the NF1 gene and copy number polymorphisms, no other changes were found. The profiles of malignant samples, however, revealed specific loss of chromosomal regions including 1p35-33, 1p21, 9p21.3, 10q25, 11q22-23, 17q11, and 20p12.2 as well as gain of 1q25, 3p26, 3q13, 5p12, 5q11.2-q14, 5q21-23, 5q31-33, 6p23-p21, 6p12, 6q15, 6q23-q24, 7p22, 7p14-p13, 7q21, 7q36, 8q22-q24, 14q22, and 17q21-q25. Copy number gains were more frequent than deletions in the MPNST samples (62% vs. 38%). The genes resident within common regions of gain were NEDL1 (7p14), AP3B1 (5q14.1), and CUL1 (7q36.1) and these were identified in >63% MPNSTs. The most frequently deleted locus encompassed CDKN2A, CDKN2B, and MTAP genes on 9p21.3 (33% cases). These genes have previously been implicated in other cancer conditions and therefore, should be considered for their therapeutic, prognostic, and diagnostic relevance in NF1 tumorigenesis.

Tissue-specific variation in DNA methylation levels along human chromosome 1.

BACKGROUND: DNA methylation is a major epigenetic modification important for regulating gene expression and suppressing spurious transcription. Most methods to scan the genome in different tissues for differentially methylated sites have focused on the methylation of CpGs in CpG islands, which are concentrations of CpGs often associated with gene promoters. RESULTS: Here, we use a methylation profiling strategy that is predominantly responsive to methylation differences outside of CpG islands. The method compares the yield from two samples of size-selected fragments generated by a methylation-sensitive restriction enzyme. We then profile nine different normal tissues from two human donors relative to spleen using a custom array of genomic clones covering the euchromatic portion of human chromosome 1 and representing 8% of the human genome. We observe gross regional differences in methylation states across chromosome 1 between tissues from the same individual, with the most striking differences detected in the comparison of cerebellum and spleen. Profiles of the same tissue from different donors are strikingly similar, as are the profiles of different lobes of the brain. Comparing our results with published gene expression levels, we find that clones exhibiting extreme ratios reflecting low relative methylation are statistically enriched for genes with high expression ratios, and vice versa, in most pairs of tissues examined. CONCLUSION: The varied patterns of methylation differences detected between tissues by our methylation profiling method reinforce the potential functional significance of regional differences in methylation levels outside of CpG islands.