RNA-Seq Experiment and Data Analysis.

With the ability to obtain several millions of reads per sample, high-throughput RNA sequencing (RNA-Seq) enables investigation of any transcriptome at a fine resolution. Not just the messenger RNA (mRNA), but a wide variety of different RNA populations (e.g., total RNA, microRNA, long ncRNA, pre-mRNA) can also be investigated using RNA-Seq. While facilitating accurate quantification of gene expression, RNA-Seq offers the opportunity to estimate abundance of isoforms and find novel transcripts and allele-specific transcripts. In this chapter, we describe a protocol to construct an RNA-Seq library for sequencing on Illumina NGS platforms and a computational pipeline to perform RNA-Seq data analysis. The protocols described in this chapter can be applied to the analysis of differential gene expression in control versus 17ß-estradiol treatment of in vivo or in vitro systems.

A novel homozygous missense mutation p.P388S in TULP1 causes protein instability and retinitis pigmentosa.

Purpose: Retinitis pigmentosa (RP) is an inherited retinal disorder that results in the degeneration of photoreceptor cells, ultimately leading to severe visual impairment. We characterized a consanguineous family from Southern India wherein a 25 year old individual presented with night blindness since childhood. The purpose of this study was to identify the causative mutation for RP in this individual as well as characterize how the mutation may ultimately affect protein function. Methods: We performed a complete ophthalmologic examination of the proband followed by exome sequencing. The likely causative mutation was identified and modeled in cultured cells, evaluating its expression, solubility (both with western blotting), subcellular distribution, (confocal microscopy), and testing whether this variant induced endoplasmic reticulum (ER) stress (quantitative PCR [qPCR] and western blotting). Results: The proband presented with generalized and parafoveal retinal pigmented epithelium (RPE) atrophy with bone spicule-like pigmentation in the midperiphery and arteriolar attenuation. Optical coherence tomography scans through the macula of both eyes showed atrophy of the outer retinal layers with loss of the ellipsoid zone, whereas the systemic examination of this individual was normal. The proband's parents and sibling were asymptomatic and had normal funduscopic examinations. We discovered a novel homozygous p.Pro388Ser mutation in the tubby-like protein 1 (TULP1) gene in the individual with RP. In cultured cells, the P388S mutation does not alter the subcellular distribution of TULP1 or induce ER stress when compared to wild-type TULP1, but instead significantly lowers protein stability as indicated with steady-state and cycloheximide-chase experiments. Conclusions: These results add to the list of known mutations in TULP1 identified in individuals with RP and suggest a possible unique pathogenic mechanism in TULP1-induced RP, which may be shared among select mutations in TULP1.

Excessive deubiquitination of NLRP3-R779C variant contributes to very-early-onset inflammatory bowel disease development.

BACKGROUND: Very-early-onset inflammatory bowel disease (VEOIBD) is a chronic inflammatory disease of the gastrointestinal tract occurring during infancy or early childhood. NOD-like receptor family, pyrin domain containing 3 (NLRP3) inflammasome has emerged as a crucial regulator of intestinal homeostasis; however, whether NLRP3 variants may modify VEOIBD risk is unknown. OBJECTIVE: We sought to investigate whether and how a rare NLRP3 variant, found in 3 patients with gastrointestinal symptoms, contributes to VEOIBD development. METHODS: Whole-exome sequencing and bioinformatic analysis were performed to screen disease-associated NLRP3 variants from a cohort of children with VEOIBD. Inflammasome activation was determined in reconstituted HEK293T human embryonic kidney cells with NLRP3 inflammasome components, doxycycline-inducible NLRP3 macrophages, as well as PBMCs and biopsies from patients with NLRP3 variants. Pathogenesis of the variants was determined using a dextran sulfate sodium-induced acute colitis model. RESULTS: We identified a dominant gain-of-function missense variant of NLRP3, encoded by rs772009059 (R779C), in 3 patients with gastrointestinal symptoms. Functional analysis revealed that R779C increased NLRP3 inflammasome activation and pyroptosis in macrophages. This was mediated by enhanced deubiquitination of NLRP3 via binding with deubiquitinases BRCC3 and JOSD2, which are highly expressed in myeloid cells. In a dextran sulfate sodium-induced acute colitis model, NLRP3-R779C in hematopoietic cells resulted in more severe colitis, which can be ameliorated via knockdown of BRCC3 or JOSD2. CONCLUSIONS: BRCC3 and JOSD2 mediate NLRP3-R779C deubiquitination, which promotes NLRP3 inflammasome activation and the risk of developing VEOIBD.

