Complementary Techniques

Identification of Reference Gene for Quantitative Gene Expression in Early-Term and Late-Term Cultured Canine Fibroblasts Derived from Ear Skin

Fibroblasts are cells that reside within the fibrous or loose connective tissues of most mammalian organs. For research purposes, fibroblasts are often subjected to long-term culture under defined conditions, during which their properties can significantly change. It is essential to understand and document these changes to obtain reliable outcomes. For the quantification of specific gene expressions, the most reliable and widely used technique is quantitative real-time polymerase chain reaction (qRT-PCR). Here, we assessed the impact of a reference gene's stability on a qRT-PCR analysis of long-term cultured canine skin fibroblasts. After successfully isolating the fibroblasts from canine skin tissues, they were cultured and evaluated for proliferation and β-galactosidase activity at different passage numbers. With extended culture, the fibroblasts showed a long doubling time and elevated β-galactosidase activity. Using three widely used algorithms, geNorm, Normfinder, and Bestkeeper, we identified HPRT1, YWHAZ, and GUSB as the most stable reference genes for both early- and late-passage fibroblasts. Conventional reference genes such as GAPDH were found to be less stable than those genes. The normalization of Vimentin by the stable genes showed statistical differences, whereas normalization by an unstable gene did not. Collectively, this study indicates that using stable reference genes is essential for accurately and reliably measuring gene expression in both early- and late-passage fibroblasts. These findings provide valuable insights into internal controls for gene expression studies and are expected to be utilized for analyzing gene expression patterns in molecular biology research.

Assessing equivalent and inverse change in genes between diverse experiments

Background: It is important to identify when two exposures impact a molecular marker (e.g., a gene's expression) in similar ways, for example, to learn that a new drug has a similar effect to an existing drug. Currently, statistically robust approaches for making comparisons of equivalence of effect sizes obtained from two independently run treatment vs. control comparisons have not been developed. Results: Here, we propose two approaches for evaluating the question of equivalence between effect sizes of two independent studies: a bootstrap test of the Equivalent Change Index (ECI), which we previously developed, and performing Two One-Sided t-Tests (TOST) on the difference in log-fold changes directly. The ECI of a gene is computed by taking the ratio of the effect size estimates obtained from the two different studies, weighted by the maximum of the two p-values and giving it a sign indicating if the effects are in the same or opposite directions, whereas TOST is a test of whether the difference in log-fold changes lies outside a region of equivalence. We used a series of simulation studies to compare the two tests on the basis of sensitivity, specificity, balanced accuracy, and F1-score. We found that TOST is not efficient for identifying equivalently changed gene expression values (F1-score = 0) because it is too conservative, while the ECI bootstrap test shows good performance (F1-score = 0.95). Furthermore, applying the ECI bootstrap test and TOST to publicly available microarray expression data from pancreatic cancer showed that, while TOST was not able to identify any equivalently or inversely changed genes, the ECI bootstrap test identified genes associated with pancreatic cancer. Additionally, when investigating publicly available RNAseq data of smoking vs. vaping, no equivalently changed genes were identified by TOST, but ECI bootstrap test identified genes associated with smoking. Conclusion: A bootstrap test of the ECI is a promising new statistical approach for determining if two diverse studies show similarity in the differential expression of genes and can help to identify genes which are similarly influenced by a specific treatment or exposure. The R package for the ECI bootstrap test is available at https://github.com/Hecate08/ECIbootstrap.

Screening Potential Reference Genes in Tuta absoluta with Real-Time Quantitative PCR Analysis under Different Experimental Conditions

Tuta absoluta is one of the most significant invasive pests affecting tomato plants worldwide. RT-qPCR has emerged as one of the most sensitive and accurate methods for detecting gene expression data. The screening of stable internal reference genes is the most critical step for studying the molecular mechanisms of environmental adaptability. The stable reference genes expressed in T. absoluta under specific experimental conditions have not yet been clarified. In this study, seven candidate reference genes (RPL27, RPS13, RPS15, EF1-α, TUB, TBP, and β-actin) and their optimal numbers were evaluated under biotic (developmental stages and adult tissues) and abiotic (insecticide, temperature, and plant VOC) conditions using four software programs. Our results identified the following reference genes and numbers as optimal: three genes (EF1-α, RPS13, and RPL27) for different developmental stages (egg, larva, pupa, unmated adult), two genes (RPS13 and TBP) for adult tissues (antenna, head, thorax, abdomen, leg), two genes (TBP and RPS13) for insecticides (Bacillus thuringiensis, chlorpyrifos, abamectin-aminomethyl, and chlorantraniliprole), two genes (RPL27 and TUB) for temperature-induced stresses (0, 25, and 40 °C), and two genes (RPS13 and TUB) for VOC-induced stresses (nonanal, α-phellandrene, and tomato leaves). Our results provide a reference for selecting appropriate reference genes for further study of the functional genes of T. absoluta under different experimental conditions.

Identification and Evaluation of Suitable Reference Genes for RT-qPCR Analysis in Hippodamia variegata (Coleoptera: Coccinellidae) Under Different Biotic and Abiotic Conditions

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is an accurate and convenient technique for quantifying expression levels of the target genes. Selection of the appropriate reference gene is of the vital importance for RT-qPCR analysis. Hippodamia variegata is one of the most important predatory natural enemies of aphids. Recently, transcriptome and genome sequencings of H. variegata facilitate the gene functional studies. However, there has been rare investigation on the detection of stably expressed reference genes in H. variegata. In the current study, by using five analytical tools (Delta Ct, geNorm, NormFinder, BestKeeper, and RefFinder), eight candidate reference genes, namely, Actin, EF1α, RPL7, RPL18, RPS23, Tubulin-α, Tubulin-β, and TufA, were evaluated under four experimental conditions including developmental stages, tissues, temperatures, and diets. As a result, a specific set of reference genes were recommended for each experimental condition. These findings will help to improve the accuracy and reliability of RT-qPCR data, and lay a foundation for further exploration on the gene function of H. variegata.

Comprehensive Assessment of Candidate Reference Genes for Gene Expression Studies Using RT-qPCR in Tamarixia radiata, a Predominant Parasitoid of Diaphorina citri

Tamarixia radiata (Waterston) is a predominant parasitoid of the Asian citrus psyllid (ACP), a destructive citrus pest and vector of huanglongbing (HLB) disease in the fields of southern China. To explore the functioning of target genes in T. radiata, the screening of specific reference genes is critical for carrying out the reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) under different experimental conditions. However, no reference gene(s) for T. radiata has yet been reported. Here, we selected seven housekeeping genes of T. radiate and evaluated their stability under the six conditions (developmental stage, sex, tissue, population, temperature, diet) by using RefFinder software, which contains four different programs (geNorm, ΔCt, BestKeeper, and NormFinder). Pairwise variation was analyzed by geNorm software to determine the optimal number of reference genes during the RT-qPCR analysis. The results reveal better reference genes for differing research foci: 18S and EF1A for the developmental stage; PRS18 and EF1A for sex, PRS18 and RPL13 for different tissues (head, thorax, abdomen); EF1A and ArgK between two populations; β-tubulin and EF1A for different temperatures (5, 15, 25, 35 °C); and ArgK and PRS18 for different feeding diets. Furthermore, when the two optimal and two most inappropriate reference genes were chosen in different temperatures and tissue treatments, respectively, the corresponding expression patterns of HSP70 (as the reporter gene) differed substantially. Our study provides, for the first time, a more comprehensive list of optimal reference genes from T. radiata for use in RT-qPCR analysis, which should prove beneficial for subsequent functional investigations of target gene(s) in this natural enemy of ACP.

Screening and evaluation of the stability of expression of reference genes in Lymantria dispar (Lepidoptera: Erebidae) using qRT-PCR

The quantitative real-time polymerase chain reaction (qRT-PCR) has rapidly become the most sensitive and accurate method for the quantitative analysis of gene expression. Normalization of gene expression to that of relatively stably expressed housekeeping genes is required to facilitate the study of gene expression and to obtain more accurate RT-PCR data. However, no studies of the stability of expression of housekeeping genes in Lymantria dispar have been reported. In the present study, BestKeeper, GeNorm and NormFinder statistical software was used to evaluate the expression of thirteen candidate reference genes in L. dispar under different conditions. The expression levels of candidate reference genes were determined for two biological factors (developmental stages and tissues) and four abiotic treatments (temperature, insecticide, CO₂ and starvation). The results showed that the best candidate reference genes in L. dispar were TUB, AK, RPS15 for developmental stages, RPL32 and GAPDH for tissues, ACTB and EF1-α for CO₂ stress, GAPDH and RPL32 for temperature stress, RPS3 and GAPDH for insecticide stress, and GAPDH and RPS3 for starvation stress. In summary, EF1-α and TUB are preferential housekeeping genes in L. dispar under various conditions. These results provide a basis for the further study of functional genes of L. dispar.

