Identification of suitable endogenous control genes for microRNA expression profiling of childhood medulloblastoma and human neural stem cells

Analysis of circulating microRNA during early gestation in Japanese black cattle

Circulating microRNAs (miRNAs) have been used as biomarkers for various diseases and physiological conditions in humans and mice; studies in domestic animals, particularly cattle, are limited. The importance of early pregnancy diagnosis (especially within the 21-d cow estrous cycle) in the livestock industry is extremely high. This study compared the circulating miRNAs in bred non-pregnant and pregnant Japanese Black cows, explored miRNAs as biomarkers for early pregnancy diagnosis, and established a measurement system that included selecting an appropriate reference miRNA and determining the effect of hemolysis on miRNA quantification in plasma. miRNA was extracted from the plasma of Japanese Black cows on day 21 after artificial insemination and subjected to a customized bovine oligonucleotide microarray for expression analysis. Differentially expressed miRNAs and reference miRNA candidates were selected and validated using reverse transcription-quantitative PCR (RT-qPCR). An appropriate endogenous reference miRNA for normalization was selected using NormFinder software. To evaluate the effect of hemolysis on miRNA quantification, hemolyzed samples were prepared using plasma from four cows in the estrous cycle and subjected to RT-qPCR. A total of 124 miRNAs were detected in bovine plasma by microarray analysis in bred non-pregnant and pregnant cows. The levels of five circulating miRNAs were significantly higher in pregnant cows than in bred non-pregnant cows, and 24 miRNAs were detected only in the pregnant group. NormFinder analysis and RT-qPCR validation showed that miR-2455 was an appropriate reference miRNA in the plasma of bred non-pregnant and pregnant Japanese Black cows, and miR-19b, miR-25, miR-29a, and miR-148a were significantly higher in the pregnant group. These four circulating miRNAs did not change during the estrous cycle and were less affected by hemolysis. In the current study, we found four miRNAs, miR-19b, miR-25, miR-29a, and miR-148a, which were present at high levels in the plasma of pregnant Japanese Black cows. Since these miRNAs are less affected by hemolysis, they may potentially be used as biomarkers for early pregnancy diagnosis in cattle.

Differences in toxicity, mitochondrial function and miRNome in human cells exposed in vitro to Cd as CdS quantum dots or ionic Cd

Cadmium is toxic to humans, although Cd-based quantum dots exerts less toxicity. Human hepatocellular carcinoma cells (HepG2) and macrophages (THP-1) were exposed to ionic Cd, Cd(II), and cadmium sulfide quantum dots (CdS QDs), and cell viability, cell integrity, Cd accumulation, mitochondrial function and miRNome profile were evaluated. Cell-type and Cd form-specific responses were found: CdS QDs affected cell viability more in HepG2 than in THP-1; respective IC₂₀ values were ∼3 and ∼50 μg ml^-1. In both cell types, Cd(II) exerted greater effects on viability. Mitochondrial membrane function in HepG2 cells was reduced 70 % with 40 μg ml^-1 CdS QDs but was totally inhibited by Cd(II) at corresponding amounts. In THP-1 cells, CdS QDs has less effect on mitochondrial function; 50 μg ml^-1 CdS QDs or equivalent Cd(II) caused 30 % reduction or total inhibition, respectively. The different in vitro effects of CdS QDs were unrelated to Cd uptake, which was greater in THP-1 cells. For both cell types, changes in the expression of miRNAs (miR-222, miR-181a, miR-142-3p, miR-15) were found with CdS QDs, which may be used as biomarkers of hazard nanomaterial exposure. The cell-specific miRNome profiles were indicative of a more conservative autophagic response in THP-1 and as apoptosis as in HepG2.

The Long-Term Effects of Early Life Stress on the Modulation of miR-19 Levels

