The miR-200c/141-ZEB2-TGFß axis is aberrant in human T-cell prolymphocytic leukemia.

T-cell prolymphocytic leukemia (T-PLL) is mostly characterized by aberrant expansion of small- to medium-sized prolymphocytes with a mature post-thymic phenotype, high aggressiveness of the disease and poor prognosis. However, T-PLL is more heterogeneous with a wide range of clinical, morphological, and molecular features, which occasionally impedes the diagnosis. We hypothesized that T-PLL consists of phenotypic and/or genotypic subgroups that may explain the heterogeneity of the disease. Multi-dimensional immuno-phenotyping and gene expression profiling did not reveal clear T-PLL subgroups, and no clear T-cell receptor a or ß CDR3 skewing was observed between different T-PLL cases. We revealed that the expression of microRNA (miRNA) is aberrant and often heterogeneous in T-PLL. We identified 35 miRNA that were aberrantly expressed in T-PLL with miR-200c/141 as the most differentially expressed cluster. High miR- 200c/141 and miR-181a/181b expression was significantly correlated with increased white blood cell counts and poor survival. Furthermore, we found that overexpression of miR-200c/141 correlated with downregulation of their targets ZEB2 and TGFßR3 and aberrant TGFß1- induced phosphorylated SMAD2 (p-SMAD2) and p-SMAD3, indicating that the TGFß pathway is affected in T-PLL. Our results thus highlight the potential role for aberrantly expressed oncogenic miRNA in T-PLL and pave the way for new therapeutic targets in this disease.

Evaluation of FRET X for single-molecule protein fingerprinting.

Single-molecule protein identification is an unrealized concept with potentially ground-breaking applications in biological research. We propose a method called FRET X (Förster Resonance Energy Transfer via DNA eXchange) fingerprinting, in which the FRET efficiency is read out between exchangeable dyes on protein-bound DNA docking strands and accumulated FRET efficiencies constitute the fingerprint for a protein. To evaluate the feasibility of this approach, we simulated fingerprints for hundreds of proteins using a coarse-grained lattice model and experimentally demonstrated FRET X fingerprinting on model peptides. Measured fingerprints are in agreement with our simulations, corroborating the validity of our modeling approach. In a simulated complex mixture of >300 human proteins of which only cysteines, lysines, and arginines were labeled, a support vector machine was able to identify constituents with 95% accuracy. We anticipate that our FRET X fingerprinting approach will form the basis of an analysis tool for targeted proteomics.

Hap10: reconstructing accurate and long polyploid haplotypes using linked reads.

BACKGROUND: Haplotype information is essential for many genetic and genomic analyses, including genotype-phenotype associations in human, animals and plants. Haplotype assembly is a method for reconstructing haplotypes from DNA sequencing reads. By the advent of new sequencing technologies, new algorithms are needed to ensure long and accurate haplotypes. While a few linked-read haplotype assembly algorithms are available for diploid genomes, to the best of our knowledge, no algorithms have yet been proposed for polyploids specifically exploiting linked reads. RESULTS: The first haplotyping algorithm designed for linked reads generated from a polyploid genome is presented, built on a typical short-read haplotyping method, SDhaP. Using the input aligned reads and called variants, the haplotype-relevant information is extracted. Next, reads with the same barcodes are combined to produce molecule-specific fragments. Then, these fragments are clustered into strongly connected components which are then used as input of a haplotype assembly core in order to estimate accurate and long haplotypes. CONCLUSIONS: Hap10 is a novel algorithm for haplotype assembly of polyploid genomes using linked reads. The performance of the algorithms is evaluated in a number of simulation scenarios and its applicability is demonstrated on a real dataset of sweet potato.

Family-Based Haplotype Estimation and Allele Dosage Correction for Polyploids Using Short Sequence Reads.

DNA sequence reads contain information about the genomic variants located on a single chromosome. By extracting and extending this information using the overlaps between the reads, the haplotypes of an individual can be obtained. Using parent-offspring relationships in a population can considerably improve the quality of the haplotypes obtained from short reads, as pedigree information can be used to correct for spurious overlaps (due to sequencing errors) and insufficient overlaps (due to short read lengths, low genomic variation and shallow coverage). We developed a novel method, PopPoly, to estimate polyploid haplotypes in an F1-population from short sequence data by taking into consideration the transmission of the haplotypes from the parents to the offspring. In addition, this information is employed to improve genotype dosage estimation and to call missing genotypes in the population. Through simulations, we compare PopPoly to other haplotyping methods and show its better performance. We evaluate PopPoly by applying it to a tetraploid potato cross at nine genomic regions involved in tuber formation.

