Search

Down syndrome (DS) is characterised by several clinical features including intellectual disability (ID) and craniofacial dysmorphisms. In 1976, Jackson and coll. identified a checklist of signs for clinical diagnosis of DS; the utility of these checklists in improving the accuracy of clinical diagnosis has been recently reaffirmed, but they have rarely been revised. The purpose of this work is to reassess the characteristic phenotypic signs and their frequencies in 233 DS subjects, following Jackson's checklist. 63.77% of the subjects showed more than 12 signs while none showed less than 5, confirming the effectiveness of Jackson's checklist for the clinical diagnosis of DS. An association between three phenotypic signs emerged, allowing us to distinguish two sub-phenotypes: Brachycephaly, short and broad Hands, short Neck (BHN), which is more frequent, and "non-BHN". The strong association of these signs might be interpreted in the context of the growth defects observed in DS children suggesting decreased cell proliferation. Lastly, cognitive assessments were investigated for 114 subjects. The lack of association between the presence of a physical sign or the number of signs present in a subject and cognitive skills disproves the stereotype that physical characteristics are predictive of degree of ID.

Background Trisomy 21 (T21) is a genetic alteration characterised by the presence of an extra full or partial human chromosome 21 (Hsa21) leading to Down syndrome (DS), the most common form of intellectual disability (ID). It is broadly agreed that the presence of extra genetic material in T21 gives origin to an altered expression of genes located on Hsa21 leading to DS phenotype. The aim of this study was to analyse T21 and normal control blood cell gene expression profiles obtained by total RNA sequencing (RNA-Seq). Results The results were elaborated by the TRAM (Transcriptome Mapper) software which generated a differential transcriptome map between human T21 and normal control blood cells providing the gene expression ratios for 17,867 loci. The obtained gene expression profiles were validated through real-time reverse transcription polymerase chain reaction (RT-PCR) assay and compared with previously published data. A post-analysis through transcriptome mapping allowed the identification of the segmental (regional) variation of the expression level across the whole genome (segment-based analysis of expression). Interestingly, the most over-expressed genes encode for interferon-induced proteins, two of them (MX1 and MX2 genes) mapping on Hsa21 (21q22.3). The altered expression of genes involved in mitochondrial translation and energy production also emerged, followed by the altered expression of genes encoding for the folate cycle enzyme, GART, and the folate transporter, SLC19A1. Conclusions The alteration of these pathways might be linked and involved in the manifestation of ID in DS.

Down Syndrome (DS) is the most common genetic alteration responsible for intellectual disability, which refers to deficits in both intellectual and adaptive functioning. According to this, individuals with Down Syndrome (DS) reach developmental milestones (e.g., sitting, walking, and babbling) in the same order as their typically developing peers, but later in life. Since developmental milestones are the first blocks on which development builds, the aims of the current study are to: (i) expand the knowledge of developmental milestone acquisition; and (ii) explore the relationship between developmental milestone acquisition and later development. For this purpose 105 children/adolescents with DS were involved in this study, divided in two groups, Preschoolers (n = 39) and School-age participants (n = 66). Information on the age of acquisition of Sitting, Walking, Babbling, and Sphincter Control was collected, together with cognitive, motor, and adaptive functioning. Sitting predicted later motor development, but, with age, it became less important in predicting motor development in everyday life. Babbling predicted later language development in older children. Finally, Sphincter Control emerged as the strongest predictor of motor, cognitive, language, and adaptive skills, with its role being more evident with increasing age. Our data suggest that the age of reaching the milestones considered in the study has an influence on successive development, a role that can be due to common neural substrates, the environment, and the developmental cascade effect.

