Screening the optimal housekeeping genes (HKGs) of placenta tissues by RNA-sequence and qRT-PCR throughout gestation in goat (Capra Hircus).

Selection of stable housekeeping genes (HKGs) is very important for accurate calculation of relative expression levels of target genes by quantitative real-time polymerase chain reaction (qRT-PCR). At present, the appropriate HKGs have not been identified in placental tissues throughout the pregnancy of the goat. In our study, 20 HKGs were tentatively selected from RNA-seq data and previous reports. The cycle threshold (Ct) of HKGs was determined by qRT-PCR in trophoblast membrane and cotyledon villus collected from 38 Dazu Black goats on gestation days of 20, 25, 30, 45, 60, 90, 120, and 150 (birth). The expression stability of the HKGs was analyzed by geNorm, Normfinder, Bestkeeper and Delta Ct algorithms, and comprehensively evaluated by ReFinder and ComprFinder. In addition, the optimal HKGs were further verified by placenta-specific genes (SPP1, VEGFA and PAG6). The 16 candidate HKGs (except POP4, TBP, RNF10, UBC) showed a qualified Ct value, less than 28. Among them, YWHAZ, EIF3K and PPIB showed the most stable expression in placental tissues during early, mid-late pregnancy and postpartum, but the least stable expression was B2M at early and mid-late stage, and PPIB at postpartum. After comprehensive analysis, RPLP0, EIF3K and YWHAZ were found to be the most stable placental HKGs throughout pregnancy. The classical HKGs, ACTB, GAPDH and 18S RNA have unstable expressions and even ranked at the bottom of the list from comprehensive index, suggesting an inappropriate for target gene normalization. Taken together, our study confirmed that YWHAZ, EIF3K, HMBS and RPLP0 may be the optimal HKGs in goat placenta at different stage of pregnancy, which provided a valuable reference of HKGs on functional gene expression detection for further research on placenta development and growth in ruminants.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors.

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

Systematic and scalable genome-wide essentiality mapping to identify nonessential genes in phages.

Phages are one of the key ecological drivers of microbial community dynamics, function, and evolution. Despite their importance in bacterial ecology and evolutionary processes, phage genes are poorly characterized, hampering their usage in a variety of biotechnological applications. Methods to characterize such genes, even those critical to the phage life cycle, are labor intensive and are generally phage specific. Here, we develop a systematic gene essentiality mapping method scalable to new phage-host combinations that facilitate the identification of nonessential genes. As a proof of concept, we use an arrayed genome-wide CRISPR interference (CRISPRi) assay to map gene essentiality landscape in the canonical coliphages λ and P1. Results from a single panel of CRISPRi probes largely recapitulate the essential gene roster determined from decades of genetic analysis for lambda and provide new insights into essential and nonessential loci in P1. We present evidence of how CRISPRi polarity can lead to false positive gene essentiality assignments and recommend caution towards interpreting CRISPRi data on gene essentiality when applied to less studied phages. Finally, we show that we can engineer phages by inserting DNA barcodes into newly identified inessential regions, which will empower processes of identification, quantification, and tracking of phages in diverse applications.

Comparative genomics hints at dispensability of multiple essential genes in two Escherichia coli L-form strains.

Despite the critical role of bacterial cell walls in maintaining cell shapes, certain environmental stressors can induce the transition of many bacterial species into a wall-deficient state called L-form. Long-term induced Escherichia coli L-forms lose their rod shape and usually hold significant mutations that affect cell division and growth. Besides this, the genetic background of L-form bacteria is still poorly understood. In the present study, the genomes of two stable L-form strains of E. coli (NC-7 and LWF+) were sequenced and their gene mutation status was determined and compared with their parental strains. Comparative genomic analysis between two L-forms reveals both unique adaptions and common mutated genes, many of which belong to essential gene categories not involved in cell wall biosynthesis, indicating that L-form genetic adaptation impacts crucial metabolic pathways. Missense variants from L-forms and Lenski's long-term evolution experiment (LTEE) were analyzed in parallel using an optimized DeepSequence pipeline to investigate predicted mutation effects (α) on protein functions. We report that the two L-form strains analyzed display a frequency of 6-10% (0% for LTEE) in mutated essential genes where the missense variants have substantial impact on protein functions (α<0.5). This indicates the emergence of different survival strategies in L-forms through changes in essential genes during adaptions to cell wall deficiency. Collectively, our results shed light on the detailed genetic background of two E. coli L-forms and pave the way for further investigations of the gene functions in L-form bacterial models.

Regulatory architecture of cell identity genes and housekeeping genes.

