Detecting, quantifying, and discriminating the mechanism of mosaic chromosomal aneuploidies using MAD-seq.

Current approaches to detect and characterize mosaic chromosomal aneuploidy are limited by sensitivity, efficiency, cost, or the need to culture cells. We describe the mosaic aneuploidy detection by massively parallel sequencing (MAD-seq) capture assay and the MADSEQ analytical approach that allow low (<10%) levels of mosaicism for chromosomal aneuploidy or regional loss of heterozygosity to be detected, assigned to a meiotic or mitotic origin, and quantified as a proportion of the cells in the sample. We show results from a multi-ethnic MAD-seq (meMAD-seq) capture design that works equally well in populations of diverse racial and ethnic origins and how the MADSEQ analytical approach can be applied to exome or whole-genome sequencing data, revealing previously unrecognized aneuploidy or copy number neutral loss of heterozygosity in samples studied by the 1000 Genomes Project, cell lines from public repositories, and one of the Illumina Platinum Genomes samples. We have made the meMAD-seq capture design and MADSEQ analytical software open for unrestricted use, with the goal that they can be applied in clinical samples to allow new insights into the unrecognized prevalence of mosaic chromosomal aneuploidy in humans and its phenotypic associations.

Altered hydroxymethylation is seen at regulatory regions in pancreatic cancer and regulates oncogenic pathways.

Transcriptional deregulation of oncogenic pathways is a hallmark of cancer and can be due to epigenetic alterations. 5-Hydroxymethylcytosine (5-hmC) is an epigenetic modification that has not been studied in pancreatic cancer. Genome-wide analysis of 5-hmC-enriched loci with hmC-seal was conducted in a cohort of low-passage pancreatic cancer cell lines, primary patient-derived xenografts, and pancreatic controls and revealed strikingly altered patterns in neoplastic tissues. Differentially hydroxymethylated regions preferentially affected known regulatory regions of the genome, specifically overlapping with known H3K4me1 enhancers. Furthermore, base pair resolution analysis of cytosine methylation and hydroxymethylation with oxidative bisulfite sequencing was conducted and correlated with chromatin accessibility by ATAC-seq and gene expression by RNA-seq in pancreatic cancer and control samples. 5-hmC was specifically enriched at open regions of chromatin, and gain of 5-hmC was correlated with up-regulation of the cognate transcripts, including many oncogenic pathways implicated in pancreatic neoplasia, such as MYC, KRAS, VEGFA, and BRD4 Specifically, BRD4 was overexpressed and acquired 5-hmC at enhancer regions in the majority of neoplastic samples. Functionally, acquisition of 5-hmC at BRD4 promoter was associated with increase in transcript expression in reporter assays and primary samples. Furthermore, blockade of BRD4 inhibited pancreatic cancer growth in vivo. In summary, redistribution of 5-hmC and preferential enrichment at oncogenic enhancers is a novel regulatory mechanism in human pancreatic cancer.

The nutritional environment determines which and how intestinal stem cells contribute to homeostasis and tumorigenesis.

Sporadic colon cancer accounts for approximately 80% of colorectal cancer (CRC) with high incidence in Western societies strongly linked to long-term dietary patterns. A unique mouse model for sporadic CRC results from feeding a purified rodent Western-style diet (NWD1) recapitulating intake for the mouse of common nutrient risk factors each at its level consumed in higher risk Western populations. This causes sporadic large and small intestinal tumors in wild-type mice at an incidence and frequency similar to that in humans. NWD1 perturbs intestinal cell maturation and Wnt signaling throughout villi and colonic crypts and decreases mouse Lgr5hi intestinal stem cell contribution to homeostasis and tumor development. Here we establish that NWD1 transcriptionally reprograms Lgr5hi cells, and that nutrients are interactive in reprogramming. Furthermore, the DNA mismatch repair pathway is elevated in Lgr5hi cells by lower vitamin D3 and/or calcium in NWD1, paralleled by reduced accumulation of relevant somatic mutations detected by single-cell exome sequencing. In compensation, NWD1 also reprograms Bmi1+ cells to function and persist as stem-like cells in mucosal homeostasis and tumor development. The data establish the key role of the nutrient environment in defining the contribution of two different stem cell populations to both mucosal homeostasis and tumorigenesis. This raises important questions regarding impact of variable human diets on which and how stem cell populations function in the human mucosa and give rise to tumors. Moreover, major differences reported in turnover of human and mouse crypt base stem cells may be linked to their very different nutrient exposures.

