Identification of differentially expressed small RNAs and prediction of target genes in Italian Large White pigs with divergent backfat deposition.

The identification of the molecular mechanisms regulating pathways associated with the potential for fat deposition in pigs can lead to the detection of key genes and markers for the genetic improvement of fat traits. Interactions of microRNAs (miRNAs) with target RNAs regulate gene expression and modulate pathway activation in cells and tissues. In pigs, miRNA discovery is far from saturation, and the knowledge of miRNA expression in backfat tissue and particularly of the impact of miRNA variations is still fragmentary. Using RNA-seq, we characterized the small RNA (sRNA) expression profiles in Italian Large White pig backfat tissue. Comparing two groups of pigs divergent for backfat deposition, we detected 31 significant differentially expressed (DE) sRNAs: 14 up-regulated (including ssc-miR-132, ssc-miR-146b, ssc-miR-221-5p, ssc-miR-365-5p and the moRNA ssc-moR-21-5p) and 17 down-regulated (including ssc-miR-136, ssc-miR-195, ssc-miR-199a-5p and ssc-miR-335). To understand the biological impact of the observed miRNA expression variations, we used the expression correlation of DE miRNA target transcripts expressed in the same samples to define a regulatory network of 193 interactions between DE miRNAs and 40 DE target transcripts showing opposite expression profiles and being involved in specific pathways. Several miRNAs and mRNAs in the network were found to be expressed from backfat-related pig QTL. These results are informative for the complex mechanisms influencing fat traits, shed light on a new aspect of the genetic regulation of fat deposition in pigs and facilitate the prospective implementation of innovative strategies of pig genetic improvement based on genomic markers.

Large-scale individual thyroid monitoring following nuclear accidents by means of non-spectrometric devices.

In order to properly respond to an emergency caused by an accident in a nuclear power plant with a spread of radionuclides in the atmosphere, we propose a field procedure to perform a large-scale individual thyroid monitoring of internal contamination due to inhalation of 131I, by means of non-spectrometric equipment, in particular dose rate meters. Specific attention is paid to the individual monitoring of children, because of the very high radiosensitivity of the child's thyroid to the carcinogenic effects of ionising radiation. The device performance was evaluated by measuring mock iodine sources provided in the Child and Adult Thyroid Monitoring After Reactor Accident (CAThyMARA) intercomparison and, just for a scintillator dose rate meter, by means of 60 s acquisitions of healthy volunteers' thyroids. All the devices showed a remarkable accuracy in quantification of equivalent 131I activity in the thyroids of persons of all ages. The selected scintillator dose rate meter showed detection limit values resulting in a maximum committed equivalent dose to thyroid HT, assuming an acute 131I inhalation occurred five days before the measurement, equal to 10 mSv (related to five-year-old children). Considering the level of HT values associated with the calculated detection limit activities, the proposed procedure has a significant sensitivity to be used for fast internally thyroid monitoring in nuclear or radiological emergencies, allowing daily monitoring a large amount of individuals.

Transcriptional profiling of subcutaneous adipose tissue in Italian Large White pigs divergent for backfat thickness.

Fat deposition is a widely studied trait in pigs because of its implications with animal growth efficiency, technological and nutritional characteristics of meat products, but the global framework of the biological and molecular processes regulating fat deposition in pigs is still incomplete. This study describes the backfat tissue transcription profile in Italian Large White pigs and reports genes differentially expressed between fat and lean animals according to RNA-seq data. The backfat transcription profile was characterised by the expression of 23 483 genes, of which 54.1% were represented by known genes. Of 63 418 expressed transcripts, about 80% were non-previously annotated isoforms. By comparing the expression level of fat vs. lean pigs, we detected 86 robust differentially expressed transcripts, 72 more highly expressed (e.g. ACP5, BCL2A1, CCR1, CD163, CD1A, EGR2, ENPP1, GPNMB, INHBB, LYZ, MSR1, OLR1, PIK3AP1, PLIN2, SPP1, SLC11A1, STC1) and 14 lower expressed (e.g. ADSSL1, CDO1, DNAJB1, HSPA1A, HSPA1B, HSPA2, HSPB8, IGFBP5, OLFML3) in fat pigs. The main functional categories enriched in differentially expressed genes were immune system process, response to stimulus, cell activation and skeletal system development, for the overexpressed genes, and unfolded protein binding and stress response, for the underexpressed genes, which included five heat shock proteins. Adipose tissue alterations and impaired stress response are linked to inflammation and, in turn, to adipose tissue secretory activity, similar to what is observed in human obesity. Our results provide the opportunity to identify biomarkers of carcass fat traits to improve the pig production chain and to identify genetic factors that regulate the observed differential expression.