Analyzing pre-symptomatic tissue to gain insights into the molecular and mechanistic origins of late-onset degenerative trinucleotide repeat disease.

How genetic defects trigger the molecular changes that cause late-onset disease is important for understanding disease progression and therapeutic development. Fuchs' endothelial corneal dystrophy (FECD) is an RNA-mediated disease caused by a trinucleotide CTG expansion in an intron within the TCF4 gene. The mutant intronic CUG RNA is present at one-two copies per cell, posing a challenge to understand how a rare RNA can cause disease. Late-onset FECD is a uniquely advantageous model for studying how RNA triggers disease because: (i) Affected tissue is routinely removed during surgery; (ii) The expanded CTG mutation is one of the most prevalent disease-causing mutations, making it possible to obtain pre-symptomatic tissue from eye bank donors to probe how gene expression changes precede disease; and (iii) The affected tissue is a homogeneous single cell monolayer, facilitating accurate transcriptome analysis. Here, we use RNA sequencing (RNAseq) to compare tissue from individuals who are pre-symptomatic (Pre_S) to tissue from patients with late stage FECD (FECD_REP). The abundance of mutant repeat intronic RNA in Pre_S and FECD_REP tissue is elevated due to increased half-life in a corneal cells. In Pre_S tissue, changes in splicing and extracellular matrix gene expression foreshadow the changes observed in advanced disease and predict the activation of the fibrosis pathway and immune system seen in late-stage patients. The absolute magnitude of splicing changes is similar in pre-symptomatic and late stage tissue. Our data identify gene candidates for early drivers of disease and biomarkers that may represent diagnostic and therapeutic targets for FECD. We conclude that changes in alternative splicing and gene expression are observable decades prior to the diagnosis of late-onset trinucleotide repeat disease.

Cord Blood CD8+ T Cells Have a Natural Propensity to Express IL-4 in a Fatty Acid Metabolism and Caspase Activation-Dependent Manner.

How T cells differentiate in the neonate may critically determine the ability of the infant to cope with infections, respond to vaccines and avert allergies. Previously, we found that naïve cord blood CD4+ T cells differentiated toward an IL-4-expressing phenotype when activated in the presence of TGF-ß and monocyte-derived inflammatory cytokines, the latter are more highly secreted by infants who developed food allergy. Here, we show that in the absence of IL-2 or IL-12, naïve cord blood CD8+ T cells have a natural propensity to differentiate into IL-4-producing non-classic TC2 cells when they are activated alone, or in the presence of TGF-ß and/or inflammatory cytokines. Mechanistically, non-classic TC2 development is associated with decreased expression of IL-2 receptor alpha (CD25) and glycolysis, and increased fatty acid metabolism and caspase-dependent cell death. Consequently, the short chain fatty acid, sodium propionate (NaPo), enhanced IL-4 expression, but exogenous IL-2 or pan-caspase inhibition prevented IL-4 expression. In children with endoscopically and histologically confirmed non-inflammatory bowel disease and non-infectious pediatric idiopathic colitis, the presence of TGF-ß, NaPo, and IL-1ß or TNF-α promoted TC2 differentiation in vitro. In vivo, colonic mucosa of children with colitis had significantly increased expression of IL-4 in CD8+ T cells compared with controls. In addition, activated caspase-3 and IL-4 were co-expressed in CD8+ T cells in the colonic mucosa of children with colitis. Thus, in the context of colonic inflammation and limited IL-2 signaling, CD8+ T cells differentiate into non-classic TC2 that may contribute to the pathology of inflammatory/allergic diseases in children.

The IFN Response in Bats Displays Distinctive IFN-Stimulated Gene Expression Kinetics with Atypical RNASEL Induction.

