Search

Long non-coding RNAs (lncRNAs) are transcripts with a length of more than 200 nucleotides that lack protein-coding capacity. LncRNAs play important roles in the regulation of genes during many biological processes, development and disease progression, such as obesity. Obesity is an increasing health concern around the world. Although multiple studies linking lncRNAs and fully differentiated adipocytes have been published, the systematic analysis of those lncRNAs involved in early preadipocyte differentiation, when fate determination decisions are made, has not been reported. In order to fill this gap, we conducted strand-specific RNA-Seq on mouse 3T3-L1 preadipocyte cells and compared the expression profiles before and after early differentiation for 2 days. We identified 82 lncRNAs that significantly changed their expression after early differentiation, 98% of which were newly discovered in association with early adipogenesis. The most remarkable lncRNAs were U90926, Wincr1, Kcnq1ot1, Malat1 and Hotairm1. Expression patterns for these identified genes were also highly correlated with markers of adipogenesis, including PPARγ, CEBPα, FABP4 and FASN. Further analysis implies that many of those altered lncRNAs might coordinately inhibit Wnt/β-catenin signaling pathway to promote preadipocyte differentiation. This work strongly suggests that lncRNAs are critical in the proliferation and differentiation of adipose tissue.

Exosomes are paracrine regulators of the tumor microenvironment and contain complex cargo. We previously reported that exosomes released from acute myeloid leukemia (AML) cells can suppress residual hematopoietic stem and progenitor cell (HSPC) function indirectly through stromal reprogramming of niche retention factors. We found that the systemic loss of hematopoietic function is also in part a consequence of AML exosome-directed microRNA (miRNA) trafficking to HSPCs. Exosomes isolated from cultured AML or the plasma from mice bearing AML xenografts exhibited enrichment of miR-150 and miR-155. HSPCs cocultured with either of these exosomes exhibited impaired clonogenicity, through the miR-150- and miR-155-mediated suppression of the translation of transcripts encoding c-MYB, a transcription factor involved in HSPC differentiation and proliferation. To discover additional miRNA targets, we captured miR-155 and its target transcripts by coimmunoprecipitation with an attenuated RNA-induced silencing complex (RISC)-trap, followed by high-throughput sequencing. This approach identified known and previously unknown miR-155 target transcripts. Integration of the miR-155 targets with information from the protein interaction database STRING revealed proteins indirectly affected by AML exosome-derived miRNA. Our findings indicate a direct effect of AML exosomes on HSPCs that, through a stroma-independent mechanism, compromises hematopoiesis. Furthermore, combining miRNA target data with protein-protein interaction data may be a broadly applicable strategy to define the effects of exosome-mediated trafficking of regulatory molecules within the tumor microenvironment.

Identifying targets is critical for understanding the biological effects of microRNA (miRNA) expression. The challenge lies in characterizing the cohort of targets for a specific miRNA, especially when targets are being actively down-regulated in miRNA– RNA-induced silencing complex (RISC)–messengerRNA (mRNA) complexes. We have developed a robust and versatile strategy called RISCtrap to stabilize and purify targets from this transient interaction. Its utility was demonstrated by determining specific high-confidence target datasets for miR-124, miR-132, and miR-181 that contained known and previously unknown transcripts. Two previously unknown miR-132 targets identified with RISCtrap, adaptor protein CT10 regulator of kinase 1 (CRK1) and tight junction-associated protein 1 (TJAP1), were shown to be endogenously regulated by miR-132 in adult mouse forebrain. The datasets, moreover, differed in the number of targets and in the types and frequency of microRNA recognition element (MRE) motifs, thus revealing a previously underappreciated level of specificity in the target sets regulated by individual miRNAs.

Combinatorial transcription codes generate the myriad of cell types during development and thus likely provide crucial insights into directed differentiation of stem cells to a specific cell type. The LIM complex composed of Isl1 and Lhx3 directs the specification of spinal motor neurons (MNs) in embryos. Here, we report that Isl1–Lhx3, a LIM-complex mimicking fusion, induces a signature of MN transcriptome and concomitantly suppresses interneuron differentiation programs, thereby serving as a potent and specific inducer of MNs in stem cells. We show that an equimolar ratio of Isl1 and Lhx3 and the LIM domain of Lhx3 are crucial for generating MNs without up-regulating interneuron genes. These led us to design Isl1–Lhx3, which maintains the desirable 1:1 ratio of Isl1 and Lhx3 and the LIM domain of Lhx3. Isl1–Lhx3 drives MN differentiation with high specificity and efficiency in the spinal cord and embryonic stem cells, bypassing the need for sonic hedgehog (Shh). RNA-seq analysis revealed that Isl1–Lhx3 induces the expression of a battery of MN genes that control various functional aspects of MNs, while suppressing key interneuron genes. Our studies uncover a highly efficient method for directed MN generation and MN gene networks. Our results also demonstrate a general strategy of using embryonic transcription complexes for producing specific cell types from stem cells.

