Origin of phenotypes: genes and transcripts

Deciphering the nexus between long non-coding RNAs and endoplasmic reticulum stress in hepatocellular carcinoma: biomarker discovery and therapeutic horizons

Hepatocellular carcinoma (HCC) remains a significant global health challenge with few effective treatment options. The dysregulation of endoplasmic reticulum (ER) stress responses has emerged as a pivotal factor in HCC progression and therapy resistance. Long non-coding RNAs (lncRNAs) play a crucial role as key epigenetic modifiers in this process. Recent research has explored how lncRNAs influence ER stress which in turn affects lncRNAs activity in HCC. We systematically analyze the current literature to highlight the regulatory roles of lncRNAs in modulating ER stress and vice versa in HCC. Our scrutinization highlights how dysregulated lncRNAs contribute to various facets of HCC, including apoptosis resistance, enhanced proliferation, invasion, and metastasis, all driven by ER stress. Moreover, we delve into the emerging paradigm of the lncRNA-miRNA-mRNA axis, elucidating it as the promising avenue for developing novel biomarkers and paving the way for more personalized treatment options in HCC. Nevertheless, we acknowledge the challenges and future directions in translating these insights into clinical practice. In conclusion, our review provides insights into the complex regulatory mechanisms governing ER stress modulation by lncRNAs in HCC.

Advances in big data and omics: Paving the way for discovery in childhood epilepsies

The insights gained from big data and omics approaches have transformed the field of childhood genetic epilepsy. With an increasing number of individuals receiving genetic testing for seizures, we are provided with an opportunity to identify clinically relevant subgroups and extract meaningful observations from this large-scale clinical data. However, the volume of data from electronic medical records and omics (e.g., genomics, transcriptomics) is so vast that standardized methods, such as the Human Phenotype Ontology, are necessary for reliable and comprehensive characterization. Here, we explore the integration of clinical and omics data, highlighting how these approaches pave the way for discovery in childhood epilepsies.

Adaptive Selection of Cis-regulatory Elements in the Han Chinese

Cis-regulatory elements have an important role in human adaptation to the living environment. However, the lag in population genomic cohort studies and epigenomic studies, hinders the research in the adaptive analysis of cis-regulatory elements in human populations. In this study, we collected 4,013 unrelated individuals and performed a comprehensive analysis of adaptive selection of genome-wide cis-regulatory elements in the Han Chinese. In total, 12.34% of genomic regions are under the influence of adaptive selection, where 1.00% of enhancers and 2.06% of promoters are under positive selection, and 0.06% of enhancers and 0.02% of promoters are under balancing selection. Gene ontology enrichment analysis of these cis-regulatory elements under adaptive selection reveals that many positive selections in the Han Chinese occur in pathways involved in cell-cell adhesion processes, and many balancing selections are related to immune processes. Two classes of adaptive cis-regulatory elements related to cell adhesion were in-depth analyzed, one is the adaptive enhancers derived from neanderthal introgression, leads to lower hyaluronidase level in skin, and brings better performance on UV-radiation resistance to the Han Chinese. Another one is the cis-regulatory elements regulating wound healing, and the results suggest the positive selection inhibits coagulation and promotes angiogenesis and wound healing in the Han Chinese. Finally, we found that many pathogenic alleles, such as risky alleles of type 2 diabetes or schizophrenia, remain in the population due to the hitchhiking effect of positive selections. Our findings will help deepen our understanding of the adaptive evolution of genome regulation in the Han Chinese.

Identification of polymorphic loci in OSMR and GHR genes and analysis of their association with growth traits in sheep

The aim of this study was to analyze the effect of OSMR and GHR genes polymorphisms on growth traits in sheep. The single nucleotide polymorphisms of sheep OSMR and GHR genes were identified by DNA sequencing technology. A total of two intronic mutations g.2443 T > C and g.170196 A > G were identified in OSMR and GHR, respectively. Correlation analysis was carried out between the obtained genotypes and the growth traits of sheep. The results showed that at the OSMR g.2443 T > C locus, the body weight, chest circumference and cannon circumference of the TT genotype sheep were significantly higher than those of the CC genotype sheep (p < .05). At the GHR g.170196 A > G locus, the body weight, body length, chest circumference and cannon circumference of the AA genotype sheep were significantly higher than those of the AG genotype and GG genotype sheep (p < .05). Moreover, the body weight of sheep of combination TTOSMR/AAGHR genotype was significantly higher than that of other combination genotypes (p < .05). Therefore, we believe that the polymorphic sites identified in the OSMR and GHR genes can be used as candidate molecular markers for breeding good traits in sheep.

Long non-coding RNAs: definitions, functions, challenges and recommendations

Genes specifying long non-coding RNAs (lncRNAs) occupy a large fraction of the genomes of complex organisms. The term 'lncRNAs' encompasses RNA polymerase I (Pol I), Pol II and Pol III transcribed RNAs, and RNAs from processed introns. The various functions of lncRNAs and their many isoforms and interleaved relationships with other genes make lncRNA classification and annotation difficult. Most lncRNAs evolve more rapidly than protein-coding sequences, are cell type specific and regulate many aspects of cell differentiation and development and other physiological processes. Many lncRNAs associate with chromatin-modifying complexes, are transcribed from enhancers and nucleate phase separation of nuclear condensates and domains, indicating an intimate link between lncRNA expression and the spatial control of gene expression during development. lncRNAs also have important roles in the cytoplasm and beyond, including in the regulation of translation, metabolism and signalling. lncRNAs often have a modular structure and are rich in repeats, which are increasingly being shown to be relevant to their function. In this Consensus Statement, we address the definition and nomenclature of lncRNAs and their conservation, expression, phenotypic visibility, structure and functions. We also discuss research challenges and provide recommendations to advance the understanding of the roles of lncRNAs in development, cell biology and disease.

Identification of TRAPPC9 and BAIAP2 Gene Polymorphisms and Their Association With Fat Deposition-Related Traits in Hu Sheep

Fat deposition is an important economic trait that is closely related to feed efficiency and carcass performance in livestock. In this study, the fat deposition-related traits of 1,293 Hu sheep were measured and descriptive statistical analysis was conducted. The results showed that the coefficient of variation of all fat deposition-related traits was higher than 24%. In addition, single nucleotide polymorphisms and the expression characteristics of TRAPPC9 (encoding trafficking protein particle complex subunit 9) and BAIAP2 (encoding brain-specific Angiogenesis inhibitor 1-associated protein 2) genes in Hu sheep were detected using PCR amplification, Sanger sequencing, KASPar genotyping, and quantitative real-time reverse transcription PCR (qRT-PCR). The associations between SNPs and fat deposition-related traits were also analyzed. Two intronic mutations, TRAPPC9 g.57654 A > G and BAIAP2 g.46061 C > T, were identified in Hu sheep. The result of association analysis showed that TRAPPC9 g.57654 A > G and BAIAP2 g.46061 C > T were both significantly associated with the weight of tail fat, tail fat relative weight (body weight), and tail fat relative weight (carcass) (P < 0.05). Comprehensive effects analysis showed that there were significant differences between the combined genotypes and tail fat and perirenal fat deposition. Moreover, qRT-PCR analysis showed that TRAPPC9 and BAIAP2 are widely expressed, and their expression levels were significantly higher in the small-tail group compared with those in the big-tail group (P < 0.01). These results provided important candidate molecular markers that could be used in strategies to reduce tail fat deposition in Hu sheep.

