High-field FMRI reveals brain activation patterns underlying saccade execution in the human superior colliculus.

BACKGROUND: The superior colliculus (SC) has been shown to play a crucial role in the initiation and coordination of eye- and head-movements. The knowledge about the function of this structure is mainly based on single-unit recordings in animals with relatively few neuroimaging studies investigating eye-movement related brain activity in humans. METHODOLOGY/PRINCIPAL FINDINGS: The present study employed high-field (7 Tesla) functional magnetic resonance imaging (fMRI) to investigate SC responses during endogenously cued saccades in humans. In response to centrally presented instructional cues, subjects either performed saccades away from (centrifugal) or towards (centripetal) the center of straight gaze or maintained fixation at the center position. Compared to central fixation, the execution of saccades elicited hemodynamic activity within a network of cortical and subcortical areas that included the SC, lateral geniculate nucleus (LGN), occipital cortex, striatum, and the pulvinar. CONCLUSIONS/SIGNIFICANCE: Activity in the SC was enhanced contralateral to the direction of the saccade (i.e., greater activity in the right as compared to left SC during leftward saccades and vice versa) during both centrifugal and centripetal saccades, thereby demonstrating that the contralateral predominance for saccade execution that has been shown to exist in animals is also present in the human SC. In addition, centrifugal saccades elicited greater activity in the SC than did centripetal saccades, while also being accompanied by an enhanced deactivation within the prefrontal default-mode network. This pattern of brain activity might reflect the reduced processing effort required to move the eyes toward as compared to away from the center of straight gaze, a position that might serve as a spatial baseline in which the retinotopic and craniotopic reference frames are aligned.

Electron tomography of cryofixed, isometrically contracting insect flight muscle reveals novel actin-myosin interactions.

BACKGROUND: Isometric muscle contraction, where force is generated without muscle shortening, is a molecular traffic jam in which the number of actin-attached motors is maximized and all states of motor action are trapped with consequently high heterogeneity. This heterogeneity is a major limitation to deciphering myosin conformational changes in situ. METHODOLOGY: We used multivariate data analysis to group repeat segments in electron tomograms of isometrically contracting insect flight muscle, mechanically monitored, rapidly frozen, freeze substituted, and thin sectioned. Improved resolution reveals the helical arrangement of F-actin subunits in the thin filament enabling an atomic model to be built into the thin filament density independent of the myosin. Actin-myosin attachments can now be assigned as weak or strong by their motor domain orientation relative to actin. Myosin attachments were quantified everywhere along the thin filament including troponin. Strong binding myosin attachments are found on only four F-actin subunits, the "target zone", situated exactly midway between successive troponin complexes. They show an axial lever arm range of 77°/12.9 nm. The lever arm azimuthal range of strong binding attachments has a highly skewed, 127° range compared with X-ray crystallographic structures. Two types of weak actin attachments are described. One type, found exclusively in the target zone, appears to represent pre-working-stroke intermediates. The other, which contacts tropomyosin rather than actin, is positioned M-ward of the target zone, i.e. the position toward which thin filaments slide during shortening. CONCLUSION: We present a model for the weak to strong transition in the myosin ATPase cycle that incorporates azimuthal movements of the motor domain on actin. Stress/strain in the S2 domain may explain azimuthal lever arm changes in the strong binding attachments. The results support previous conclusions that the weak attachments preceding force generation are very different from strong binding attachments.

An approach to the study of G-protein-coupled receptor kinases: an in vitro-purified membrane assay reveals differential receptor specificity and regulation by G beta gamma subunits.

Phosphorylation of GTP-binding-regulatory (G)-protein-coupled receptors by specific G-protein-coupled receptor kinases (GRKs) is a major mechanism responsible for agonist-mediated desensitization of signal transduction processes. However, to date, studies of the specificity of these enzymes have been hampered by the difficulty of preparing the purified and reconstituted receptor preparations required as substrates. Here we describe an approach that obviates this problem by utilizing highly purified membrane preparations from Sf9 and 293 cells overexpressing G-protein-coupled receptors. We use this technique to demonstrate specificity of several GRKs with respect to both receptor substrates and the enhancing effects of G-protein beta gamma subunits on phosphorylation. Enriched membrane preparations of the beta 2- and alpha 2-C2-adrenergic receptors (ARs, where alpha 2-C2-AR refers to the AR whose gene is located on human chromosome 2) prepared by sucrose density gradient centrifugation from Sf9 or 293 cells contain the receptor at 100-300 pmol/mg of protein and serve as efficient substrates for agonist-dependent phosphorylation by beta-AR kinase 1 (GRK2), beta-AR kinase 2 (GRK3), or GRK5. Stoichiometries of agonist-mediated phosphorylation of the receptors by GRK2 (beta-AR kinase 1), in the absence and presence of G beta gamma, are 1 and 3 mol/mol, respectively. The rate of phosphorylation of the membrane receptors is 3 times faster than that of purified and reconstituted receptors. While phosphorylation of the beta 2-AR by GRK2, -3, and -5 is similar, the activity of GRK2 and -3 is enhanced by G beta gamma whereas that of GRK5 is not. In contrast, whereas GRK2 and -3 efficiently phosphorylate alpha 2-C2-AR, GRK5 is quite weak. The availability of a simple direct phosphorylation assay applicable to any cloned G-protein-coupled receptor should greatly facilitate elucidation of the mechanisms of regulation of these receptors by the expanding family of GRKs.