Polymorphisms associated with gene coding regions are useful tools for use as molecular markers, for example in marker-assisted breeding efforts. In particular, single nucleotide polymorphisms (SNPs) have emerged as a preferred marker for high-throughput genotyping studies. The objective of this study was to develop a set of SNPs for European hazelnut (Corylus avellana) by exploiting the exceptional depth and breadth of transcriptome sequencing made possible by high-throughput RNA sequencing (RNA-seq) using the Illumina Genome Analyzer platform. Bioinformatics tools were used to mine SNPs from a database of EST sequences derived from de novo assembly of Illumina RNA-seq reads representing hazelnut transcripts from multiple accessions. This resulted in the identification of 5,398 SNPs. This set of SNP markers for hazelnut provides a new resource for genetic mapping of important agronomic traits and will be useful markers for aligning future physical and genetic maps.

Alternative splicing can enhance transcriptome plasticity and proteome diversity. In plants, alternative splicing can be manifested at different developmental stages, and is frequently associated with specific tissue types or environmental conditions such as abiotic stress. We mapped the Arabidopsis transcriptome at single-base resolution using the Illumina platform for ultrahigh-throughput RNA sequencing (RNA-seq). Deep transcriptome sequencing confirmed a majority of annotated introns and identified thousands of novel alternatively spliced mRNA isoforms. Our analysis suggests that at least ∼42% of intron-containing genes in Arabidopsis are alternatively spliced; this is significantly higher than previous estimates based on cDNA/expressed sequence tag sequencing. Random validation confirmed that novel splice isoforms empirically predicted by RNA-seq can be detected in vivo. Novel introns detected by RNA-seq were substantially enriched in nonconsensus terminal dinucleotide splice signals. Alternative isoforms with premature termination codons (PTCs) comprised the majority of alternatively spliced transcripts. Using an example of an essential circadian clock gene, we show that intron retention can generate relatively abundant PTC+ isoforms and that this specific event is highly conserved among diverse plant species. Alternatively spliced PTC+ isoforms can be potentially targeted for degradation by the nonsense mediated mRNA decay (NMD) surveillance machinery or regulate the level of functional transcripts by the mechanism of regulated unproductive splicing and translation (RUST). We demonstrate that the relative ratios of the PTC+ and reference isoforms for several key regulatory genes can be considerably shifted under abiotic stress treatments. Taken together, our results suggest that like in animals, NMD and RUST may be widespread in plants and may play important roles in regulating gene expression.

While microRNAs have emerged as an important component of gene regulatory networks, it remains unclear how microRNAs collaborate with transcription factors in the gene networks that determines neuronal cell fate. Here we show that in the developing spinal cord, the expression of miR-218 is directly upregulated by the Isl1-Lhx3 complex, which drives motor neuron fate. Inhibition of miR-218 suppresses the generation of motor neurons in both chick neural tube and mouse embryonic stem cells, suggesting that miR-218 plays a crucial role in motor neuron differentiation. Results from unbiased RISC-trap screens, in vivo reporter assays and overexpression studies indicated that miR-218 directly represses transcripts that promote developmental programs for interneurons. In addition, we found that miR-218 activity is required for Isl1-Lhx3 to effectively induce motor neurons and suppress interneuron fates. Together our results reveal an essential role of miR-218 as a downstream effector of the Isl1-Lhx3 complex in establishing motor neuron identity.

