Survey of sugar beet (Beta vulgaris L.) hAT transposons and MITE-like hATpin derivatives

Genome-wide analyses of repetitive DNA suggest a significant impact particularly of transposable elements on genome size and evolution of virtually all eukaryotic organisms. In this study, we analyzed the abundance and diversity of the hAT transposon superfamily of the sugar beet (B. vulgaris) genome, using molecular, bioinformatic and cytogenetic approaches. We identified 81 transposase-coding sequences, three of which are part of structurally intact but nonfunctional hAT transposons (BvhAT), in a B. vulgaris BAC library as well as in whole genome sequencing-derived data sets. Additionally, 116 complete and 497 truncated non-autonomous BvhAT derivatives lacking the transposase gene were in silico-detected. The 116 complete derivatives were subdivided into four BvhATpin groups each characterized by a distinct terminal inverted repeat motif. Both BvhAT and BvhATpin transposons are specific for species of the genus Beta and closely related species, showing a localization on B. vulgaris chromosomes predominantely in euchromatic regions. The lack of any BvhAT transposase function together with the high degree of degeneration observed for the BvhAT and the BvhATpin genomic fraction contrasts with the abundance and activity of autonomous and non-autonomous hAT transposons revealed in other plant species. This indicates a possible genus-specific structural and functional repression of the hAT transposon superfamily during Beta diversification and evolution.

Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia

Here we perform whole-exome sequencing of samples from 105 individuals with chronic lymphocytic leukemia (CLL), the most frequent leukemia in adults in Western countries. We found 1,246 somatic mutations potentially affecting gene function and identified 78 genes with predicted functional alterations in more than one tumor sample. Among these genes, SF3B1, encoding a subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP), is somatically mutated in 9.7% of affected individuals. Further analysis in 279 individuals with CLL showed that SF3B1 mutations were associated with faster disease progression and poor overall survival. This work provides the first comprehensive catalog of somatic mutations in CLL with relevant clinical correlates and defines a large set of new genes that may drive the development of this common form of leukemia. The results reinforce the idea that targeting several well-known genetic pathways, including mRNA splicing, could be useful in the treatment of CLL and other malignancies.

Role of medium- and short-chain L-3-hydroxyacyl-CoA dehydrogenase in the regulation of body weight and thermogenesis

Dysregulation of fatty acid oxidation plays a pivotal role in the pathophysiology of obesity and insulin resistance. Medium- and short-chain-3-hydroxyacyl-coenzyme A (CoA) dehydrogenase (SCHAD) (gene name, hadh) catalyze the third reaction of the mitochondrial β-oxidation cascade, the oxidation of 3-hydroxyacyl-CoA to 3-ketoacyl-CoA, for medium- and short-chain fatty acids. We identified hadh as a putative obesity gene by comparison of two genome-wide scans, a quantitative trait locus analysis previously performed in the polygenic obese New Zealand obese mouse and an earlier described small interfering RNA-mediated mutagenesis in Caenorhabditis elegans. In the present study, we show that mice lacking SCHAD (hadh(-/-)) displayed a lower body weight and a reduced fat mass in comparison with hadh(+/+) mice under high-fat diet conditions, presumably due to an impaired fuel efficiency, the loss of acylcarnitines via the urine, and increased body temperature. Food intake, total energy expenditure, and locomotor activity were not altered in knockout mice. Hadh(-/-) mice exhibited normal fat tolerance at 20 C. However, during cold exposure, knockout mice were unable to clear triglycerides from the plasma and to maintain their normal body temperature, indicating that SCHAD plays an important role in adaptive thermogenesis. Blood glucose concentrations in the fasted and postprandial state were significantly lower in hadh(-/-) mice, whereas insulin levels were elevated. Accordingly, insulin secretion in response to glucose and glucose plus palmitate was elevated in isolated islets of knockout mice. Therefore, our data indicate that SCHAD is involved in thermogenesis, in the maintenance of body weight, and in the regulation of nutrient-stimulated insulin secretion.

Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and genome analyzer systems

BACKGROUND

The generation and analysis of high-throughput sequencing data are becoming a major component of many studies in molecular biology and medical research. Illumina's Genome Analyzer (GA) and HiSeq instruments are currently the most widely used sequencing devices. Here, we comprehensively evaluate properties of genomic HiSeq and GAIIx data derived from two plant genomes and one virus, with read lengths of 95 to 150 bases.

RESULTS

