Validation of differential gene expression in muscle engineered from rat groin adipose tissue by quantitative real-time PCR

Optimizations for identifying reference genes in bone and cartilage bioengineering

BACKGROUND

Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) remains one of the best-established techniques to assess gene expression patterns. However, appropriate reference gene(s) selection remains a critical and challenging subject in which inappropriate reference gene selction can distort results leading to false interpretations. To date, mixed opinions still exist in how to choose the most optimal reference gene sets in accodrance to the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guideline. Therefore, the purpose of this study was to investigate which schemes were the most feasible for the identification of reference genes in a bone and cartilage bioengineering experimental setting. In this study, rat bone mesenchymal stem cells (rBMSCs), skeletal muscle tissue and adipose tissue were utilized, undergoing either chondrogenic or osteogenic induction, to investigate the optimal reference gene set identification scheme that would subsequently ensure stable and accurate interpretation of gene expression in bone and cartilage bioengineering.

RESULTS

The stability and pairwise variance of eight candidate reference genes were analyzed using geNorm. The V_0.15- vs. V_min-based normalization scheme in rBMSCs had no significant effect on the eventual normalization of target genes. In terms of the muscle tissue, the results of the correlation of NF values between the V_0.15 and V_min schemes and the variance of target genes expression levels generated by these two schemes showed that different schemes do indeed have a significant effect on the eventual normalization of target genes. Three selection schemes were adopted in terms of the adipose tissue, including the three optimal reference genes (Opt₃), V_0.20 and V_min schemes, and the analysis of NF values with eventual normalization of target genes showed that the different selection schemes also have a significant effect on the eventual normalization of target genes.

CONCLUSIONS

Based on these results, the proposed cut-off value of Vn/n + 1 under 0.15, according to the geNorm algorithm, should be considered with caution. For cell only experiments, at least rBMSCs, a Vn/n + 1 under 0.15 is sufficient in RT-qPCR studies. However, when using certain tissue types such as skeletal muscle and adipose tissue the minimum Vn/n + 1 should be used instead as this provides a far superior mode of generating accurate gene expression results. We thus recommended that when the stability and variation of a candidate reference genes in a specific study is unclear the minimum Vn/n + 1 should always be used as this ensures the best and most accurate gene expression value is achieved during RT-qPCR assays.

Large-Volume Vascularized Muscle Grafts Engineered From Groin Adipose Tissue in Perfusion Bioreactor Culture

BACKGROUND

Muscle tissue engineering still remains a major challenge. An axial vascular pedicle and a perfusion bioreactor are necessary for the development and maintenance of a large-volume engineered muscle tissue to provide circulation within the construct. This study aimed to determine whether large-volume vascularized muscle-like constructs could be made from rat groin adipose tissue in a perfusion bioreactor.

METHODS

Epigastric adipofascial flaps based on the inferior superficial epigastric vessels were elevated bilaterally in male Lewis rats and connected to the bioreactor. The system was run using a cable pump and filled with myogenic differentiation medium in the perfusion bioreactor for 1, 3, 5, or 7 weeks. The resulting tissue constructs were characterized with respect to the morphology and muscle-related expression of genes and proteins.

RESULTS

The histological examination demonstrated intact muscle-like tissue fibers; myogenesis was verified by the expression of myosin, MADS box transcription enhancer factor 2 D, desmin-a disintegrin and metalloproteinase domain (ADAM) 12-and M-cadherin using reverse transcription-polymerase chain reaction. Western blot analysis for desmin, MyoD1, N-cadherin, and ADAM12 was performed to verify the myogenic phenotype of the extracted differentiated tissue and prove the formation of muscle-like constructs.

CONCLUSIONS

A large-volume vascularized muscle tissue could be engineered in a perfusion bioreactor. The resulting tissue had muscle-like histological features and expressed muscle-related genes and proteins, indicating that the trans-differentiation of adipose tissue into muscle tissue occurred.

Identification of Suitable Reference Genes for Quantitative Gene Expression Analysis in Innervated and Denervated Adipose Tissue from Cafeteria Diet-Fed Rats