Bats host a large number of zoonotic viruses, including several viruses that are highly pathogenic to other mammals. The mechanisms underlying this rich viral diversity are unknown, but they may be linked to unique immunological features that allow bats to act as asymptomatic viral reservoirs. Vertebrates respond to viral infection by inducing IFNs, which trigger antiviral defenses through IFN-stimulated gene (ISG) expression. Although the IFN system of several bats is characterized at the genomic level, less is known about bat IFN-mediated transcriptional responses. In this article, we show that IFN signaling in bat cells from the black flying fox (Pteropus alecto) consists of conserved and unique ISG expression profiles. In IFN-stimulated cells, bat ISGs comprise two unique temporal subclusters with similar early induction kinetics but distinct late-phase declines. In contrast, human ISGs lack this decline phase and remained elevated for longer periods. Notably, in unstimulated cells, bat ISGs were expressed more highly than their human counterparts. We also found that the antiviral effector 2-5A-dependent endoribonuclease, which is not an ISG in humans, is highly IFN inducible in black flying fox cells and contributes to cell-intrinsic control of viral infection. These studies reveal distinctive innate immune features that may underlie a unique virus-host relationship in bats.

A comprehensively characterized cell line panel highly representative of clinical ovarian high-grade serous carcinomas.

Recent literature suggests that most widely used ovarian cancer (OVCA) cell models do not recapitulate the molecular features of clinical tumors. To address this limitation, we generated 18 cell lines and 3 corresponding patient-derived xenografts predominantly from high-grade serous carcinoma (HGSOC) peritoneal effusions. Comprehensive genomic characterization and comparison of each model to its parental tumor demonstrated a high degree of molecular similarity. Our characterization included whole exome-sequencing and copy number profiling for cell lines, xenografts, and matched non-malignant tissues, and DNA methylation, gene expression, and spectral karyotyping for a subset of specimens. Compared to the Cancer Genome Atlas (TCGA), our models more closely resembled HGSOC than any other tumor type, justifying their validity as OVCA models. Our meticulously characterized models provide a crucial resource for the OVCA research community that will advance translational findings and ultimately lead to clinical applications.

Prostaglandin dehydrogenase is a target for successful induction of cervical ripening.

The cervix represents a formidable structural barrier for successful induction of labor. Approximately 10% of pregnancies undergo induction of cervical ripening and labor with prostaglandin (PG) E2 or PGE analogs, often requiring many hours of hospitalization and monitoring. On the other hand, preterm cervical ripening in the second trimester predicts preterm birth. The regulatory mechanisms of this paradoxical function of the cervix are unknown. Here, we show that PGE2 uses cell-specific EP2 receptor-mediated increases in Ca2+ to dephosphorylate and translocate histone deacetylase 4 (HDAC4) to the nucleus for repression of 15-hydroxy prostaglandin dehydrogenase (15-PGDH). The crucial role of 15-PGDH in cervical ripening was confirmed in vivo. Although PGE2 or 15-PGDH inhibitor alone did not alter gestational length, treatment with 15-PGDH inhibitor + PGE2 or metabolism-resistant dimethyl-PGE2 resulted in preterm cervical ripening and delivery in mice. The ability of PGE2 to selectively autoamplify its own synthesis in stromal cells by signaling transcriptional repression of 15-PGDH elucidates long sought-after molecular mechanisms that govern PG action in the cervix. This report details unique mechanisms of action in the cervix and serves as a catalyst for (i) the use of 15-PGDH inhibitors to initiate or amplify low-dose PGE2-mediated cervical ripening or (ii) EP2 receptor antagonists, HDAC4 inhibitors, and 15-PGDH activators to prevent preterm cervical ripening and preterm birth.

Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp. tritici isolates revealed by the comparative gene co-expression network and genome analyses.

BACKGROUND: Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt). RESULTS: The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub-network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors. CONCLUSIONS: Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta.

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Suppression of the SWI/SNF Component Arid1a Promotes Mammalian Regeneration.

Mammals have partially lost the extensive regenerative capabilities of some vertebrates, possibly as a result of chromatin-remodeling mechanisms that enforce terminal differentiation. Here, we show that deleting the SWI/SNF component Arid1a substantially improves mammalian regeneration. Arid1a expression is suppressed in regenerating tissues, and genetic deletion of Arid1a increases tissue repair following an array of injuries. Arid1a deficiency in the liver increases proliferation, reduces tissue damage and fibrosis, and improves organ function following surgical resection and chemical injuries. Hepatocyte-specific deletion is also sufficient to increase proliferation and regeneration without excessive overgrowth, and global Arid1a disruption potentiates soft tissue healing in the ear. We show that Arid1a loss reprograms chromatin to restrict promoter access by transcription factors such as C/ebpα, which enforces differentiation, and E2F4, which suppresses cell-cycle re-entry. Thus, epigenetic reprogramming mediated by deletion of a single gene improves mammalian regeneration and suggests strategies to promote tissue repair after injury.

