Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

Altered gene expression and DNA damage in peripheral blood cells from Friedreich's ataxia patients: cellular model of pathology.

The neurodegenerative disease Friedreich's ataxia (FRDA) is the most common autosomal-recessively inherited ataxia and is caused by a GAA triplet repeat expansion in the first intron of the frataxin gene. In this disease, transcription of frataxin, a mitochondrial protein involved in iron homeostasis, is impaired, resulting in a significant reduction in mRNA and protein levels. Global gene expression analysis was performed in peripheral blood samples from FRDA patients as compared to controls, which suggested altered expression patterns pertaining to genotoxic stress. We then confirmed the presence of genotoxic DNA damage by using a gene-specific quantitative PCR assay and discovered an increase in both mitochondrial and nuclear DNA damage in the blood of these patients (p<0.0001, respectively). Additionally, frataxin mRNA levels correlated with age of onset of disease and displayed unique sets of gene alterations involved in immune response, oxidative phosphorylation, and protein synthesis. Many of the key pathways observed by transcription profiling were downregulated, and we believe these data suggest that patients with prolonged frataxin deficiency undergo a systemic survival response to chronic genotoxic stress and consequent DNA damage detectable in blood. In conclusion, our results yield insight into the nature and progression of FRDA, as well as possible therapeutic approaches. Furthermore, the identification of potential biomarkers, including the DNA damage found in peripheral blood, may have predictive value in future clinical trials.

Diagnosis of partial body radiation exposure in mice using peripheral blood gene expression profiles.

In the event of a terrorist-mediated attack in the United States using radiological or improvised nuclear weapons, it is expected that hundreds of thousands of people could be exposed to life-threatening levels of ionizing radiation. We have recently shown that genome-wide expression analysis of the peripheral blood (PB) can generate gene expression profiles that can predict radiation exposure and distinguish the dose level of exposure following total body irradiation (TBI). However, in the event a radiation-mass casualty scenario, many victims will have heterogeneous exposure due to partial shielding and it is unknown whether PB gene expression profiles would be useful in predicting the status of partially irradiated individuals. Here, we identified gene expression profiles in the PB that were characteristic of anterior hemibody-, posterior hemibody- and single limb-irradiation at 0.5 Gy, 2 Gy and 10 Gy in C57Bl6 mice. These PB signatures predicted the radiation status of partially irradiated mice with a high level of accuracy (range 79-100%) compared to non-irradiated mice. Interestingly, PB signatures of partial body irradiation were poorly predictive of radiation status by site of injury (range 16-43%), suggesting that the PB molecular response to partial body irradiation was anatomic site specific. Importantly, PB gene signatures generated from TBI-treated mice failed completely to predict the radiation status of partially irradiated animals or non-irradiated controls. These data demonstrate that partial body irradiation, even to a single limb, generates a characteristic PB signature of radiation injury and thus may necessitate the use of multiple signatures, both partial body and total body, to accurately assess the status of an individual exposed to radiation.

Pax3 expression enhances PDGF-B-induced brainstem gliomagenesis and characterizes a subset of brainstem glioma.

High-grade Brainstem Glioma (BSG), also known as Diffuse Intrinsic Pontine Glioma (DIPG), is an incurable pediatric brain cancer. Increasing evidence supports the existence of regional differences in gliomagenesis such that BSG is considered a distinct disease from glioma of the cerebral cortex (CG). In an effort to elucidate unique characteristics of BSG, we conducted expression analysis of mouse PDGF-B-driven BSG and CG initiated in Nestin progenitor cells and identified a short list of expression changes specific to the brainstem gliomagenesis process, including abnormal upregulation of paired box 3 (Pax3). In the neonatal mouse brain, Pax3 expression marks a subset of brainstem progenitor cells, while it is absent from the cerebral cortex, mirroring its regional expression in glioma. Ectopic expression of Pax3 in normal brainstem progenitors in vitro shows that Pax3 inhibits apoptosis. Pax3-induced inhibition of apoptosis is p53-dependent, however, and in the absence of p53, Pax3 promotes proliferation of brainstem progenitors. In vivo, Pax3 enhances PDGF-B-driven gliomagenesis by shortening tumor latency and increasing tumor penetrance and grade, in a region-specific manner, while loss of Pax3 function extends survival of PDGF-B-driven;p53-deficient BSG-bearing mice by 33%. Importantly, Pax3 is regionally expressed in human glioma as well, with high PAX3 mRNA characterizing 40% of human BSG, revealing a subset of tumors that significantly associates with PDGFRA alterations, amplifications of cell cycle regulatory genes, and is exclusive of ACVR1 mutations. Collectively, these data suggest that regional Pax3 expression not only marks a novel subset of BSG but also contributes to PDGF-B-induced brainstem gliomagenesis.

