Description: Comparative functional genomics studies the evolution of biological processes by analyzing functional data, such as gene expression profiles, across species. A major challenge is to compare profiles collected in a complex phylogeny. Here, we present Arboretum, a novel scalable computational algorithm that integrates expression data from multiple species with species and gene phylogenies to infer modules of coexpressed genes in extant species and their evolutionary histories. We also develop new, generally applicable measures of conservation and divergence in gene regulatory modules to assess the impact of changes in gene content and expression on module evolution. We used Arboretum to study the evolution of the transcriptional response to heat shock in eight species of Ascomycota fungi and to reconstruct modules of the ancestral environmental stress response (ESR). We found substantial conservation in the stress response across species and in the reconstructed components of the ancestral ESR modules. The greatest divergence was in the most induced stress, primarily through module expansion. The divergence of the heat stress response exceeds that observed in the response to glucose depletion in the same species. Arboretum and its associated analyses provide a comprehensive framework to systematically study regulatory evolution of condition-specific responses.

Description: Artemisinin-derived monomers and dimers inhibit human cytomegalovirus (CMV) replication in human foreskin fibroblasts (HFFs). The monomer artesunate (AS) inhibits CMV at micromolar concentrations, while dimers inhibit CMV replication at nanomolar concentrations, without increased toxicity in HFFs. We report on the variable anti-CMV activity of AS compared to the consistent and reproducible CMV inhibition by dimer 606 and ganciclovir (GCV). Investigation of this phenomenon revealed that the anti-CMV activity of AS correlated with HFFs synchronized to the G 0 /G 1 stage of the cell cycle. In contact-inhibited serum-starved HFFs or cells arrested at early/late G 1 with specific checkpoint regulators, AS and dimer 606 efficiently inhibited CMV replication. However, in cycling HFFs, in which CMV replication was productive, virus inhibition by AS was significantly reduced, but inhibition by dimer 606 and GCV was maintained. Cell cycle analysis in noninfected HFFs revealed that AS induced early G 1 arrest, while dimer 606 partially blocked cell cycle progression. In infected HFFs, AS and dimer 606 prevented the progression of cell cycle toward the G 1 /S checkpoint. AS reduced the expression of cyclin-dependent kinases (CDK) 2, 4, and 6 in noninfected cycling HFFs, while the effect of dimer 606 on these CDKs was moderate. Neither compound affected CDK expression in noninfected contact-inhibited HFFs. In CMV-infected cells, AS activity correlated with reduced CDK2 levels. CMV inhibition by AS and dimer 606 also correlated with hypophosphorylation (activity) of the retinoblastoma protein (pRb). AS activity was strongly associated with pRb hypophosphorylation, while its reduced anti-CMV activity was marked by pRb phosphorylation. Roscovitine, a CDK2 inhibitor, antagonized the anti-CMV activities of AS and dimer 606. These data suggest that cell cycle modulation through CDKs and pRb might play a role in the anti-CMV activities of artemisinins. Proteins involved in this modulation may be identified and targeted for CMV inhibition.

Description: Chromatin boundary elements (CBEs) are widely distributed in the genome and mediate formation of chromatin loops, but their roles in gene regulation remain poorly understood. The complex expression pattern of the Drosophila homeotic gene Sex combs reduced ( Scr ) is directed by an unusually long regulatory sequence harboring diverse cis elements and an intervening neighbor gene fushi tarazu ( ftz ). Here we report the presence of a multitude of CBEs in the Scr regulatory region. Selective and dynamic pairing among these CBEs mediates developmentally regulated chromatin loops. In particular, the SF1 boundary plays a central role in organizing two subsets of chromatin loops: one subset encloses ftz , limiting its access by the surrounding Scr enhancers and compartmentalizing distinct histone modifications, and the other subset subdivides the Scr regulatory sequences into independent enhancer access domains. We show that these CBEs exhibit diverse enhancer-blocking activities that vary in strength and tissue distribution. Tandem pairing of SF1 and SF2, two strong CBEs that flank the ftz domain, allows the distal enhancers to bypass their block in transgenic Drosophila , providing a mechanism for the endogenous Scr enhancer to circumvent the ftz domain. Our study demonstrates how an endogenous CBE network, centrally orchestrated by SF1, could remodel the genomic environment to facilitate gene regulation during development.

