A Homozygous Missense Variant in INPP5E Associated with Joubert Syndrome and Related Disorders

Joubert syndrome and related disorders (JSRD; ORPHA 140874) is a complex set of neurodevelopmental disorders with multiple organ involvement. JSRD is a type of ciliopathy which is caused by the presence of defective primary cilia in an individual. JSRD is commonly inherited in an autosomal recessive pattern, and more than 23 genes are known to be associated with JSRD. We report a novel homozygous mutation identified in the INPP5E gene, c.1303C>T, which leads to a change of an amino acid from arginine to tryptophan at residue 435 in the protein chain. In silico analysis indicates that p.Arg435Trp substitution affects the functionality of the protein product of the gene. Our result adds to the growing body of evidences that underlines the clinical utility of next-generation sequencing in the diagnosis of a genetic disorder when clinical features are inconclusive.

Genome Wide Methylome Alterations in Lung Cancer

Aberrant cytosine 5-methylation underlies many deregulated elements of cancer. Among paired non-small cell lung cancers (NSCLC), we sought to profile DNA 5-methyl-cytosine features which may underlie genome-wide deregulation. In one of the more dense interrogations of the methylome, we sampled 1.2 million CpG sites from twenty-four NSCLC tumor (T)-non-tumor (NT) pairs using a methylation-sensitive restriction enzyme- based HELP-microarray assay. We found 225,350 differentially methylated (DM) sites in adenocarcinomas versus adjacent non-tumor tissue that vary in frequency across genomic compartment, particularly notable in gene bodies (GB; p<2.2E-16). Further, when DM was coupled to differential transcriptome (DE) in the same samples, 37,056 differential loci in adenocarcinoma emerged. Approximately 90% of the DM-DE relationships were non-canonical; for example, promoter DM associated with DE in the same direction. Of the canonical changes noted, promoter (PR) DM loci with reciprocal changes in expression in adenocarcinomas included HBEGF, AGER, PTPRM, DPT, CST1, MELK; DM GB loci with concordant changes in expression included FOXM1, FERMT1, SLC7A5, and FAP genes. IPA analyses showed adenocarcinoma-specific promoter DMxDE overlay identified familiar lung cancer nodes [tP53, Akt] as well as less familiar nodes [HBEGF, NQO1, GRK5, VWF, HPGD, CDH5, CTNNAL1, PTPN13, DACH1, SMAD6, LAMA3, AR]. The unique findings from this study include the discovery of numerous candidate The unique findings from this study include the discovery of numerous candidate methylation sites in both PR and GB regions not previously identified in NSCLC, and many non-canonical relationships to gene expression. These DNA methylation features could potentially be developed as risk or diagnostic biomarkers, or as candidate targets for newer methylation locus-targeted preventive or therapeutic agents.

Comparison of newly diagnosed and relapsed patients with acute promyelocytic leukemia treated with arsenic trioxide: insight into mechanisms of resistance

There is limited data on the clinical, cellular and molecular changes in relapsed acute promyeloytic leukemia (RAPL) in comparison with newly diagnosed cases (NAPL). We undertook a prospective study to compare NAPL and RAPL patients treated with arsenic trioxide (ATO) based regimens. 98 NAPL and 28 RAPL were enrolled in this study. RAPL patients had a significantly lower WBC count and higher platelet count at diagnosis. IC bleeds was significantly lower in RAPL cases (P=0.022). The ability of malignant promyelocytes to concentrate ATO intracellularly and their in-vitro IC50 to ATO was not significantly different between the two groups. Targeted NGS revealed PML B2 domain mutations in 4 (15.38%) of the RAPL subset and none were associated with secondary resistance to ATO. A microarray GEP revealed 1744 genes were 2 fold and above differentially expressed between the two groups. The most prominent differentially regulated pathways were cell adhesion (n=92), cell survival (n=50), immune regulation (n=74) and stem cell regulation (n=51). Consistent with the GEP data, immunophenotyping revealed significantly increased CD34 expression (P=0.001) in RAPL cases and there was in-vitro evidence of significant microenvironment mediated innate resistance (EM-DR) to ATO. Resistance and relapse following treatment with ATO is probably multi-factorial, mutations in PML B2 domain while seen only in RAPL may not be the major clinically relevant cause of subsequent relapses. In RAPL additional factors such as expansion of the leukemia initiating compartment along with EM-DR may contribute significantly to relapse following treatment with ATO based regimens.

