Dynamics of dendritic cell maturation are identified through a novel filtering strategy applied to biological time-course microarray replicates

Long tracks of homozygosity predict the severity of alcohol use disorders in an American Indian population

Runs of homozygosity (ROH) arise when an individual inherits two copies of the same haplotype segment. While ROH are ubiquitous across human populations, Native populations-with shared parental ancestry arising from isolation and endogamy-can carry a substantial enrichment for ROH. We have been investigating genetic and environmental risk factors for alcohol use disorders (AUD) in a group of American Indians (AI) who have higher rates of AUD than the general U. S. population. Here we explore whether ROH might be associated with incidence and severity of AUD in this admixed AI population (n = 742) that live on geographically contiguous reservations, using low-coverage whole genome sequences. We have found that the genomic regions in the ROH that were identified in this population had significantly elevated American Indian heritage compared with the rest of the genome. Increased ROH abundance and ROH burden are likely risk factors for AUD severity in this AI population, especially in those diagnosed with severe and moderate AUD. The association between ROH and AUD was mostly driven by ROH of moderate lengths between 1 and 2 Mb. An ROH island on chromosome 1p32.3 and a rare ROH pool on chromosome 3p12.3 were found to be significantly associated with AUD severity. They contain genes involved in lipid metabolism, oxidative stress and inflammatory responses; and OSBPL9 was found to reside on the consensus part of the ROH island. These data demonstrate that ROH are associated with risk for AUD severity in this AI population.

An in silico model of cytotoxic T-lymphocyte activation in the lymph node following short peptide vaccination

Tumour immunotherapy is dependent upon activation and expansion of tumour-targetting immune cells, known as cytotoxic T-lymphocytes (CTLs). Cancer vaccines developed in the past have had limited success and the mechanisms resulting in failure are not well characterized. To elucidate these mechanisms, we developed a human-parametrized, in silico, agent-based model of vaccination-driven CTL activation within a clinical short-peptide vaccination context. The simulations predict a sharp transition in the probability of CTL activation, which occurs with variation in the separation rate (or off-rate) of tumour-specific immune response-inducing peptides (cognate antigen) from the major histocompatibility class I (MHC-I) receptors of dendritic cells (DCs) originally at the vaccination site. For peptides with MHC-I off-rates beyond this transition, it is predicted that no vaccination strategy will lead to successful expansion of CTLs. For slower off-rates, below the transition, the probability of CTL activation becomes sensitive to the numbers of DCs and T cells that interact subsequent to DC migration to the draining lymph node of the vaccination site. Thus, the off-rate is a key determinant of vaccine design.

Quantitative proteomic changes in LPS-activated monocyte-derived dendritic cells: A SWATH-MS study

Dendritic cells are key immune cells that respond to pathogens and co-ordinate many innate and adaptive immune responses. Quantitative mass spectrometry using Sequential Window Acquisition of all THeoretical fragment-ion spectra-Mass Spectrometry (SWATH-MS) was performed here to determine the global alterations in monocyte-derived dendritic cells (moDCs) in response to stimulation with lipopolysaccharide (LPS). A moDC library of 4,666 proteins was generated and proteins were quantified at 0, 6 and 24 h post-LPS stimulation using SWATH-MS. At 6 h and 24 h post-LPS exposure, the relative abundance of 227 and 282 proteins was statistically significantly altered (p-value ≤ 0.05), respectively. Functional annotation of proteins exhibiting significant changes in expression between the various time points led to the identification of clusters of proteins implicated in distinct cellular processes including interferon and interleukin signalling, endocytosis, the ER-phagosome pathway and antigen-presentation. In SWATH-MS major histocompatibility complex (MHC) class I proteins were highly upregulated at 24 h, whilst MHC class II proteins exhibited comparatively fewer changes over this period. This study provides new detailed insight into the global proteomic changes that occur in moDCs during antigen processing and presentation and further demonstrates the potential of SWATH-MS for the quantitative study of proteins involved in cellular processes.

Identification of Transcriptional and Receptor Networks That Control Root Responses to Ethylene

Transcriptomic analyses with high temporal resolution provide substantial new insight into hormonal response networks. This study identified the kinetics of genome-wide transcript abundance changes in response to elevated levels of the plant hormone ethylene in roots from light-grown Arabidopsis (Arabidopsis thaliana) seedlings, which were overlaid on time-matched developmental changes. Functional annotation of clusters of transcripts with similar temporal patterns revealed rapidly induced clusters with known ethylene function and more slowly regulated clusters with novel predicted functions linked to root development. In contrast to studies with dark-grown seedlings, where the canonical ethylene response transcription factor, EIN3, is central to ethylene-mediated development, the roots of ein3 and eil1 single and double mutants still respond to ethylene in light-grown seedlings. Additionally, a subset of these clusters of ethylene-responsive transcripts were enriched in targets of EIN3 and ERFs. These results are consistent with EIN3-independent developmental and transcriptional changes in light-grown roots. Examination of single and multiple gain-of-function and loss-of-function receptor mutants revealed that, of the five ethylene receptors, ETR1 controls lateral root and root hair initiation and elongation and the synthesis of other receptors. These results provide new insight into the transcriptional and developmental responses to ethylene in light-grown seedlings.