RNA-Seq Experiment and Data Analysis.

With the ability to obtain tens of millions of reads, high-throughput messenger RNA sequencing (RNA-Seq) data offers the possibility of estimating abundance of isoforms and finding novel transcripts. In this chapter, we describe a protocol to construct an RNA-Seq library for sequencing on Illumina NGS platforms, and a computational pipeline to perform RNA-Seq data analysis. The protocols described in this chapter can be applied to the analysis of differential gene expression in control versus 17ß-estradiol treatment of in vivo or in vitro systems.

Hypoxia fate mapping identifies cycling cardiomyocytes in the adult heart.

Although the adult mammalian heart is incapable of meaningful functional recovery following substantial cardiomyocyte loss, it is now clear that modest cardiomyocyte turnover occurs in adult mouse and human hearts, mediated primarily by proliferation of pre-existing cardiomyocytes. However, fate mapping of these cycling cardiomyocytes has not been possible thus far owing to the lack of identifiable genetic markers. In several organs, stem or progenitor cells reside in relatively hypoxic microenvironments where the stabilization of the hypoxia-inducible factor 1 alpha (Hif-1α) subunit is critical for their maintenance and function. Here we report fate mapping of hypoxic cells and their progenies by generating a transgenic mouse expressing a chimaeric protein in which the oxygen-dependent degradation (ODD) domain of Hif-1α is fused to the tamoxifen-inducible CreERT2 recombinase. In mice bearing the creERT2-ODD transgene driven by either the ubiquitous CAG promoter or the cardiomyocyte-specific α myosin heavy chain promoter, we identify a rare population of hypoxic cardiomyocytes that display characteristics of proliferative neonatal cardiomyocytes, such as smaller size, mononucleation and lower oxidative DNA damage. Notably, these hypoxic cardiomyocytes contributed widely to new cardiomyocyte formation in the adult heart. These results indicate that hypoxia signalling is an important hallmark of cycling cardiomyocytes, and suggest that hypoxia fate mapping can be a powerful tool for identifying cycling cells in adult mammals.

Comparative analysis of syntenic genes in grass genomes reveals accelerated rates of gene structure and coding sequence evolution in polyploid wheat.

Cycles of whole-genome duplication (WGD) and diploidization are hallmarks of eukaryotic genome evolution and speciation. Polyploid wheat (Triticum aestivum) has had a massive increase in genome size largely due to recent WGDs. How these processes may impact the dynamics of gene evolution was studied by comparing the patterns of gene structure changes, alternative splicing (AS), and codon substitution rates among wheat and model grass genomes. In orthologous gene sets, significantly more acquired and lost exonic sequences were detected in wheat than in model grasses. In wheat, 35% of these gene structure rearrangements resulted in frame-shift mutations and premature termination codons. An increased codon mutation rate in the wheat lineage compared with Brachypodium distachyon was found for 17% of orthologs. The discovery of premature termination codons in 38% of expressed genes was consistent with ongoing pseudogenization of the wheat genome. The rates of AS within the individual wheat subgenomes (21%-25%) were similar to diploid plants. However, we uncovered a high level of AS pattern divergence between the duplicated homeologous copies of genes. Our results are consistent with the accelerated accumulation of AS isoforms, nonsynonymous mutations, and gene structure rearrangements in the wheat lineage, likely due to genetic redundancy created by WGDs. Whereas these processes mostly contribute to the degeneration of a duplicated genome and its diploidization, they have the potential to facilitate the origin of new functional variations, which, upon selection in the evolutionary lineage, may play an important role in the origin of novel traits.

H3K9 methylation is a barrier during somatic cell reprogramming into iPSCs.

