A C-terminal motif found in the beta2-adrenergic receptor, P2Y1 receptor and cystic fibrosis transmembrane conductance regulator determines binding to the Na+/H+ exchanger regulatory factor family of PDZ proteins.

The Na+/H+ exchanger regulatory factor (NHERF) binds to the tail of the beta2-adrenergic receptor and plays a role in adrenergic regulation of Na+/H+ exchange. NHERF contains two PDZ domains, the first of which is required for its interaction with the beta2 receptor. Mutagenesis studies of the beta2 receptor tail revealed that the optimal C-terminal motif for binding to the first PDZ domain of NHERF is D-S/T-x-L, a motif distinct from those recognized by other PDZ domains. The first PDZ domain of NHERF-2, a protein that is 52% identical to NHERF and also known as E3KARP, SIP-1, and TKA-1, exhibits binding preferences very similar to those of the first PDZ domain of NHERF. The delineation of the preferred binding motif for the first PDZ domain of the NHERF family of proteins allows for predictions for other proteins that may interact with NHERF or NHERF-2. For example, as would be predicted from the beta2 receptor tail mutagenesis studies, NHERF binds to the tail of the purinergic P2Y1 receptor, a seven-transmembrane receptor with an intracellular C-terminal tail ending in D-T-S-L. NHERF also binds to the tail of the cystic fibrosis transmembrane conductance regulator, which ends in D-T-R-L. Because the preferred binding motif of the first PDZ domain of the NHERF family of proteins is found at the C termini of a variety of intracellular proteins, NHERF and NHERF-2 may be multifunctional adaptor proteins involved in many previously unsuspected aspects of intracellular signaling.

Evidence-ranked motif identification.

cERMIT is a computationally efficient motif discovery tool based on analyzing genome-wide quantitative regulatory evidence. Instead of pre-selecting promising candidate sequences, it utilizes information across all sequence regions to search for high-scoring motifs. We apply cERMIT on a range of direct binding and overexpression datasets; it substantially outperforms state-of-the-art approaches on curated ChIP-chip datasets, and easily scales to current mammalian ChIP-seq experiments with data on thousands of non-coding regions.

HDAC6 and Ubp-M BUZ domains recognize specific C-terminal sequences of proteins.

The BUZ/Znf-UBP domain is a protein module found in the cytoplasmic deacetylase HDAC6, E3 ubiquitin ligase BRAP2/IMP, and a subfamily of ubiquitin-specific proteases. Although several BUZ domains have been shown to bind ubiquitin with high affinity by recognizing its C-terminal sequence (RLRGG-COOH), it is currently unknown whether the interaction is sequence-specific or whether the BUZ domains are capable of binding to proteins other than ubiquitin. In this work, the BUZ domains of HDAC6 and Ubp-M were subjected to screening against a one-bead-one-compound (OBOC) peptide library that exhibited random peptide sequences with free C-termini. Sequence analysis of the selected binding peptides as well as alanine scanning studies revealed that the BUZ domains require a C-terminal Gly-Gly motif for binding. At the more N-terminal positions, the two BUZ domains have distinct sequence specificities, allowing them to bind to different peptides and/or proteins. A database search of the human proteome on the basis of the BUZ domain specificities identified 11 and 24 potential partner proteins for Ubp-M and HDAC6 BUZ domains, respectively. Peptides corresponding to the C-terminal sequences of four of the predicted binding partners (FBXO11, histone H4, PTOV1, and FAT10) were synthesized and tested for binding to the BUZ domains by fluorescence polarization. All four peptides bound to the HDAC6 BUZ domain with low micromolar K(D) values and less tightly to the Ubp-M BUZ domain. Finally, in vitro pull-down assays showed that the Ubp-M BUZ domain was capable of binding to the histone H3-histone H4 tetramer protein complex. Our results suggest that BUZ domains are sequence-specific protein-binding modules, with each BUZ domain potentially binding to a different subset of proteins.

Determinants of protection among HIV‐exposed seronegative persons: an overview.

Both clinical experience and a growing medical literature indicate that some persons who have been exposed to human immunodeficiency virus (HIV) infection remain uninfected. Although in some instances this may represent good fortune, cohorts of uninfected persons have been reported who are considered at high risk for infection. In these cohorts a variety of characteristics have been proposed as mediating protection, but to date only the 32–base pair deletion in the chemokine (C‐C motif) receptor 5 gene, which results in complete failure of cell surface expression of this coreceptor, has been associated with high‐level protection from HIV infection. With this in mind, there are probably many other factors that may individually or in combination provide some level of protection from acquisition of HIV infection. Because some of these factors are probably incompletely protective or inconsistently active, identifying them with confidence will be difficult. Nonetheless, clarifying the determinants of protection against HIV infection is a high priority that will require careful selection of high‐risk uninfected cohorts, who should undergo targeted studies of plausible mediators and broad screening for unexpected determinants of protection.

Bayesian variable selection in structured high-dimensional covariate spaces with applications in genomics

