Zinc metabolism and its role in immunity status in subjects with trisomy 21: chromosomal dosage effect

Introduction

Trisomy 21 (T21), which causes Down syndrome (DS), is the most common chromosomal aneuploidy in humankind and includes different clinical comorbidities, among which the alteration of the immune system has a heavy impact on patient's lives. A molecule with an important role in immune response is zinc and it is known that its concentration is significantly lower in children with T21. Different hypotheses were made about this metabolic alteration and one of the reasons might be the overexpression of superoxide dismutase 1 (SOD1) gene, as zinc is part of the SOD1 active enzymatic center.

Methods

The aim of our work is to explore if there is a linear correlation between zinc level and immune cell levels measured in a total of 217 blood samples from subjects with T21. Furthermore, transcriptome map analyses were performed using Transcriptome Mapper (TRAM) software to investigate whether a difference in gene expression is detectable between subjects with T21 and euploid control group in tissues and cells involved in the immune response such as lymphoblastoid cells, thymus and white blood cells.

Results

Our results have confirmed the literature data stating that the blood zinc level in subjects with T21 is lower compared to the general population; in addition, we report that the T21/control zinc concentration ratio is 2:3, consistent with a chromosomal dosage effect due to the presence of three copies of chromosome 21. The transcriptome map analyses showed an alteration of some gene's expression which might explain low levels of zinc in the blood.

Discussion

Our data suggest that zinc level is not associated with the levels of immunity cells or proteins analyzed themselves and rather the main role of this ion might be played in altering immune cell function.

Machine learning based analysis for intellectual disability in Down syndrome

Down syndrome (DS) or trisomy 21 is the most common genetic cause of intellectual disability (ID), but a pathogenic mechanism has not been identified yet. Studying a complex and not monogenic condition such as DS, a clear correlation between cause and effect might be difficult to find through classical analysis methods, thus different approaches need to be used. The increased availability of big data has made the use of artificial intelligence (AI) and in particular machine learning (ML) in the medical field possible. The purpose of this work is the application of ML techniques to provide an analysis of clinical records obtained from subjects with DS and study their association with ID. We have applied two tree-based ML models (random forest and gradient boosting machine) to the research question: how to identify key features likely associated with ID in DS. We analyzed 109 features (or variables) in 106 DS subjects. The outcome of the analysis was the age equivalent (AE) score as indicator of intellectual functioning, impaired in ID. We applied several methods to configure the models: feature selection through Boruta framework to minimize random correlation; data augmentation to overcome the issue of a small dataset; age effect mitigation to take into account the chronological age of the subjects. The results show that ML algorithms can be applied with good accuracy to identify variables likely involved in cognitive impairment in DS. In particular, we show how random forest and gradient boosting machine produce results with low error (MSE <0.12) and an acceptable R2 (0.70 and 0.93). Interestingly, the ranking of the variables point to several features of interest related to hearing, gastrointestinal alterations, thyroid state, immune system and vitamin B12 that can be considered with particular attention for improving care pathways for people with DS. In conclusion, ML-based model may assist researchers in identifying key features likely correlated with ID in DS, and ultimately, may improve research efforts focused on the identification of possible therapeutic targets and new care pathways. We believe this study can be the basis for further testing/validating of our algorithms with multiple and larger datasets.

Partial trisomy 21 with or without highly restricted Down syndrome critical region (HR-DSCR): report of two new cases and reanalysis of the genotype-phenotype association

BACKGROUND

Down syndrome (DS) is caused by the presence of an extra copy of full or partial human chromosome 21 (Hsa21). Partial (segmental) trisomy 21 (PT21) is the duplication of only a delimited region of Hsa21 and can be associated or not to DS: the study of PT21 cases is an invaluable model for addressing genotype-phenotype correlation in DS. Previous works reported systematic reanalyses of 132 subjects with PT21 and allowed the identification of a 34-kb highly restricted DS critical region (HR-DSCR) as the minimal region whose duplication is shared by all PT21 subjects diagnosed with DS.

METHODS

We report clinical data and cytogenetic analysis of two children with PT21, one with DS and the other without DS. Moreover, we performed a systematic bibliographic search for any new PT21 report.

RESULTS

Clinical and cytogenetic analyses of the two PT21 children have been reported: in Case 1 the duplication involves the whole long arm of Hsa21, except for the last 2.7 Mb, which are deleted as a consequence of an isodicentric 21: the HR-DSCR is within the duplicated regions and the child is diagnosed with DS. In Case 2 the duplication involves 7.1 Mb of distal 21q22, with a deletion of 2.1 Mb of proximal 20p, as a consequence of an unbalanced translocation: the HR-DSCR is not duplicated and the child presents with psychomotor development delay but no clinical signs of DS. Furthermore, two PT21 reports recently published (named Case 3 and 4) have been discussed: Case 3 has DS diagnosis, nearly full trisomy for Hsa21 and a monosomy for the 21q22.3 region. Case 4 is a baby without DS and a 0.56-Mb duplication of 21q22.3. Genotype-phenotype correlation confirmed the presence of three copies of the HR-DSCR in all DS subjects and two copies in all non-DS individuals.

