Allison Piovesan

An estimation of the number of cells in the human body

Abstract Background: All living organisms are made of individual and identifiable cells, whose number, together with their size and type, ultimately defines the structure and functions of an organism. While the total cell number of lower organisms is often known, it has not yet been defined in higher organisms. In particular, the reported total cell number of a human being ranges between 1012 and 1016 and it is widely mentioned without a proper reference. Aim: To study and discuss the theoretical issue of the total number of cells that compose the standard human adult organism. Subjects and methods: A systematic calculation of the total cell number of the whole human body and of the single organs was carried out using bibliographical and/or mathematical approaches. Results: A current estimation of human total cell number calculated for a variety of organs and cell types is presented. These partial data correspond to a total number of 3.72 × 1013. Conclusions: Knowing the total cell number of the human body as well as of individual organs is important from a cultural, biological, medical and comparative modelling point of view. The presented cell count could be a starting point for a common effort to complete the total calculation.

Extensive microRNA-mediated crosstalk between lncRNAs and mRNAs in mouse embryonic stem cells

Recently, a handful of intergenic long noncoding RNAs (lncRNAs) have been shown to compete with mRNAs for binding to miRNAs and to contribute to development and disease. Beyond these reports, little is yet known of the extent and functional consequences of miRNA-mediated regulation of mRNA levels by lncRNAs. To gain further insight into lncRNA-mRNA miRNA-mediated crosstalk, we reanalyzed transcriptome-wide changes induced by the targeted knockdown of over 100 lncRNA transcripts in mouse embryonic stem cells (mESCs). We predicted that, on average, almost one-fifth of the transcript level changes induced by lncRNAs are dependent on miRNAs that are highly abundant in mESCs. We validated these findings experimentally by temporally profiling transcriptome-wide changes in gene expression following the loss of miRNA biogenesis in mESCs. Following the depletion of miRNAs, we found that >50% of lncRNAs and their miRNA-dependent mRNA targets were up-regulated coordinately, consistent with their interaction being miRNA-mediated. These lncRNAs are preferentially located in the cytoplasm, and the response elements for miRNAs they share with their targets have been preserved in mammals by purifying selection. Lastly, miRNA-dependent mRNA targets of each lncRNA tended to share common biological functions. Post-transcriptional miRNA-mediated crosstalk between lncRNAs and mRNA, in mESCs, is thus surprisingly prevalent, conserved in mammals, and likely to contribute to critical developmental processes.

Systematic reanalysis of partial trisomy 21 cases with or without Down syndrome suggests a small region on 21q22.13 as critical to the phenotype

A ‘Down Syndrome critical region’ (DSCR) sufficient to induce the most constant phenotypes of Down syndrome (DS) had been identified by studying partial (segmental) trisomy 21 (PT21) as an interval of 0.6–8.3 Mb within human chromosome 21 (Hsa21), although its existence was later questioned. We propose an innovative, systematic reanalysis of all described PT21 cases (from 1973 to 2015). In particular, we built an integrated, comparative map from 125 cases with or without DS fulfilling stringent cytogenetic and clinical criteria. The map allowed to define or exclude as candidates for DS fine Hsa21 sequence intervals, also integrating duplication copy number variants (CNVs) data. A highly restricted DSCR (HR-DSCR) of only 34 kb on distal 21q22.13 has been identified as the minimal region whose duplication is shared by all DS subjects and is absent in all non-DS subjects. Also being spared by any duplication CNV in healthy subjects, HR-DSCR is proposed as a candidate for the typical DS features, the intellectual disability and some facial phenotypes. HR-DSCR contains no known gene and has relevant homology only to the chimpanzee genome. Searching for HR-DSCR functional loci might become a priority for understanding the fundamental genotype-phenotype relationships in DS.

A quantitative transcriptome reference map of the normal human brain.

We performed an innovative systematic meta-analysis of 60 gene expression profiles of whole normal human brain, to provide a quantitative transcriptome reference map of it, i.e. a reference typical value of expression for each of the 39,250 known, mapped and 26,026 uncharacterized (unmapped) transcripts. To this aim, we used the software named Transcriptome Mapper (TRAM), which is able to generate transcriptome maps based on gene expression data from multiple sources. We also analyzed differential expression by comparing the brain transcriptome with those derived from human foetal brain gene expression, from a pool of human tissues (except the brain) and from the two normal human brain regions cerebellum and cerebral cortex, which are two of the main regions severely affected when cognitive impairment occurs, as happens in the case of trisomy 21. Data were downloaded from microarray databases, processed and analyzed using TRAM software and validated in vitro by assaying gene expression through several magnitude orders by ‘real-time’ reverse transcription polymerase chain reaction (RT-PCR). The excellent agreement between in silico and experimental data suggested that our transcriptome maps may be a useful quantitative reference benchmark for gene expression studies related to the human brain. Furthermore, our analysis yielded biological insights about those genes which have an intrinsic over-/under-expression in the brain, in addition offering a basis for the regional analysis of gene expression. This could be useful for the study of chromosomal alterations associated to cognitive impairment, such as trisomy 21, the most common genetic cause of intellectual disability.

