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

Screening and identification of serum exosomal protein ZNF587B in liquid biopsy for ovarian cancer diagnosis

Addressing the critical challenge of early ovarian cancer (OC) detection, our study focuses on identifying novel biomarkers by analyzing preoperative peripheral blood exosomes from high-grade serous ovarian cancer (HGSC) patients and healthy controls. Utilizing high-performance liquid chromatography-mass spectrometry-based quantitative proteomics, we isolated and analyzed peripheral blood exosomes to identify differentially expressed proteins (DEPs). This comprehensive analysis, supported by gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) database assessments, revealed 28 proteins with decreased abundance and 33 with increased abundance in HGSC patients compared to controls. Notably, Zinc Finger Protein 587B (ZNF587B) exhibited a significant reduction in abundance, confirmed by decreased mRNA and protein levels in HGSC and normal ovarian tissues, consistent with omes exosomal protein expression levels. Immunohistochemical staining further confirmed reduced ZNF587B protein levels in HGSC tissues. The significant correlation between ZNF587B expression levels and tumor stage underscores its potential as a valuable biomarker for early liquid biopsy screening of OC. Our findings suggest ZNF587B plays a crucial role in early HGSC detection, highlighting the importance of further research to validate its clinical utility and improve ovarian cancer patient outcomes.

Acute Myeloid Leukemia Expresses a Specific Group of Olfactory Receptors

Acute myeloid leukemia (AML) is the most common form of acute leukemia in adults, with a 5-year overall survival rate of approximately 30%. Despite recent advances in therapeutic options, relapse remains the leading cause of death and poor survival outcomes. New drugs benefit specific small subgroups of patients with actionable therapeutic targets. Thus, finding new targets with greater applicability should be pursued. Olfactory receptors (ORs) are seven transmembrane G-protein coupled receptors preferentially expressed in sensory neurons with a critical role in recognizing odorant molecules. Recent studies have revealed ectopic expression and putative function of ORs in nonolfactory tissues and pathologies, including AML. Here, we investigated OR expression in 151 AML samples, 6400 samples of 15 other cancer types, and 11,200 samples of 51 types of healthy tissues. First, we identified 19 ORs with a distinct and major expression pattern in AML, which were experimentally validated by RT-PCR in an independent set of 13 AML samples, 13 healthy donors, and 8 leukemia cell lines. We also identified an OR signature with prognostic potential for AML patients. Finally, we found cancer-related genes coexpressed with the ORs in the AML samples. In summary, we conducted an extensive study to identify ORs that can be used as novel biomarkers for the diagnosis of AML and as potential drug targets.

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.

High Expression of DC-STAMP Gene Predicts Adverse Outcomes in AML

Acute myeloid leukemia (AML) is a genetically heterogeneous hematological malignancy with poor prognosis. We explored the RNA sequence data and clinical information of AML patients from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) database to search for the core molecule for prognosis. The DC-STAMP expression was significantly higher in AML patients, which was linked to old age, unfavorable cytogenetic risk, and death (all p < 0.05). Furthermore, it was revealed that high DC-STAMP expression was an independent unfavorable factor for overall survival (OS) by univariate analysis [hazard ratio (HR): 2.683; 95% confidence interval (CI): 1.723–4.178; p < 0.001] and multivariate analysis (HR: 1.733; 95% CI: 1.079–2.781; p = 0.023). The concordance index (C-index 0.734, 95% CI: 0.706–0.762), calibration curves, and decision curve analysis showed the certain predictive accuracy of a nomogram model based on multivariate analysis for OS. In addition, we found that the differentially expressed gene (DEG) enrichment pathways of high- and low-DC-STAMP expression group enrichment pathways were focused on channel activity and platelet alpha granule by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), while gene set enrichment analysis (GSEA) pathways were mainly involved in mTORC1 signaling and TNF-α signaling via the NF-kB pathway. Moreover, a protein–protein interaction (PPI) network demonstrated that DC-STAMP interacted with two hub genes (PPBP and PF4), which were highly regulated and associated with poor survival. Finally, high DC-STAMP expression showed a significantly positive correlation with four immune cell [NK CD56 (dim) cells, macrophages, cytotoxic cells, and CD8 (+) T cells] infiltration and high level of immune checkpoint genes (PDCD1, CD274, CTLA-4, and TIGIT). Therefore, our results suggest that high expression of DC-STAMP predicts adverse outcomes for AML patients.

The Novel Zinc Finger Protein 587B Gene, ZNF587B, Regulates Cell Proliferation and Metastasis in Ovarian Cancer Cells in vivo and in vitro

Background

The zinc finger protein 587B (ZNF587B) is a novel cisplatin-sensitive gene that was identified in our previous research by using a genome-scale CRISPR-Cas9 knockout library in ovarian cancer (OC) cell lines. ZNF587B belongs to the C2H2-type zinc finger protein (ZFP) family. Many ZFP protein could inhibit tumor development and malignancy. However, the function of ZNF587B remains unknown.

Methods

Quantitative PCR (qPCR) was utilized to compare ZNF587B mRNA expression levels in OC and normal ovarian cell lines. The small interfering RNA (siRNA) and full-length ZNF587B eukaryotic expression plasmid were constructed and transfected into OC cells later. Colony formation, 5-ethynyl-2′-deoxyuridine (EdU) assay, transwell assay, and xenograft experiment were conducted to evaluate the effect of ZNF587B on OC cells.