Wheat glutamine synthetase TaGSr-4B is a candidate gene for a QTL of thousand grain weight on chromosome 4B

KEY MESSAGE

Glutamine synthetase TaGSr-4B is a candidate gene for a QTL of thousand grain weight on 4B, and the gene marker is ready for wheat breeding. A QTL for thousand grain weight (TGW) in wheat was previously mapped on chromosome 4B in a DH population of Westonia × Kauz. For identifying the candidate genes of the QTL, wheat 90 K SNP array was used to saturate the existing linkage map, and four field trials plus one glasshouse experiment over five locations were conducted to refine the QTL. Three nitrogen levels were applied to two of those field trials, resulting in a TGW phenotype data set from nine environments. A robust TGW QTL cluster including 773 genes was detected in six environments with the highest LOD value of 13.4. Based on differentiate gene expression within the QTL cluster in an RNAseq data of Westonia and Kauz during grain filling, a glutamine synthesis gene (GS: TaGSr-4B) was selected as a potential candidate gene for the QTL. A SNP on the promoter region between Westonia and Kauz was used to develop a cleaved amplified polymorphic marker for TaGSr-4B gene mapping and QTL reanalysing. As results, TGW QTL appeared in seven environments, and in four out of seven environments, the TGW QTL were localized on the TaGSr-4B locus and showed significant contributions to the phenotype. Based on the marker, two allele groups of Westonia and Kauz formed showed significant differences on TGW in eight environments. In agreement with the roles of GS genes on nitrogen and carbon remobilizations, TaGSr-4B is likely the candidate gene of the TGW QTL on 4B and the TaGSr-4B gene marker is ready for wheat breeding.

Extent and complexity of RNA processing in honey bee queen and worker caste development

The distinct honeybee (Apis mellifera) worker and queen castes have become a model for the study of genomic mechanisms of phenotypic plasticity. Here we performed a nanopore-based direct RNA sequencing with exceptionally long reads to compare the mRNA transcripts between queen and workers at three points during their larval development. We found thousands of significantly differentially expressed transcript isoforms (DEIs) between queen and worker larvae. These DEIs were formatted by a flexible splicing system. We showed that poly(A) tails participated in this caste differentiation by negatively regulating the expression of DEIs. Hundreds of isoforms uniquely expressed in either queens or workers during their larval development, and isoforms were expressed at different points in queen and worker larval development demonstrating a dynamic relationship between isoform expression and developmental mechanisms. These findings show the full complexity of RNA processing and transcript expression in honey bee phenotypic plasticity.

•

Honeybee caste differentiation has a complexity of RNA processing

•

Isoforms differentially express between queens and workers during larval development

•

Isoforms are formatted by a flexible alternative splicing system

•

Poly(A) tails are negatively correlated with isoform expression

MicroRNAs and their delivery in diabetic fibrosis

The global prevalence of diabetes mellitus was estimated to be 463 million people in 2019 and is predicted to rise to 700 million by 2045. The associated financial and societal costs of this burgeoning epidemic demand an understanding of the pathology of this disease, and its complications, that will inform treatment to enable improved patient outcomes. Nearly two decades after the sequencing of the human genome, the significance of noncoding RNA expression is still being assessed. The family of functional noncoding RNAs known as microRNAs regulates the expression of most genes encoded by the human genome. Altered microRNA expression profiles have been observed both in diabetes and in diabetic complications. These transcripts therefore have significant potential and novelty as targets for therapy, therapeutic agents and biomarkers.

The PI3K/Akt Pathway: Emerging Roles in Skin Homeostasis and a Group of Non-Malignant Skin Disorders

The phosphatidylinositol 3-kinase (PI3K)/protein kinase B (Akt) signaling pathway regulates cell proliferation, differentiation, and migration, along with angiogenesis and metabolism. Additionally, it could mediate skin development and homeostasis. There is much evidence to suggest that dysregulation of PI3K/Akt pathway is frequently associated with several human cutaneous malignancies like malignant melanoma (MM), basal cell carcinoma (BCC), and cutaneous squamous cell carcinoma (SCC), as well as their poor outcomes. Nevertheless, emerging roles of PI3K/Akt pathway cascade in a group of common non-malignant skin disorders including acne and psoriasis, among others, have been recognized. The enhanced understanding of dysfunction of PI3K/Akt pathway in patients with these non-malignant disorders has offered a solid foundation for the progress of updated therapeutic targets. This article reviews the latest advances in the roles of PI3K/Akt pathway and their targets in the skin homeostasis and progression of a wide range of non-malignant skin disorders and describes the current progress in preclinical and clinical researches on the involvement of PI3K/Akt pathway targeted therapies.

HRT Atlas v1.0 database: redefining human and mouse housekeeping genes and candidate reference transcripts by mining massive RNA-seq datasets

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.

Terminus-Associated Non-coding RNAs: Trash or Treasure?

3′ untranslated regions (3′ UTRs) of protein-coding genes are well known for their important roles in determining the fate of mRNAs in diverse processes, including trafficking, stabilization, translation, and RNA–protein interactions. However, non-coding RNAs (ncRNAs) scattered around 3′ termini of the protein-coding genes, here referred to as terminus-associated non-coding RNAs (TANRs), have not attracted wide attention in RNA research. Indeed, whether TANRs are transcriptional noise, degraded mRNA products, alternative 3′ UTRs, or functional molecules has remained unclear for a long time. As a new category of ncRNAs, TANRs are widespread, abundant, and conserved in diverse eukaryotes. The biogenesis of TANRs mainly follows the same promoter model, the RNA-dependent RNA polymerase activity-dependent model, or the independent promoter model. Functional studies of TANRs suggested that they are significantly involved in the versatile regulation of gene expression. For instance, at the transcriptional level, they can lead to transcriptional interference, induce the formation of gene loops, and participate in transcriptional termination. Furthermore, at the posttranscriptional level, they can act as microRNA sponges, and guide cleavage or modification of target RNAs. Here, we review current knowledge of the potential role of TANRs in the modulation of gene expression. In this review, we comprehensively summarize the current state of knowledge about TANRs, and discuss TANR nomenclature, relation to ncRNAs, cross-talk biogenesis pathways and potential functions. We further outline directions of future studies of TANRs, to promote investigations of this emerging and enigmatic category of RNA.

Genetic architecture of complex traits and disease risk predictors

Genomic prediction of complex human traits (e.g., height, cognitive ability, bone density) and disease risks (e.g., breast cancer, diabetes, heart disease, atrial fibrillation) has advanced considerably in recent years. Using data from the UK Biobank, predictors have been constructed using penalized algorithms that favor sparsity: i.e., which use as few genetic variants as possible. We analyze the specific genetic variants (SNPs) utilized in these predictors, which can vary from dozens to as many as thirty thousand. We find that the fraction of SNPs in or near genic regions varies widely by phenotype. For the majority of disease conditions studied, a large amount of the variance is accounted for by SNPs outside of coding regions. The state of these SNPs cannot be determined from exome-sequencing data. This suggests that exome data alone will miss much of the heritability for these traits-i.e., existing PRS cannot be computed from exome data alone. We also study the fraction of SNPs and of variance that is in common between pairs of predictors. The DNA regions used in disease risk predictors so far constructed seem to be largely disjoint (with a few interesting exceptions), suggesting that individual genetic disease risks are largely uncorrelated. It seems possible in theory for an individual to be a low-risk outlier in all conditions simultaneously.

Metabolic Engineering Strategies in Diatoms Reveal Unique Phenotypes and Genetic Configurations With Implications for Algal Genetics and Synthetic Biology