This work investigates the role of metabolite levels in the intellectual impairment of subjects with Down syndrome (DS). Homocysteine, folate, vitamin B12, uric acid (UA), creatinine levels and MTHFR C677T genotype were analyzed in 147 subjects with DS. For 77 subjects, metabolite levels were correlated with cognitive tests. Griffiths-III test was administered to 28 subjects (3.08–6.16 years) and WPPSI-III test was administered to 49 subjects (7.08–16.08 years). Significant correlations were found among some metabolite levels and between homocysteine levels and MTHFR C677T genotype. Moreover, homocysteine, UA and creatinine levels resulted increased with age. We did not find any correlation between metabolites and cognitive test score in the younger group. Homocysteine showed statistically significant correlation with WPPSI-III subtest scores when its level is ≥ 7.35 µmol/L, remaining correlated in higher thresholds only for non-verbal area scores. Vitamin B12 showed correlations with all WPPSI-III subtest scores when its level is

Trisomy 21 (Down syndrome, DS) is the main human genetic cause of intellectual disability (ID). Lejeune hypothesized that DS could be considered a metabolic disease, and we found that subjects with DS have a specific plasma and urinary metabolomic profile. In this work we confirmed the alteration of mitochondrial metabolism in DS and also investigated if metabolite levels are related to cognitive aspects of DS. We analyzed the metabolomic profiles of plasma samples from 129 subjects with DS and 46 healthy control (CTRL) subjects by 1H Nuclear Magnetic Resonance (NMR). Multivariate analysis of the NMR metabolomic profiles showed a clear discrimination (up to 94% accuracy) between the two groups. The univariate analysis revealed a significant alteration in 7 metabolites out of 28 assigned unambiguously. Correlations among the metabolite levels in DS and CTRL groups were separately investigated and statistically significant relationships appeared. On the contrary, statistically significant correlations among the NMR-detectable part of DS plasma metabolome and the different intelligence quotient ranges obtained by Griffiths-III or WPPSI-III tests were not found. Even if metabolic imbalance provides a clear discrimination between DS and CTRL groups, it appears that the investigated metabolomic profiles cannot be associated with the degree of ID.

Among Down syndrome (DS) children, 40-50% have congenital heart disease (CHD). Although trisomy 21 is not sufficient to cause CHD, three copies of at least part of chromosome 21 (Hsa21) increases the risk for CHD. In order to establish a genotype-phenotype correlation for CHD in DS, we built an integrated Hsa21 map of all described partial trisomy 21 (PT21) cases with sufficient indications regarding presence or absence of CHD (n=107), focusing on DS PT21 cases. We suggest a DS CHD candidate region on 21q22.2 (0.96Mb), being shared by most PT21 cases with CHD and containing three known protein-coding genes (DSCAM, BACE2, PLAC4) and four known non-coding RNAs (DSCAM-AS1, DSCAM-IT1, LINC00323, MIR3197). The characterization of a DS CHD candidate region provides a useful approach to identify specific genes contributing to the pathology and to orient further investigations and possibly more effective therapy in relation to the multifactorial pathogenesis of CHD.

This article contains further data and information from our published manuscript [1]. We aim to identify significant transcriptome alterations of total normal human thyroid vs. histologically normal thyroid adjacent to papillary thyroid carcinoma. We performed a systematic meta-analysis of all the available gene expression profiles for the whole organ also collecting gene expression data for the normal thyroid adjacent to papillary thyroid carcinoma. A differential quantitative transcriptome reference map was generated by using TRAM (Transcriptome Mapper) software able to combine, normalize and integrate a total of 35 datasets from total normal thyroid and 40 datasets from histologically normal thyroid adjacent to papillary thyroid carcinoma from different sources. This analysis identified genes and genome segments that significantly discriminated the two groups of samples. Differentially expressed genes were grouped and enrichment function analyses were performed identifying the main features of the differentially expressed genes between total normal thyroid and histologically normal thyroid adjacent to papillary thyroid carcinoma. The search for housekeeping genes retrieved 414 loci.