Results

ZNF587B was downregulated by approximately 43% and 17% in the OC cell lines SKOV3 and A2780, respectively, compared with that in the normal ovarian cell line IOSE80. Overexpression of ZNF587B reduced cell proliferation, colony formation, migration, and invasion, which could be reversed by knockdown of ZNF587B via siRNA. Xenograft experiments also confirmed that ZNF587B could suppress tumor growth. Survival data of OC patients in the SurvExpress database showed that with respect to overall survival, low-risk patients grouped by the prognostic index had a higher expression of ZNF587B and a better prognosis than high-risk group (HR = 1.77, 95% CI: 0.55–0.70, p = 0.023). Moreover, overexpression of ZNF587B promoted OC cells apoptosis when pretreated with cisplatin.

Conclusion

ZNF587B is a novel potential tumor suppressor of OC and may be a therapeutic target for OC.

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.

Transient leukemia of Down syndrome

Childhood leukemia is mostly a "developmental accident" during fetal hematopoiesis and may require multiple prenatal and postnatal "hits". The World Health Organization defines transient leukemia of Down syndrome (DS) as increased peripheral blood blasts in neonates with DS and classifies this type of leukemia as a separate entity. Although it was shown that DS predisposes children to myeloid leukemia, neither the nature of the predisposition nor the associated genetic lesions have been defined. Acute myeloid leukemia of DS is a unique disease characterized by a long pre-leukemic, myelodysplastic phase, unusual chromosomal findings and a high cure rate. In the present manuscript, we present a comprehensive review of the literature about clinical and biological findings of transient leukemia of DS (TL-DS) and link them with the genetic discoveries in the field. We address the manuscript to the pediatric generalist and especially to the next generation of pediatric hematologists.

Analysis of a nanoparticle‑enriched fraction of plasma reveals miRNA candidates for Down syndrome pathogenesis

Down syndrome (DS) is caused by the presence of part or all of a third copy of chromosome 21. DS is associated with several phenotypes, including intellectual disability, congenital heart disease, childhood leukemia and immune defects. Specific microRNAs (miRNAs/miR) have been described to be associated with DS, although none of them so far have been unequivocally linked to the pathology. The present study focuses to the best of our knowledge for the first time on the miRNAs contained in nanosized RNA carriers circulating in the blood. Fractions enriched in nanosized RNA-carriers were separated from the plasma of young participants with DS and their non-trisomic siblings and miRNAs were extracted. A microarray-based analysis on a small cohort of samples led to the identification of the three most abundant miRNAs, namely miR-16-5p, miR-99b-5p and miR-144-3p. These miRNAs were then profiled for 15 pairs of DS and non-trisomic sibling couples by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Results identified a clear differential expression trend of these miRNAs in DS with respect to their non-trisomic siblings and gene ontology analysis pointed to their potential role in a number of typical DS features, including 'nervous system development', 'neuronal cell body' and certain forms of 'leukemia'. Finally, these expression levels were associated with certain typical quantitative and qualitative clinical features of DS. These results contribute to the efforts in defining the DS-associated pathogenic mechanisms and emphasize the importance of properly stratifying the miRNA fluid vehicles in order to probe biomolecules that are otherwise hidden and/or not accessible to (standard) analysis.

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.

Sex-Specific Transcriptome Differences in Substantia Nigra Tissue: A Meta-Analysis of Parkinson's Disease Data

Parkinson’s disease (PD) is one of the most common progressive neurodegenerative diseases. Clinical and epidemiological studies indicate that sex differences, as well as genetic components and ageing, can influence the prevalence, age at onset and symptomatology of PD. This study undertook a systematic meta-analysis of substantia nigra microarray data using the Transcriptome Mapper (TRAM) software to integrate and normalize a total of 10 suitable datasets from multiple sources. Four different analyses were performed according to default parameters, to better define the segments differentially expressed between PD patients and healthy controls, when comparing men and women data sets. The results suggest a possible regulation of specific sex-biased systems in PD susceptibility. TRAM software allowed us to highlight the different activation of some genomic regions and loci involved in molecular pathways related to neurodegeneration and neuroinflammatory mechanisms.

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.

Common Subcluster Mining in Microarray Data for Molecular Biomarker Discovery

Molecular biomarkers can be potential facilitators for detection of cancer at early stage which is otherwise difficult through conventional biomarkers. Gene expression data from microarray experiments on both normal and diseased cell samples provide enormous scope to explore genetic relations of disease using computational techniques. Varied patterns of expressions of thousands of genes at different cell conditions along with inherent experimental error make the task of isolating disease related genes challenging. In this paper, we present a data mining method, common subcluster mining (CSM), to discover highly perturbed genes under diseased condition from differential expression patterns. The method builds heap through superposing near centroid clusters from gene expression data of normal samples and extracts its core part. It, thus, isolates genes exhibiting the most stable state across normal samples and constitute a reference set for each centroid. It performs the same operation on datasets from corresponding diseased samples and isolates the genes showing drastic changes in their expression patterns. The method thus finds the disease-sensitive genesets when applied to datasets of lung cancer, prostrate cancer, pancreatic cancer, breast cancer, leukemia and pulmonary arterial hypertension. In majority of the cases, few new genes are found over and above some previously reported ones. Genes with distinct deviations in diseased samples are prospective candidates for molecular biomarkers of the respective disease.

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.

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.

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.