Transcription factors MYOCD, SRF, Mesp1 and SMARCD3 enhance the cardio-inducing effect of GATA4, TBX5, and MEF2C during direct cellular reprogramming.

Transient overexpression of defined combinations of master regulator genes can effectively induce cellular reprogramming: the acquisition of an alternative predicted phenotype from a differentiated cell lineage. This can be of particular importance in cardiac regenerative medicine wherein the heart lacks the capacity to heal itself, but simultaneously contains a large pool of fibroblasts. In this study we determined the cardio-inducing capacity of ten transcription factors to actuate cellular reprogramming of mouse embryonic fibroblasts into cardiomyocyte-like cells. Overexpression of transcription factors MYOCD and SRF alone or in conjunction with Mesp1 and SMARCD3 enhanced the basal but necessary cardio-inducing effect of the previously reported GATA4, TBX5, and MEF2C. In particular, combinations of five or seven transcription factors enhanced the activation of cardiac reporter vectors, and induced an upregulation of cardiac-specific genes. Global gene expression analysis also demonstrated a significantly greater cardio-inducing effect when the transcription factors MYOCD and SRF were used. Detection of cross-striated cells was highly dependent on the cell culture conditions and was enhanced by the addition of valproic acid and JAK inhibitor. Although we detected Ca(2+) transient oscillations in the reprogrammed cells, we did not detect significant changes in resting membrane potential or spontaneously contracting cells. This study further elucidates the cardio-inducing effect of the transcriptional networks involved in cardiac cellular reprogramming, contributing to the ongoing rational design of a robust protocol required for cardiac regenerative therapies.

A host transcriptional signature for presymptomatic detection of infection in humans exposed to influenza H1N1 or H3N2.

There is great potential for host-based gene expression analysis to impact the early diagnosis of infectious diseases. In particular, the influenza pandemic of 2009 highlighted the challenges and limitations of traditional pathogen-based testing for suspected upper respiratory viral infection. We inoculated human volunteers with either influenza A (A/Brisbane/59/2007 (H1N1) or A/Wisconsin/67/2005 (H3N2)), and assayed the peripheral blood transcriptome every 8 hours for 7 days. Of 41 inoculated volunteers, 18 (44%) developed symptomatic infection. Using unbiased sparse latent factor regression analysis, we generated a gene signature (or factor) for symptomatic influenza capable of detecting 94% of infected cases. This gene signature is detectable as early as 29 hours post-exposure and achieves maximal accuracy on average 43 hours (p = 0.003, H1N1) and 38 hours (p-value = 0.005, H3N2) before peak clinical symptoms. In order to test the relevance of these findings in naturally acquired disease, a composite influenza A signature built from these challenge studies was applied to Emergency Department patients where it discriminates between swine-origin influenza A/H1N1 (2009) infected and non-infected individuals with 92% accuracy. The host genomic response to Influenza infection is robust and may provide the means for detection before typical clinical symptoms are apparent.

Integrative Genomics Reveals a Role for GNA13 in Lymphomagenesis

Lymphomas comprise a diverse group of malignancies derived from immune cells. High throughput sequencing has recently emerged as a powerful and versatile method for analysis of the cancer genome and transcriptome. As these data continue to emerge, the crucial work lies in sorting through the wealth of information to hone in on the critical aspects that will give us a better understanding of biology and new insight for how to treat disease. Finding the important signals within these large data sets is one of the major challenges of next generation sequencing.

In this dissertation, I have developed several complementary strategies to describe the genetic underpinnings of lymphomas. I begin with developing a better method for RNA sequencing that enables strand-specific total RNA sequencing and alternative splicing profiling in the same analysis. I then combine this RNA sequencing technique with whole exome sequencing to better understand the global landscape of aberrations in these diseases. Finally, I use traditional cell and molecular biology techniques to define the consequences of major genetic alterations in lymphoma.