Gene regulation and chromosome architecture are intimately linked. Genes with prominent roles in cell identity are often regulated by clusters of enhancer elements. By contrast, a recent study shows housekeeping genes are often regulated through clustering of promoters. We discuss here new regulatory insights for these two types of genes.

YWHAZ and TBP are potential reference gene candidates for qPCR analysis of response to radiation therapy in colorectal cancer.

The expression profiles of conventional reference genes (RGs), including ACTB and GAPDH, used in quantitative real-time PCR (qPCR), vary depending on tissue types and environmental conditions. We searched for suitable RGs for qPCR to determine the response to radiotherapy in colorectal cancer (CRC) cell lines, organoids, and patient-derived tissues. Ten CRC cell lines (Caco-2, COLO 205, DLD-1, HCT116, HCT-15, HT-29, RKO, SW1116, SW480, and SW620) and organoids were selected and irradiated with 2, 10 or 21 grays (Gy) based on the previous related studies conducted over the last decade. The expression stability of 14 housekeeping genes (HKGs; ACTB, B2M, G6PD, GAPDH, GUSB, HMBS, HPRT1, IPO8, PGK1, PPIA, TBP, TFRC, UBC, and YWHAZ) after irradiation was evaluated using RefFinder using raw quantification cycle (Cq) values obtained from samples before and after irradiation. The expression stability of HKGs were also evaluated for paired fresh frozen tissues or formalin-fixed, paraffin-embedded samples obtained from CRC patients before and after chemoradiotherapy. The expression of YWHAZ and TBP encoding 14-3-3-zeta protein and TATA-binding protein were more stable than the other 12 HKGs in CRC cell lines, organoids, and patient-derived tissues after irradiation. The findings suggest that YWHAZ and TBP are potential RG candidates for normalizing qPCR results in CRC radiotherapy experiments.

Regulatory architecture of housekeeping genes is driven by promoter assemblies.

Genes that are key to cell identity are generally regulated by cell-type-specific enhancer elements bound by transcription factors, some of which facilitate looping to distant gene promoters. In contrast, genes that encode housekeeping functions, whose regulation is essential for normal cell metabolism and growth, generally lack interactions with distal enhancers. We find that Ronin (Thap11) assembles multiple promoters of housekeeping and metabolic genes to regulate gene expression. This behavior is analogous to how enhancers are brought together with promoters to regulate cell identity genes. Thus, Ronin-dependent promoter assemblies provide a mechanism to explain why housekeeping genes can forgo distal enhancer elements and why Ronin is important for cellular metabolism and growth control. We propose that clustering of regulatory elements is a mechanism common to cell identity and housekeeping genes but is accomplished by different factors binding distinct control elements to establish enhancer-promoter or promoter-promoter interactions, respectively.

Identifying essential genes in genome-scale metabolic models of consensus molecular subtypes of colorectal cancer.

Identifying essential targets in the genome-scale metabolic networks of cancer cells is a time-consuming process. The present study proposed a fuzzy hierarchical optimization framework for identifying essential genes, metabolites and reactions. On the basis of four objectives, the present study developed a framework for identifying essential targets that lead to cancer cell death and evaluating metabolic flux perturbations in normal cells that have been caused by cancer treatment. Through fuzzy set theory, a multiobjective optimization problem was converted into a trilevel maximizing decision-making (MDM) problem. We applied nested hybrid differential evolution to solve the trilevel MDM problem to identify essential targets in genome-scale metabolic models for five consensus molecular subtypes (CMSs) of colorectal cancer. We used various media to identify essential targets for each CMS and discovered that most targets affected all five CMSs and that some genes were CMS-specific. We obtained experimental data on the lethality of cancer cell lines from the DepMap database to validate the identified essential genes. The results reveal that most of the identified essential genes were compatible with the colorectal cancer cell lines obtained from DepMap and that these genes, with the exception of EBP, LSS, and SLC7A6, could generate a high level of cell death when knocked out. The identified essential genes were mostly involved in cholesterol biosynthesis, nucleotide metabolisms, and the glycerophospholipid biosynthetic pathway. The genes involved in the cholesterol biosynthetic pathway were also revealed to be determinable, if a cholesterol uptake reaction was not induced when the cells were in the culture medium. However, the genes involved in the cholesterol biosynthetic pathway became non-essential if such a reaction was induced. Furthermore, the essential gene CRLS1 was revealed as a medium-independent target for all CMSs.

Using RNA-seq to identify suitable housekeeping genes for hypoxia studies in human adipose-derived stem cells.

