Genome-wide screen reveals Rab12 GTPase as a critical activator of Parkinson's disease-linked LRRK2 kinase.

Activating mutations in the leucine-rich repeat kinase 2 (LRRK2) cause Parkinson's disease. LRRK2 phosphorylates a subset of Rab GTPases, particularly Rab10 and Rab8A, and we showed previously that these phosphoRabs play an important role in LRRK2 membrane recruitment and activation (Vides et al., 2022). To learn more about LRRK2 pathway regulation, we carried out an unbiased, CRISPR-based genome-wide screen to identify modifiers of cellular phosphoRab10 levels. A flow cytometry assay was developed to detect changes in phosphoRab10 levels in pools of mouse NIH-3T3 cells harboring unique CRISPR guide sequences. Multiple negative and positive regulators were identified; surprisingly, knockout of the Rab12 gene was especially effective in decreasing phosphoRab10 levels in multiple cell types and knockout mouse tissues. Rab-driven increases in phosphoRab10 were specific for Rab12, LRRK2-dependent and PPM1H phosphatase-reversible, and did not require Rab12 phosphorylation; they were seen with wild type and pathogenic G2019S and R1441C LRRK2. As expected for a protein that regulates LRRK2 activity, Rab12 also influenced primary cilia formation. AlphaFold modeling revealed a novel Rab12 binding site in the LRRK2 Armadillo domain, and we show that residues predicted to be essential for Rab12 interaction at this site influence phosphoRab10 and phosphoRab12 levels in a manner distinct from Rab29 activation of LRRK2. Our data show that Rab12 binding to a new site in the LRRK2 Armadillo domain activates LRRK2 kinase for Rab phosphorylation and could serve as a new therapeutic target for a novel class of LRRK2 inhibitors that do not target the kinase domain.

Global mapping of RNA-chromatin contacts reveals a proximity-dominated connectivity model for ncRNA-gene interactions.

Non-coding RNAs (ncRNAs) are transcribed throughout the genome and provide regulatory inputs to gene expression through their interaction with chromatin. Yet, the genomic targets and functions of most ncRNAs are unknown. Here we use chromatin-associated RNA sequencing (ChAR-seq) to map the global network of ncRNA interactions with chromatin in human embryonic stem cells and the dynamic changes in interactions during differentiation into definitive endoderm. We uncover general principles governing the organization of the RNA-chromatin interactome, demonstrating that nearly all ncRNAs exclusively interact with genes in close three-dimensional proximity to their locus and provide a model predicting the interactome. We uncover RNAs that interact with many loci across the genome and unveil thousands of unannotated RNAs that dynamically interact with chromatin. By relating the dynamics of the interactome to changes in gene expression, we demonstrate that activation or repression of individual genes is unlikely to be controlled by a single ncRNA.

A feed-forward pathway drives LRRK2 kinase membrane recruitment and activation.

Activating mutations in the leucine-rich repeat kinase 2 (LRRK2) cause Parkinson's disease, and previously we showed that activated LRRK2 phosphorylates a subset of Rab GTPases (Steger et al., 2017). Moreover, Golgi-associated Rab29 can recruit LRRK2 to the surface of the Golgi and activate it there for both auto- and Rab substrate phosphorylation. Here, we define the precise Rab29 binding region of the LRRK2 Armadillo domain between residues 360-450 and show that this domain, termed 'site #1,' can also bind additional LRRK2 substrates, Rab8A and Rab10. Moreover, we identify a distinct, N-terminal, higher-affinity interaction interface between LRRK2 phosphorylated Rab8 and Rab10 termed 'site #2' that can retain LRRK2 on membranes in cells to catalyze multiple, subsequent phosphorylation events. Kinase inhibitor washout experiments demonstrate that rapid recovery of kinase activity in cells depends on the ability of LRRK2 to associate with phosphorylated Rab proteins, and phosphorylated Rab8A stimulates LRRK2 phosphorylation of Rab10 in vitro. Reconstitution of purified LRRK2 recruitment onto planar lipid bilayers decorated with Rab10 protein demonstrates cooperative association of only active LRRK2 with phospho-Rab10-containing membrane surfaces. These experiments reveal a feed-forward pathway that provides spatial control and membrane activation of LRRK2 kinase activity.