The induction of pluripotent stem cells (iPSCs) by defined factors is poorly understood stepwise. Here, we show that histone H3 lysine 9 (H3K9) methylation is the primary epigenetic determinant for the intermediate pre-iPSC state, and its removal leads to fully reprogrammed iPSCs. We generated a panel of stable pre-iPSCs that exhibit pluripotent properties but do not activate the core pluripotency network, although they remain sensitive to vitamin C for conversion into iPSCs. Bone morphogenetic proteins (BMPs) were subsequently identified in serum as critical signaling molecules in arresting reprogramming at the pre-iPSC state. Mechanistically, we identified H3K9 methyltransferases as downstream targets of BMPs and showed that they function with their corresponding demethylases as the on/off switch for the pre-iPSC fate by regulating H3K9 methylation status at the core pluripotency loci. Our results not only establish pre-iPSCs as an epigenetically stable signpost along the reprogramming road map, but they also provide mechanistic insights into the epigenetic reprogramming of cell fate.

Homoeolog-specific transcriptional bias in allopolyploid wheat.

BACKGROUND: Interaction between parental genomes is accompanied by global changes in gene expression which, eventually, contributes to growth vigor and the broader phenotypic diversity of allopolyploid species. In order to gain a better understanding of the effects of allopolyploidization on the regulation of diverged gene networks, we performed a genome-wide analysis of homoeolog-specific gene expression in re-synthesized allohexaploid wheat created by the hybridization of a tetraploid derivative of hexaploid wheat with the diploid ancestor of the wheat D genome Ae. tauschii. RESULTS: Affymetrix wheat genome arrays were used for both the discovery of divergent homoeolog-specific mutations and analysis of homoeolog-specific gene expression in re-synthesized allohexaploid wheat. More than 34,000 detectable parent-specific features (PSF) distributed across the wheat genome were used to assess AB genome (could not differentiate A and B genome contributions) and D genome parental expression in the allopolyploid transcriptome. In re-synthesized polyploid 81% of PSFs detected mid-parent levels of gene expression, and only 19% of PSFs showed the evidence of non-additive expression. Non-additive expression in both AB and D genomes was strongly biased toward up-regulation of parental type of gene expression with only 6% and 11% of genes, respectively, being down-regulated. Of all the non-additive gene expression, 84% can be explained by differences in the parental genotypes used to make the allopolyploid. Homoeolog-specific co-regulation of several functional gene categories was found, particularly genes involved in photosynthesis and protein biosynthesis in wheat. CONCLUSIONS: Here, we have demonstrated that the establishment of interactions between the diverged regulatory networks in allopolyploids is accompanied by massive homoeolog-specific up- and down-regulation of gene expression. This study provides insights into interactions between homoeologous genomes and their role in growth vigor, development, and fertility of allopolyploid species.

Autographa californica multiple nucleopolyhedrovirus 38K is a novel nucleocapsid protein that interacts with VP1054, VP39, VP80, and itself.

It has been shown that the Autographa californica multiple nucleopolyhedrovirus (AcMNPV) 38K (ac98) is required for nucleocapsid assembly. However, the exact role of 38K in nucleocapsid assembly remains unknown. In the present study, we investigated the relationship between 38K and the nucleocapsid. Western blotting using polyclonal antibodies raised against 38K revealed that 38K was expressed in the late phase of infection in AcMNPV-infected Spodoptera frugiperda cells and copurified with budded virus (BV) and occlusion-derived virus (ODV). Biochemical fractionation of BV and ODV into the nucleocapsid and envelope components followed by Western blotting showed that 38K was associated with the nucleocapsids. Immunoelectron microscopic analysis revealed that 38K was specifically localized to the nucleocapsids in infected cells and appeared to be distributed over the cylindrical capsid sheath of nucleocapsid. Yeast two-hybrid assays were performed to examine potential interactions between 38K and nine known nucleocapsid shell-associated proteins (PP78/83, PCNA, VP1054, FP25, VLF-1, VP39, BV/ODV-C42, VP80, and P24), three non-nucleocapsid shell-associated proteins (P6.9, PP31, and BV/ODV-E26), and itself. The results revealed that 38K interacted with the nucleocapsid proteins VP1054, VP39, VP80, and 38K itself. These interactions were confirmed by coimmunoprecipitation assays in vivo. These data demonstrate that 38K is a novel nucleocapsid protein and provide a rationale for why 38K is essential for nucleocapsid assembly.