24,189 are all the possible non-synonymous amino acid changes potentially affecting the human mitochondrial DNA. Only a tiny subset was functionally evaluated with certainty so far, while the pathogenicity of the vast majority was only assessed in-silico by software predictors. Since these tools proved to be rather incongruent, we have designed and implemented APOGEE, a machine-learning algorithm that outperforms all existing prediction methods in estimating the harmfulness of mitochondrial non-synonymous genome variations. We provide a detailed description of the underlying algorithm, of the selected and manually curated training and test sets of variants, as well as of its classification ability.

Multigenerational diabetes of adulthood is a mostly overlooked entity, simplistically lumped into the large pool of type 2 diabetes. The general aim of our research in the past few years is to unravel the genetic causes of this form of diabetes. Identifying among families with multigenerational diabetes those who carry mutations in known monogenic diabetes genes is the first step to then allow us to concentrate on remaining pedigrees in which to unravel new diabetes genes. Targeted next-generation sequencing of 27 monogenic diabetes genes was carried out in 55 family probands and identified mutations verified among their relatives by Sanger sequencing. Nine variants (in eight probands) survived our filtering/prioritization strategy. After likelihood of causality assessment by established guidelines, six variants were classified as “pathogenetic/likely pathogenetic” and two as “of uncertain significance.” Combining present results with our previous data on the six genes causing the most common forms of maturity-onset diabetes of the young allows us to infer that 23.6% of families with multigenerational diabetes of adulthood carry mutations in known monogenic diabetes genes. Our findings indicate that the genetic background of hyperglycemia is unrecognized in the vast majority of families with multigenerational diabetes of adulthood. These families now become the object of further research aimed at unraveling new diabetes genes.

In living organisms, biological clocks regulate 24 h (circadian) molecular, physiological, and behavioral rhythms to maintain homeostasis and synchrony with predictable environmental changes, in particular with those induced by Earth’s rotation on its axis. Harmonics of these circadian rhythms having periods of 8 and 12 h (ultradian) have been documented in several species. In mouse liver, harmonics of the 24-h period of gene transcription hallmarked genes oscillating with a frequency two or three times faster than circadian periodicity. Many of these harmonic transcripts enriched pathways regulating responses to environmental stress and coinciding preferentially with subjective dawn and dusk. At this time, the evolutionary history of genes with rhythmic expression is still poorly known and the role of length-of-day changes due to Earth’s rotation speed decrease over the last four billion years is totally ignored. We hypothesized that ultradian and stress anticipatory genes would be more evolutionarily conserved than circadian genes and background non-oscillating genes. To investigate this issue, we performed broad computational analyses of genes/proteins oscillating at different frequency ranges across several species and showed that ultradian genes/proteins, especially those oscillating with a 12-h periodicity, are more likely to be of ancient origin and essential in mice. In summary, our results show that genes with ultradian transcriptional patterns are more likely to be phylogenetically conserved and associated with the primeval and inevitable dawn/dusk transitions.

In several fields of research, molecular dynamics simulation techniques are exploited to evaluate the temporal motion of particles constituting water, ions, small molecules, macromolecules, or more complex systems over time. These techniques are considered difficult to setup, computationally demanding and require high specialization and scientific skills. Moreover, they need specialized computing infrastructures to run faster and make the simulation of big systems feasible. Here, we have simulated 3 systems of increasing sizes on scientific- and gaming-enabled graphic processing unit (GPU) cards with Amber, GROMACS, and NAMD and measured their performance accounting also for the market prices of the GPU cards where they were run on.

Altered functioning of the biological clock is involved in cancer onset and progression. MicroRNAs (miRNAs) interact with the clock genes modulating the function of genetically encoded molecular clockworks. Collaborative interactions may take place within the coding-noncoding RNA regulatory networks. We aimed to evaluate the cross-talk among miRNAs and clock genes in colorectal cancer (CRC). We performed an integrative analysis of miRNA-miRNA and miRNA-mRNA interactions on high-throughput molecular profiling of matched human CRC tissue and non-tumor mucosa, pinpointing core clock genes and their targeting miRNAs. Data obtained in silico were validated in CRC patients and human colon cancer cell lines. In silico we found severe alterations of clock gene–related coding-noncoding RNA regulatory networks in tumor tissues, which were later corroborated by the analysis of human CRC specimens and experiments performed in vitro. In conclusion, specific miRNAs target and regulate the transcription/translation of clock genes and clock gene-related miRNA-miRNA as well as mRNA-miRNA interactions are altered in colorectal cancer. Exploration of the interplay between specific miRNAs and genes, which are critically involved in the functioning of the biological clock, provides a better understanding of the importance of the miRNA-clock genes axis and its derangement in colorectal cancer.

Alterations in the balance of mRNA and microRNA (miRNA) expression profiles contribute to the onset and development of colorectal cancer. The regulatory functions of individual miRNA-gene pairs are widely acknowledged, but group effects are largely unexplored. We performed an integrative analysis of mRNA–miRNA and miRNA–miRNA interactions using high-throughput mRNA and miRNA expression profiles obtained from matched specimens of human colorectal cancer tissue and adjacent non-tumorous mucosa. This investigation resulted in a hypernetwork-based model, whose functional backbone was fulfilled by tight micro-societies of miRNAs. These proved to modulate several genes that are known to control a set of significantly enriched cancer-enhancer and cancer-protection biological processes, and that an array of upstream regulatory analyses demonstrated to be dependent on miR-145, a cell cycle and MAPK signaling cascade master regulator. In conclusion, we reveal miRNA-gene clusters and gene families with close functional relationships and highlight the role of miR-145 as potent upstream regulator of a complex RNA–RNA crosstalk, which mechanistically modulates several signaling pathways and regulatory circuits that when deranged are relevant to the changes occurring in colorectal carcinogenesis.

Oral squamous cell carcinoma (OSCC) is the most common oral and pharyngeal cancer, and is responsible of approximately 3% of cancers in men and 2% in women in the Western World, with increasing incidence rates in developing countries. Early detection by screening is necessary to prevent fatal disease because early, curable lesions are rarely symptomatic. The overall 5-yr survival rate is approximately 50% when surgery, radiation, or both are employed as treatment options, but lymph node involvement greatly influences this estimate, by decreasing the survival rate by about 50%. Here, we aimed at finding genetic signatures associated with lymph node metastasis in OSCC patients. We addressed this issue by whole transcriptome analysis through microarray expression profiling of a set of OSSC specimens of patients without lymph node involvement (10 patients, mean age ± SD 61.2±13.8, male 7, female 3) and with lymph node involvement (11 patients, mean age ± SD 62.1±15.1, male 8, female 3). We evidenced a gene expression signature associated to muscle contraction-related genes in specimens obtained from OSCC patients with lymph node involvement. This gene signature suggests the presence of myofibroblasts in tumor stoma of patients with lymph node involvement and emphasizes the decisive role played by myofibroblasts probably through their secretome in determining OSCC invasiveness.

Molecular dynamics (MD) simulation allows one to predict the time evolution of a system of interacting particles. It is widely used in physics, chemistry and biology to address specific questions about the structural properties and dynamical mechanisms of model systems. MD earned a great success in genome research, as it proved to be beneficial in sorting pathogenic from neutral genomic mutations. Considering their computational requirements, simulations are commonly performed on HPC computing devices, which are generally expensive and hard to administer. However, variables like the software tool used for modeling and simulation or the size of the molecule under investigation might make one hardware type or configuration more advantageous than another or even make the commodity hardware definitely suitable for MD studies. This work aims to shed lights on this aspect.