We provide quantifications and evidence for GC bias, error rates, error sequence context, effects of quality filtering, and the reliability of quality values. By combining different filtering criteria we reduced error rates 7-fold at the expense of discarding 12.5% of alignable bases. While overall error rates are low in HiSeq data we observed regions of accumulated wrong base calls. Only 3% of all error positions accounted for 24.7% of all substitution errors. Analyzing the forward and reverse strands separately revealed error rates of up to 18.7%. Insertions and deletions occurred at very low rates on average but increased to up to 2% in homopolymers. A positive correlation between read coverage and GC content was found depending on the GC content range.

CONCLUSIONS

The errors and biases we report have implications for the use and the interpretation of Illumina sequencing data. GAIIx and HiSeq data sets show slightly different error profiles. Quality filtering is essential to minimize downstream analysis artifacts. Supporting previous recommendations, the strand-specificity provides a criterion to distinguish sequencing errors from low abundance polymorphisms.

Partial short-read sequencing of a highly inbred Iberian pig and genomics inference thereof

Despite dramatic reduction in sequencing costs with the advent of next generation sequencing technologies, obtaining a complete mammalian genome sequence at sufficient depth is still costly. An alternative is partial sequencing. Here, we have sequenced a reduced representation library of an Iberian sow from the Guadyerbas strain, a highly inbred strain that has been used in numerous QTL studies because of its extreme phenotypic characteristics. Using the Illumina Genome Analyzer II (San Diego, CA, USA), we resequenced ∼ 1% of the genome with average 4 × depth, identifying 68,778 polymorphisms. Of these, 55,457 were putative fixed differences with respect to the assembly, based on the genome of a Duroc pig, and 13,321 were heterozygous positions within Guadyerbas. Despite being highly inbred, the estimate of heterozygosity within Guadyerbas was ∼ 0.78 kb(-1) in autosomes, after correcting for low depth. Nucleotide variability was consistently higher at the telomeric regions than on the rest of the chromosome, likely a result of increased recombination rates. Further, variability was 50% lower in the X-chromosome than in autosomes, which may be explained by a recent bottleneck or by selection. We divided the whole genome in 500 kb windows and we analyzed overrepresented gene ontology terms in regions of low and high variability. Multi organism process, pigmentation and cell killing were overrepresented in high variability regions and metabolic process ontology, within low variability regions. Further, a genome wide Hudson-Kreitman-Aguadé test was carried out per window; overall, variability was in agreement with neutral expectations.

p53 Gene repair with zinc finger nucleases optimised by yeast 1-hybrid and validated by Solexa sequencing

The tumor suppressor gene p53 is mutated or deleted in over 50% of human tumors. As functional p53 plays a pivotal role in protecting against cancer development, several strategies for restoring wild-type (wt) p53 function have been investigated. In this study, we applied an approach using gene repair with zinc finger nucleases (ZFNs). We adapted a commercially-available yeast one-hybrid (Y1H) selection kit to allow rapid building and optimization of 4-finger constructs from randomized PCR libraries. We thus generated novel functional zinc finger nucleases against two DNA sites in the human p53 gene, near cancer mutation ‘hotspots’. The ZFNs were first validated using in vitro cleavage assays and in vivo episomal gene repair assays in HEK293T cells. Subsequently, the ZFNs were used to restore wt-p53 status in the SF268 human cancer cell line, via ZFN-induced homologous recombination. The frequency of gene repair and mutation by non-homologous end-joining was then ascertained in several cancer cell lines, using a deep sequencing strategy. Our Y1H system facilitates the generation and optimisation of novel, sequence-specific four- to six-finger peptides, and the p53-specific ZFN described here can be used to mutate or repair p53 in genomic loci.

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in cancer.

Expression analysis of proline rich 15 (Prr15) in mouse and human gastrointestinal tumors

Proline rich 15 (Prr15), which encodes a protein of unknown function, is expressed almost exclusively in postmitotic cells both during fetal development and in adult tissues, such as the intestinal epithelium and the testis. To determine if this specific expression is lost in intestinal neoplasias, we examined Prr15 expression by in situ hybridization (ISH) on mouse intestinal tumors caused by different gene mutations, and on human colorectal cancer (CRC) samples. Prr15/PRR15 expression was consistently observed in mouse gastrointestinal (GI) tumors caused by mutations in the Apc gene, as well as in several advanced stage human CRCs. In contrast, no Prr15 expression was detected in intestinal tumors derived from mice carrying mutations in the Smad3, Smad4, or Cdkn1b genes. These findings, combined with the fact that a majority of sporadic human CRCs carry APC mutations, strongly suggest that the expression of Prr15/PRR15 in mouse and human GI tumors is linked, directly or indirectly, to the absence of the APC protein or, more generally, to the disruption of the Wnt signaling pathway.

