Genomic Signatures for Sedimentary Microbial Utilization of Phytoplankton Detritus in a Fast-Flowing Estuary.

In fast-flowing, river-dominated estuaries, "hotspots" of microbial biogeochemical cycling can be found within areas of extended water retention. Lateral bays located off of the North and South channels of the Columbia River estuary are proposed to be such hotspots. Previous metagenomic studies on water samples indicated that these regions function both as sources and sinks of biogenic particles, with potential to impact organic matter fluxes in the estuary. To extend this work, we analyzed 11 sediment metagenomes from three disparate bays: the freshwater Cathlamet Bay, and the brackish Youngs Bay and more saline Baker Bay located nearer the mouth to the south and north of the main channel, respectively. Samples were collected from upper layers of sediments in August of 2011 and 2013 for DNA extraction and metagenome sequencing. All metagenomes were dominated by bacterial sequences, although diatom sequences as high as 26% of the total annotated sequences were observed in the higher salinity samples. Unsupervised 2D hierarchical clustering analysis resulted in the eleven metagenome samples clustered into four groups by microbial taxonomic composition, with Bacteroides, diatom, and phage levels driving most of the grouping. Results of functional gene clustering further indicated that diatom bloom degradation stage (early vs. late) was an important factor. While the Flavobacteriia and Cytophagia classes were well represented in metagenomes containing abundant diatoms, taxa from the Bacteroidia class, along with certain members of the Sphingobacteriia class, were particularly abundant in metagenomes representing later stages of diatom decomposition. In contrast, the sediment metagenomes with low relative abundance of diatom and Bacteroidetes sequences appeared to have a metabolic potential biased toward microbial growth under nutrient limitation. While differences in water salinity clearly also influenced the microbial community composition and metabolic potential, our results highlight a central role for allochthonous labile organic matter (i.e., diatom detritus), in shaping bacterial taxonomic and functional properties in the Columbia River estuary lateral bay sediments. These results suggest that in fast-flowing, river-dominated estuaries, sediment microbial communities in areas of extended water retention, such as the lateral bays, may contribute disproportionately to estuarine organic matter degradation and recycling.

Detection of sub-microscopic blood levels of Plasmodium falciparum using Tandem Oligonucleotide Repeat Cascade Amplification (TORCA) assay with an attomolar detection limit.

Tandem Oligonucleotide Repeat Cascade Amplification (TORCA) based on signal rather than target amplification under isothermal conditions was developed for nucleic acid assays. The initial signal was generated by hybridization of single stranded DNA targets to immobilized recognition probes followed by hybrid cleavage with specific restriction endonuclease (REase), and release of trigger oligonucleotides (Tr1). The signal amplification chamber contained two bead types carrying single-stranded amplification probes and two amplification REases. The probes consisted of multiple tandem repeats of either Tr1 or another trigger Tr2, with the tandem-Tr1 anchored to the beads through the antisense Tr2 linker and vice versa. Addition of the recognition reaction solution and Tr1 hybridization to the anti-Tr1 linkers started cleavage and release of additional Tr1 and Tr2, resulting in exponential signal amplification. The cleavage cascade also released horseradish peroxidase (HRP) pre-attached to the amplification probes, and the resultant signal was measured colorimetrically. A TORCA assay was developed for detection of Plasmodium falciparum parasites in blood. It had the detection limit in the attomolar concentration range, successfully detecting sub-microscopic P. falciparum infections at less than 0.75 infected erythrocytes per microliter. Further TORCA optimization will likely produce the quantitative isothermal alternative to PCR at a fraction of its cost.

Host-specific response to HCV infection in the chimeric SCID-beige/Alb-uPA mouse model: role of the innate antiviral immune response.

The severe combined immunodeficiency disorder (SCID)-beige/albumin (Alb)-urokinase plasminogen activator (uPA) mouse containing a human-mouse chimeric liver is currently the only small animal model capable of supporting hepatitis C virus (HCV) infection. This model was utilized to characterize the host transcriptional response to HCV infection. The purpose of these studies was to investigate the genetic component of the host response to HCV infection and also to distinguish virus-induced gene expression changes from adaptive HCV-specific immune-mediated effects. Gene expression profiles from HCV-infected mice were also compared to those from HCV-infected patients. Analyses of the gene expression data demonstrate that host factors regulate the response to HCV infection, including the nature of the innate antiviral immune response. They also indicate that HCV mediates gene expression changes, including regulation of lipid metabolism genes, which have the potential to be directly cytopathic, indicating that liver pathology may not be exclusively mediated by HCV-specific adaptive immune responses. This effect appears to be inversely related to the activation of the innate antiviral immune response. In summary, the nature of the initial interferon response to HCV infection may determine the extent of viral-mediated effects on host gene expression.

