Variation In Transcript Abundance Research Articles

Distinct mechanisms of PTEN inactivation in dogs and humans highlight convergent molecular events that drive cell division in the pathogenesis of osteosarcoma

A hallmark of osteosarcoma in both human and canine tumors is somatic fragmentation and rearrangement of chromosome structure which leads to recurrent increases and decreases in DNA copy number. The PTEN gene has been implicated as an important tumor suppressor in osteosarcoma via forward genetic screens. Here, we analyzed copy number changes, promoter methylation and transcriptomes to better understand the role of PTEN in canine and human osteosarcoma. Reduction in PTEN copy number was observed in 23 of 95 (25%) of the canine tumors examined leading to corresponding decreases in PTEN transcript levels from RNA-Seq samples. Unexpectedly, canine tumors with an intact PTEN locus had higher levels of PTEN transcripts than human tumors. This variation in transcript abundance was used to evaluate the role of PTEN in osteosarcoma biology. Decreased PTEN copy number and transcript level was observed in - and likely an important driver of - increases in cell cycle transcripts in four independent canine transcriptional datasets. In human osteosarcoma, homozygous copy number loss was not observed, instead increased methylation of the PTEN promoter was associated with increased cell cycle transcripts. Somatic modification of PTEN, either by homozygous deletion in dogs or by promoter methylation in humans, is clinically relevant to osteosarcoma, because the cell cycle related transcripts are associated with patient outcomes. The PTEN gene is part of a syntenic rearrangement unique to the canine genome, making it susceptible to somatic loss of both copies of distal chromosome 26 which also includes the FAS death receptor. Significance StatementPTEN function is abrogated by different mechanisms in canine and human osteosarcoma tumors leading to uncontrolled cell cycling. Somatic loss of this canine specific syntenic region may help explain why the canine genome appears to be uniquely susceptible to osteosarcoma. Syntenic arrangement, in the context of copy number change, may lead to synergistic interactions that in turn modify species specific cancer risk. Comparative models of tumorigenesis may utilize different driver mechanisms.

The role of non-additive gene action on gene expression variation in plant domestication

BackgroundPlant domestication is a remarkable example of rapid phenotypic transformation of polygenic traits, such as organ size. Evidence from a handful of study cases suggests this transformation is due to gene regulatory changes that result in non-additive phenotypes. Employing data from published genetic crosses, we estimated the role of non-additive gene action in the modulation of transcriptional landscapes in three domesticated plants: maize, sunflower, and chili pepper. Using A. thaliana, we assessed the correlation between gene regulatory network (GRN) connectivity properties, transcript abundance variation, and gene action. Finally, we investigated the propagation of non-additive gene action in GRNs.ResultsWe compared crosses between domesticated plants and their wild relatives to a set of control crosses that included a pair of subspecies evolving under natural selection and a set of inbred lines evolving under domestication. We found abundance differences on a higher portion of transcripts in crosses between domesticated-wild plants relative to the control crosses. These transcripts showed non-additive gene action more often in crosses of domesticated-wild plants than in our control crosses. This pattern was strong for genes associated with cell cycle and cell fate determination, which control organ size. We found weak but significant negative correlations between the number of targets of trans-acting genes (Out-degree) and both the magnitude of transcript abundance difference a well as the absolute degree of dominance. Likewise, we found that the number of regulators that control a gene’s expression (In-degree) is weakly but negatively correlated with the magnitude of transcript abundance differences. We observed that dominant-recessive gene action is highly propagable through GRNs. Finally, we found that transgressive gene action is driven by trans-acting regulators showing additive gene action.ConclusionsOur study highlights the role of non-additive gene action on modulating domestication-related traits, such as organ size via regulatory divergence. We propose that GRNs are shaped by regulatory changes at genes with modest connectivity, which reduces the effects of antagonistic pleiotropy. Finally, we provide empirical evidence of the propagation of non-additive gene action in GRNs, which suggests a transcriptional epistatic model for the control of polygenic traits, such as organ size.

Effect of Different Anticoagulant Agents on Immune-Related Genes in Leukocytes Isolated from the Whole-Blood of Holstein Cows

