Improved Denaturation of Small RNA Duplexes and Its Application for Northern Blotting.

Small RNAs (sRNAs) are short (18-30 nucleotide) noncoding RNA molecules, which control gene expression and pathogen response in eukaryotes. They are associated with and guide nucleases to target nucleic acids by nucleotide base pairing. We found that current techniques for small RNA detection are adversely affected by the presence of complementary RNA. Thus we established FDF-PAGE (fully denaturing formaldehyde polyacrylamide gel electrophoresis), which dramatically improves denaturation efficiency and subsequently the detection of sequestered sRNAs.

A strand specific real-time RT-PCR method for the targeted detection of the three species (vRNA, cRNA and mRNA) of infectious salmon anaemia virus (ISAV) replicative RNA.

Three species of viral-derived RNA (vRNA, cRNA and mRNA) are produced during an infectious salmon anaemia virus (ISAV) infection. Conventional real-time RT-PCR (RT-qPCR) targeting ISAV segment 8 provides a very sensitive method for the detection of ISAV RNA, however it does not differentiate between these three individual RNA species. In this study, strand-specific tagged primers have been utilised in the RT reaction to specifically produce cDNA corresponding to each of the 3 viral RNA types produced from ISAV segment 8 for the subsequent detection by real-time PCR. The RNA species-specific assay was successfully used to specifically distinguish synthetic T7-produced RNA transcripts representing the 3 species of ISAV RNA at levels up to approximately 10(5)-fold higher than the other types. In addition, the method was applied to investigate the production of segment 8 RNA in time-course tissue culture experiments performed at optimal (15°C), sub-optimal (20°C) and inadequate (25°C) temperatures for replication or in the presence of a chemical inhibitor to vary the RNA populations and investigate its effectiveness. Variation in RNA production was observed between the optimal and sub-optimal temperatures and in the presence of the chemical inhibitor. Production of all RNA species was completely inhibited at 25°C indicating the potential usefulness of the assay as a tool in the understanding of ISAV replication and transcription dynamics.

CFTR regulation of epithelial sodium channel.

Cystic fibrosis (CF) is a lethal genetic disorder, characterized by both clinical and genetic complexities, and arises as a result of mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. The gene encodes a Cl(-) channel belonging to the ABC (ATP Binding Cassette) superfamily of transporters. The members of this superfamily use ATP hydrolysis to fulfill their function as active transporters. So far, CFTR is the only member of this family to function as a cAMP-activated Cl(-) channel. Intense research following the cloning of the CFTR gene has extended the role of the CFTR beyond that of a Cl(-) channel. One of the best recognized, yet still controversial, functions of the CFTR is its ability to modulate the functioning of other transporters. The modulation of epithelial Na(+) channel (ENaC) function serves as a prime example of regulatory function of the CFTR. In this chapter, we will briefly describe an integrated protocol consisting of biochemical and electrophysiological approaches to study the regulation of ENaC by CFTR.

Comprehensive gene expression analysis in human periodontal ligaments of the mandibular third molars performing vertical movement and the maxillary second premolars with occlusal contact.

OBJECTIVES: The periodontal ligament (PDL) is thought to be an important tissue in vertical movement during tooth eruption, but the precise molecular mechanism is not known. Thereto, comprehensive gene expression was analyzed in human PDL of mandibular third molars performing vertical movement and maxillary second premolars with occlusal contact. DESIGN: The expression profile of 9,243 genes in the PDL of one subject was compared between vertically moving third molars and second premolars with occlusal contact by DNA microarray. RESULTS: The expression of 27 genes showed more than a 10-fold difference between third molars and second premolars. The expression of CALB1 (encoding calbindin 1), CYP26A1 (encoding cytochrome P450, family 26, subfamily A, polypeptide 1), SPOCK3 (encoding testican-3), CCK (encoding cholecystokinin) and SCRG1 (encoding scrapie responsive protein 1) was more than 30-fold higher in PDLs of the third molars than the second premolars. CALB1 is reported to increase at the pressure side of PDL during experimental orthodontic tooth movement in rats. Interestingly, in this study, CALB1 expression showed the largest difference. In contrast, CRCT1 (encoding cysteine-rich C-terminal 1), SPRP3 (encoding small proline-rich protein 3), IL8 (encoding interleukin 8) and MMP12 (encoding matrix metalloproteinase 12) showed more than 100-fold higher expression in PDLs of the second premolars than the third molars. CONCLUSION: The present comprehensive gene expression in PDLs provides new insights into the molecular mechanism during the vertical tooth movement.