PPIA, HPRT1, and YWHAZ Genes Are Suitable for Normalization of mRNA Expression in Long-Term Expanded Human Mesenchymal Stem Cells

Long-term expansion of mesenchymal stem cells (MSCs) under defined culture conditions is necessary in human stem cell therapy. However, it alters the characteristics of MSCs. Since quantitative real time polymerase chain reaction (qRT-PCR) is widely used as one of the key analytical methods for comparative characterization, the validation of reference genes (RGs) for normalization under each experimental condition is important to achieve reliable qRT-PCR results. Therefore, the most stable RGs for long-term expanded bone marrow- and umbilical cord blood-derived MSCs (BM-MSCs and UCB-MSCs) under defined culture conditions for up to 20 passages were evaluated. The more apparent alterations in characteristics such as differentiation capacity, proliferation, senescence, and the expression of RGs were noted in BM-MSCs than UCB-MSCs during long-term expansion. The RG validation programs (GeNorm and NormFinder) suggested that PPIA, HPRT1, and YWHAZ were suitable for normalization in qRT-PCR regardless of MSC types and long-term culture expansion, and the traditional RGs (ACTB and GAPDH) were less stable in long-term expanded MSCs. In addition, the use of these RGs for normalization of OCT4 expression in long-term expanded BM-MSCs showed that a less stable RG (GAPDH) showed contrasting data compared to other RGs. Therefore, the use of RGs such as PPIA, HPRT1, and YWHAZ for normalization in qRT-PCR experiments is highly recommended for long-term expanded MSCs to generate accurate and reliable data.

Evaluation of the expression stability of β-actin under bacterial infection in Macrobrachium nipponense

The selection of a suitable reference gene is an important prerequisite for the precise analysis of target gene expression by real-time quantitative PCR (qPCR). The present study aims to explore the expression pattern of the Macrobrachium nipponense (M. nipponense) β-actin gene under Aeromonas hydrophila bacterial infection conditions. The complete sequence of the β-actin gene from M. nipponense was cloned by PCR. Identified and named β-actin genes were searched in the NCBI database, and the characteristics of the β-actin gene were analyzed using bioinformatics methods. The expression profiles of β-actin under stresses challenged by bacteria after 3, 6, 12, 24 and 48 h were investigated by measuring Ct values by qPCR. The prokaryotic expression vector pET-30a-actin was constructed by PCR and recombinant DNA techniques. Fused protein was induced by IPTG in the transformed Escherichia coli BL21 (DE3). Recombinant rActin was purified by nickel column. The bioinformatics analysis result revealed that the deduced protein encoded by the β-actin gene from M. nipponense had the highest homology with other prawns in the homologous assay (99%). The phylogenetic tree indicates that the β-actin from M. nipponense and other crustaceans have a single cluster. The qPCR results revealed that a stable expression of β-actin was observed in response to the A. hydrophila challenge for 3-48 h, and the Ct value was 22 ± 1.5. β-actin was ranked as a stable gene after the bacterial challenge, which was selected as the appropriate reference gene in M. nipponense.

Quantification of gene expression while taking into account RNA alternative splicing

Gene expression has been widely used in functional genomics research; however, the gene expressions quantified with different methods have been frequently inconsistent, thus challenging the conclusions from such research. Here we have addressed this issue, while taking into account RNA alternative splicing. We found that when a gene was subjected to RNA alternative splicing, it was impossible or difficult to properly quantify the expression of a transcript of the gene or its overall expression using quantitative real-time PCR (qPCR), Northern hybridization, microarray, or serial analysis of gene expression. Shot-gun RNA-seq was the most proper to quantify the expression of a transcript or a gene in such cases. Moreover, the expressions of individual transcripts quantified by shot-gun RNA-seq were highly reproducible (r = 0.90-0.98) between individuals. Therefore, shot-gun or full-length RNA-seq should be the method of choice to properly quantify the expression of a transcript or a gene.

Ectopic expression of GhCOBL9A, a cotton glycosyl-phosphatidyl inositol-anchored protein encoding gene, promotes cell elongation, thickening and increased plant biomass in transgenic Arabidopsis

Cellulose is a major component of plant cell walls and is necessary for plant morphogenesis and biomass. COBL (COBRA-Like) proteins have been shown to be key regulators in the orientation of cell expansion and cellulose crystallinity status. To clarify the role of a cotton COBL gene, GhCOBL9A, we conducted the ectopic expression and functional analysis in Arabidopsis. Previous study showed that GhCOBL9A was preferentially expressed during secondary cell wall biosynthesis in cotton fibers, and showed a significant co-expression pattern with cellulose synthase genes. Here, we detected that overexpression of GhCOBL9A induced the up-regulation of genes related to cellulose synthesis and enhanced the cellulose deposition. As a result, GhCOBL9A transgenic plants displayed increased hypocotyl and root lengths in early development, and cell wall thickening at the SCW stage. Notably, overexpression of GhCOBL9A led to an erect, robust-stature phenotype and brought higher biomass in mature plants. In addition, overexpression of GhCOBL9A in Arabidopsis AtCOBL4 mutants, a paralogous gene of GhCOBL9A, also led to a stronger growth potential, but the Atcobl4 mutant phenotype could not be rescued, implying the functional divergence of GhCOBL9A and AtCOBL4 paralogs. Taken together, these results suggest that overexpression of GhCOBL9A contributes to plant cell elongation and thickening, and increased biomass, which provides references for further utilizing GhCOBL9A to improve yield and quality traits in cotton and other species.

Genome-Wide Association Studies Reveal Genetic Variation and Candidate Genes of Drought Stress Related Traits in Cotton (Gossypium hirsutum L.)

Cotton is an important industrial crop worldwide and upland cotton (Gossypium hirsutum L.) is most widely cultivated in the world. Due to ever-increasing water deficit, drought stress brings a major threat to cotton production. Thus, it is important to reveal the genetic basis under drought stress and develop drought tolerant cotton cultivars. To address this issue, in present study, 319 upland cotton accessions were genotyped by 55,060 single nucleotide polymorphisms (SNPs) from high-density CottonSNP80K array and phenotyped nine drought tolerance related traits. The two datasets were used to identify quantitative trait nucleotides (QTNs) for the above nine traits using multi-locus random-SNP-effect mixed linear model method. As a result, a total of 20 QTNs distributed on 16 chromosomes were found to be significantly associated with six drought tolerance related traits. Of the 1,326 genes around the 20 QTNs, 205 were induced after drought stress treatment, and 46 were further mapped to Gene ontology (GO) term "response to stress." Taken genome-wide association study (GWAS) analysis, RNA-seq data and qRT-PCR verification, four genes, RD2 encoding a response to desiccation 2 protein, HAT22 encoding a homeobox-leucine zipper protein, PIP2 encoding a plasma membrane intrinsic protein 2, and PP2C encoding a protein phosphatase 2C, were proposed to be potentially important for drought tolerance in cotton. These results will deepen our understanding of the genetic basis of drought stress tolerance in cotton and provide candidate markers to accelerate the development of drought-tolerant cotton cultivars.

Identification and validation of reference genes for gene expression analysis in Aphidius gifuensis (Hymenoptera: Aphidiidae)

Reference genes have been utilized in estimating gene expression levels using quantitative reverse transcriptase-quantitative polymerase chain reaction (qRT-PCR) analysis. Aphidius gifuensis Ashmaed is one of the most widely used biological control agents for aphids. The biological properties of this species have been studied in detail, and current investigations are focused on elucidating the regulatory mechanisms in its host However, the appropriate reference genes for target gene expression studies have not been identified. In this study, the expression profiles of 12 candidate reference genes were evaluated under different experimental conditions(development stage, sex, tissue type, diet) by using dedicated algorithms, including geNorm, Normfinder, BestKeeper, and ΔCt. In addition, RefFinder was used to rank the overall stability of the candidate genes. Finally, we recommend three optimal reference genes for the normalization of qRT-PCR data in the presence of specific variables, which include ACTB, RPL13, and PPI for different developmental stages; RPS18, ACTB, and RPL13 for sexes; RPL13, PRII3, and RPS18 in different tissue types; and RPL13, RPL27, and ACTB in diverse diets. The present study has identified optimal reference genes that could be used in estimating the expression levels of specific genes under these conditions following the Minimum Information for publication of Quantitative real-time PCR Experiments (MIQE) guidelines, which would facilitate in advancements in functional genomics research on A. gifuensis.