Rare Risk Variants Identification by Identity-by-Descent Mapping and Whole-Exome Sequencing Implicates Neuronal Development Pathways in Schizophrenia and Bipolar Disorder.

Schizophrenia (SCZ) and bipolar disorder (BPD) are highly heritable disorders with an estimated co-heritability of 68%. Hundreds of common alleles have been implicated, but recently a role for rare, high-penetrant variants has been also suggested in both disorders. This study investigated a familial cohort of SCZ and BPD patients from a closed population sample, where the high recurrence of the disorders and the homogenous genetic background indicate a possible enrichment in rare risk alleles. A total of 230 subjects (161 cases, 22 unaffected relatives, and 47 controls) were genetically investigated through an innovative strategy that integrates identity-by-descent (IBD) mapping and whole-exome sequencing (WES). IBD analysis allowed to track high-risk haplotypes (IBDrisk) shared exclusively by multiple patients from different families and possibly carrying the most penetrant alleles. A total of 444 non-synonymous sequence variants, of which 137 disruptive, were identified in IBDrisk haplotypes by WES. Interestingly, gene sets previously implicated in SCZ (i.e., post-synaptic density (PSD) proteins, voltage-gated calcium channels (VGCCs), and fragile X mental retardation protein (FMRP) targets) were found significantly enriched in genes carrying IBDrisk variants. Further, IBDrisk variants were preferentially affecting genes involved in the extracellular matrix (ECM) biology and axon guidance processes which appeared to be functionally connected in the pathway-derived meta-network analysis. Results thus confirm rare risk variants as key factors in SCZ and BPD pathogenesis and highlight a role for the development of neuronal connectivity in the etiology of both disorders.

CircRNAs in hematopoiesis and hematological malignancies.

Cell states in hematopoiesis are controlled by master regulators and by complex circuits of a growing family of RNA species impacting cell phenotype maintenance and plasticity. Circular RNAs (circRNAs) are rapidly gaining the status of particularly stable transcriptome members with distinctive qualities. RNA-seq identified thousands of circRNAs with developmental stage- and tissue-specific expression corroborating earlier suggestions that circular isoforms are a natural feature of the cell expression program. CircRNAs are abundantly expressed also in the hematopoietic compartment. There are a number of studies on circRNAs in blood cells, a specific overview is however lacking. In this review we first present current insight in circRNA biogenesis discussing the relevance for hematopoiesis of the highly interleaved processes of splicing and circRNA biogenesis. Regarding molecular functions circRNAs modulate host gene expression, but also compete for binding of microRNAs, RNA-binding proteins or translation initiation and participate in regulatory circuits. We examine circRNA expression in the hematopoietic compartment and in hematologic malignancies and review the recent breakthrough study that identified pathogenic circRNAs derived from leukemia fusion genes. CircRNA high and regulated expression in blood cell types indicate that further studies are warranted to inform the position of these regulators in normal and malignant hematopoiesis.

A data-driven network model of primary myelofibrosis: transcriptional and post-transcriptional alterations in CD34+ cells.