Integrated differential transcriptome maps of Acute Megakaryoblastic Leukemia (AMKL) in children with or without Down Syndrome (DS).

BackgroundThe incidence of Acute Megakaryoblastic Leukemia (AMKL) is 500-fold higher in children with Down Syndrome (DS) compared with non-DS children, but the relevance of trisomy 21 as a specific background of AMKL in DS is still an open issue. Several Authors have determined gene expression profiles by microarray analysis in DS and/or non-DS AMKL. Due to the rarity of AMKL, these studies were typically limited to a small group of samples.MethodsWe generated integrated quantitative transcriptome maps by systematic meta-analysis from any available gene expression profile dataset related to AMKL in pediatric age. This task has been accomplished using a tool recently described by us for the generation and the analysis of quantitative transcriptome maps, TRAM (Transcriptome Mapper), which allows effective integration of data obtained from different experimenters, experimental platforms and data sources. This allowed us to explore gene expression changes involved in transition from normal megakaryocytes (MK, n=19) to DS (n=43) or non-DS (n=45) AMKL blasts, including the analysis of Transient Myeloproliferative Disorder (TMD, n=20), a pre-leukemia condition.ResultsWe propose a biological model of the transcriptome depicting progressive changes from MK to TMD and then to DS AMKL. The data indicate the repression of genes involved in MK differentiation, in particular the cluster on chromosome 4 including PF4 (platelet factor 4) and PPBP (pro-platelet basic protein); the gene for the mitogen-activated protein kinase MAP3K10 and the thrombopoietin receptor gene MPL. Moreover, comparing both DS and non-DS AMKL with MK, we identified three potential clinical markers of progression to AMKL: TMEM241 (transmembrane protein 241) was the most over-expressed single gene, while APOC2 (apolipoprotein C-II) and ZNF587B (zinc finger protein 587B) appear to be the most discriminant markers of progression, specifically to DS AMKL. Finally, the chromosome 21 (chr21) genes resulted to be the most over-expressed in DS and non-DS AMKL, as well as in TMD, pointing out a key role of chr21 genes in differentiating AMKL from MK.ConclusionsOur study presents an integrated original model of the DS AMLK transcriptome, providing the identification of genes relevant for its pathophysiology which can potentially be new clinical markers.

A quantitative transcriptome reference map of the normal human hippocampus

We performed an innovative systematic meta‐analysis of 41 gene expression profiles of normal human hippocampus to provide a quantitative transcriptome reference map of it, i.e. a reference typical value of expression for each of the 30,739 known mapped and the 16,258 uncharacterized (unmapped) transcripts. For this aim, we used the software called TRAM (Transcriptome Mapper), which is able to generate transcriptome maps based on gene expression data from multiple sources. We also analyzed differential expression by comparing the hippocampus with the whole brain transcriptome map to identify a typical expression pattern of this subregion compared with the whole organ. Finally, due to the fact that the hippocampus is one of the main brain region to be severely affected in trisomy 21 (the best known genetic cause of intellectual disability), a particular attention was paid to the expression of chromosome 21 (chr21) genes. Data were downloaded from microarray databases, processed, and analyzed using TRAM software. Among the main findings, the most over‐expressed loci in the hippocampus are the expressed sequence tag cluster Hs.732685 and the member of the calmodulin gene family CALM2. The tubulin folding cofactor B (TBCB) gene is the best gene at behaving like a housekeeping gene. The hippocampus vs. the whole brain differential transcriptome map shows the over‐expression of LINC00114, a long non‐coding RNA mapped on chr21. The hippocampus transcriptome map was validated in vitro by assaying gene expression through several magnitude orders by “Real‐Time” reverse transcription polymerase chain reaction (RT‐PCR). The highly significant agreement between in silico and experimental data suggested that our transcriptome map may be a useful quantitative reference benchmark for gene expression studies related to human hippocampus. Furthermore, our analysis yielded biological insights about those genes that have an intrinsic over‐/under‐expression in the hippocampus.

Identification of minimal eukaryotic introns through GeneBase, a user-friendly tool for parsing the NCBI Gene databank