CombiMatrix oligonucleotide arrays: genotyping and gene expression assays employing electrochemical detection.

Electrochemical detection has been developed and assay performances studied for the CombiMatrix oligonucleotide microarray platform that contains 12,544 individually addressable microelectrodes (features) in a semiconductor matrix. The approach is based on the detection of redox active chemistries (such as horseradish peroxidase (HRP) and the associated substrate TMB) proximal to specific microarray electrodes. First, microarray probes are hybridized to biotin-labeled targets, second, the HRP-streptavidin conjugate binds to biotin, and enzymatic oxidation of the electron donor substrate then occurs. The detection current is generated due to electro-reduction of the HRP reaction product, and it is measured with the CombiMatrix ElectraSense Reader. Performance of the ElectraSense platform has been characterized using gene expression and genotyping assays to analyze: (i) signal to concentration dependence, (ii) assay resolution, (iii) coefficients of variation, (CV) and (iv) array-to-array reproducibility and data correlation. The ElectraSense platform was also compared to the standard fluorescent detection, and good consistency was observed between these two different detection techniques. A lower detection limit of 0.75 pM was obtained for ElectraSense as compared to the detection limit of 1.5 pM obtained for fluorescent detection. Thus, the ElectraSense platform has been used to develop nucleic acid assays for highly accurate genotyping of a variety of pathogens including bio-threat agents (such as Bacillus anthracis, Yersinia pestis, and other microorganisms including Escherichia coli, Bacillus subtilis, etc.) and common pathogens of the respiratory tract (e.g. influenza A virus).

Microbial players and processes involved in phytoplankton bloom utilization in the water column of a fast-flowing, river-dominated estuary.

Fueled by seasonal phytoplankton blooms, the Columbia River estuary is a natural bioreactor for organic matter transformations. Prior metagenome analyses indicated high abundances of diverse Bacteroidetes taxa in estuarine samples containing phytoplankton. To examine the hypothesis that Bacteroidetes taxa have important roles in phytoplankton turnover, we further analyzed metagenomes from water collected along a salinity gradient at 0, 5, 15, 25, and 33 PSU during bloom events. Size fractions were obtained by using a 3-µm prefilter and 0.2-µm collection filter. Although this approach targeted bacteria by removing comparatively large eukaryotic cells, the metagenome from the ES-5 sample (5 PSU) nevertheless contained an abundance of diatom DNA. Biogeochemical measurements and prior studies indicated that this finding resulted from the leakage of cellular material due to freshwater diatom lysis at low salinity. Relative to the other metagenomes, the bacterial fraction of ES-5 was dramatically depleted of genes annotated as Bacteroidetes and lysogenic bacteriophages, but was overrepresented in DNA of protists and Myxococcales bacterivores. We suggest the following equally plausible scenarios for the microbial response to phytoplankton lysis: (1) Bacteroidetes depletion in the free-living fraction may at least in part be caused by their attachment to fluvial diatoms as the latter are lysed upon contact with low-salinity estuarine waters; (2) diatom particle colonization is likely followed by rapid bacterial growth and lytic phage infection, resulting in depletion of lysogenic bacteriophages and host bacteria; and (3) the subsequent availability of labile organic matter attracted both grazers and predators to feed in this estuarine biogeochemical "hotspot," which may have additionally depleted Bacteroidetes populations. These results represent the first detailed molecular analysis of the microbial response to phytoplankton lysis at the freshwater-brackish water interface in the fast-flowing Columbia River estuary.

Application of functional genomics to the chimeric mouse model of HCV infection: optimization of microarray protocols and genomics analysis.