Through this analysis, I find recurrent silencing mutations in the G alpha binding protein GNA13 and associated focal adhesion proteins. I aim to describe how loss-of-function mutations in GNA13 can be oncogenic in the context of germinal center B cell biology. Using in vitro techniques including liquid chromatography-mass spectrometry and knockdown and overexpression of genes in B cell lymphoma cell lines, I determine protein binding partners and downstream effectors of GNA13. I also develop a transgenic mouse model to study the role of GNA13 in the germinal center in vivo to determine effects of GNA13 deletion on germinal center structure and cell migration.

Thus, I have developed complementary approaches that span the spectrum from discovery to context-dependent gene models that afford a better understanding of the biological function of aberrant events and ultimately result in a better understanding of disease.

Recovery from an acute infection in C. elegans requires the GATA transcription factor ELT-2.

The mechanisms involved in the recognition of microbial pathogens and activation of the immune system have been extensively studied. However, the mechanisms involved in the recovery phase of an infection are incompletely characterized at both the cellular and physiological levels. Here, we establish a Caenorhabditis elegans-Salmonella enterica model of acute infection and antibiotic treatment for studying biological changes during the resolution phase of an infection. Using whole genome expression profiles of acutely infected animals, we found that genes that are markers of innate immunity are down-regulated upon recovery, while genes involved in xenobiotic detoxification, redox regulation, and cellular homeostasis are up-regulated. In silico analyses demonstrated that genes altered during recovery from infection were transcriptionally regulated by conserved transcription factors, including GATA/ELT-2, FOXO/DAF-16, and Nrf/SKN-1. Finally, we found that recovery from an acute bacterial infection is dependent on ELT-2 activity.

Hybrid error correction and de novo assembly of single-molecule sequencing reads.

Single-molecule sequencing instruments can generate multikilobase sequences with the potential to greatly improve genome and transcriptome assembly. However, the error rates of single-molecule reads are high, which has limited their use thus far to resequencing bacteria. To address this limitation, we introduce a correction algorithm and assembly strategy that uses short, high-fidelity sequences to correct the error in single-molecule sequences. We demonstrate the utility of this approach on reads generated by a PacBio RS instrument from phage, prokaryotic and eukaryotic whole genomes, including the previously unsequenced genome of the parrot Melopsittacus undulatus, as well as for RNA-Seq reads of the corn (Zea mays) transcriptome. Our long-read correction achieves >99.9% base-call accuracy, leading to substantially better assemblies than current sequencing strategies: in the best example, the median contig size was quintupled relative to high-coverage, second-generation assemblies. Greater gains are predicted if read lengths continue to increase, including the prospect of single-contig bacterial chromosome assembly.

Mitochondrial mutations in adenoid cystic carcinoma of the salivary glands.

BACKGROUND: The MitoChip v2.0 resequencing array is an array-based technique allowing for accurate and complete sequencing of the mitochondrial genome. No studies have investigated mitochondrial mutation in salivary gland adenoid cystic carcinomas. METHODOLOGY: The entire mitochondrial genome of 22 salivary gland adenoid cystic carcinomas (ACC) of salivary glands and matched leukocyte DNA was sequenced to determine the frequency and distribution of mitochondrial mutations in ACC tumors. PRINCIPAL FINDINGS: Seventeen of 22 ACCs (77%) carried mitochondrial mutations, ranging in number from 1 to 37 mutations. A disproportionate number of mutations occurred in the D-loop. Twelve of 17 tumors (70.6%) carried mutations resulting in amino acid changes of translated proteins. Nine of 17 tumors (52.9%) with a mutation carried an amino acid changing mutation in the nicotinamide adenine dinucleotide dehydrogenase (NADH) complex. CONCLUSIONS/SIGNIFICANCE: Mitochondrial mutation is frequent in salivary ACCs. The high incidence of amino acid changing mutations implicates alterations in aerobic respiration in ACC carcinogenesis. D-loop mutations are of unclear significance, but may be associated with alterations in transcription or replication.