Our previous studies show that cafeteria diet increases body adiposity, plasma insulin levels, and sympathetic activity to brown adipose tissue (BAT) and white adipose tissue (WAT) of Wistar rats, leading to rapid and progressive changes in the metabolic profile. The identification of suitable reference genes that are not affected by the experimental conditions is a critical step in accurate normalization of the reverse transcription quantitative real-time PCR (qRT-PCR), a commonly used assay to elucidate changes in the gene expression profile. In the present study, the effects of the cafeteria diet and sympathetic innervation on the gene expression of adrenoceptor beta 3 (Adrb3) from BAT and WAT were assessed using one of the most stable and one of the least stable genes as normalizers. Rats were fed the cafeteria diet and on the 17th day, interscapular BAT or retroperitoneal WAT was denervated and, 7 days after surgery, the contralateral innervated tissue was used as control. Ten reference genes were evaluated (18S, B2m, Actb, CypA, Gapdh, Hprt1, Rpl32, Tbp, Ubc, and Ywhaz) and ranked according to their stability using the following algorithms: geNorm, NormFinder, BestKeeper, and comparative delta threshold cycle (ΔC _t ) method. According to the algorithms employed, the normalization of Adrb3 expression by the least stable genes produced opposite results compared with the most stable genes and literature data. In cafeteria and control diet-fed rats, the three most stable genes were Hprt1, Tbp, and Rpl32 for interscapular BAT and Tbp, B2m, and Hprt1 for retroperitoneal WAT, while the least stable genes were 18S, Actb, and Gapdh for both tissues.

Selection of Reference Genes for Expression Studies of Xenobiotic Adaptation in Tetranychus urticae

Quantitative real-time PCR (qRT-PCR) is an extensively used, high-throughput method to analyze transcriptional expression of genes of interest. An appropriate normalization strategy with reliable reference genes is required for calculating gene expression across diverse experimental conditions. In this study, we aim to identify the most stable reference genes for expression studies of xenobiotic adaptation in Tetranychus urticae, an extremely polyphagous herbivore causing significant yield reduction of agriculture. We chose eight commonly used housekeeping genes as candidates. The qRT-PCR expression data for these genes were evaluated from seven populations: a susceptible and three acaricide resistant populations feeding on lima beans, and three other susceptible populations which had been shifted host from lima beans to three other plant species. The stability of the candidate reference genes was then assessed using four different algorithms (comparative ΔCt method, geNorm, NormFinder, and BestKeeper). Additionally, we used an online web-based tool (RefFinder) to assign an overall final rank for each candidate gene. Our study found that CycA and Rp49 are best for investigating gene expression in acaricide susceptible and resistant populations. GAPDH, Rp49, and Rpl18 are best for host plant shift studies. And GAPDH and Rp49 were the most stable reference genes when investigating gene expression under changes in both experimental conditions. These results will facilitate research in revealing molecular mechanisms underlying the xenobiotic adaptation of this notorious agricultural pest.

Identification and de novo sequencing of housekeeping genes appropriate for gene expression analyses in farmed maraena whitefish (Coregonus maraena) during crowding stress

Maraena whitefish (Coregonus maraena; synonym Coregonus lavaretus f. balticus) is a high-quality food fish in the Southern Baltic Sea belonging to the group of salmonid fishes. Coregonus sp. is successfully kept in aquaculture throughout northern Europe (e.g. in Finland, Germany, Russia) and North America. In this regard, the molecular and immunological characterisation of stress response in maraena whitefish contributes to the development of robust and fast-growing maraena whitefish breeding strains for aquaculture. Thus, in the present study, the potential housekeeping genes beta actin (ACTB), elongation factor 1 alpha (EEF1A1), glyceraldehydes-3-phosphate dehydrogenase (GAPDH), ribosomal protein 9 (RPL9), ribosomal protein 32 (RPL32) and ribosomal protein S20 (RPS20) were de novo sequenced and tested concerning their applicability as reference genes in quantitative real-time PCR (qPCR) in maraena whitefish under different stocking densities. For this purpose, tissue samples of liver, kidney, gills, head kidney, skin, adipose tissue, heart and dorsal fin were investigated. qPCR data were analysed with Normfinder tool to determine gene expression stability. DNA sequencing exposed transcribed paralogous EEF1A1A and EEF1A1B genes differing in their putative protein structure. Normfinder analysis revealed RPL9 and RPL32 as most stable, GAPDH and ACTB as least stable genes for qPCR analyses, respectively. This is the first study that provides a subset of seven de novo sequenced housekeeping genes usable as reference genes in studies of stress response in maraena whitefish.