Comprehensive evaluation of candidate reference genes for gene expression studies in Lysiphlebia japonica (Hymenoptera: Aphidiidae) using RT-qPCR

Lysiphlebia japonica (Ashmead) is a predominant parasitoid of cotton-melon aphids in the fields of northern China with a proven ability to effectively control cotton aphid populations in early summer. For accurate normalization of gene expression in L. japonica using quantitative reverse transcriptase-polymerase chain reaction (RT-qPCR), reference genes with stable gene expression patterns are essential. However, no appropriate reference genes is L. japonica have been investigated to date. In the present study, 12 selected housekeeping genes from L. japonica were cloned. We evaluated the stability of these genes under various experimental treatments by RT-qPCR using four independent (geNorm, NormFinder, BestKeeper and Delta Ct) and one comparative (RefFinder) algorithm. We identified genes showing the most stable levels of expression: DIMT, 18S rRNA, and RPL13 during different stages; AK, RPL13, and TBP among sexes; EF1A, PPI, and RPL27 in different tissues, and EF1A, RPL13, and PPI in adults fed on different diets. Moreover, the expression profile of a target gene (odorant receptor 1, OR1) studied during the developmental stages confirms the reliability of the chosen selected reference genes. This study provides for the first time a comprehensive list of suitable reference genes for gene expression studies in L. japonica and will benefit subsequent genomics and functional genomics research on this natural enemy.

Stable reference genes in granulosa cells of bovine dominant follicles during follicular growth, FSH stimulation and maternal aging

The aim of the present study was to determine a set of reference genes in granulosa cells of dominant follicles that are suitable for relative gene expression analyses during maternal and follicular aging. Granulosa cells of growing and preovulatory dominant follicles were collected from aged and young cows (maternal aging study) and from FSH-stimulated follicles developing under different durations of FSH treatment (follicular aging study). The mRNA levels of the two commonly used reference genes (GAPDH, ACTB) and four novel genes (UBE2D2, EIF2B2, SF3A1, RNF20) were analysed using cycle threshold values. Results revealed that mRNA levels of GAPDH, ACTB, EIF2B2, RNF20, SF3A1 and UBE2D2 were similar (P>0.05) between dominant follicle type, age and among follicles obtained after FSH-stimulation, but differed (P=0.005) due to mRNA processing (i.e. with versus without amplification). The stability of reference genes was analysed using GeNorm, DeltaCT and NormFinder programs and comprehensive ranking order was determined using RefFinder. The mRNA levels of GAPDH and ACTB were less stable than those of UBE2D2 and EIF2B2. The geometric mean of multiple genes (UBE2D2, EIF2B2, GAPDH and SF3A1) is a more appropriate reference control than the use of a single reference gene to compare relative gene expression among dominant and FSH-stimulated follicles during maternal and/or follicular aging studies.

Granulosa cell biomarkers to predict pregnancy in ART: pieces to solve the puzzle

Cumulus and mural granulosa cells of the ovarian follicle surround and interact with the developing oocyte. These follicular cells reflect the oocyte's overall health and may indicate subsequent developmental competence of embryos. Biomarkers of granulosa cells associated with individual oocytes could potentially be used in assisted reproduction to indicate which embryos have the best chance of implanting in the uterus and completing gestation. In this review, we have performed a comprehensive assessment of the recent literature for human cumulus and mural granulosa cell mRNA biomarkers as they relate to pregnancy and live birth. A critical discussion of variables affecting granulosa gene expression profiles for in vitro fertilization patients, including patient demographics and ovarian stimulation regimens, is presented. Although studies with microarray data were evaluated, this synopsis focuses on expressed genes that have been validated by quantitative RT-PCR. Furthermore, we summarize the current published data that support or refute identified granulosa expressed genes as potential biomarkers of embryos that give rise to ongoing pregnancy and live birth. Finally, we review studies that offer predictive models for embryo selection for uterine transfer based on biomarkers that show differential gene expression.

Tracking vaginal, anal and oral infection in a mouse papillomavirus infection model

Noninvasive and practical techniques to longitudinally track viral infection are sought after in clinical practice. We report a proof-of-principle study to monitor the viral DNA copy number using a newly established mouse papillomavirus (MmuPV1) mucosal infection model. We hypothesized that viral presence could be identified and quantified by collecting lavage samples from cervicovaginal, anal and oral sites. Nude mice infected at these sites with infectious MmuPV1 were tracked for up to 23 weeks starting at 6 weeks post-infection. Viral DNA copy number was determined by SYBR Green Q-PCR analysis. In addition, we tracked viral DNA load through three complete oestrous cycles to pinpoint whether there was a correlation between the DNA load and the four stages of the oestrous cycle. Our results showed that high viral DNA copy number was reproducibly detected from both anal and cervicovaginal lavage samples. The infection and disease progression were further confirmed by histology, cytology, in situ hybridization, immunohistochemistry and transmission electron microscopy. Interestingly, the viral copy number fluctuated over the oestrous cycle, with the highest level at the oestrus stage, implying that multiple sampling might be necessary to provide a reliable diagnosis. Virus DNA was detected in oral lavage samples at a later time after infection. Lower viral DNA load was found in oral samples when compared with those in anal and vaginal tracts. To our knowledge, our study is the first in vivo study to sequentially monitor papillomavirus infection from mucosal anal, oral and vaginal tracts in a preclinical model.

Human amniotic fluid stem cells labeled with up-conversion nanoparticles for imaging-monitored repairing of acute lung injury

Human amniotic fluid stem (hAFS) cells have generated a great deal of excitement in cell-based therapies and regenerative medicine. Here, we examined the effect of hAFS cells labeled with dual-polymer-coated UCNP-PEG-PEI nanoparticles in a murine model of acute lung injury (ALI). We observed hAFS cells migration to the lung using highly sensitive in vivo upconversion luminescence (UCL) imaging. We demonstrated that hAFS cells remained viable and retained their ability to differentiate even after UCNP-PEG-PEI labeling. More importantly, hAFS cells displayed remarkable positive effects on ALI-damaged lung tissue repair compared with mouse bone marrow mesenchymal stem cells (mBMSCs), which include recovery of the integrity of alveolar-capillary membrane, attenuation of transepithelial leukocyte and neutrophil migration, and down-regulation of proinflammatory cytokine and chemokine expression. Our work highlights a promising role for imaging-guided hAFS cell-based therapy in ALI.

Gene expression profiling of selenophosphate synthetase 2 knockdown in Drosophila melanogaster

Selenium (Se) is an important trace element for many organisms and is incorporated into selenoproteins as selenocysteine (Sec). In eukaryotes, selenophosphate synthetase SPS2 is essential for Sec biosynthesis. In recent years, genetic disruptions of both Sec biosynthesis genes and selenoprotein genes have been investigated in different animal models, which provide important clues for understanding the Se metabolism and function in these organisms. However, a systematic study on the knockdown of SPS2 has not been performed in vivo. Herein, we conducted microarray experiments to study the transcriptome of fruit flies with knockdown of SPS2 in larval and adult stages. Several hundred differentially expressed genes were identified in each stage. In spite that the expression levels of other Sec biosynthesis genes and selenoprotein genes were not significantly changed, it is possible that selenoprotein translation might be reduced without impacting the mRNA level. Functional enrichment and network-based analyses revealed that although different sets of differentially expressed genes were obtained in each stage, they were both significantly enriched in the carbohydrate metabolism and redox processes. Furthermore, protein-protein interaction (PPI)-based network clustering analysis implied that several hub genes detected in the top modules, such as Nimrod C1 and regucalcin, could be considered as key regulators that are responsible for the complex responses caused by SPS2 knockdown. Overall, our data provide new insights into the relationship between Se utilization and several fundamental cellular processes as well as diseases.

Genome-wide analysis of CrRLK1L gene family in Gossypium and identification of candidate CrRLK1L genes related to fiber development