Homogeneous fluorescence assays for RNA diagnosis by pyrene-conjugated 2'-O-methyloligoribonucleotides.

We developed a bispyrene-conjugated 2'-O-methyloligoribonucleotide as an RNA-specific RNA-probe. The probe hybridized with the complementary RNA, greatly enhancing fluorescence and discriminating RNA from DNA. The assay was carried out in homogeneous aqueous media without removing the unbound probe from the detection solution. This homogeneous fluorescence assay also discriminated mismatch sequences in the target RNA. These pyrene probes could possess high potential to detect RNA in biological specimens simply.

Reagent preparation and storage for amplification of microarray hybridization targets with a fully automated system.

The advent of automated systems for gene expression profiling has accentuated the need for the development of convenient and cost-effective methods for reagent preparation. We have developed a method for the preparation and storage of pre-aliquoted cocktail plates that contain all reagents required for amplification of nucleic acid by reverse transcription and in vitro transcription reactions. Plates can be stored at -80 degrees C for at least 1 month and kept in a hotel at 4 degrees C for at least 24 h prior to use. Microarray data quality generated from these pre-aliquoted reagent plates is not statistically different between cRNA amplified with stored cocktails and cRNA amplified with freshly prepared cocktails. Deployment of pre-aliquoted, stored cocktail plates in a fully automated system not only increases the throughput of amplifying cRNA targets from thousands of RNA samples, but could also considerably reduce reagent costs and potentially improve process robustness.

Signal transducer and activator of transcription-6 (STAT6) is a constitutively expressed survival factor in human prostate cancer.

BACKGROUND: Signal transducer and activator of transcription (STAT)-6 is a member of the STAT family of latent transcription factors. In this investigation, we examined STAT6 expression in clinical prostate cancer tissue specimen and determined its role in prostate cell proliferation and migration. METHODS: STAT6 expression in cell lines and tissues was analyzed by RT-PCR, IHC and/or immunoblot analyses. Down-regulation of STAT6 expression was achieved by STAT6 siRNA and its effect on cell migration and apoptosis was measured. RESULTS: STAT6 is highly expressed in the fibromuscular stroma of prostate cancer specimens. STAT6 is also expressed in the malignant epithelial layer and prostate intraepithelial neoplasia (PIN). STAT6 expression was significantly correlated with high histological grades of prostate cancer as well as with tumor size. Our data indicate deregulated STAT6 mRNA and protein expression in prostate cancer cells with high levels in the non-cancerous HPV 18C-1 and cancerous DU145 cell lines and low levels in PC3 and LNCaP cells. Phosphorylated STAT6 was expressed in all three cancer cell lines DU145, PC3, and LNCaP. Down-regulation of STAT6 using siRNA leads to the induction of early apoptosis in DU145 cells and inhibits migration of these cells. Significant reduction in cell viability and transcriptional down-regulation of the anti-apoptotic protein Bcl-X(L) was observed followed by STAT6 down-regulation in DU145 cells. Interestingly STAT6 also regulates transcription of 15-lipoxygenase-1 gene in DU145 cells. CONCLUSIONS: Our data suggest that STAT6 is a survival factor in prostate cancer and regulates the genetic transcriptional program that is responsible for prostate cancer progression.

Molecular detection of SS18-SSX fusion gene transcripts by cRNA in situ hybridization in synovial sarcoma using formalin-fixed, paraffin-embedded tumor tissue specimens.