We have developed GeneBase, a full parser of the National Center for Biotechnology Information (NCBI) Gene database, which generates a fully structured local database with an intuitive user-friendly graphic interface for personal computers. Features of all the annotated eukaryotic genes are accessible through three main software tables, including for each entry details such as the gene summary, the gene exon/intron structure and the specific Gene Ontology attributions. The structuring of the data, the creation of additional calculation fields and the integration with nucleotide sequences allow users to make many types of comparisons and calculations that are useful for data retrieval and analysis. We provide an original example analysis of the existing introns across all the available species, through which the classic biological problem of the ‘minimal intron’ may find a solution using available data. Based on all currently available data, we can define the shortest known eukaryotic GT-AG intron length, setting the physical limit at the 30 base pair intron belonging to the human MST1L gene. This ‘model intron’ will shed light on the minimal requirement elements of recognition used for conventional splicing functioning. Remarkably, this size is indeed consistent with the sum of the splicing consensus sequence lengths.

Integrated Transcriptome Map Highlights Structural and Functional Aspects of the Normal Human Heart

A systematic meta‐analysis of the available gene expression profiling datasets for the whole normal human heart generated a quantitative transcriptome reference map of this organ. Transcriptome Mapper (TRAM) software integrated 32 gene expression profile datasets from different sources returning a reference value of expression for each of the 43,360 known, mapped transcripts assayed by any of the experimental platforms used in this regard. Main findings include the visualization at the gene and chromosomal levels of the classical description of the basic histology and physiology of the heart, the identification of suitable housekeeping reference genes, the analysis of stoichiometry of gene products, and the focusing on chromosome 21 genes, which are present in one excess copy in Down syndrome subjects, presenting cardiovascular defects in 30–40% of cases. Independent in vitro validation showed an excellent correlation coefficient (r = 0.98) with the in silico data. Remarkably, heart/non‐cardiac tissue expression ratio may also be used to anticipate that effects of mutations will most probably affect or not the heart. The quantitative reference global portrait of gene expression in the whole normal human heart illustrates the structural and functional aspects of the whole organ and is a general model to understand the mechanisms underlying heart pathophysiology. J. Cell. Physiol. 232: 759–770, 2017.

GeneBase 1.1: a tool to summarize data from NCBI gene datasets and its application to an update of human gene statistics

We release GeneBase 1.1, a local tool with a graphical interface useful for parsing, structuring and indexing data from the National Center for Biotechnology Information (NCBI) Gene data bank. Compared to its predecessor GeneBase (1.0), GeneBase 1.1 now allows dynamic calculation and summarization in terms of median, mean, standard deviation and total for many quantitative parameters associated with genes, gene transcripts and gene features (exons, introns, coding sequences, untranslated regions). GeneBase 1.1 thus offers the opportunity to perform analyses of the main gene structure parameters also following the search for any set of genes with the desired characteristics, allowing unique functionalities not provided by the NCBI Gene itself. In order to show the potential of our tool for local parsing, structuring and dynamic summarizing of publicly available databases for data retrieval, analysis and testing of biological hypotheses, we provide as a sample application a revised set of statistics for human nuclear genes, gene transcripts and gene features. In contrast with previous estimations strongly underestimating the length of human genes, a ‘mean’ human protein-coding gene is 67 kbp long, has eleven 309 bp long exons and ten 6355 bp long introns. Median, mean and extreme values are provided for many other features offering an updated reference source for human genome studies, data useful to set parameters for bioinformatic tools and interesting clues to the biomedical meaning of the gene features themselves. Database URL: http://apollo11.isto.unibo.it/software/

Characterization of human gene locus CYYR1: a complex multi-transcript system.

Cysteine/tyrosine-rich 1 (CYYR1) is a gene we previously identified on human chromosome 21 starting from an in-depth bioinformatics analysis of chromosome 21 segment 40/105 (21q21.3), where no coding region had previously been predicted. CYYR1 was initially characterized as a four-exon gene, whose brain-derived cDNA sequencing predicts a 154-amino acid product. In this study we provide, with in silico and in vitro analyses, the first detailed description of the human CYYR1 locus. The analysis of this locus revealed that it is composed of a multi-transcript system, which includes at least seven CYYR1 alternative spliced isoforms and a new CYYR1 antisense gene (named CYYR1-AS1). In particular, we cloned, for the first time, the following isoforms: CYYR1-1,2,3,4b and CYYR1-1,2,3b, which present a different 3′ transcribed region, with a consequent different carboxy-terminus of the predicted proteins; CYYR1-1,2,4 lacks exon 3; CYYR1-1,2,2bis,3,4 presents an additional exon between exon 2 and exon 3; CYYR1-1b,2,3,4 presents a different 5′ untranslated region when compared to CYYR1. The complexity of the locus is enriched by the presence of an antisense transcript. We have cloned a long transcript overlapping with CYYR1 as an antisense RNA, probably a non-coding RNA. Expression analysis performed in different normal tissues, tumour cell lines as well as in trisomy 21 and euploid fibroblasts has confirmed a quantitative and qualitative variability in the expression pattern of the multi-transcript locus, suggesting a possible role in complex diseases that should be further investigated.