BACKGROUND: Many model systems of human viral disease involve human-mouse chimeric tissue. One such system is the recently developed SCID-beige/Alb-uPA mouse model of hepatitis C virus (HCV) infection which involves a human-mouse chimeric liver. The use of functional genomics to study HCV infection in these chimeric tissues is complicated by the potential cross-hybridization of mouse mRNA on human oligonucleotide microarrays. To identify genes affected by mouse liver mRNA hybridization, mRNA from identical human liver samples labeled with either Cy3 or Cy5 was compared in the presence and absence of known amounts of mouse liver mRNA labeled in only one dye. RESULTS: The results indicate that hybridization of mouse mRNA to the corresponding human gene probe on Agilent Human 22 K oligonucleotide microarray does occur. The number of genes affected by such cross-hybridization was subsequently reduced to approximately 300 genes both by increasing the hybridization temperature and using liver samples which contain at least 80% human tissue. In addition, Real Time quantitative RT-PCR using human specific probes was shown to be a valid method to verify the expression level in human cells of known cross-hybridizing genes. CONCLUSION: The identification of genes affected by cross-hybridization of mouse liver RNA on human oligonucleotide microarrays makes it feasible to use functional genomics approaches to study the chimeric SCID-beige/Alb-uPA mouse model of HCV infection. This approach used to study cross-species hybridization on oligonucleotide microarrays can be adapted to other chimeric systems of viral disease to facilitate selective analysis of human gene expression.

Identification of novel tumor markers in hepatitis C virus-associated hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is a common primary cancer associated frequently with hepatitis C virus (HCV). To gain insight into the molecular mechanisms of hepatocarcinogenesis, and to identify potential HCC markers, we performed cDNA microarray analysis on surgical liver samples from 20 HCV-infected patients. RNA from individual tumors was compared with RNA isolated from adjacent nontumor tissue that was cirrhotic in all of the cases. Gene expression changes related to cirrhosis were filtered out using experiments in which pooled RNA from HCV-infected cirrhotic liver without tumors was compared with pooled RNA from normal liver. Expression of approximately 13,600 genes was analyzed using the advanced analysis tools of the Rosetta Resolver System. This analysis revealed a set of 50 potential HCC marker genes, which were up-regulated in the majority of the tumors analyzed, much more widely than common clinical markers such as cell proliferation-related genes. This HCC marker set contained several cancer-related genes, including serine/threonine kinase 15 (STK15), which has been implicated in chromosome segregation abnormalities but which has not been linked previously with liver cancer. In addition, a set of genes encoding secreted or plasma proteins was identified, including plasma glutamate carboxypeptidase (PGCP) and two secreted phospholipases A2 (PLA2G13 and PLA2G7). These genes may provide potential HCC serological markers because of their strong up-regulation in more than half of the tumors analyzed. Thus, high throughput methods coupled with high-order statistical analyses may result in the development of new diagnostic tools for liver malignancies.

Metagenomic evidence for reciprocal particle exchange between the mainstem estuary and lateral bay sediments of the lower Columbia River.

Lateral bays of the lower Columbia River estuary are areas of enhanced water retention that influence net ecosystem metabolism through activities of their diverse microbial communities. Metagenomic characterization of sediment microbiota from three disparate sites in two brackish lateral bays (Baker and Youngs) produced ∼100 Gbp of DNA sequence data analyzed subsequently for predicted SSU rRNA and peptide-coding genes. The metagenomes were dominated by Bacteria. A large component of Eukaryota was present in Youngs Bay samples, i.e., the inner bay sediment was enriched with the invasive New Zealand mudsnail, Potamopyrgus antipodarum, known for high ammonia production. The metagenome was also highly enriched with an archaeal ammonia oxidizer closely related to Nitrosoarchaeum limnia. Combined analysis of sequences and continuous, high-resolution time series of biogeochemical data from fixed and mobile platforms revealed the importance of large-scale reciprocal particle exchanges between the mainstem estuarine water column and lateral bay sediments. Deposition of marine diatom particles in sediments near Youngs Bay mouth was associated with a dramatic enrichment of Bacteroidetes (58% of total Bacteria) and corresponding genes involved in phytoplankton polysaccharide degradation. The Baker Bay sediment metagenome contained abundant Archaea, including diverse methanogens, as well as functional genes for methylotrophy and taxonomic markers for syntrophic bacteria, suggesting that active methane cycling occurs at this location. Our previous work showed enrichments of similar anaerobic taxa in particulate matter of the mainstem estuarine water column. In total, our results identify the lateral bays as both sources and sinks of biogenic particles significantly impacting microbial community composition and biogeochemical activities in the estuary.

Restriction Cascade Exponential Amplification (RCEA) assay with an attomolar detection limit: a novel, highly specific, isothermal alternative to qPCR.