High-throughput identification of genetic markers using representational oligonucleotide microarray analysis

We describe a novel approach for high-throughput development of genetic markers using representational oligonucleotide microarray analysis. We test the performance of the method in sugar beet (Beta vulgaris L.) as a model for crop plants with little sequence information available. Genomic representations of both parents of a mapping population were hybridized on microarrays containing in total 146,554 custom made oligonucleotides based on sugar beet bacterial artificial chromosome (BAC) end sequences and expressed sequence tags (ESTs). Oligonucleotides showing a signal with one parental line only, were selected as potential marker candidates and placed onto an array, designed for genotyping of 184 F(2) individuals from the mapping population. Utilizing known co-dominant anchor markers we obtained 511 new dominant markers (392 derived from BAC end sequences, and 119 from ESTs) distributed over all nine sugar beet linkage groups and calculated genetic maps. Further improvements for large-scale application of the approach are discussed and its feasibility for the cost-effective and flexible generation of genetic markers is presented.

Strand-specific deep sequencing of the transcriptome

Several studies support that antisense-mediated regulation may affect a large proportion of genes. Using the Illumina next-generation sequencing platform, we developed DSSS (direct strand specific sequencing), a strand-specific protocol for transcriptome sequencing. We tested DSSS with RNA from two samples, prokaryotic (Mycoplasma pneumoniae) as well as eukaryotic (Mus musculus), and obtained data containing strand-specific information, using single-read and paired-end sequencing. We validated our results by comparison with a strand-specific tiling array data set for strain M129 of the simple prokaryote M. pneumoniae, and by quantitative PCR (qPCR). The results of DSSS were very well supported by the results from tiling arrays and qPCR. Moreover, DSSS provided higher dynamic range and single-base resolution, thus enabling efficient antisense detection and the precise mapping of transcription start sites and untranslated regions. DSSS data for mouse confirmed strand specificity of the protocol and the general applicability of the approach to studying eukaryotic transcription. We propose DSSS as a simple and efficient strategy for strand-specific transcriptome sequencing and as a tool for genome annotation exploiting the increased read lengths that next-generation sequencing technology now is capable to deliver.

Global transcriptomic analysis of murine embryonic stem cell-derived brachyury(+) (T) cells

Brachyury(+) mesodermal cell population with purity over 79% was obtained from differentiating brachyury embryonic stem cells (ESC) generated with brachyury promoter driven enhanced green fluorescent protein and puromycin-N-acetyltransferase. A comprehensive transcriptomic analysis of brachyury(+) cells enriched with puromycin application from 6-day-old embryoid bodies (EBs), 6-day-old control EBs and undifferentiated ESCs led to identification of 1573 uniquely up-regulated and 1549 uniquely down-regulated transcripts in brachyury(+) cells. Furthermore, transcripts up-regulated in brachyury(+) cells have overrepresented the Gene Ontology annotations (cell differentiation, blood vessel morphogenesis, striated muscle development, placenta development and cell motility) and Kyoto Encyclopedia of Genes and Genomes pathway annotations (mitogen-activated protein kinase signaling and transforming growth factor beta signaling). Transcripts representing Larp2 and Ankrd34b are notably up-regulated in brachyury(+) cells. Knockdown of Larp2 resulted in a significantly down-regulation BMP-2 expression, and knockdown of Ankrd34b resulted in alteration of NF-H, PPARγ and PECAM1 expression. The elucidation of transcriptomic signatures of ESCs-derived brachyury(+) cells will contribute toward defining the genetic and cellular identities of presumptive mesodermal cells. Furthermore, there is a possible involvement of Larp2 in the regulation of the late mesodermal marker BMP-2. Ankrd34b might be a positive regulator of neurogenesis and a negative regulator of adipogenesis.

Cloning of mouse ojoplano, a reticular cytoplasmic protein expressed during embryonic development