Lung cancer transcriptomes refined with laser capture microdissection

We evaluated the importance of tumor cell selection for generating gene signatures in non-small cell lung cancer. Tumor and nontumor tissue from macroscopically dissected (Macro) surgical specimens (31 pairs from 32 subjects) was homogenized, extracted, amplified, and hybridized to microarrays. Adjacent scout sections were histologically mapped; sets of approximately 1000 tumor cells and nontumor cells (alveolar or bronchial) were procured by laser capture microdissection (LCM). Within histological strata, LCM and Macro specimens exhibited approximately 67% to 80% nonoverlap in differentially expressed (DE) genes. In a representative subset, LCM uniquely identified 300 DE genes in tumor versus nontumor specimens, largely attributable to cell selection; 382 DE genes were common to Macro, Macro with preamplification, and LCM platforms. RT-qPCR validation in a 33-gene subset was confirmatory (ρ = 0.789 to 0.964, P = 0.0013 to 0.0028). Pathway analysis of LCM data suggested alterations in known cancer pathways (cell growth, death, movement, cycle, and signaling components), among others (eg, immune, inflammatory). A unique nine-gene LCM signature had higher tumor-nontumor discriminatory accuracy (100%) than the corresponding Macro signature (87%). Comparison with Cancer Genome Atlas data sets (based on homogenized Macro tissue) revealed both substantial overlap and important differences from LCM specimen results. Thus, cell selection via LCM enhances expression profiling precision, and confirms both known and under-appreciated lung cancer genes and pathways.

Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation

Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.

Cryptococcus neoformans var. grubii is a major human pathogen responsible for deadly meningoencephalitis in immunocompromised patients. Here, we report the sequencing and annotation of its genome. Evidence for extensive intron splicing, antisense transcription, non-coding RNAs, and alternative polyadenylation indicates the potential for highly intricate regulation of gene expression in this opportunistic pathogen. In addition, detailed molecular, genetic, and genomic studies were performed to characterize structural features of the genome, including centromeres and origins of replication. Finally, the phenotypic and genome re-sequencing analysis of a collection of isolates of the reference H99 strain resulting from laboratory passage revealed that microevolutionary processes during in vitro culturing of pathogenic fungi can impact virulence.

An integrated pipeline for next generation sequencing and annotation of the complete mitochondrial genome of the giant intestinal fluke, Fasciolopsis buski (Lankester, 1857) Looss, 1899