Objective A well-known limit of genome browsers is that the large amount of genome and gene data is not organized in the form of a searchable database, hampering full management of numerical data and free calculations. Due to the continuous increase of data deposited in genomic repositories, their content revision and analysis is recommended. Using GeneBase, a software with a graphical interface able to import and elaborate National Center for Biotechnology Information (NCBI) Gene database entries, we provide tabulated spreadsheets updated to 2019 about human nuclear protein-coding gene data set ready to be used for any type of analysis about genes, transcripts and gene organization. Results Comparison with previous reports reveals substantial change in the number of known nuclear protein-coding genes (now 19,116), the protein-coding non-redundant transcriptome space [now 59,281,518 base pair (bp), 10.1% increase], the number of exons (now 562,164, 36.2% increase) due to a relevant increase of the RNA isoforms recorded. Other parameters such as gene, exon or intron mean and extreme length appear to have reached a stability that is unlikely to be substantially modified by human genome data updates, at least regarding protein-coding genes. Finally, we confirm that there are no human introns shorter than 30 bp.

Objective Basic parameters commonly used to describe genomes including length, weight and relative guanine-cytosine (GC) content are widely cited in absence of a primary source. By using updated data and original software we determined these values to the best of our knowledge as standard reference for the whole human nuclear genome, for each chromosome and for mitochondrial DNA. We also devised a method to calculate the relative GC content in the whole messenger RNA sequence set and in transcriptomes by multiplying the GC content of each gene by its mean expression level. Results The male nuclear diploid genome extends for 6.27 Gigabase pairs (Gbp), is 205.00 cm (cm) long and weighs 6.41 picograms (pg). Female values are 6.37 Gbp, 208.23 cm, 6.51 pg. The individual variability and the implication for the DNA informational density in terms of bits/volume were discussed. The genomic GC content is 40.9%. Following analysis in different transcriptomes and species, we showed that the greatest deviation was observed in the pathological condition analysed (trisomy 21 leukaemic cells) and in Caenorhabditis elegans. Our results may represent a solid basis for further investigation on human structural and functional genomics while also providing a framework for other genome comparative analysis. Electronic supplementary material The online version of this article (10.1186/s13104-019-4137-z) contains supplementary material, which is available to authorized users.

Background: 5,10-Methylentetrahydrofolate reductase (MTHFR) C677T polymorphism is one of the most studied genetic variations in the human genome. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) is one of the most used techniques to characterize the point mutations in genomic sequences because of its suitability and low cost. The most widely used method for the MTHFR C677T polymorphism characterization was developed by Frosst et al. (1995) but appears to have some technical limitations. The aim of this study was to propose a novel PCR-RFLP method for the detection of this polymorphism. Methods: In order to retrieve all published articles possibly describing any PCR-RFLP methods useful to analyze MTHFR C677T polymorphism, we performed systematic queries on PubMed, using a combination of Boolean operators (AND/OR) and MeSH terms. Amplify software was used in order to design a new primer pair following the optimal standard criteria. Primer-BLAST software was used to check primer pair's biological specificity. Results: The analysis of previous literature showed that PCR-RFLP method remains the most used technique. None of the 108 primer pairs described was ideal with regard to main accepted primer pair biochemical technical parameters. The new primer pair amplifies a DNA-fragment of 513 base pair (bp) that, in the presence of the polymorphism, is cut by HinfI enzyme in two pieces of 146bp and 367bp and clearly visible on 2% agarose gel. The level of expertise and the materials required are minimal and the protocol takes one day to carry out. Conclusion: Our original PCR-RFLP strategy, specifically designed to make the analysis optimal with respect to PCR primers and gel analysis, fits the ideal criteria compared to the widely used strategy by Frosst et al (1995) as well as any other PCR-RFLP strategies proposed for MTHFR C677T polymorphism genotyping to date.