Description: Human La protein is known to be an essential host factor for translation and replication of hepatitis C virus (HCV) RNA. Previously, we have demonstrated that residues responsible for interaction of human La protein with the HCV internal ribosomal entry site (IRES) around the initiator AUG within stem-loop IV form a β-turn in the RNA recognition motif (RRM) structure. In this study, sequence alignment and mutagenesis suggest that the HCV RNA-interacting β-turn is conserved only in humans and chimpanzees, the species primarily known to be infected by HCV. A 7-mer peptide corresponding to the HCV RNA-interacting region of human La inhibits HCV translation, whereas another peptide corresponding to the mouse La sequence was unable to do so. Furthermore, IRES-mediated translation was found to be significantly high in the presence of recombinant human La protein in vitro in rabbit reticulocyte lysate. We observed enhanced replication with HCV subgenomic and full-length replicons upon overexpression of either human La protein or a chimeric mouse La protein harboring a human La β-turn sequence in mouse cells. Taken together, our results raise the possibility of creating an immunocompetent HCV mouse model using human-specific cell entry factors and a humanized form of La protein. IMPORTANCE Hepatitis C virus is known to infect only humans and chimpanzees under natural conditions. This has prevented the development of a small-animal model, which is important for development of new antiviral drugs. Although a number of human-specific proteins are responsible for this species selectivity and some of these proteins—-mostly entry factors—-have been identified, full multiplication of the virus in mouse cells is still not possible. In this study, we show that a turn in the human La protein that is responsible for the interaction with the viral RNA is highly specific for the human sequence. Replacement of the corresponding mouse sequence with the human sequence allows the mouse La to behave like its human counterpart and support viral growth in the mouse cell efficiently. This observation, in combination with previously identified cell entry factors, should open up the possibility of creating a mouse model of hepatitis C.

Description: Changes in chromatin state play important roles in cell fate transitions. Current computational approaches to analyze chromatin modifications across multiple cell types do not model how the cell types are related on a lineage or over time . To overcome this limitation, we developed a method called Chromatin Module INference on Trees (CMINT), a probabilistic clustering approach to systematically capture chromatin state dynamics across multiple cell types. Compared to existing approaches, CMINT can handle complex lineage topologies, capture higher quality clusters, and reliably detect chromatin transitions between cell types. We applied CMINT to gain novel insights in two complex processes: reprogramming to induced pluripotent stem cells (iPSCs) and hematopoiesis. In reprogramming, chromatin changes could occur without large gene expression changes, different combinations of activating marks were associated with specific reprogramming factors, there was an order of acquisition of chromatin marks at pluripotency loci, and multivalent states (comprising previously undetermined combinations of activating and repressive histone modifications) were enriched for CTCF. In the hematopoietic system, we defined critical decision points in the lineage tree, identified regulatory elements that were enriched in cell-type–specific regions, and found that the underlying chromatin state was achieved by specific erasure of preexisting chromatin marks in the precursor cell or by de novo assembly. Our method provides a systematic approach to model the dynamics of chromatin state to provide novel insights into the relationships among cell types in diverse cell-fate specification processes.

Description: Porphyromonas gingivalis is a keystone pathogen that contributes to periodontal pathogenesis by disrupting host-microbe homeostasis and promoting dysbiosis. The virulence of P. gingivalis likely reflects an alteration in the lipid A composition of its lipopolysaccharide (LPS) from the penta-acylated ( Pg LPS 1690 ) to the tetra-acylated ( Pg LPS 1435/1449 ) form. Mast cells play an important role in periodontitis, but the mechanisms of their activation and regulation remain unknown. The expression of epithelium- and neutrophil-derived host defense peptides (HDPs) (LL-37 and human β-defensin-3), which activate mast cells via Mas-related G protein-coupled receptor X2 (MRGPRX2), is increased in periodontitis. We found that MRGPRX2-expressing mast cells are present in normal gingiva and that their numbers are elevated in patients with chronic periodontitis. Furthermore, HDPs stimulated degranulation in a human mast cell line (LAD2) and in RBL-2H3 cells stably expressing MRGPRX2 (RBL-MRGPRX2). Pg LPS 1690 caused substantial inhibition of HDP-induced mast cell degranulation, but Pg LPS 1435/1449 had no effect. A fluorescently labeled HDP (FAM-LL-37) bound to RBL-MRGPRX2 cells, and Pg LPS 1690 inhibited this binding, but Pg LPS 1435/1449 had no effect. These findings suggest that low-level inflammation induced by HDP/MRGPRX2-mediated mast cell degranulation contributes to gingival homeostasis but that sustained inflammation due to elevated levels of both HDPs and MRGPRX2-expressing mast cells promotes periodontal disease. Furthermore, differential regulation of HDP-induced mast cell degranulation by Pg LPS 1690 and Pg LPS 1435/1449 may contribute to the modulation of disease progression.