Helminths include both parasitic nematodes (roundworms) and platyhelminths (trematode and cestode flatworms) that are abundant, and are of clinical importance. The genetic characterization of parasitic flatworms using advanced molecular tools is central to the diagnosis and control of infections. Although the nuclear genome houses suitable genetic markers (e.g., in ribosomal (r) DNA) for species identification and molecular characterization, the mitochondrial (mt) genome consistently provides a rich source of novel markers for informative systematics and epidemiological studies. In the last decade, there have been some important advances in mtDNA genomics of helminths, especially lung flukes, liver flukes and intestinal flukes. Fasciolopsis buski, often called the giant intestinal fluke, is one of the largest digenean trematodes infecting humans and found primarily in Asia, in particular the Indian subcontinent. Next-generation sequencing (NGS) technologies now provide opportunities for high throughput sequencing, assembly and annotation within a short span of time. Herein, we describe a high-throughput sequencing and bioinformatics pipeline for mt genomics for F. buski that emphasizes the utility of short read NGS platforms such as Ion Torrent and Illumina in successfully sequencing and assembling the mt genome using innovative approaches for PCR primer design as well as assembly. We took advantage of our NGS whole genome sequence data (unpublished so far) for F. buski and its comparison with available data for the Fasciola hepatica mtDNA as the reference genome for design of precise and specific primers for amplification of mt genome sequences from F. buski. A long-range PCR was carried out to create an NGS library enriched in mt DNA sequences. Two different NGS platforms were employed for complete sequencing, assembly and annotation of the F. buski mt genome. The complete mt genome sequences of the intestinal fluke comprise 14,118 bp and is thus the shortest trematode mitochondrial genome sequenced to date. The noncoding control regions are separated into two parts by the tRNA-Gly gene and don’t contain either tandem repeats or secondary structures, which are typical for trematode control regions. The gene content and arrangement are identical to that of F. hepatica. The F. buski mtDNA genome has a close resemblance with F. hepatica and has a similar gene order tallying with that of other trematodes. The mtDNA for the intestinal fluke is reported herein for the first time by our group that would help investigate Fasciolidae taxonomy and systematics with the aid of mtDNA NGS data. More so, it would serve as a resource for comparative mitochondrial genomics and systematic studies of trematode parasites.

Next generation sequencing and de novo transcriptome analysis of Costus pictus D. Don, a non-model plant with potent anti-diabetic properties

Background

Phyto-remedies for diabetic control are popular among patients with Type II Diabetes mellitus (DM), in addition to other diabetic control measures. A number of plant species are known to possess diabetic control properties. Costus pictus D. Don is popularly known as “Insulin Plant” in Southern India whose leaves have been reported to increase insulin pools in blood plasma. Next Generation Sequencing is employed as a powerful tool for identifying molecular signatures in the transcriptome related to physiological functions of plant tissues. We sequenced the leaf transcriptome of C. pictus using Illumina reversible dye terminator sequencing technology and used combination of bioinformatics tools for identifying transcripts related to anti-diabetic properties of C. pictus.

Results

A total of 55,006 transcripts were identified, of which 69.15% transcripts could be annotated. We identified transcripts related to pathways of bixin biosynthesis and geraniol and geranial biosynthesis as major transcripts from the class of isoprenoid secondary metabolites and validated the presence of putative norbixin methyltransferase, a precursor of Bixin. The transcripts encoding these terpenoids are known to be Peroxisome Proliferator-Activated Receptor (PPAR) agonists and anti-glycation agents. Sequential extraction and High Performance Liquid Chromatography (HPLC) confirmed the presence of bixin in C. pictus methanolic extracts. Another significant transcript identified in relation to anti-diabetic, anti-obesity and immuno-modulation is of Abscisic Acid biosynthetic pathway. We also report many other transcripts for the biosynthesis of antitumor, anti-oxidant and antimicrobial metabolites of C. pictus leaves.

Conclusion

Solid molecular signatures (transcripts related to bixin, abscisic acid, and geranial and geraniol biosynthesis) for the anti-diabetic properties of C. pictus leaves and vital clues related to the other phytochemical functions like antitumor, anti-oxidant, immuno-modulatory, anti-microbial and anti-malarial properties through the secondary metabolite pathway annotations are reported. The data provided will be of immense help to researchers working in the treatment of DM using herbal therapies.

Partial promoter substitutions generating transcriptional sentinels of diverse signaling pathways in embryonic stem cells and mice

Extracellular signals in development, physiology, homeostasis and disease often act by regulating transcription. Herein we describe a general method and specific resources for determining where and when such signaling occurs in live animals and for systematically comparing the timing and extent of different signals in different cellular contexts. We used recombinase-mediated cassette exchange (RMCE) to test the effect of successively deleting conserved genomic regions of the ubiquitously active Rosa26 promoter and substituting the deleted regions for regulatory sequences that respond to diverse extracellular signals. We thereby created an allelic series of embryonic stem cells and mice, each containing a signal-responsive sentinel with different fluorescent reporters that respond with sensitivity and specificity to retinoic acids, bone morphogenic proteins, activin A, Wnts or Notch, and that can be adapted to any pathway that acts via DNA elements.