Down Syndrome (DS) is the most frequent form of intellectual disability (ID) of genetic origin, whose main features include craniofacial dysmorphisms and cardiovascular defects. In 1959, Lejeune and coll. described an extra copy of chromosome 21 (Hsa21) in children with DS (trisomy 21, or T21). We first review how different biological mechanisms may lead to the gain of genetic material of Hsa21 in the cells, originating from different combinations of genetic conditions, including a free or translocated extra copy of Hsa21, distributed in all cells or only in a part of them (mosaicism), with a complete or partial representation of the Hsa21 long arm (21q). Although it is broadly agreed that the DS phenotype originates from the altered expression of the genes located on Hsa21, its molecular pathogenesis is still unknown. We therefore illustrate how recent genomic science may be useful in the elucidation of the genotype-phenotype relationship in DS.

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.

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. © 2016 Wiley Periodicals, Inc.

The incomplete determination of the mRNA 5' end sequence may lead to the incorrect assignment of the first AUG codon and to errors in the prediction of the encoded protein product. Due to the significance of the mouse as a model organism in biomedical research, we performed a systematic identification of coding regions at the 5' end of all known mouse mRNAs, using an automated expressed sequence tag (EST)-based approach which we have previously described. By parsing almost 4 million BLAT alignments we found 351 mouse loci, out of 20,221 analyzed, in which an extension of the mRNA 5' coding region was identified. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for Apc2 and Mknk2 cDNAs. We also generated a list of 16,330 mouse mRNAs where the presence of an in-frame stop codon upstream of the known start codon indicates completeness of the coding sequence at 5' end in the current form. Systematic searches in the main mouse genome databases and genome browsers showed that 82 % of our results are original and have not been identified by their annotation pipelines. Moreover, the same information is not easily derivable from RNA-Seq data, due to short sequence length and laboriousness in building full-length transcript structures. In conclusion, our results improve the determination of full-length 5' coding sequences and might be useful in order to reduce errors when studying mouse gene structure and function in biomedical research.

Trisomy 21 causes Down syndrome (DS), the most common human genetic disorder and the leading genetic cause of intellectual disability. The alteration of one-carbon metabolism was described as the possible metabolic cause of the intellectual disability development in subjects with DS. One of the biochemical pathways involved in the one-carbon group transfer is the folate cycle. The cytotoxic drug methotrexate (MTX) is a folic acid (FA) analogue which inhibits the activity of dihydrofolate reductase enzyme involved in the one-carbon metabolic cycle. Trisomy 21 cells are more sensitive to the MTX effect than euploid cells, and in 1986 Jerome Lejeune and Coll. demonstrated that MTX was twice as toxic in trisomy 21 lymphocytes than in control cells. In the present work, the rescue effect on MTX toxicity mediated by FA and some of its derivatives, tetrahydrofolate (THF), 5-formyl-THF, and 5-methyl-THF, in both normal and trisomy 21 skin fibroblast cells, was evaluated. A statistically significant rescue effect was obtained by 5-formyl-THF, 5-methyl-THF, and their combination, administered together with MTX. In conclusion, trisomy 21 fibroblast cell lines showed a good response to the rescue effects of 5-formyl-THF and 5-methyl-THF on the MTX toxicity almost as normal cell lines.

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. © 2015 Wiley Periodicals, Inc.

Background The thyroid is the earliest endocrine structure to appear during human development, and thyroid hormones are necessary for proper organism development, in particular for the nervous system and heart, normal growth and skeletal maturation. To date a quantitative, validated transcriptional atlas of the whole normal human thyroid does not exist and the availability of a detailed expression map might be an excellent occasion to investigate the many features of the thyroid transcriptome. Results We present a view at the molecular level of the normal human thyroid histology and physiology obtained by a systematic meta-analysis of all the available gene expression profiles for the whole organ. A quantitative transcriptome reference map was generated by using the TRAM (Transcriptome Mapper) software able to combine, normalize and integrate a total of 35 suitable datasets from different sources thus providing a typical reference expression value for each of the 27,275 known, mapped transcripts obtained. The experimental in vitro validation of data was performed by “Real-Time” reverse transcription polymerase chain reaction showing an excellent correlation coefficient (r = 0.93) with data obtained in silico. Conclusions Our study provides a quantitative global reference portrait of gene expression in the normal human thyroid and highlights differential expression between normal human thyroid and a pool of non-thyroid tissues useful for modeling correlations between thyroidal gene expression and specific thyroid functions and diseases. The experimental in vitro validation supports the possible usefulness of the human thyroid transcriptome map as a reference for molecular studies of the physiology and pathology of this organ. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-4049-z) contains supplementary material, which is available to authorized users.