Metabolic programming and PDHK1 control CD4+ T cell subsets and inflammation.

Activation of CD4+ T cells results in rapid proliferation and differentiation into effector and regulatory subsets. CD4+ effector T cell (Teff) (Th1 and Th17) and Treg subsets are metabolically distinct, yet the specific metabolic differences that modify T cell populations are uncertain. Here, we evaluated CD4+ T cell populations in murine models and determined that inflammatory Teffs maintain high expression of glycolytic genes and rely on high glycolytic rates, while Tregs are oxidative and require mitochondrial electron transport to proliferate, differentiate, and survive. Metabolic profiling revealed that pyruvate dehydrogenase (PDH) is a key bifurcation point between T cell glycolytic and oxidative metabolism. PDH function is inhibited by PDH kinases (PDHKs). PDHK1 was expressed in Th17 cells, but not Th1 cells, and at low levels in Tregs, and inhibition or knockdown of PDHK1 selectively suppressed Th17 cells and increased Tregs. This alteration in the CD4+ T cell populations was mediated in part through ROS, as N-acetyl cysteine (NAC) treatment restored Th17 cell generation. Moreover, inhibition of PDHK1 modulated immunity and protected animals against experimental autoimmune encephalomyelitis, decreasing Th17 cells and increasing Tregs. Together, these data show that CD4+ subsets utilize and require distinct metabolic programs that can be targeted to control specific T cell populations in autoimmune and inflammatory diseases.

Nutritional control of mRNA isoform expression during developmental arrest and recovery in C. elegans.

Nutrient availability profoundly influences gene expression. Many animal genes encode multiple transcript isoforms, yet the effect of nutrient availability on transcript isoform expression has not been studied in genome-wide fashion. When Caenorhabditis elegans larvae hatch without food, they arrest development in the first larval stage (L1 arrest). Starved larvae can survive L1 arrest for weeks, but growth and post-embryonic development are rapidly initiated in response to feeding. We used RNA-seq to characterize the transcriptome during L1 arrest and over time after feeding. Twenty-seven percent of detectable protein-coding genes were differentially expressed during recovery from L1 arrest, with the majority of changes initiating within the first hour, demonstrating widespread, acute effects of nutrient availability on gene expression. We used two independent approaches to track expression of individual exons and mRNA isoforms, and we connected changes in expression to functional consequences by mining a variety of databases. These two approaches identified an overlapping set of genes with alternative isoform expression, and they converged on common functional patterns. Genes affecting mRNA splicing and translation are regulated by alternative isoform expression, revealing post-transcriptional consequences of nutrient availability on gene regulation. We also found that phosphorylation sites are often alternatively expressed, revealing a common mode by which alternative isoform expression modifies protein function and signal transduction. Our results detail rich changes in C. elegans gene expression as larvae initiate growth and post-embryonic development, and they provide an excellent resource for ongoing investigation of transcriptional regulation and developmental physiology.

The evolutionary history of ferns inferred from 25 low-copy nuclear genes.

UNLABELLED: • PREMISE OF THE STUDY: Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data.• METHODS: Here we take a curated phylogenomics approach to infer the first broad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35877 bp), along with rigorous alignment, orthology inference and model selection.• KEY RESULTS: Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we find strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ∼200 MYA. Finally, our species-tree inferences are congruent with analyses of concatenated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data.• CONCLUSIONS: Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.

High-coverage sequencing and annotated assemblies of the budgerigar genome.