Transcriptomic Analysis of the Host Response and Innate Resilience to Enterotoxigenic Escherichia coli Infection in Humans.

BACKGROUND: Enterotoxigenic Escherichia coli (ETEC) is a globally prevalent cause of diarrhea. Though usually self-limited, it can be severe and debilitating. Little is known about the host transcriptional response to infection. We report the first gene expression analysis of the human host response to experimental challenge with ETEC. METHODS: We challenged 30 healthy adults with an unattenuated ETEC strain, and collected serial blood samples shortly after inoculation and daily for 8 days. We performed gene expression analysis on whole peripheral blood RNA samples from subjects in whom severe symptoms developed (n = 6) and a subset of those who remained asymptomatic (n = 6) despite shedding. RESULTS: Compared with baseline, symptomatic subjects demonstrated significantly different expression of 406 genes highlighting increased immune response and decreased protein synthesis. Compared with asymptomatic subjects, symptomatic subjects differentially expressed 254 genes primarily associated with immune response. This comparison also revealed 29 genes differentially expressed between groups at baseline, suggesting innate resilience to infection. Drug repositioning analysis identified several drug classes with potential utility in augmenting immune response or mitigating symptoms. CONCLUSIONS: There are statistically significant and biologically plausible differences in host gene expression induced by ETEC infection. Differential baseline expression of some genes may indicate resilience to infection.

Altered ultrasonic vocalization and impaired learning and memory in Angelman syndrome mouse model with a large maternal deletion from Ube3a to Gabrb3.

Angelman syndrome (AS) is a neurobehavioral disorder associated with mental retardation, absence of language development, characteristic electroencephalography (EEG) abnormalities and epilepsy, happy disposition, movement or balance disorders, and autistic behaviors. The molecular defects underlying AS are heterogeneous, including large maternal deletions of chromosome 15q11-q13 (70%), paternal uniparental disomy (UPD) of chromosome 15 (5%), imprinting mutations (rare), and mutations in the E6-AP ubiquitin ligase gene UBE3A (15%). Although patients with UBE3A mutations have a wide spectrum of neurological phenotypes, their features are usually milder than AS patients with deletions of 15q11-q13. Using a chromosomal engineering strategy, we generated mutant mice with a 1.6-Mb chromosomal deletion from Ube3a to Gabrb3, which inactivated the Ube3a and Gabrb3 genes and deleted the Atp10a gene. Homozygous deletion mutant mice died in the perinatal period due to a cleft palate resulting from the null mutation in Gabrb3 gene. Mice with a maternal deletion (m-/p+) were viable and did not have any obvious developmental defects. Expression analysis of the maternal and paternal deletion mice confirmed that the Ube3a gene is maternally expressed in brain, and showed that the Atp10a and Gabrb3 genes are biallelically expressed in all brain sub-regions studied. Maternal (m-/p+), but not paternal (m+/p-), deletion mice had increased spontaneous seizure activity and abnormal EEG. Extensive behavioral analyses revealed significant impairment in motor function, learning and memory tasks, and anxiety-related measures assayed in the light-dark box in maternal deletion but not paternal deletion mice. Ultrasonic vocalization (USV) recording in newborns revealed that maternal deletion pups emitted significantly more USVs than wild-type littermates. The increased USV in maternal deletion mice suggests abnormal signaling behavior between mothers and pups that may reflect abnormal communication behaviors in human AS patients. Thus, mutant mice with a maternal deletion from Ube3a to Gabrb3 provide an AS mouse model that is molecularly more similar to the contiguous gene deletion form of AS in humans than mice with Ube3a mutation alone. These mice will be valuable for future comparative studies to mice with maternal deficiency of Ube3a alone.

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.

Engineering Highly-functional, Self-regenerative Skeletal Muscle Tissues with Enhanced Vascularization and Survival in Vivo

Tissue engineering of biomimetic skeletal muscle may lead to development of new therapies for myogenic repair and generation of improved in vitro models for studies of muscle function, regeneration, and disease. For the optimal therapeutic and in vitro results, engineered muscle should recreate the force-generating and regenerative capacities of native muscle, enabled respectively by its two main cellular constituents, the mature myofibers and satellite cells (SCs). Still, after 20 years of research, engineered muscle tissues fall short of mimicking contractile function and self-repair capacity of native skeletal muscle. To overcome this limitation, we set the thesis goals to: 1) generate a highly functional, self-regenerative engineered skeletal muscle and 2) explore mechanisms governing its formation and regeneration in vitro and survival and vascularization in vivo.