Our laboratory has reported data suggesting a role for T4 phage gene 32 single-stranded DNA-binding protein in organizing a complex of deoxyribonucleotide-synthesizing enzymes at the replication fork. In this article we examined the effects of gene 32 ablation on the association of these enzymes with DNA-protein complexes. These experiments showed several deoxyribonucleotide-synthesizing enzymes to be present in DNA-protein complexes, with some of these associations being dependent on gene 32 protein. To further understand the role of gp32, we created amber mutations at codons 24 and 204 of gene 32, which encodes a 301-residue protein. We used the newly created mutants along with several experimental approaches--DNA-cellulose chromatography, immunoprecipitation, optical biosensor analysis and glutathione-S-transferase pulldowns--to identify relevant protein-protein and protein-DNA interactions. These experiments identified several proteins whose interactions with DNA depend on the presence of intact gp32, notably thymidylate synthase, dihydrofolate (DHF) reductase, ribonucleotide reductase (RNR) and Escherichia coli nucleoside diphosphate (NDP) kinase, and they also demonstrated direct associations between gp32 and RNR and NDP kinase, but not dCMP hydroxymethylase, deoxyribonucleoside monophosphate kinase, or DHF reductase. Taken together, the results support the hypothesis that the gene 32 protein helps to recruit enzymes of deoxyribonucleoside triphosphates synthesis to DNA replication sites.

The histone H3-lysine-4 methyltransferase mixed-lineage leukemia 3 (MLL3) and its closest homolog, MLL4 (aka KMT2D), belong to two homologous transcriptional coactivator complexes, named MLL3 and MLL4 complexes, respectively. MLL3 plays crucial roles in multiple metabolic processes. However, the physiological roles of MLL4 in metabolism and the relationship between MLL3 and MLL4 in metabolic gene regulation are unclear. To address these issues, we analyzed the phenotypes of newly generated MLL4 mutant mice, along with MLL3 mutant and MLL3;MLL4 compound mutant mice. We also performed comparative genome-wide transcriptome analyses in livers of MLL3, MLL4, and MLL3;MLL4 mutant mice. These analyses revealed that MLL3 and MLL4 complexes are key epigenetic regulators of common metabolic processes and the hepatic circadian clock. Subsequent mechanistic analyses uncovered that MLL3/4 complexes function as pivotal coactivators of the circadian transcription factors (TFs), retinoid-related orphan receptor (ROR)-α and -γ, in the hepatic circadian clock. Consistent with disturbed hepatic clock gene expression in MLL4 mutant mice, we found that rhythmic fluctuation of hepatic and serum bile acid (BA) levels over the circadian cycle is abolished in MLL4 mutant mice. Our analyses also demonstrate that MLL4 primarily impinges on hepatic BA production among several regulatory pathways to control BA homeostasis. Together, our results provide strong in vivo support for important roles of both MLL3 and MLL4 in similar metabolic pathways. Conclusion: Both MLL3 and MLL4 complexes act as major epigenetic regulators of diverse metabolic processes (including circadian control of bile acid homeostasis) and as critical transcriptional coactivators of the circadian TFs, RORs. (Hepatology 2015;61:1012–1023)

While microRNAs have emerged as an important component of gene regulatory networks, it remains unclear how microRNAs collaborate with transcription factors in the gene networks that determines neuronal cell fate. Here we show that in the developing spinal cord, the expression of miR-218 is directly upregulated by the Isl1–Lhx3 complex, which drives motor neuron fate. Inhibition of miR-218 suppresses the generation of motor neurons in both chick neural tube and mouse embryonic stem cells, suggesting that miR-218 plays a crucial role in motor neuron differentiation. Results from unbiased RISC-trap screens, in vivo reporter assays and overexpression studies indicated that miR-218 directly represses transcripts that promote developmental programs for interneurons. In addition, we found that miR-218 activity is required for Isl1–Lhx3 to effectively induce motor neurons and suppress interneuron fates. Together our results reveal an essential role of miR-218 as a downstream effector of the Isl1–Lhx3 complex in establishing motor neuron identity. microRNAs have emerged as important components of numerous gene regulatory networks. Here the authors demonstrate that miR-218 is an essential component of the gene regulatory network that controls motor neuron fate specification in the developing spinal cord.