Diatoms are photosynthetic microeukaryotes that dominate phytoplankton populations and have increasing applicability in biotechnology. Uncovering their complex biology and elevating strains to commercial standards depends heavily on robust genetic engineering tools. However, engineering microalgal genomes predominantly relies on random integration of transgenes into nuclear DNA, often resulting in detrimental “position-effects” such as transgene silencing, integration into transcriptionally-inactive regions, and endogenous sequence disruption. With the recent development of extrachromosomal transgene expression via independent episomes, it is timely to investigate both strategies at the phenotypic and genomic level. Here, we engineered the model diatom Phaeodactylum tricornutum to produce the high-value heterologous monoterpenoid geraniol, which, besides applications as fragrance and insect repellent, is a key intermediate of high-value pharmaceuticals. Using high-throughput phenotyping we confirmed the suitability of episomes for synthetic biology applications and identified superior geraniol-yielding strains following random integration. We used third generation long-read sequencing technology to generate a complete analysis of all transgene integration events including their genomic locations and arrangements associated with high-performing strains at a genome-wide scale with subchromosomal detail, never before reported in any microalga. This revealed very large, highly concatenated insertion islands, offering profound implications on diatom functional genetics and next generation genome editing technologies, and is key for developing more precise genome engineering approaches in diatoms, including possible genomic safe harbour locations to support high transgene expression for targeted integration approaches. Furthermore, we have demonstrated that exogenous DNA is not integrated inadvertently into the nuclear genome of extrachromosomal-expression clones, an important characterisation of this novel engineering approach that paves the road to synthetic biology applications.

Innate Immune Response Regulation by the Human RNASET2 Tumor Suppressor Gene

The link between cancer development or progression and immune system dysregulation has long been established. Virtually every cell type belonging to both the innate and adaptive immune system has been reported to be involved in a complex interplay that might culminate into either a pro- or anti-tumorigenic response. Among the cellular components of the innate immune system, cells belonging to the monocyte/macrophage lineage have been consistently shown to play a key role in the tumorigenic process. The most advanced human tumors are reported to be strongly infiltrated with Tumor-Associated Macrophages (TAMs) endowed with the ability to contribute to tumor growth and dissemination. However, given their widely acknowledged functional plasticity, macrophages can display anti-tumor properties as well. Based on these premises, experimental approaches to promote the in vivo macrophage shift from pro-tumor to anti-tumor phenotype represent one of the most promising research field aimed at developing immune system-mediated tumor suppressive therapies. In this context, the human RNASET2 oncosuppressor gene has emerged as a potential tool for macrophage-mediated tumor suppression. A growing body of experimental evidence has been reported to suggest a role for this gene in the regulation of macrophage activity in both in vitro and in vivo experimental models. Moreover, several recent reports suggest a role for this gene in a broad range of cell types involved in immune response, pointing at RNASET2 as a putative regulator of several functional features within the immune system.

Divergent and convergent evolution of housekeeping genes in human-pig lineage

Housekeeping genes are ubiquitously expressed and maintain basic cellular functions across tissue/cell type conditions. The present study aimed to develop a set of pig housekeeping genes and compare the structure, evolution and function of housekeeping genes in the human–pig lineage. By using RNA sequencing data, we identified 3,136 pig housekeeping genes. Compared with human housekeeping genes, we found that pig housekeeping genes were longer and subjected to slightly weaker purifying selection pressure and faster neutral evolution. Common housekeeping genes, shared by the two species, achieve stronger purifying selection than species-specific genes. However, pig- and human-specific housekeeping genes have similar functions. Some species-specific housekeeping genes have evolved independently to form similar protein active sites or structure, such as the classical catalytic serine–histidine–aspartate triad, implying that they have converged for maintaining the basic cellular function, which allows them to adapt to the environment. Human and pig housekeeping genes have varied structures and gene lists, but they have converged to maintain basic cellular functions essential for the existence of a cell, regardless of its specific role in the species. The results of our study shed light on the evolutionary dynamics of housekeeping genes.

Modeling Transcriptional Rewiring in Neutrophils Through the Course of Treated Juvenile Idiopathic Arthritis

Neutrophils in children with the polyarticular form of juvenile idiopathic arthritis (JIA) display abnormal transcriptional patterns linked to fundamental metabolic derangements. In this study, we sought to determine the effects of therapy on mRNA and miRNA expression networks in polyarticular JIA. Using exon and miRNA microarrays, we studied children with untreated active JIA (ADU, n = 35), children with active disease on therapy with methotrexate ± etanercept (ADT, n = 26), and children with inactive disease also on therapy (ID, n = 14). We compared the results to findings from healthy control children (HC, n = 35). We found substantial re-ordering of mRNA and miRNA expression networks after the initiation of therapy. Each disease state was associated with a distinct transcriptional profile, with the ADT state differing the most from HC, and ID more strongly resembling HC. Changes at the mRNA level were mirrored in changes in miRNA expression patterns. The analysis of the expression dynamics from differentially expressed genes across three disease states indicated that therapeutic response is a complex process. This process does not simply involve genes slowly correcting in a linear fashion over time. Computational modeling of miRNA and transcription factor (TF) co-regulatory networks demonstrated that combinational regulation of miRNA and TF might play an important role in dynamic transcriptome changes.

Classification of Transcription Boundary-Associated RNAs (TBARs) in Animals and Plants

There is increasing evidence suggesting the contribution of non-coding RNAs (ncRNAs) to the phenotypic and physiological complexity of organisms. A novel ncRNA species has been identified near the transcription boundaries of protein-coding genes in eukaryotes, bacteria, and archaea. This review provides a detailed description of these transcription boundary-associated RNAs (TBARs), including their classification. Based on their genomic distribution, TBARs are divided into two major groups: promoter-associated RNAs (PARs) and terminus-associated RNAs (TARs). Depending on the sequence length, each group is further classified into long RNA species (>200 nt) and small RNA species (<200 nt). According to these rules of TBAR classification, divergent ncRNAs with confusing nomenclatures, such as promoter upstream transcripts (PROMPTs), upstream antisense RNAs (uaRNAs), stable unannotated transcripts (SUTs), cryptic unstable transcripts (CUTs), upstream non-coding transcripts (UNTs), transcription start site-associated RNAs (TSSaRNAs), transcription initiation RNAs (tiRNAs), and transcription termination site-associated RNAs (TTSaRNAs), were assigned to specific classes. Although the biogenesis pathways of PARs and TARs have not yet been clearly elucidated, previous studies indicate that some of the PARs have originated either through divergent transcription or via RNA polymerase pausing. Intriguing findings regarding the functional implications of the TBARs such as the long-range “gene looping” model, which explains their role in the transcriptional regulation of protein-coding genes, are also discussed. Altogether, this review provides a comprehensive overview of the current research status of TBARs, which will promote further investigations in this research area.

Conserved noncoding transcription and core promoter regulatory code in early Drosophila development

Multicellular development is driven by regulatory programs that orchestrate the transcription of protein-coding and noncoding genes. To decipher this genomic regulatory code, and to investigate the developmental relevance of noncoding transcription, we compared genome-wide promoter activity throughout embryogenesis in 5 Drosophila species. Core promoters, generally not thought to play a significant regulatory role, in fact impart restrictions on the developmental timing of gene expression on a global scale. We propose a hierarchical regulatory model in which core promoters define broad windows of opportunity for expression, by defining a range of transcription factors from which they can receive regulatory inputs. This two-tiered mechanism globally orchestrates developmental gene expression, including extremely widespread noncoding transcription. The sequence and expression specificity of noncoding RNA promoters are evolutionarily conserved, implying biological relevance. Overall, this work introduces a hierarchical model for developmental gene regulation, and reveals a major role for noncoding transcription in animal development.