Members of the CrRLK1L family, a subgroup of the receptor-like kinase (RLK) gene family, are thought to act as sensors for the integrity of the cell wall and regulators of polar elongation. To better understand the various functions in fiber development, we conducted genome-wide identification and characterization analyses of CrRLK1L family in cotton. Here 44, 40, and 79 CrRLK1L genes were identified from three cotton species: diploid G. raimondii (D5), diploid G. arboreum (A2), and tetraploid G. hirsutum TM-1 (AD1), respectively. The 44 CrRLK1Ls in G. raimondii were anchored to the 12 chromosomes unevenly and were classified into six groups (I-VI), with group II and group IV being further divided into two subgroups (groups IIa and IIb, and IVa and IVb, respectively). These CrRLK1Ls displayed a highly regular pattern of developmental and spatial regulation in cotton. Using the transcriptome data of five chromosomal segment introgression lines (CSILs) and the physical integration of CrRLK1Ls with the quantitative trait loci (QTLs) related to fiber quality traits, we revealed that six CrRLK1L genes were highly associated with fiber development. This study brings new insights into the integrated genome-wide identification of CrRLK1Ls in cotton and provides references for the genetic improvement of cotton fiber.

Development and utilization of a custom PCR array workflow: analysis of gene expression in mycoplasma genitalium and guinea pig (Cavia porcellus)

Transcriptome analysis is a powerful tool for evaluating molecular pathways central to maturation of specific biological processes and disease states. Recently, PCR-based arrays have supplemented microarray and RNA-seq methodologies for studying changes in gene expression levels. PCR arrays are a more cost efficient alternative, however commercially available assemblies are generally limited to only a few more widely researched species (e.g., rat, human, and mouse). Consequently, the investigation of emerging or under-studied species is hindered until such assays are created. To address this need, we present data documenting the success of a developed workflow with enhanced potential to create and validate novel RT-PCR arrays for underrepresented species with whole or partial genome annotation. Utilizing this enhanced workflow, we have achieved a success rate of 80 % for first-round designs for over 400 primer pairs. Of these, ~160 distinct targets were sequence confirmed. Proof of concept studies using two unique arrays, one targeting the pathogenic bacterium Mycoplasma genitalium and the other specific for the guinea pig (Cavia porcellus), allowed us to identify significant (P < 0.05) changes in mRNA expression validated by subsequent qPCR. This flexible and adaptable platform provides a valuable and cost-effective alternative for gene expression analysis.

TqPCR: A Touchdown qPCR Assay with Significantly Improved Detection Sensitivity and Amplification Efficiency of SYBR Green qPCR

The advent of fluorescence-based quantitative real-time PCR (qPCR) has revolutionized the quantification of gene expression analysis in many fields, including life sciences, agriculture, forensic science, molecular diagnostics, and medicine. While SYBR Green-based qPCR is the most commonly-used platform due to its inexpensive nature and robust chemistry, quantifying the expression of genes with low abundance or RNA samples extracted from highly restricted or limited sources can be challenging because the detection sensitivity of SYBR Green-based qPCR is limited. Here, we develop a novel and effective touchdown qPCR (TqPCR) protocol by incorporating a 4-cycle touchdown stage prior to the quantification amplification stage. Using the same cDNA templates, we find that TqPCR can reduce the average Cq values for Gapdh, Rps13, and Hprt1 reference genes by 4.45, 5.47, and 4.94 cycles, respectively, when compared with conventional qPCR; the overall average Cq value reduction for the three reference genes together is 4.95. We further find that TqPCR can improve PCR amplification efficiency and thus increase detection sensitivity. When the quantification of Wnt3A-induced target gene expression in mesenchymal stem cells is analyzed, we find that, while both conventional qPCR and TqPCR can detect the up-regulation of the relatively abundant target Axin2, only TqPCR can detect the up-regulation of the lowly-expressed targets Oct4 and Gbx2. Finally, we demonstrate that the MRQ2 and MRQ3 primer pairs derived from mouse reference gene Tbp can be used to validate the RNA/cDNA integrity of qPCR samples. Taken together, our results strongly suggest that TqPCR may increase detection sensitivity and PCR amplification efficiency. Overall, TqPCR should be advantageous over conventional qPCR in expression quantification, especially when the transcripts of interest are lowly expressed, and/or the availability of total RNA is highly restricted or limited.

Glycosylation-related genes in NS0 cells are insensitive to moderately elevated ammonium concentrations

NS0 and Chinese hamster ovary (CHO) cell lines are used to produce recombinant proteins for human therapeutics; however, ammonium accumulation can negatively impact cell growth, recombinant protein production, and protein glycosylation. To improve product quality and decrease costs, the relationship between ammonium and protein glycosylation needs to be elucidated. While ammonium has been shown to adversely affect glycosylation-related gene expression in CHO cells, NS0 studies have not been performed. Therefore, this study sought to determine if glycosylation in NS0 cells were ammonium-sensitive at the gene expression level. Using a DNA microarray that contained mouse glycosylation-related and housekeeping genes, these genes were analyzed in response to various culture conditions - elevated ammonium, elevated salt, and elevated ammonium with proline. Surprisingly, no significant differences in gene expression levels were observed between the control and these conditions. Further, the elevated ammonium cultures were analyzed using real-time quantitative reverse transcriptase PCR (qRT-PCR) for key glycosylation genes, and the qRT-PCR results corroborated the DNA microarray results, demonstrating that NS0 cells are ammonium-insensitive at the gene expression level. Since NS0 are known to have elevated nucleotide sugar pools under ammonium stress, and none of the genes directly responsible for these metabolic pools were changed, consequently cellular control at the translational or substrate-level must be responsible for the universally observed decreased glycosylation quality under elevated ammonium.

Reference genes for expression studies in hypoxia and hyperglycemia models in human umbilical vein endothelial cells

Human umbilical vein endothelial cell (HUVEC)-based gene expression studies performed under hypoxia and/or hyperglycemia show huge potential for modeling endothelial cell response in cardiovascular disease and diabetes. However, such studies require reference genes that are stable across the whole range of experimental conditions. These reference genes have not been comprehensively defined to date. We applied human genome-wide microarrays and quantitative real-time PCR (qRT-PCR) on RNA obtained from primary HUVEC cultures that were incubated for 24 hr either in euglycemic or in hyperglycemic conditions and then subjected to short-term CoCl₂-induced hypoxia for 1, 3, or 12 hr. Using whole-transcript arrays, we selected 10 commonly used reference genes with no significant expression variation across eight different conditions. These genes were ranked using NormFinder software according to their stability values. Consequently, five genes were selected for validation by qRT-PCR. These were ribosomal protein large P0 (RPLP0), transferrin receptor (TFRC), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), β-glucuronidase (GUSB), and β-actin (ACTB). All five genes displayed stable expression under hyperglycemia. However, only RPLP0 and TFRC genes were stable under hypoxia up to 12 hr. Under hyperglycemia combined with hypoxia up to 12 hr, the expression of RPLP0, TFRC, GUSB, and ACTB genes remained unchanged. Our findings strongly confirm that RPLP0 and TFRC are the most suitable reference genes for HUVEC gene expression experiments subjected to hypoxia and/or hyperglycemia for the given experimental conditions. We provide further evidence that even commonly known references genes require experimental validation for all conditions involved.

SecDF as part of the Sec-translocase facilitates efficient secretion of Bacillus cereus toxins and cell wall-associated proteins

The aim of this study was to explore the role of SecDF in protein secretion in Bacillus cereus ATCC 14579 by in-depth characterization of a markerless secDF knock out mutant. Deletion of secDF resulted in pleiotropic effects characterized by a moderately slower growth rate, aberrant cell morphology, enhanced susceptibility to xenobiotics, reduced virulence and motility. Most toxins, including food poisoning-associated enterotoxins Nhe, Hbl, and cytotoxin K, as well as phospholipase C were less abundant in the secretome of the ΔsecDF mutant as determined by label-free mass spectrometry. Global transcriptome studies revealed profound transcriptional changes upon deletion of secDF indicating cell envelope stress. Interestingly, the addition of glucose enhanced the described phenotypes. This study shows that SecDF is an important part of the Sec-translocase mediating efficient secretion of virulence factors in the Gram-positive opportunistic pathogen B. cereus, and further supports the notion that B. cereus enterotoxins are secreted by the Sec-system.

Analysis of microRNAs in patients with systemic lupus erythematosus, using Solexa deep sequencing

OBJECTIVES

Our aim in this study was to identify and examine the differential expression of microRNAs in patients with systemic lupus erythematosus (SLE).

METHODS

We employed high-quality, high-throughput Solexa sequencing to clone and identify microRNAs in SLE patients and a control group.

RESULTS

From the sequencing data, we identified numerous microRNAs displaying significantly different levels of expression in patients with SLE and in healthy controls. The 212 and 199 microRNAs were upregulated and downregulated, respectively. Only 61 novel microRNAs exhibited significantly different levels of expression in the two groups. The target genes of the novel microRNAs identified in the SLE group were found to have cell metabolism functions. We also analyzed the chromosomal locations of the microRNAs with high level of expression between the two groups. A profile comparison revealed that the majority of transcripts were expressed at a similar level. The functional classes of the most abundant microRNAs were equally represented on each chromosome.