Abstract A38: Whole methylome explorations of paired lung tumor and non-tumor clinical samples

Background: Hypermethylation in specific candidate gene promoters has been found during progressive lung carcinogenesis. To explore common methylation events on a genome-wide scale in lung cancer, we analyzed the methylation profiles of paired NSCLC tumor and far adjacent non-tumor samples using the HELP-microarray assay, which yields information on 1.2 million fragments throughout the genome.

Methods: The HELP (HpaII tiny fragment Enriched by Ligation mediated PCR) assay is based on the generation of restriction enzyme libraries generated by methylation sensitive (HpaII) and methylation insensitive (MspI) isoschizomers, in its second generation as a microarray platform. The assay lends itself to low starting amounts of DNA (3 ug) and robust assessment of methylation status by comparing ratios of HpaII- generated-fragments to MspI- generated fragments co-hybridized to a Nimbelgen custom high-density microarrays. The CCGG sites were weighted if neighboring CCGGs were methylated in same direction. Here, 24 pairs of tumor and adjacent non-tumor samples were analyzed using the HELP assay.

Summary of results: At p = 5E-6, we identified 26,138 differentially methylated fragments (corresponding to 2 CpG sites each) in tumor versus non-tumor. The overall trend was consistent with genome-wide hypomethylation and locus specific hypermethylation (localized to CG-island containing promoters). We could identify both known and novel regions of the genome as well as specific gene-promoters that are hypermethylated in tumor versus non-tumor.

Region # loci # signif loci T Hypomethyl T Hypermethyl

Promoter 151,568 576 69% 31%

Gene Body 551,628 9,817 62% 38%

Intergenic 540,473 15,745 97% 3%

Conclusion: An interrogation of methylation status of 1.2 million loci throughout the genome in paired lung tumor/non-tumor specimens reveals many more differential methylation events in gene bodies and intergenic regions than in promoters. That said, many previously unreported differentially-methylated gene promoters were identified. We are able to discover individual methylation events common across different clinical specimens. Based on a set of priors, we have narrowed down promoter-specific hypermethylation events for further validation using tagged Bisulfite Genomic Sequencing. We are also working on the integration of methylome date with other genome-wide epigenetic and expression data from the same clinical samples.

[Funding source: NCI 1RC1 CA145422-01; 1K24-CA139054-01]

Abstract 4820: Whole methylome explorations of paired tumor and non-tumor clinical samples using HELP

Introduction: Gene silencing by promoter methylation is common in carcinogenesis. In lung cancer, hypermethylation of specific gene promoters has been found to occur during the earliest histologic stages of the cancer. Thus far surveys of methylation in lung tumors have focused on specific gene panels (5 – 10 genes). A few large scale studies have reported selection of longer gene lists or focused on larger sets of CpG islands.

Methods: We have adapted an approach to survey genome-wide methylation events in tumor and adjacent non-tumor tissue from clinical lung resection samples using the HELP (HpaII tiny fragment Enrichment by Ligation mediated PCR) assay. This assay, based on the generation of restriction enzyme libraries generated by methylation sensitive (HpaII) and methylation insensitive (MspI) isoschizomers yields information on 1.2 million HpaII fragments throughout the genome. Each HpaII fragment corresponds to two CpG sites. The assay lends itself to low starting amounts of DNA (3 ug) and robust assessment of methylation status by comparing ratios of HpaII- generated- fragments to MspI- generated fragments co-hybridized to a Nimbelgen custom high-density microarray.

Summary of results: In a pilot set of three paired tumor and non-tumor tissue samples, we identified 1504 fragments (corresponding to two CpG sites each) that were hypermethylated in tumor versus non-tumor in all three pairs. Of these, 165 fragments were located within gene promoters. In the opposite direction, 2907 fragments were found to be hypermethylated in non-tumor versus tumor. Of these, 81 were located within gene promoters. We could identify both known and novel regions of the genome as well as specific gene-promoters that are hypermethylated in tumor versus non-tumor. We will discuss selected methylation events in further detail.