p21 Restricts HIV-1 in Monocyte-Derived Dendritic Cells through the Reduction of Deoxynucleoside Triphosphate Biosynthesis and Regulation of SAMHD1 Antiviral Activity

HIV-1 infection of noncycling cells, such as dendritic cells (DCs), is impaired due to limited availability of deoxynucleoside triphosphates (dNTPs), which are needed for HIV-1 reverse transcription. The levels of dNTPs are tightly regulated during the cell cycle and depend on the balance between dNTP biosynthesis and degradation. SAMHD1 potently blocks HIV-1 replication in DCs, although the underlying mechanism is still unclear. SAMHD1 has been reported to be able to degrade dNTPs and viral nucleic acids, which may both hamper HIV-1 reverse transcription. The relative contribution of these activities may differ in cycling and noncycling cells. Here, we show that inhibition of HIV-1 replication in monocyte-derived DCs (MDDCs) is associated with an increased expression of p21cip1/waf, a cell cycle regulator that is involved in the differentiation and maturation of DCs. Induction of p21 in MDDCs decreases the pool of dNTPs and increases the antiviral active isoform of SAMHD1. Although both processes are complementary in inhibiting HIV-1 replication, the antiviral activity of SAMHD1 in our primary cell model appears to be, at least partially, independent of its dNTPase activity. The reduction in the pool of dNTPs in MDDCs appears rather mostly due to a p21-mediated suppression of several enzymes involved in dNTP synthesis (i.e., RNR2, TYMS, and TK-1). These results are important to better understand the interplay between HIV-1 and DCs and may inform the design of new therapeutic approaches to decrease viral dissemination and improve immune responses against HIV-1.IMPORTANCE DCs play a key role in the induction of immune responses against HIV. However, HIV has evolved ways to exploit these cells, facilitating immune evasion and virus dissemination. We have found that the expression of p21, a cyclin-dependent kinase inhibitor involved in cell cycle regulation and monocyte differentiation and maturation, potentially can contribute to the inhibition of HIV-1 replication in monocyte-derived DCs through multiple mechanisms. p21 decreased the size of the intracellular dNTP pool. In parallel, p21 prevented SAMHD1 phosphorylation and promoted SAMHD1 dNTPase-independent antiviral activity. Thus, induction of p21 resulted in conditions that allowed the effective inhibition of HIV-1 replication through complementary mechanisms. Overall, p21 appears to be a key regulator of HIV infection in myeloid cells.

Impact of the Type I Interferon Receptor on the Global Gene Expression Program During the Course of Dendritic Cell Maturation Induced by Polyinosinic Polycytidylic Acid

Dendritic cell (DC) maturation involves widespread changes in cellular function and gene expression. The regulatory role of IFNAR in the program of DC maturation remains incompletely defined. Thus, the time evolution impact of IFNAR on this process was evaluated. Changes in DC phenotype, function, and gene expression induced by poly I:C were measured in wild-type and IFNAR(-/-) DC at 9 time points over 24 h. Temporal gene expression profiles were filtered on consistency and response magnitude across replicates. The number of genes whose expression was altered by poly I:C treatment was greatly reduced in IFNAR(-/-) DC, including the majority of the downregulated gene expression program previously observed in wild-type (WT) DC. Furthermore, the number of genes upregulated was almost equal between WT and IFNAR(-/-) DC, yet the identities of those genes were distinct. Integrating these data with protein-protein interaction data revealed several novel subnetworks active during maturation, including nucleotide synthesis, metabolism, and repair. A subnetwork associated with redox activity was uniquely identified in IFNAR(-/-) DC. Overall, temporal gene expression and network analyses identified many genes regulated by the type I interferon response and revealed previously unidentified aspects of the DC maturation process.

A kinetic analysis of the auxin transcriptome reveals cell wall remodeling proteins that modulate lateral root development in Arabidopsis

To identify gene products that participate in auxin-dependent lateral root formation, a high temporal resolution, genome-wide transcript abundance analysis was performed with auxin-treated Arabidopsis thaliana roots. Data analysis identified 1246 transcripts that were consistently regulated by indole-3-acetic acid (IAA), partitioning into 60 clusters with distinct response kinetics. We identified rapidly induced clusters containing auxin-response functional annotations and clusters exhibiting delayed induction linked to cell division temporally correlated with lateral root induction. Several clusters were enriched with genes encoding proteins involved in cell wall modification, opening the possibility for understanding mechanistic details of cell structural changes that result in root formation following auxin treatment. Mutants with insertions in 72 genes annotated with a cell wall remodeling function were examined for alterations in IAA-regulated root growth and development. This reverse-genetic screen yielded eight mutants with root phenotypes. Detailed characterization of seedlings with mutations in cellulase3/glycosylhydrolase9b3 and leucine rich extensin2, genes not normally linked to auxin response, revealed defects in the early and late stages of lateral root development, respectively. The genes identified here using kinetic insight into expression changes lay the foundation for mechanistic understanding of auxin-mediated cell wall remodeling as an essential feature of lateral root development.