From telomere to telomere: The transcriptional and epigenetic state of human repeat elements.

Mobile elements and repetitive genomic regions are sources of lineage-specific genomic innovation and uniquely fingerprint individual genomes. Comprehensive analyses of such repeat elements, including those found in more complex regions of the genome, require a complete, linear genome assembly. We present a de novo repeat discovery and annotation of the T2T-CHM13 human reference genome. We identified previously unknown satellite arrays, expanded the catalog of variants and families for repeats and mobile elements, characterized classes of complex composite repeats, and located retroelement transduction events. We detected nascent transcription and delineated CpG methylation profiles to define the structure of transcriptionally active retroelements in humans, including those in centromeres. These data expand our insight into the diversity, distribution, and evolution of repetitive regions that have shaped the human genome.

Identification and characterization of centromeric sequences in Xenopus laevis.

Centromeres play an essential function in cell division by specifying the site of kinetochore formation on each chromosome for mitotic spindle attachment. Centromeres are defined epigenetically by the histone H3 variant Centromere Protein A (Cenpa). Cenpa nucleosomes maintain the centromere by designating the site for new Cenpa assembly after dilution by replication. Vertebrate centromeres assemble on tandem arrays of repetitive sequences, but the function of repeat DNA in centromere formation has been challenging to dissect due to the difficulty in manipulating centromeres in cells. Xenopus laevis egg extracts assemble centromeres in vitro, providing a system for studying centromeric DNA functions. However, centromeric sequences in Xenopus laevis have not been extensively characterized. In this study, we combine Cenpa ChIP-seq with a k-mer based analysis approach to identify the Xenopus laevis centromere repeat sequences. By in situ hybridization, we show that Xenopus laevis centromeres contain diverse repeat sequences, and we map the centromere position on each Xenopus laevis chromosome using the distribution of centromere-enriched k-mers. Our identification of Xenopus laevis centromere sequences enables previously unapproachable centromere genomic studies. Our approach should be broadly applicable for the analysis of centromere and other repetitive sequences in any organism.

Mapping Transcriptome-Wide and Genome-Wide RNA-DNA Contacts with Chromatin-Associated RNA Sequencing (ChAR-seq).

RNAs play key roles in the cell as molecular intermediates for protein synthesis and as regulators of nuclear processes such as splicing, posttranscriptional regulation, or chromatin remodeling. Various classes of non-coding RNAs, including long non-coding RNAs (lncRNAs), can bind chromatin either directly or via interaction with chromatin binding proteins. It has been proposed that lncRNAs regulate cell-state-specific genes by coordinating the locus-dependent activity of chromatin-modifying complexes. Yet, the vast majority of lncRNAs have unknown functions, and we know little about the specific loci they regulate. A key step toward understanding chromatin regulation by RNAs is to map the genomic loci with which every nuclear RNA interacts and, reciprocally, to identify all RNAs that target a given locus. Our ability to generate such data has been limited, until recently, by the lack of methods to probe the genomic localization of more than a few RNAs at a time. Here, we describe a protocol for ChAR-seq, an RNA-DNA proximity ligation method that maps the binding loci for thousands of RNAs at once and without the need for specific RNA or DNA probe sequences. The ChAR-seq approach generates chimeric RNA-DNA molecules in situ and then converts those chimeras to DNA for next-generation sequencing. Using ChAR-seq we detect many types of chromatin-associated RNA, both coding and non-coding. Understanding the RNA-DNA interactome and its changes during differentiation or disease with ChAR-seq will likely provide key insights into chromatin and RNA biology.

Chromatin-Associated RNA Sequencing (ChAR-seq).