BACKGROUND: Hypoxic culture conditions have been used to study the impact of oxygen deprivation has on gene expression in a number of disease models. However, hypoxia response elements present in the promoter regions of some commonly used housekeeping genes, such as GAPDH and PGK1, can confound the relative gene expression analysis. Thus, there is ongoing debate as to which housekeeping gene is appropriate for studies investigating hypoxia-induced cell responses. Specifically, there is still contradicting information for which housekeeping genes are stable in hypoxia cultures of mesenchymal stem cells. In this study, candidate housekeeping genes curated from the literature were matched to RNAseq data of normoxic and hypoxic human adipose-derived stem cell cultures to determine if gene expression was modulated by hypoxia or not. Expression levels of selected candidates were used to calculate coefficient of variation. Then, accounting for the mean coefficient of variation, and normalised log twofold change, genes were ranked and shortlisted, before validating with qRT-PCR. Housekeeping gene suitability were then determined using GeNorm, NormFinder, BestKeeper, comparative[Formula: see text], RefFinder, and the Livak method. RESULTS: Gene expression levels of 78 candidate genes identified in the literature were analysed in the RNAseq dataset generated from hADSC cultured under Nx and Hx conditions. From the dataset, 15 candidates with coefficient of variation ≤ 0.15 were identified, where differential expression analysis results further shortlisted 8 genes with least variation in expression levels. The top 4 housekeeping gene candidates, ALAS1, RRP1, GUSB, and POLR2B, were chosen for qRT-PCR validation. Additionally, 18S, a ribosomal RNA commonly used as housekeeping gene but not detected in the RNAseq method, was added to the list of housekeeping gene candidates to validate. From qRT-PCR results, 18S and RRP1 were determined to be stably expressed in cells cultured under hypoxic conditions. CONCLUSIONS: We have demonstrated that 18S and RRP1 are suitable housekeeping genes for use in hypoxia studies with human adipose-derived stem cell and should be used in combination. Additionally, these data shown that the commonly used GAPDH and PGK1 are not suitable housekeeping genes for investigations into the effect of hypoxia in human adipose-derived stem cell.

Essential genes: a cross-species perspective.

Protein coding genes exhibit different degrees of intolerance to loss-of-function variation. The most intolerant genes, whose function is essential for cell or/and organism survival, inform on fundamental biological processes related to cell proliferation and organism development and provide a window on the molecular mechanisms of human disease. Here we present a brief overview of the resources and knowledge gathered around gene essentiality, from cancer cell lines to model organisms to human development. We outline the implications of using different sources of evidence and definitions to determine which genes are essential and highlight how information on the essentiality status of a gene can inform novel disease gene discovery and therapeutic target identification.

Identification of essential genes in Coxiella burnetii.

Coxiella burnetii is an intracellular pathogen responsible for causing Q fever in humans, a disease with varied presentations ranging from a mild flu-like sickness to a debilitating illness that can result in endocarditis. The intracellular lifestyle of C. burnetii is unique, residing in an acidic phagolysosome-like compartment within host cells. An understanding of the core molecular biology of C. burnetii will greatly increase our understanding of C. burnetii growth, survival and pathogenesis. We used transposon-directed insertion site sequencing (TraDIS) to reveal C. burnetii Nine Mile Phase II genes fundamental for growth and in vitro survival. Screening a transposon library containing >10 000 unique transposon mutants revealed 512 predicted essential genes. Essential routes of synthesis were identified for the mevalonate pathway, as well as peptidoglycan and biotin synthesis. Some essential genes identified (e.g. predicted type IV secretion system effector genes) are typically considered to be associated with C. burnetii virulence, a caveat concerning the axenic media used in the study. Investigation into the conservation of the essential genes identified revealed that 78 % are conserved across all C. burnetii strains sequenced to date, which probably play critical functions. This is the first report of a whole genome transposon screen in C. burnetii that has been undertaken for the identification of essential genes.

Evaluation of housekeeping genes as references for quantitative real-time PCR analysis of European eel, Anguilla anguilla.