Bortezomib-induced miRNAs direct epigenetic silencing of locus genes and trigger apoptosis in leukemia

MicroRNAs (miRNAs) have been suggested to repress transcription via binding the 3′-untranslated regions of mRNAs. However, the involvement and details of miRNA-mediated epigenetic regulation, particularly in targeting genomic DNA and mediating epigenetic regulation, remain largely uninvestigated. In the present study, transcription factor CCAAT/enhancer binding protein delta (CEBPD) was responsive to the anticancer drug bortezomib, a clinical and highly selective drug for leukemia treatment, and contributed to bortezomib-induced cell death. Interestingly, following the identification of CEBPD-induced miRNAs, we found that miR-744, miR-3154 and miR-3162 could target CpG islands in the 5′-flanking region of the CEBPD gene. We previously demonstrated that the Yin Yang 1 (YY1)/polycomb group (PcG) protein/DNA methyltransferase (DNMT) complex is important for CCAAT/enhancer binding protein delta (CEBPD) gene inactivation; we further found that Argonaute 2 (Ago2) interacts with YY1 and binds to the CEBPD promoter. The miRNA/Ago2/YY1/PcG group protein/DNMT complex linked the inactivation of CEBPD and genes adjacent to its 5′-flanking region, including protein kinase DNA-activated catalytic polypeptide (PRKDC), minichromosome maintenance-deficient 4 (MCM4) and ubiquitin-conjugating enzyme E2 variant 2 (UBE2V2), upon bortezomib treatment. Moreover, we revealed that miRNA binding is necessary for YY1/PcG group protein/DNMT complex-mediated epigenetic gene silencing and is associated with bortezomib-induced methylation on genomic DNA. The present study successfully characterized the interactions of the miRNA/Ago2/YY1/PcG group protein/DNMT complex and provided new insights for miRNA-mediated epigenetic regulation in bortezomib-induced leukemic cell arrest and cell death.

LA ONTOGENIA Y LA EVOLUCIÓN DESDE LA PERSPECTIVA DE LA TEORÍA DE LOS SISTEMAS DE DESARROLLO (TSD)

La teoría de los sistemas de desarrollo (TSD) pretende realizar una síntesis conceptual que vincule el desarrollo ontogenético con la evolución. Sus antecedentes pueden ser encontrados básicamente en los trabajos de Waddington y de Bertalanffy quienes aportaron las bases de la canalización del desarrollo y la teoría de sistemas biológicos, respectivamente. El objetivo de este artículo es realizar un análisis conceptual preliminar de la TSD y reflexionar acerca de los aportes potenciales de la TSD como marco teórico para la biología del desarrollo en particular y la biología evolutiva en general. Para ello, se tendrán en cuenta algunos de los conceptos y propuestas que componen este marco y se trabajará sobre datos secundarios obtenidos de la bibliografía. Se concluye que la TSD: 1-logra argumentar en contra de la visión gen-centrista respecto de las explicaciones que pretenden justificar el desarrollo biológico y evolutivo; 2- argumenta de manera coherente a favor del rol de la epigenética en la ontogenia y la evolución; 3- en relación con lo anterior el rol de la selección natural se restringe a un segundo plano; 4- propone que la dicotomía naturaleza/cultura debe ser superada; y 5- constituye un posible programa de investigación metodológico compuesto de una diversidad de hipótesis y teorías no necesariamente relacionadas que pueden ser corroboradas de manera relativamente independientes del resto de la red teórica.   

Molecular Resources from Transcriptomes in the Brassicaceae Family

The rapidly falling costs and the increasing availability of large DNA sequence data sets facilitate the fast and affordable mining of large molecular markers data sets for comprehensive evolutionary studies. The Brassicaceae (mustards) are an important species-rich family in the plant kingdom with taxa distributed worldwide and a complex evolutionary history. We performed Simple Sequence Repeats (SSRs) mining using de novo assembled transcriptomes from 19 species across the Brassicaceae in order to study SSR evolution and provide comprehensive sets of molecular markers for genetic studies within the family. Moreover, we selected the genus Cochlearia to test the transferability and polymorphism of these markers among species. Additionally, we annotated Cochlearia pyrenaica transcriptome in order to identify the position of each of the mined SSRs. While we introduce a new set of tools that will further enable evolutionary studies across the Brassicaceae, we also discuss some broader aspects of SSR evolution. Overall, we developed 2012 ready-to-use SSR markers with their respective primers in 19 Brassicaceae species and a high quality annotated transcriptome for C. pyrenaica. As indicated by our transferability test with the genus Cochlearia these SSRs are transferable to species within the genus increasing exponentially the number of targeted species. Also, our polymorphism results showed substantial levels of variability for these markers. Finally, despite its complex evolutionary history, SSR evolution across the Brassicaceae family is highly conserved and we found no deviation from patterns reported in other Angiosperms.

The Evolving Definition of the Term "Gene"

This paper presents a history of the changing meanings of the term "gene," over more than a century, and a discussion of why this word, so crucial to genetics, needs redefinition today. In this account, the first two phases of 20th century genetics are designated the "classical" and the "neoclassical" periods, and the current molecular-genetic era the "modern period." While the first two stages generated increasing clarity about the nature of the gene, the present period features complexity and confusion. Initially, the term "gene" was coined to denote an abstract "unit of inheritance," to which no specific material attributes were assigned. As the classical and neoclassical periods unfolded, the term became more concrete, first as a dimensionless point on a chromosome, then as a linear segment within a chromosome, and finally as a linear segment in the DNA molecule that encodes a polypeptide chain. This last definition, from the early 1960s, remains the one employed today, but developments since the 1970s have undermined its generality. Indeed, they raise questions about both the utility of the concept of a basic "unit of inheritance" and the long implicit belief that genes are autonomous agents. Here, we review findings that have made the classic molecular definition obsolete and propose a new one based on contemporary knowledge.

Between-species differences in gene copy number are enriched among functions critical for adaptive evolution in Arabidopsis halleri

Background

Gene copy number divergence between species is a form of genetic polymorphism that contributes significantly to both genome size and phenotypic variation. In plants, copy number expansions of single genes were implicated in cultivar- or species-specific tolerance of high levels of soil boron, aluminium or calamine-type heavy metals, respectively. Arabidopsis halleri is a zinc- and cadmium-hyperaccumulating extremophile species capable of growing on heavy-metal contaminated, toxic soils. In contrast, its non-accumulating sister species A. lyrata and the closely related reference model species A. thaliana exhibit merely basal metal tolerance.

Results

For a genome-wide assessment of the role of copy number divergence (CND) in lineage-specific environmental adaptation, we conducted cross-species array comparative genome hybridizations of three plant species and developed a global signal scaling procedure to adjust for sequence divergence. In A. halleri, transition metal homeostasis functions are enriched twofold among the genes detected as copy number expanded. Moreover, biotic stress functions including mostly disease Resistance (R) gene-related genes are enriched twofold among genes detected as copy number reduced, when compared to the abundance of these functions among all genes.

Conclusions

