Jeong Hyeon Choi

The ecoresponsive genome of Daphnia pulex

The Daphnia genome reveals a multitude of genes and shows adaptation through gene family expansions. We describe the draft genome of the microcrustacean Daphnia pulex, which is only 200 megabases and contains at least 30,907 genes. The high gene count is a consequence of an elevated rate of gene duplication resulting in tandem gene clusters. More than a third of Daphnia’s genes have no detectable homologs in any other available proteome, and the most amplified gene families are specific to the Daphnia lineage. The coexpansion of gene families interacting within metabolic pathways suggests that the maintenance of duplicated genes is not random, and the analysis of gene expression under different environmental conditions reveals that numerous paralogs acquire divergent expression patterns soon after duplication. Daphnia-specific genes, including many additional loci within sequenced regions that are otherwise devoid of annotations, are the most responsive genes to ecological challenges.

Gut pathology and responses to the microsporidium Nosema ceranae in the honey bee Apis mellifera.

The microsporidium Nosema ceranae is a newly prevalent parasite of the European honey bee (Apis mellifera). Although this parasite is presently spreading across the world into its novel host, the mechanisms by it which affects the bees and how bees respond are not well understood. We therefore performed an extensive characterization of the parasite effects at the molecular level by using genetic and biochemical tools. The transcriptome modifications at the midgut level were characterized seven days post-infection with tiling microarrays. Then we tested the bee midgut response to infection by measuring activity of antioxidant and detoxification enzymes (superoxide dismutases, glutathione peroxidases, glutathione reductase, and glutathione-S-transferase). At the gene-expression level, the bee midgut responded to N. ceranae infection by an increase in oxidative stress concurrent with the generation of antioxidant enzymes, defense and protective response specifically observed in the gut of mammals and insects. However, at the enzymatic level, the protective response was not confirmed, with only glutathione-S-transferase exhibiting a higher activity in infected bees. The oxidative stress was associated with a higher transcription of sugar transporter in the gut. Finally, a dramatic effect of the microsporidia infection was the inhibition of genes involved in the homeostasis and renewal of intestinal tissues (Wnt signaling pathway), a phenomenon that was confirmed at the histological level. This tissue degeneration and prevention of gut epithelium renewal may explain early bee death. In conclusion, our integrated approach not only gives new insights into the pathological effects of N. ceranae and the bee gut response, but also demonstrate that the honey bee gut is an interesting model system for studying host defense responses.

Linking DNA Methyltransferases to Epigenetic Marks and Nucleosome Structure Genome-wide in Human Tumor Cells

DNA methylation, mediated by the combined action of three DNA methyltransferases (DNMT1, DNMT3A, and DNMT3B), is essential for mammalian development and is a major contributor to cellular transformation. To elucidate how DNA methylation is targeted, we mapped the genome-wide localization of all DNMTs and methylation, and examined the relationships among these markers, histone modifications, and nucleosome structure in a pluripotent human tumor cell line in its undifferentiated and differentiated states. Our findings reveal a strong link between DNMTs and transcribed loci, and that DNA methylation is not a simple sum of DNMT localization patterns. A comparison of the epigenomes of normal and cancerous stem cells, and pluripotent and differentiated states shows that the presence of at least two DNMTs is strongly associated with loci targeted for DNA hypermethylation. Taken together, these results shed important light on the determinants of DNA methylation and how it may become disrupted in cancer cells.

Function and Evolution of DNA Methylation in Nasonia vitripennis

The parasitoid wasp Nasonia vitripennis is an emerging genetic model for functional analysis of DNA methylation. Here, we characterize genome-wide methylation at a base-pair resolution, and compare these results to gene expression across five developmental stages and to methylation patterns reported in other insects. An accurate assessment of DNA methylation across the genome is accomplished using bisulfite sequencing of adult females from a highly inbred line. One-third of genes show extensive methylation over the gene body, yet methylated DNA is not found in non-coding regions and rarely in transposons. Methylated genes occur in small clusters across the genome. Methylation demarcates exon-intron boundaries, with elevated levels over exons, primarily in the 5′ regions of genes. It is also elevated near the sites of translational initiation and termination, with reduced levels in 5′ and 3′ UTRs. Methylated genes have higher median expression levels and lower expression variation across development stages than non-methylated genes. There is no difference in frequency of differential splicing between methylated and non-methylated genes, and as yet no established role for methylation in regulating alternative splicing in Nasonia. Phylogenetic comparisons indicate that many genes maintain methylation status across long evolutionary time scales. Nasonia methylated genes are more likely to be conserved in insects, but even those that are not conserved show broader expression across development than comparable non-methylated genes. Finally, examination of duplicated genes shows that those paralogs that have lost methylation in the Nasonia lineage following gene duplication evolve more rapidly, show decreased median expression levels, and increased specialization in expression across development. Methylation of Nasonia genes signals constitutive transcription across developmental stages, whereas non-methylated genes show more dynamic developmental expression patterns. We speculate that loss of methylation may result in increased developmental specialization in evolution and acquisition of methylation may lead to broader constitutive expression.