Genome-wide analysis of alternative splicing during dendritic cell response to a bacterial challenge

The immune system relies on the plasticity of its components to produce appropriate responses to frequent environmental challenges. Dendritic cells (DCs) are critical initiators of innate immunity and orchestrate the later and more specific adaptive immunity. The generation of diversity in transcriptional programs is central for effective immune responses. Alternative splicing is widely considered a key generator of transcriptional and proteomic complexity, but its role has been rarely addressed systematically in immune cells. Here we used splicing-sensitive arrays to assess genome-wide gene- and exon-level expression profiles in human DCs in response to a bacterial challenge. We find widespread alternative splicing events and splicing factor transcriptional signatures induced by an E. coli challenge to human DCs. Alternative splicing acts in concert with transcriptional modulation, but these two mechanisms of gene regulation affect primarily distinct functional gene groups. Alternative splicing is likely to have an important role in DC immunobiology because it affects genes known to be involved in DC development, endocytosis, antigen presentation and cell cycle arrest.

Disease progression and phasic changes in gene expression in a mouse model of osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis and has multiple risk factors including joint injury. The purpose of this study was to characterize the histologic development of OA in a mouse model where OA is induced by destabilization of the medial meniscus (DMM model) and to identify genes regulated during different stages of the disease, using RNA isolated from the joint “organ” and analyzed using microarrays. Histologic changes seen in OA, including articular cartilage lesions and osteophytes, were present in the medial tibial plateaus of the DMM knees beginning at the earliest (2 week) time point and became progressively more severe by 16 weeks. 427 probe sets (371 genes) from the microarrays passed consistency and significance filters. There was an initial up-regulation at 2 and 4 weeks of genes involved in morphogenesis, differentiation, and development, including growth factor and matrix genes, as well as transcription factors including Atf2, Creb3l1, and Erg. Most genes were off or down-regulated at 8 weeks with the most highly down-regulated genes involved in cell division and the cytoskeleton. Gene expression increased at 16 weeks, in particular extracellular matrix genes including Prelp, Col3a1 and fibromodulin. Immunostaining revealed the presence of these three proteins in cartilage and soft tissues including ligaments as well as in the fibrocartilage covering osteophytes. The results support a phasic development of OA with early matrix remodeling and transcriptional activity followed by a more quiescent period that is not maintained. This implies that the response to an OA intervention will depend on the timing of the intervention. The quiescent period at 8 weeks may be due to the maturation of the osteophytes which are thought to temporarily stabilize the joint.

Continuous cotemporal probabilistic modeling of systems biology networks from sparse data

Modeling of biological networks is a difficult endeavor, but exploration of this problem is essential for understanding the systems behavior of biological processes. In this contribution, developed for sparse data, we present a new continuous Bayesian graphical learning algorithm to cotemporally model proteins in signaling networks and genes in transcriptional regulatory networks. In this continuous Bayesian algorithm, the correlation matrix is singular because the number of time points is less than the number of biological entities (genes or proteins). A suitable restriction on the degree of the graph's vertices is applied and a Metropolis-Hastings algorithm is guided by a BIC-based posterior probability score. Ten independent and diverse runs of the algorithm are conducted, so that the probability space is properly well-explored. Diagnostics to test the applicability of the algorithm to the specific data sets are developed; this is a major benefit of the methodology. This novel algorithm is applied to two time course experimental data sets: 1) protein modification data identifying a potential signaling network in chondrocytes, and 2) gene expression data identifying the transcriptional regulatory network underlying dendritic cell maturation. This method gives high estimated posterior probabilities to many of the proteins' directed edges that are predicted by the literature; for the gene study, the method gives high posterior probabilities to many of the literature-predicted sibling edges. In simulations, the method gives substantially higher estimated posterior probabilities for true edges and true subnetworks than for their false counterparts.

SC²ATmd: a tool for integration of the figure of merit with cluster analysis for gene expression data

Standard and Consensus Clustering Analysis Tool for Microarray Data (SC²ATmd) is a MATLAB-implemented application specifically designed for the exploration of microarray gene expression data via clustering. Implementation of two versions of the clustering validation method figure of merit allows for performance comparisons between different clustering algorithms, and tailors the cluster analysis process to the varying characteristics of each dataset. Along with standard clustering algorithms this application also offers a consensus clustering method that can generate reproducible clusters across replicate experiments or different clustering algorithms. This application was designed specifically for the analysis of gene expression data, but may be used with any numerical data as long as it is in the right format.

AVAILABILITY

SC²ATmd may be freely downloaded from http://www.compbiosci.wfu.edu/tools.htm.