TriPoly: haplotype estimation for polyploids using sequencing data of related individuals.

Motivation: Knowledge of haplotypes, i.e. phased and ordered marker alleles on a chromosome, is essential to answer many questions in genetics and genomics. By generating short pieces of DNA sequence, high-throughput modern sequencing technologies make estimation of haplotypes possible for single individuals. In polyploids, however, haplotype estimation methods usually require deep coverage to achieve sufficient accuracy. This often renders sequencing-based approaches too costly to be applied to large populations needed in studies of Quantitative Trait Loci. Results: We propose a novel haplotype estimation method for polyploids, TriPoly, that combines sequencing data with Mendelian inheritance rules to infer haplotypes in parent-offspring trios. Using realistic simulations of both short and long-read sequencing data for banana (Musa acuminata) and potato (Solanum tuberosum) trios, we show that TriPoly yields more accurate progeny haplotypes at low coverages compared to existing methods that work on single individuals. We also apply TriPoly to phase Single Nucleotide Polymorphisms on chromosome 5 for a family of tetraploid potato with 2 parents and 37 offspring sequenced with an RNA capture approach. We show that TriPoly haplotype estimates differ from those of the other methods mainly in regions with imperfect sequencing or mapping difficulties, as it does not rely solely on sequence reads and aims to avoid phasings that are not likely to have been passed from the parents to the offspring. Availability and implementation: TriPoly has been implemented in Python 3.5.2 (also compatible with Python 2.7.3 and higher) and can be freely downloaded at https://github.com/EhsanMotazedi/TriPoly. Supplementary information: Supplementary data are available at Bioinformatics online.

Exploiting next-generation sequencing to solve the haplotyping puzzle in polyploids: a simulation study.

Haplotypes are the units of inheritance in an organism, and many genetic analyses depend on their precise determination. Methods for haplotyping single individuals use the phasing information available in next-generation sequencing reads, by matching overlapping single-nucleotide polymorphisms while penalizing post hoc nucleotide corrections made. Haplotyping diploids is relatively easy, but the complexity of the problem increases drastically for polyploid genomes, which are found in both model organisms and in economically relevant plant and animal species. Although a number of tools are available for haplotyping polyploids, the effects of the genomic makeup and the sequencing strategy followed on the accuracy of these methods have hitherto not been thoroughly evaluated.We developed the simulation pipeline haplosim to evaluate the performance of three haplotype estimation algorithms for polyploids: HapCompass, HapTree and SDhaP, in settings varying in sequencing approach, ploidy levels and genomic diversity, using tetraploid potato as the model. Our results show that sequencing depth is the major determinant of haplotype estimation quality, that 1 kb PacBio circular consensus sequencing reads and Illumina reads with large insert-sizes are competitive and that all methods fail to produce good haplotypes when ploidy levels increase. Comparing the three methods, HapTree produces the most accurate estimates, but also consumes the most resources. There is clearly room for improvement in polyploid haplotyping algorithms.

Dysregulated signaling, proliferation and apoptosis impact on the pathogenesis of TCRγδ+ T cell large granular lymphocyte leukemia.