The establishment of correct neurotransmitter characteristics is an essential step of neuronal fate specification in CNS development. However, very little is known about how a battery of genes involved in the determination of a specific type of chemical-driven neurotransmission is coordinately regulated during vertebrate development. Here, we investigated the gene regulatory networks that specify the cholinergic neuronal fates in the spinal cord and forebrain, specifically, spinal motor neurons (MNs) and forebrain cholinergic neurons (FCNs). Conditional inactivation of Isl1, a LIM homeodomain factor expressed in both differentiating MNs and FCNs, led to a drastic loss of cholinergic neurons in the developing spinal cord and forebrain. We found that Isl1 forms two related, but distinct types of complexes, the Isl1-Lhx3-hexamer in MNs and the Isl1-Lhx8-hexamer in FCNs. Interestingly, our genome-wide ChIP-seq analysis revealed that the Isl1-Lhx3-hexamer binds to a suite of cholinergic pathway genes encoding the core constituents of the cholinergic neurotransmission system, such as acetylcholine synthesizing enzymes and transporters. Consistently, the Isl1-Lhx3-hexamer directly coordinated upregulation of cholinergic pathways genes in embryonic spinal cord. Similarly, in the developing forebrain, the Isl1-Lhx8-hexamer was recruited to the cholinergic gene battery and promoted cholinergic gene expression. Furthermore, the expression of the Isl1-Lhx8-complex enabled the acquisition of cholinergic fate in embryonic stem cell-derived neurons. Together, our studies show a shared molecular mechanism that determines the cholinergic neuronal fate in the spinal cord and forebrain, and uncover an important gene regulatory mechanism that directs a specific neurotransmitter identity in vertebrate CNS development., Author Summary Neurons utilize various chemicals to transmit signals to a target cell. Distinct types of neurons in the spinal cord and forebrain, collectively termed cholinergic neurons, utilize the same chemical, acetylcholine, for signal transmission. These neurons play critical roles in controlling locomotion and cognition. In this study, we have found that the Isl1 gene orchestrates the process to generate cholinergic neurons in the spinal cord and forebrain. Isl1 forms two different types of multi-protein complexes in the spinal cord and forebrain. Both complexes bind the same genomic regions in a group of genes critical for cholinergic signal transmission, and promote their simultaneous expression. These cholinergic genes include enzymes that synthesize acetylcholine and proteins required to package acetylcholine into vesicles. The Isl1-containing multi-protein complexes were able to trigger the generation of cholinergic neurons in embryonic stem cells and neural stem cells. Our study reveals crucial mechanisms to coordinate the expression of genes in the same biological pathway in different cell types. Furthermore, it suggests a new strategy to produce cholinergic neurons from stem cells.

The motor neuron (MN)–hexamer complex consisting of LIM homeobox 3, Islet-1, and nuclear LIM interactor is a key determinant of motor neuron specification and differentiation. To gain insights into the transcriptional network in motor neuron development, we performed a genome-wide ChIP-sequencing analysis and found that the MN–hexamer directly regulates a wide array of motor neuron genes by binding to the HxRE (hexamer response element) shared among the target genes. Interestingly, STAT3-binding motif is highly enriched in the MN–hexamer–bound peaks in addition to the HxRE. We also found that a transcriptionally active form of STAT3 is expressed in embryonic motor neurons and that STAT3 associates with the MN–hexamer, enhancing the transcriptional activity of the MN–hexamer in an upstream signal-dependent manner. Correspondingly, STAT3 was needed for motor neuron differentiation in the developing spinal cord. Together, our studies uncover crucial gene regulatory mechanisms that couple MN–hexamer and STAT-activating extracellular signals to promote motor neuron differentiation in vertebrate spinal cord.