CONCLUSION

We identified novel and known microRNAs significantly enhancing our understanding of the microRNA expression profiles of SLE patients. These data also provide insight into the function of microRNAs in SLE and provide new strategies for future therapies.

CXCL5 as a potential novel prognostic factor in early stage non-small cell lung cancer: results of a study of expression levels of 23 genes

As the current staging system is imprecise for estimating prognosis of early stage non-small cell lung cancer (NSCLC), it is important to identify other methods for selecting high-risk patients after failed surgical treatment. The aim of the study was to evaluate the expression of 23 genes as putative prognostic markers in early stage NSCLC. The study was performed on 109 pairs of tumor and matched unaffected lung tissue surgical specimens taken from stage I and II NSCLC patients. We evaluated the mRNA level of 23 genes using the real-time PCR method. The difference in the expression between the tumor and normal tissue for each gene was analyzed using a general linear model. The influence of gene expression on survival was analyzed by using the proportional hazards model. Eighteen out of the 23 genes showed statistically significant differences in expression between the tumor and non-tumor tissue. For 12 genes (ITGB1, ITGB3, CXCL1, CXCL8, CXCL9, CXCL10, CXCL11, CXCR3, CXCR4, TNF, CHKA, AGFG1, and CTC1), the expression was lower, and for six genes (ITGA5, IL8, IL6, CXCL2, CXCL3, and CXCL12), it was higher in the tumor tissue as compared to the matched normal tissue. Expression changes were more pronounced in squamous cell carcinomas than in adenocarcinomas or large cell carcinomas. Of all the analyzed genes, only CXCL5 was found to statistically significantly (p = 0.04) influence both overall and disease-free survival. Among the 23 genes previously suggested to be relevant for early staged NSCLC patients’ postoperative outcome, only CXCL5 showed a statistically significant prognostic effect.

Selection and evaluation of reference genes for expression analysis using qRT-PCR in the beet armyworm Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae)

Quantitative real-time PCR (qRT-PCR) is a reliable and reproducible technique for measuring and evaluating changes in gene expression. The most common method for analyzing qRT-PCR data is to normalize mRNA levels of target genes to internal reference genes. Evaluating and selecting stable reference genes on a case-by-case basis is critical. The present study aimed to facilitate gene expression studies by identifying the most suitable reference genes for normalization of mRNA expression in qRT-PCR analysis of the beet armyworm Spodoptera exigua (Lepidoptera: Noctuidae). For this purpose, three software tools (geNorm, NormFinder and BestKeeper) were used to investigate 10 candidate reference genes in nine developmental stages and five different tissues (epidermis, head, midgut, fat body and hemolymph) in three larval physiological stages (molting, feeding and wandering stages) of, S. exigua. With the exception of 18S ribosomal RNA (18S), all other candidate genes evaluated, β-actin1(ACT1), β-actin2 (ACT2), elongation factor1(EF1), elongation factor 2 (EF2), Glyceralde hyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L10 (L10), ribosomal protein L17A (L17A), superoxide dismutase (SOD), α-tubulin (TUB),proved to be acceptable reference genes. However, their suitability partly differed between physiological stages and different tissues. L10, EF2 and L17A ranked highest in all tissue sample sets. SOD, ACT2, GAPDH, EF1 and ACT1 were stably expressed in all developmental stage sample sets; ACT2, ACT1 and L10 for larvae sample sets; GAPDH, ACT1 and ACT2 for pupae and adults; SOD and L17A for males; and EF2 and SOD for females. The expression stability of genes varied in different conditions. The findings provided here demonstrated, with a few exceptions, the suitability of most of the 10 reference genes tested in tissues and life developmental stages. Overall, this study emphasizes the importance of validating reference genes for qRT-PCR analysis in S. exigua.

Keratinocyte growth factor gene delivery via mesenchymal stem cells protects against lipopolysaccharide-induced acute lung injury in mice

Acute lung injury (ALI) and acute respiratory distress syndrome (ARDS) are associated with high morbidity and mortality, and have no specific therapy. Keratinocyte growth factor (KGF) is a critical factor for pulmonary epithelial repair and acts via the stimulation of epithelial cell proliferation. Mesenchymal stem cells (MSCs) have been proved as good therapeutic vectors. Thus, we hypothesized that MSC-based KGF gene therapy would have beneficial effects on lipopolysaccharide(LPS)-induced lung injury. After two hours of intratracheal LPS administration to induce lung injury, mice received saline, MSCs alone, empty vector-engineered MSCs (MSCs-vec) or KGF-engineered MSCs (MSCs-kgf) via the tail vein. The MSCs-kgf could be detected in the recipient lungs and the level of KGF expression significantly increased in the MSCs-kgf mice. The MSC-mediated administration of KGF not only improved pulmonary microvascular permeability but also mediated a down-regulation of proinflammatory responses (reducing IL-1β and TNF-α) and an up-regulation of anti-inflammatory responses (increasing cytokine IL-10). Furthermore, the total severity scores of lung injury were significantly reduced in the MSCs-kgf group compared with the other three groups. The underlying mechanism of the protective effect of KGF on ALI may be attributed to the promotion of type II lung epithelial cell proliferation and the enhancement of surfactant synthesis. These findings suggest that MSCs-based KGF gene therapy may be a promising strategy for ALI treatment.

Compound screening and transcriptional profiling in human primary keratinocytes: a brief guideline

Cultured human primary keratinocytes constitute suitable targets for in-depth evaluation of the proliferative or differentiative potential of compounds. There is, however, a double-edged and intrinsically inseparable transition from biological activity to cytotoxicity for any agent under investigation. For that reason, we here first of all present an established protocol for the isolation, cultivation, and analysis of primary foreskin-derived keratinocytes. Taking calcitriol as example, we then reveal how a straightforward photometric cell culture assay can be exploited to assess overall cell viability in response to increasing compound doses. With predetermined cellular cytotoxicity at hand, physiologically meaningful (sub-toxic) compound concentrations for subsequent stimulation of cells can be readily selected, and, in doing so, differentially expressed genes with biological significance can be reliably identified.

Modeling real-time PCR kinetics: Richards reparametrized equation for quantitative estimation of European hake (Merluccius merluccius)

Real-time PCR is the most sensitive method for detection and precise quantification of specific DNA sequences, but it is not usually applied as a quantitative method in seafood. In general, benchmark techniques, mainly cycle threshold (Ct), are the routine method for quantitative estimations, but they are not the most precise approaches for a standard assay. In the present work, amplification data from European hake (Merluccius merluccius) DNA samples were accurately modeled by three sigmoid reparametrized equations, where the lag phase parameter (λc) from the Richards equation with four parameters was demonstrated to be the perfect substitute for Ct for PCR quantification. The concentrations of primers and probes were subsequently optimized by means of that selected kinetic parameter. Finally, the linear correlation among DNA concentration and λc was also confirmed.

Deciphering the calcitriol-induced transcriptomic response in keratinocytes: presentation of novel target genes

Numerous studies to date have been aimed at unraveling the large suite of calcitriol (1α,25-dihydroxyvitamin D(3)) response genes in diverse tissues including skin, where this hormone is involved in regulating keratinocyte proliferation, differentiation, permeability barrier formation, innate immunity promotion, antimicrobial peptide production, and wound healing. However, the various approaches differ considerably in probed cell types, scale, throughput, and statistical reliability and do, of note, not reveal much overlap. To further expand our knowledge on presently elusive targets and characterize the extent of fragmentation of existing datasets, we have performed whole-transcriptome microarray examinations of calcitriol-treated human primary keratinocytes. Out of 28,869 genes investigated, we uncovered 86 differentially expressed (67 upregulated and 19 downregulated) candidates that were functionally clustered into five annotation categories: response to wounding, protease inhibition, secondary metabolite biosynthesis, cellular migration, and amine biosynthetic processes. A complementary RTq-PCR study of 78 nominees selected thereof demonstrated significant differential expression of 55 genes (48 upregulated and seven downregulated) within biological replicates. Our hit list contains nine previously authenticated targets (16.36%, proof of concept) and 46 novel genes (83.6%) that have not yet been explicitly described as being differentially regulated within human primary keratinocytes. Direct vitamin D receptor response element predictions within the regulatory promoter regions of 50 of the RTq-PCR-validated targets agreed with known biological functionality and corroborated our stringent data validation pipeline. Altogether, our results indicate the value of continuing these kinds of gene expression studies, which contribute to an enhanced comprehension of calcitriol-mediated processes that may be dysregulated in human skin pathophysiology.