SS18-SSX fusion genes resulting from a chromosomal translocation t(X;18)(p11.2;q11.2) are a genetic hallmark of synovial sarcoma. Although such cytogenetic or molecular aberrations have mostly been detected by fluorescence in situ hybridization or reverse transcription-polymerase chain reaction, the expression of SS18-SSX has been poorly investigated at a cellular or tissue level. In this study, biotinylated tyramide (BT)-based in situ hybridization (ISH) was performed to detect SS18-SSX transcripts using formalin-fixed, paraffin-embedded tissues from 15 synovial sarcomas. Digoxigenin-labeled cRNA probes flanking the fusion points of SS18-SSX1 and SS18-SSX2 were generated by in vitro transcription, and hybridized signals were detected by a streptavidin-biotin complex method after chemical enhancement with BT. The localizations of signals were compared with the immunohistochemical expressions of epithelial or neuroectodermal markers and those of cell adhesion including cytokeratins (CAM5.2, AE1/AE3, CK7), epithelial membrane antigen, E-cadherin, beta-catenin, c-erbB-2 (HER2/neu), CD56, and claudin-1. The ISH signals of the SS18-SSX transcripts were identified in 13 synovial sarcomas, and their fusion types correlated with those determined by reverse transcription-polymerase chain reaction. In biphasic tumors, the ISH signals tended to localize to epithelial areas, whereas spindle-cell areas or monophasic fibrous tumors showed a less intense or focal expression pattern. Notably, the expression patterns of AE1/AE3, CK7, and c-erbB-2 often colocalized with the ISH signals (7 of 11 cases positive for each marker). Our results suggest that BT-based ISH can be used as a molecular technique for the detection of SS18-SSX using formalin-fixed, paraffin-embedded tissues.

Using luminescent nanoparticles as staining probes for Affymetrix GeneChips.

Microarray technology provides efficient access to genetic information using miniaturized, high-density arrays of DNA probes. We investigated the application of luminescent nanoparticles as probes for Affymetrix GeneChips detection without the need for signal amplification. Our goal is to investigate the feasibility of using luminescent nanoparticles as probes in a commercial microarray system without changing its configurations. With the present imaging modality and existing optical excitation and detection systems of the Affymetrix GeneChips, our early results indicate that nanoparticles not only can be used for GeneChip labeling but also are superior to the traditional fluorescent protein streptavidin-phycoerythrin (SAPE). The advantage of the particles lies in a simplified staining procedure, higher photobleaching threshold, and enhanced luminescence signal. The nanoparticles can be used for detection of low-abundance targets without any amplification step. A concentration detection limit of 50 fM has been achieved. This work demonstrates the feasibility of using luminescent nanoparticles as probes for commercial microarray systems, making them less costly, more reproducible, and potentially quantitative.

Detection of cRNA hybridized on a DNA chip using a tetrakis-acridinyl peptide cassette, consisting of TAP and partially self-complementary oligonucleotide, d[A18(TA)51].

A tetrakis-acridinyl peptide (TAP) cassette, consisting of a double-stranded region of alternating AT sequence bound to TAP and a single stranded overhanging sequence of continuous dA, was prepared by mixing TAP with d[A18(TA)51]. A TAP cassette could be applied to the fluorometric detection of hybridized DNA on the DNA chip, which was prepared by stamping a 45-meric DNA probe onto a gold-coated plastic chip using a high-precision spotter developed at RIKEN. Spots on the DNA chip were imaged by a CCD camera after hybridization with 65-meric target single-stranded DNAs carrying a continuous dA20 sequence (dA tail) on the DNA chip after treatment with a TAP cassette. Their fluorescence intensity on the DNA chip showed a good linear correlation with the concentration of the target DNAs in the range from 10 pM to 1 nM. Fluorescence of their spots derived from the TAP cassette remaining on the surface of the DNA chip through the dA tail of the hybridized target DNA. Furthermore, the TAP cassette could be successfully applied to the quantitative detection of complementary RNAs (cRNAs) prepared from rat brain with reverse transcription and in vitro transcription.

Microarray analysis of cardiovascular diseases.

Microarray analysis is a powerful technique for high-throughput, global transcriptonomic profiling of gene expression. It holds great promise for analyzing the genetic and molecular bases of cardiovascular diseases and various other complex diseases and permits the analysis of thousands of genes simultaneously, both in diseased and nondiseased tissues and/or cell lines. Microarrays or microchips are made by depositing spots of DNA or oligonucleotides representing thousands of genes on a solid support such as a coated glass surface, and can allow the comparison of gene expression patterns in any two samples. Total RNA is isolated from the tissue or cells of interest, converted to cDNA and then cRNA labeled with biotin, and hybridized to the chips. Hybridization signals are then quantified and compared among different samples. We used oligonucleotide microarrays to obtain an unbiased assessment of expression levels of thousands of genes simultaneously in normal and diseased coronary arteries. Fifty-six genes showed differential expression in atherosclerotic coronary artery tissues, and 49 of them represent new linked genes for coronary artery disease. These studies can generate novel hypotheses relating to the pathologies of disease and further studies with animal models, molecular biology, cell biology, and biochemistry will validate these hypotheses and provide novel insights into the pathogenesis of disease.