We consider the problem of variable selection in regression modeling in high-dimensional spaces where there is known structure among the covariates. This is an unconventional variable selection problem for two reasons: (1) The dimension of the covariate space is comparable, and often much larger, than the number of subjects in the study, and (2) the covariate space is highly structured, and in some cases it is desirable to incorporate this structural information in to the model building process. We approach this problem through the Bayesian variable selection framework, where we assume that the covariates lie on an undirected graph and formulate an Ising prior on the model space for incorporating structural information. Certain computational and statistical problems arise that are unique to such high-dimensional, structured settings, the most interesting being the phenomenon of phase transitions. We propose theoretical and computational schemes to mitigate these problems. We illustrate our methods on two different graph structures: the linear chain and the regular graph of degree k. Finally, we use our methods to study a specific application in genomics: the modeling of transcription factor binding sites in DNA sequences. © 2010 American Statistical Association.

Ubiquitylation of p53 by the APC/C inhibitor Trim39.

Tripartite motif 39 (Trim39) is a RING domain-containing E3 ubiquitin ligase able to inhibit the anaphase-promoting complex (APC/C) directly. Through analysis of Trim39 function in p53-positive and p53-negative cells, we have found, surprisingly, that p53-positive cells lacking Trim39 could not traverse the G1/S transition. This effect did not result from disinhibition of the APC/C. Moreover, although Trim39 loss inhibited etoposide-induced apoptosis in p53-negative cells, apoptosis was enhanced by Trim39 knockdown in p53-positive cells. Furthermore, we show here that the Trim39 can directly bind and ubiquitylate p53 in vitro and in vivo, leading to p53 degradation. Depletion of Trim39 significantly increased p53 protein levels and cell growth retardation in multiple cell lines. We found that the relative importance of Trim39 and the well-characterized p53-directed E3 ligase, murine double minute 2 (MDM2), varied between cell types. In cells that were relatively insensitive to the MDM2 inhibitor, nutlin-3a, apoptosis could be markedly enhanced by siRNA directed against Trim39. As such, Trim39 may serve as a potential therapeutic target in tumors with WT p53 when MDM2 inhibition is insufficient to elevate p53 levels and apoptosis.

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.

Nutritional Control of L1 Arrest and Recovery in Caenorhabditis elegans by Insulin-like Peptides and Signaling

Animals must coordinate development with fluctuating nutrient availability. Nutrient availability governs post-embryonic development in Caenorhabditis elegans: larvae that hatch in the absence of food do not initiate post-embryonic development but enter "L1 arrest" (or "L1 diapause") and can survive starvation for weeks, while rapidly resume normal development once get fed. Insulin-like signaling (IIS) has been shown to be a key regulator of L1 arrest and recovery. However, the C. elegans genome encodes 40 insulin-like peptides (ILPs), and it is unknown which peptides participate in nutritional control of L1 arrest and recovery. Work in other contexts has identified putative receptor agonists and antagonists, but the extent of specificity versus redundancy is unclear beyond this distinction.

We measured mRNA expression dynamics with high temporal resolution for all 40 insulin-like genes during entry into and recovery from L1 arrest. Nutrient availability influences expression of the majority of insulin-like genes, with variable dynamics suggesting complex regulation. We identified 13 candidate agonists and 8 candidate antagonists based on expression in response to nutrient availability. We selected ten candidate agonists (daf-28, ins-3, ins-4, ins-5, ins-6, ins-7, ins-9, ins-26, ins-33 and ins-35) for further characterization in L1 stage larvae. We used destabilized reporter genes to determine spatial expression patterns. Expression of candidate agonists was largely overlapping in L1 stage larvae, suggesting a role of the intestine, chemosensory neurons ASI and ASJ, and the interneuron PVT in systemic control of L1 development. Transcriptional regulation of candidate agonists was most significant in the intestine, as if nutrient uptake was a more important influence on transcription than sensory perception. Scanning in the 5' upstream promoter region of these 40 ILPs, We found that transcription factor PQM-1 and GATA putative binding sites are depleted in the promoter region of antagonists. A novel motif was also found to be over-represented in ILPs.

Phenotypic analysis of single and compound deletion mutants did not reveal effects on L1 recovery/developmental dynamics, though simultaneous disruption of ins-4 and daf-28 extended survival of L1 arrest without enhancing thermal tolerance, while overexpression of ins-4, ins-6 or daf-28 shortened L1 survival. Simultaneous disruption of several ILPs showed a temperature independent, transient dauer phenotype. These results revealed the relative redundancy and specificity among agonistic ILPs.

TGF- β and steroid hormone (SH) signaling have been reported to control the dauer formation along with IIS. Our preliminary results suggest they may also mediate the IIS control of L1 arrest and recovery, as the expression of several key components of TGF-β and SH signaling pathway genes are negatively regulated by DAF-16, and loss-of-function of these genes partially represses daf-16 null phenotype in L1 arrest, and causes a retardation in L1 development.