The miR-183/Taok1 target pair is implicated in cochlear responses to acoustic trauma

Acoustic trauma, one of the leading causes of sensorineural hearing loss, induces sensory hair cell damage in the cochlea. Identifying the molecular mechanisms involved in regulating sensory hair cell death is critical towards developing effective treatments for preventing hair cell damage. Recently, microRNAs (miRNAs) have been shown to participate in the regulatory mechanisms of inner ear development and homeostasis. However, their involvement in cochlear sensory cell degeneration following acoustic trauma is unknown. Here, we profiled the expression pattern of miRNAs in the cochlear sensory epithelium, defined miRNA responses to acoustic overstimulation, and explored potential mRNA targets of miRNAs that may be responsible for the stress responses of the cochlea. Expression analysis of miRNAs in the cochlear sensory epithelium revealed constitutive expression of 176 miRNAs, many of which have not been previously reported in cochlear tissue. Exposure to intense noise caused significant threshold shift and apoptotic activity in the cochleae. Gene expression analysis of noise-traumatized cochleae revealed time-dependent transcriptional changes in the expression of miRNAs. Target prediction analysis revealed potential target genes of the significantly downregulated miRNAs, many of which had cell death- and apoptosis-related functions. Verification of the predicted targets revealed a significant upregulation of Taok1, a target of miRNA-183. Moreover, inhibition of miR-183 with morpholino antisense oligos in cochlear organotypic cultures revealed a negative correlation between the expression levels of miR-183 and Taok1, suggesting the presence of a miR-183/Taok1 target pair. Together, miRNA profiling as well as the target analysis and validation suggest the involvement of miRNAs in the regulation of the degenerative process of the cochlea following acoustic overstimulation. The miR-183/Taok1 target pair is likely to play a role in this regulatory process.

What is beyond a qRT-PCR study on mesenchymal stem cell differentiation properties: how to choose the most reliable housekeeping genes

In the last years, mesenchymal stem cells (MSCs) have been identified as an attractive cell population in regenerative medicine. In view of future therapeutic applications, the study of specific differentiation-related gene expression is a pivotal prerequisite to define the most appropriate MSC source for clinical translation. In this context, it is crucial to use stable housekeeping genes (HGs) for normalization of qRT-PCR to obtain validated and comparable results. By our knowledge, an exhaustive validation study of HGs comparing MSCs from different sources under various differentiation conditions is still missing. In this pivotal study, we compared the expression levels of 12 genes (ACTB, Β2M, EF1alpha, GAPDH, GUSB, PPIA, RPL13A, RPLP0, TBP, UBC, YWHAZ and 18S rRNA) to assess their suitability as HGs in MSCs during adipogenic, osteogenic and chondrogenic differentiation. We demonstrated that many of the most popular HGs including 18S rRNA, B2M and ACTB were inadequate for normalization, whereas TBP/YWHAZ/GUSB were frequently identified among the best performers. Moreover, we showed the dramatic effects of suboptimal HGs choice on the quantification of cell differentiation markers, thus interfering with a reliable comparison of the lineage potential properties among various MSCs. Thus, in the emerging field of regenerative medicine, the identification of the most appropriate MSC source and cell line is so crucial for the treatment of patients that being inaccurate in the first step of the stem cell characterization can bring important consequences for the patients and for the promising potential of stem cell therapy.

Gene profiles during vestibular compensation in rats after unilateral labyrinthectomy

OBJECTIVES

The aims of this study were to determine changes in gene expression in the chronic state of vestibular compensation by microarray analysis and to validate the asymmetrical levels of gene expression in the ipsilateral and contralateral vestibular nucleus complexes (VNCs).

METHODS

Microarray analysis was used to examine the expression of genes up-regulated or down-regulated in the ipsilateral VNC at 1 and 7 days after unilateral labyrinthectomy. Up-regulated or down-regulated gene expression in the ipsilateral and contralateral VNCs was then validated by reverse transcriptase polymerase chain reaction at 1, 7, 14, and 28 days after labyrinthectomy.

RESULTS

The genes down-regulated at 1 day after labyrinthectomy and up-regulated at 7 days after labyrinthectomy as determined by microarray analysis and reverse transcriptase polymerase chain reaction were zinc finger protein 307, zinc metallopeptidase, P34, calcitonin receptor, insulin-like growth factor binding protein 5, GATA binding protein 3, and CD151. Expression of zinc finger protein 307, zinc metallopeptidase, P34, and calcitonin receptor was up-regulated even after 7 days in the contralateral VNC of rats that had labyrinthectomy.

CONCLUSIONS

This study demonstrated changes in gene expression in rats during the chronic phase of vestibular compensation after unilateral labyrinthectomy and provided profiles of these changes in gene expression.

Mesenchymal stem cells relieve fibrosis of Schistosoma japonicum-induced mouse liver injury

Mesenchymal stem cells (MSCs) have gained popularity for their potential as seed cells to treat various human diseases, including pathogenic infections. Schistosoma japonicum (S. japonicum) infection is characterized by formation of parasite egg granulomas and host liver fibrosis. MSCs have been proposed as useful treatments of S. japonicum infection, but the efficacy and underlying mechanisms remain unknown. Herein, we report that MSCs were able to ameliorate S. japonicum-induced liver injury in vivo and this effect was enhanced by combining MSCs with conventional drug praziquantel (PZQ). Kunming strains of mice were infected with S. japonicum and treated with vehicle, MSCs, PZQ or PZQ + MSCs. MSC treatment not only prolonged the survival time of infected mice but reduced egg granuloma diameter and decreased the concentrations of serum transforming growth factor-β1 and hyaluronic acid. MSC treatment also inhibited collagen deposition and reduced the expression of collagen type 3, α-smooth muscle actin and vimentin in infected mouse liver tissues. Collectively, our findings suggest that MSC treatment represents a novel therapeutic approach for S. japonicum-induced liver injury and fibrosis.

Discrimination of genotoxic and non-genotoxic hepatocarcinogens by statistical analysis based on gene expression profiling in the mouse liver as determined by quantitative real-time PCR

The general aim of the present study is to discriminate between mouse genotoxic and non-genotoxic hepatocarcinogens via selected gene expression patterns in the liver as analyzed by quantitative real-time PCR (qPCR) and statistical analysis. qPCR was conducted on liver samples from groups of 5 male, 9-week-old B6C3F(1) mice, at 4 and 48h following a single intraperitoneal administration of chemicals. We quantified 35 genes selected from our previous DNA microarray studies using 12 different chemicals: 8 genotoxic hepatocarcinogens (2-acetylaminofluorene, 2,4-diaminotoluene, diisopropanolnitrosamine, 4-dimethylaminoazobenzene, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, N-nitrosomorpholine, quinoline and urethane) and 4 non-genotoxic hepatocarcinogens (1,4-dichlorobenzene, dichlorodiphenyltrichloroethane, di(2-ethylhexyl)phthalate and furan). A considerable number of genes exhibited significant changes in their gene expression ratios (experimental group/control group) analyzed statistically by the Dunnett's test and Welch's t-test. Finally, we distinguished between the genotoxic and non-genotoxic hepatocarcinogens by statistical analysis using principal component analysis (PCA) of the gene expression profiles for 7 genes (Btg2, Ccnf, Ccng1, Lpr1, Mbd1, Phlda3 and Tubb2c) at 4h and for 12 genes (Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Gdf15, Lrp1, Mbd1, Phlda3, Plk2 and Tubb2c) at 48h. Seven major biological processes were extracted from the gene ontology analysis: apoptosis, the cell cycle, cell proliferation, DNA damage, DNA repair, oncogenes and tumor suppression. The major, biologically relevant gene pathway suggested was the DNA damage response pathway, resulting from signal transduction by a p53-class mediator leading to the induction of apoptosis. Eight genes (Aen, Bax, Btg2, Ccng1, Cdkn1a, Gdf15, Phlda3 and Plk2) that are directly associated with Trp53 contributed to the PCA. The current findings demonstrate a successful discrimination between genotoxic and non-genotoxic hepatocarcinogens, using qPCR and PCA, on 12 genes associated with a Trp53-mediated signaling pathway for DNA damage response at 4 and 48 h after a single administration of chemicals.

Selection of housekeeping genes for normalization of RT-PCR in hypoxic neural stem cells of rat in vitro