The DIURNAL project ( http://diurnal.cgrb.oregonstate.edu/ ) provides a graphical interface for mining and viewing diurnal and circadian microarray data for Arabidopsis thaliana, poplar, and rice. The database is searchable and provides access to several user-friendly Web-based data-mining tools with easy-to-understand output. The associated tools include HAYSTACK ( http://haystack.cgrb.oregonstate.edu/ ) and ELEMENT ( http://element.cgrb.oregonstate.edu/ ). HAYSTACK is a model-based pattern-matching algorithm for identifying genes that are coexpressed and potentially coregulated. HAYSTACK can be used to analyze virtually any large-scale microarray data set and provides an alternative method for clustering microarray data from any experimental system by grouping together genes whose expression patterns match the same or similar user-defined patterns. ELEMENT is a Web-based program for identifying potential cis-regulatory elements in the promoters of coregulated genes in Arabidopsis, poplar, and rice. Together, DIURNAL, HAYSTACK, and ELEMENT can be used to facilitate cross-species comparisons among the plant species supported and to accelerate functional genomics efforts in the laboratory.

Despite representing only a small fraction of hippocampal granule cells, adult-generated newborn granule cells have been implicated in learning and memory (Aimone et al., 2011). Newborn granule cells undergo functional maturation and circuit integration over a period of weeks. However, it is difficult to assess the accompanying gene expression profiles in vivo with high spatial and temporal resolution using traditional methods. Here we used a novel method [“thiouracil (TU) tagging”] to map the profiles of nascent mRNAs in mouse immature newborn granule cells compared with mature granule cells. We targeted a nonmammalian uracil salvage enzyme, uracil phosphoribosyltransferase, to newborn neurons and mature granule cells using retroviral and lentiviral constructs, respectively. Subsequent injection of 4-TU tagged nascent RNAs for analysis by RNA sequencing. Several hundred genes were significantly enhanced in the retroviral dataset compared with the lentiviral dataset. We compared a selection of the enriched genes with steady-state levels of mRNAs using quantitative PCR. Ontology analysis revealed distinct patterns of nascent mRNA expression, with newly generated immature neurons showing enhanced expression for genes involved in synaptic function, and neural differentiation and development, as well as genes not previously associated with granule cell maturation. Surprisingly, the nascent mRNAs enriched in mature cells were related to energy homeostasis and metabolism, presumably indicative of the increased energy demands of synaptic transmission and their complex dendritic architecture. The high spatial and temporal resolution of our modified TU-tagging method provides a foundation for comparison with steady-state RNA analyses by traditional transcriptomic approaches in defining the functional roles of newborn neurons.

Our laboratory has reported data suggesting a role for T4 phage gene 32 single-stranded DNA-binding protein in organizing a complex of deoxyribonucleotide-synthesizing enzymes at the replication fork. In this article we examined the effects of gene 32 ablation on the association of these enzymes with DNA-protein complexes. These experiments showed several deoxyribonucleotide-synthesizing enzymes to be present in DNA-protein complexes, with some of these associations being dependent on gene 32 protein. To further understand the role of gp32, we created amber mutations at codons 24 and 204 of gene 32, which encodes a 301-residue protein. We used the newly created mutants along with several experimental approaches--DNA-cellulose chromatography, immunoprecipitation, optical biosensor analysis and glutathione-S-transferase pulldowns--to identify relevant protein-protein and protein-DNA interactions. These experiments identified several proteins whose interactions with DNA depend on the presence of intact gp32, notably thymidylate synthase, dihydrofolate (DHF) reductase, ribonucleotide reductase (RNR) and Escherichia coli nucleoside diphosphate (NDP) kinase, and they also demonstrated direct associations between gp32 and RNR and NDP kinase, but not dCMP hydroxymethylase, deoxyribonucleoside monophosphate kinase, or DHF reductase. Taken together, the results support the hypothesis that the gene 32 protein helps to recruit enzymes of deoxyribonucleoside triphosphates synthesis to DNA replication sites.

In both prokaryotic and eukaryotic organisms, nucleoside diphosphate kinase is a multifunctional protein, with well defined functions in ribo- and deoxyribonucleoside triphosphate biosynthesis and more recently described functions in genetic and metabolic regulation, signal transduction, and DNA repair. This paper concerns two unusual properties of nucleoside diphosphate (NDP) kinase from Escherichia coli: 1) its ability to interact specifically with enzymes encoded by the virulent bacteriophage T4 and 2) its roles in regulating metabolism of the host cell. By means of optical biosensor analysis, fluorescence spectroscopy, immunoprecipitation, and glutathione S-transferase pull-down assays, we have shown that E. coli NDP kinase interacts directly with T4 thymidylate synthase, aerobic ribonucleotide reductase, dCTPase-dUTPase, gene 32 single-strand DNA-binding protein, and deoxycytidylate hydroxymethylase. The interactions with ribonucleotide reductase and with gp32 are enhanced by nucleoside triphosphates, suggesting that the integrity of the T4 dNTP synthetase complex in vivo is influenced by the composition of the nucleotide pool. The other investigations in this work stem from the unexpected finding that E. coli NDP kinase is dispensable for successful T4 phage infection, and they deal with two observations suggesting that the NDP kinase protein plays a genetic role in regulating metabolism of the host cell: 1) the elevation of CTP synthetase activity in an ndk mutant, in which the structural gene for NDP kinase is disrupted, and 2) the apparent ability of NDP kinase to suppress anaerobic growth in a pyruvate kinase-negative E. coli mutant. Our data indicate that the regulatory roles are metabolic, not genetic, in nature.