MicroRNAs (miRNAs), one of the major small non-coding RNA classes, have been proposed as regulatory molecules in neurodevelopment and stress response. Although alterations in miRNAs profiles have been implicated in several psychiatric and neurodevelopmental disorders, the contribution of individual miRNAs in brain development and function is still unknown. Recent studies have identified miR-19 as a key regulator of brain trajectories, since it drives the differentiation of neural stem cells into mature neurons. However, no findings are available on how vulnerability factors for these disorders, such as early life stress (ELS), can modulate the expression of miR-19 and its target genes. To reach our aim, we investigated miR-19 modulation in human hippocampal progenitor stem cells (HPCs) treated with cortisol during 3 days of proliferation and harvested immediately after the end of the treatment or after 20 days of differentiation into mature neurons. We also analyzed the long-term expression changes of miR-19 and of its validated target genes, involved in neurodevelopment and inflammation, in the hippocampus of adult rats exposed or not to prenatal stress (PNS). Interestingly, we observed a significant downregulation of miR-19 levels both in proliferating (FC = −1.59, p-value = 0.022 for miR-19a; FC = −1.79, p-value = 0.016 for miR-19b) as well as differentiated HPCs (FC = −1.28, p-value = 0.065 for miR-19a; FC = −1.75, p-value = 0.047 for miR-19b) treated with cortisol. Similarly, we found a long-term decrease of miR-19 levels in the hippocampus of adult PNS rats (FC = −1.35, p-value = 0.025 for miR-19a; FC = −1.43, p-value = 0.032 for miR-19b). Among all the validated target genes, we observed a significant increase of NRCAM (FC = 1.20, p-value = 0.027), IL4R (FC = 1.26, p-value = 0.046), and RAPGEF2 (FC = 1.23, p-value = 0.020).We suggest that ELS can cause a long-term downregulation of miR-19 levels, which may be responsible of alterations in neurodevelopmental pathways and in immune/inflammatory processes, leading to an enhanced risk for mental disorders later in life. Intervention strategies targeting miR-19 may prevent alterations in these pathways, reducing the ELS-related effects.

Glycine-Histidine-Lysine (GHK) Alleviates Neuronal Apoptosis Due to Intracerebral Hemorrhage via the miR-339-5p/VEGFA Pathway

Glycine-histidine-lysine (GHK) is a human tripeptide that enhances wound healing, exerts neuroprotective effects against neurodegenerative disease, and improves tissue regeneration. This study examined whether GHK can alleviate injury due to intracerebral hemorrhage (ICH). Briefly, adult Wistar rats in GHK pretreatment groups were injected with GHK (1 or 10 mg/kg, i.p.) every 24 h for 3 days. Water content and intact neurons were detected in the rats 3 days after ICH, and the neurological deficit scores were examined in the rats at 4, 24, 72, and 168 h after ICH. Apoptosis was evaluated via caspase-3 immunohistochemistry, Nissl staining, and TUNEL assay. We also examined the effect of GHK on the expression of related proteins in SH-SY5Y cells via Western blotting. The expression of miR-339-5p was examined via real-time polymerase chain reaction analyses. GHK improved neurological deficits, reduced water content in the brain and inhibited neuronal apoptosis in ICH rats. It also prevented the apoptosis of SH-SY5Y cells with hemin treatment. Furthermore, GHK downregulated miR-339-5p expression, and overexpression of miR-339-5p partially reversed the anti-apoptotic effects of GHK in SH-SY5Y cells. Our findings suggest that the p38 MAPK pathway is involved in the GHK-induced downregulation of miR-339-5p, and that the miR-339-5p/VEGFA axis plays a role in preventing neuronal apoptosis following ICH injury. These findings indicate that GHK may represent a novel therapeutic strategy for ICH.

Impact of miRNA-mRNA Profiling and Their Correlation on Medulloblastoma Tumorigenesis

Medulloblastoma (MB) is a clinically challenging, childhood brain tumor with a diverse genetic makeup and differential miRNA profile. Aiming to identify deregulated miRNAs in MB, the miRNA expression profile of human MB samples was compared to that of normal cerebellar tissues. As a result, 8 upregulated and 64 downregulated miRNAs were identified in MB samples. Although various algorithms have been developed to predict the interaction between miRNA-mRNA pairs, the complexity and fidelity of miRNA-mRNA remain a concern. Therefore, to identify the signatures of miRNA-mRNA interactions essential for MB pathogenesis, miRNA profiling, RNA sequencing, and ingenuity pathway analysis (IPA) were performed in the same primary human MB samples. Further, when miR-217 was inhibited, a significant upregulation of predicted target genes SIRT1, ROBO1, FOXO3, and SMAD7 in HDMB03 cells was observed, confirming the validity of our approach. Functional analysis revealed that the inhibition of miR-217 in HDMB03 cells suppresses colony formation, migration, invasion, promoted apoptosis, and arrested cell population in S phase, indicating that manipulation of miR-217 may have a therapeutic potential for MB patients. Therefore, our study provides an essential platform for future investigations of specific miRNAs responsible for MB pathogenesis.