Conclusion: The HELP assay is a powerful approach to interrogate methylation status of individual HpaII fragments across the whole genome, offering a large sample size of 1.2 million HpaII fragments. We are able to discover individual methylation events common across different clinical specimens and are currently expanding our sample set, and planning correlations with transcriptome-wide studies.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4820. doi:10.1158/1538-7445.AM2011-4820

Identification and functional characterization of cis-regulatory elements in the apicomplexan parasite Toxoplasma gondii

Mining of genomic sequence data of the apicomplexan parasite Toxoplasma gondii identifies putative cis-regulatory elements using a de novo approach.

Background

Toxoplasma gondii is a member of the phylum Apicomplexa, which consists entirely of parasitic organisms that cause several diseases of veterinary and human importance. Fundamental mechanisms of gene regulation in this group of protistan parasites remain largely uncharacterized. Owing to their medical and veterinary importance, genome sequences are available for several apicomplexan parasites. Their genome sequences reveal an apparent paucity of known transcription factors and the absence of canonical cis-regulatory elements. We have approached the question of gene regulation from a sequence perspective by mining the genomic sequence data to identify putative cis-regulatory elements using a de novo approach.

Results

We have identified putative cis-regulatory elements present upstream of functionally related groups of genes and subsequently characterized the function of some of these conserved elements using reporter assays in the parasite. We show a sequence-specific role in gene-expression for seven out of eight identified elements.

Conclusions

This work demonstrates the power of pure sequence analysis in the absence of expression data or a priori knowledge of regulatory elements in eukaryotic organisms with compact genomes.

Identification of putative cis-regulatory elements in Cryptosporidium parvum by de novo pattern finding

Background

Cryptosporidium parvum is a unicellular eukaryote in the phylum Apicomplexa. It is an obligate intracellular parasite that causes diarrhea and is a significant AIDS-related pathogen. Cryptosporidium parvum is not amenable to long-term laboratory cultivation or classical molecular genetic analysis. The parasite exhibits a complex life cycle, a broad host range, and fundamental mechanisms of gene regulation remain unknown. We have used data from the recently sequenced genome of this organism to uncover clues about gene regulation in C. parvum. We have applied two pattern finding algorithms MEME and AlignACE to identify conserved, over-represented motifs in the 5' upstream regions of genes in C. parvum. To support our findings, we have established comparative real-time -PCR expression profiles for the groups of genes examined computationally.

Results

We find that groups of genes that share a function or belong to a common pathway share upstream motifs. Different motifs are conserved upstream of different groups of genes. Comparative real-time PCR studies show co-expression of genes within each group (in sub-sets) during the life cycle of the parasite, suggesting co-regulation of these genes may be driven by the use of conserved upstream motifs.

Conclusion

This is one of the first attempts to characterize cis-regulatory elements in the absence of any previously characterized elements and with very limited expression data (seven genes only). Using de novo pattern finding algorithms, we have identified specific DNA motifs that are conserved upstream of genes belonging to the same metabolic pathway or gene family. We have demonstrated the co-expression of these genes (often in subsets) using comparative real-time-PCR experiments thus establishing evidence for these conserved motifs as putative cis-regulatory elements. Given the lack of prior information concerning expression patterns and organization of promoters in C. parvum we present one of the first investigations of gene regulation in this important human pathogen.

Pair of unusual GCN5 histone acetyltransferases and ADA2 homologues in the protozoan parasite Toxoplasma gondii