Down syndrome (DS) is due to the presence of an extra full or partial chromosome 21 (Hsa21). The identification of genes contributing to DS pathogenesis could be the key to any rational therapy of the associated intellectual disability. We aim at generating quantitative transcriptome maps in DS integrating all gene expression profile datasets available for any cell type or tissue, to obtain a complete model of the transcriptome in terms of both expression values for each gene and segmental trend of gene expression along each chromosome. We used the TRAM (Transcriptome Mapper) software for this meta-analysis, comparing transcript expression levels and profiles between DS and normal brain, lymphoblastoid cell lines, blood cells, fibroblasts, thymus and induced pluripotent stem cells, respectively. TRAM combined, normalized, and integrated datasets from different sources and across diverse experimental platforms. The main output was a linear expression value that may be used as a reference for each of up to 37,181 mapped transcripts analyzed, related to both known genes and expression sequence tag (EST) clusters. An independent example in vitro validation of fibroblast transcriptome map data was performed through “Real-Time” reverse transcription polymerase chain reaction showing an excellent correlation coefficient (r = 0.93, p < 0.0001) with data obtained in silico. The availability of linear expression values for each gene allowed the testing of the gene dosage hypothesis of the expected 3:2 DS/normal ratio for Hsa21 as well as other human genes in DS, in addition to listing genes differentially expressed with statistical significance. Although a fraction of Hsa21 genes escapes dosage effects, Hsa21 genes are selectively over-expressed in DS samples compared to genes from other chromosomes, reflecting a decisive role in the pathogenesis of the syndrome. Finally, the analysis of chromosomal segments reveals a high prevalence of Hsa21 over-expressed segments over the other genomic regions, suggesting, in particular, a specific region on Hsa21 that appears to be frequently over-expressed (21q22). Our complete datasets are released as a new framework to investigate transcription in DS for individual genes as well as chromosomal segments in different cell types and tissues.

The ideal reference, or control, gene for the study of gene expression in a given organism should be expressed at a medium‑high level for easy detection, should be expressed at a constant/stable level throughout different cell types and within the same cell type undergoing different treatments, and should maintain these features through as many different tissues of the organism. From a biological point of view, these theoretical requirements of an ideal reference gene appear to be best suited to housekeeping (HK) genes. Recent advancements in the quality and completeness of human expression microarray data and in their statistical analysis may provide new clues toward the quantitative standardization of human gene expression studies in biology and medicine, both cross‑ and within‑tissue. The systematic approach used by the present study is based on the Transcriptome Mapper tool and exploits the automated reassignment of probes to corresponding genes, intra‑ and inter‑sample normalization, elaboration and representation of gene expression values in linear form within an indexed and searchable database with a graphical interface recording quantitative levels of expression, expression variability and cross‑tissue width of expression for more than 31,000 transcripts. The present study conducted a meta‑analysis of a pool of 646 expression profile data sets from 54 different human tissues and identified actin γ 1 as the HK gene that best fits the combination of all the traditional criteria to be used as a reference gene for general use; two ribosomal protein genes, RPS18 and RPS27, and one aquaporin gene, POM121 transmembrane nucleporin C, were also identified. The present study provided a list of tissue‑ and organ‑specific genes that may be most suited for the following individual tissues/organs: Adipose tissue, bone marrow, brain, heart, kidney, liver, lung, ovary, skeletal muscle and testis; and also provides in these cases a representative, quantitative portrait of the relative, typical gene‑expression profile in the form of searchable database tables.

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. © 2016 Wiley Periodicals, Inc.

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.