Eels are important aquaculture species for which an increasing number of reference genes are being identified and applied. In this study, five housekeeping genes [RPL7 (ribosomal protein L7), 18 S (18 S ribosomal RNA), EF1A (elongation factor 1α), ACTB (ß-actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase)] were chosen to evaluate their reliability as reference genes for quantitative real-time PCR (qPCR) for the study of Anguilla anguilla. The expression of the selected genes in different eel tissues was determined using qPCR at different growth stages or upon challenge by Anguillid herpesvirus (AngHV), and the expression levels of these genes were then compared and evaluated using the geNorm and NormFinder algorithms. Then, RefFinder was used to comprehensively rank the examined housekeeping genes. Interestingly, the expression of the evaluated housekeeping genes exhibited tissue-dependent and treatment-dependent variations. In different growth periods A. anguilla tissues, the most stable genes were the following: ACTB in mucus; 18 S in skin and kidney; RPL7 in muscle, gill, intestine and brain; EF1A in heart and liver; and GAPDH in spleen. In contrast, in AngHV-challenged A. anguilla tissues, the most stable genes were the following: 18 S in mucus; RPL7 in skin, gill, heart, spleen, kidney and intestine; EF1A in muscle and liver; and ACTB in brain. Further comparison analysis indicated that the expression of RPL7 and EF1A was stable in multiple A. anguilla tissues in different growth periods and in eels challenged by AngHV. Nonetheless, the expression level of GAPDH in eel tissues was lower, and it was unstable in several tissues. These results indicated that the selection of reference genes for qPCR analysis in A. anguilla should be made in accordance with experimental parameters, and both RPL7 and EF1A could be used as reference genes for qPCR study of A. anguilla at different growth stages or upon challenge by AngHV. The reference genes identified in this study could improve the accuracy of qPCR data and facilitate further studies aimed at understanding the biology of eels.

A bacterial pan-genome makes gene essentiality strain-dependent and evolvable.

Many bacterial species are represented by a pan-genome, whose genetic repertoire far outstrips that of any single bacterial genome. Here we investigate how a bacterial pan-genome might influence gene essentiality and whether essential genes that are initially critical for the survival of an organism can evolve to become non-essential. By using Transposon insertion sequencing (Tn-seq), whole-genome sequencing and RNA-seq on a set of 36 clinical Streptococcus pneumoniae strains representative of >68% of the species' pan-genome, we identify a species-wide 'essentialome' that can be subdivided into universal, core strain-specific and accessory essential genes. By employing 'forced-evolution experiments', we show that specific genetic changes allow bacteria to bypass essentiality. Moreover, by untangling several genetic mechanisms, we show that gene essentiality can be highly influenced by and/or be dependent on: (1) the composition of the accessory genome, (2) the accumulation of toxic intermediates, (3) functional redundancy, (4) efficient recycling of critical metabolites and (5) pathway rewiring. While this functional characterization underscores the evolvability potential of many essential genes, we also show that genes with differential essentiality remain important antimicrobial drug target candidates, as their inactivation almost always has a severe fitness cost in vivo.

Naturally occurring Neisseria gonorrhoeae can have large deletions in housekeeping gene abcZ, making them untypable with multilocus sequence typing.

The abcZ gene is an essential housekeeping gene in all the Neisseria species. It is one of the seven genes used for multilocus sequence typing (MLST) this genus. It encodes the cytosolic component of an ATP-binding cassette (ABC) transporter complex of unknown function. We report here the finding of a strain of Neisseria gonorrhoeae with a 485 base pair deletion in the 5' region of the abcZ gene that truncates the protein product from 636 amino acids to 89 amino acids. A second open reading frame (ORF), encoding the latter 388 amino acids of the abcZ gene, was predicted downstream. The deletion will affect MLST profiling; interrogation of genomic sequences from PubMLST revealed that this isolate is not an anomaly. Deletions in abcZ were identified in 256 Neisseria genomes, roughly 0.6% of isolates. Furthermore, these deletions could leave the abcZ gene in a pseudogenized state. Our strain, isolated from a patient with symptoms of gonorrheal infection, nevertheless behaved normal in terms of growth and in vitro phenotypic properties.

Gene essentiality and variability: What is the link? A within- and between-species perspective.

A gene is considered essential when its loss of function results in a deleterious phenotype, hugely reducing the organism's viability or fitness. However, the link between the essentiality of a gene and its degree of polymorphism is unclear. In this review, we show that there is a place for a certain degree of variability, even for essential genes. We first study the role of infectious diseases in the prevalence of genetic disorders among humans: balancing selection has selected harmful variants due to the selective pressure of pathogens because the heterozygous carrier can resist them. Then we show that the environment can induce adaptation of species by rapidly evolving genes. We also study the role of positive selection on speciation, particularly upon genes of the immune, reproductive and nervous systems. Finally, we highlight the role of regulatory sequences in changes in morphology between species and adaptation to the environment within species.

Exploring conditional gene essentiality through systems genetics approaches in yeast.

An essential gene encodes for a cellular function that is required for viability. Although viability is a straightforward phenotype to analyze in yeast, defining a gene as essential is not always trivial. Gene essentiality has generally been studied in specific laboratory strains and under standard growth conditions, however, essentiality can vary across species, strains, and environments. Recent systematic studies of gene essentiality revealed that two sets of essential genes exist: core essential genes that are always required for viability and conditional essential genes that vary in essentiality in different genetic and environmental contexts. Here, we review recent advances made in the systematic analysis of gene essentiality in yeast and discuss the properties that distinguish core from context-dependent essential genes.