Our results provide genome-wide support for a link between evolutionary adaptation and CND in A. halleri as shown previously for Heavy metal ATPase4. Moreover our results support the hypothesis that elemental defences, which result from the hyperaccumulation of toxic metals, allow the reduction of classical defences against biotic stress as a trade-off.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3319-5) contains supplementary material, which is available to authorized users.

What is an RNA? A top layer for RNA classification

Every ribonucleic acid begins its cellular life as a transcript. If the transcript or its processing product has a function it should be regarded an RNA. Nonfunctional transcripts, by-products from processing, degradation intermediates, even those originating from (functional) RNAs, and non-functional products of transcriptional gene regulation accomplished via the act of transcription, as well as stochastic (co)transcripts could simply be addressed as transcripts (class 0). The copious functional RNAs (class I), often maturing after one or more processing steps, already are systematized into ever expanding sub-classifications ranging from micro RNAs to rRNAs. Established sub-classifications addressing a wide functional diversity remain unaffected. mRNAs (class II) are distinct from any other RNA by virtue of their potential to be translated into (poly)peptide(s) on ribosomes. We are not proposing a novel RNA classification, but wish to add a basic concept with existing terminology (transcript, RNA, and mRNA) that should serve as an additional framework for carefully delineating RNA function from an avalanche of RNA sequencing data. At the same time, this top level hierarchical model should illuminate important principles of RNA evolution and biology thus heightening our awareness that in biology boundaries and categorizations are typically fuzzy.

LncRNAs: emerging players in gene regulation and disease pathogenesis

The advent of next-generation sequencing has demonstrated that eukaryotic genomes are extremely complex than what were previously thought. Recent studies revealed that in addition to protein-coding genes, nonprotein-coding genes have allocated a large fraction of the genome. Long noncoding RNA (lncRNA) genes are classified as nonprotein-coding genes, serving as a molecular signal, decoy, guide and scaffold. They were suggested to play important roles in chromatin states, epigenetic and posttranscriptional regulation of genes. Aberrant expression of lncRNAs and changes in their structure are associated with a wide spectrum of diseases ranging from different types of cancer and neurodegeneration to ?-thalassaemia. The purpose of this study was to summarize the current progress in understanding the genomic bases and origin of lncRNAs. Moreover, this study focusses on the diverse functions of lncRNAs in normal cells as well as various types of disease to illustrate the potential impacts of lncRNAs on diverse biological processes and their therapeutic significance.

Mammalian Synthetic Biology: Time for Big MACs

The enabling technologies of synthetic biology are opening up new opportunities for engineering and enhancement of mammalian cells. This will stimulate diverse applications in many life science sectors such as regenerative medicine, development of biosensing cell lines, therapeutic protein production, and generation of new synthetic genetic regulatory circuits. Harnessing the full potential of these new engineering-based approaches requires the design and assembly of large DNA constructs-potentially up to chromosome scale-and the effective delivery of these large DNA payloads to the host cell. Random integration of large transgenes, encoding therapeutic proteins or genetic circuits into host chromosomes, has several drawbacks such as risks of insertional mutagenesis, lack of control over transgene copy-number and position-specific effects; these can compromise the intended functioning of genetic circuits. The development of a system orthogonal to the endogenous genome is therefore beneficial. Mammalian artificial chromosomes (MACs) are functional, add-on chromosomal elements, which behave as normal chromosomes-being replicating and portioned to daughter cells at each cell division. They are deployed as useful gene expression vectors as they remain independent from the host genome. MACs are maintained as a single-copy and can accommodate multiple gene expression cassettes of, in theory, unlimited DNA size (MACs up to 10 megabases have been constructed). MACs therefore enabled control over ectopic gene expression and represent an excellent platform to rapidly prototype and characterize novel synthetic gene circuits without recourse to engineering the host genome. This review describes the obstacles synthetic biologists face when working with mammalian systems and how the development of improved MACs can overcome these-particularly given the spectacular advances in DNA synthesis and assembly that are fuelling this research area.

MicroRNA in Control of Gene Expression: An Overview of Nuclear Functions

The finding that small non-coding RNAs (ncRNAs) are able to control gene expression in a sequence specific manner has had a massive impact on biology. Recent improvements in high throughput sequencing and computational prediction methods have allowed the discovery and classification of several types of ncRNAs. Based on their precursor structures, biogenesis pathways and modes of action, ncRNAs are classified as small interfering RNAs (siRNAs), microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs), endogenous small interfering RNAs (endo-siRNAs or esiRNAs), promoter associate RNAs (pRNAs), small nucleolar RNAs (snoRNAs) and sno-derived RNAs. Among these, miRNAs appear as important cytoplasmic regulators of gene expression. miRNAs act as post-transcriptional regulators of their messenger RNA (mRNA) targets via mRNA degradation and/or translational repression. However, it is becoming evident that miRNAs also have specific nuclear functions. Among these, the most studied and debated activity is the miRNA-guided transcriptional control of gene expression. Although available data detail quite precisely the effectors of this activity, the mechanisms by which miRNAs identify their gene targets to control transcription are still a matter of debate. Here, we focus on nuclear functions of miRNAs and on alternative mechanisms of target recognition, at the promoter lavel, by miRNAs in carrying out transcriptional gene silencing.

Long non-coding RNAs in human early embryonic development and their potential in ART

BACKGROUND

Human long non-coding RNAs (lncRNAs) are an emerging category of transcripts with increasingly documented functional roles during development. LncRNAs and roles during human early embryo development have recently begun to be unravelled.

OBJECTIVE AND RATIONALE

This review summarizes the most recent knowledge on lncRNAs and focuses on their expression patterns and role during early human embryo development and in pluripotent stem cells (PSCs). Public mRNA sequencing (mRNA-seq) data were used to illustrate these expression signatures.

SEARCH METHODS

The PubMed and EMBASE databases were first interrogated using specific terms, such as 'lncRNAs', to get an extensive overview on lncRNAs up to February 2016, and then using 'human lncRNAs' and 'embryo', 'development', or 'PSCs' to focus on lncRNAs involved in human embryo development or in PSC.Recently published RNA-seq data from human oocytes and pre-implantation embryos (including single-cell data), PSC and a panel of normal and malignant adult tissues were used to describe the specific expression patterns of some lncRNAs in early human embryos.

OUTCOMES

The existence and the crucial role of lncRNAs in many important biological phenomena in each branch of the life tree are now well documented. The number of identified lncRNAs is rapidly increasing and has already outnumbered that of protein-coding genes. Unlike small non-coding RNAs, a variety of mechanisms of action have been proposed for lncRNAs. The functional role of lncRNAs has been demonstrated in many biological and developmental processes, including cell pluripotency induction, X-inactivation or gene imprinting. Analysis of RNA-seq data highlights that lncRNA abundance changes significantly during human early embryonic development. This suggests that lncRNAs could represent candidate biomarkers for developing non-invasive tests for oocyte or embryo quality. Finally, some of these lncRNAs are also expressed in human cancer tissues, suggesting that reactivation of an embryonic lncRNA program may contribute to human malignancies.

WIDER IMPLICATIONS

LncRNAs are emerging potential key players in gene expression regulation. Analysis of RNA-seq data from human pre-implantation embryos identified lncRNA signatures that are specific to this critical step. We anticipate that further studies will show that these new transcripts are major regulators of embryo development. These findings might also be used to develop new tests/treatments for improving the pregnancy success rate in IVF procedures or for regenerative medicine applications involving PSC.

Complexity and Specificity of the Neutrophil Transcriptomes in Juvenile Idiopathic Arthritis