By studying myogenic progenitors isolated from neonatal rats, we first discovered advantages of using an adherent cell fraction for engineering of skeletal muscles with robust structure and function and the formation of a SC pool. Specifically, when synergized with dynamic culture conditions, the use of adherent cells yielded muscle constructs capable of replicating the contractile output of native neonatal muscle, generating >40 mN/mm2 of specific force. Moreover, tissue structure and cellular heterogeneity of engineered muscle constructs closely resembled those of native muscle, consisting of aligned, striated myofibers embedded in a matrix of basal lamina proteins and SCs that resided in native-like niches. Importantly, we identified rapid formation of myofibers early during engineered muscle culture as a critical condition leading to SC homing and conversion to a quiescent, non-proliferative state. The SCs retained natural regenerative capacity and activated, proliferated, and differentiated to rebuild damaged myofibers and recover contractile function within 10 days after the muscle was injured by cardiotoxin (CTX). The resulting regenerative response was directly dependent on the abundance of SCs in the engineered muscle that we varied by expanding starting cell population under different levels of basic fibroblast growth factor (bFGF), an inhibitor of myogenic differentiation. Using a dorsal skinfold window chamber model in nude mice, we further demonstrated that within 2 weeks after implantation, initially avascular engineered muscle underwent robust vascularization and perfusion and exhibited improved structure and contractile function beyond what was achievable in vitro.

To enhance translational value of our approach, we transitioned to use of adult rat myogenic cells, but found that despite similar function to that of neonatal constructs, adult-derived muscle lacked regenerative capacity. Using a novel platform for live monitoring of calcium transients during construct culture, we rapidly screened for potential enhancers of regeneration to establish that many known pro-regenerative soluble factors were ineffective in stimulating in vitro engineered muscle recovery from CTX injury. This led us to introduce bone marrow-derived macrophages (BMDMs), an established non-myogenic contributor to muscle repair, to the adult-derived constructs and to demonstrate remarkable recovery of force generation (>80%) and muscle mass (>70%) following CTX injury. Mechanistically, while similar patterns of early SC activation and proliferation upon injury were observed in engineered muscles with and without BMDMs, a significant decrease in injury-induced apoptosis occurred only in the presence of BMDMs. The importance of preventing apoptosis was further demonstrated by showing that application of caspase inhibitor (Q-VD-OPh) yielded myofiber regrowth and functional recovery post-injury. Gene expression analysis suggested muscle-secreted tumor necrosis factor-α (TNFα) as a potential inducer of apoptosis as common for muscle degeneration in diseases and aging in vivo. Finally, we showed that BMDM incorporation in engineered muscle enhanced its growth, angiogenesis, and function following implantation in the dorsal window chambers in nude mice.

In summary, this thesis describes novel strategies to engineer highly contractile and regenerative skeletal muscle tissues starting from neonatal or adult rat myogenic cells. We find that age-dependent differences of myogenic cells distinctly affect the self-repair capacity but not contractile function of engineered muscle. Adult, but not neonatal, myogenic progenitors appear to require co-culture with other cells, such as bone marrow-derived macrophages, to allow robust muscle regeneration in vitro and rapid vascularization in vivo. Regarding the established roles of immune system cells in the repair of various muscle and non-muscle tissues, we expect that our work will stimulate the future applications of immune cells as pro-regenerative or anti-inflammatory constituents of engineered tissue grafts. Furthermore, we expect that rodent studies in this thesis will inspire successful engineering of biomimetic human muscle tissues for use in regenerative therapy and drug discovery applications.

The genomic analysis of erythrocyte microRNA expression in sickle cell diseases.

BACKGROUND: Since mature erythrocytes are terminally differentiated cells without nuclei and organelles, it is commonly thought that they do not contain nucleic acids. In this study, we have re-examined this issue by analyzing the transcriptome of a purified population of human mature erythrocytes from individuals with normal hemoglobin (HbAA) and homozygous sickle cell disease (HbSS). METHODS AND FINDINGS: Using a combination of microarray analysis, real-time RT-PCR and Northern blots, we found that mature erythrocytes, while lacking ribosomal and large-sized RNAs, contain abundant and diverse microRNAs. MicroRNA expression of erythrocytes was different from that of reticulocytes and leukocytes, and contributed the majority of the microRNA expression in whole blood. When we used microRNA microarrays to analyze erythrocytes from HbAA and HbSS individuals, we noted a dramatic difference in their microRNA expression pattern. We found that miR-320 played an important role for the down-regulation of its target gene, CD71 during reticulocyte terminal differentiation. Further investigation revealed that poor expression of miR-320 in HbSS cells was associated with their defective downregulation CD71 during terminal differentiation. CONCLUSIONS: In summary, we have discovered significant microRNA expression in human mature erythrocytes, which is dramatically altered in HbSS erythrocytes and their defect in terminal differentiation. Thus, the global analysis of microRNA expression in circulating erythrocytes can provide mechanistic insights into the disease phenotypes of erythrocyte diseases.

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.