Droplet Digital PCR versus qPCR for gene expression analysis with low abundant targets: from variable nonsense to publication quality data

Quantitative PCR (qPCR) has become the gold standard technique to measure cDNA and gDNA levels but the resulting data can be highly variable, artifactual and non-reproducible without appropriate verification and validation of both samples and primers. The root cause of poor quality data is typically associated with inadequate dilution of residual protein and chemical contaminants that variably inhibit Taq polymerase and primer annealing. The most susceptible, frustrating and often most interesting samples are those containing low abundant targets with small expression differences of 2-fold or lower. Here, Droplet Digital PCR (ddPCR) and qPCR platforms were directly compared for gene expression analysis using low amounts of purified, synthetic DNA in well characterized samples under identical reaction conditions. We conclude that for sample/target combinations with low levels of nucleic acids (Cq ≥ 29) and/or variable amounts of chemical and protein contaminants, ddPCR technology will produce more precise, reproducible and statistically significant results required for publication quality data. A stepwise methodology is also described to choose between these complimentary technologies to obtain the best results for any experiment.

A Proteogenomic Approach to Understanding MYC Function in Metastatic Medulloblastoma Tumors

Brain tumors are the leading cause of cancer-related deaths in children, and medulloblastoma is the most prevalent malignant childhood/pediatric brain tumor. Providing effective treatment for these cancers, with minimal damage to the still-developing brain, remains one of the greatest challenges faced by clinicians. Understanding the diverse events driving tumor formation, maintenance, progression, and recurrence is necessary for identifying novel targeted therapeutics and improving survival of patients with this disease. Genomic copy number alteration data, together with clinical studies, identifies c-MYC amplification as an important risk factor associated with the most aggressive forms of medulloblastoma with marked metastatic potential. Yet despite this, very little is known regarding the impact of such genomic abnormalities upon the functional biology of the tumor cell. We discuss here how recent advances in quantitative proteomic techniques are now providing new insights into the functional biology of these aggressive tumors, as illustrated by the use of proteomics to bridge the gap between the genotype and phenotype in the case of c-MYC-amplified/associated medulloblastoma. These integrated proteogenomic approaches now provide a new platform for understanding cancer biology by providing a functional context to frame genomic abnormalities.

Comparison of different normalization strategies for the analysis of glomerular microRNAs in IgA nephropathy

Small nucleolar RNAs (snoRNAs) have been used for normalization in glomerular microRNA (miRNA) quantification without confirmation of validity. Our aim was to identify glomerular reference miRNAs in IgA nephropathy. We compared miRNAs in human paraffin-embedded renal biopsies from patients with cellular-crescentic IgA-GN (n = 5; crescentic IgA-GN) and non-crescentic IgA-GN (n = 5; IgA-GN) to mild interstitial nephritis without glomerular abnormalities (controls, n = 5). Laser-microdissected glomeruli were used for expression profiling of 762 miRNAs by low-density TaqMan arrays (cards A and B). The comparison of different normalization methods (GeNormPlus, NormFinder, global mean and snoRNAs) in crescentic IgA-GN, IgA-GN and controls yielded similar results. However, levels of significance and the range of relative expression differed. In median, two normalization methods demonstrated similar results. GeNormPlus and NormFinder gave different top ranked reference miRNAs. Stability ranking for snoRNAs varied between cards A and B. In conclusion, we suggest the geometric mean of the most stable reference miRNAs found in GeNormPlus (miR-26b-5p), NormFinder (miR-28-5p) and snoRNAs (RNU44) as reference. It should be considered that significant differences could be missed using one particular normalization method. As a starting point for glomerular miRNA studies in IgA nephropathy we provide a library of miRNAs.

Analysis of Deregulated microRNAs and Their Target Genes in Gastric Cancer

BACKGROUND