Anticoagulants, such as ethylenediaminetetraacetic acid (EDTA), sodium citrate (Na-citrate), or heparin are normally used in hematological clinical tests to prevent coagulation. Although anticoagulants are fundamental for the correct application of clinical tests, they produce adverse effects in different fields, such as those involving specific molecular techniques; for instance, quantitative real time polymerase chain reactions (qPCR) and gene expression evaluation. For this reason, the aim of this study was to evaluate the expression of 14 genes in leukocytes that were isolated from the blood of Holstein cows, and which were collected in Li-heparin, K-EDTA, or Na-citrate tubes; then, they were analyzed using qPCR. Only the SDHA gene showed a significant dependence (p ≤ 0.05) on the anticoagulant that was used with the lowest expression; this was observed in Na-Citrate after being compared with Li-heparin and K-EDTA (p < 0.05). Although a variation in transcript abundance with the three anticoagulants was observed in almost all the investigated genes, the relative abundance levels were not statistically significant. In conclusion, the qPCR results were not influenced by the presence of the anticoagulant; thus, we had the opportunity to choose the test tube that was used in the experiment without interfering effects impacting the gene expression levels caused by the anticoagulant.

Natural variation in the maternal and zygotic mRNA complements of the early embryo in Drosophila melanogaster

BackgroundMaternal gene products supplied to the egg during oogenesis drive the earliest events of development in all metazoans. After the initial stages of embryogenesis, maternal transcripts are degraded as zygotic transcription is activated; this is known as the maternal to zygotic transition (MZT). Recently, it has been shown that the expression of maternal and zygotic transcripts have evolved in the Drosophila genus over the course of 50 million years. However, the extent of natural variation of maternal and zygotic transcripts within a species has yet to be determined. We asked how the maternal and zygotic pools of mRNA vary within and between populations of D. melanogaster. In order to maximize sampling of genetic diversity, African lines of D. melanogaster originating from Zambia as well as DGRP lines originating from North America were chosen for transcriptomic analysis.ResultsGenerally, we find that maternal transcripts are more highly conserved, and zygotic transcripts evolve at a higher rate. We find that there is more within-population variation in transcript abundance than between populations and that expression variation is highest post- MZT between African lines.ConclusionsDetermining the natural variation of gene expression surrounding the MZT in natural populations of D. melanogaster gives insight into the extent of how a tightly regulated process may vary within a species, the extent of developmental constraint at both stages and on both the maternal and zygotic genomes, and reveals expression changes allowing this species to adapt as it spread across the world.

Genetic regulation and variation of expression of miRNA and mRNA transcripts in fetal muscle tissue in the context of sex, dam and variable fetal weight

BackgroundImpaired skeletal muscle growth in utero can result in reduced birth weight and pathogenesis of intrauterine growth restriction. Fetal and placental growth is influenced by many factors including genetic, epigenetic and environmental factors. In fact, the sex and genotype of the fetus itself, as well as the mother providing it with a suitable environment, influence the growth of the fetus. Hence, our goal was to decipher and elucidate the molecular pathways of developmental processes mediated by miRNAs and mRNAs in fetal muscle tissue in the context of sex, dam, and fetal weight. Therefore, we analyse the variation of miRNA and mRNA expression in relation to these factors. In addition, the coincidence of genetic regulation of these mRNAs and miRNAs, as revealed by expression quantitative trait loci (eQTL) analyses, with sex-, mother- and weight-associated expression was investigated.MethodsA three-generation pig F2 population (n = 118) based on reciprocal crossing of German Landrace (DL) and Pietrain (Pi) was used. Genotype information and transcriptomic data (mRNA and miRNA) from longissimus dorsi muscle (LDM) of pig fetuses sampled at 63 days post-conception (dpc) were used for eQTL analyses.ResultsThe transcript abundances of 13, 853, and 275 probe-sets were influenced by sex, dam and fetal weight at 63 dpc, respectively (FDR < 5%). Most of significant transcripts affected by sex were located on the sex chromosomes including KDM6A and ANOS1 or autosomes including ANKS1B, LOC100155138 and miR-153. The fetal muscle transcripts associated with fetal weight indicated clearer metabolic directions than maternally influenced fetal muscle transcripts. Moreover, coincidence of genetic regulation (eQTL) and variation in transcript abundance due to sex, dam and fetal weight were identified.ConclusionsIntegrating information on eQTL, sex-, dam- and weight-associated differential expression and QTL for fetal weight allowed us to identify molecular pathways and shed light on the basic biological processes associated with differential muscle development in males and females, with implications for adaptive fetal programming.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