Mosaic epigenetic dysregulation of ectodermal cells in autism spectrum disorder.

DNA mutational events are increasingly being identified in autism spectrum disorder (ASD), but the potential additional role of dysregulation of the epigenome in the pathogenesis of the condition remains unclear. The epigenome is of interest as a possible mediator of environmental effects during development, encoding a cellular memory reflected by altered function of progeny cells. Advanced maternal age (AMA) is associated with an increased risk of having a child with ASD for reasons that are not understood. To explore whether AMA involves covert aneuploidy or epigenetic dysregulation leading to ASD in the offspring, we tested a homogeneous ectodermal cell type from 47 individuals with ASD compared with 48 typically developing (TD) controls born to mothers of ≥35 years, using a quantitative genome-wide DNA methylation assay. We show that DNA methylation patterns are dysregulated in ectodermal cells in these individuals, having accounted for confounding effects due to subject age, sex and ancestral haplotype. We did not find mosaic aneuploidy or copy number variability to occur at differentially-methylated regions in these subjects. Of note, the loci with distinctive DNA methylation were found at genes expressed in the brain and encoding protein products significantly enriched for interactions with those produced by known ASD-causing genes, representing a perturbation by epigenomic dysregulation of the same networks compromised by DNA mutational mechanisms. The results indicate the presence of a mosaic subpopulation of epigenetically-dysregulated, ectodermally-derived cells in subjects with ASD. The epigenetic dysregulation observed in these ASD subjects born to older mothers may be associated with aging parental gametes, environmental influences during embryogenesis or could be the consequence of mutations of the chromatin regulatory genes increasingly implicated in ASD. The results indicate that epigenetic dysregulatory mechanisms may complement and interact with DNA mutations in the pathogenesis of the disorder.

Myeloma is characterized by stage-specific alterations in DNA methylation that occur early during myelomagenesis.

Epigenetic changes play important roles in carcinogenesis and influence initial steps in neoplastic transformation by altering genome stability and regulating gene expression. To characterize epigenomic changes during the transformation of normal plasma cells to myeloma, we modified the HpaII tiny fragment enrichment by ligation-mediated PCR assay to work with small numbers of purified primary marrow plasma cells. The nano-HpaII tiny fragment enrichment by ligation-mediated PCR assay was used to analyze the methylome of CD138(+) cells from 56 subjects representing premalignant (monoclonal gammopathy of uncertain significance), early, and advanced stages of myeloma, as well as healthy controls. Plasma cells from premalignant and early stages of myeloma were characterized by striking, widespread hypomethylation. Gene-specific hypermethylation was seen to occur in the advanced stages, and cell lines representative of relapsed cases were found to be sensitive to decitabine. Aberrant demethylation in monoclonal gammopathy of uncertain significance occurred primarily in CpG islands, whereas differentially methylated loci in cases of myeloma occurred predominantly outside of CpG islands and affected distinct sets of gene pathways, demonstrating qualitative epigenetic differences between premalignant and malignant stages. Examination of the methylation machinery revealed that the methyltransferase, DNMT3A, was aberrantly hypermethylated and underexpressed, but not mutated in myeloma. DNMT3A underexpression was also associated with adverse overall survival in a large cohort of patients, providing insights into genesis of hypomethylation in myeloma. These results demonstrate widespread, stage-specific epigenetic changes during myelomagenesis and suggest that early demethylation can be a potential contributor to genome instability seen in myeloma. We also identify DNMT3A expression as a novel prognostic biomarker and suggest that relapsed cases can be therapeutically targeted by hypomethylating agents.

High resolution methylome analysis reveals widespread functional hypomethylation during adult human erythropoiesis.