MicroRNAs (miRNAs) are widely studied non-coding RNAs that modulate gene expression. MiRNAs are deregulated in different tumors including gastric cancer (GC) and have potential diagnostic and prognostic implications. The aim of our study was to determine miRNA profile in GC tissues, followed by evaluation of deregulated miRNAs in plasma of GC patients. Using available databases and bioinformatics methods we also aimed to evaluate potential target genes of confirmed differentially expressed miRNA and validate these findings in GC tissues.

METHODS

The study included 51 GC patients and 51 controls. Initially, we screened miRNA expression profile in 13 tissue samples of GC and 12 normal gastric tissues with TaqMan low density array (TLDA). In the second stage, differentially expressed miRNAs were validated in a replication cohort using qRT-PCR in tissue and plasma samples. Subsequently, we analyzed potential target genes of deregulated miRNAs using bioinformatics approach, determined their expression in GC tissues and performed correlation analysis with targeting miRNAs.

RESULTS

Profiling with TLDA revealed 15 deregulated miRNAs in GC tissues compared to normal gastric mucosa. Replication analysis confirmed that miR-148a-3p, miR-204-5p, miR-223-3p and miR-375 were consistently deregulated in GC tissues. Analysis of GC patients' plasma samples showed significant down-regulation of miR-148a-3p, miR-375 and up-regulation of miR-223-3p compared to healthy subjects. Further, using bioinformatic tools we identified targets of replicated miRNAs and performed disease-associated gene enrichment analysis. Ultimately, we evaluated potential target gene BCL2 and DNMT3B expression by qRT-PCR in GC tissue, which correlated with targeting miRNA expression.

CONCLUSIONS

Our study revealed miRNA profile in GC tissues and showed that miR-148a-3p, miR-223-3p and miR-375 are deregulated in GC plasma samples, but these circulating miRNAs showed relatively weak diagnostic performance as sole biomarkers. Target gene analysis demonstrated that BCL2 and DNMT3B expression in GC tissue correlated with their targeting miRNA expression.

Stability Assessment of Candidate Reference Genes in Urine Sediment of Prostate Cancer Patients for miRNA Applications

We aimed at assessing the stability of candidate reference genes in urine sediments of men subjected to digital rectal examination for suspected prostate cancer (PCa). Two microRNAs (miR-191 and miR-25) and 1 small nucleolar RNA (SNORD48) were assayed in 35 post-DRE urine sediments of men with PCa and in 26 subjects with histologically confirmed benign prostatic hyperplasia (BPH). The stability of candidate reference genes was assessed through BestKeeper algorithm and equivalence test. miR-200b and miR-452 were used to test for the effect of normalization on target genes. Our results proved miR-191 to be the most stable gene, showing the lowest degree of variation and the highest stability value. miR-25 and SNORD48 values fell beyond the cutoff of acceptability. In conclusion, we recommend the use of miR-191 for normalization purposes in post-DRE urine sediments.

Validation of miRNA genes suitable as reference genes in qPCR analyses of miRNA gene expression in Atlantic salmon (Salmo salar)

Background

MicroRNAs (miRNAs) are an abundant class of endogenous small RNA molecules that downregulate gene expression at the post-transcriptional level. They play important roles by regulating genes that control multiple biological processes, and recent years there has been an increased interest in studying miRNA genes and miRNA gene expression. The most common method applied to study gene expression of single genes is quantitative PCR (qPCR). However, before expression of mature miRNAs can be studied robust qPCR methods (miRNA-qPCR) must be developed. This includes identification and validation of suitable reference genes. We are particularly interested in Atlantic salmon (Salmo salar). This is an economically important aquaculture species, but no reference genes dedicated for use in miRNA-qPCR methods has been validated for this species. Our aim was, therefore, to identify suitable reference genes for miRNA-qPCR methods in Salmo salar.

Results

We used a systematic approach where we utilized similar studies in other species, some biological criteria, results from deep sequencing of small RNAs and, finally, experimental validation of candidate reference genes by qPCR to identify the most suitable reference genes. Ssa-miR-25-3p was identified as most suitable single reference gene. The best combinations of two reference genes were ssa-miR-25-3p and ssa-miR-455-5p. These two genes were constitutively and stably expressed across many different tissues. Furthermore, infectious salmon anaemia did not seem to affect their expression levels. These genes were amplified with high specificity, good efficiency and the qPCR assays showed a good linearity when applying a simple cybergreen miRNA-PCR method using miRNA gene specific forward primers.

Conclusions