BackgroundGenome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”ResultsWe experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.ConclusionWe functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Sources of transcription variation in Plasmodium falciparum

Variation in transcript abundance can contribute to both short-term environmental response and long-term evolutionary adaptation. Most studies are designed to assess differences in mean transcription levels and do not consider other potentially important and confounding sources of transcriptional variation. Detailed quantification of variation sources will improve our ability to detect and identify the mechanisms that contribute to genome-wide transcription changes that underpin adaptive responses. To quantify innate levels of expression variation, we measured mRNA levels for more than 5000 genes in the malaria parasite, Plasmodium falciparum, among clones derived from two parasite strains across biologically and experimentally replicated batches. Using a mixed effects model, we partitioned the total variation among four sources—between strain, within strain, environmental batch effects, and stochastic noise. We found 646 genes with significant variation attributable to at least one of these sources. These genes were categorized by their predominant variation source and further examined using gene ontology enrichment analysis to associate function with each source of variation. Genes with environmental batch effect and within strain transcript variation may contribute to phenotypic plasticity, while genes with between strain variation may contribute to adaptive responses and processes that lead to parasite strain-specific survival under varied conditions.

TU12. COPY NUMBER VARIATION (CNV) IN AUTISM SPECTRUM DISORDERS (ASD) AND CO-OCCURRING EPILEPSY

Recent research has shown expression of clock genes in peripheral tissue explants, targeting multiple pathways leading to the entrainment of circadian rhythms. Temporal variations are not solely regulated by photoperiod, but factors such as maternal feed availability can entrain fetal circadian clock. Currently, a paucity of information exists for clock gene expression and short-term temporal transcript abundance in the bovine placenta, which is essential for proper offspring development. Therefore, the objective of this study was to determine the effect of early to mid-gestational nutrient restriction on clock genes, angiogenic factors, and nutrient sensing genes mRNA transcript abundance in placental explants during a 24 h period. Placentomes from adequately fed and nutrient restricted heifers were collected via Cesarean section at day 180 of gestation; separated into caruncular and cotyledonary tissue and placed in culture media for a 24 h period. The mRNA transcript abundance of clock genes (ARNTL, CRY1, and PER2), angiogenic factors (HIF1A and VEGFA), and nutrient sensing genes (NAMPT and NR3C1) was determined every 4 h. Clock genes were expressed in caruncular and cotyledonary explant tissue. The caruncular explant transcript abundance of the clock genes was not influenced by time (P > 0.05); while ARNTL abundance decreased over time in the cotyledon explant (P < 0.05). A main effect of time was observed for HIF1A, VEGFA, and NR3C1 in the caruncular tissue (P < 0.05). Although, angiogenic factors and nutrient sensing genes in cotyledonary tissue displayed evident temporal variation in transcript abundance (P < 0.05). Nutrient restriction did not alter (P > 0.15) mRNA transcript abundance of clock genes, angiogenic factors, or nutrient sensing genes in either caruncular or cotyledonary tissue. Interestingly, these data may indicate limited transmission and synchronization of maternal and fetal temporal variations in transcript abundance. These findings demonstrate that multiple timepoint collections are needed in future studies due to the innate existence of temporal oscillations observed in the bovine placenta.

EORNA, a barley gene and transcript abundance database

A high-quality, barley gene reference transcript dataset (BaRTv1.0), was used to quantify gene and transcript abundances from 22 RNA-seq experiments, covering 843 separate samples. Using the abundance data we developed a Barley Expression Database (EORNA*) to underpin a visualisation tool that displays comparative gene and transcript abundance data on demand as transcripts per million (TPM) across all samples and all the genes. EORNA provides gene and transcript models for all of the transcripts contained in BaRTV1.0, and these can be conveniently identified through either BaRT or HORVU gene names, or by direct BLAST of query sequences. Browsing the quantification data reveals cultivar, tissue and condition specific gene expression and shows changes in the proportions of individual transcripts that have arisen via alternative splicing. TPM values can be easily extracted to allow users to determine the statistical significance of observed transcript abundance variation among samples or perform meta analyses on multiple RNA-seq experiments. * Eòrna is the Scottish Gaelic word for Barley.

