A case of Ververi-Brady syndrome due to QRICH1 loss of function and the literature review.

Ververi-Brady syndrome (VBS), first reported in 2018, is characterized by intellectual disability, speech delay, and mild dysmorphic facial features. VBS has been linked to de novo loss-of-function variants in the glutamine-rich protein 1 (QRICH1) on chromosome 3p21 and was reported until lately in only five individuals. Four additional cases have just been described substantiating the notion that children with VBS are mildly dysmorphic, mildly to moderately intellectually disabled, have linear growth shortage, are picky eaters, and have notable attention and social behavioral deficits. We describe a new patient and review the clinical and genetic information, on all previously reported VBS cases. The child here reported is noted for maladaptive behavior, sensory hypersensitivity, and slow linear growth. He is mainly hyperactive, distractible, impulsive, and inattentive. His speech, initially delayed, is fair and his verbal comprehension age adequate.

Using Synthetic Mouse Spike-In Transcripts to Evaluate RNA-Seq Analysis Tools.

One of the key applications of next-generation sequencing (NGS) technologies is RNA-Seq for transcriptome genome-wide analysis. Although multiple studies have evaluated and benchmarked RNA-Seq tools dedicated to gene level analysis, few studies have assessed their effectiveness on the transcript-isoform level. Alternative splicing is a naturally occurring phenomenon in eukaryotes, significantly increasing the biodiversity of proteins that can be encoded by the genome. The aim of this study was to assess and compare the ability of the bioinformatics approaches and tools to assemble, quantify and detect differentially expressed transcripts using RNA-Seq data, in a controlled experiment. To this end, in vitro synthesized mouse spike-in control transcripts were added to the total RNA of differentiating mouse embryonic bodies, and their expression patterns were measured. This novel approach was used to assess the accuracy of the tools, as established by comparing the observed results versus the results expected of the mouse controlled spiked-in transcripts. We found that detection of differential expression at the gene level is adequate, yet on the transcript-isoform level, all tools tested lacked accuracy and precision.

VAV1 and BAFF, via NFκB pathway, are genetic risk factors for myasthenia gravis.

OBJECTIVE: To identify novel genetic loci that predispose to early-onset myasthenia gravis (EOMG) applying a two-stage association study, exploration, and replication strategy. METHODS: Thirty-four loci and one confirmation loci, human leukocyte antigen (HLA)-DRA, were selected as candidate genes by team members of groups involved in different research aspects of MG. In the exploration step, these candidate genes were genotyped in 384 EOMG and 384 matched controls and significant difference in allele frequency were found in eight genes. In the replication step, eight candidate genes and one confirmation loci were genotyped in 1177 EOMG patients and 814 controls, from nine European centres. RESULTS: ALLELE FREQUENCY DIFFERENCES WERE FOUND IN FOUR NOVEL LOCI: CD86, AKAP12, VAV1, B-cell activating factor (BAFF), and tumor necrosis factor-alpha (TNF-α), and these differences were consistent in all nine cohorts. Haplotype trend test supported the differences in allele frequencies between cases and controls. In addition, allele frequency difference in female versus male patients at HLA-DRA and TNF-α loci were observed. INTERPRETATION: The genetic associations to EOMG outside the HLA complex are novel and of interest as VAV1 is a key signal transducer essential for T- and B-cell activation, and BAFF is a cytokine that plays important roles in the proliferation and differentiation of B-cells. Moreover, we noted striking epistasis between the predisposing VAV1 and BAFF haplotypes; they conferred a greater risk in combination than alone. These, and CD86, share the same signaling pathway, namely nuclear factor-kappaB (NFκB), thus implicating dysregulation of proinflammatory signaling in predisposition to EOMG.

Differences in microRNA detection levels are technology and sequence dependent.

Identification and quantification of small RNAs are challenging because of their short length, high sequence similarities within microRNA (miRNA) families, and the existence of miRNA isoforms and O-methyl 3' modifications. In this study, the detection performance of three high-throughput commercial platforms, Agilent and Affymetrix microarrays and Illumina next-generation sequencing, was systematically and comprehensively compared. The ability to detect miRNAs was shown to depend strongly on the platform and on miRNA modifications and sequence. Using synthetic transcripts, including mature, precursor, and O-methyl-modified miRNAs spiked into human RNA, a large intensity variation in all spiked-in miRNAs and a reduced capacity in detecting O-methyl-modified miRNAs were observed between the tested platforms. In addition, endogenous human miRNA expression levels were assessed across the platforms. Detected miRNA expression levels were not consistent between platforms. Although biases in miRNA detection were previously described, here the end-point result, i.e., detection intensity, of these biases was investigated on multiple platforms in a controlled fashion. A detailed exploration of a large number of attributes, including base composition, sequence structure, and isoform miRNA attributes, suggests their impact on miRNA expression detection level. This study provides a basis for understanding the attributes that should be considered to adjust platform-dependent detection biases.