The receptor kinase family: primary structure of rhodopsin kinase reveals similarities to the beta-adrenergic receptor kinase.

Light-dependent deactivation of rhodopsin as well as homologous desensitization of beta-adrenergic receptors involves receptor phosphorylation that is mediated by the highly specific protein kinases rhodopsin kinase (RK) and beta-adrenergic receptor kinase (beta ARK), respectively. We report here the cloning of a complementary DNA for RK. The deduced amino acid sequence shows a high degree of homology to beta ARK. In a phylogenetic tree constructed by comparing the catalytic domains of several protein kinases, RK and beta ARK are located on a branch close to, but separate from the cyclic nucleotide-dependent protein kinase and protein kinase C subfamilies. From the common structural features we conclude that both RK and beta ARK are members of a newly delineated gene family of guanine nucleotide-binding protein (G protein)-coupled receptor kinases that may function in diverse pathways to regulate the function of such receptors.

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.

Integrative Genomics Reveals a Role for GNA13 in Lymphomagenesis

Lymphomas comprise a diverse group of malignancies derived from immune cells. High throughput sequencing has recently emerged as a powerful and versatile method for analysis of the cancer genome and transcriptome. As these data continue to emerge, the crucial work lies in sorting through the wealth of information to hone in on the critical aspects that will give us a better understanding of biology and new insight for how to treat disease. Finding the important signals within these large data sets is one of the major challenges of next generation sequencing.

In this dissertation, I have developed several complementary strategies to describe the genetic underpinnings of lymphomas. I begin with developing a better method for RNA sequencing that enables strand-specific total RNA sequencing and alternative splicing profiling in the same analysis. I then combine this RNA sequencing technique with whole exome sequencing to better understand the global landscape of aberrations in these diseases. Finally, I use traditional cell and molecular biology techniques to define the consequences of major genetic alterations in lymphoma.

Through this analysis, I find recurrent silencing mutations in the G alpha binding protein GNA13 and associated focal adhesion proteins. I aim to describe how loss-of-function mutations in GNA13 can be oncogenic in the context of germinal center B cell biology. Using in vitro techniques including liquid chromatography-mass spectrometry and knockdown and overexpression of genes in B cell lymphoma cell lines, I determine protein binding partners and downstream effectors of GNA13. I also develop a transgenic mouse model to study the role of GNA13 in the germinal center in vivo to determine effects of GNA13 deletion on germinal center structure and cell migration.

Thus, I have developed complementary approaches that span the spectrum from discovery to context-dependent gene models that afford a better understanding of the biological function of aberrant events and ultimately result in a better understanding of disease.

Conformational kinetics reveals affinities of protein conformational states.

Most biological reactions rely on interplay between binding and changes in both macromolecular structure and dynamics. Practical understanding of this interplay requires detection of critical intermediates and determination of their binding and conformational characteristics. However, many of these species are only transiently present and they have often been overlooked in mechanistic studies of reactions that couple binding to conformational change. We monitored the kinetics of ligand-induced conformational changes in a small protein using six different ligands. We analyzed the kinetic data to simultaneously determine both binding affinities for the conformational states and the rate constants of conformational change. The approach we used is sufficiently robust to determine the affinities of three conformational states and detect even modest differences in the protein's affinities for relatively similar ligands. Ligand binding favors higher-affinity conformational states by increasing forward conformational rate constants and/or decreasing reverse conformational rate constants. The amounts by which forward rate constants increase and reverse rate constants decrease are proportional to the ratio of affinities of the conformational states. We also show that both the affinity ratio and another parameter, which quantifies the changes in conformational rate constants upon ligand binding, are strong determinants of the mechanism (conformational selection and/or induced fit) of molecular recognition. Our results highlight the utility of analyzing the kinetics of conformational changes to determine affinities that cannot be determined from equilibrium experiments. Most importantly, they demonstrate an inextricable link between conformational dynamics and the binding affinities of conformational states.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Live Imaging of Host-Parasite Interactions in a Zebrafish Infection Model Reveals Cryptococcal Determinants of Virulence and Central Nervous System Invasion.