An alternative to qPCR was developed for nucleic acid assays, involving signal rather than target amplification. The new technology, Restriction Cascade Exponential Amplification (RCEA), relies on specific cleavage of probe-target hybrids by restriction endonucleases (REase). Two mutant REases for amplification (Ramp), S17C BamHI and K249C EcoRI, were conjugated to oligonucleotides, and immobilized on a solid surface. The signal generation was based on: (i) hybridization of a target DNA to a Ramp-oligonucleotide probe conjugate, followed by (ii) specific cleavage of the probe-target hybrid using a non-immobilized recognition REase. The amount of Ramp released into solution upon cleavage was proportionate to the DNA target amount. Signal amplification was achieved through catalysis, by the free Ramp, of a restriction cascade containing additional oligonucleotide-conjugated Ramp and horseradish peroxidase (HRP). Colorimetric quantification of free HRP indicated that the RCEA achieved a detection limit of 10 aM (10(-17) M) target concentration, or approximately 200 molecules, comparable to the sensitivity of qPCR-based assays. The RCEA assay had high specificity, it was insensitive to non-specific binding, and detected target sequences in the presence of foreign DNA. RCEA is an inexpensive isothermal assay that allows coupling of the restriction cascade signal amplification with any DNA target of interest.

Metagenomic insights into particles and their associated microbiota in a coastal margin ecosystem.

Our previously published research was one of the pioneering studies on the use of metagenomics to directly compare taxonomic and metabolic properties of aquatic microorganisms from different filter size-fractions. We compared size-fractionated water samples representing free-living and particle-attached communities from four diverse habitats in the Columbia River coastal margin, analyzing 12 metagenomes consisting of >5 million sequence reads (>1.6 Gbp). With predicted peptide and rRNA data we evaluated eukaryotic, bacterial and archaeal populations across size fractions and related their properties to attached and free-living lifestyles, and their potential roles in carbon and nutrient cycling. In this focused review, we expand our discussion on the use of high-throughput sequence data to relate microbial community structure and function to the origin, fate and transport of particulate organic matter (POM) in coastal margins. We additionally discuss the potential impact of the priming effect on organic matter cycling at the land-ocean interface, and build a case for the importance, in particle-rich estuaries and coastal margin waters, of microbial activities in low-oxygen microzones within particle interiors.

A new restriction endonuclease-based method for highly-specific detection of DNA targets from methicillin-resistant Staphylococcus aureus.

PCR multiplexing has proven to be challenging, and thus has provided limited means for pathogen genotyping. We developed a new approach for analysis of PCR amplicons based on restriction endonuclease digestion. The first stage of the restriction enzyme assay is hybridization of a target DNA to immobilized complementary oligonucleotide probes that carry a molecular marker, horseradish peroxidase (HRP). At the second stage, a target-specific restriction enzyme is added, cleaving the target-probe duplex at the corresponding restriction site and releasing the HRP marker into solution, where it is quantified colorimetrically. The assay was tested for detection of the methicillin-resistant Staphylococcus aureus (MRSA) pathogen, using the mecA gene as a target. Calibration curves indicated that the limit of detection for both target oligonucleotide and PCR amplicon was approximately 1 nM. Sequences of target oligonucleotides were altered to demonstrate that (i) any mutation of the restriction site reduced the signal to zero; (ii) double and triple point mutations of sequences flanking the restriction site reduced restriction to 50-80% of the positive control; and (iii) a minimum of a 16-bp target-probe dsDNA hybrid was required for significant cleavage. Further experiments showed that the assay could detect the mecA amplicon from an unpurified PCR mixture with detection limits similar to those with standard fluorescence-based qPCR. Furthermore, addition of a large excess of heterologous genomic DNA did not affect amplicon detection. Specificity of the assay is very high because it involves two biorecognition steps. The proposed assay is low-cost and can be completed in less than 1 hour. Thus, we have demonstrated an efficient new approach for pathogen detection and amplicon genotyping in conjunction with various end-point and qPCR applications. The restriction enzyme assay may also be used for parallel analysis of multiple different amplicons from the same unpurified mixture in broad-range PCR applications.

Contrasting genomic properties of free-living and particle-attached microbial assemblages within a coastal ecosystem.