Epigenetic polymorphism and the stochastic formation of differentially methylated regions in normal and cancerous tissues.

DNA methylation has been comprehensively profiled in normal and cancer cells, but the dynamics that form, maintain and reprogram differentially methylated regions remain enigmatic. Here, we show that methylation patterns within populations of cells from individual somatic tissues are heterogeneous and polymorphic. Using in vitro evolution of immortalized fibroblasts for over 300 generations, we track the dynamics of polymorphic methylation at regions developing significant differential methylation on average. The data indicate that changes in population-averaged methylation occur through a stochastic process that generates a stream of local and uncorrelated methylation aberrations. Despite the stochastic nature of the process, nearly deterministic epigenetic remodeling emerges on average at loci that lose or gain resistance to methylation accumulation. Changes in the susceptibility to methylation accumulation are correlated with changes in histone modification and CTCF occupancy. Characterizing epigenomic polymorphism within cell populations is therefore critical to understanding methylation dynamics in normal and cancer cells.

Expression profiling and biochemical analysis suggest stress response as a potential mechanism inhibiting proliferation of polyamine-depleted cells.

Polyamines are small organic polycations that are absolutely required for cell growth and proliferation; yet the basis for this requirement is mostly unknown. Here, we combined a genome-wide expression profiling with biochemical analysis to reveal the molecular basis for inhibited proliferation of polyamine-depleted cells. Transcriptional responses accompanying growth arrest establishment in polyamine-depleted cells or growth resumption following polyamine replenishment were monitored and compared. Changes in the expression of genes related to various fundamental cellular processes were established. Analysis of mirror-symmetric expression patterns around the G(1)-arrest point identified a set of genes representing a stress-response signature. Indeed, complementary biochemical analysis demonstrated activation of the PKR-like endoplasmic reticulum kinase arm of the unfolded protein response and of the stress-induced p38 MAPK. These changes were accompanied by induction of key growth-inhibitory factors such as p21 and Gadd45a and reduced expression of various cyclins, most profoundly cyclin D1, setting the basis for the halted proliferation. However, although the induced stress response could arrest growth, polyamine depletion also inhibited proliferation of PKR-like endoplasmic reticulum kinase and p38α-deficient cells and of cells harboring a nonphosphorylatable mutant eIF2α (S51A), suggesting that additional yet unidentified mechanisms might inhibit proliferation of polyamine-depleted cells. Despite lengthy persistence of the stress and activation of apoptotic signaling, polyamine-depleted cells remained viable, apparently due to induced expression of protective genes and development of autophagy.

The H3K27 demethylase Utx regulates somatic and germ cell epigenetic reprogramming.

Induced pluripotent stem cells (iPSCs) can be derived from somatic cells by ectopic expression of different transcription factors, classically Oct4 (also known as Pou5f1), Sox2, Klf4 and Myc (abbreviated as OSKM). This process is accompanied by genome-wide epigenetic changes, but how these chromatin modifications are biochemically determined requires further investigation. Here we show in mice and humans that the histone H3 methylated Lys 27 (H3K27) demethylase Utx (also known as Kdm6a) regulates the efficient induction, rather than maintenance, of pluripotency. Murine embryonic stem cells lacking Utx can execute lineage commitment and contribute to adult chimaeric animals; however, somatic cells lacking Utx fail to robustly reprogram back to the ground state of pluripotency. Utx directly partners with OSK reprogramming factors and uses its histone demethylase catalytic activity to facilitate iPSC formation. Genomic analysis indicates that Utx depletion results in aberrant dynamics of H3K27me3 repressive chromatin demethylation in somatic cells undergoing reprogramming. The latter directly hampers the derepression of potent pluripotency promoting gene modules (including Sall1, Sall4 and Utf1), which can cooperatively substitute for exogenous OSK supplementation in iPSC formation. Remarkably, Utx safeguards the timely execution of H3K27me3 demethylation observed in embryonic day 10.5-11 primordial germ cells (PGCs), and Utx-deficient PGCs show cell-autonomous aberrant epigenetic reprogramming dynamics during their embryonic maturation in vivo. Subsequently, this disrupts PGC development by embryonic day 12.5, and leads to diminished germline transmission in mouse chimaeras generated from Utx-knockout pluripotent cells. Thus, we identify Utx as a novel mediator with distinct functions during the re-establishment of pluripotency and germ cell development. Furthermore, our findings highlight the principle that molecular regulators mediating loss of repressive chromatin during in vivo germ cell reprogramming can be co-opted during in vitro reprogramming towards ground state pluripotency.