What are housekeeping genes?

The concept of "housekeeping gene" has been used for four decades but remains loosely defined. Housekeeping genes are commonly described as "essential for cellular existence regardless of their specific function in the tissue or organism", and "stably expressed irrespective of tissue type, developmental stage, cell cycle state, or external signal". However, experimental support for the tenet that gene essentiality is linked to stable expression across cell types, conditions, and organisms has been limited. Here we use genome-scale functional genomic screens together with bulk and single-cell sequencing technologies to test this link and optimize a quantitative and experimentally validated definition of housekeeping gene. Using the optimized definition, we identify, characterize, and provide as resources, housekeeping gene lists extracted from several human datasets, and 10 other animal species that include primates, chicken, and C. elegans. We find that stably expressed genes are not necessarily essential, and that the individual genes that are essential and stably expressed can considerably differ across organisms; yet the pathways enriched among these genes are conserved. Further, the level of conservation of housekeeping genes across the analyzed organisms captures their taxonomic groups, showing evolutionary relevance for our definition. Therefore, we present a quantitative and experimentally supported definition of housekeeping genes that can contribute to better understanding of their unique biological and evolutionary characteristics.

Evidence that conserved essential genes are enriched for pro-longevity factors.

At the cellular level, many aspects of aging are conserved across species. This has been demonstrated by numerous studies in simple model organisms like Saccharomyces cerevisiae, Caenorhabdits elegans, and Drosophila melanogaster. Because most genetic screens examine loss of function mutations or decreased expression of genes through reverse genetics, essential genes have often been overlooked as potential modulators of the aging process. By taking the approach of increasing the expression level of a subset of conserved essential genes, we found that 21% of these genes resulted in increased replicative lifespan in S. cerevisiae. This is greater than the ~ 3.5% of genes found to affect lifespan upon deletion, suggesting that activation of essential genes may have a relatively disproportionate effect on increasing lifespan. The results of our experiments demonstrate that essential gene overexpression is a rich, relatively unexplored means of increasing eukaryotic lifespan.

A novel system for the generation of baculoviruses mutant for an essential gene.

BACKGROUND: Currently, to delete an essential gene from a baculovirus genome, a cell line stably expressing the gene to be knocked-out should be first generated, which is time-consuming. Alternatively, essential genes can be deleted in E. coli using the λ Red recombination system, which requires an electroporation system. Here, based on homologous recombination in insect cells, we develop an alternative efficient system that requires neither generation of a cell line nor an electroporation system. METHODS AND RESULTS: Using puc19-based inverse PCR, a transfer vector for deleting BmNPV orf92 (Bm92, an essential gene) was efficiently constructed. A copy of Bm92 was introduced into the polyhedrin locus of BmNPV bacmid. The transfer vector was then co-transfected into BmN cell with the modified bacmid to enable homologous recombination at the Bm92 locus. An agarose-free approach was developed for the purification of Bm92-disrupted bacmid viruses in insect cells. Subsequently, BmN cells were co-infected with purified Bm92-disrupted bacmid viruses and unmodified bacmid viruses to allow recombination at the Tn7 insertion site between the two viruses. Finally, bacmid DNA extracted from BmN cells was transformed into chemically-treated competent DH10B cells, and blue colonies containing Bm92-disrupted bacmid were selected using PCR. CONCLUSIONS: For its efficiency and convenience, the system has great potential to be used for the generation of baculovirus knockout mutants.

Transposon insertional mutagenesis of diverse yeast strains suggests coordinated gene essentiality polymorphisms.

Due to epistasis, the same mutation can have drastically different phenotypic consequences in different individuals. This phenomenon is pertinent to precision medicine as well as antimicrobial drug development, but its general characteristics are largely unknown. We approach this question by genome-wide assessment of gene essentiality polymorphism in 16 Saccharomyces cerevisiae strains using transposon insertional mutagenesis. Essentiality polymorphism is observed for 9.8% of genes, most of which have had repeated essentiality switches in evolution. Genes exhibiting essentiality polymorphism lean toward having intermediate numbers of genetic and protein interactions. Gene essentiality changes tend to occur concordantly among components of the same protein complex or metabolic pathway and among a group of over 100 mitochondrial proteins, revealing molecular machines or functional modules as units of gene essentiality variation. Most essential genes tolerate transposon insertions consistently among strains in one or more coding segments, delineating nonessential regions within essential genes.