NIH projects such as ENCODE and Roadmap Epigenomics have revealed surprising complexity in the transcriptomes of mammalian cells. In this study, we explored transcriptional complexity in human neutrophils, cells generally regarded as nonspecific in their functions and responses. We studied distinct human disease phenotypes and found that, at the gene, gene isoform, and miRNA level, neutrophils exhibit considerable specificity in their transcriptomes. Thus, even cells whose responses are considered non-specific show tailoring of their transcriptional repertoire toward specific physiologic or pathologic contexts. We also found that miRNAs had a global impact on neutrophil transcriptome and are associated with innate immunity in juvenile idiopathic arthritis (JIA). These findings have important implications for our understanding of the link between genes, non-coding transcripts and disease phenotypes.

Mining microsatellite markers from public expressed sequence tags databases for the study of threatened plants

BACKGROUND

Simple Sequence Repeats (SSRs) are widely used in population genetic studies but their classical development is costly and time-consuming. The ever-increasing available DNA datasets generated by high-throughput techniques offer an inexpensive alternative for SSRs discovery. Expressed Sequence Tags (ESTs) have been widely used as SSR source for plants of economic relevance but their application to non-model species is still modest.

METHODS

Here, we explored the use of publicly available ESTs (GenBank at the National Center for Biotechnology Information-NCBI) for SSRs development in non-model plants, focusing on genera listed by the International Union for the Conservation of Nature (IUCN). We also search two model genera with fully annotated genomes for EST-SSRs, Arabidopsis and Oryza, and used them as controls for genome distribution analyses. Overall, we downloaded 16 031 555 sequences for 258 plant genera which were mined for SSRsand their primers with the help of QDD1. Genome distribution analyses in Oryza and Arabidopsis were done by blasting the sequences with SSR against the Oryza sativa and Arabidopsis thaliana reference genomes implemented in the Basal Local Alignment Tool (BLAST) of the NCBI website. Finally, we performed an empirical test to determine the performance of our EST-SSRs in a few individuals from four species of two eudicot genera, Trifolium and Centaurea.

RESULTS

We explored a total of 14 498 726 EST sequences from the dbEST database (NCBI) in 257 plant genera from the IUCN Red List. We identify a very large number (17 102) of ready-to-test EST-SSRs in most plant genera (193) at no cost. Overall, dinucleotide and trinucleotide repeats were the prevalent types but the abundance of the various types of repeat differed between taxonomic groups. Control genomes revealed that trinucleotide repeats were mostly located in coding regions while dinucleotide repeats were largely associated with untranslated regions. Our results from the empirical test revealed considerable amplification success and transferability between congenerics.

CONCLUSIONS

The present work represents the first large-scale study developing SSRs by utilizing publicly accessible EST databases in threatened plants. Here we provide a very large number of ready-to-test EST-SSR (17 102) for 193 genera. The cross-species transferability suggests that the number of possible target species would be large. Since trinucleotide repeats are abundant and mainly linked to exons they might be useful in evolutionary and conservation studies. Altogether, our study highly supports the use of EST databases as an extremely affordable and fast alternative for SSR developing in threatened plants.

Disease-Associated Single-Nucleotide Polymorphisms From Noncoding Regions in Juvenile Idiopathic Arthritis Are Located Within or Adjacent to Functional Genomic Elements of Human Neutrophils and CD4+ T Cells

OBJECTIVE

Juvenile idiopathic arthritis (JIA) is considered a complex disease in which the environment interacts with inherited genes to produce a phenotype that shows broad interindividual variability. Twenty-four regions of genetic risk for JIA were identified in a recent genome-wide association study (GWAS); however, as is typical of the results of GWAS, most of the regions of genetic risk (22 of 24) were in noncoding regions of the genome. This study was undertaken to identify functional elements (other than genes) that might be located within the regions of genetic risk.

METHODS

We used paired-end RNA sequencing to identify noncoding RNAs (ncRNAs) located within 5 kb of disease-associated single-nucleotide polymorphisms (SNPs). In addition, we used chromatin immunoprecipitation (ChIP) followed by sequencing to identify epigenetic marks associated with enhancer function (H3K4me1 and H3K27ac) in human neutrophils to determine whether enhancer-associated histone marks were enriched in the linkage disequilibrium blocks that encompassed the 22 SNPs identified in the GWAS.

RESULTS

In human neutrophils, we identified H3K4me1 and/or H3K27ac marks in 15 of the 22 regions previously identified as risk loci for JIA. In CD4+ T cells, 18 regions had H3K4me1 and/or H3K27ac marks. In addition, we identified ncRNA transcripts at the rs4705862 and rs6894249 loci in human neutrophils.

CONCLUSION

Much of the genetic risk for JIA lies within or adjacent to regions of neutrophil and CD4+ T cell genomes that carry epigenetic marks associated with enhancer function and/or ncRNA transcripts. These findings are consistent with the hypothesis that JIA is fundamentally a disorder of gene regulation that includes both the innate and the adaptive immune system. Elucidating the specific roles of these noncoding elements within leukocyte genomes will be critical to our understanding of JIA pathogenesis.

RNA sequencing from human neutrophils reveals distinct transcriptional differences associated with chronic inflammatory states

BACKGROUND

The transcriptional complexity of mammalian cells suggests that they have broad abilities to respond to specific environmental stimuli and physiologic contexts. These abilities were not apparent a priori from the structure of mammalian genomes, but have been identified through detailed transcriptome analyses. In this study, we examined the transcriptomes of cells of the innate immune system, human neutrophils, using RNA sequencing (RNAseq).

METHODS

We sequenced poly-A RNA from nine individual samples corresponding to specific phenotypes: three children with active, untreated juvenile idiopathic arthritis (JIA)(AD), three children with the same disease whose disease was inactive on medication (CRM), and three children with cystic fibrosis (CF).

RESULTS

We demonstrate that transcriptomes of neutrophils, typically considered non-specific in their responses and functions, display considerable specificity in their transcriptional repertoires dependent on the pathologic context, and included genes, gene isoforms, and long non-coding RNA transcripts. Furthermore, despite the small sample numbers, these findings demonstrate the potential of RNAseq approaches to biomarker development in rheumatic diseases.

CONCLUSIONS

These data demonstrate the capacity of cells previously considered non-specific in function to adapt their transcriptomes to specific biologic contexts. These data also provide insight into previously unrecognized pathological pathways and show considerable promise for elucidating disease and disease-state specific regulatory networks.

Regulation of alternative gene transcription in the striatum in response to antidepressant drugs

The mechanisms that control the selection of transcription initiation and termination sites in response to pharmacological stimulation of neuronal cells are poorly understood. We used next-generation sequencing and bioinformatics to construct a genome-wide inventory of protein-coding and non-coding transcripts altered by antidepressant treatment. We analyzed available ChIP-seq data to identify mechanisms that control drug-inducible expression of alternative gene variants in the brain. We identified 153 transcripts of various biotypes regulated in the mouse striatum in response to tranylcypromine or mianserin (at a 0.1% FDR threshold). Five drug-responsive gene patterns are enriched in protein-coding variants (77%), regulated by different sets of transcriptional factors (including SRF/CREB1 and GR/CTCF) and expressed in separate cellular compartments of the brain. We found that alterations mediated by proximal promoters in neurons are more specific in the selection of regulated transcriptional isoforms compared with enhancer-dependent alterations in glia. The change in transcriptional programs, from housekeeping to inducible, provides cells with the resource of functionally distinct proteins. We conclude that the regulation of drug-induced brain plasticity may occur at the level of transcripts rather than genes. The expression of specific isoforms in response to antidepressants may constitute a trigger that initiates the long-lasting effects of these drugs.