TCRγδ+ T-LGL leukemia is a rare form of chronic mature T cell disorders in elderly, which is generally characterized by a persistently enlarged CD3+CD57+TCRγδ+ large granular lymphocyte population in the peripheral blood with a monoclonal phenotype. Clinically, the disease is heterogeneous, most patients being largely asymptomatic, although neutropenia, fatigue and B symptoms and underlying diseases such as autoimmune diseases or malignancies are also often observed. The etiology of TCRγδ+ T-LGL proliferations is largely unknown. Here, we aimed to investigate underlying molecular mechanisms of these rare proliferations by performing gene expression profiling of TCRγδ+ T-LGL versus normal TCRγδ+ T cell subsets. From our initial microarray dataset we observed that TCRγδ+ T-LGL leukemia forms a separate group when compared with different healthy control TCRγδ+ T cell subsets, correlating best with the healthy TemRA subset. The lowest correlation was seen with the naive subset. Based on specific comparison between healthy control cells and TCRγδ+ T-LGL leukemia cells we observed up-regulation of survival, proliferation and hematopoietic system related genes, with a remarkable down-regulation of apoptotic pathway genes. RQ-PCR validation of important genes representative for the dataset, including apoptosis (XIAP, CASP1, BCLAF1 and CFLAR), proliferation/development (ID3) and inflammation (CD28, CCR7, CX3CR1 and IFNG) processes largely confirmed the dysregulation in proliferation and apoptosis. Based on these expression data we conclude that TCRγδ+ T-LGL leukemia is likely the result of an underlying aberrant molecular mechanisms leading to increased proliferation and reduced apoptosis.

Circulating Human CD27-IgA+ Memory B Cells Recognize Bacteria with Polyreactive Igs.

The vast majority of IgA production occurs in mucosal tissue following T cell-dependent and T cell-independent Ag responses. To study the nature of each of these responses, we analyzed the gene-expression and Ig-reactivity profiles of T cell-dependent CD27(+)IgA(+) and T cell-independent CD27(-)IgA(+) circulating memory B cells. Gene-expression profiles of IgA(+) subsets were highly similar to each other and to IgG(+) memory B cell subsets, with typical upregulation of activation markers and downregulation of inhibitory receptors. However, we identified the mucosa-associated CCR9 and RUNX2 genes to be specifically upregulated in CD27(-)IgA(+) B cells. We also found that CD27(-)IgA(+) B cells expressed Abs with distinct Ig repertoire and reactivity compared with those from CD27(+)IgA(+) B cells. Indeed, Abs from CD27(-)IgA(+) B cells were weakly mutated, often used Igλ chain, and were enriched in polyreactive clones recognizing various bacterial species. Hence, T cell-independent IgA responses are likely involved in the maintenance of gut homeostasis through the production of polyreactive mutated IgA Abs with cross-reactive anti-commensal reactivity.

B-cell prolymphocytic leukemia: a specific subgroup of mantle cell lymphoma.

B-cell prolymphocytic leukemia (B-PLL) is a rare mature B-cell malignancy that may be hard to distinguish from mantle cell lymphoma (MCL) and chronic lymphocytic leukemia (CLL). B-PLL cases with a t(11;14) were redefined as MCL in the World Health Organization 2008 classification. We evaluated 13 B-PLL patients [7 being t(11;14)-positive (B-PLL+) and 6 negative (B-PLL-)] and compared them with MCL and CLL patients. EuroFlow-based immunophenotyping showed significant overlap between B-PLL+ and B-PLL-, as well as between B-PLL and MCL, whereas CLL clustered separately. Immunogenotyping showed specific IGHV gene usage partly resembling MCL. Gene expression profiling showed no separation between B-PLL+ and B-PLL- but identified 3 subgroups. One B-PLL subgroup clustered close to CLL and another subgroup clustered with leukemic MCL; both were associated with prolonged survival. A third subgroup clustered close to nodal MCL and was associated with short survival. Gene expression profiles of both B-PLL+ and B-PLL- showed best resemblance with normal immunoglobulin M-only B-cells. Our data confirm that B-PLL+ is highly comparable to MCL, indicate that B-PLL- also may be considered as a specific subgroup of MCL, and suggest that B-PLL is part of a spectrum, ranging from CLL-like B-PLL, to leukemic MCL-like B-PLL, to nodal MCL-like B-PLL.

Genome duplication and mutations in ACE2 cause multicellular, fast-sedimenting phenotypes in evolved Saccharomyces cerevisiae.