Fragmentation of Small-Cell Lung Cancer Regulatory States in Heterotypic Microenvironments.

Small-cell lung cancers derive from pulmonary neuroendocrine cells, which have stem-like properties to reprogram into other cell types upon lung injury. It is difficult to uncouple transcriptional plasticity of these transformed cells from genetic changes that evolve in primary tumors or secondary metastases. Profiling of single cells is also problematic if the required sample dissociation activates injury-like signaling and reprogramming. Here we defined cell-state heterogeneities in situ through laser capture microdissection-based 10-cell transcriptomics coupled with stochastic-profiling fluctuation analysis. In labeled cells from a small-cell lung cancer mouse model initiated by neuroendocrine deletion of Rb1-Trp53, variations in transcript abundance revealed cell-to-cell differences in regulatory state in vitro and in vivo. Fluctuating transcripts in spheroid culture were partly shared among Rb1-Trp53-null models, and heterogeneities increased considerably when cells were delivered intravenously to colonize the liver. Colonization of immunocompromised animals drove a fractional appearance of alveolar type II-like markers and poised cells for paracrine stimulation from immune cells and hepatocytes. Immunocompetency further exaggerated the fragmentation of tumor states in the liver, yielding mixed stromal signatures evident in bulk sequencing from autochthonous tumors and metastases. Dozens of transcript heterogeneities recurred irrespective of biological context; their mapped orthologs brought together observations of murine and human small-cell lung cancer. Candidate heterogeneities recurrent in the liver also stratified primary human tumors into discrete groups not readily explained by molecular subtype but with prognostic relevance. These data suggest that heterotypic interactions in the liver and lung are an accelerant for intratumor heterogeneity in small-cell lung cancer. SIGNIFICANCE: These findings demonstrate that the single-cell regulatory heterogeneity of small-cell lung cancer becomes increasingly elaborate in the liver, a common metastatic site for the disease.See related articles by Singh and colleagues, p. 1840 and Sutcliffe and colleagues, p. 1868.

Temporal transcript abundance of clock genes, angiogenic factors and nutrient sensing genes in bovine placental explants

Recent research has shown expression of clock genes in peripheral tissue explants, targeting multiple pathways leading to the entrainment of circadian rhythms. Temporal variations are not solely regulated by photoperiod, but factors such as maternal feed availability can entrain fetal circadian clock. Currently, a paucity of information exists for clock gene expression and short-term temporal transcript abundance in the bovine placenta, which is essential for proper offspring development. Therefore, the objective of this study was to determine the effect of early to mid-gestational nutrient restriction on clock genes, angiogenic factors, and nutrient sensing genes mRNA transcript abundance in placental explants during a 24 h period. Placentomes from adequately fed and nutrient restricted heifers were collected via Cesarean section at day 180 of gestation; separated into caruncular and cotyledonary tissue and placed in culture media for a 24 h period. The mRNA transcript abundance of clock genes (ARNTL, CRY1, and PER2), angiogenic factors (HIF1A and VEGFA), and nutrient sensing genes (NAMPT and NR3C1) was determined every 4 h. Clock genes were expressed in caruncular and cotyledonary explant tissue. The caruncular explant transcript abundance of the clock genes was not influenced by time (P > 0.05); while ARNTL abundance decreased over time in the cotyledon explant (P < 0.05). A main effect of time was observed for HIF1A, VEGFA, and NR3C1 in the caruncular tissue (P < 0.05). Although, angiogenic factors and nutrient sensing genes in cotyledonary tissue displayed evident temporal variation in transcript abundance (P < 0.05). Nutrient restriction did not alter (P > 0.15) mRNA transcript abundance of clock genes, angiogenic factors, or nutrient sensing genes in either caruncular or cotyledonary tissue. Interestingly, these data may indicate limited transmission and synchronization of maternal and fetal temporal variations in transcript abundance. These findings demonstrate that multiple timepoint collections are needed in future studies due to the innate existence of temporal oscillations observed in the bovine placenta.