Does every transcript originate from a gene?

Outdated gene definitions favored regions corresponding to mature messenger RNAs, in particular, the open reading frame. In eukaryotes, the intergenic space was widely regarded nonfunctional and devoid of RNA transcription. Original concepts were based on the assumption that RNA expression was restricted to known protein-coding genes and a few so-called structural RNA genes, such as ribosomal RNAs or transfer RNAs. With the discovery of introns and, more recently, sensitive techniques for monitoring genome-wide transcription, this view had to be substantially modified. Tiling microarrays and RNA deep sequencing revealed myriads of transcripts, which cover almost entire genomes. The tremendous complexity of non-protein-coding RNA transcription has to be integrated into novel gene definitions. Despite an ever-growing list of functional RNAs, questions concerning the mass of identified transcripts are under dispute. Here, we examined genome-wide transcription from various angles, including evolutionary considerations, and suggest, in analogy to novel alternative splice variants that do not persist, that the vast majority of transcripts represent raw material for potential, albeit rare, exaptation events.

Transcriptomic landscape of lncRNAs in inflammatory bowel disease

BACKGROUND

Inflammatory bowel disease (IBD) is a complex multi-factorial inflammatory disease with Crohn's disease (CD) and ulcerative colitis (UC) being the two most common forms. A number of transcriptional profiling studies have provided compelling evidence that describe the role of protein-coding genes and microRNAs in modulating the immune responses in IBD.

METHODS

In the present study, we performed a genome-wide transcriptome profiling of lncRNAs and protein-coding genes in 96 colon pinch biopsies (inflamed and non-inflamed) extracted from multiple colonic locations from 45 patients (CD = 13, UC = 20, controls = 12) using an expression microarray platform.

RESULTS