We have identified suitable reference genes for miRNA-qPCR in Atlantic salmon. These results will greatly facilitate further studies on miRNA genes in this species. The reference genes identified are conserved genes that are identical in their mature sequence in many aquaculture species. Therefore, they may also be suitable as reference genes in other teleosts. Finally, the systematic approach used in our study successfully identified suitable reference genes, suggesting that this may be a useful strategy to apply in similar validation studies in other aquaculture species.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-945) contains supplementary material, which is available to authorized users.

Identification of circulating microRNAs for the differential diagnosis of Parkinson's disease and Multiple System Atrophy

Background: Parkinson's disease (PD) is a progressive neurodegenerative disorder which may be misdiagnosed with atypical conditions such as Multiple System Atrophy (MSA), due to overlapping clinical features. MicroRNAs (miRNAs) are small non-coding RNAs with a key role in post-transcriptional gene regulation. We hypothesized that identification of a distinct set of circulating miRNAs (cmiRNAs) could distinguish patients affected by PD from MSA and healthy individuals. Results. Using TaqMan Low Density Array technology, we analyzed 754 miRNAs and found 9 cmiRNAs differentially expressed in PD and MSA patients compared to healthy controls. We also validated a set of 4 differentially expressed cmiRNAs in PD and MSA patients vs. controls. More specifically, miR-339-5p was downregulated, whereas miR-223^*, miR-324-3p, and mir-24 were upregulated in both diseases. We found cmiRNAs specifically deregulated in PD (downregulation of miR-30c and miR-148b) and in MSA (upregulation of miR-148b). Finally, comparing MSA and PD, we identified 3 upregulated cmiRNAs in MSA serum (miR-24, miR-34b, miR-148b). Conclusions. Our results suggest that cmiRNA signatures discriminate PD from MSA patients and healthy controls and may be considered specific, non-invasive biomarkers for differential diagnosis.

Identification and validation of miRNAs as endogenous controls for RQ-PCR in blood specimens for breast cancer studies

Introduction

A prerequisite to accurate interpretation of RQ-PCR data is robust data normalization. A commonly used method is to compare the cycle threshold (C_T) of target miRNAs with those of a stably expressed endogenous (EC) miRNA(s) from the same sample. Despite the large number of studies reporting miRNA expression patterns, comparatively few appropriate ECs have been reported thus far. The purpose of this study was to identify stably expressed miRNAs with which to normalize RQ-PCR data derived from human blood specimens.

Methods

MiRNA profiling of approximately 380 miRNAs was performed on RNA derived from blood specimens from 10 women with breast cancer and 10 matched controls. Analysis of mean expression values across the dataset (GME) identified stably expressed candidates. Additional candidates were selected from the literature and analyzed by the geNorm algorithm. Further validation of three candidate ECs by RQ-PCR was performed in a larger cohort (n = 40 cancer, n = 20 control) was performed, including analysis by geNorm and NormFinder algorithms.

Results

Microarray screening identified 10 candidate ECs with expression patterns closest to the global mean. Geometric averaging of candidate ECs from the literature using geNorm identified miR-425 as the most stably expressed miRNA. MiR-425 and miR-16 were the best combination, achieving the lowest V-value of 0.185. Further validation by RQ-PCR confirmed that miR-16 and miR-425 were the most stably expressed ECs overall. Their combined use to normalize expression data enabled the detection of altered target miRNA expression that reliably differentiated between cancers and controls in human blood specimens.

Conclusion

This study identified that the combined use of 2 miRNAs, (miR-16 and miR-425) to normalize RQ-PCR data generated more reliable results than using either miRNA alone, or use of U6. Further investigation into suitable ECs for use in miRNA RQ-PCR studies is warranted.

Identification and validation of reference genes for qPCR detection of serum microRNAs in colorectal adenocarcinoma patients