Expression profiling in Bemisia tabaci under insecticide treatment: indicating the necessity for custom reference gene selection

Finding a suitable reference gene is the key for qRT-PCR analysis. However, none of the reference gene discovered thus far can be utilized universally under various biotic and abiotic experimental conditions. In this study, we further examine the stability of candidate reference genes under a single abiotic factor, insecticide treatment. After being exposed to eight commercially available insecticides, which belong to five different classes, the expression profiles of eight housekeeping genes in the sweetpotato whitefly, Bemisia tabaci, one of the most invasive and destructive pests in the world, were investigated using qRT-PCR analysis. In summary, elongation factor 1α (EF1α), α-tubulin (TUB1α) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) were identified as the most stable reference genes under the insecticide treatment. The initial assessment of candidate reference genes was further validated with the expression of two target genes, a P450 (Cyp6cm1) and a glutathione S-transferase (GST). However, ranking of reference genes varied substantially among intra- and inter-classes of insecticides. These combined data strongly suggested the necessity of conducting custom reference gene selection designed for each and every experimental condition, even when examining the same abiotic or biotic factor.

Selection of valid reference genes for expression studies of hepatic cell lines under IFN-α treatment

The proper selection of reference genes to normalize the quantitative real-time PCR (RT-qPCR) results under particular experimental conditions is crucial for validation of the gene quantification data. Herein, using SYBR green RT-qPCR, five reference genes (GAPDH, ACTB, HMBS, HPRT-1 and TBP) were evaluated to determine the most stable reference genes in hepatic cell lines (Huh-7 and HepG(2)) under IFN-α treatment conditions. Analyses by geNorm program ranked GAPDH and HPRT-1 in Huh-7 and that of ACTB and HMBS in HepG(2) cells as the most stable reference genes under IFN-α treatment. While, same reference gene pairs were ranked by NormFinder program in Huh-7 cells, GAPDH was assessed as the most stable gene in HepG(2) group by this program, implying the importance of the employed algorithm in comparative interpretation of the data. Finally, cumulative analyses by one-way ANOVA, geNorm and NormFinder programs indicated that use of two reference genes (HMBS and GAPDH) in Huh-7 and three (HMBS, ACTB and GAPDH) in HepG(2) cells would greatly improve the normalization of the RT-qPCR data under IFN-α. Data presented in this paper will aid the selection of the most stable reference genes in RT-qPCR studies on evaluation of hepatic viral proteins and IFN pathway.

Identification of condition-specific reference genes from microarray data for locusts exposed to hypobaric hypoxia

Real-time quantitative polymerase chain reaction (qPCR) is a routine and robust approach for measuring gene expression. The stability of reference genes in qPCR is crucial for the accurate quantification of gene expression. To provide reliable reference genes for studying the transcriptional responses of locust muscles to hypobaric hypoxia, we first examined the gene expression stability of the frequently used housekeeping genes 18S, GAPDH, and β-actin. However, the expression of these three housekeeping genes was influenced by hypobaric hypoxia. Consequently, we identified five novel candidate reference genes from the locust microarray data. The gene expression stability of the five candidates, together with the three classical housekeeping genes, were evaluated using two distinct algorithms implemented in geNorm and NormFinder. GeNorm identified Ach (acetyl-CoA hydrolase) and Pgp (phosphoglycolate phosphatase-like) as the most stable genes and NormFinder further distinguished Ach as the most stable one. The validity of Ach as a reference gene was confirmed through comparison with 18S. This study exemplifies the necessity of validating reference genes before their application and the feasibility of identifying condition-specific reference genes from large-scale gene expression data.