Laboratory evolution of the yeast Saccharomyces cerevisiae in bioreactor batch cultures yielded variants that grow as multicellular, fast-sedimenting clusters. Knowledge of the molecular basis of this phenomenon may contribute to the understanding of natural evolution of multicellularity and to manipulating cell sedimentation in laboratory and industrial applications of S. cerevisiae. Multicellular, fast-sedimenting lineages obtained from a haploid S. cerevisiae strain in two independent evolution experiments were analyzed by whole genome resequencing. The two evolved cell lines showed different frameshift mutations in a stretch of eight adenosines in ACE2, which encodes a transcriptional regulator involved in cell cycle control and mother-daughter cell separation. Introduction of the two ace2 mutant alleles into the haploid parental strain led to slow-sedimenting cell clusters that consisted of just a few cells, thus representing only a partial reconstruction of the evolved phenotype. In addition to single-nucleotide mutations, a whole-genome duplication event had occurred in both evolved multicellular strains. Construction of a diploid reference strain with two mutant ace2 alleles led to complete reconstruction of the multicellular-fast sedimenting phenotype. This study shows that whole-genome duplication and a frameshift mutation in ACE2 are sufficient to generate a fast-sedimenting, multicellular phenotype in S. cerevisiae. The nature of the ace2 mutations and their occurrence in two independent evolution experiments encompassing fewer than 500 generations of selective growth suggest that switching between unicellular and multicellular phenotypes may be relevant for competitiveness of S. cerevisiae in natural environments.

Topology of molecular interaction networks.

Molecular interactions are often represented as network models which have become the common language of many areas of biology. Graphs serve as convenient mathematical representations of network models and have themselves become objects of study. Their topology has been intensively researched over the last decade after evidence was found that they share underlying design principles with many other types of networks.Initial studies suggested that molecular interaction network topology is related to biological function and evolution. However, further whole-network analyses did not lead to a unified view on what this relation may look like, with conclusions highly dependent on the type of molecular interactions considered and the metrics used to study them. It is unclear whether global network topology drives function, as suggested by some researchers, or whether it is simply a byproduct of evolution or even an artefact of representing complex molecular interaction networks as graphs.Nevertheless, network biology has progressed significantly over the last years. We review the literature, focusing on two major developments. First, realizing that molecular interaction networks can be naturally decomposed into subsystems (such as modules and pathways), topology is increasingly studied locally rather than globally. Second, there is a move from a descriptive approach to a predictive one: rather than correlating biological network topology to generic properties such as robustness, it is used to predict specific functions or phenotypes.Taken together, this change in focus from globally descriptive to locally predictive points to new avenues of research. In particular, multi-scale approaches are developments promising to drive the study of molecular interaction networks further.

GATA3 expression is decreased in psoriasis and during epidermal regeneration; induction by narrow-band UVB and IL-4.

Psoriasis is characterized by hyperproliferation of keratinocytes and by infiltration of activated Th1 and Th17 cells in the (epi)dermis. By expression microarray, we previously found the GATA3 transcription factor significantly downregulated in lesional psoriatic skin. Since GATA3 serves as a key switch in both epidermal and T helper cell differentiation, we investigated its function in psoriasis. Because psoriatic skin inflammation shares many characteristics of epidermal regeneration during wound healing, we also studied GATA3 expression under such conditions.Psoriatic lesional skin showed decreased GATA3 mRNA and protein expression compared to non-lesional skin. GATA3 expression was also markedly decreased in inflamed skin of mice with a psoriasiform dermatitis induced with imiquimod. Tape-stripping of non-lesional skin of patients with psoriasis, a standardized psoriasis-triggering and skin regeneration-inducing technique, reduced the expression of GATA3. In wounded skin of mice, low GATA3 mRNA and protein expression was detected. Taken together, GATA3 expression is downregulated under regenerative and inflammatory hyperproliferative skin conditions. GATA3 expression could be re-induced by successful narrow-band UVB treatment of both human psoriasis and imiquimod-induced psoriasiform dermatitis in mice. The prototypic Th2 cytokine IL-4 was the only cytokine capable of inducing GATA3 in skin explants from healthy donors. Based on these findings we argue that GATA3 serves as a key regulator in psoriatic inflammation, keratinocyte hyperproliferation and skin barrier dysfunction.

Effective treatment of psoriasis with narrow-band UVB phototherapy is linked to suppression of the IFN and Th17 pathways.