CONCLUSIONS

The results presented here are fully consistent with the hypothesis that the HR-DSCR is critically associated with DS diagnosis. No exception to this pathogenetic model was found. Further studies are needed to detect genetic determinants likely located in the HR-DSCR and possibly responsible for core DS features, in particular intellectual disability.

One-carbon pathway metabolites are altered in the plasma of subjects with Down syndrome: Relation to chromosomal dosage

Introduction

Down syndrome (DS) is the most common chromosomal disorder and it is caused by trisomy of chromosome 21 (Hsa21). Subjects with DS show a large heterogeneity of phenotypes and the most constant clinical features present are typical facies and intellectual disability (ID). Several studies demonstrated that trisomy 21 causes an alteration in the metabolic profile, involving among all the one-carbon cycle.

Methods

We performed enzyme-linked immunosorbent assays (ELISAs) to identify the concentration of 5 different intermediates of the one-carbon cycle in plasma samples obtained from a total of 164 subjects with DS compared to 54 euploid subjects. We investigated: tetrahydrofolate (THF; DS n = 108, control n = 41), 5-methyltetrahydrofolate (5-methyl-THF; DS n = 140, control n = 34), 5-formyltetrahydrofolate (5-formyl-THF; DS n = 80, control n = 21), S-adenosyl-homocysteine (SAH; DS n = 94, control n = 20) and S-adenosyl-methionine (SAM; DS n = 24, control n = 15).

Results

Results highlight specific alterations of THF with a median concentration ratio DS/control of 2:3, a decrease of a necessary molecule perfectly consistent with a chromosomal dosage effect. Moreover, SAM and SAH show a ratio DS/control of 1.82:1 and 3.6:1, respectively.

Discussion

The relevance of these results for the biology of intelligence and its impairment in trisomy 21 is discussed, leading to the final proposal of 5-methyl-THF as the best candidate for a clinical trial aimed at restoring the dysregulation of one-carbon cycle in trisomy 21, possibly improving cognitive skills of subjects with DS.

Structural Characterization of the Highly Restricted Down Syndrome Critical Region on 21q22.13: New KCNJ6 and DSCR4 Transcript Isoforms

Down syndrome (DS) is caused by trisomy of chromosome 21 and it is the most common genetic cause of intellectual disability (ID) in humans. Subjects with DS show a typical phenotype marked by facial dysmorphisms and ID. Partial trisomy 21 (PT21) is a rare genotype characterized by the duplication of a delimited chromosome 21 (Hsa21) portion and it may or may not be associated with DS diagnosis. The highly restricted Down syndrome critical region (HR-DSCR) is a region of Hsa21 present in three copies in all individuals with PT21 and a diagnosis of DS. This region, located on distal 21q22.13, is 34 kbp long and does not include characterized genes. The HR-DSCR is annotated as an intergenic region between KCNJ6-201 transcript encoding for potassium inwardly rectifying channel subfamily J member 6 and DSCR4-201 transcript encoding Down syndrome critical region 4. Two transcripts recently identified by massive RNA-sequencing (RNA-Seq) and automatically annotated on Ensembl database reveal that the HR-DSCR seems to be partially crossed by KCNJ6-202 and DSCR4-202 isoforms. KCNJ6-202 shares the coding sequence with KCNJ6-201 which is involved in many physiological processes, including heart rate in cardiac cells and circuit activity in neuronal cells. DSCR4-202 transcript has the first two exons in common with DSCR4-201, the only experimentally verified gene uniquely present in Hominidae. In this study, we performed in silico and in vitro analyses of the HR-DSCR. Bioinformatic data, obtained using Sequence Read Archive (SRA) and SRA-BLAST software, were confirmed by Reverse Transcription-Polymerase Chain Reaction (RT-PCR) and Sanger sequencing on a panel of human tissues. Our data demonstrate that the HR-DSCR cannot be defined as an intergenic region. Further studies are needed to investigate the functional role of the new transcripts, likely involved in DS phenotypes.

Is the Age of Developmental Milestones a Predictor for Future Development in Down Syndrome?

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.

The transcriptome profile of human trisomy 21 blood cells

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.

One-carbon pathway and cognitive skills in children with Down syndrome

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 < 442 pg/mL. The relevance of the present findings is the detection of a specific metabolite threshold related with a better or worse cognitive score, suggesting that vitamin B12 and homocysteine may have a role in cognitive development in children with DS.

Plasma metabolome and cognitive skills in Down syndrome

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.