Description: Among emerging non- albicans Candida species, Candida parapsilosis is of particular concern as a cause of nosocomial bloodstream infections in neonatal and intensive care unit patients. While fluconazole and echinocandins are considered effective treatments for such infections, recent reports of fluconazole and echinocandin resistance in C. parapsilosis indicate a growing problem. The present study describes a novel mechanism of antifungal resistance in this organism affecting susceptibility to azole and echinocandin antifungals in a clinical isolate obtained from a patient with prosthetic valve endocarditis. Transcriptome analysis indicated differential expression of several genes in the resistant isolate, including upregulation of ergosterol biosynthesis pathway genes ERG2 , ERG5 , ERG6 , ERG11 , ERG24 , ERG25 , and UPC2 . Whole-genome sequencing revealed that the resistant isolate possessed an ERG3 mutation resulting in a G111R amino acid substitution. Sterol profiles indicated a reduction in sterol desaturase activity as a result of this mutation. Replacement of both mutant alleles in the resistant isolate with the susceptible isolate's allele restored wild-type susceptibility to all azoles and echinocandins tested. Disruption of ERG3 in the susceptible and resistant isolates resulted in a loss of sterol desaturase activity, high-level azole resistance, and an echinocandin-intermediate to -resistant phenotype. While disruption of ERG3 in C. albicans resulted in azole resistance, echinocandin MICs, while elevated, remained within the susceptible range. This work demonstrates that the G111R substitution in Erg3 is wholly responsible for the altered azole and echinocandin susceptibilities observed in this C. parapsilosis isolate and is the first report of an ERG3 mutation influencing susceptibility to the echinocandins.

Description: Gaining insights on gene regulation from large-scale functional data sets is a grand challenge in systems biology. In this article, we develop and apply methods for transcriptional regulatory network inference from diverse functional genomics data sets and demonstrate their value for gene function and gene expression prediction. We formulate the network inference problem in a machine-learning framework and use both supervised and unsupervised methods to predict regulatory edges by integrating transcription factor (TF) binding, evolutionarily conserved sequence motifs, gene expression, and chromatin modification data sets as input features. Applying these methods to Drosophila melanogaster, we predict ~300,000 regulatory edges in a network of ~600 TFs and 12,000 target genes. We validate our predictions using known regulatory interactions, gene functional annotations, tissue-specific expression, protein–protein interactions, and three-dimensional maps of chromosome conformation. We use the inferred network to identify putative functions for hundreds of previously uncharacterized genes, including many in nervous system development, which are independently confirmed based on their tissue-specific expression patterns. Last, we use the regulatory network to predict target gene expression levels as a function of TF expression, and find significantly higher predictive power for integrative networks than for motif or ChIP-based networks. Our work reveals the complementarity between physical evidence of regulatory interactions (TF binding, motif conservation) and functional evidence (coordinated expression or chromatin patterns) and demonstrates the power of data integration for network inference and studies of gene regulation at the systems level.

Description: The increased awareness of the role of environmental matrices in enteric disease transmission has resulted in the need for rapid, field-based methods for fecal indicator bacteria and pathogen detection. Evidence of the specificity of β-glucuronidase-based assays for detection of Escherichia coli from environmental matrices relevant to enteric pathogen transmission in developing countries, such as hands, soils, and surfaces, is limited. In this study, we quantify the false-positive rate of a β-glucuronidase-based E. coli detection assay (Colilert) for two environmental reservoirs in Bangladeshi households (hands and soils) and three fecal composite sources (cattle, chicken, and humans). We investigate whether or not the isolation source of E. coli influences phenotypic and genotypic characteristics. Phenotypic characteristics include results of biochemical assays provided by the API-20E test; genotypic characteristics include the Clermont phylogroup and the presence of enteric and/or environmental indicator genes sfmH , rfaI , and fucK . Our findings demonstrate no statistically significant difference in the false-positive rate of Colilert for environmental compared to enteric samples. E. coli isolates from all source types are genetically diverse, representing six of the seven phylogroups, and there is no difference in relative frequency of phylogroups between enteric and environmental samples. We conclude that Colilert, and likely other β-glucuronidase-based assays, is appropriate for detection of E. coli on hands and in soils with low false-positive rates. Furthermore, E. coli isolated from hands and soils in Bangladeshi households are diverse and indistinguishable from cattle, chicken, and human fecal isolates, using traditional biochemical assays and phylogrouping.