Much of the morbidity and mortality caused by Acute Myeloid Leukemia (AML) is secondary to failure of normal hematopoiesis. Several recent studies indicate that this is not due to "overcrowding" within the bone marrow, suggesting other mechanisms account for the loss of bone marrow regenerative function in the leukemic niche. We recently reported that AML patient leukemia cells produce exosomes – small, membrane-enclosed extracellular vesicles – that are enriched in microRNA (miRNA) content and alter gene expression, cytokine secretion, and homing upon entry into neighboring stromal and progenitor cells. In order to examine the consequences of exosome trafficking in vivo, we conducted a series of experiments in a xenograft model. Here, we show that introduction of either AML cells or isolated AML-derived exosomes into the marrow of immunodeficient mice causes suppression of marrow hematopoietic stem and progenitor cell (HSPC) clonogenicity and a dramatic increase in peripheral circulating HSPC that coincided with a marked decrease in expression of Scf and Cxcl12 in murine stromal cells. Hypothesizing that transferred miRNA mediates the changes observed in recipient cells, we evaluated the miRNA content of two AML cell lines and their exosomes using microarrays paired with qRT-PCR. These experiments revealed stark differences between the total miRNA produced by AML cells and the selective incorporation of miRNA in exosomes. Combining AML exosome-enriched miRNA identified in our screen with a review of miRNA previously identified as significant in AML biology, we selected a panel of highly enriched miRNA, including miR-155 (100-fold enriched in exosomes) for further study. In order to mechanistically identify the mRNA targets of these miRNA, we used the uniquely powerful RISC-Trap method (Cambronne et al, PNAS, 2012), beginning with miR-155. This assay revealed 131 likely targets of miR-155, a set which we compared with several other target detection methods previously attempted with this miRNA, finding both significant overlap and several promising new putative mRNA targets. Of these 131 RISC-Trap targets, 98 were predicted by one or more of six tested target prediction algorithms (miRWalk). We compared the miR-155 targets identified by RISC-Trap to predicted targets for several other miRNA in the panel selected from our microarray studies, identifying both miRNA potentially responsible for the expression changes seen in stroma as well as several potential targets common to multiple miRNA. Using bioinformatics analysis of direct targets of our miRNA of interest and their common interacting partners, we were able to identify a network of targets and cellular processes regulated by AML exosomes, including regulation of transcription (p53, SOX9, CEBPB) and apoptosis (ESR1, TRAF2, CHEK2). These networks may provide insight into the mechanisms by which AML impairs hematopoiesis, both directly, through effects on HSPC, and indirectly, through effects on stroma. Ongoing RISC-Trap studies of other miRNA provide a forward genetic screen to evaluate the paracrine targets of exosomal miRNA in a mechanistic way, in specific niche cell types. By identifying the wider networks in addition to the individual contributors of exosomal miRNA in the AML microenvironment, we hope to uncover new non-leukemia targeted “niche” therapies to relieve hematopoietic suppression in AML patients. Disclosures Cambronne: Clontech/Takara Bio: commercial licensing of the RISC-trap technology Patents & Royalties.

Organellar DNA sequences are widely used in evolutionary and population genetic studies, however, the conservative nature of chloroplast gene and genome evolution often limits phylogenetic resolution and statistical power. To gain maximal access to the historical record contained within chloroplast genomes, we have adapted multiplex sequencing-by-synthesis (MSBS) to simultaneously sequence multiple genomes using the Illumina Genome Analyzer. We PCR-amplified approximately 120 kb plastomes from eight species (seven Pinus, one Picea) in 35 reactions. Pooled products were ligated to modified adapters that included 3 bp indexing tags and samples were multiplexed at four genomes per lane. Tagged microreads were assembled by de novo and reference-guided assembly methods, using previously published Pinus plastomes as surrogate references. Assemblies for these eight genomes are estimated at 88-94% complete, with an average sequence depth of 55x to 186x. Mononucleotide repeats interrupt contig assembly with increasing repeat length, and we estimate that the limit for their assembly is 16 bp. Comparisons to 37 kb of Sanger sequence show a validated error rate of 0.056%, and conspicuous errors are evident from the assembly process. This efficient sequencing approach yields high-quality draft genomes and should have immediate applicability to genomes with comparable complexity.