The known difficulty in obtaining the actual full length, complete sequence of a messenger RNA (mRNA) may lead to the erroneous determination of its coding sequence at the 5' region (5' end mRNA artifact), and consequently to the wrong assignment of the translation start codon, leading to the inaccurate prediction of the encoded polypeptide at its amino terminus. Among the known human genes whose study was affected by this artifact, we can include disco interacting protein 2 homolog A (DIP2A; KIAA0184), Down syndrome critical region 1 (DSCR1), SON DNA binding protein (SON), trefoil factor 3 (TFF3) and URB1 ribosome biogenesis 1 homolog (URB1; KIAA0539) on chromosome 21, as well as receptor for activated C kinase 1 (RACK1, also known as GNB2L1), glutaminyl‑tRNA synthetase (QARS) and tyrosyl-DNA phosphodiesterase 2 (TDP2) along with another 474 loci, including interleukin 16 (IL16). In this review, we discuss the causes of this issue, its quantitative incidence in biomedical research, the consequences in biology and medicine, and the possible solutions for obtaining the actual amino acid sequence of proteins in the post-genomics era.

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/

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.

The term "5´ end mRNA artifact" refers to the incorrect assignment of the first AUG codon in an mRNA, due to the incomplete determination of its 5´ end sequence (Casadei et al., 2003). Since the '70s, the amino acid sequence of gene products has been routinely deduced from the nucleotide sequence of the relative cloned cDNA (DNA complementary to mRNA), according to rules for recognition of the start codon (first-AUG rule, optimal sequence context) and the genetic code (Kozak, 2002). All standard methods for the cloning of cDNA are affected by a potential inability to effectively clone the 5´ region of mRNA. This is due to the reverse transcriptase failure to extend first-strand cDNA along the full length of the mRNA template toward its 5´ end (Sambrook, 2001). The identification of a more complete mRNA 5´ end could reveal an additional upstream AUG – in-frame with the previously determined one – thus extending the predicted amino terminus sequence of the product and avoiding subsequent relevant errors in the experimental study of the relative cDNA (Casadei et al., 2003). The continuous incorporation of information derived from individual and large-scale cDNA sequencing projects, including those specifically designed to characterize mRNA 5´ end (Carninci et al., 1996; Suzuki et al., 2000; Porcel et al., 2004), in the last few years led to continuous improvement of completeness of mRNA reference sequences (e.g., RefSeq), and also to the corresponding protein coding sequences. However, genome browsers do not appear to systematically extract useful information from the ever-increasing vast quantity of EST (expressed sequence tag) data. To date, EST data remain invaluable due to significantly longer continuous RNA sequences they may provide in comparison with the very short fragments typically deposited in current high-throughput nucleotide sequencing databases. We previously used individual EST-based gene model refinement by classic in silico sequence analysis to revise the mRNA sequence of 109 human chromosome 21 protein-coding genes (Casadei et al., 2003). The success of this approach encouraged us to develop a piece of software ("5'_ORF_Extender" software) in order to automate the steps that were previously performed manually, applying it to the Danio rerio (zebrafish) genome (Frabetti et al., 2007). In the present work, we present a modified strategy able to analyze the much more numerous human sequences. Firstly, we fully revised the software algorithm by using pre-computed coordinates of the UCSC-downloaded RefSeqs and ESTs genome alignment data and specific UCSC-downloaded EST sequence entries. Furthermore, we adopted an original quality filter which was able to test if each single EST candidate with sequence information of possible use for extending a known mRNA, was attributed to the same locus of that mRNA by an updated, complete and embedded version of UniGene. Lastly, we automated data summarization for an analyzed genome. Following these improvements, parsing more than 7 million BLAT alignment, 5'_ORF_Extender 2.0 recognized a total of 477 loci, out of the 18,665 human loci represented in the mRNA reference set, as bona fide candidates for extension. Proof-of-concept confirmation was obtained by in vitro cloning and sequencing for GNB2L1 (guanine nucleotide binding protein (G protein), beta polypeptide 2-like 1), QARS (glutaminyl-tRNA synthetase) and TDP2 (tyrosyl-DNA phosphodiesterase 2) cDNAs, and the consequences for the functional studies of these loci are discussed. In addition, we generated a list of 20,775 human mRNAs in which the presence of an in-frame stop codon upstream of the known start codon indicates completeness of the coding sequence at 5´ in the current form. Bibliografia: R. Casadei et al., mRNA 5’ region sequence incompleteness: a potential source of systematic errors in translation initiation codon assignment in human mRNAs, Gene 321 (2003) 185–193. M. Kozak, Pushing the limits of the scanning mechanism for initiation of translation, Gene 99 (2002) 1–34. P. Carninci et al., High-efficiency fulllength cDNA cloning by biotinylated CAP trapper, Genomics 37 (1996) 327–336. F. Frabetti et al., Systematic analysis of mRNA 5’ coding sequence incompleteness in Danio rerio: an automated EST-based approach, Biol. Direct 2 (2007) 34.

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.