Differentiation of hematopoietic stem cells to red cells requires coordinated expression of numerous erythroid genes and is characterized by nuclear condensation and extrusion during terminal development. To understand the regulatory mechanisms governing these widespread phenotypic changes, we conducted a high resolution methylomic and transcriptomic analysis of six major stages of human erythroid differentiation. We observed widespread epigenetic differences between early and late stages of erythropoiesis with progressive loss of methylation being the dominant change during differentiation. Gene bodies, intergenic regions, and CpG shores were preferentially demethylated during erythropoiesis. Epigenetic changes at transcription factor binding sites correlated significantly with changes in gene expression and were enriched for binding motifs for SCL, MYB, GATA, and other factors not previously implicated in erythropoiesis. Demethylation at gene promoters was associated with increased expression of genes, whereas epigenetic changes at gene bodies correlated inversely with gene expression. Important gene networks encoding erythrocyte membrane proteins, surface receptors, and heme synthesis proteins were found to be regulated by DNA methylation. Furthermore, integrative analysis enabled us to identify novel, potential regulatory areas of the genome as evident by epigenetic changes in a predicted PU.1 binding site in intron 1 of the GATA1 gene. This intronic site was found to be conserved across species and was validated to be a novel PU.1 binding site by quantitative ChIP in erythroid cells. Altogether, our study provides a comprehensive analysis of methylomic and transcriptomic changes during erythroid differentiation and demonstrates that human terminal erythropoiesis is surprisingly associated with hypomethylation of the genome.

Direct, genome-wide assessment of DNA mutations in single cells.

DNA mutations are the inevitable consequences of errors that arise during replication and repair of DNA damage. Because of their random and infrequent occurrence, quantification and characterization of DNA mutations in the genome of somatic cells has been difficult. Random, low-abundance mutations are currently inaccessible by standard high-throughput sequencing approaches because they cannot be distinguished from sequencing errors. One way to circumvent this problem and simultaneously account for the mutational heterogeneity within tissues is whole genome sequencing of a representative number of single cells. Here, we show elevated mutation levels in single cells from Drosophila melanogaster S2 and mouse embryonic fibroblast populations after treatment with the powerful mutagen N-ethyl-N-nitrosourea. This method can be applied as a direct measure of exposure to mutagenic agents and for assessing genotypic heterogeneity within tissues or cell populations.

An optimized approach for multiplexing single-nuclear ATAC-seq using oligonucleotide-conjugated antibodies.

BACKGROUND: Single-cell technologies to analyze transcription and chromatin structure have been widely used in many research areas to reveal the functions and molecular properties of cells at single-cell resolution. Sample multiplexing techniques are valuable when performing single-cell analysis, reducing technical variation and permitting cost efficiencies. Several commercially available methods have been used in many scRNA-seq studies. On the other hand, while several methods have been published, multiplexing techniques for single nuclear assay for transposase-accessible chromatin (snATAC)-seq assays remain under development. We developed a simple nucleus hashing method using oligonucleotide-conjugated antibodies recognizing nuclear pore complex proteins, NuHash, to perform snATAC-seq library preparations by multiplexing. RESULTS: We performed multiplexing snATAC-seq analyses on a mixture of human and mouse cell samples (two samples, 2-plex, and four samples, 4-plex) using NuHash. The analyses on nuclei with at least 10,000 read counts showed that the demultiplexing accuracy of NuHash was high, and only ten out of 9144 nuclei (2-plex) and 150 of 12,208 nuclei (4-plex) had discordant classifications between NuHash demultiplexing and discrimination using reference genome alignments. The differential open chromatin region (OCR) analysis between female and male samples revealed that male-specific OCRs were enriched in chromosome Y (four out of nine). We also found that five female-specific OCRs (20 OCRs) were on chromosome X. A comparative analysis between snATAC-seq and deeply sequenced bulk ATAC-seq on the same samples revealed that the bulk ATAC-seq signal intensity was positively correlated with the number of cell clusters detected in snATAC-seq. Moreover, when we categorized snATAC-seq peaks based on the number of cell clusters in which the peak was present, we observed different distributions over different genomic features between the groups. This result suggests that the peak intensities of bulk ATAC-seq can be used to identify different types of functional loci. CONCLUSIONS: Our multiplexing method using oligo-conjugated anti-nuclear pore complex proteins, NuHash, permits high-accuracy demultiplexing of samples. The NuHash protocol is straightforward, works on frozen samples, and requires no modifications for snATAC-seq library preparation.

Dampened activity of E2F1-DP and Myb-MuvB transcription factors in Drosophila endocycling cells.