RAPSearch: a fast protein similarity search tool for short reads.

BackgroundNext Generation Sequencing (NGS) is producing enormous corpuses of short DNA reads, affecting emerging fields like metagenomics. Protein similarity search--a key step to achieve annotation of protein-coding genes in these short reads, and identification of their biological functions--faces daunting challenges because of the very sizes of the short read datasets.ResultsWe developed a fast protein similarity search tool RAPSearch that utilizes a reduced amino acid alphabet and suffix array to detect seeds of flexible length. For short reads (translated in 6 frames) we tested, RAPSearch achieved ~20-90 times speedup as compared to BLASTX. RAPSearch missed only a small fraction (~1.3-3.2%) of BLASTX similarity hits, but it also discovered additional homologous proteins (~0.3-2.1%) that BLASTX missed. By contrast, BLAT, a tool that is even slightly faster than RAPSearch, had significant loss of sensitivity as compared to RAPSearch and BLAST.ConclusionsRAPSearch is implemented as open-source software and is accessible at http://omics.informatics.indiana.edu/mg/RAPSearch. It enables faster protein similarity search. The application of RAPSearch in metageomics has also been demonstrated.

DNMT1 is essential for mammary and cancer stem cell maintenance and tumorigenesis

Mammary stem/progenitor cells (MaSCs) maintain self-renewal of the mammary epithelium during puberty and pregnancy. DNA methylation provides a potential epigenetic mechanism for maintaining cellular memory during self-renewal. Although DNA methyltransferases (DNMTs) are dispensable for embryonic stem cell maintenance, their role in maintaining MaSCs and cancer stem cells (CSCs) in constantly replenishing mammary epithelium is unclear. Here we show that DNMT1 is indispensable for MaSC maintenance. Furthermore, we find that DNMT1 expression is elevated in mammary tumors, and mammary gland-specific DNMT1 deletion protects mice from mammary tumorigenesis by limiting the CSC pool. Through genome-scale methylation studies, we identify ISL1 as a direct DNMT1 target, hypermethylated and downregulated in mammary tumors and CSCs. DNMT inhibition or ISL1 expression in breast cancer cells limits CSC population. Altogether, our studies uncover an essential role for DNMT1 in MaSC and CSC maintenance and identify DNMT1-ISL1 axis as a potential therapeutic target for breast cancer treatment.

Genome-wide DNA methylation analysis reveals novel epigenetic changes in chronic lymphocytic leukemia

We conducted a genome-wide DNA methylation analysis in CD19+ B-cells from chronic lymphocytic leukemia (CLL) patients and normal control samples using reduced representation bisulfite sequencing (RRBS). The methylation status of 1.8–2.3 million CpGs in the CLL genome was determined; about 45% of these CpGs were located in more than 23,000 CpG islands (CGIs). While global CpG methylation was similar between CLL and normal B-cells, 1764 gene promoters were identified as being differentially methylated in at least one CLL sample when compared with normal B-cell samples. Nineteen percent of the differentially methylated genes were involved in transcriptional regulation. Aberrant hypermethylation was found in all HOX gene clusters and a significant number of WNT signaling pathway genes. Hypomethylation occurred more frequently in the gene body including introns, exons, and 3′-UTRs in CLL. The NFATc1 P2 promoter and first intron was found to be hypomethylated and correlated with upregulation of both NFATc1 RNA and protein expression levels in CLL suggesting that an epigenetic mechanism is involved in the constitutive activation of NFAT activity in CLL cells. This comprehensive DNA methylation analysis will further our understanding of the epigenetic contribution to cellular dysfunction in CLL.

A garter snake transcriptome: pyrosequencing, de novo assembly, and sex-specific differences