Gene expression analysis under various conditions using real-time reverse transcription polymerase chain reaction (RT-PCR) needs reliable control genes. Housekeeping genes are commonly used as the control. However, no validated housekeeping gene is available for study of hypoxic neural stem cell culture. To choose appropriate internal control genes, the expression of eight commonly used housekeeping genes was examined in rat neural stem cell model to find one or more stably expressed genes under hypoxic/ischemic conditions. Two genes, HPRT and RPL13A were identified as the most confidential housekeeping genes in this research by geNorm and NormFinder softwares. As a groundwork, the most stable housekeeping genes for neural stem cells under hypoxic/ischemic conditions are initially investigated and validated in this experiment, which might provide a better understanding for the gene expression study in ischemic and necrotic neural stem cell cultures or in ischemic diseases of the central nervous system (CNS).

GATA depletion impacts insulin-like growth factor 1 mRNA and protein levels in luteinizing porcine granulosa cells

GATA4 and GATA6 are zinc-finger transcription factors that regulate specific genes involved in steroidogenesis. Using RNA interference (RNAi)-mediated reduction of GATA4 and/or GATA6 with microarray analysis, we aimed to identify novel GATA target genes in luteinizing porcine granulosa cells under vehicle- and cAMP-treated conditions. Microarray analysis identified IGF1 mRNA to be cAMP- and GATA-responsive, and real-time PCR demonstrated that the cAMP-induced increase in IGF1 mRNA was reduced under conditions of GATA6 depletion and GATA4 plus GATA6 depletion, but not GATA4 depletion. Insulin-like growth factor 1 protein levels in media were also decreased by GATA6 or GATA4 plus GATA6 reduction. IGFBP2 and IGFBP4 mRNAs were increased and IGFBP5 mRNA decreased with vehicle and cAMP treatment under GATA4 plus GATA6 RNAi conditions. GATA6 reduction alone increased basal IGFBP4 and decreased IGFBP5 with both vehicle and cAMP, and GATA4 reduction alone lowered cAMP IGFBP5 levels with cAMP. No changes in IGFBP3 mRNA were observed with GATA reduction relative to the control RNAi condition. Levels of insulin-like growth factor binding proteins 2-5 in media as assessed by Western ligand blotting were not altered by GATA reduction. Electromobility gel shift assays with two GATA-containing oligonucleotides of the IGF1 5'-regulatory region showed GATA4 and GATA6 could bind the more proximal GATA-B site. These studies indicate that although GATA4 and GATA6 can bind the porcine IGF1 5'-region, GATA6 is functionally most important for cAMP-stimulated mRNA levels. Using microarray analysis, we identified other mRNAs that were altered by GATA-reduced conditions, including ALDH1, DIO2, and EDNRB. Our findings further support GATA as a coordinator of endocrine/paracrine/autocrine signals in the ovary.

Transcriptomic profile indicative of immunotoxic exposure: in vitro studies in peripheral blood mononuclear cells

Investigating the immunotoxic effects of exposure to chemicals usually comprises evaluation of weight and histopathology of lymphoid tissues, various lymphocyte parameters in the circulation, and immune function. Immunotoxicity assessment is time consuming in humans or requires a high number of animals, making it expensive. Furthermore, reducing the use of animals in research is an important ethical and political issue. Immunotoxicogenomics represents a novel approach to investigate immunotoxicity able of overcoming these limitations. The current research, embedded in the European Union project NewGeneris, aimed to retrieve gene expression profiles that are indicative of exposure to immunotoxicants. To this end, whole-genome gene expression was investigated in human peripheral blood mononuclear cells in response to in vitro exposure to a range of immunotoxic chemicals (4-hydroxy-2-nonenal, aflatoxin B1, benzo[a]pyrene, deoxynivalenol, ethanol, malondialdehyde, polychlorinated biphenyl 153, and 2,3,7,8-tetrachlorodibenzo-p-dioxin) and nonimmunotoxic chemicals (acrylamide, dimethylnitrosamine, 2-amino-3-methyl-3H-imidazo[4,5-F]quinoline, and 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine). Using Agilent oligonucleotide microarrays, whole-genome gene expression profiles were generated, which were analyzed using Genedata's Expressionist software. Using Recursive Feature Elimination and Support Vector Machine, a set of 48 genes was identified that distinguishes the immunotoxic from the nonimmunotoxic compounds. Analysis for enrichment of biological processes showed the gene set to be highly biologically and immunologically relevant. We conclude that we have identified a promising transcriptomic profile indicative of immunotoxic exposure.

Microarray analysis of relative gene expression stability for selection of internal reference genes in the rhesus macaque brain

Background

Normalization of gene expression data refers to the comparison of expression values using reference standards that are consistent across all conditions of an experiment. In PCR studies, genes designated as "housekeeping genes" have been used as internal reference genes under the assumption that their expression is stable and independent of experimental conditions. However, verification of this assumption is rarely performed. Here we assess the use of gene microarray analysis to facilitate selection of internal reference sequences with higher expression stability across experimental conditions than can be expected using traditional selection methods.

We recently demonstrated that relative gene expression from qRT-PCR data normalized using GAPDH, ALG9 and RPL13A expression values mirrored relative expression using quantile normalization in Robust Multichip Analysis (RMA) on the Affymetrix^® GeneChip^® rhesus Macaque Genome Array.

Having shown that qRT-PCR and Affymetrix^® GeneChip^® data from the same hormone replacement therapy (HRT) study yielded concordant results, we used quantile-normalized gene microarray data to identify the most stably expressed among probe sets for prospective internal reference genes across three brain regions from the HRT study and an additional study of normally menstruating rhesus macaques (cycle study). Gene selection was limited to 575 previously published human "housekeeping" genes. Twelve animals were used per study, and three brain regions were analyzed from each animal. Gene expression stabilities were determined using geNorm, NormFinder and BestKeeper software packages.

Results

Sequences co-annotated for ribosomal protein S27a (RPS27A), and ubiquitin were among the most stably expressed under all conditions and selection criteria used for both studies. Higher annotation quality on the human GeneChip^® facilitated more targeted analysis than could be accomplished using the rhesus GeneChip^®. In the cycle study, multiple probe sets annotated for actin, gamma 1 (ACTG1) showed high signal intensity and were among the most stably expressed.

Conclusions

Using gene microarray analysis, we identified genes showing high expression stability under various sex-steroid environments in different regions of the rhesus macaque brain. Use of quantile-normalized microarray gene expression values represents an improvement over traditional methods of selecting internal reference genes for PCR analysis.

Fine ambient air particulate matter exposure induces molecular alterations associated with vascular disease progression within plaques of atherosclerotic susceptible mice

Epidemiology studies have reported associations between increased mortality and morbidity with exposure to particulate air pollution, particularly within individuals with preexisting cardiovascular disease (CVD). Clinical and toxicological studies have provided evidence that exposure to ambient air particulate matter (PM) impacts CVD by increasing plaque size. It is unclear whether PM-induced increased plaque size is associated with molecular disease progression. This study examines molecular profiles within plaques recovered from ApoE(-/-) mice exposed to concentrated ambient air particles (CAPs) to determine whether pulmonary deposition of PM contributes to molecular alterations leading to vascular disease progression. Laser capture microdissection was used to recover atherosclerotic plaques from ApoE(-/-) male mice exposed daily for 5 mo to filtered air or CAPs. Alterations in mRNA expression was assessed in microdissected plaques of CAPs-exposed and air controls using the Affymetrix microarray platform. Bioinformatic analysis indicated alterations in 611 genes: 395 genes downregulated and 216 genes upregulated. Gene ontology revealed CAPs-induced changes to inflammation, proliferation, cell cycle, hematological system, and cardiovascular pathways. Quantitative reverse-transcription polymerase chain reaction (qRT-PCR) verified microarray data also revealing gene expression alterations undetected by the microarray analysis, i.e., decreased expression of alpha-actin for smooth muscle cells, and increased expression of the macrophage marker Cd68 and of beta-actin. Comparison of CAPs-induced gene expression profiles demonstrated consistency with previously published gene expression profiles in the ApoE(-/-) mouse model and humans associated with plaque progression. These results indicate that exposure to fine PM induces molecular alterations associated with vascular disease progression and provides insight into potential biological pathways responsible for this effect.

Kidney injury molecule-1 expression in rat proximal tubule after treatment with segment-specific nephrotoxicants: a tool for early screening of potential kidney toxicity