A comparison of methods accounting for batch effects in differential expression analysis of UMI count based single cell RNA sequencing

Accounting for batch effects, especially latent batch effects, in differential expression (DE) analysis is critical for identifying true biological effects. Single-cell RNA sequencing (scRNA-seq) is a powerful tool for quantifying cell-to-cell variation in transcript abundance and characterizing cellular dynamics. Although many scRNA-seq DE analysis methods accommodate known batch variables, their performance has not been systematically evaluated. Moreover, the challenge of accounting for latent batch variables in scRNA-seq DE analysis is largely unmet. In contrast, many methods have been developed to account for batch variables (either known or latent) in other high-dimensional data, especially bulk RNA-seq. We extensively evaluate 11 methods for batch variables in different scRNA-seq DE analysis scenarios, with a primary focus on latent batch variables. We demonstrate that for known batch variables, incorporating them as covariates into a regression model outperformed approaches using a batch-corrected matrix. For latent batches, fixed effects models have inflated FDRs, whereas aggregation-based methods and mixed effects models have significant power loss. Surrogate variable based methods generally control the FDR well while achieving good power with small group effects. However, their performance (except that of SVA) deteriorated substantially in scenarios involving large group effects and/or group label impurity. In these settings, SVA achieves relatively good performance despite an occasionally inflated FDR (up to 0.2). Finally we make the following recommendations for scRNA-seq DE analysis: 1) incorporate known batch variables instead of using batch-corrected data; and 2) employ SVA for latent batch correction. However, better methods are still needed to fully unleash the power of scRNA-seq.

Identification of differentially expressed genes of Haloxylon ammodendronin response to salinity stress

Haloxylon ammodendron (C.A. Mey.), an endangered desert tree with excellent drought and salinity tolerance, provides a unique genotype to characterize and understand the tolerance mechanisms. In this study, four RNA-Seq libraries were constructed and sequenced from H. ammodendron under salinity stress. Total 12 027 differentially expressed genes (DEGs) were identified, in which 4 023, 3 517, 4 487 genes were differentially expressed under light salinity stress (200 mM NaCl), moderate salinity stress (400 mM NaCl), and severe salinity stress (800 mM NaCl), respectively. The up-regulated DEGs included several transcription factors (e.g., MYB and bHLH), hormone-related genes (e.g., cytochrome P450), protein kinases (e.g., Atpk2-Atpk19 like), and genes involved in carbon metabolism (e.g., UDP glycosyltransferase), osmotic regulation (e.g., proline transporter), and ubiquitin proteasome system (e.g., ubiquitin-conjugating enzymes). Heat shock proteins were identified as positive regulators of salinity tolerance in H. ammodendron. The expression patterns of 13 DEGs verified by real-time quantitative PCR were identically consistent with the variations in transcript abundance identified by RNA-Seq. Our results provide new insights into molecular mechanism of H. ammodendron in response to salinity stress.

Faster-X evolution of gene expression is driven by recessive adaptive cis-regulatory variation in Drosophila.

The hemizygosity of the X (Z) chromosome fully exposes the fitness effects of mutations on that chromosome and has evolutionary consequences on the relative rates of evolution of X and autosomes. Specifically, several population genetics models predict increased rates of evolution in X-linked loci relative to autosomal loci. This prediction of faster-X evolution has been evaluated and confirmed for both protein coding sequences and gene expression. In the case of faster-X evolution for gene expression divergence, it is often assumed that variation in 5' noncoding sequences is associated with variation in transcript abundance between species but a formal, genomewide test of this hypothesis is still missing. Here, I use whole genome sequence data in Drosophila yakuba and D.santomea to evaluate this hypothesis and report positive correlations between sequence divergence at 5' noncoding sequences and gene expression divergence. I also examine polymorphism and divergence in 9,279 noncoding sequences located at the 5' end of annotated genes and detected multiple signals of positive selection. Notably, I used the traditional synonymous sites as neutral reference to test for adaptive evolution, but I also used bases 8-30 of introns <65bp, which have been proposed to be a better neutral choice. X-linked genes with high degree of male-biased expression show the most extreme adaptive pattern at 5' noncoding regions, in agreement with faster-X evolution for gene expression divergence and a higher incidence of positively selected recessive mutations.