BACKGROUND: Parrots belong to a group of behaviorally advanced vertebrates and have an advanced ability of vocal learning relative to other vocal-learning birds. They can imitate human speech, synchronize their body movements to a rhythmic beat, and understand complex concepts of referential meaning to sounds. However, little is known about the genetics of these traits. Elucidating the genetic bases would require whole genome sequencing and a robust assembly of a parrot genome. FINDINGS: We present a genomic resource for the budgerigar, an Australian Parakeet (Melopsittacus undulatus) -- the most widely studied parrot species in neuroscience and behavior. We present genomic sequence data that includes over 300× raw read coverage from multiple sequencing technologies and chromosome optical maps from a single male animal. The reads and optical maps were used to create three hybrid assemblies representing some of the largest genomic scaffolds to date for a bird; two of which were annotated based on similarities to reference sets of non-redundant human, zebra finch and chicken proteins, and budgerigar transcriptome sequence assemblies. The sequence reads for this project were in part generated and used for both the Assemblathon 2 competition and the first de novo assembly of a giga-scale vertebrate genome utilizing PacBio single-molecule sequencing. CONCLUSIONS: Across several quality metrics, these budgerigar assemblies are comparable to or better than the chicken and zebra finch genome assemblies built from traditional Sanger sequencing reads, and are sufficient to analyze regions that are difficult to sequence and assemble, including those not yet assembled in prior bird genomes, and promoter regions of genes differentially regulated in vocal learning brain regions. This work provides valuable data and material for genome technology development and for investigating the genomics of complex behavioral traits.

Three crocodilian genomes reveal ancestral patterns of evolution among archosaurs.

To provide context for the diversification of archosaurs--the group that includes crocodilians, dinosaurs, and birds--we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.

Sea Urchin Body Plan Development and Evolution: An Integrative Transcriptomic Approach

My dissertation work integrates comparative transcriptomics and functional analyses to investigate gene expression changes underlying two significant aspects of sea urchin evolution and development: the dramatic developmental changes associated with an ecologically significant shift in life history strategy and the development of the unusual radial body plan of adult sea urchins.

In Chapter 2, I investigate evolutionary changes in gene expression underlying the switch from feeding (planktotrophic) to nonfeeding (lecithotrophic) development in sea urchins. In order to identify these changes, I used Illumina RNA-seq to measure expression dynamics across 7 developmental stages in three sea urchin species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph Heliocidaris tuberculata, and an outgroup planktotroph Lytechinus variegatus. My analyses draw on a well-characterized developmental gene regulatory network (GRN) in sea urchins to understand how the ancestral planktotrophic developmental program was altered during the evolution of lecithotrophic development. My results suggest that changes in gene expression profiles occurred more frequently across the transcriptome during the evolution of lecithotrophy than during the persistence of planktotrophy. These changes were even more pronounced within the GRN than across the transcriptome as a whole, and occurred in each network territory (skeletogenic, endomesoderm and ectoderm). I found evidence for both conservation and divergence of regulatory interactions in the network, as well as significant changes in the expression of genes with known roles in larval skeletogenesis, which is dramatically altered in lecithotrophs. I further explored network dynamics between species using coexpression analyses, which allowed me to identify novel players likely involved in sea urchin neurogenesis and endoderm patterning.

In Chapter 3, I investigate developmental changes in gene expression underlying radial body plan development and metamorphosis in H. erythrogramma. Using Illumina RNA-seq, I measured gene expression profiles across larval, metamorphic, and post-metamorphic life cycle phases. My results present a high-resolution view of gene expression dynamics during the complex transition from pre- to post-metamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases.

Collectively, my investigations provide an important foundation for future, empirical studies to investigate the functional role of gene expression change in the evolution of developmental differences between species and also for the generation of the unusual radial body plan of sea urchins.

A computational analysis of antisense off-targets in prokaryotic organisms