GCN5 is a histone acetyltransferase (HAT) essential for development in mammals and critical to stress responses in yeast. The protozoan parasite Toxoplasma gondii is a serious opportunistic pathogen. The study of epigenetics and gene expression in this ancient eukaryote has pharmacological relevance and may facilitate the understanding of these processes in higher eukaryotes. Here we show that the disruption of T. gondii GCN5 yields viable parasites, which were subsequently employed in a proteomics study to identify gene products affected by its loss. Promoter analysis of these TgGCN5-dependent genes, which were mostly parasite specific, reveals a conserved T-rich element. The loss of TgGCN5 does not attenuate virulence in an in vivo mouse model. We also discovered that T. gondii is the only invertebrate reported to date possessing a second GCN5 (TgGCN5-B). TgGCN5-B harbors a strikingly divergent N-terminal domain required for nuclear localization. Despite high homology between the HAT domains, the two TgGCN5s exhibit differing substrate specificities. In contrast to TgGCN5-A, which exclusively targets lysine 18 of H3, TgGCN5-B acetylates multiple lysines in the H3 tail. We also identify two ADA2 homologues that interact differently with the TgGCN5s. TgGCN5-B has the potential to compensate for TgGCN5-A, which probably arose from a gene duplication unique to T. gondii. Our work reveals an unexpected complexity in the GCN5 machinery of this primitive eukaryote.

Phylogenomic evidence supports past endosymbiosis, intracellular and horizontal gene transfer in Cryptosporidium parvum

An analysis of Cryptosporidium parvum genes of likely endosymbiont or prokaryotic origin supports the hypothesis that C. arvum evolved from a plastid-containing lineage.

Background

The apicomplexan parasite Cryptosporidium parvum is an emerging pathogen capable of causing illness in humans and other animals and death in immunocompromised individuals. No effective treatment is available and the genome sequence has recently been completed. This parasite differs from other apicomplexans in its lack of a plastid organelle, the apicoplast. Gene transfer, either intracellular from an endosymbiont/donor organelle or horizontal from another organism, can provide evidence of a previous endosymbiotic relationship and/or alter the genetic repertoire of the host organism. Given the importance of gene transfers in eukaryotic evolution and the potential implications for chemotherapy, it is important to identify the complement of transferred genes in Cryptosporidium.

Results

We have identified 31 genes of likely plastid/endosymbiont (n = 7) or prokaryotic (n = 24) origin using a phylogenomic approach. The findings support the hypothesis that Cryptosporidium evolved from a plastid-containing lineage and subsequently lost its apicoplast during evolution. Expression analyses of candidate genes of algal and eubacterial origin show that these genes are expressed and developmentally regulated during the life cycle of C. parvum.

Conclusions

Cryptosporidium is the recipient of a large number of transferred genes, many of which are not shared by other apicomplexan parasites. Genes transferred from distant phylogenetic sources, such as eubacteria, may be potential targets for therapeutic drugs owing to their phylogenetic distance or the lack of homologs in the host. The successful integration and expression of the transferred genes in this genome has changed the genetic and metabolic repertoire of the parasite.

A first glimpse into the pattern and scale of gene transfer in Apicomplexa

Reports of plant-like and bacterial-like genes for a number of parasitic organisms, most notably those within the Apicomplexa and Kinetoplastida, have appeared in the literature over the last few years. Among the apicomplexan organisms, following discovery of the apicomplexan plastid (apicoplast), the discovery of plant-like genes was less surprising although the extent of transfer and the relationship of transferred genes to the apicoplast remained unclear. We used new genome sequence data to begin a systematic examination of the extent and origin of transferred genes in the Apicomplexa combined with a phylogenomic approach to detect potential gene transfers in four apicomplexan genomes. We have detected genes of algal nuclear, chloroplast (cyanobacterial) and proteobacterial origin. Plant-like genes were detected in species not currently harbouring a plastid (e.g. Cryptosporidium parvum) and putatively transferred genes were detected that appear to be unrelated to the function of the apicoplast. While the mechanism of acquisition for many of the identified genes is not certain, it appears that some were most likely acquired via intracellular gene transfer from an algal endosymbiont while others may have been acquired via horizontal gene transfer.