Automated target preparation for microarray-based gene expression analysis.

DNA microarrays have rapidly evolved toward a platform for massively paralleled gene expression analysis. Despite its widespread use, the technology has been criticized to be vulnerable to technical variability. Addressing this issue, recent comparative, interplatform, and interlaboratory studies have revealed that, given defined procedures for "wet lab" experiments and data processing, a satisfactory reproducibility and little experimental variability can be achieved. In view of these advances in standardization, the requirement for uniform sample preparation becomes evident, especially if a microarray platform is used as a facility, i.e., by different users working in the laboratory. While one option to reduce technical variability is to dedicate one laboratory technician to all microarray studies, we have decided to automate the entire RNA sample preparation implementing a liquid handling system coupled to a thermocycler and a microtiter plate reader. Indeed, automated RNA sample preparation prior to chip analysis enables (1) the reduction of experimentally caused result variability, (2) the separation of (important) biological variability from (undesired) experimental variation, and (3) interstudy comparison of gene expression results. Our robotic platform can process up to 24 samples in parallel, using an automated sample preparation method that produces high-quality biotin-labeled cRNA ready to be hybridized on Affymetrix GeneChips. The results show that the technical interexperiment variation is less pronounced than with manually prepared samples. Moreover, experiments using the same starting material showed that the automated process yields a good reproducibility between samples.

Gene expression profiling of primary cutaneous melanoma and clinical outcome.

BACKGROUND: Gene expression profiling data for human primary cutaneous melanomas are scarce because of the lack of retrospective collections of frozen tumors. To identify differentially expressed genes that may be involved in melanoma progression and prognosis, we investigated the relationship between gene expression profiles and clinical outcome in a cohort of patients with primary melanoma. METHODS: Labeled complementary RNA (cRNA) from each tissue sample was hybridized to a pangenomic 44K 60-mer oligonucleotide microarray. Class comparison and class prediction analyses were performed to identify genes whose expression in primary melanomas was associated with 4-year distant metastasis-free survival among 58 patients with at least 4 years of follow-up, distant metastasis, or death. Results were validated immunohistochemically at the protein level in 176 independent primary melanomas from patients with a median clinical follow-up of 8.5 years. Survival was analyzed with a Cox multivariable model and stratified log-rank test. All statistical tests were two-sided. RESULTS: We identified 254 genes that were associated with distant metastasis-free survival of patients with primary melanoma. These 254 genes include genes involved in activating DNA replication origins, such as minichromosome maintenance genes and geminin. Twenty-three of these genes were studied at the protein level; expression of five (MCM4, P = .002; MCM3, P = .030; MCM6, P = .004; KPNA2, P = .021; and geminin, P = .004) was statistically significantly associated with overall survival in the validation set. In a multivariable Cox model adjusted for tumor thickness, ulceration, age, and sex, expression of MCM4 (hazard ratio [HR] of death = 4.04, 95% confidence interval [CI] = 1.39 to 11.76; P = .010) and MCM6 (HR of death = 7.42, 95% CI = 1.99 to 27.64; P = .003) proteins was still statistically significantly associated with overall survival. CONCLUSION: We identified 254 genes whose expression was associated with metastatic dissemination of cutaneous melanomas. These genes may shed light on the molecular mechanisms underlying poor prognosis in melanoma patients.

Gene chip analyses reveal differential genetic responses to iron deficiency in rat duodenum and jejunum.