The endocycle is a variant cell cycle comprised of alternating gap (G) and DNA synthesis (S) phases (endoreplication) without mitosis (M), which results in DNA polyploidy and large cell size. Endocycles occur widely in nature, but much remains to be learned about the regulation of this modified cell cycle. Here, we compared gene expression profiles of mitotic cycling larval brain and disc cells with the endocycling cells of fat body and salivary gland of the Drosophila larva. The results indicated that many genes that are positively regulated by the heterodimeric E2F1-DP or Myb-MuvB complex transcription factors are expressed at lower levels in endocycling cells. Many of these target genes have functions in M phase, suggesting that dampened E2F1 and Myb activity promote endocycles. Many other E2F1 target genes that are required for DNA replication were also repressed in endocycling cells, an unexpected result given that these cells must duplicate up to thousands of genome copies during each S phase. For some EF2-regulated genes, the lower level of mRNA in endocycling cells resulted in lower protein concentration, whereas for other genes it did not, suggesting a contribution of post-transcriptional regulation. Both knockdown and overexpression of E2F1-DP and Myb-MuvB impaired endocycles, indicating that transcriptional activation and repression must be balanced. Our data suggest that dampened transcriptional activation by E2F1-DP and Myb-MuvB is important to repress mitosis and coordinate the endocycle transcriptional and protein stability oscillators.

ETS1 is a novel transcriptional regulator of adult T-cell leukemia/lymphoma of North American descent.

Adult T-cell leukemia/lymphoma (ATLL) is an aggressive T-cell lymphoma associated with the human T-cell lymphotropic virus type 1 virus endemic in regions including Japan, the Caribbean islands, and Latin America. Although progress has been made to understand the disease, survival outcomes with current standard therapy remain extremely poor particularly in acute ATLL, underlying the need for better understanding of its biology and identification of novel therapeutic targets. Recently, it was demonstrated that ATLL of North American-descendent patients (NA-ATLL) is both clinically and molecularly distinct from Japanese-descendent (J-ATLL), with inferior prognosis and higher incidence of epigenetic-targeting mutations compared with J-ATLL. In this study, combined chromatin accessibility and transcriptomic profiling were used to further understand the key transcriptional regulators of NA-ATLL compared with J-ATLL. The ETS1 motif was found to be enriched in chromatin regions that were differentially open in NA-ATLL, whereas the AP1/IRF4 motifs were enriched in chromatin regions more open in J-ATLL. ETS1 expression was markedly elevated in NA-ATLL in both cell line and primary tumor samples, and knockdown of ETS1 in NA-ATLL cells resulted in inhibition of cell growth. CCR4, a previously identified oncogenic factor in ATLL, was found to be a direct ETS1 transcriptional target in NA-ATLL. As such, ETS1 provides an alternate mechanism to enhance CCR4 expression/activity in NA-ATLL, even in the absence of activating CCR4 mutations (CCR4 mutations were identified in 4 of 9 NA-ATLL cases). Taken together, this study identifies ETS1 as a novel dominant oncogenic transcriptional regulator in NA-ATLL.

Genomic insights into host and parasite interactions during intracellular infection by Toxoplasma gondii.

To gain insights into the molecular interactions of an intracellular pathogen and its host cell, we studied the gene expression and chromatin states of human fibroblasts infected with the Apicomplexan parasite Toxoplasma gondii. We show a striking activation of host cell genes that regulate a number of cellular processes, some of which are protective of the host cell, others likely to be advantageous to the pathogen. The simultaneous capture of host and parasite genomic information allowed us to gain insights into the regulation of the T. gondii genome. We show how chromatin accessibility and transcriptional profiling together permit novel annotation of the parasite's genome, including more accurate mapping of known genes and the identification of new genes and cis-regulatory elements. Motif analysis reveals not only the known T. gondii AP2 transcription factor-binding site but also a previously-undiscovered candidate TATA box-containing motif at one-quarter of promoters. By inferring the transcription factor and upstream cell signaling responses involved in the host cell, we can use genomic information to gain insights into T. gondii's perturbation of host cell physiology. Our resulting model builds on previously-described human host cell signalling responses to T. gondii infection, linked to induction of specific transcription factors, some of which appear to be solely protective of the host cell, others of which appear to be co-opted by the pathogen to enhance its own survival.