RNA is a fundamental component of chromatin. Noncoding RNAs (ncRNAs) can associate with chromatin to influence gene expression and chromatin state; many also act at long distances from their transcriptional origin. Yet we know almost nothing about the functions or sites of action for most ncRNAs. Current methods to identify sites of RNA interaction with the genome are limited to the study of a single RNA at a time. Here we describe a protocol for ChAR-seq, a strategy to identify all chromatin-associated RNAs and map their DNA contacts genome-wide. In ChAR-seq, proximity ligation of RNA and DNA to a linker molecule is used to construct a chimeric RNA-DNA molecule that is converted to DNA for sequencing. In a single assay, ChAR-seq can discover de novo chromatin interactions of distinct RNAs, including nascent transcripts, splicing RNAs, and long noncoding RNAs (lncRNAs). Resulting "maps" of genome-bound RNAs should provide new insights into RNA biology. © 2019 by John Wiley & Sons, Inc.

Measurement of Mesoscale Conformational Dynamics of Freely Diffusing Molecules with Tracking FCS.

Few techniques are suited to probe the structure and dynamics of molecular complexes at the mesoscale level (∼100-1000 nm). We have developed a single-molecule technique that uses tracking fluorescence correlation spectroscopy (tFCS) to probe the conformation and dynamics of mesoscale molecular assemblies. tFCS measures the distance fluctuations between two fluorescently labeled sites within an untethered, freely diffusing biomolecule. To achieve subdiffraction spatial resolution, we developed a feedback scheme that allows us to maintain the molecule at an optimal position within the laser intensity gradient for fluorescence correlation spectroscopy. We characterized tFCS spatial sensitivity by measuring the Brownian end-to-end dynamics of DNA molecules as short as 1000 bp. We demonstrate that tFCS detects changes in the compaction of reconstituted nucleosome arrays and can assay transient protein-mediated interactions between distant sites in an individual DNA molecule. Our measurements highlight the applicability of tFCS to a wide variety of biochemical processes involving mesoscale conformational dynamics.

Dentate gyrus mossy cells control spontaneous convulsive seizures and spatial memory.

Temporal lobe epilepsy (TLE) is characterized by debilitating, recurring seizures and an increased risk for cognitive deficits. Mossy cells (MCs) are key neurons in the hippocampal excitatory circuit, and the partial loss of MCs is a major hallmark of TLE. We investigated how MCs contribute to spontaneous ictal activity and to spatial contextual memory in a mouse model of TLE with hippocampal sclerosis, using a combination of optogenetic, electrophysiological, and behavioral approaches. In chronically epileptic mice, real-time optogenetic modulation of MCs during spontaneous hippocampal seizures controlled the progression of activity from an electrographic to convulsive seizure. Decreased MC activity is sufficient to impede encoding of spatial context, recapitulating observed cognitive deficits in chronically epileptic mice.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway.

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg2+ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

Quantitative tests of a reconstitution model for RNA folding thermodynamics and kinetics.

Decades of study of the architecture and function of structured RNAs have led to the perspective that RNA tertiary structure is modular, made of locally stable domains that retain their structure across RNAs. We formalize a hypothesis inspired by this modularity-that RNA folding thermodynamics and kinetics can be quantitatively predicted from separable energetic contributions of the individual components of a complex RNA. This reconstitution hypothesis considers RNA tertiary folding in terms of ΔGalign, the probability of aligning tertiary contact partners, and ΔGtert, the favorable energetic contribution from the formation of tertiary contacts in an aligned state. This hypothesis predicts that changes in the alignment of tertiary contacts from different connecting helices and junctions (ΔGHJH) or from changes in the electrostatic environment (ΔG+/-) will not affect the energetic perturbation from a mutation in a tertiary contact (ΔΔGtert). Consistent with these predictions, single-molecule FRET measurements of folding of model RNAs revealed constant ΔΔGtert values for mutations in a tertiary contact embedded in different structural contexts and under different electrostatic conditions. The kinetic effects of these mutations provide further support for modular behavior of RNA elements and suggest that tertiary mutations may be used to identify rate-limiting steps and dissect folding and assembly pathways for complex RNAs. Overall, our model and results are foundational for a predictive understanding of RNA folding that will allow manipulation of RNA folding thermodynamics and kinetics. Conversely, the approaches herein can identify cases where an independent, additive model cannot be applied and so require additional investigation.