Previous studies revealed novel genetic changes in the duodenal mucosa of iron-deprived rats during postnatal development. These observations are now extended to compare the genetic response to iron deficiency in the duodenum versus jejunum of 12-wk-old rats. cRNA samples were prepared from the duodenal and jejunal mucosa of three groups each of control and iron-deficient rats and hybridized with RAE 230A and 230B gene chips (Affymetrix). Stringent data reduction strategies were employed. Results showed that several genes were similarly induced in both gut segments, including DMT1, Dcytb, transferrin receptor 1, heme oxygenase 1, metallothionein, the Menkes copper ATPase (ATP7A), tripartitie motif protein 27, and the sodium-dependent vitamin C transporter. However, a subset of genes showed regulation in only one or the other gut segment. In duodenum only, gastrokine 1, trefoil factor 1 and claudin 2 were induced by iron-deficiency. Other genes previously identified were only regulated in the duodenum. Overall, these studies demonstrate similarities and distinct differences in the genetic response to iron deprivation in the duodenum versus jejunum and provide evidence that more distal gut segments also may play a role in increasing iron absorption in iron-deficiency anemia.

Gene chip analyses reveal differential genetic responses to iron deficiency in rat duodenum and jejunum

Previous studies revealed novel genetic changes in the duodenal mucosa of iron-deprived rats during post-natal development. These observations are now extended to compare the genetic response to iron deficiency in the duodenum versus jejunum of 12-wk-old rats. cRNA samples were prepared from the duodenal and jejunal mucosa of three groups each of control and iron-deficient rats and hybridized with RAE 230A and 230B gene chips (Affymetrix). Stringent data reduction strategies were employed. Results showed that several genes were similarly induced in both gut segments, including DMT1, Dcytb, transferrin receptor 1, heme oxygenase 1, metallothionein, the Menkes copper ATPase (ATP7A), tripartitie motif protein 27, and the sodium-dependent vitamin C transporter. However, a subset of genes showed regulation in only one or the other gut segment. In duodenum only, gastrokine 1, trefoil factor 1 and claudin 2 were induced by iron-deficiency. Other genes previously identified were only regulated in the duodenum. Overall, these studies demonstrate similarities and distinct differences in the genetic response to iron deprivation in the duodenum versus jejunum and provide evidence that more distal gut segments also may play a role in increasing iron absorption in iron-deficiency anemia.

Gene expression profiling of primary tumor cell populations using laser capture microdissection, RNA transcript amplification, and GeneChip microarrays.

Gene expression profiling from microdissected cell populations is a powerful approach to explore molecular processes involved in development and solid tumor biology. In this chapter, we detail robust and validated methods for tissue preparation and isolation of high-quality RNA from microdissected cell populations. A protocol is also provided for linear transcript amplification using as little as 10 ng of total RNA to produce labeled cRNA targets for hybridization to GeneChip high-density oligonucleotide microarrays. Particular emphasis is placed on troubleshooting each technical step in the protocol and measures of quality assurance for both RNA isolation and resulting microarray data.

Neuropeptide Y Y1 receptor mRNA in rodent brain: distribution and colocalization with melanocortin-4 receptor.

The central neuropeptide Y (NPY) Y1 receptor (Y1-R) system has been implicated in feeding, endocrine, and autonomic regulation. In the present study, we systematically examined the brain distribution of Y1-R mRNA in rodents by using radioisotopic in situ hybridization histochemistry (ISHH) with a novel sensitive cRNA probe. Within the rat hypothalamus, Y1-R-specific hybridization was observed in the anteroventral periventricular, ventromedial preoptic, suprachiasmatic, paraventricular (PVH), dorsomedial, ventromedial, arcuate, and mamillary nuclei. In the rat, Y1-R mRNA expression was also seen in the subfornical organ, anterior hypothalamic area, dorsal hypothalamic area, and in the lateral hypothalamic area. In addition, Y1-R hybridization was evident in several extrahypothalamic forebrain and hindbrain sites involved in feeding and/or autonomic regulation in the rat. A similar distribution pattern of Y1-R mRNA was observed in the mouse brain. Moreover, by using a transgenic mouse line expressing green fluorescent protein under the control of the melanocortin-4 receptor (MC4-R) promoter, we observed Y1-R mRNA expression in MC4-R-positive cells in several brain sites such as the PVH and central nucleus of the amygdala. Additionally, dual-label ISHH demonstrated that hypophysiotropic PVH cells coexpress Y1-R and pro-thyrotropin-releasing hormone mRNAs in the rat. These observations are consistent with the proposed roles of the central NPY/Y1-R system in energy homeostasis.

Direct labeling of RNA with multiple biotins allows sensitive expression profiling of acute leukemia class predictor genes.