ASXL1 mutations are associated with distinct epigenomic alterations that lead to sensitivity to venetoclax and azacytidine.

The BCL2-inhibitor, Venetoclax (VEN), has shown significant anti-leukemic efficacy in combination with the DNMT-inhibitor, Azacytidine (AZA). To explore the mechanisms underlying the selective sensitivity of mutant leukemia cells to VEN and AZA, we used cell-based isogenic models containing a common leukemia-associated mutation in the epigenetic regulator ASXL1. KBM5 cells with CRISPR/Cas9-mediated correction of the ASXL1G710X mutation showed reduced leukemic growth, increased myeloid differentiation, and decreased HOXA and BCL2 gene expression in vitro compared to uncorrected KBM5 cells. Increased expression of the anti-apoptotic gene, BCL2, was also observed in bone marrow CD34+ cells from ASXL1 mutant MDS patients compared to CD34+ cells from wild-type MDS cases. ATAC-sequencing demonstrated open chromatin at the BCL2 promoter in the ASXL1 mutant KBM5 cells. BH3 profiling demonstrated increased dependence of mutant cells on BCL2. Upon treatment with VEN, mutant cells demonstrated increased growth inhibition. In addition, genome-wide methylome analysis of primary MDS samples and isogenic cell lines demonstrated increased gene-body methylation in ASXL1 mutant cells, with consequently increased sensitivity to AZA. These data mechanistically link the common leukemia-associated mutation ASXL1 to enhanced sensitivity to VEN and AZA via epigenetic upregulation of BCL2 expression and widespread alterations in DNA methylation.

A Cellular Stress Response Induced by the CRISPR-dCas9 Activation System Is Not Heritable Through Cell Divisions.

The CRISPR-Cas9 system can be modified to perform "epigenetic editing" by utilizing the catalytically inactive (dead) Cas9 (dCas9) to recruit regulatory proteins to specific genomic locations. In prior studies, epigenetic editing with multimers of the transactivator VP16 and guide RNAs (gRNAs) was found to cause adverse cellular responses. These side effects may confound studies inducing new cellular properties, especially if the cellular responses are maintained through cell divisions-an epigenetic regulatory property. Here, we show how distinct components of this CRISPR-dCas9 activation system, particularly dCas9 with untargeted gRNAs, upregulate genes associated with transcriptional stress, defense response, and regulation of cell death. Our results highlight a previously undetected acute stress response to CRISPR-dCas9 components in human cells, which is transient and not maintained through cell divisions.

Transformation of oats and its application to improving osmotic stress tolerance.

Oat (Avena sativa L.), a worldwide temperate cereal crop, is deficient in tolerance to osmotic stress due to drought and/or salinity. To genetically transform the available commercial oat cultivars, a genotype-independent and efficient regeneration system from shoot apical meristems was developed using four oat cultivars: Prairie, Porter, Ogle, and Pacer. All these oat cultivars generated a genotype-independent in vitro differentiated multiple shoots from shoot apical meristems at a high frequency. Using this system, three oat cultivars were genetically co-transformed with pBY520 (containing hva1 and bar) and pAct1-D (containing gus) using biolistic trade mark bombardment. Transgenic plants were selected and regenerated using herbicide resistance and GUS as a marker. Molecular and biochemical analyses of putative transgenic plants confirmed the co-integration of hva1 and bar genes with a frequency of 100%, and 61.6% of the transgenic plants carried all three genes (hva1, bar and gus). Further analyses of R0, R1, and R2 progenies confirmed stable integration, expression, and Mendalian inheritance for all transgenes. Histochemical analysis of GUS protein in transgenic plants showed a high level of GUS expression in vascular tissues and in the pollen grains of mature flowers. Immunochemical analysis of transgenic plants indicated a constitutive expression of hva1 at all developmental stages. However, the level of HVA1 was higher during the early seedling stages. The characteristic of HVA1 expression for osmotic tolerance in transgenic oat progeny was analyzed in vitro as well as in vivo. Transgenic plants exhibited significantly (P<0.05) increased tolerance to stress conditions than non-transgenic control plants. The symptoms of wilting or death of leaves as observed in 80% of non-transgenic plants due to osmotic stress was delayed and detected only in less than 10% of trans-genic plants. These observations confirmed the characteristic of HVA1 protein as providing or enhancing the osmotic tolerance in transgenic plants against salinity and possible water-deficiency stress conditions.