microRNAs (miRNAs) are relevant in the pathogenesis of primary myelofibrosis (PMF) but our understanding is limited to specific target genes and the overall systemic scenario islacking. By both knowledge-based and ab initio approaches for comparative analysis of CD34+ cells of PMF patients and healthy controls, we identified the deregulated pathways involving miRNAs and genes and new transcriptional and post-transcriptional regulatory circuits in PMF cells. These converge in a unique and integrated cellular process, in which the role of specific miRNAs is to wire, co-regulate and allow a fine crosstalk between the involved processes. The PMF pathway includes Akt signaling, linked to Rho GTPases, CDC42, PLD2, PTEN crosstalk with the hypoxia response and Calcium-linked cellular processes connected to cyclic AMP signaling. Nested on the depicted transcriptional scenario, predicted circuits are reported, opening new hypotheses. Links between miRNAs (miR-106a-5p, miR-20b-5p, miR-20a-5p, miR-17-5p, miR-19b-3p and let-7d-5p) and key transcription factors (MYCN, ATF, CEBPA, REL, IRF and FOXJ2) and their common target genes tantalizingly suggest new path to approach the disease. The study provides a global overview of transcriptional and post-transcriptional deregulations in PMF, and, unifying consolidated and predicted data, could be helpful to identify new combinatorial therapeutic strategy. Interactive PMF network model: http://compgen.bio.unipd.it/pmf-net/.

DHPLC analysis of the MECP2 gene in Italian Rett patients.

Rett Syndrome (RTT) is an X-linked dominant neurodevelopmental disorder, which almost exclusively affects girls, with an estimated prevalence of one in 10,000-15,000 female births. Mutations in the methyl CpG binding protein 2 gene (MECP2) have been identified in roughly 75% of classical Rett girls. The vast majority of Rett cases (99%) are sporadic in origin, and are due to de novo mutations. We collected DNA samples from 50 Italian classical Rett girls, and screened the MECP2 coding region for mutations by denaturing high-performance liquid chromatography (DHPLC) and subsequent direct sequencing. DHPLC is a recently developed method for mutation screening which identifies heteroduplexes formed in DNA samples containing mismatches between wild type and mutant DNA strands, combining high sensitivity, reduced cost per run, and high throughput. In our series, 19 different de novo MECP2 mutations, eight of which were previously unreported, were found in 35 out of 50 Rett girls (70%). Seven recurrent mutations were characterized in a total of 22 unrelated cases. Initial DHPLC screening allowed the identification of 17 out of 19 different mutations (90%); after optimal conditions were established, this figure increased to 100%, with all recurrent MECP2 mutations generating a characteristic chromatographic profile. Detailed clinical data were available for 27 out of 35 mutation carrying Rett girls. Milder disease was detectable in patients carrying nonsense mutation as compared to patients carrying missense mutations, although this difference was not statistically significant (P = 0.077).

Disease genes and intracellular protein networks.

By a computational approach we reconstructed genomic transcriptional profiles of 19 different adult human tissues, based on information on activity of 27,924 genes obtained from unbiased UniGene cDNA libraries. In each considered tissue, a small number of genes resulted highly expressed or "tissue specific." Distribution of gene expression levels in a tissue appears to follow a power law, thus suggesting a correspondence between transcriptional profile and "scale-free" topology of protein networks. The expression of 737 genes involved in Mendelian diseases was analyzed, compared with a large reference set of known human genes. Disease genes resulted significantly more expressed than expected. The possible correspondence of their products to important nodes of intracellular protein network is suggested. Auto-organization of the protein network, its stability in time in the differentiated state, and relationships with the degree of genetic variability at genome level are discussed.

A novel resource for the study of genes expressed in the adult human retina.