Dysregulation of NIPBL leads to impaired RUNX1 expression and haematopoietic defects

The transcription factor RUNX1, a pivotal regulator of HSCs and haematopoiesis, is a frequent target of chromosomal translocations, point mutations or altered gene/protein dosage. These modifications lead or contribute to the development of myelodysplasia, leukaemia or platelet disorders. A better understanding of how regulatory elements contribute to fine‐tune the RUNX1 expression in haematopoietic tissues could improve our knowledge of the mechanisms responsible for normal haematopoiesis and malignancy insurgence. The cohesin RAD21 was reported to be a regulator of RUNX1 expression in the human myeloid HL60 cell line and during primitive haematopoiesis in zebrafish. In our study, we demonstrate that another cohesin, NIPBL, exerts positive regulation of RUNX1 in three different contexts in which RUNX1 displays important functions: in megakaryocytes derived from healthy donors, in bone marrow samples obtained from adult patients with acute myeloid leukaemia and during zebrafish haematopoiesis. In this model, we demonstrate that alterations in the zebrafish orthologue nipblb reduce runx1 expression with consequent defects in its erythroid and myeloid targets such as gata1a and spi1b in an opposite way to rad21. Thus, also in the absence of RUNX1 translocation or mutations, additional factors such as defects in the expression of NIPBL might induce haematological diseases.

Partial trisomy 21 map: Ten cases further supporting the highly restricted Down syndrome critical region (HR-DSCR) on human chromosome 21

BACKGROUND

Down syndrome (DS) is characterized by the presence of an extra full or partial human chromosome 21 (Hsa21). An invaluable model to define genotype-phenotype correlations in DS is the study of the extremely rare cases of partial (segmental) trisomy 21 (PT21), the duplication of only a delimited region of Hsa21 associated or not to DS. A systematic retrospective reanalysis of 125 PT21 cases described up to 2015 allowed the creation of the most comprehensive PT21 map and the identification of a 34-kb highly restricted DS critical region (HR-DSCR) as the minimal region whose duplication is shared by all PT21 subjects diagnosed with DS. We reanalyzed at higher resolution three cases previously published and we accurately searched for any new PT21 reports in order to verify whether HR-DSCR limits could prospectively be confirmed and possibly refined.

METHODS

Hsa21 partial duplications of three PT21 subjects were refined by adding array-based comparative genomic hybridization data. Seven newly described PT21 cases fulfilling stringent cytogenetic and clinical criteria have been incorporated into the PT21 integrated map.

RESULTS

The PT21 map now integrates fine structure of Hsa21 sequence intervals of 132 subjects onto a common framework fully consistent with the presence of a duplicated HR-DSCR, on distal 21q22.13 sub-band, only in DS subjects and not in non-DS individuals. No documented exception to the HR-DSCR model was found.

CONCLUSIONS

The findings presented here further support the association of the HR-DSCR with the diagnosis of DS, representing an unbiased validation of the original model. Further studies are needed to identify and characterize genetic determinants presumably located in the HR-DSCR and functionally associated to the critical manifestations of DS.

Human protein-coding genes and gene feature statistics in 2019

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.

Dataset of differential gene expression between total normal human thyroid and histologically normal thyroid adjacent to papillary thyroid carcinoma

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.

MTHFR C677T polymorphism analysis: A simple, effective restriction enzyme-based method improving previous protocols

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 Hinf I enzyme in two pieces of 146 bp and 367 bp 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.

On the length, weight and GC content of the human genome

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.

Integrated Quantitative Transcriptome Maps of Human Trisomy 21 Tissues and Cells

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.

Plasma and urinary metabolomic profiles of Down syndrome correlate with alteration of mitochondrial metabolism

Down syndrome (DS) is caused by the presence of a supernumerary copy of the human chromosome 21 (Hsa21) and is the most frequent genetic cause of intellectual disability (ID). Key traits of DS are the distinctive facies and cognitive impairment. We conducted for the first time an analysis of the Nuclear Magnetic Resonance (NMR)-detectable part of the metabolome in plasma and urine samples, studying 67 subjects with DS and 29 normal subjects as controls selected among DS siblings. Multivariate analysis of the NMR metabolomic profiles showed a clear discrimination (up to of 80% accuracy) between the DS and the control groups. The univariate analysis of plasma and urine revealed a significant alteration for some interesting metabolites. Remarkably, most of the altered concentrations were consistent with the 3:2 gene dosage model, suggesting effects caused by the presence of three copies of Hsa21 rather than two: DS/normal ratio in plasma was 1.23 (pyruvate), 1.47 (succinate), 1.39 (fumarate), 1.33 (lactate), 1.4 (formate). Several significantly altered metabolites are produced at the beginning or during the Krebs cycle. Accounting for sex, age and fasting state did not significantly affect the main result of both multivariate and univariate analysis.

