1
|
Age-Invariant Genes: Multi-Tissue Identification and Characterization of Murine Reference Genes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.04.09.588721. [PMID: 38645168 PMCID: PMC11030416 DOI: 10.1101/2024.04.09.588721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/23/2024]
Abstract
Studies of the aging transcriptome focus on genes that change with age. But what can we learn from age-invariant genes-those that remain unchanged throughout the aging process? These genes also have a practical application: they serve as reference genes (often called housekeeping genes) in expression studies. Reference genes have mostly been identified and validated in young organisms, and no systematic investigation has been done across the lifespan. Here, we build upon a common pipeline for identifying reference genes in RNA-seq datasets to identify age-invariant genes across seventeen C57BL/6 mouse tissues (brain, lung, bone marrow, muscle, white blood cells, heart, small intestine, kidney, liver, pancreas, skin, brown, gonadal, marrow, and subcutaneous adipose tissue) spanning 1 to 21+ months of age. We identify 9 pan-tissue age-invariant genes and many tissue-specific age-invariant genes. These genes are stable across the lifespan and are validated in independent bulk RNA-seq datasets and RT-qPCR. We find age-invariant genes have shorter transcripts on average and are enriched for CpG islands. Interestingly, pathway enrichment analysis for age-invariant genes identifies an overrepresentation of molecular functions associated with some, but not all, hallmarks of aging. Thus, though hallmarks of aging typically involve changes in cell maintenance mechanisms, select genes associated with these hallmarks resist fluctuations in expression with age. Finally, our analysis concludes no classical reference gene is appropriate for aging studies in all tissues. Instead, we provide tissue-specific and pan-tissue genes for assays utilizing reference gene normalization (i.e., RT-qPCR) that can be applied to animals across the lifespan.
Collapse
|
2
|
Selection and validation of novel stable reference genes for qPCR analysis in EMT and MET. Exp Cell Res 2023; 428:113619. [PMID: 37146958 DOI: 10.1016/j.yexcr.2023.113619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2023] [Revised: 04/18/2023] [Accepted: 04/29/2023] [Indexed: 05/07/2023]
Abstract
Quantitative real-time polymerase chain reaction is a powerful tool for quantifying gene expression. The relative quantification relies on normalizing the data to reference genes or internal controls not modulated by the experimental conditions. The most widely used internal controls occasionally show changed expression patterns in different experimental settings, such as the mesenchymal to epithelial transition. Thus, identifying appropriate internal controls is of utmost importance. We analyzed multiple RNA-Seq datasets using a combination of statistical approaches such as percent relative range and coefficient of variance to define a list of candidate internal control genes, which was then validated experimentally and by using in silico analyses as well. We identified a group of genes as strong internal control candidates with high stability compared to the classical ones. We also presented evidence for the superiority of the percent relative range method for calculating expression stability in data sets with larger sample sizes. We used multiple methods to analyze data collected from several RNA-Seq datasets; we identified Rbm17 and Katna1 as the most stable reference genes in EMT/MET studies. The percent relative range approach surpasses other methods when analyzing datasets of larger sample sizes.
Collapse
|
3
|
A Multi-Tissue Gene Expression Atlas of Water Buffalo ( Bubalus bubalis) Reveals Transcriptome Conservation between Buffalo and Cattle. Genes (Basel) 2023; 14:genes14040890. [PMID: 37107649 PMCID: PMC10137413 DOI: 10.3390/genes14040890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 04/04/2023] [Accepted: 04/07/2023] [Indexed: 04/29/2023] Open
Abstract
We generated 73 transcriptomic data of water buffalo, which were integrated with publicly available data in this species, yielding a large dataset of 355 samples representing 20 major tissue categories. We established a multi-tissue gene expression atlas of water buffalo. Furthermore, by comparing them with 4866 cattle transcriptomic data from the cattle genotype-tissue expression atlas (CattleGTEx), we found that the transcriptomes of the two species exhibited conservation in their overall gene expression patterns, tissue-specific gene expression and house-keeping gene expression. We further identified conserved and divergent expression genes between the two species, with the largest number of differentially expressed genes found in the skin, which may be related to structural and functional differences in the skin of the two species. This work provides a source of functional annotation of the buffalo genome and lays the foundations for future genetic and evolutionary studies in water buffalo.
Collapse
|
4
|
Comparative analysis of novel and common reference genes in adult tissues of the mussel Mytilus galloprovincialis. BMC Genomics 2022; 23:349. [PMID: 35524177 PMCID: PMC9077915 DOI: 10.1186/s12864-022-08553-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2021] [Accepted: 04/14/2022] [Indexed: 11/18/2022] Open
Abstract
Background Real-time quantitative PCR is a widely used method for gene expression analyses in various organisms. Its accuracy mainly relies on the correct selection of reference genes. Any experimental plan involving real-time PCR needs to evaluate the characteristics of the samples to be examined and the relative stability of reference genes. Most studies in mollusks rely on reference genes commonly used in vertebrates. Results In this study, we focused on the transcriptome of the bivalve mollusk Mytilus galloprovincialis in physiological state to identify suitable reference genes in several adult tissues. Candidate genes with highly stable expression across 51 RNA-seq datasets from multiple tissues were selected through genome-wide bioinformatics analysis. This approach led to the identification of three genes (Rpl14, Rpl32 and Rpl34), whose suitability was evaluated together with 7 other reference genes commonly reported in literature (Act, Cyp-A, Ef1α, Gapdh, 18S, 28S and Rps4). The stability analyses performed with geNorm, NormFinder and Bestkeeper identified specific either single or pairs of genes suitable as references for gene expression analyses in specific tissues and revealed the Act/Cyp-A pair as the most appropriate to analyze gene expression across different tissues. Conclusion Mytilus galloprovincialis is a model system increasingly used in ecotoxicology and molecular studies. Our transcriptome-wide approach represents the first comprehensive investigation aimed at the identification of suitable reference genes for expression studies in this species. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08553-1.
Collapse
|
5
|
Transcriptional atlas analysis from multiple tissues reveals the expression specificity patterns in beef cattle. BMC Biol 2022; 20:79. [PMID: 35351103 PMCID: PMC8966188 DOI: 10.1186/s12915-022-01269-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Accepted: 03/03/2022] [Indexed: 11/12/2022] Open
Abstract
Background A comprehensive analysis of gene expression profiling across tissues can provide necessary information for an in-depth understanding of their biological functions. We performed a large-scale gene expression analysis and generated a high-resolution atlas of the transcriptome in beef cattle. Results Our transcriptome atlas was generated from 135 bovine tissues in adult beef cattle, covering 51 tissue types of major organ systems (e.g., muscular system, digestive system, immune system, reproductive system). Approximately 94.76% of sequencing reads were successfully mapped to the reference genome assembly ARS-UCD1.2. We detected a total of 60,488 transcripts, and 32% of them were not reported before. We identified 2654 housekeeping genes (HKGs) and 477 tissue-specific genes (TSGs) across tissues. Using weighted gene co-expression network analysis, we obtained 24 modules with 237 hub genes (HUBGs). Functional enrichment analysis showed that HKGs mainly maintain the basic biological activities of cells, while TSGs were involved in tissue differentiation and specific physiological processes. HKGs in bovine tissues were more conserved in terms of expression pattern as compared to TSGs and HUBGs among multiple species. Finally, we obtained a subset of tissue-specific differentially expressed genes (DEGs) between beef and dairy cattle and several functional pathways, which may be involved in production and health traits. Conclusions We generated a large-scale gene expression atlas across the major tissues in beef cattle, providing valuable information for enhancing genome assembly and annotation. HKGs, TSGs, and HUBGs further contribute to better understanding the biology and evolution of multiple tissues in cattle. DEGs between beef and dairy cattle also fill in the knowledge gaps about differential transcriptome regulation of bovine tissues underlying economically important traits. Supplementary Information The online version contains supplementary material available at 10.1186/s12915-022-01269-4.
Collapse
|
6
|
Understanding the constitutive presentation of MHC class I immunopeptidomes in primary tissues. iScience 2022; 25:103768. [PMID: 35141507 PMCID: PMC8810409 DOI: 10.1016/j.isci.2022.103768] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Revised: 06/15/2021] [Accepted: 01/11/2022] [Indexed: 12/20/2022] Open
Abstract
Understanding the molecular principles that govern the composition of the MHC-I immunopeptidome across different primary tissues is fundamentally important to predict how T cells respond in different contexts in vivo. Here, we performed a global analysis of the MHC-I immunopeptidome from 29 to 19 primary human and mouse tissues, respectively. First, we observed that different HLA-A, HLA-B, and HLA-C allotypes do not contribute evenly to the global composition of the MHC-I immunopeptidome across multiple human tissues. Second, we found that tissue-specific and housekeeping MHC-I peptides share very distinct properties. Third, we discovered that proteins that are evolutionarily hyperconserved represent the primary source of the MHC-I immunopeptidome at the organism-wide scale. Fourth, we uncovered new components of the antigen processing and presentation network, including the carboxypeptidases CPE, CNDP1/2, and CPVL. Together, this study opens up new avenues toward a system-wide understanding of antigen presentation in vivo across mammalian species. Tissue-specific and housekeeping MHC class I peptides share distinct properties HLA-A, HLA-B, and HLA-C allotypes contribute very unevenly to the pool of class I peptides MHC-I immunopeptidomes are represented by evolutionarily conserved proteins An extended antigen processing and presentation pathway is uncovered
Collapse
|
7
|
Unraveling the Spatiotemporal Distribution of VPS13A in the Mouse Brain. Int J Mol Sci 2021; 22:ijms222313018. [PMID: 34884823 PMCID: PMC8657609 DOI: 10.3390/ijms222313018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Revised: 11/22/2021] [Accepted: 11/25/2021] [Indexed: 11/27/2022] Open
Abstract
Loss-of-function mutations in the human vacuolar protein sorting the 13 homolog A (VPS13A) gene cause Chorea-acanthocytosis (ChAc), with selective degeneration of the striatum as the main neuropathologic feature. Very little is known about the VPS13A expression in the brain. The main objective of this work was to assess, for the first time, the spatiotemporal distribution of VPS13A in the mouse brain. We found VPS13A expression present in neurons already in the embryonic stage, with stable levels until adulthood. VPS13A mRNA and protein distributions were similar in the adult mouse brain. We found a widespread VPS13A distribution, with the strongest expression profiles in the pons, hippocampus, and cerebellum. Interestingly, expression was weak in the basal ganglia. VPS13A staining was positive in glutamatergic, GABAergic, and cholinergic neurons, but rarely in glial cells. At the cellular level, VPS13A was mainly located in the soma and neurites, co-localizing with both the endoplasmic reticulum and mitochondria. However, it was not enriched in dendritic spines or the synaptosomal fraction of cortical neurons. In vivo pharmacological modulation of the glutamatergic, dopaminergic or cholinergic systems did not modulate VPS13A concentration in the hippocampus, cerebral cortex, or striatum. These results indicate that VPS13A has remarkable stability in neuronal cells. Understanding the distinct expression pattern of VPS13A can provide relevant information to unravel pathophysiological hallmarks of ChAc.
Collapse
|
8
|
Identification of Reference Genes for Expression Studies in the Whole-Blood from Three Cattle Breeds under Two States of Livestock Weather Safety. Animals (Basel) 2021; 11:ani11113073. [PMID: 34827805 PMCID: PMC8614315 DOI: 10.3390/ani11113073] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2021] [Revised: 09/08/2021] [Accepted: 09/14/2021] [Indexed: 01/03/2023] Open
Abstract
Simple Summary Reductions in the fertility, body weight, and growth rate of cattle across the world are associated with the global warming phenomenon. Developing optimal management strategies is an important aspect of breeding programs for different breeds. Blood tissue undergoes dramatic physiological and metabolic changes during heat stress conditions, which involves the expression and regulation of a great number of genes across species. Real-time quantitative PCR (qPCR) is a method for the rapid and reliable quantification of mRNA transcription. Reference genes are used to normalize mRNA levels between different samples. Thus, the selection of high-quality reference genes is necessary for the interpretation of data generated by real-time PCR. Abstract Real-time PCR is widely used to study the relative abundance of mRNA due to its specificity, sensitivity, and repeatability quantification. However, relative quantification requires a reference gene, which should be stable in its expression, showing lower variation by experimental conditions or tissues. The aim of this study was to evaluate the stability of the expression of five commonly used reference genes (actb, ywhaz, b2m, sdha, and 18s rRNA) at different physiological stages (alert and emergency) in three different cattle breeds. In this study, five genes (actb, ywhaz, b2m, sdha, and 18s rRNA) were selected as candidate reference genes for expression studies in the whole blood from three cattle breeds (Romosinuano, Gyr, and Brahman) under heat stress conditions. The transcription stability of the candidate reference genes was evaluated using geNorm and NormFinder. The results showed that actb, 18SrRNA, and b2m expression were the most stable reference genes for whole blood of Gyr and Brahman breeds under two states of livestock weather safety (alert and emergency). Meanwhile, actb, b2m, and ywhaz were the most stable reference genes for the Romosinuano breed.
Collapse
|
9
|
HRT Atlas v1.0 database: redefining human and mouse housekeeping genes and candidate reference transcripts by mining massive RNA-seq datasets. Nucleic Acids Res 2021; 49:D947-D955. [PMID: 32663312 PMCID: PMC7778946 DOI: 10.1093/nar/gkaa609] [Citation(s) in RCA: 88] [Impact Index Per Article: 29.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2020] [Accepted: 07/08/2020] [Indexed: 12/18/2022] Open
Abstract
Housekeeping (HK) genes are constitutively expressed genes that are required for the maintenance of basic cellular functions. Despite their importance in the calibration of gene expression, as well as the understanding of many genomic and evolutionary features, important discrepancies have been observed in studies that previously identified these genes. Here, we present Housekeeping and Reference Transcript Atlas (HRT Atlas v1.0, www.housekeeping.unicamp.br) a web-based database which addresses some of the previously observed limitations in the identification of these genes, and offers a more accurate database of human and mouse HK genes and transcripts. The database was generated by mining massive human and mouse RNA-seq data sets, including 11 281 and 507 high-quality RNA-seq samples from 52 human non-disease tissues/cells and 14 healthy tissues/cells of C57BL/6 wild type mouse, respectively. User can visualize the expression and download lists of 2158 human HK transcripts from 2176 HK genes and 3024 mouse HK transcripts from 3277 mouse HK genes. HRT Atlas also offers the most stable and suitable tissue selective candidate reference transcripts for normalization of qPCR experiments. Specific primers and predicted modifiers of gene expression for some of these HK transcripts are also proposed. HRT Atlas has also been integrated with a regulatory elements resource from Epiregio server.
Collapse
|
10
|
Selection and validation of optimal endogenous reference genes for analysis of quantitative PCR in four tissues pathologically associated with Kidney-yang deficiency syndrome following influenza A infection. Exp Ther Med 2020; 20:244. [PMID: 33178342 PMCID: PMC7651868 DOI: 10.3892/etm.2020.9374] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2019] [Accepted: 07/22/2020] [Indexed: 11/06/2022] Open
Abstract
Kidney-yang deficiency syndrome (KYDS) infected with the influenza virus is a suitable model to imitate a population at high-risk to influenza infection with a high rate of morbidity and mortality. However, the specific molecular mechanisms underlying this disease remain unclear. A stable reference gene is essential as an internal control for molecular biology research of this condition. Reverse-transcription-quantitative PCR (RT-qPCR) is considered an extremely sensitive technique used for absolute and relative quantification of target genes transcript levels. To accurately estimate the relative expression of genes in cells from mice with KYDS in response to infection with influenza A virus subtype H1N1 (A/H1N1) virus using RT-qPCR, it is necessary to identify suitable reference genes. In the present study, analysis of 10 reference genes (Act-β, β2m, GAPDH, Gusβ, Tubα, Grcc10, Eif4h, Rnf187, Nedd8 and Ywhae) was performed across a set of 4 tissue types: Lung; heart; liver; and kidney. KYDS mice were inoculated with A/H1N1 virus or a mock control. For analysis, geNorm, BestKeeper, NormFinder, and Bio-Rad Maestro™ statistical programs were used for the estimation of the stability of the reference genes. The results were authenticated through extended experimental settings using a group of 10 samples, parallel to 3 additional innate immune system-associated genes of the host, TLR3, TLR7 and RIG-I, which were also analyzed using the same algorithms. From the 4 algorithms, taking into account the joint analyses of the ranking order outputs, the 2 genes Ywhae and Nedd8 were identified to be the most stable for mice with KYDS following infection with A/H1N1 virus. In contrast, the least stable genes in all 4 tissues were GAPDH and β2m. These results may affect the choice of reference genes in future studies that use RT-qPCR analysis of target genes in experimental conditions, such as mice with KYDS infected with influenza A virus.
Collapse
|
11
|
A Preliminary Study on the Characteristics of microRNAs in Ovarian Stroma and Follicles of Chuanzhong Black Goat during Estrus. Genes (Basel) 2020; 11:genes11090970. [PMID: 32825655 PMCID: PMC7564575 DOI: 10.3390/genes11090970] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Revised: 08/16/2020] [Accepted: 08/19/2020] [Indexed: 02/08/2023] Open
Abstract
microRNAs (miRNAs) play a significant role in ovarian follicular maturity, but miRNA expression patterns in ovarian stroma (OS), large follicles (LF), and small follicles (SF) have been rarely explored. We herein aimed to identify miRNAs, their target genes and signaling pathways, as well as their interaction networks in OS, LF, and SF of Chuanzhong black goats at the estrus phase using small RNA-sequencing. We found that the miRNA expression profiles of LF and SF were more similar than those of OS—32, 16, and 29 differentially expressed miRNAs were identified in OS vs. LF, OS vs. SF, and LF vs. SF, respectively. Analyses of functional enrichment and the miRNA-targeted gene interaction network suggested that miR-182 (SMC3), miR-122 (SGO1), and miR-206 (AURKA) were involved in ovarian organogenesis and hormone secretion by oocyte meiosis. Furthermore, miR-202-5p (EREG) and miR-485-3p (FLT3) were involved in follicular maturation through the MAPK signaling pathway, and miR-2404 (BMP7 and CDKN1C) played a key role in follicular development through the TGF-β signaling pathway and cell cycle; nevertheless, further research is warranted. To our knowledge, this is the first study to investigate miRNA expression patterns in OS, LF, and SF of Chuanzhong black goats during estrus. Our findings provide a theoretical basis to elucidate the role of miRNAs in follicular maturation. These key miRNAs might provide candidate biomarkers for the diagnosis of follicular maturation and will assist in developing new therapeutic targets for female goat infertility.
Collapse
|
12
|
Abstract
Approximately 5% of all cases of papillary thyroid cancer (PTC) are inherited. However, the susceptibility gene(s) for nonsyndromic familial PTC (FPTC) remain unclear. We performed whole genome sequencing of peripheral blood DNA samples from two affected family members with PTC. CHEK2 transcript expression and the protein levels of CHK2 and p53 were evaluated in the thyroid tissues from two affected members of the kindred and sporadic PTC cases. The entire CHEK2 coding sequence was examined by Sanger sequencing in blood DNA samples from 242 sporadic PTC patients. We identified a novel heterozygous germline mutation in CHEK2 (c.417C→A) that was detected in all available affected members of a kindred with FPTC. This variant was found in only 1 out of 264,200 persons in the Genome Aggregation Database and the NHLBI Trans-Omics for Precision Medicine program. The CHEK2 c.417C→A variant introduces a premature termination codon (Y139X). We found reduced CHK2 protein expression in tumor samples from the two patients who carried the variant as compared with sporadic cases without the variant. The Y139X loss-of-function variant led to reduced p53 phosphorylation and decreased p53 protein level. In addition, two rare missense variants (R180C and H371Y) in CHEK2 were identified in 5 (2%) of 242 patients with sporadic PTC. Our findings suggest that the CHEK2 Y139X variant may be associated with FPTC.
Collapse
|
13
|
An Efficient Exogenous Gene Insertion Site in CHO Cells with High Transcription Level to Enhance AID-Induced Mutation. Biotechnol J 2020; 15:e1900313. [PMID: 31975519 DOI: 10.1002/biot.201900313] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Revised: 12/18/2019] [Indexed: 11/10/2022]
Abstract
Antibodies have been extensively used for the purpose of scientific research, clinical diagnosis, and therapy. Combination of in vitro somatic hypermutation and mammalian cell surface display has been an efficient technology for antibody or other proteins optimization, in which the efficiency of activation-induced cytidine deaminase (AID) mutations in genes is one of the most important key factors. Gene transcriptional level has been found to be positively proportional to AID-induced mutation frequency. Thus, construction of the cell clone bearing a gene of interest (GOI) with high transcription level can increase AID-induced mutations. In this study, a retargetable gene cassette is inserted onto predetermined chromosome site (ywhae gene site) which is among the genes with the highest as well as stable transcription, and is found that one subsite is suitable to be retargeted for efficient protein display in Chinese hamster ovary (CHO) cells. The resultant cell clone (T31) has higher and more stable transcription/expression than CHO-puro clone which was previously established through the strategy of random insertion followed by a high-throughput selection. It also possesses a significantly higher mutation frequency to GOI than CHO-puro cells; thus, it is a better clone for the in vitro improvement of antibody affinity, and probably other properties.
Collapse
|
14
|
Comprehensive Profiles of mRNAs and miRNAs Reveal Molecular Characteristics of Multiple Organ Physiologies and Development in Pigs. Front Genet 2019; 10:756. [PMID: 31552085 PMCID: PMC6737989 DOI: 10.3389/fgene.2019.00756] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2019] [Accepted: 07/17/2019] [Indexed: 12/13/2022] Open
Abstract
The pig (Sus scrofa) is not only an important livestock animal but also widely used as a biomedical model. However, the understanding of the molecular characteristics of organs and of the developmental skeletal muscle of the pig is severely limited. Here, we performed a comprehensive transcriptome profiling of mRNAs and miRNAs across nine tissues and three skeletal muscle developmental stages in the Guizhou miniature pig. The reproductive organs (ovary and testis) had greater transcriptome complexity and activity than other tissues, and the highest transcriptome similarity was between skeletal muscle and heart (R = 0.79). We identified 1,819 mRNAs and 96 miRNAs to be tissue-specific in nine organs. Testis had the largest number of tissue-specific mRNAs (992) and miRNAs (40). Only 15 genes and two miRNAs were specifically expressed in skeletal muscle and fat, respectively. During postnatal skeletal muscle development, the mRNAs associated with focal adhesion, Notch signaling, protein digestion, and absorption pathways were up-regulated from D0 to D30 and then down-regulated from D30 and D240, while genes with opposing expression patterns were significantly enriched in the oxidative phosphorylation and proteasome pathways. The miRNAs mainly regulated genes associated with insulin, Wnt, fatty acid biosynthesis, Notch, MAPK, TGF-beta, insulin secretion, ECM-receptor interaction, focal adhesion, and calcium signaling pathways. We also identified 37 new miRNA-mRNA interaction pairs involved in skeletal muscle development. Overall, our data not only provide a rich resource for understanding pig organ physiology and development but also aid the study of the molecular functions of mRNA and miRNA in mammals.
Collapse
|
15
|
Defining housekeeping genes suitable for RNA-seq analysis of the human allograft kidney biopsy tissue. BMC Med Genomics 2019; 12:86. [PMID: 31208411 PMCID: PMC6580566 DOI: 10.1186/s12920-019-0538-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2018] [Accepted: 05/24/2019] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND RNA-seq is poised to play a major role in the management of kidney transplant patients. Rigorous definition of housekeeping genes (HKG) is essential for further progress in this field. Using single genes or a limited set HKG is inherently problematic since their expression might be altered by specific diseases in the patients being studied. METHODS To generate a HKG set specific for kidney transplantation, we performed RNA-sequencing from renal allograft biopsies collected in a variety of clinical settings. Various normalization methods were applied to identify transcripts that had a coefficient of variation of expression that was below the 2nd percentile across all samples, and the corresponding genes were designated as housekeeping genes. Comparison with transcriptomic data from the Gene Expression Omnibus (GEO) database, pathway analysis and molecular biological functions were utilized to validate the housekeeping genes set. RESULTS We have developed a bioinformatics solution to this problem by using nine different normalization methods to derive large HKG gene sets from a RNA-seq data set of 47,611 transcripts derived from 30 biopsies. These biopsies were collected in a variety of clinical settings, including normal function, acute rejection, interstitial nephritis, interstitial fibrosis/tubular atrophy and polyomavirus nephropathy. Transcripts with coefficient of variation below the 2nd percentile were designated as HKG, and validated by showing their virtual absence in diseased allograft derived transcriptomic data sets available in the GEO. Pathway analysis indicated a role for these genes in maintenance of cell morphology, pyrimidine metabolism, and intracellular protein signaling. CONCLUSIONS Utilization of these objectively defined HKG data sets will guard against errors resulting from focusing on individual genes like 18S RNA, actin & tubulin, which do not maintain constant expression across the known spectrum of renal allograft pathology.
Collapse
|
16
|
Systematic identification and validation of the reference genes from 60 RNA-Seq libraries in the scallop Mizuhopecten yessoensis. BMC Genomics 2019; 20:288. [PMID: 30975074 PMCID: PMC6460854 DOI: 10.1186/s12864-019-5661-x] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2018] [Accepted: 03/31/2019] [Indexed: 12/15/2022] Open
Abstract
Background Reverse transcription quantitative PCR (RT-qPCR) is widely used for gene expression analysis in various organisms. Its accuracy largely relies on the stability of reference genes, making reference gene selection a vital step in RT-qPCR experiments. However, previous studies in mollusks only focused on the reference genes widely used in vertebrates. Results In this study, we conducted the transcriptome-wide identification of reference genes in the bivalve mollusk Mizuhopecten yessoensis based on 60 transcriptomes covering early development, adult tissues and gonadal development. A total of 964, 1210 and 2097 candidate reference genes were identified, respectively, resulting in a core set of 568 genes. Functional enrichment analysis showed that these genes are significantly overrepresented in Gene Ontology (GO) terms or Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to ribosomes, energy production, etc. Six genes (RS23, EF1A, NDUS4, SELR1, EIF3F, and OLA1) were selected from the candidate genes for RT-qPCR validation, together with 6 commonly used reference genes (ACT, CYTC, HEL, EF1B, GAPDH and RPL16). Stability analyses using geNorm, NormFinder and the comparative delta-Ct method revealed that the new candidate reference genes are more stable than the traditionally used genes, and ACT and CYTC are not recommended under either of the three circumstances. There was a significant correlation between the Ct of RT-qPCR and the log2(TPM) of RNA-Seq data (Ct = − 0.94 log2(TPM) + 29.67, R2 = 0.73), making it easy to estimate the Ct values from transcriptome data prior to RT-qPCR experiments. Conclusion Our study represents the first transcriptome-wide identification of reference genes for early development, adult tissues, and gonadal development in the Yesso scallop and will benefit gene expression studies in other bivalve mollusks. Electronic supplementary material The online version of this article (10.1186/s12864-019-5661-x) contains supplementary material, which is available to authorized users.
Collapse
|
17
|
RNAi-Based Identification of Gene-Specific Nuclear Cofactor Networks Regulating Interleukin-1 Target Genes. Front Immunol 2018; 9:775. [PMID: 29755455 PMCID: PMC5934416 DOI: 10.3389/fimmu.2018.00775] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2017] [Accepted: 03/28/2018] [Indexed: 12/22/2022] Open
Abstract
The potent proinflammatory cytokine interleukin (IL)-1 triggers gene expression through the NF-κB signaling pathway. Here, we investigated the cofactor requirements of strongly regulated IL-1 target genes whose expression is impaired in p65 NF-κB-deficient murine embryonic fibroblasts. By two independent small-hairpin (sh)RNA screens, we examined 170 genes annotated to encode nuclear cofactors for their role in Cxcl2 mRNA expression and identified 22 factors that modulated basal or IL-1-inducible Cxcl2 levels. The functions of 16 of these factors were validated for Cxcl2 and further analyzed for their role in regulation of 10 additional IL-1 target genes by RT-qPCR. These data reveal that each inducible gene has its own (quantitative) requirement of cofactors to maintain basal levels and to respond to IL-1. Twelve factors (Epc1, H2afz, Kdm2b, Kdm6a, Mbd3, Mta2, Phf21a, Ruvbl1, Sin3b, Suv420h1, Taf1, and Ube3a) have not been previously implicated in inflammatory cytokine functions. Bioinformatics analysis indicates that they are components of complex nuclear protein networks that regulate chromatin functions and gene transcription. Collectively, these data suggest that downstream from the essential NF-κB signal each cytokine-inducible target gene has further subtle requirements for individual sets of nuclear cofactors that shape its transcriptional activation profile.
Collapse
|