PURPOSE: To reconstruct the transcriptional profile of the human adult retina and the genomic map of the genes expressed in this tissue. METHODS: Original software was used for the retrieval and analysis of records from UniGene (http://www.ncbi.nlm.nih. gov/UniGene/) pertaining to selected cDNA libraries from adult human retina. RESULTS: The 4974 genes reported so far to be expressed in retina were included in a catalog available on the Internet. For each entry, an estimation of the level of expression of the corresponding gene in the retina was provided. A high-resolution genomic map of the human retina was built up by inclusion of 3152 genes showing a precise and unique map assignment. The correspondence was established between 53 gene-orphan retinal diseases and clusters of genes expressed in the retina. CONCLUSIONS: The in silico reconstruction of the transcriptional profile of the adult human retina provides preliminary information on the pattern of genomic expression in this tissue. The chromosomal location of many retinal genes, combined with their expression data, should speed up the identification of genes involved in retinal diseases.

Environmental radioactivity analyses in Italy following the Fukushima Dai-ichi nuclear accident.

Following the Fukushima power plants accident on the 11th March 2011, the radioactivity monitoring programme at the Italian ENEA research centres was activated in order to detect the possible new input of radionuclides through atmospheric transport and precipitation. Measurements of (131)I and (134,137)Cs were carried out on atmospheric particulate, atmospheric deposition, seawater and mussels and sheep milk. In the daily samples of air particulate, (131)I was detectable between March 28 and April 12, with extremely low concentrations (<1 mBq m(-3); the detection limit for (131)I was ~0.2 mBq m(-3)) while Cs isotopes were always below the detection limit (<0.2 mBq m(-3)). The two main episodes of (131)I atmospheric deposition were registered in La Spezia research centre, around March 28 and April 15, reaching values of 17.8 ± 1.1 and 8.0 ± 2.5 Bq m(-2) respectively; maximum values of (134)Cs and (137)Cs were 0.11 ± 0.03 and 0.17 ± 0.02 Bq m(-2), respectively, detected in Brasimone research centre in April (reference date April 15). Mussels and seawater were collected in the Gulf of La Spezia: only mussels after the main (131)I deposition, on March 28, contained a measurable, although very small, amount of (131)I (0.18 ± 0.05 Bq kg(-1), detection limit (131)I = 0.03 Bq kg(-1) wet weight - soft parts). The (131)I was also detected in sheep milk in Rome (Casaccia research centre) until May 5, showing a maximum concentration of 4.9 ± 0.4 Bq L(-1). As for other European Countries for which data are available, activity levels remain of no concern for public health.

Age estimation in subadult Egyptian remains.

In anthropological analyses of past populations, it is very important to be able to accurately reconstruct the palaeodemographic profile in order to interpret infant mortality as an indicator of the environmental, social and cultural conditions. There are various methods to evaluate the age of immature individuals but some of these methods are strongly influenced by the different rates of skeletal development observed in populations from various geographical areas and/or from various time periods, as well as between the sexes. Clearly, there is a need for adopting a method of estimation of age at death, which will be the one most suitable for analysing the particular skeletal sample. In this study we investigated subadults from the Egyptian osteological collection housed in the Museum of Anthropology and Ethnography of the University of Turin. For each individual, the age at death was estimated based on the degree of eruption and mineralisation of the teeth. Then the estimated age at death was correlated with the measurements of the long bones and ilium. We showed that greater regularity and constancy of rates of skeletal growth could be assessed with measurements, alternative to using maximum length of diaphysis. Moreover, using alternative characters, it was possible to markedly increase the number of individuals whose age at death could be estimated. Our study also showed the need to use a reference sample consistent with the sample being analysed and, which was derived from similar biological-environmental context. Therefore, our proposed method can be used for the estimation of age at death in pre/protohistorical populations from the Mediterranean region.

A computational reconstruction of the adult human heart transcriptional profile.

The reconstruction of the transcriptional profile of the adult human heart was attempted, by applying a bioinformatic and computational approach to UniGene data. A catalogue of 2077 expressed genes was produced. Over 1000 entries of the catalogue corresponded to putative novel genes. Highly expressed genes accounted for about 20% of the total. Almost all genes expressed in adult heart resulted to be active in at least one additional tissue and about 90% were found in over five additional tissues. A genomic map of 1364 genes expressed in heart, which also indicated chromosomal location, was produced, which could be conveniently used for the discovery of the determinants of gene-orphan heart diseases and for the detection of clusters of highly expressed genes. The catalogue and the genomic map of genes expressed in adult human heart are available on Internet at the sites: http://telethon.bio.unipd.it/GETProfiles/heart and http://telethon.bio.unipd.it/GETMaps/heart.

Detecting differentially expressed genes in multiple tag sampling experiments: comparative evaluation of statistical tests.

The comparison of several statistical methods currently used for detection of differentially expressed genes was attempted both by a simulation approach and by the analysis of data sets of human expressed sequence tags, obtained from UniGene. In the simulated mixed case, mimicking a situation close to reality, the general chi(2) test was unexpectedly the most efficient in multiple tag sampling experiments, especially when dealing with variations affecting weakly expressed genes. On the other hand, Audic and Claverie's method proved the most efficient for detecting differences in gene expression when dealing with pairwise comparisons. By applying the above methods on UniGene-based data sets concerning two human kidney tumours compared with normal kidney tissue, three novel genes overexpressed in these tumours were identified. Software and additional information on statistical methodologies, simulation approach and data are available at http://telethon.bio.unipd.it/bioinfo/IDEG6/.

The human adult skeletal muscle transcriptional profile reconstructed by a novel computational approach.

By applying a novel software tool, information on 4080 UniGene clusters was retrieved from three adult human skeletal muscle cDNA libraries, which were selected for being neither normalized nor subtracted. Reconstruction of a transcriptional profile of the corresponding tissue was attempted by a computational approach, classifying each transcript according to its level of expression. About 25% of the transcripts accounted for about 80% of the detected transcriptional activity, whereas most genes showed a low level of expression. This in silico transcriptional profile was then compared with data obtained by a SAGE study. A fairly good agreement between the two methods was observed. About 400 genes, highly expressed in skeletal muscle or putatively skeletal muscle-specific, may represent the minimal set of genes needed to determine the tissue specificity. These genes could be used as a convenient reference to monitor major changes in the transcriptional profile of adult human skeletal muscle in response to different physiological or pathological conditions, thus providing a framework for designing DNA microarrays and initiating biological studies.

Characterization of thin layers of glucose oxidase.

Glucose Oxidase (GOD) has been covalently bound to functionalized glass cover slips. The surface density of immobilized GOD molecules was measured by a method based on the amperometric determination of Flavin Adenine Dinucleotide (FAD). Atomic Force Microscopy (AFM) images, obtained in aqueous solution for the covalently bound enzyme, show a monomolecular layer of the enzyme on a functionalized glass surface. The catalytic constants were measured for the immobilized GOD and compared with those of the free enzyme.

Differential expression of genes coding for ribosomal proteins in different human tissues.

MOTIVATION: To perform a computational and statistical study on a large set of gene expression data pertaining six adult human tissues (brain, liver, skeletal muscle, ovary, retina and uterus) for analyzing the expression of ribosomal protein genes. RESULTS: Unexpectedly, in each of the considered tissues large variations in the expression of ribosomal protein genes were observed. Moreover, when comparing the expression levels of 89 ribosomal protein genes in six different tissues, 13 genes appeared differentially expressed among tissues. AVAILABILITY: The expression data of the ribosomal protein genes together with supplementary material (complete transcriptional profiles of the considered human tissues) are freely available at the site GETProfiles (http://telethon.bio.unipd.it/GETProfiles/). CONTACT: danieli@bio.unipd.it

Chromosomal localization of four MAPK signaling cascade genes: MEK1, MEK3, MEK4 and MEKK5.

The mitogen-activated protein kinase (MAPK) signaling cascade is one of the most important mechanisms for the cytoplasmic transduction of extracellular signals. We report the chromosomal localization of the human MEK1, MEK3, MEK4 and MEKK5 genes, involved in the MAPK cascade. Using radiation hybrid mapping, MEK1 was assigned to chromosome 15q22.1 --> q22.33, MEK3 to chromosome 17q11.2, MEK4 to chromosome 17p12, and MEKK5 to chromosome 6q22.33.