UNLABELLED: The human fungal pathogen Cryptococcus neoformans is capable of infecting a broad range of hosts, from invertebrates like amoebas and nematodes to standard vertebrate models such as mice and rabbits. Here we have taken advantage of a zebrafish model to investigate host-pathogen interactions of Cryptococcus with the zebrafish innate immune system, which shares a highly conserved framework with that of mammals. Through live-imaging observations and genetic knockdown, we establish that macrophages are the primary immune cells responsible for responding to and containing acute cryptococcal infections. By interrogating survival and cryptococcal burden following infection with a panel of Cryptococcus mutants, we find that virulence factors initially identified as important in causing disease in mice are also necessary for pathogenesis in zebrafish larvae. Live imaging of the cranial blood vessels of infected larvae reveals that C. neoformans is able to penetrate the zebrafish brain following intravenous infection. By studying a C. neoformans FNX1 gene mutant, we find that blood-brain barrier invasion is dependent on a known cryptococcal invasion-promoting pathway previously identified in a murine model of central nervous system invasion. The zebrafish-C. neoformans platform provides a visually and genetically accessible vertebrate model system for cryptococcal pathogenesis with many of the advantages of small invertebrates. This model is well suited for higher-throughput screening of mutants, mechanistic dissection of cryptococcal pathogenesis in live animals, and use in the evaluation of therapeutic agents. IMPORTANCE: Cryptococcus neoformans is an important opportunistic pathogen that is estimated to be responsible for more than 600,000 deaths worldwide annually. Existing mammalian models of cryptococcal pathogenesis are costly, and the analysis of important pathogenic processes such as meningitis is laborious and remains a challenge to visualize. Conversely, although invertebrate models of cryptococcal infection allow high-throughput assays, they fail to replicate the anatomical complexity found in vertebrates and, specifically, cryptococcal stages of disease. Here we have utilized larval zebrafish as a platform that overcomes many of these limitations. We demonstrate that the pathogenesis of C. neoformans infection in zebrafish involves factors identical to those in mammalian and invertebrate infections. We then utilize the live-imaging capacity of zebrafish larvae to follow the progression of cryptococcal infection in real time and establish a relevant model of the critical central nervous system infection phase of disease in a nonmammalian model.

The Physarum polycephalum Genome Reveals Extensive Use of Prokaryotic Two-Component and Metazoan-Type Tyrosine Kinase Signaling.

Physarum polycephalum is a well-studied microbial eukaryote with unique experimental attributes relative to other experimental model organisms. It has a sophisticated life cycle with several distinct stages including amoebal, flagellated, and plasmodial cells. It is unusual in switching between open and closed mitosis according to specific life-cycle stages. Here we present the analysis of the genome of this enigmatic and important model organism and compare it with closely related species. The genome is littered with simple and complex repeats and the coding regions are frequently interrupted by introns with a mean size of 100 bases. Complemented with extensive transcriptome data, we define approximately 31,000 gene loci, providing unexpected insights into early eukaryote evolution. We describe extensive use of histidine kinase-based two-component systems and tyrosine kinase signaling, the presence of bacterial and plant type photoreceptors (phytochromes, cryptochrome, and phototropin) and of plant-type pentatricopeptide repeat proteins, as well as metabolic pathways, and a cell cycle control system typically found in more complex eukaryotes. Our analysis characterizes P. polycephalum as a prototypical eukaryote with features attributed to the last common ancestor of Amorphea, that is, the Amoebozoa and Opisthokonts. Specifically, the presence of tyrosine kinases in Acanthamoeba and Physarum as representatives of two distantly related subdivisions of Amoebozoa argues against the later emergence of tyrosine kinase signaling in the opisthokont lineage and also against the acquisition by horizontal gene transfer.

Single-cell functional analysis of parathyroid adenomas reveals distinct classes of calcium sensing behaviour in primary hyperparathyroidism.

Primary hyperparathyroidism (PHPT) is a common endocrine neoplastic disorder caused by a failure of calcium sensing secondary to tumour development in one or more of the parathyroid glands. Parathyroid adenomas are comprised of distinct cellular subpopulations of variable clonal status that exhibit differing degrees of calcium responsiveness. To gain a clearer understanding of the relationship among cellular identity, tumour composition and clinical biochemistry in PHPT, we developed a novel single cell platform for quantitative evaluation of calcium sensing behaviour in freshly resected human parathyroid tumour cells. Live-cell intracellular calcium flux was visualized through Fluo-4-AM epifluorescence, followed by in situ immunofluorescence detection of the calcium sensing receptor (CASR), a central component in the extracellular calcium signalling pathway. The reactivity of individual parathyroid tumour cells to extracellular calcium stimulus was highly variable, with discrete kinetic response patterns observed both between and among parathyroid tumour samples. CASR abundance was not an obligate determinant of calcium responsiveness. Calcium EC50 values from a series of parathyroid adenomas revealed that the tumours segregated into two distinct categories. One group manifested a mean EC50 of 2.40 mM (95% CI: 2.37-2.41), closely aligned to the established normal range. The second group was less responsive to calcium stimulus, with a mean EC50 of 3.61 mM (95% CI: 3.45-3.95). This binary distribution indicates the existence of a previously unappreciated biochemical sub-classification of PHPT tumours, possibly reflecting distinct etiological mechanisms. Recognition of quantitative differences in calcium sensing could have important implications for the clinical management of PHPT.