© 2014 .The adoption of antisense gene silencing as a novel disinfectant for prokaryotic organisms is hindered by poor silencing efficiencies. Few studies have considered the effects of off-targets on silencing efficiencies, especially in prokaryotic organisms. In this computational study, a novel algorithm was developed that determined and sorted the number of off-targets as a function of alignment length in Escherichia coli K-12 MG1655 and Mycobacterium tuberculosis H37Rv. The mean number of off-targets per a single location was calculated to be 14.1. ±. 13.3 and 36.1. ±. 58.5 for the genomes of E. coli K-12 MG1655 and M. tuberculosis H37Rv, respectively. Furthermore, when the entire transcriptome was analyzed, it was found that there was no general gene location that could be targeted to minimize or maximize the number of off-targets. In an effort to determine the effects of off-targets on silencing efficiencies, previously published studies were used. Analyses with acpP, ino1, and marORAB revealed a statistically significant relationship between the number of short alignment length off-targets hybrids and the efficacy of the antisense gene silencing, suggesting that the minimization of off-targets may be beneficial for antisense gene silencing in prokaryotic organisms.

The Physarum polycephalum Genome Reveals Extensive Use of Prokaryotic Two-Component and Metazoan-Type Tyrosine Kinase Signaling.

Physarum polycephalum is a well-studied microbial eukaryote with unique experimental attributes relative to other experimental model organisms. It has a sophisticated life cycle with several distinct stages including amoebal, flagellated, and plasmodial cells. It is unusual in switching between open and closed mitosis according to specific life-cycle stages. Here we present the analysis of the genome of this enigmatic and important model organism and compare it with closely related species. The genome is littered with simple and complex repeats and the coding regions are frequently interrupted by introns with a mean size of 100 bases. Complemented with extensive transcriptome data, we define approximately 31,000 gene loci, providing unexpected insights into early eukaryote evolution. We describe extensive use of histidine kinase-based two-component systems and tyrosine kinase signaling, the presence of bacterial and plant type photoreceptors (phytochromes, cryptochrome, and phototropin) and of plant-type pentatricopeptide repeat proteins, as well as metabolic pathways, and a cell cycle control system typically found in more complex eukaryotes. Our analysis characterizes P. polycephalum as a prototypical eukaryote with features attributed to the last common ancestor of Amorphea, that is, the Amoebozoa and Opisthokonts. Specifically, the presence of tyrosine kinases in Acanthamoeba and Physarum as representatives of two distantly related subdivisions of Amoebozoa argues against the later emergence of tyrosine kinase signaling in the opisthokont lineage and also against the acquisition by horizontal gene transfer.

Characterization of the human and mouse ETV1/ER81 transcription factor genes: Role of the two alternatively spliced isoforms in the human

The Ets transcription factors of the PEA3 group - E1AF/PEA3, ETV1/ER81 and ERM - are almost identical in the ETS DNA-binding and the transcriptional acidic domains. To accelerate our understanding of the molecular basis of putative diseases linked to ETV1 such as Ewing's sarcoma we characterized the human ETV1 and the mouse ER81 genes. We showed that these genes are both encoded by 13 exons in more than 90 kbp genomic DNA, and that the classical acceptor and donor splicing sites are present in each junction except for the 5' donor site of intron 9 where GT is replaced by TT. The genomic organization of the ETS and acidic domains in the human ETV1 and mouse ER81 (localized to chromosome 12) genes is similar to that observed in human ERM and human E1AF/PEA3 genes. Moreover, as in human ERM and human E1AF/PEA3 genes, a first untranslated exon is upstream from the first methionine, and the mouse ER81 gene transcription is regulated by a 1.8 kbp of genomic DNA upstream from this exon. In human, the alternative splicing of the ETV1 gene leads to the presence (ETV1α) or the absence (ETV1β) of exon 5 encoding the C-terminal part of the transcriptional acidic domain, but without affecting the alpha helix previously described as crucial for transactivation. We demonstrated here that the truncated isoform (human ETV1β) and the full-length isoform (human ETV1α) bind similarly specific DNA Ets binding sites. Moreover, they both activate transcription similarly through the PKA-transduction pathway, so suggesting that this alternative splicing is not crucial for the function of this protein as a transcription factor. The comparison of human ETV1α and human ETV1β expression in the same tissues, such as the adrenal gland or the bladder, showed no clear-cut differences. Altogether, these data open a new avenue of investigation leading to a better understanding of the functional role of this transcription factor.