Serum microRNAs (miRNAs) have become a highlighted research hotspot, especially for their great potential as a novel promising non-invasive biomarker in cancer diagnosis. The most frequently used approach for serum miRNAs detection is quantitative real time polymerase chain reaction (qPCR). In order to obtain reliable qPCR data of miRNAs expression, the use of reference genes as endogenous control is undoubtly necessary. However, no systematic evaluation and validation of reference genes for normalizing qPCR analysis of serum miRNAs has been reported in colorectal adenocarcinoma. We firstly profiled pooled serum of colorectal adenocarcinoma, colorectal adenoma and healthy controls and selected a list of 13 miRNAs as candidate reference genes. U6 snRNA (U6) and above-mentioned 13 miRNAs were included in further confirmation by qPCR. As a result, 5 miRNAs (miR-151a-3p, miR-4446-3p, miR-221-3p, miR-93-5p and miR-3184-3p) were not detected in all samples and 2 miRNAs (miR-197-3p and miR-26a-5p) were relatively low with median Cq more than 35, and were excluded from further stability analysis. Then variable stability of other 6 miRNAs (miR-103b, miR-484, miR-16-5p, miR-3615, miR-18a-3p and miR-191-5p) and U6 were evaluated using two algorithms: geNorm and NormFinder which both identified miR-191-5p as the most stably expressed reference gene and selected miR-191-5p and U6 as the most stable pair of reference genes. After validating in an independent large cohorts and selecting miR-92a-3p as target miRNA to evaluate the effect of reference gene, we propose that combination of miR-191-5p and U6 could be used as reference genes for serum microRNAs qPCR data in colorectal adenocarcinoma, colorectal adenoma and healthy controls.

Similarity on neural stem cells and brain tumor stem cells in transgenic brain tumor mouse models

Although it is believed that glioma is derived from brain tumor stem cells, the source and molecular signal pathways of these cells are still unclear. In this study, we used stable doxycycline-inducible transgenic mouse brain tumor models (c-myc⁺/SV40Tag⁺/Tet-on⁺) to explore the malignant trans-formation potential of neural stem cells by observing the differences of neural stem cells and brain tumor stem cells in the tumor models. Results showed that chromosome instability occurred in brain tumor stem cells. The numbers of cytolysosomes and autophagosomes in brain tumor stem cells and induced neural stem cells were lower and the proliferative activity was obviously stronger than that in normal neural stem cells. Normal neural stem cells could differentiate into glial fibrillary acidic protein-positive and microtubule associated protein-2-positive cells, which were also negative for nestin. However, glial fibrillary acidic protein/nestin, microtubule associated protein-2/nestin, and glial fibrillary acidic protein/microtubule associated protein-2 double-positive cells were found in induced neural stem cells and brain tumor stem cells. Results indicate that induced neural stem cells are similar to brain tumor stem cells, and are possibly the source of brain tumor stem cells.

Identification of optimal reference genes for RT-qPCR in the rat hypothalamus and intestine for the study of obesity

BACKGROUND

Obesity has a complicated metabolic pathology, and defining the underlying mechanisms of obesity requires integrative studies with molecular end points. Real-time quantitative PCR (RT-qPCR) is a powerful tool that has been widely utilized. However, the importance of using carefully validated reference genes in RT-qPCR seems to have been overlooked in obesity-related research. The objective of this study was to select a set of reference genes with stable expressions to be used for RT-qPCR normalization in rats under fasted vs re-fed and chow vs high-fat diet (HFD) conditions.

DESIGN

Male long-Evans rats were treated under four conditions: chow/fasted, chow/re-fed, HFD/fasted and HFD/re-fed. Expression stabilities of 13 candidate reference genes were evaluated in the rat hypothalamus, duodenum, jejunum and ileum using the ReFinder software program. The optimal number of reference genes needed for RT-qPCR analyses was determined using geNorm.

RESULTS

Using geNorm analysis, we found that it was sufficient to use the two most stably expressed genes as references in RT-qPCR analyses for each tissue under specific experimental conditions. B2M and RPLP0 in the hypothalamus, RPS18 and HMBS in the duodenum, RPLP2 and RPLP0 in the jejunum and RPS18 and YWHAZ in the ileum were the most suitable pairs for a normalization study when the four aforementioned experimental conditions were considered.

CONCLUSIONS

Our study demonstrates that gene expression levels of reference genes commonly used in obesity-related studies, such as ACTB or RPS18, are altered by changes in acute or chronic energy status. These findings underline the importance of using reference genes that are stable in expression across experimental conditions when studying the rat hypothalamus and intestine, because these tissues have an integral role in the regulation of energy homeostasis. It is our hope that this study will raise awareness among obesity researchers on the essential need for reference gene validation in gene expression studies.