BackgroundThe reptiles, characterized by both diversity and unique evolutionary adaptations, provide a comprehensive system for comparative studies of metabolism, physiology, and development. However, molecular resources for ectothermic reptiles are severely limited, hampering our ability to study the genetic basis for many evolutionarily important traits such as metabolic plasticity, extreme longevity, limblessness, venom, and freeze tolerance. Here we use massively parallel sequencing (454 GS-FLX Titanium) to generate a transcriptome of the western terrestrial garter snake (Thamnophis elegans) with two goals in mind. First, we develop a molecular resource for an ectothermic reptile; and second, we use these sex-specific transcriptomes to identify differences in the presence of expressed transcripts and potential genes of evolutionary interest.ResultsUsing sex-specific pools of RNA (one pool for females, one pool for males) representing 7 tissue types and 35 diverse individuals, we produced 1.24 million sequence reads, which averaged 366 bp in length after cleaning. Assembly of the cleaned reads from both sexes with NEWBLER and MIRA resulted in 96,379 contigs containing 87% of the cleaned reads. Over 34% of these contigs and 13% of the singletons were annotated based on homology to previously identified proteins. From these homology assignments, additional clustering, and ORF predictions, we estimate that this transcriptome contains ~13,000 unique genes that were previously identified in other species and over 66,000 transcripts from unidentified protein-coding genes. Furthermore, we use a graph-clustering method to identify contigs linked by NEWBLER-split reads that represent divergent alleles, gene duplications, and alternatively spliced transcripts. Beyond gene identification, we identified 95,295 SNPs and 31,651 INDELs. From these sex-specific transcriptomes, we identified 190 genes that were only present in the mRNA sequenced from one of the sexes (84 female-specific, 106 male-specific), and many highly variable genes of evolutionary interest.ConclusionsThis is the first large-scale, multi-organ transcriptome for an ectothermic reptile. This resource provides the most comprehensive set of EST sequences available for an individual ectothermic reptile species, increasing the number of snake ESTs 50-fold. We have identified genes that appear to be under evolutionary selection and those that are sex-specific. This resource will assist studies on gene expression and comparative genomics, and will facilitate the study of evolutionarily important traits at the molecular level.

De novo identification of LTR retrotransposons in eukaryotic genomes

BackgroundLTR retrotransposons are a class of mobile genetic elements containing two similar long terminal repeats (LTRs). Currently, LTR retrotransposons are annotated in eukaryotic genomes mainly through the conventional homology searching approach. Hence, it is limited to annotating known elements.ResultsIn this paper, we report a de novo computational method that can identify new LTR retrotransposons without relying on a library of known elements. Specifically, our method identifies intact LTR retrotransposons by using an approximate string matching technique and protein domain analysis. In addition, it identifies partially deleted or solo LTRs using profile Hidden Markov Models (pHMMs). As a result, this method can de novo identify all types of LTR retrotransposons. We tested this method on the two pairs of eukaryotic genomes, C. elegans vs. C. briggsae and D. melanogaster vs. D. pseudoobscura. LTR retrotransposons in C. elegans and D. melanogaster have been intensively studied using conventional annotation methods. Comparing with previous work, we identified new intact LTR retroelements and new putative families, which may imply that there may still be new retroelements that are left to be discovered even in well-studied organisms. To assess the sensitivity and accuracy of our method, we compared our results with a previously published method, LTR_STRUC, which predominantly identifies full-length LTR retrotransposons. In summary, both methods identified comparable number of intact LTR retroelements. But our method can identify nearly all known elements in C. elegans, while LTR_STRUCT missed about 1/3 of them. Our method also identified more known LTR retroelements than LTR_STRUCT in the D. melanogaster genome. We also identified some LTR retroelements in the other two genomes, C. briggsae and D. pseudoobscura, which have not been completely finished. In contrast, the conventional method failed to identify those elements. Finally, the phylogenetic and chromosomal distributions of the identified elements are discussed.ConclusionWe report a novel method for de novo identification of LTR retrotransposons in eukaryotic genomes with favorable performance over the existing methods.

Characteristics of the genome of Arsenophonus nasoniae, son-killer bacterium of the wasp Nasonia

We report the properties of a draft genome sequence of the bacterium Arsenophonus nasoniae, son‐killer bacterium of Nasonia vitripennis. The genome sequence data from this study are the first for a male‐killing bacterium, and represent a microorganism that is unusual compared with other sequenced symbionts, in having routine vertical and horizontal transmission, two alternating hosts, and being culturable on cell‐free media. The resulting sequence totals c. 3.5 Mbp and is annotated to contain 3332 predicted open reading frames (ORFs). Therefore, Arsenophonus represents a relatively large genome for an insect symbiont. The annotated ORF set suggests that the microbe is capable of a broad array of metabolic functions, well beyond those found for reproductive parasite genomes sequenced to date and more akin to horizontally transmitted and secondary symbionts. We also find evidence of genetic transfer from Wolbachia symbionts, and phage exchange with other gammaproteobacterial symbionts. These findings reflect the complex biology of a bacterium that is able to live, invade and survive multiple host environments while resisting immune responses.