In summary, my dissertation study focused on the IIS, characterized the dynamics and sites of ILPs expression in response to nutrient availability, revealed the function of specific agonistic ILPs in L1 arrest, and suggested potential cross-regulation among IIS, TGF-β signaling and SH signaling in controlling L1 arrest and recovery. These findings provide insights into how post-embryonic development is governed by insulin-like signaling and nutrient availability.

Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex.

The human neocortex differs from that of other great apes in several notable regards, including altered cell cycle, prolonged corticogenesis, and increased size [1-5]. Although these evolutionary changes most likely contributed to the origin of distinctively human cognitive faculties, their genetic basis remains almost entirely unknown. Highly conserved non-coding regions showing rapid sequence changes along the human lineage are candidate loci for the development and evolution of uniquely human traits. Several studies have identified human-accelerated enhancers [6-14], but none have linked an expression difference to a specific organismal trait. Here we report the discovery of a human-accelerated regulatory enhancer (HARE5) of FZD8, a receptor of the Wnt pathway implicated in brain development and size [15, 16]. Using transgenic mice, we demonstrate dramatic differences in human and chimpanzee HARE5 activity, with human HARE5 driving early and robust expression at the onset of corticogenesis. Similar to HARE5 activity, FZD8 is expressed in neural progenitors of the developing neocortex [17-19]. Chromosome conformation capture assays reveal that HARE5 physically and specifically contacts the core Fzd8 promoter in the mouse embryonic neocortex. To assess the phenotypic consequences of HARE5 activity, we generated transgenic mice in which Fzd8 expression is under control of orthologous enhancers (Pt-HARE5::Fzd8 and Hs-HARE5::Fzd8). In comparison to Pt-HARE5::Fzd8, Hs-HARE5::Fzd8 mice showed marked acceleration of neural progenitor cell cycle and increased brain size. Changes in HARE5 function unique to humans thus alter the cell-cycle dynamics of a critical population of stem cells during corticogenesis and may underlie some distinctive anatomical features of the human brain.

Human Non-neutralizing HIV-1 Envelope Monoclonal Antibodies Limit the Number of Founder Viruses during SHIV Mucosal Infection in Rhesus Macaques.

HIV-1 mucosal transmission begins with virus or virus-infected cells moving through mucus across mucosal epithelium to infect CD4+ T cells. Although broadly neutralizing antibodies (bnAbs) are the type of HIV-1 antibodies that are most likely protective, they are not induced with current vaccine candidates. In contrast, antibodies that do not neutralize primary HIV-1 strains in the TZM-bl infection assay are readily induced by current vaccine candidates and have also been implicated as secondary correlates of decreased HIV-1 risk in the RV144 vaccine efficacy trial. Here, we have studied the capacity of anti-Env monoclonal antibodies (mAbs) against either the immunodominant region of gp41 (7B2 IgG1), the first constant region of gp120 (A32 IgG1), or the third variable loop (V3) of gp120 (CH22 IgG1) to modulate in vivo rectal mucosal transmission of a high-dose simian-human immunodeficiency virus (SHIV-BaL) in rhesus macaques. 7B2 IgG1 or A32 IgG1, each containing mutations to enhance Fc function, was administered passively to rhesus macaques but afforded no protection against productive clinical infection while the positive control antibody CH22 IgG1 prevented infection in 4 of 6 animals. Enumeration of transmitted/founder (T/F) viruses revealed that passive infusion of each of the three antibodies significantly reduced the number of T/F genomes. Thus, some antibodies that bind HIV-1 Env but fail to neutralize virus in traditional neutralization assays may limit the number of T/F viruses involved in transmission without leading to enhancement of viral infection. For one of these mAbs, gp41 mAb 7B2, we provide the first co-crystal structure in complex with a common cyclical loop motif demonstrated to be critical for infection by other retroviruses.

Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

ARFGAP1 promotes the formation of COPI vesicles, suggesting function as a component of the coat.

The role of GTPase-activating protein (GAP) that deactivates ADP-ribosylation factor 1 (ARF1) during the formation of coat protein I (COPI) vesicles has been unclear. GAP is originally thought to antagonize vesicle formation by triggering uncoating, but later studies suggest that GAP promotes cargo sorting, a process that occurs during vesicle formation. Recent models have attempted to reconcile these seemingly contradictory roles by suggesting that cargo proteins suppress GAP activity during vesicle formation, but whether GAP truly antagonizes coat recruitment in this process has not been assessed directly. We have reconstituted the formation of COPI vesicles by incubating Golgi membrane with purified soluble components, and find that ARFGAP1 in the presence of GTP promotes vesicle formation and cargo sorting. Moreover, the presence of GTPgammaS not only blocks vesicle uncoating but also vesicle formation by preventing the proper recruitment of GAP to nascent vesicles. Elucidating how GAP functions in vesicle formation, we find that the level of GAP on the reconstituted vesicles is at least as abundant as COPI and that GAP binds directly to the dilysine motif of cargo proteins. Collectively, these findings suggest that ARFGAP1 promotes vesicle formation by functioning as a component of the COPI coat.