Single Molecule Fluorescence &lt;em&gt;In Situ&lt;/em&gt; Hybridization (smFISH) Analysis in Budding Yeast Vegetative Growth and Meiosis

Single molecule fluorescence in situ hybridization (smFISH) is a powerful technique to study gene expression in single cells due to its ability to detect and count individual RNA molecules. Complementary to deep sequencing-based methods, smFISH provides information about the cell-to-cell variation in transcript abundance and the subcellular localization of a given RNA. Recently, we have used smFISH to study the expression of the gene NDC80 during meiosis in budding yeast, in which two transcript isoforms exist and the short transcript isoform has its entire sequence shared with the long isoform. To confidently identify each transcript isoform, we optimized known smFISH protocols and obtained high consistency and quality of smFISH data for the samples acquired during budding yeast meiosis. Here, we describe this optimized protocol, the criteria that we use to determine whether high quality of smFISH data is obtained, and some tips for implementing this protocol in other yeast strains and growth conditions.

Single Molecule Fluorescence &lt;em&gt;In Situ&lt;/em&gt; Hybridization (smFISH) Analysis in Budding Yeast Vegetative Growth and Meiosis

Single molecule fluorescence in situ hybridization (smFISH) is a powerful technique to study gene expression in single cells due to its ability to detect and count individual RNA molecules. Complementary to deep sequencing-based methods, smFISH provides information about the cell-to-cell variation in transcript abundance and the subcellular localization of a given RNA. Recently, we have used smFISH to study the expression of the gene NDC80 during meiosis in budding yeast, in which two transcript isoforms exist and the short transcript isoform has its entire sequence shared with the long isoform. To confidently identify each transcript isoform, we optimized known smFISH protocols and obtained high consistency and quality of smFISH data for the samples acquired during budding yeast meiosis. Here, we describe this optimized protocol, the criteria that we use to determine whether high quality of smFISH data is obtained, and some tips for implementing this protocol in other yeast strains and growth conditions.

Multiple Transcript Properties Related to Translation Affect mRNA Degradation Rates in Saccharomyces cerevisiae.

Degradation of mRNA contributes to variation in transcript abundance. Studies of individual mRNAs have shown that both cis and trans factors affect mRNA degradation rates. However, the factors underlying transcriptome-wide variation in mRNA degradation rates are poorly understood. We investigated the contribution of different transcript properties to transcriptome-wide degradation rate variation in the budding yeast, Saccharomyces cerevisiae, using multiple regression analysis. We find that multiple transcript properties are significantly associated with variation in mRNA degradation rates, and that a model incorporating these properties explains ∼50% of the genome-wide variance. Predictors of mRNA degradation rates include transcript length, ribosome density, biased codon usage, and GC content of the third position in codons. To experimentally validate these factors, we studied individual transcripts expressed from identical promoters. We find that decreasing ribosome density by mutating the first translational start site of a transcript increases its degradation rate. Using coding sequence variants of green fluorescent protein (GFP) that differ only at synonymous sites, we show that increased GC content of the third position of codons results in decreased rates of mRNA degradation. Thus, in steady-state conditions, a large fraction of genome-wide variation in mRNA degradation rates is determined by inherent properties of transcripts, many of which are related to translation, rather than specific regulatory mechanisms.

Genome-wide identification and characterization of the NF-Y gene family in grape (vitis vinifera L.).