Signaling to transcription networks in the neuronal retrograde injury response.

Retrograde signaling from axon to soma activates intrinsic regeneration mechanisms in lesioned peripheral sensory neurons; however, the links between axonal injury signaling and the cell body response are not well understood. Here, we used phosphoproteomics and microarrays to implicate approximately 900 phosphoproteins in retrograde injury signaling in rat sciatic nerve axons in vivo and approximately 4500 transcripts in the in vivo response to injury in the dorsal root ganglia. Computational analyses of these data sets identified approximately 400 redundant axonal signaling networks connected to 39 transcription factors implicated in the sensory neuron response to axonal injury. Experimental perturbation of individual overrepresented signaling hub proteins, including Abl, AKT, p38, and protein kinase C, affected neurite outgrowth in sensory neurons. Paradoxically, however, combined perturbation of Abl together with other hub proteins had a reduced effect relative to perturbation of individual proteins. Our data indicate that nerve injury responses are controlled by multiple regulatory components, and suggest that network redundancies provide robustness to the injury response.

A critical role for ceramide synthase 2 in liver homeostasis: II. insights into molecular changes leading to hepatopathy.

We have generated a mouse that cannot synthesize very long acyl chain (C22-C24) ceramides (Pewzner-Jung, Y., Park, H., Laviad, E. L., Silva, L. C., Lahiri, S., Stiban, J., Erez-Roman, R., Brugger, B., Sachsenheimer, T., Wieland, F. T., Prieto, M., Merrill, A. H., and Futerman, A. H. (2010) J. Biol. Chem. 285, 10902-10910) due to ablation of ceramide synthase 2 (CerS2). As a result, significant changes were observed in the sphingolipid profile of livers from these mice, including elevated C16-ceramide and sphinganine levels. We now examine the functional consequences of these changes. CerS2 null mice develop severe nonzonal hepatopathy from about 30 days of age, the age at which CerS2 expression peaks in wild type mice, and display increased rates of hepatocyte apoptosis and proliferation. In older mice there is extensive and pronounced hepatocellular anisocytosis with widespread formation of nodules of regenerative hepatocellular hyperplasia. Progressive hepatomegaly and noninvasive hepatocellular carcinoma are also seen from approximately 10 months of age. Even though CerS2 is found at equally high mRNA levels in kidney and liver, there are no changes in renal function and no pathological changes in the kidney. High throughput analysis of RNA expression in liver revealed up-regulation of genes associated with cell cycle regulation, protein transport, cell-cell interactions and apoptosis, and down-regulation of genes associated with intermediary metabolism, such as lipid and steroid metabolism, adipocyte signaling, and amino acid metabolism. In addition, levels of the cell cycle regulator, the cyclin dependent-kinase inhibitor p21(WAF1/CIP1), were highly elevated, which occurs by at least two mechanisms, one of which may involve p53. We propose a functional rationale for the synthesis of sphingolipids with very long acyl chains in liver homeostasis and in cell physiology.

Expoldb: expression linked polymorphism database with inbuilt tools for analysis of expression and simple repeats.

BACKGROUND: Quantitative variation in gene expression has been proposed to underlie phenotypic variation among human individuals. A facilitating step towards understanding the basis for gene expression variability is associating genome wide transcription patterns with potential cis modifiers of gene expression. DESCRIPTION: EXPOLDB, a novel Database, is a new effort addressing this need by providing information on gene expression levels variability across individuals, as well as the presence and features of potentially polymorphic (TG/CA)n repeats. EXPOLDB thus enables associating transcription levels with the presence and length of (TG/CA)n repeats. One of the unique features of this database is the display of expression data for 5 pairs of monozygotic twins, which allows identification of genes whose variability in expression, are influenced by non-genetic factors including environment. In addition to queries by gene name, EXPOLDB allows for queries by a pathway name. Users can also upload their list of HGNC (HUGO (The Human Genome Organisation) Gene Nomenclature Committee) symbols for interrogating expression patterns. The online application 'SimRep' can be used to find simple repeats in a given nucleotide sequence. To help illustrate primary applications, case examples of Housekeeping genes and the RUNX gene family, as well as one example of glycolytic pathway genes are provided. CONCLUSION: The uniqueness of EXPOLDB is in facilitating the association of genome wide transcription variations with the presence and type of polymorphic repeats while offering the feature for identifying genes whose expression variability are influenced by non genetic factors including environment. In addition, the database allows comprehensive querying including functional information on biochemical pathways of the human genes. EXPOLDB can be accessed at http://expoldb.igib.res.in/expol.