Ojoplano (Opo) is a morphogenetic gene playing an important role during embryogenesis in medaka. This report focuses on the identification and characterization of the mouse Opo gene. We examined Opo expression by whole-mount in situ hybridization and in situ hybridization on sagittal sections during mouse embryogenesis. First expression in whole-mounts was detected at Theiler stages 15-17 (E 9.5-10.5dpc) as a spotted specific staining in migrating neural crest cells and in placodal structures. A complex expression pattern was observed in Theiler stage 22-23 (E 14.5dpc) in sagittal sections, including expression in skeletal structures (skull, vertebrae, ribs, bones of the locomotor system), in the nasal region, the heart and the eye. Fusion proteins revealed the localization of OPO within the cytoplasm with a reticular distribution that largely overlapped with the endoplasmic reticulum. Opo shows homology to human transcripts linked to a hereditary craniofacial malformation, orofacial cleft 1 (OFC1). The expression of mouse Opo in neural crest derivatives and skull elements further supports this link.

The FunGenES database: a genomics resource for mouse embryonic stem cell differentiation

Embryonic stem (ES) cells have high self-renewal capacity and the potential to differentiate into a large variety of cell types. To investigate gene networks operating in pluripotent ES cells and their derivatives, the “Functional Genomics in Embryonic Stem Cells” consortium (FunGenES) has analyzed the transcriptome of mouse ES cells in eleven diverse settings representing sixty-seven experimental conditions. To better illustrate gene expression profiles in mouse ES cells, we have organized the results in an interactive database with a number of features and tools. Specifically, we have generated clusters of transcripts that behave the same way under the entire spectrum of the sixty-seven experimental conditions; we have assembled genes in groups according to their time of expression during successive days of ES cell differentiation; we have included expression profiles of specific gene classes such as transcription regulatory factors and Expressed Sequence Tags; transcripts have been arranged in “Expression Waves” and juxtaposed to genes with opposite or complementary expression patterns; we have designed search engines to display the expression profile of any transcript during ES cell differentiation; gene expression data have been organized in animated graphs of KEGG signaling and metabolic pathways; and finally, we have incorporated advanced functional annotations for individual genes or gene clusters of interest and links to microarray and genomic resources. The FunGenES database provides a comprehensive resource for studies into the biology of ES cells.

ojoplano-mediated basal constriction is essential for optic cup morphogenesis

Although the vertebrate retina is a well-studied paradigm for organogenesis, the morphogenetic mechanisms that carve the architecture of the vertebrate optic cup remain largely unknown. Understanding how the hemispheric shape of an eye is formed requires addressing the fundamental problem of how individual cell behaviour is coordinated to direct epithelial morphogenesis. Here, we analyze the role of ojoplano (opo), an uncharacterized gene whose human ortholog is associated with orofacial clefting syndrome, in the morphogenesis of epithelial tissues. Most notably,when opo is mutated in medaka fish, optic cup folding is impaired. We characterize optic cup morphogenesis in vivo and determine at the cellular level how opo affects this process. opo encodes a developmentally regulated transmembrane protein that localizes to compartments of the secretory pathway and to basal end-feet of the neuroepithelial precursors. We show that Opo regulates the polarized localization of focal adhesion components to the basal cell surface. Furthermore, tissue-specific interference with integrin-adhesive function impairs optic cup folding,resembling the ocular phenotype observed in opo mutants. We propose a model of retinal morphogenesis whereby opo-mediated formation of focal contacts is required to transmit the mechanical tensions that drive the macroscopic folding of the vertebrate optic cup.

Construction and characterization of a sugar beet (Beta vulgaris) fosmid library

A sugar beet (Beta vulgaris) fosmid library from the doubled haploid accession KWS2320 encompassing 115 200 independent clones was constructed and characterized. The average insert size of the fosmid library was determined by pulsed field gel electrophoresis to be 39 kbp on average, thus representing 5.9-fold coverage of the sugar beet genome (758 Mbp). PCR screening of plate pools with primer pairs against nine sugar beet genes supported the insert size estimation. BLAST searches with 2951 fosmid end-sequences originating from 1510 clones (1536 clones attempted) revealed little contamination with organellar DNA (2.1% chloroplast DNA, 0.3% mitochondrial DNA). The sugar beet fosmid library will be integrated in the presently ongoing efforts to determine the sequence of the sugar beet genome. Fosmids will be publicly available in the format of plate pools and individual clones.

Haplotype divergence in Beta vulgaris and microsynteny with sequenced plant genomes