Xenopus laevis M18BP1 Directly Binds Existing CENP-A Nucleosomes to Promote Centromeric Chromatin Assembly.

Vertebrate centromeres are epigenetically defined by nucleosomes containing the histone H3 variant, CENP-A. CENP-A nucleosome assembly requires the three-protein Mis18 complex (Mis18α, Mis18ß, and M18BP1) that recruits the CENP-A chaperone HJURP to centromeres, but how the Mis18 complex recognizes centromeric chromatin is unknown. Using Xenopus egg extract, we show that direct, cell-cycle-regulated binding of M18BP1 to CENP-A nucleosomes recruits the Mis18 complex to interphase centromeres to promote new CENP-A nucleosome assembly. We demonstrate that Xenopus M18BP1 binds CENP-A nucleosomes using a motif that is widely conserved except in mammals. The M18BP1 motif resembles a CENP-A nucleosome binding motif in CENP-C, and we show that CENP-C competes with M18BP1 for CENP-A nucleosome binding at centromeres. We show that both CENP-C and M18BP1 recruit HJURP to centromeres for new CENP-A assembly. This study defines cellular mechanisms for recruiting CENP-A assembly factors to existing CENP-A nucleosomes for the epigenetic inheritance of centromeres.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway.

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

Protein flexibility is required for vesicle tethering at the Golgi.

The Golgi is decorated with coiled-coil proteins that may extend long distances to help vesicles find their targets. GCC185 is a trans Golgi-associated protein that captures vesicles inbound from late endosomes. Although predicted to be relatively rigid and highly extended, we show that flexibility in a central region is required for GCC185's ability to function in a vesicle tethering cycle. Proximity ligation experiments show that that GCC185's N-and C-termini are within <40 nm of each other on the Golgi. In physiological buffers without fixatives, atomic force microscopy reveals that GCC185 is shorter than predicted, and its flexibility is due to a central bubble that represents local unwinding of specific sequences. Moreover, 85% of the N-termini are splayed, and the splayed N-terminus can capture transport vesicles in vitro. These unexpected features support a model in which GCC185 collapses onto the Golgi surface, perhaps by binding to Rab GTPases, to mediate vesicle tethering.

The structure of DNA overstretched from the 5'5' ends differs from the structure of DNA overstretched from the 3'3' ends.

It has been suggested that the structure that results when double-stranded DNA (dsDNA) is pulled from the 3'3' ends differs from that which results when it is pulled from the 5'5' ends. In this work, we demonstrate, using lambda phage dsDNA, that the overstretched states do indeed show different properties, suggesting that they correspond to different structures. For 3'3' pulling versus 5'5' pulling, the following differences are observed: (i) the forces at which half of the molecules in the ensemble have made a complete force-induced transition to single stranded DNA are 141 +/- 3 pN and 122 +/- 4 pN, respectively; (ii) the extension vs. force curve for overstretched DNA has a marked change in slope at 127 +/- 3 pN for 3'3' and 110 +/- 3 pN for 5'5'; (iii) the hysteresis (H) in the extension vs. force curves at 150 mM NaCl is 0.3 +/- 0.8 pN microm for 3'3' versus 13 +/- 8 pN for 5'5'; and (iv) 3'3' and 5'5' molecules show different changes in hysteresis due to interactions with beta-cyclodextrin, a molecule that is known to form stable host-guest complexes with rotated base pairs, and glyoxal that is known to bind stably to unpaired bases. These differences and additional findings are well-accommodated by the corresponding structures predicted on theoretical grounds.