Misidentification of MLL3 and other mutations in cancer due to highly homologous genomic regions.

The MLL3 gene has been shown to be recurrently mutated in many malignancies including in families with acute myeloid leukemia. We demonstrate that many MLL3 variant calls made by exome sequencing are false positives due to misalignment to homologous regions, including a region on chr21, and can only be validated by long-range PCR. Numerous other recurrently mutated genes reported in COSMIC and TCGA databases have pseudogenes and cannot also be validated by conventional short read-based sequencing approaches. Genome-wide identification of pseudogene regions demonstrates that frequency of these homologous regions is increased with sequencing read lengths below 200 bps. To enable identification of poor quality sequencing variants in prospective studies, we generated novel genome-wide maps of regions with poor mappability that can be used in variant calling algorithms. Taken together, our findings reveal that pseudogene regions are a source of false-positive mutations in cancers.

Detection of circulating tumor DNA in patients with osteosarcoma.

Identification and quantification of somatic alterations in plasma-derived, circulating tumor DNA (ctDNA) is gaining traction as a non-invasive and cost effective method of disease monitoring in cancer patients, particularly to evaluate response to treatment and monitor for disease recurrence. To our knowledge, genetic analysis of ctDNA in osteosarcoma has not yet been studied. To determine whether somatic alterations can be detected in ctDNA and perhaps applied to patient management in this disease, we collected germline, tumor, and serial plasma samples from pediatric, adolescent, and young adult patients with osteosarcoma and used targeted Next Generation Sequencing (NGS) to identify somatic single nucleotide variants (SNV), insertions and deletions (INDELS), and structural variants (SV) in 7 genes commonly mutated in osteosarcoma. We demonstrate that patient-specific somatic alterations identified through comparison of tumor-germline pairs can be detected and quantified in cell-free DNA of osteosarcoma patients.

RNASeq in C. elegans Following Manganese Exposure.

Manganese is a metal that is required for optimal biological functioning of organisms. Absorption, cellular import and export, and excretion of manganese are all tightly regulated. While some genes involved in regulation, such as DMT-1 and ferroportin, are known, it is presumed that many more are involved and as yet unknown. Excessive exposure to manganese, usually in industrial settings such as mining or welding, can lead to neurotoxicity and a condition known as manganism that closely resembles Parkinson's disease. Elucidating transcriptional changes following manganese exposure could lead to the development of biomarkers for exposure. This unit presents a protocol for RNA sequencing in the worm Caenorhabditis elegans to assay for transcriptional changes following exposure to manganese. This protocol is adaptable to any environmental exposure in C. elegans. The protocol results in counts of gene transcripts in control versus exposed conditions and a ranked list of differentially expressed genes for further study.

The Complete Genome Sequence of the Emerging Pathogen Mycobacterium haemophilum Explains Its Unique Culture Requirements.

UNLABELLED: Mycobacterium haemophilum is an emerging pathogen associated with a variety of clinical syndromes, most commonly skin infections in immunocompromised individuals. M. haemophilum exhibits a unique requirement for iron supplementation to support its growth in culture, but the basis for this property and how it may shape pathogenesis is unclear. Using a combination of Illumina, PacBio, and Sanger sequencing, the complete genome sequence of M. haemophilum was determined. Guided by this sequence, experiments were performed to define the basis for the unique growth requirements of M. haemophilum. We found that M. haemophilum, unlike many other mycobacteria, is unable to synthesize iron-binding siderophores known as mycobactins or to utilize ferri-mycobactins to support growth. These differences correlate with the absence of genes associated with mycobactin synthesis, secretion, and uptake. In agreement with the ability of heme to promote growth, we identified genes encoding heme uptake machinery. Consistent with its propensity to infect the skin, we show at the whole-genome level the genetic closeness of M. haemophilum with Mycobacterium leprae, an organism which cannot be cultivated in vitro, and we identify genes uniquely shared by these organisms. Finally, we identify means to express foreign genes in M. haemophilum. These data explain the unique culture requirements for this important pathogen, provide a foundation upon which the genome sequence can be exploited to improve diagnostics and therapeutics, and suggest use of M. haemophilum as a tool to elucidate functions of genes shared with M. leprae. IMPORTANCE: Mycobacterium haemophilum is an emerging pathogen with an unknown natural reservoir that exhibits unique requirements for iron supplementation to grow in vitro. Understanding the basis for this iron requirement is important because it is fundamental to isolation of the organism from clinical samples and environmental sources. Defining the molecular basis for M. haemophilium's growth requirements will also shed new light on mycobacterial strategies to acquire iron and can be exploited to define how differences in such strategies influence pathogenesis. Here, through a combination of sequencing and experimental approaches, we explain the basis for the iron requirement. We further demonstrate the genetic closeness of M. haemophilum and Mycobacterium leprae, the causative agent of leprosy which cannot be cultured in vitro, and we demonstrate methods to genetically manipulate M. haemophilum. These findings pave the way for the use of M. haemophilum as a model to elucidate functions of genes shared with M. leprae.