A molecular view of the normal human thyroid structure and function reconstructed from its reference transcriptome map

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.

Systematic identification of human housekeeping genes possibly useful as references in gene expression studies

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.

Difficulty in obtaining the complete mRNA coding sequence at 5' region (5' end mRNA artifact): Causes, consequences in biology and medicine and possible solutions for obtaining the actual amino acid sequence of proteins (Review)

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.

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/

Meta-Analysis of Parkinson's Disease Transcriptome Data Using TRAM Software: Whole Substantia Nigra Tissue and Single Dopamine Neuron Differential Gene Expression

The understanding of the genetic basis of the Parkinson's disease (PD) and the correlation between genotype and phenotype has revolutionized our knowledge about the pathogenetic mechanisms of neurodegeneration, opening up exciting new therapeutic and neuroprotective perspectives. Genomic knowledge of PD is still in its early stages and can provide a good start for studies of the molecular mechanisms that underlie the gene expression variations and the epigenetic mechanisms that may contribute to the complex and characteristic phenotype of PD. In this study we used the software TRAM (Transcriptome Mapper) to analyse publicly available microarray data of a total of 151 PD patients and 130 healthy controls substantia nigra (SN) samples, to identify chromosomal segments and gene loci differential expression. In particular, we separately analyzed PD patients and controls data from post-mortem snap-frozen SN whole tissue and from laser microdissected midbrain dopamine (DA) neurons, to better characterize the specific DA neuronal expression profile associated with the late-stage Parkinson's condition. The default "Map" mode analysis resulted in 10 significantly over/under-expressed segments, mapping on 8 different chromosomes for SN whole tissue and in 4 segments mapping on 4 different chromosomes for DA neurons. In conclusion, TRAM software allowed us to confirm the deregulation of some genomic regions and loci involved in key molecular pathways related to neurodegeneration, as well as to provide new insights about genes and non-coding RNA transcripts not yet associated with the disease.

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

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.

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.

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.

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

BACKGROUND

The 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.

METHODS

We 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.

RESULTS

We 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.

CONCLUSIONS

Our 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.

An estimation of the number of cells in the human body

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 10(12) and 10(16) 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 × 10(13).

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.

Universal tight correlation of codon bias and pool of RNA codons (codonome): The genome is optimized to allow any distribution of gene expression values in the transcriptome from bacteria to humans

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.

Genome-scale analysis of human mRNA 5' coding sequences based on expressed sequence tag (EST) database

The "5' end mRNA artifact" issue refers to the incorrect assignment of the first AUG codon in an mRNA, due to the incomplete determination of its 5' end sequence. We performed a systematic identification of coding regions at the 5' end of all human known mRNAs, using an automated expressed sequence tag (EST)-based approach. Following parsing of more than 7 million BLAT alignments, we found 477 human loci, out of 18,665 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 GNB2L1, QARS and TDP2 cDNAs, and the consequences for the functional studies of these loci are discussed. We also generated a list of 20,775 human mRNAs where 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.

Complexity of bidirectional transcription and alternative splicing at human RCAN3 locus

Human RCAN3 (regulator of calcineurin 3) belongs to the human RCAN gene family.

In this study we provide, with in silico and in vitro analyses, the first detailed description of the human multi-transcript RCAN3 locus. Its analysis revealed that it is composed of a multigene system that includes at least 21 RCAN3 alternative spliced isoforms (16 of them identified here for the first time) and a new RCAN3 antisense gene (RCAN3AS). In particular, we cloned RCAN3-1,3,4,5 (lacking exon 2), RCAN3-1a,2,3,4,5, RCAN3-1a,3,4,5, RCAN3-1b,2,3,4,5, RCAN3-1c,2,3,4,5, RCAN3-1c,2,4,5 and RCAN3-1c,3,4,5, isoforms that present a different 5′ untranslated region when compared to RCAN3. Moreover, in order to verify the possible 5′ incompleteness of previously identified cDNA isoforms with the reference exon 1, ten more alternative isoforms were retrieved. Bioinformatic searches allowed us to identify RCAN3AS, which overlaps in part with exon 1a, on the opposite strand, for which four different RCAN3AS isoforms were cloned.

In order to analyze the different expression patterns of RCAN3 alternative first exons and of RCAN3AS mRNA isoforms, RT-PCR was performed in 17 human tissues.

Finally, analyses of RCAN3 and RCAN3AS genomic sequences were performed to identify possible promoter regions, to examine donor and acceptor splice sequences and to compare evolutionary conservation, in particular of alternative exon 1 or 1c - exon 2 junctions in different species.

The description of its number of transcripts, of their expression patterns and of their regulatory regions can be important to clarify the functions of RCAN3 gene in different pathways and cellular processes.