Direct labeling of RNA is an expedient method for labeling large quantities (e.g. micrograms) of target RNA for microarray analysis. We have developed an efficient labeling system that uses T4 RNA ligase to attach a 3'-biotinylated donor molecule to target RNA. Microarray analyses indicate that directly labeled RNA is uniformly labeled, has higher signal intensity than comparable labeling methods and achieves high transcript detection sensitivity. The labeled donor molecule we have developed allows the attachment of multiple biotins, which increases target signal intensity up to 30%. We have used this direct-labeling method to detect previously discovered class predictor genes for two types of cancer: acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). In order to test the sensitivity of direct RNA labeling, we analyzed the AML and ALL expression profiles for predictor genes that were previously found to show elevated expression in the disease state. Direct labeling of AML poly(A) RNA detects 90% of the class predictor genes that are detected by the IVT-based target amplification method used to discover the genes. These results indicate that the detection sensitivity, simplicity (single tube reaction) and speed (2 h) of this direct labeling protocol may be ideal for diagnostic applications that do not require target amplification.

Establishment of a real-time RT-PCR for the determination of absolute amounts of IGF-I and IGF-II gene expression in liver and extrahepatic sites of the tilapia.

We developed a one-tube two-temperature real-time RT-PCR that allows to absolutely quantify the gene expression of hormones using the standard curve method. As our research focuses on the expression of the insulin-like growth factors (IGFs) in bony fish, we established the technique for IGF-I and IGF-II using the tilapia (Oreochromis niloticus) as model species. As approach, we used primer extension adding a T7 phage polymerase promoter (21 nt) to the 5' end of the antisense primers. This procedure avoids the disadvantages arising from plasmids. Total RNA extracted from liver was subjected to conventional RT-PCR to create templates for in vitro transcription of IGF-I and IGF-II cRNA. Correct template sizes including the T7 promoter were verified (IGF-I: 91 nt; IGF-II: 94 nt). The PCR products were used to create IGF-I and IGF-II cRNAs which were quantified in dot blot by comparison with defined amounts of standardised kanamycin mRNA. Standardised threshold cycle (Ct) values for IGF-I and IGF-II mRNA were achieved by real-time RT-PCR and used to create standard curves. To allow sample normalisation the standard curve was also established for beta-actin as internal calibrator (template: 86 nt), and validation experiments were performed demonstrating similar amplification efficiencies for target and reference genes. Based on the standard curves, the absolute amounts of IGF-I and IGF-II mRNA were determined for liver (IGF-I: 8.90+/-1.90 pg/microg total RNA, IGF-II: 3.59+/-0.98 pg/microg total RNA) and extrahepatic sites, such as heart, kidney, intestine, spleen, gills, gonad, and brain considering the different lengths of cRNAs and mRNAs by correction factors. The reliability of the method was confirmed in additional experiments. The amplification of descending dilutions of cRNA and total liver RNA resulted in parallel slopes of the amplification curves. Furthermore, amplification plots of the standard cRNA and the IGF-I and IGF-II mRNAs showed signals starting at the expected Ct values. Thus, the one-tube RT-PCR described here is highly sensitive (detection level approximately 2 pg/microg total RNA) and allows precise absolute quantification. The method is rapid as there are neither separate reverse transcriptions nor post-amplification steps, and can be executed with low risk of contamination. Therefore, it will be helpful when investigating gene expression in any species and tissue whenever absolute levels are of concern.

Graphical technique for identifying a monotonic variance stabilizing transformation for absolute gene intensity signals.

BACKGROUND: The usefulness of log2 transformation for cDNA microarray data has led to its widespread application to Affymetrix data. For Affymetrix data, where absolute intensities are indicative of number of transcripts, there is a systematic relationship between variance and magnitude of measurements. Application of the log2 transformation expands the scale of genes with low intensities while compressing the scale of genes with higher intensities thus reversing the mean by variance relationship. The usefulness of these transformations needs to be examined. RESULTS: Using an Affymetrix GeneChip dataset, problems associated with applying the log2 transformation to absolute intensity data are demonstrated. Use of the spread-versus-level plot to identify an appropriate variance stabilizing transformation is presented. For the data presented, the spread-versus-level plot identified a power transformation that successfully stabilized the variance of probe set summaries. CONCLUSION: The spread-versus-level plot is helpful to identify transformations for variance stabilization. This is robust against outliers and avoids assumption of models and maximizations.