Dose-response expression of kidney injury molecule-1 (KIM-1) gene in kidney cortex and its correlation with morphology and traditional biomarkers of nephrotoxicity (plasma creatinine and blood urea nitrogen, BUN) or segment-specific marker of proximal tubule injury (kidney glutamine synthetase, GSK) were studied in male rats treated with proximal tubule segment-specific nephrotoxicants. These included hexachloro-1:3-butadiene (HCBD, S(3) segment-specific), potassium dichromate (chromate, S(1)-S(2) segment-specific), and cephaloridine (Cph, S(2) segment-specific). Rats were treated with a single intraperitoneal (ip) injection of HCBD 25, 50, and 100 mg/kg, subcutaneous (sc) injection of chromate 8, 12.5, and 25 mg/kg; or ip injection of Cph 250, 500, and 1,000 mg/kg. KIM-1 gene showed a dose-dependent up-regulation induced by all segment-specific nephrotoxicants. Interestingly, magnitude of the up-regulation reflected the severity of microscopic tubular changes (degeneration, necrosis, and regeneration). Even low-severity microscopic observations were evidenced by significant gene expression changes. Furthermore, KIM-1 showed significant up-regulation even in the absence of morphological changes. In contrast, traditional and specific markers demonstrated low sensitivity or specificity. In conclusion, this study suggested KIM-1 as a sensitive molecular marker of different levels of tubular injury, and it is likely to represent a potential tool for early screening of nephrotoxicants.

Comparison of lung cancer cell lines representing four histopathological subtypes with gene expression profiling using quantitative real-time PCR

BACKGROUND

Lung cancers are the most common type of human malignancy and are intractable. Lung cancers are generally classified into four histopathological subtypes: adenocarcinoma (AD), squamous cell carcinoma (SQ), large cell carcinoma (LC), and small cell carcinoma (SC). Molecular biological characterization of these subtypes has been performed mainly using DNA microarrays. In this study, we compared the gene expression profiles of these four subtypes using twelve human lung cancer cell lines and the more reliable quantitative real-time PCR (qPCR).

RESULTS

We selected 100 genes from public DNA microarray data and examined them by DNA microarray analysis in eight test cell lines (A549, ABC-1, EBC-1, LK-2, LU65, LU99, STC 1, RERF-LC-MA) and a normal control lung cell line (MRC-9). From this, we extracted 19 candidate genes. We quantified the expression of the 19 genes and a housekeeping gene, GAPDH, with qPCR, using the same eight cell lines plus four additional validation lung cancer cell lines (RERF-LC-MS, LC-1/sq, 86-2, and MS-1-L). Finally, we characterized the four subtypes of lung cancer cell lines using principal component analysis (PCA) of gene expression profiling for 12 of the 19 genes (AMY2A, CDH1, FOXG1, IGSF3, ISL1, MALL, PLAU, RAB25, S100P, SLCO4A1, STMN1, and TGM2). The combined PCA and gene pathway analyses suggested that these genes were related to cell adhesion, growth, and invasion. S100P in AD cells and CDH1 in AD and SQ cells were identified as candidate markers of these lung cancer subtypes based on their upregulation and the results of PCA analysis. Immunohistochemistry for S100P and RAB25 was closely correlated to gene expression.

CONCLUSIONS

These results show that the four subtypes, represented by 12 lung cancer cell lines, were well characterized using qPCR and PCA for the 12 genes examined. Certain genes, in particular S100P and CDH1, may be especially important for distinguishing the different subtypes. Our results confirm that qPCR and PCA analysis provide a useful tool for characterizing cancer cell subtypes, and we discuss the possible clinical applications of this approach.

Customized PCR-array analysis informed by gene-chip microarray and biological hypothesis reveals pathways involved in lung inflammatory response to titanium dioxide in pregnancy

Validation of gene-chip microarray results is one of the challenges in genomic studies. The successful use of a custom-designed 96-well polymerase chain reaction (PCR) array to study the unexpected inflammatory effect of environmental titanium dioxide (TiO2) particles on the lungs of pregnant mice, with similar results not seen in control mice, is reported. In our approach, selection of candidate genes for the custom PCR array was informed by prior gene-chip microarray profiling. Results demonstrated multiple upregulation of genes in the lungs of pregnant but not control mice produced by TiO2 exposure. Customized PCR array is a flexible tool that offers the ability to combine the "blind" genome-wide scan with a hypothesis-driven approach, by including both the "candidate" genes for validation positively identified by the microarray and biologically relevant "suspects" that failed to be found in the genomic data. Compared to conventional gene-by-gene qPCR or manufacturer-preset pathway kits, this technique provides a cost-effective and time-saving method of analysis and allows for a strong, easily detectable signal. Genes with confirmed differential expression were further used for pathway analysis and indicated involvement in several biologically relevant pathways including allergy mediator signaling in dendritic cells. Finally, an analytical network was created that will inform further mechanistic studies. The dual purpose of the work was to demonstrate that the novel custom PCR array is a convenient approach to validate the microarray results, and to obtain biologically significant data on TiO2-induced inflammation by following the PCR array with pathway analysis, which provided feasible hypotheses to support future experimental studies.

Validation of reference genes for quantitative expression analysis by real-time RT-PCR in Saccharomyces cerevisiae

Background

Real-time RT-PCR is the recommended method for quantitative gene expression analysis. A compulsory step is the selection of good reference genes for normalization. A few genes often referred to as HouseKeeping Genes (HSK), such as ACT1, RDN18 or PDA1 are among the most commonly used, as their expression is assumed to remain unchanged over a wide range of conditions. Since this assumption is very unlikely, a geometric averaging of multiple, carefully selected internal control genes is now strongly recommended for normalization to avoid this problem of expression variation of single reference genes. The aim of this work was to search for a set of reference genes for reliable gene expression analysis in Saccharomyces cerevisiae.

Results

From public microarray datasets, we selected potential reference genes whose expression remained apparently invariable during long-term growth on glucose. Using the algorithm geNorm, ALG9, TAF10, TFC1 and UBC6 turned out to be genes whose expression remained stable, independent of the growth conditions and the strain backgrounds tested in this study. We then showed that the geometric averaging of any subset of three genes among the six most stable genes resulted in very similar normalized data, which contrasted with inconsistent results among various biological samples when the normalization was performed with ACT1. Normalization with multiple selected genes was therefore applied to transcriptional analysis of genes involved in glycogen metabolism. We determined an induction ratio of 100-fold for GPH1 and 20-fold for GSY2 between the exponential phase and the diauxic shift on glucose. There was no induction of these two genes at this transition phase on galactose, although in both cases, the kinetics of glycogen accumulation was similar. In contrast, SGA1 expression was independent of the carbon source and increased by 3-fold in stationary phase.

Conclusion

In this work, we provided a set of genes that are suitable reference genes for quantitative gene expression analysis by real-time RT-PCR in yeast biological samples covering a large panel of physiological states. In contrast, we invalidated and discourage the use of ACT1 as well as other commonly used reference genes (PDA1, TDH3, RDN18, etc) as internal controls for quantitative gene expression analysis in yeast.

Glycogenome expression dynamics during mouse C2C12 myoblast differentiation suggests a sequential reorganization of membrane glycoconjugates

Background

Several global transcriptomic and proteomic approaches have been applied in order to obtain new molecular insights on skeletal myogenesis, but none has generated any specific data on glycogenome expression, and thus on the role of glycan structures in this process, despite the involvement of glycoconjugates in various biological events including differentiation and development. In the present study, a quantitative real-time RT-PCR technology was used to profile the dynamic expression of 375 glycogenes during the differentiation of C2C12 myoblasts into myotubes.

Results

Of the 276 genes expressed, 95 exhibited altered mRNA expression when C2C12 cells differentiated and 37 displayed more than 4-fold up- or down-regulations. Principal Component Analysis and Hierarchical Component Analysis of the expression dynamics identified three groups of coordinately and sequentially regulated genes. The first group included 12 down-regulated genes, the second group four genes with an expression peak at 24 h of differentiation, and the last 21 up-regulated genes. These genes mainly encode cell adhesion molecules and key enzymes involved in the biosynthesis of glycosaminoglycans and glycolipids (neolactoseries, lactoseries and ganglioseries), providing a clearer indication of how the plasma membrane and extracellular matrix may be modified prior to cell fusion. In particular, an increase in the quantity of ganglioside G_M3 at the cell surface of myoblasts is suggestive of its potential role during the initial steps of myogenic differentiation.

Conclusion

For the first time, these results provide a broad description of the expression dynamics of glycogenes during C2C12 differentiation. Among the 37 highly deregulated glycogenes, 29 had never been associated with myogenesis. Their biological functions suggest new roles for glycans in skeletal myogenesis.