The Columbia River (CR) is a powerful economic and environmental driver in the US Pacific Northwest. Microbial communities in the water column were analyzed from four diverse habitats: (1) an estuarine turbidity maximum (ETM), (2) a chlorophyll maximum of the river plume, (3) an upwelling-associated hypoxic zone, and (4) the deep ocean bottom. Three size fractions, 0.1-0.8, 0.8-3, and 3-200 µm were collected for each habitat in August 2007, and used for DNA isolation and 454 sequencing, resulting in 12 metagenomes of >5 million reads (>1.6 Gbp). To characterize the dominant microorganisms and metabolisms contributing to coastal biogeochemistry, we used predicted peptide and rRNA data. The 3- and 0.8-µm metagenomes, representing particulate fractions, were taxonomically diverse across habitats. The 3-µm size fractions contained a high abundance of eukaryota with diatoms dominating the hypoxic water and plume, while cryptophytes were more abundant in the ETM. The 0.1-µm metagenomes represented mainly free-living bacteria and archaea. The most abundant archaeal hits were observed in the deep ocean and hypoxic water (19% of prokaryotic peptides in the 0.1-µm metagenomes), and were homologous to Nitrosopumilus maritimus (ammonia-oxidizing Thaumarchaeota). Bacteria dominated metagenomes of all samples. In the euphotic zone (estuary, plume and hypoxic ocean), the most abundant bacterial taxa (≥40% of prokaryotic peptides) represented aerobic photoheterotrophs. In contrast, the low-oxygen, deep water metagenome was enriched with sequences for strict and facultative anaerobes. Interestingly, many of the same anaerobic bacterial families were enriched in the 3-µm size fraction of the ETM (2-10X more abundant relative to the 0.1-µm metagenome), indicating possible formation of anoxic microniches within particles. Results from this study provide a metagenome perspective on ecosystem-scale metabolism in an upwelling-influenced river-dominated coastal margin.

Development of real-time assays for impedance-based detection of microbial double-stranded DNA targets: optimization and data analysis.

A real-time, label free assay was developed for microbial detection, utilizing double-stranded DNA targets and employing the next generation of an impedimetric sensor array platform designed by Sharp Laboratories of America (SLA). Real-time curves of the impedimetric signal response were obtained at fixed frequency and voltage for target binding to oligonucleotide probes attached to the sensor array surface. Kinetic parameters of these curves were analyzed by the integrated data analysis package for signal quantification. Non-specific binding presented a major challenge for assay development, and required assay optimization. For this, differences were maximized between binding curve kinetic parameters for probes binding to complementary targets versus non-target controls. Variables manipulated for assay optimization included target concentration, hybridization temperature, buffer concentration, and the use of surfactants. Our results showed that (i) different target-probe combinations required optimization of specific sets of variables; (ii) for each assay condition, the optimum range was relatively narrow, and had to be determined empirically; and (iii) outside of the optimum range, the assay could not distinguish between specific and non-specific binding. For each target-probe combination evaluated, conditions resulting in good separation between specific and non-specific binding signals were established, generating high confidence in the SLA impedimetric dsDNA assay results.

Seasonal changes in bacterial and archaeal gene expression patterns across salinity gradients in the Columbia River coastal margin.

Through their metabolic activities, microbial populations mediate the impact of high gradient regions on ecological function and productivity of the highly dynamic Columbia River coastal margin (CRCM). A 2226-probe oligonucleotide DNA microarray was developed to investigate expression patterns for microbial genes involved in nitrogen and carbon metabolism in the CRCM. Initial experiments with the environmental microarrays were directed toward validation of the platform and yielded high reproducibility in multiple tests. Bioinformatic and experimental validation also indicated that >85% of the microarray probes were specific for their corresponding target genes and for a few homologs within the same microbial family. The validated probe set was used to query gene expression responses by microbial assemblages to environmental variability. Sixty-four samples from the river, estuary, plume, and adjacent ocean were collected in different seasons and analyzed to correlate the measured variability in chemical, physical and biological water parameters to differences in global gene expression profiles. The method produced robust seasonal profiles corresponding to pre-freshet spring (April) and late summer (August). Overall relative gene expression was high in both seasons and was consistent with high microbial abundance measured by total RNA, heterotrophic bacterial production, and chlorophyll a. Both seasonal patterns involved large numbers of genes that were highly expressed relative to background, yet each produced very different gene expression profiles. April patterns revealed high differential gene expression in the coastal margin samples (estuary, plume and adjacent ocean) relative to freshwater, while little differential gene expression was observed along the river-to-ocean transition in August. Microbial gene expression profiles appeared to relate, in part, to seasonal differences in nutrient availability and potential resource competition. Furthermore, our results suggest that highly-active particle-attached microbiota in the Columbia River water column may perform dissimilatory nitrate reduction (both dentrification and DNRA) within anoxic particle microniches.