Narrow-band ultraviolet-B (NB-UVB) phototherapy is an effective treatment for psoriasis. The molecular mechanisms underlying its efficacy are incompletely understood. To identify NB-UVB-induced molecular pathways that may account for its anti-inflammatory efficacy, gene expression profiling was performed using epidermal RNA from lesional and nonlesional skin from patients with psoriasis undergoing NB-UVB therapy. Downregulation of Th17 signaling pathway was observed during NB-UVB therapy in psoriatic epidermis. Strong inhibition of the Th17 pathway by UVB was confirmed in an ex vivo organ culture system by demonstrating reduced signal transducer and activator of transcription 3 (STAT3) phosphorylation and ß-defensin-2 production. These results were further substantiated by demonstrating that NB-UVB inhibited the Th17-dependent psoriasis-like dermatitis in mice. Other pathways affected by NB-UVB therapy include the IFN signaling pathway, epidermal differentiation, and other well-known therapeutic targets in psoriasis, such as the glucocorticoid, vitamin D, peroxisome proliferator-activated receptor, and IL-4 signaling pathways. In conclusion, clinical improvement of psoriasis by NB-UVB is linked to suppression of Th17 and type I and type II IFN signaling pathways, which are critical in the pathogenesis of the disease. Our results show that clinically effective NB-UVB therapy is based on suppression of a broad range of important molecular pathways in psoriatic skin.

Shifts in growth strategies reflect tradeoffs in cellular economics.

The growth rate-dependent regulation of cell size, ribosomal content, and metabolic efficiency follows a common pattern in unicellular organisms: with increasing growth rates, cell size and ribosomal content increase and a shift to energetically inefficient metabolism takes place. The latter two phenomena are also observed in fast growing tumour cells and cell lines. These patterns suggest a fundamental principle of design. In biology such designs can often be understood as the result of the optimization of fitness. Here we show that in basic models of self-replicating systems these patterns are the consequence of maximizing the growth rate. Whereas most models of cellular growth consider a part of physiology, for instance only metabolism, the approach presented here integrates several subsystems to a complete self-replicating system. Such models can yield fundamentally different optimal strategies. In particular, it is shown how the shift in metabolic efficiency originates from a tradeoff between investments in enzyme synthesis and metabolic yields for alternative catabolic pathways. The models elucidate how the optimization of growth by natural selection shapes growth strategies.

DNA microarray studies of hematopoietic subpopulations.

In recent years, one of the most quickly incorporated methods in biomedical research has been microarray technology. Microarrays have been designed for the purpose of genotyping (e.g., SNP analysis), expression analysis (mRNA, miRNA, exon arrays), chromatin immunoprecipitations (ChIP-on-chip technology), and DNA sequencing. In this chapter we will focus on the application of DNA microarrays in gene expression analysis of mRNA. This technology allows for the simultaneous analysis of transcription patterns of literally thousands of genes, making it a very powerful approach. Hematopoietic sub-populations are relatively easy to obtain by using fluorescence-activated cell sorting or magnetic bead cell separation. Over the last decade, the combination of these purification techniques with microarray analysis has resulted in an enormous gain of knowledge of blood lineages and their development.

Insertional mutagenesis combined with acquired somatic mutations causes leukemogenesis following gene therapy of SCID-X1 patients.

X-linked SCID (SCID-X1) is amenable to correction by gene therapy using conventional gammaretroviral vectors. Here, we describe the occurrence of clonal T cell acute lymphoblastic leukemia (T-ALL) promoted by insertional mutagenesis in a completed gene therapy trial of 10 SCID-X1 patients. Integration of the vector in an antisense orientation 35 kb upstream of the protooncogene LIM domain only 2 (LMO2) caused overexpression of LMO2 in the leukemic clone. However, leukemogenesis was likely precipitated by the acquisition of other genetic abnormalities unrelated to vector insertion, including a gain-of-function mutation in NOTCH1, deletion of the tumor suppressor gene locus cyclin-dependent kinase 2A (CDKN2A), and translocation of the TCR-beta region to the STIL-TAL1 locus. These findings highlight a general toxicity of endogenous gammaretroviral enhancer elements and also identify a combinatorial process during leukemic evolution that will be important for risk stratification and for future protocol design.

Vector integration is nonrandom and clustered and influences the fate of lymphopoiesis in SCID-X1 gene therapy.