Down Syndrome (DS) is the most frequent human chromosomal disorder. Main symptoms include intellectual disability (ID), cardiovascular defects and craniofacial dysmorphisms. Despite ID being measured by a test of symbolic logic skills, it is common for children with DS to arouse a climate of affective intensity greater than the norm. In 1959, Jérôme Lejeune (1926-1994) and coll. described an additional chromosome 21 (Hsa21) in children with DS (trisomy 21), giving origin to the field of medical genetics. Remarkably, the discovery of trisomy 21 had relevant social consequences for the affected children, in that their parents were no longer suspected to be alcoholics or infected with syphilis. Although it is broadly agreed that the DS phenotype originates from the altered expression of the genes located on Hsa21, its molecular pathogenesis is still unknown. To date, no therapy is recognized and recommended by guidelines as being effective in improving the cognitive abilities of persons with DS. The aim of this article is to categorize main therapeutical approaches or pathways to new approaches reported in the biomedical literature, to extract critical methodological points from the works of Lejeune and then to propose a new research project aimed to generate and integrate clinical, biochemical, genetic and bioinformatic data in order to identify novel therapeutic targets for this form of trisomy. We show here that nearly all the current lines of research were pursued, theorized or foreseen by Lejeune, and that central points of his method remain current: positive hypothesis about the existence of a solution, envision of systematic investigation of cell machinery, anchoring of clinical and biochemical finding to the chromosome physical map, and continuing clinical observation of the affected children. We therefore propose a project aimed at producing both experimentally and by meta-analysis state-of-the-art maps and databases related to clinical/phenotype, cytogenetics, exome, transcriptome, methylome, molecular biology, metabolome and mutations data. The primary expected outcome of this research project is the identification of a restricted list of strong candidate genes and mechanisms for ID in persons with DS in order to devise new rational therapeutic approaches.

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.

Codon bias is the phenomenon in which distinct synonymous codons are used with different frequencies. We define here the “codonome value” as the total number of codons present across all the expressed mRNAs in a given biological condition. We have developed the “CODONOME” software, which calculates the codon bias and, following integration with a gene expression profile, estimates the actual frequency of each codon at the transcriptome level (codonome bias) of a given tissue. Systematic analysis across different human tissues and multiple species shows a surprisingly tight correlation between the codon bias and the codonome bias. An aneuploidy and cancer condition such as that of Down Syndrome-related acute megakaryoblastic leukemia (DS-AMKL), does not appear to alter this relationship. The law of correlation between codon bias and codonome emerges as a property of the distribution and range of the number, sequence and expression level of the genes in a genome.