BackgroundNuclear factor Y (NF-Y) transcription factor is composed of three distinct subunits: NF-YA, NF-YB and NF-YC. Many members of NF-Y family have been reported to be key regulators in plant development, phytohormone signaling and drought tolerance. However, the function of the NF-Y family is less known in grape (Vitis vinifera L.).ResultsA total of 34 grape NF-Y genes that distributed unevenly on grape (V. vinifera) chromosomes were identified in this study. Phylogenetic analysis was performed to predict functional similarities between Arabidopsis thaliana and grape NF-Y genes. Comparison of the structures of grape NF-Y genes (VvNF-Ys) revealed their functional conservation and alteration. Furthermore, we investigated the expression profiles of VvNF-Ys in response to various stresses, phytohormone treatments, and in leaves and grape berries with various sugar contents at different developmental stages. The relationship between VvNF-Y transcript levels and sugar content was examined to select candidates for exogenous sugar treatments. Quantitative real-time PCR (qPCR) indicated that many VvNF-Ys responded to different sugar stimuli with variations in transcript abundance. qPCR and publicly available microarray data suggest that VvNF-Ys exhibit distinct expression patterns in different grape organs and developmental stages, and a number of VvNF-Ys may participate in responses to multiple abiotic and biotic stresses, phytohormone treatments and sugar accumulation or metabolism.ConclusionsIn this study, we characterized 34 VvNF-Ys based on their distributions on chromosomes, gene structures, phylogenetic relationship with Arabidopsis NF-Y genes, and their expression patterns. The potential roles of VvNF-Ys in sugar accumulation or metabolism were also investigated. Altogether, the data provide significant insights on VvNF-Ys, and lay foundations for further functional studies of NF-Y genes in grape.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-2989-3) contains supplementary material, which is available to authorized users.

Haplotype and diplotype analyses of variation in ERCC5 transcription cis-regulation in normal bronchial epithelial cells.

Excision repair cross-complementation group 5 (ERCC5) gene plays an important role in nucleotide excision repair, and dysregulation of ERCC5 is associated with increased lung cancer risk. Haplotype and diplotype analyses were conducted in normal bronchial epithelial cells (NBEC) to better understand mechanisms responsible for interindividual variation in transcript abundance regulation of ERCC5 We determined genotypes at putative ERCC5 cis-regulatory SNPs (cis-rSNP) rs751402 and rs2296147, and marker SNPs rs1047768 and rs17655. ERCC5 allele-specific transcript abundance was assessed by a recently developed targeted sequencing method. Syntenic relationships among alleles at rs751402, rs2296147, and rs1047768 were assessed by allele-specific PCR followed by Sanger sequencing. We then assessed association of ERCC5 allele-specific expression at rs1047768 with haplotype and diplotype structure at cis-rSNPs rs751402 and rs2296147. Genotype analysis revealed significantly (P < 0.005) higher interindividual variation in allelic ratios in cDNA samples relative to matched gDNA samples at both rs1047768 and rs17655. By diplotype analysis, mean expression was higher at the rs1047768 alleles syntenic with rs2296147 T allele compared with rs2296147 C allele. Furthermore, mean expression was lower at rs17655 C allele, which is syntenic with G allele at a linked SNP rs873601 (D' = 0.95). These data support the conclusions that in NBEC, T allele at SNP rs2296147 upregulates ERCC5, variation at rs751402 does not alter ERCC5 regulation, and that C allele at SNP rs17655 downregulates ERCC5 Variation in ERCC5 transcript abundance associated with allelic variation at these SNPs could result in variation in NER function in NBEC and lung cancer risk.

Transcriptional regulation of topology modulators and transcription regulators of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb) is a formidable pathogen which has the ability to survive the hostile environment of the host by evading the host defense system. The re-configuration of its transcriptional and metabolic process allows the pathogen to confront the adverse environment within the host macrophages. The factors that assist the transcription and modulate the DNA topology would have to play a key role in the regulation of global gene expression of the organism. How transcription of these essential housekeeping genes alters in response to growth conditions and environmental stress has not been addressed together in a set of experimental conditions in Mtb. Now, we have mapped the transcription start sites (TSS) and promoters of several genes that play a central role in the regulation of DNA topology and transcription in Mtb. Using in vivo reporter assays, we validated the activity of the identified promoter elements in different growth conditions. The variation in transcript abundance of these essential genes was also analyzed in growth phase-dependent manner. These data provide the first glimpse into the specific adaptive changes in the expression of genes involved in transcription and DNA topology modulation in Mtb.