Identification of a specific gene expression pattern associated with HCV-induced pathogenesis in HCV- and HCV/HIV-infected individuals.

Gene expression profiling was performed on liver biopsies from 28 patients (12 HCV and 16 HCV/HIV infected) in an attempt to understand the mechanisms of HCV liver disease in the presence and absence of HIV coinfection. The data were compared with clinical observations and a gene expression database obtained for transplant HCV-infected samples. This is the first report of functional genomics being used to compare intrahepatic gene expression profiles of HCV- and HCV/HIV-infected individuals. Significantly, the intrahepatic global gene expression profiles do not differ between HCV- and HCV/HIV-infected individuals. However, a subset of patients was identified who share a specific pattern of gene expression, termed the enhanced gene expression (EGE) pattern. Specifically, the EGE (+) patients show a dramatic decreased expression of multiple genes associated with the FAS-apoptosis pathway and increased expression of lymphocyte adhesion molecules and lymphocyte-specific genes. The EGE (+) patients also have partially impaired Type I and II IFN-mediated antiviral responses, including a lack of induction of the anti-fibrogenic cytokine IFN-gamma. Importantly, the pattern of gene expression observed in EGE (+) patients has similarities to patients who developed fibrosis within 1 year of receiving a liver transplant.

Gene expression patterns that correlate with hepatitis C and early progression to fibrosis in liver transplant recipients.

BACKGROUND & AIMS: Liver transplant recipients infected with hepatitis C virus (HCV) develop recurrent hepatitis soon after transplantation and, in some cases, progress to fibrosis within the first 2 years. Our goals were to identify molecular processes influencing the liver disease progression and to find potential gene markers of early fibrosis. METHODS: We performed gene expression profiling on serial liver biopsy specimens obtained from 13 (11 infected and 2 uninfected) transplant recipients within the first year after transplantation at 0, 3, 6, and 12 months. The data were compared with clinical observations and with a gene expression database obtained for 55 nontransplant HCV-infected and uninfected liver samples. RESULTS: We identified several specific gene expression patterns. The first pattern was unique for the transplant recipients regardless of their infection status. The corresponding genes encoded stress response proteins and blood proteins involved in coagulation that were differentially expressed in response to posttransplantation graft recovery. The second pattern was specific to HCV-infected samples and included up-regulation of genes encoding components of the interferon-mediated antiviral response and immune system (antigen presentation, cytotoxic response). This up-regulation pattern was absent or suppressed in the patients who developed early fibrosis, indicating that the disease progression might result from an impaired liver response to infection. Finally, we identified gene expression patterns that were specific for 12-month biopsy specimens in all 4 HCV-infected patients who developed early fibrosis. CONCLUSIONS: The identified gene expression patterns may prove useful for diagnostic and prognostic applications in HCV-infected patients, including predicting early progression to fibrosis.

Hepatitis C virus and liver disease: global transcriptional profiling and identification of potential markers.

Microarray analysis of RNA from hepatitis C virus (HCV)-infected cirrhotic livers was performed to identify a gene expression signature of liver disease. The expression levels of approximately 13600 genes were analyzed using surgical material and core biopsy specimens from HCV-infected cirrhotic liver explants in comparison with reference samples of normal nondiseased liver. In addition, normal liver samples were compared with each other to determine normal physiologic variation in gene expression. A set of genes, including some associated with stress, acute-phase immune response, and hepatic stellate cell activation, had variable expression levels in normal livers. These genes were subtracted from the sets of genes differentially expressed in cirrhotic livers. To exclude cancer-related genes from our marker sets, we subtracted genes that also were expressed differentially in hepatocellular carcinomas. The resultant HCV- and liver disease-associated gene set provided a molecular portrait of several processes occurring in the HCV-infected liver. It included (1). genes expressed in activated lymphocytes infiltrating the cirrhotic liver, and activated liver macrophages; (2). genes involved in remodeling of extracellular matrix-cell and cell-cell interactions associated with cytoskeleton rearrangements; (3). genes related to the anti-apoptotic pathway of Bcl-2 signaling; and (4). genes involved with the interferon response and virus-host interactions. In conclusion, our microarray analysis identified several potential gene markers of HCV-associated liver disease and contributed to our rapidly expanding database of experiments describing HCV pathogenesis.