Recent reports have challenged the notion that retroviruses and retroviral vectors integrate randomly into the host genome. These reports pointed to a strong bias toward integration in and near gene coding regions and, for gammaretroviral vectors, around transcription start sites. Here, we report the results obtained from a large-scale mapping of 572 retroviral integration sites (RISs) isolated from cells of 9 patients with X-linked SCID (SCID-X1) treated with a retrovirus-based gene therapy protocol. Our data showed that two-thirds of insertions occurred in or very near to genes, of which more than half were highly expressed in CD34(+) progenitor cells. Strikingly, one-fourth of all integrations were clustered as common integration sites (CISs). The highly significant incidence of CISs in circulating T cells and the nature of their locations indicate that insertion in many gene loci has an influence on cell engraftment, survival, and proliferation. Beyond the observed cases of insertional mutagenesis in 3 patients, these data help to elucidate the relationship between vector insertion and long-term in vivo selection of transduced cells in human patients with SCID-X1.

Gammaretrovirus-mediated correction of SCID-X1 is associated with skewed vector integration site distribution in vivo.

We treated 10 children with X-linked SCID (SCID-X1) using gammaretrovirus-mediated gene transfer. Those with sufficient follow-up were found to have recovered substantial immunity in the absence of any serious adverse events up to 5 years after treatment. To determine the influence of vector integration on lymphoid reconstitution, we compared retroviral integration sites (RISs) from peripheral blood CD3(+) T lymphocytes of 5 patients taken between 9 and 30 months after transplantation with transduced CD34(+) progenitor cells derived from 1 further patient and 1 healthy donor. Integration occurred preferentially in gene regions on either side of transcription start sites, was clustered, and correlated with the expression level in CD34(+) progenitors during transduction. In contrast to those in CD34(+) cells, RISs recovered from engrafted CD3(+) T cells were significantly overrepresented within or near genes encoding proteins with kinase or transferase activity or involved in phosphorus metabolism. Although gross patterns of gene expression were unchanged in transduced cells, the divergence of RIS target frequency between transduced progenitor cells and post-thymic T lymphocytes indicates that vector integration influences cell survival, engraftment, or proliferation.

Retroviral vector insertion sites associated with dominant hematopoietic clones mark "stemness" pathways.

Evidence from model organisms and clinical trials reveals that the random insertion of retrovirus-based vectors in the genome of long-term repopulating hematopoietic cells may increase self-renewal or initiate malignant transformation. Clonal dominance of nonmalignant cells is a particularly interesting phenotype as it may be caused by the dysregulation of genes that affect self-renewal and competitive fitness. We have accumulated 280 retrovirus vector insertion sites (RVISs) from murine long-term studies resulting in benign or malignant clonal dominance. RVISs (22.5%) are located in or near (up to 100 kb [kilobase]) to known proto-oncogenes, 49.6% in signaling genes, and 27.9% in other or unknown genes. The resulting insertional dominance database (IDDb) shows substantial overlaps with the transcriptome of hematopoietic stem/progenitor cells and the retrovirus-tagged cancer gene database (RTCGD). RVISs preferentially marked genes with high expression in hematopoietic stem/progenitor cells, and Gene Ontology revealed an overrepresentation of genes associated with cell-cycle control, apoptosis signaling, and transcriptional regulation, including major "stemness" pathways. The IDDb forms a powerful resource for the identification of genes that stimulate or transform hematopoietic stem/progenitor cells and is an important reference for vector biosafety studies in human gene therapy.

Gene therapy: is IL2RG oncogenic in T-cell development?

The gene IL2RG encodes the gamma-chain of the interleukin-2 receptor and is mutated in patients with X-linked severe combined immune deficiency (X-SCID). Woods et al. report the development of thymus tumours in a mouse model of X-SCID after correction by lentiviral overexpression of IL2RG and claim that these were caused by IL2RG itself. Here we find that retroviral overexpression of IL2RG in human CD34+ cells has no effect on T-cell development, whereas overexpression of the T-cell acute lymphoblastic leukaemia (T-ALL) oncogene LMO2 leads to severe abnormalities. Retroviral expression of IL2RG may therefore not be directly oncogenic--rather, the restoration of normal signalling by the interleukin-7 receptor to X-SCID precursor cells allows progression of T-cell development to stages that are permissive for the pro-leukaemic effects of ectopic LMO2.