Assessing natural variations in gene expression in humans by comparing with monozygotic twins using microarrays.

Quantitative variation in gene expression in humans is the outcome of various factors, including differences in genetic background, gender, age, and environment. However, the extent of the influence of these factors on gene expression is not clear. We attempted to address this issue by carrying out gene expression profiling in blood leukocytes with 13 individuals (including 5 pairs of monozygotic twins) on 10,000 genes using HG-U95Av2 oligonucleotide microarrays. The proportion of differentially expressed genes between monozygotic twins was low (up to 1.76%). Most of the variations belonged to the least variable category. These genes, exhibiting "random variations," did not show clear preference to any functional class, although "signaling and communication" and "immune and related functions" generally topped the list. The extent of variation in gene expression increased in comparisons between unrelated individuals (up to 14.13%). Most of the genes (89%) exhibiting random variations in twins also varied in expression in unrelated individuals. As with twins, signaling and communication topped the list, and substantial variations were observed in all three categories: least variable, moderately variable, and most variable. An important outcome of this study was that the housekeeping genes were nearly insensitive to random variations but appeared to be more susceptible to genetic differences. However, the highly expressed housekeeping genes exhibited low variation and appeared to be insensitive to all known factors. Gene expression profiling in monozygotic twins can provide useful data for the assessment of natural variation in gene expression in humans.

Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification.

MOTIVATION: Genes are often characterized dichotomously as either housekeeping or single-tissue specific. We conjectured that crucial functional information resides in genes with midrange profiles of expression. RESULTS: To obtain such novel information genome-wide, we have determined the mRNA expression levels for one of the largest hitherto analyzed set of 62 839 probesets in 12 representative normal human tissues. Indeed, when using a newly defined graded tissue specificity index tau, valued between 0 for housekeeping genes and 1 for tissue-specific genes, genes with midrange profiles having 0.15< tau<0.85 were found to constitute >50% of all expression patterns. We developed a binary classification, indicating for every gene the I(B) tissues in which it is overly expressed, and the 12-I(B) tissues in which it shows low expression. The 85 dominant midrange patterns with I(B)=2-11 were found to be bimodally distributed, and to contribute most significantly to the definition of tissue specification dendrograms. Our analyses provide a novel route to infer expression profiles for presumed ancestral nodes in the tissue dendrogram. Such definition has uncovered an unsuspected correlation, whereby de novo enhancement and diminution of gene expression go hand in hand. These findings highlight the importance of gene suppression events, with implications to the course of tissue specification in ontogeny and phylogeny. AVAILABILITY: All data and analyses are publically available at the GeneNote website, http://genecards.weizmann.ac.il/genenote/ and, GEO accession GSE803. CONTACT: doron.lancet@weizmann.ac.il SUPPLEMENTARY INFORMATION: Four tables available at the above site.

Prediction of high risk Ewing's sarcoma by gene expression profiling.

Ewing's sarcoma (ES) is the second most common primary malignant bone tumor in children and adolescents. Currently accepted clinical prognostic factors fail to classify ES patients' risk to relapse at diagnosis. We aimed to find a new strategy to distinguish between poor and good prognosis ES patients already at diagnosis. We analysed the gene expression profiles of 14 primary tumor specimens and six metastases from ES patients, using oligonucleotide microarray analysis. The over-expression of two genes was validated by quantitative PCR using the LightCycler system. We identified two distinct gene expression signatures distinguishing high-risk ES patients that are likely to progress from low-risk ES patients with a favorable prognosis of long-term progression-free survival. The microarray-based classification was superior to currently used prognostic parameters. Over-expressed genes in the poor prognosis patients included genes regulating the cell cycle and genes associated with invasion and metastasis, while among the downregulated genes were tumor suppressor genes and inducers of apoptosis. Our results indicate the existence of a specific gene expression signature of outcome in ES already at diagnosis, and provide a strategy to select patients who would benefit from risk-adapted improved therapy.

GeneNote: whole genome expression profiles in normal human tissues.

A novel data set, GeneNote (Gene Normal Tissue Expression), was produced to portray complete gene expression profiles in healthy human tissues using the Affymetrix GeneChip HG-U95 set, which includes 62 839 probe-sets. The hybridization intensities of two replicates were processed and analyzed to yield the complete transcriptome for twelve human tissues. Abundant novel information on tissue specificity provides a baseline for past and future expression studies related to diseases. The data is posted in GeneNote (http://genecards.weizmann.ac.il/genenote/), a widely used compendium of human genes (http://bioinfo.weizmann.ac.il/genecards).