Transcriptional and metabolic signatures of Arabidopsis responses to chewing damage by an insect herbivore and bacterial infection and the consequences of their interaction.

Plants use multiple interacting signaling systems to identify and respond to biotic stresses. Although it is often assumed that there is specificity in signaling responses to specific pests, this is rarely examined outside of the gene-for-gene relationships of plant-pathogen interactions. In this study, we first compared early events in gene expression and later events in metabolite profiles of Arabidopsis thaliana following attack by either the caterpillar Spodoptera exigua or avirulent (DC3000 avrRpm1) Pseudomonas syringae pv. tomato at three time points. Transcriptional responses of the plant to caterpillar feeding were rapid, occurring within 1 h of feeding, and then decreased at 6 and 24 h. In contrast, plant response to the pathogen was undetectable at 1 h but grew larger and more significant at 6 and 24 h. There was a surprisingly large amount of overlap in jasmonate and salicylate signaling in responses to the insect and pathogen, including levels of gene expression and individual hormones. The caterpillar and pathogen treatments induced different patterns of expression of glucosinolate biosynthesis genes and levels of glucosinolates. This suggests that when specific responses develop, their regulation is complex and best understood by characterizing expression of many genes and metabolites. We then examined the effect of feeding by the caterpillar Spodoptera exigua on Arabidopsis susceptibility to virulent (DC3000) and avirulent (DC3000 avrRpm1) P. syringae pv. tomato, and found that caterpillar feeding enhanced Arabidopsis resistance to the avirulent pathogen and lowered resistance to the virulent strain. We conclude that efforts to improve plant resistance to bacterial pathogens are likely to influence resistance to insects and vice versa. Studies explicitly comparing plant responses to multiple stresses, including the role of elicitors at early time points, are critical to understanding how plants organize responses in natural settings.

RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders.

Genome-wide expression analysis using next generation sequencing (RNA-Seq) provides an opportunity for in-depth molecular profiling of fundamental biological processes, such as cellular differentiation and malignant transformation. Differentiating human neurons derived from induced pluripotent stem cells (iPSCs) provide an ideal system for RNA-Seq since defective neurogenesis caused by abnormalities in transcription factors, DNA methylation, and chromatin modifiers lie at the heart of some neuropsychiatric disorders. As a preliminary step towards applying next generation sequencing using neurons derived from patient-specific iPSCs, we have carried out an RNA-Seq analysis on control human neurons. Dramatic changes in the expression of coding genes, long non-coding RNAs (lncRNAs), pseudogenes, and splice isoforms were seen during the transition from pluripotent stem cells to early differentiating neurons. A number of genes that undergo radical changes in expression during this transition include candidates for schizophrenia (SZ), bipolar disorder (BD) and autism spectrum disorders (ASD) that function as transcription factors and chromatin modifiers, such as POU3F2 and ZNF804A, and genes coding for cell adhesion proteins implicated in these conditions including NRXN1 and NLGN1. In addition, a number of novel lncRNAs were found to undergo dramatic changes in expression, one of which is HOTAIRM1, a regulator of several HOXA genes during myelopoiesis. The increase we observed in differentiating neurons suggests a role in neurogenesis as well. Finally, several lncRNAs that map near SNPs associated with SZ in genome wide association studies also increase during neuronal differentiation, suggesting that these novel transcripts may be abnormally regulated in a subgroup of patients.