We characterized two overlapping sugar beet (Beta vulgaris) bacterial artificial chromosome (BAC) clones representing different haplotypes. A total of 254 kbp of the genomic sequence was determined, of which the two BACs share 92 kbp. Eleven of 15 genes discovered in the sequenced interval locate to the overlap region. The haplotypes differ in exons by 1% (nucleotide level) and in non-coding regions by 9% (6% mismatches, 3% gaps; alignable regions only). Large indels or high sequence divergence comprised 11% of either sequence. Of such indels, 68 and 45%, respectively, could be attributed to haplotype-specific integration of transposable elements. We identified novel repeat candidates by comparing the two BAC sequences to a set of genomic sugar beet sequences. Synteny was found with Arabidopsis chromosome 1 (At1), At2 and At4, Medicago chromosome 7, Vitis chromosome 15 and paralogous regions on poplar chromosomes II and XIV.

Neuronal functions, feeding behavior, and energy balance in Slc2a3+/- mice

Homozygous deletion of the gene of the neuronal glucose transporter GLUT3 (Slc2a3) in mice results in embryonic lethality, whereas heterozygotes (Slc2a3+/-) are viable. Here, we describe the characterization of heterozygous mice with regard to neuronal function, glucose homeostasis, and, since GLUT3 might be a component of the neuronal glucose-sensing mechanism, food intake and energy balance. Levels of GLUT3 mRNA and protein in brain were reduced by 50% in Slc2a3+/- mice. Electrographic features examined by electroencephalographic recordings give evidence for slightly but significantly enhanced cerebrocortical activity in Slc2a3+/- mice. In addition, Slc2a3+/- mice were slightly more sensitive to an acoustic startle stimulus (elevated startle amplitude and reduced prepulse inhibition). However, systemic behavioral testing revealed no other functional abnormalities, e.g., in coordination, reflexes, motor abilities, anxiety, learning, and memory. Furthermore, no differences in body weight, blood glucose, and insulin levels were detected between wild-type and Slc2a3+/- littermates. Food intake as monitored randomly or after intracerebroventricular administration of 2-deoxyglucose or d-glucose, or food choice for carbohydrates/fat was not affected in Slc2a3+/- mice. Taken together, our data indicate that, in contrast to Slc2a1, a single allele of Slc2a3 is sufficient for maintenance of neuronal energy supply, motor abilities, learning and memory, and feeding behavior.

WDR55 is a nucleolar modulator of ribosomal RNA synthesis, cell cycle progression, and teleost organ development

The thymus is a vertebrate-specific organ where T lymphocytes are generated. Genetic programs that lead to thymus development are incompletely understood. We previously screened ethylnitrosourea-induced medaka mutants for recessive defects in thymus development. Here we report that one of those mutants is caused by a missense mutation in a gene encoding the previously uncharacterized protein WDR55 carrying the tryptophan-aspartate-repeat motif. We find that WDR55 is a novel nucleolar protein involved in the production of ribosomal RNA (rRNA). Defects in WDR55 cause aberrant accumulation of rRNA intermediates and cell cycle arrest. A mutation in WDR55 in zebrafish also leads to analogous defects in thymus development, whereas WDR55-null mice are lethal before implantation. These results indicate that WDR55 is a nuclear modulator of rRNA synthesis, cell cycle progression, and embryonic organogenesis including teleost thymus development.

Higher-order genome organization in platypus and chicken sperm and repositioning of sex chromosomes during mammalian evolution

In mammals, chromosomes occupy defined positions in sperm, whereas previous work in chicken showed random chromosome distribution. Monotremes (platypus and echidnas) are the most basal group of living mammals. They have elongated sperm like chicken and a complex sex chromosome system with homology to chicken sex chromosomes. We used platypus and chicken genomic clones to investigate genome organization in sperm. In chicken sperm, about half of the chromosomes investigated are organized non-randomly, whereas in platypus chromosome organization in sperm is almost entirely non-random. The use of genomic clones allowed us to determine chromosome orientation and chromatin compaction in sperm. We found that in both species chromosomes maintain orientation of chromosomes in sperm independent of random or non-random positioning along the sperm nucleus. The distance of loci correlated with the total length of sperm nuclei, suggesting that chromatin extension depends on sperm elongation. In platypus, most sex chromosomes cluster in the posterior region of the sperm nucleus, presumably the result of postmeiotic association of sex chromosomes. Chicken and platypus autosomes sharing homology with the human X chromosome located centrally in both species suggesting that this is the ancestral position. This suggests that in some therian mammals a more anterior position of the X chromosome has evolved independently.