In our study, we identified widespread dysregulation of lncRNAs and protein-coding genes in both inflamed and non-inflamed CD and UC compared to the healthy controls. In cases of inflamed CD and UC, we identified 438 and 745 differentially expressed lncRNAs, respectively, while in cases of the non-inflamed CD and UC, we identified 12 and 19 differentially expressed lncRNAs, respectively. We also observed significant enrichment (P-value <0.001, Pearson's Chi-squared test) for 96 differentially expressed lncRNAs and 154 protein-coding genes within the IBD susceptibility loci. Furthermore, we found strong positive expression correlations for the intersecting and cis-neighboring differentially expressed IBD loci-associated lncRNA-protein-coding gene pairs. The functional annotation analysis of differentially expressed genes revealed their involvement in the immune response, pro-inflammatory cytokine activity and MHC protein complex.

CONCLUSIONS

The lncRNA expression profiling in both inflamed and non-inflamed CD and UC successfully stratified IBD patients from the healthy controls. Taken together, the identified lncRNA transcriptional signature along with clinically relevant parameters suggest their potential as biomarkers in IBD.

Repetitive element signature-based visualization, distance computation, and classification of 1766 microbial genomes

The genomes of living organisms are populated with pleomorphic repetitive elements (REs) of varying densities. Our hypothesis that genomic RE landscapes are species/strain/individual-specific was implemented into the Genome Signature Imaging system to visualize and compute the RE-based signatures of any genome. Following the occurrence profiling of 5-nucleotide REs/words, the information from top-50 frequency words was transformed into a genome-specific signature and visualized as Genome Signature Images (GSIs), using a CMYK scheme. An algorithm for computing distances among GSIs was formulated using the GSIs' variables (word identity, frequency, and frequency order). The utility of the GSI-distance computation system was demonstrated with control genomes. GSI-based computation of genome-relatedness among 1766 microbes (117 archaea and 1649 bacteria) identified their clustering patterns; although the majority paralleled the established classification, some did not. The Genome Signature Imaging system, with its visualization and distance computation functions, enables genome-scale evolutionary studies involving numerous genomes with varying sizes.

Long Non-coding RNA

Rapid development in genome-wide transcriptional analyses has led to the discovery of a large number of non-coding transcripts, also called long non-coding RNA (lncRNA). LncRNAs harbor biological activities including regulation of protein-coding gene expression at epigenetic, transcriptional and post-transcriptional levels. They also take a part in various physiological and pathological processes, participating in cell development, immunity, disease processes and oncogenesis.

Here I discuss and summarize, current knowledge about lncRNA origin, function and involvement in human disease.

Complex tissue-specific patterns and distribution of multiple RAGE splice variants in different mammals

The receptor for advanced glycosylation end products (RAGE) is a multiligand receptor involved in diverse cell signaling pathways. Previous studies show that this gene expresses several splice variants in human, mouse, and dog. Alternative splicing (AS) plays an important role in expanding transcriptomic and proteomic diversity, and it has been related to disease. AS is also one of the main evolutionary mechanisms in mammalian genomes. However, limited information is available regarding the AS of RAGE in a wide context of mammalian tissues. In this study, we examined in detail the different RAGE mRNAs generated by AS from six mammals, including two primates (human and monkey), two artiodactyla (cow and pig), and two rodentia (mouse and rat) in 6–18 different tissues including fetal, adult, and tumor. By nested reverse transcription-polymerase chain reaction (RT-PCR) we identified a high number of splice variants including noncoding transcripts and predicted coding ones with different potential protein modifications affecting mainly the transmembrane and ligand-binding domains that could influence their biological function. However, analysis of RNA-seq data enabled detecting only the most abundant splice variants. More than 80% of the detected RT-PCR variants (87 of 101 transcripts) are novel (different exon/intron structure to the previously described ones), and interestingly, 20–60% of the total transcripts (depending on the species) are noncoding ones that present tissue specificity. Our results suggest that RAGE undergoes extensive AS in mammals, with different expression patterns among adult, fetal, and tumor tissues. Moreover, most splice variants seem to be species specific, especially the noncoding variants, with only two (canonical human Tv1-RAGE, and human N-truncated or Tv10-RAGE) conserved among the six different species. This could indicate a special evolution pattern of this gene at mRNA level.

Fam40b is required for lineage commitment of murine embryonic stem cells

FAM40B (STRIP2) is a member of the striatin-interacting phosphatase and kinase (STRIPAK) complex that is involved in the regulation of various processes such as cell proliferation and differentiation. Its role for differentiation processes in embryonic stem cells (ESCs) is till now completely unknown. Short hairpin RNA (shRNA)-mediated silencing of Fam40b expression in ESCs and differentiating embryoid bodies (EBs) led to perturbed differentiation to embryonic germ layers and their derivatives including a complete abrogation of cardiomyogenesis. Pluripotency factors such as Nanog, Oct4 and Sox2 as well as epigenetic factors such as histone acetyltransferase type B (HAT1) and DNA (cytosine-5)-methyltransferase 3-β (Dnmt3b) were highly upregulated in Fam40b knockdown EBs as compared with control and scrambled EBs. To examine the relevance of Fam40b for development in vivo, Fam40b was knocked down in developing zebrafish. Morpholino-mediated knockdown of Fam40b led to severe abnormalities of the cardiovascular system, including an impaired expression of ventricular myosin heavy chain (vmhc) and of cardiac myosin light chain 2 (cmlc2) in the heart. We identified the gene product of Fam40b in ESCs as a perinuclear and nucleolar protein with a molecular weight of 96 kDa. We conclude that the expression of Fam40b is essential for the lineage commitment of murine embryonic stem cells (mESCs) into differentiated somatic cells via mechanisms involving pluripotency and epigenetic networks.

The rise of regulatory RNA

Discoveries over the past decade portend a paradigm shift in molecular biology. Evidence suggests that RNA is not only functional as a messenger between DNA and protein but also involved in the regulation of genome organization and gene expression, which is increasingly elaborate in complex organisms. Regulatory RNA seems to operate at many levels; in particular, it plays an important part in the epigenetic processes that control differentiation and development. These discoveries suggest a central role for RNA in human evolution and ontogeny. Here, we review the emergence of the previously unsuspected world of regulatory RNA from a historical perspective.

Choice of transcripts and software has a large effect on variant annotation

Background

Variant annotation is a crucial step in the analysis of genome sequencing data. Functional annotation results can have a strong influence on the ultimate conclusions of disease studies. Incorrect or incomplete annotations can cause researchers both to overlook potentially disease-relevant DNA variants and to dilute interesting variants in a pool of false positives. Researchers are aware of these issues in general, but the extent of the dependency of final results on the choice of transcripts and software used for annotation has not been quantified in detail.

Methods

This paper quantifies the extent of differences in annotation of 80 million variants from a whole-genome sequencing study. We compare results using the RefSeq and Ensembl transcript sets as the basis for variant annotation with the software Annovar, and also compare the results from two annotation software packages, Annovar and VEP (Ensembl’s Variant Effect Predictor), when using Ensembl transcripts.

Results

We found only 44% agreement in annotations for putative loss-of-function variants when using the RefSeq and Ensembl transcript sets as the basis for annotation with Annovar. The rate of matching annotations for loss-of-function and nonsynonymous variants combined was 79% and for all exonic variants it was 83%. When comparing results from Annovar and VEP using Ensembl transcripts, matching annotations were seen for only 65% of loss-of-function variants and 87% of all exonic variants, with splicing variants revealed as the category with the greatest discrepancy. Using these comparisons, we characterised the types of apparent errors made by Annovar and VEP and discuss their impact on the analysis of DNA variants in genome sequencing studies.

Conclusions

Variant annotation is not yet a solved problem. Choice of transcript set can have a large effect on the ultimate variant annotations obtained in a whole-genome sequencing study. Choice of annotation software can also have a substantial effect. The annotation step in the analysis of a genome sequencing study must therefore be considered carefully, and a conscious choice made as to which transcript set and software are used for annotation.

Stem cell systems informatics for advanced clinical biodiagnostics: tracing molecular signatures from bench to bedside

Development of innovative high throughput technologies has enabled a variety of molecular landscapes to be interrogated with an unprecedented degree of detail. Emergence of next generation nucleotide sequencing methods, advanced proteomic techniques, and metabolic profiling approaches continue to produce a wealth of biological data that captures molecular frameworks underlying phenotype. The advent of these novel technologies has significant translational applications, as investigators can now explore molecular underpinnings of developmental states with a high degree of resolution. Application of these leading-edge techniques to patient samples has been successfully used to unmask nuanced molecular details of disease vs healthy tissue, which may provide novel targets for palliative intervention. To enhance such approaches, concomitant development of algorithms to reprogram differentiated cells in order to recapitulate pluripotent capacity offers a distinct advantage to advancing diagnostic methodology. Bioinformatic deconvolution of several “-omic” layers extracted from reprogrammed patient cells, could, in principle, provide a means by which the evolution of individual pathology can be developmentally monitored. Significant logistic challenges face current implementation of this novel paradigm of patient treatment and care, however, several of these limitations have been successfully addressed through continuous development of cutting edge in silico archiving and processing methods. Comprehensive elucidation of genomic, transcriptomic, proteomic, and metabolomic networks that define normal and pathological states, in combination with reprogrammed patient cells are thus poised to become high value resources in modern diagnosis and prognosis of patient disease.

Functional transcriptomics in the post-ENCODE era

The last decade has seen tremendous effort committed to the annotation of the human genome sequence, most notably perhaps in the form of the ENCODE project. One of the major findings of ENCODE, and other genome analysis projects, is that the human transcriptome is far larger and more complex than previously thought. This complexity manifests, for example, as alternative splicing within protein-coding genes, as well as in the discovery of thousands of long noncoding RNAs. It is also possible that significant numbers of human transcripts have not yet been described by annotation projects, while existing transcript models are frequently incomplete. The question as to what proportion of this complexity is truly functional remains open, however, and this ambiguity presents a serious challenge to genome scientists. In this article, we will discuss the current state of human transcriptome annotation, drawing on our experience gained in generating the GENCODE gene annotation set. We highlight the gaps in our knowledge of transcript functionality that remain, and consider the potential computational and experimental strategies that can be used to help close them. We propose that an understanding of the true overlap between transcriptional complexity and functionality will not be gained in the short term. However, significant steps toward obtaining this knowledge can now be taken by using an integrated strategy, combining all of the experimental resources at our disposal.

Understanding the regulatory and transcriptional complexity of the genome through structure

An expansive functionality and complexity has been ascribed to the majority of the human genome that was unanticipated at the outset of the draft sequence and assembly a decade ago. We are now faced with the challenge of integrating and interpreting this complexity in order to achieve a coherent view of genome biology. We argue that the linear representation of the genome exacerbates this complexity and an understanding of its three-dimensional structure is central to interpreting the regulatory and transcriptional architecture of the genome. Chromatin conformation capture techniques and high-resolution microscopy have afforded an emergent global view of genome structure within the nucleus. Chromosomes fold into complex, territorialized three-dimensional domains in concert with specialized subnuclear bodies that harbor concentrations of transcription and splicing machinery. The signature of these folds is retained within the layered regulatory landscapes annotated by chromatin immunoprecipitation, and we propose that genome contacts are reflected in the organization and expression of interweaved networks of overlapping coding and noncoding transcripts. This pervasive impact of genome structure favors a preeminent role for the nucleoskeleton and RNA in regulating gene expression by organizing these folds and contacts. Accordingly, we propose that the local and global three-dimensional structure of the genome provides a consistent, integrated, and intuitive framework for interpreting and understanding the regulatory and transcriptional complexity of the human genome.

RIG-I and MDA-5 detection of viral RNA-dependent RNA polymerase activity restricts positive-strand RNA virus replication

Type I interferons (IFN) are important for antiviral responses. Melanoma differentiation-associated gene 5 (MDA-5) and retinoic acid-induced gene I (RIG-I) proteins detect cytosolic double-stranded RNA (dsRNA) or 5'-triphosphate (5'-ppp) RNA and mediate IFN production. Cytosolic 5'-ppp RNA and dsRNA are generated during viral RNA replication and transcription by viral RNA replicases [RNA-dependent RNA polymerases (RdRp)]. Here, we show that the Semliki Forest virus (SFV) RNA replicase can induce IFN-β independently of viral RNA replication and transcription. The SFV replicase converts host cell RNA into 5'-ppp dsRNA and induces IFN-β through the RIG-I and MDA-5 pathways. Inactivation of the SFV replicase RdRp activity prevents IFN-β induction. These IFN-inducing modified host cell RNAs are abundantly produced during both wild-type SFV and its non-pathogenic mutant infection. Furthermore, in contrast to the wild-type SFV replicase a non-pathogenic mutant replicase triggers increased IFN-β production, which leads to a shutdown of virus replication. These results suggest that host cells can restrict RNA virus replication by detecting the products of unspecific viral replicase RdRp activity.