A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Bedside Back to Bench: Building Bridges between Basic and Clinical Genomic Research.

Genome sequencing has revolutionized the diagnosis of genetic diseases. Close collaborations between basic scientists and clinical genomicists are now needed to link genetic variants with disease causation. To facilitate such collaborations, we recommend prioritizing clinically relevant genes for functional studies, developing reference variant-phenotype databases, adopting phenotype description standards, and promoting data sharing.

A micropeptide encoded by a putative long noncoding RNA regulates muscle performance.

Functional micropeptides can be concealed within RNAs that appear to be noncoding. We discovered a conserved micropeptide, which we named myoregulin (MLN), encoded by a skeletal muscle-specific RNA annotated as a putative long noncoding RNA. MLN shares structural and functional similarity with phospholamban (PLN) and sarcolipin (SLN), which inhibit SERCA, the membrane pump that controls muscle relaxation by regulating Ca(2+) uptake into the sarcoplasmic reticulum (SR). MLN interacts directly with SERCA and impedes Ca(2+) uptake into the SR. In contrast to PLN and SLN, which are expressed in cardiac and slow skeletal muscle in mice, MLN is robustly expressed in all skeletal muscle. Genetic deletion of MLN in mice enhances Ca(2+) handling in skeletal muscle and improves exercise performance. These findings identify MLN as an important regulator of skeletal muscle physiology and highlight the possibility that additional micropeptides are encoded in the many RNAs currently annotated as noncoding.

Osmosensor-mediated control of Ca2+ spiking in pollen germination.

Higher plants survive terrestrial water deficiency and fluctuation by arresting cellular activities (dehydration) and resuscitating processes (rehydration). However, how plants monitor water availability during rehydration is unknown. Although increases in hypo-osmolarity-induced cytosolic Ca2+ concentration (HOSCA) have long been postulated to be the mechanism for sensing hypo-osmolarity in rehydration1,2, the molecular basis remains unknown. Because osmolarity triggers membrane tension and the osmosensing specificity of osmosensing channels can only be determined in vivo3-5, these channels have been classified as a subtype of mechanosensors. Here we identify bona fide cell surface hypo-osmosensors in Arabidopsis and find that pollen Ca2+ spiking is controlled directly by water through these hypo-osmosensors-that is, Ca2+ spiking is the second messenger for water status. We developed a functional expression screen in Escherichia coli for hypo-osmosensitive channels and identified OSCA2.1, a member of the hyperosmolarity-gated calcium-permeable channel (OSCA) family of proteins6. We screened single and high-order OSCA mutants, and observed that the osca2.1/osca2.2 double-knockout mutant was impaired in pollen germination and HOSCA. OSCA2.1 and OSCA2.2 function as hypo-osmosensitive Ca2+-permeable channels in planta and in HEK293 cells. Decreasing osmolarity of the medium enhanced pollen Ca2+ oscillations, which were mediated by OSCA2.1 and OSCA2.2 and required for germination. OSCA2.1 and OSCA2.2 convert extracellular water status into Ca2+ spiking in pollen and may serve as essential hypo-osmosensors for tracking rehydration in plants.

A myocardin-adjacent lncRNA balances SRF-dependent gene transcription in the heart.

Myocardin, a potent coactivator of serum response factor (SRF), competes with ternary complex factor (TCF) proteins for SRF binding to balance opposing mitogenic and myogenic gene programs in cardiac and smooth muscle. Here we identify a cardiac lncRNA transcribed adjacent to myocardin, named CARDINAL, which antagonizes SRF-dependent mitogenic gene transcription in the heart. CARDINAL-deficient mice show ectopic TCF/SRF-dependent mitogenic gene expression and decreased cardiac contractility in response to age and ischemic stress. CARDINAL forms a nuclear complex with SRF and inhibits TCF-mediated transactivation of the promitogenic gene c-fos, suggesting CARDINAL functions as an RNA cofactor for SRF in the heart.

Will variants of uncertain significance still exist in 2030?

In 2020, the National Human Genome Research Institute (NHGRI) made ten "bold predictions," including that "the clinical relevance of all encountered genomic variants will be readily predictable, rendering the diagnostic designation 'variant of uncertain significance (VUS)' obsolete." We discuss the prospects for this prediction, arguing that many, if not most, VUS in coding regions will be resolved by 2030. We outline a confluence of recent changes making this possible, especially advances in the standards for variant classification that better leverage diverse types of evidence, improvements in computational variant effect predictor performance, scalable multiplexed assays of variant effect capable of saturating the genome, and data-sharing efforts that will maximize the information gained from each new individual sequenced and variant interpreted. We suggest that clinicians and researchers can realize a future where VUSs have largely been eliminated, in line with the NHGRI's bold prediction. The length of time taken to reach this future, and thus whether we are able to achieve the goal of largely eliminating VUSs by 2030, is largely a consequence of the choices made now and in the next few years. We believe that investing in eliminating VUSs is worthwhile, since their predominance remains one of the biggest challenges to precision genomic medicine.

A Combined Approach Reveals a Regulatory Mechanism Coupling Src's Kinase Activity, Localization, and Phosphotransferase-Independent Functions.

Multiple layers of regulation modulate the activity and localization of protein kinases. However, many details of kinase regulation remain incompletely understood. Here, we apply saturation mutagenesis and a chemical genetic method for allosterically modulating kinase global conformation to Src kinase, providing insight into known regulatory mechanisms and revealing a previously undiscovered interaction between Src's SH4 and catalytic domains. Abrogation of this interaction increased phosphotransferase activity, promoted membrane association, and provoked phosphotransferase-independent alterations in cell morphology. Thus, Src's SH4 domain serves as an intramolecular regulator coupling catalytic activity, global conformation, and localization, as well as mediating a phosphotransferase-independent function. Sequence conservation suggests that the SH4 domain regulatory interaction exists in other Src-family kinases. Our combined approach's ability to reveal a regulatory mechanism in one of the best-studied kinases suggests that it could be applied broadly to provide insight into kinase structure, regulation, and function.

Limits of economy and fidelity for programmable assembly of size-controlled triply periodic polyhedra.

We propose and investigate an extension of the Caspar-Klug symmetry principles for viral capsid assembly to the programmable assembly of size-controlled triply periodic polyhedra, discrete variants of the Primitive, Diamond, and Gyroid cubic minimal surfaces. Inspired by a recent class of programmable DNA origami colloids, we demonstrate that the economy of design in these crystalline assemblies-in terms of the growth of the number of distinct particle species required with the increased size-scale (e.g., periodicity)-is comparable to viral shells. We further test the role of geometric specificity in these assemblies via dynamical assembly simulations, which show that conditions for simultaneously efficient and high-fidelity assembly require an intermediate degree of flexibility of local angles and lengths in programmed assembly. Off-target misassembly occurs via incorporation of a variant of disclination defects, generalized to the case of hyperbolic crystals. The possibility of these topological defects is a direct consequence of the very same symmetry principles that underlie the economical design, exposing a basic tradeoff between design economy and fidelity of programmable, size controlled assembly.

Minimizing higher-order aggregation maximizes iron mobilization by small molecules.

The natural product hinokitiol mobilizes iron across lipid bilayers at low concentrations and restores hemoglobinization in iron transporter protein-deficient systems. But hinokitiol fails to similarly mobilize iron at higher concentrations, limiting its uses in chemical biology and medicine. Here we show that at higher concentrations, hinokitiol3:Fe(III) complexes form large, higher-order aggregates, leading to loss of transmembrane iron mobilization. Guided by this understanding and systematic structure-function studies enabled by modular synthesis, we identified FeM-1269, which minimally aggregates and dose-dependently mobilizes iron across lipid bilayers even at very high concentrations. In contrast to hinokitiol, FeM-1269 is also well-tolerated in animals at high doses for extended periods of time. In a mouse model of anemia of inflammation, FeM-1269 increases serum iron, transferrin saturation, hemoglobin and hematocrit. This rationally developed iron-mobilizing small molecule has enhanced potential as a molecular prosthetic for understanding and potentially treating iron transporter deficiencies.

Hydrogen peroxide sensor HPCA1 is an LRR receptor kinase in Arabidopsis.

Hydrogen peroxide (H2O2) is a major reactive oxygen species in unicellular and multicellular organisms, and is produced extracellularly in response to external stresses and internal cues1-4. H2O2 enters cells through aquaporin membrane proteins and covalently modifies cytoplasmic proteins to regulate signalling and cellular processes. However, whether sensors for H2O2 also exist on the cell surface remains unknown. In plant cells, H2O2 triggers an influx of Ca2+ ions, which is thought to be involved in H2O2 sensing and signalling. Here, by using forward genetic screens based on Ca2+ imaging, we isolated hydrogen-peroxide-induced Ca2+ increases (hpca) mutants in Arabidopsis, and identified HPCA1 as a leucine-rich-repeat receptor kinase belonging to a previously uncharacterized subfamily that features two extra pairs of cysteine residues in the extracellular domain. HPCA1 is localized to the plasma membrane and is activated by H2O2 via covalent modification of extracellular cysteine residues, which leads to autophosphorylation of HPCA1. HPCA1 mediates H2O2-induced activation of Ca2+ channels in guard cells and is required for stomatal closure. Our findings help to identify how the perception of extracellular H2O2 is integrated with responses to various external stresses and internal cues in plants, and have implications for the design of crops with enhanced fitness.

Antigen perception in T cells by long-term Erk and NFAT signaling dynamics.

Immune system threat detection hinges on T cells' ability to perceive varying peptide-major histocompatibility complex (pMHC) antigens. As the Erk and NFAT pathways link T cell receptor engagement to gene regulation, their signaling dynamics may convey information about pMHC inputs. To test this idea, we developed a dual reporter mouse strain and a quantitative imaging assay that, together, enable simultaneous monitoring of Erk and NFAT dynamics in live T cells over day-long timescales as they respond to varying pMHC inputs. Both pathways initially activate uniformly across various pMHC inputs but diverge only over longer (9+ h) timescales, enabling independent encoding of pMHC affinity and dose. These late signaling dynamics are decoded via multiple temporal and combinatorial mechanisms to generate pMHC-specific transcriptional responses. Our findings underscore the importance of long timescale signaling dynamics in antigen perception and establish a framework for understanding T cell responses under diverse contexts.

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Measuring Pharmacogene Variant Function at Scale Using Multiplexed Assays.

As costs of next-generation sequencing decrease, identification of genetic variants has far outpaced our ability to understand their functional consequences. This lack of understanding is a central challenge to a key promise of pharmacogenomics: using genetic information to guide drug selection and dosing. Recently developed multiplexed assays of variant effect enable experimental measurement of the function of thousands of variants simultaneously. Here, we describe multiplexed assays that have been performed on nearly 25,000 variants in eight key pharmacogenes (ADRB2, CYP2C9, CYP2C19, NUDT15, SLCO1B1, TMPT, VKORC1, and the LDLR promoter), discuss advances in experimental design, and explore key challenges that must be overcome to maximize the utility of multiplexed functional data.

Pacybara: accurate long-read sequencing for barcoded mutagenized allelic libraries.

MOTIVATION: Long-read sequencing technologies, an attractive solution for many applications, often suffer from higher error rates. Alignment of multiple reads can improve base-calling accuracy, but some applications, e.g. sequencing mutagenized libraries where multiple distinct clones differ by one or few variants, require the use of barcodes or unique molecular identifiers. Unfortunately, sequencing errors can interfere with correct barcode identification, and a given barcode sequence may be linked to multiple independent clones within a given library. RESULTS: Here we focus on the target application of sequencing mutagenized libraries in the context of multiplexed assays of variant effects (MAVEs). MAVEs are increasingly used to create comprehensive genotype-phenotype maps that can aid clinical variant interpretation. Many MAVE methods use long-read sequencing of barcoded mutant libraries for accurate association of barcode with genotype. Existing long-read sequencing pipelines do not account for inaccurate sequencing or nonunique barcodes. Here, we describe Pacybara, which handles these issues by clustering long reads based on the similarities of (error-prone) barcodes while also detecting barcodes that have been associated with multiple genotypes. Pacybara also detects recombinant (chimeric) clones and reduces false positive indel calls. In three example applications, we show that Pacybara identifies and correctly resolves these issues. AVAILABILITY AND IMPLEMENTATION: Pacybara, freely available at https://github.com/rothlab/pacybara, is implemented using R, Python, and bash for Linux. It runs on GNU/Linux HPC clusters via Slurm, PBS, or GridEngine schedulers. A single-machine simplex version is also available.

A chemically controlled Cas9 switch enables temporal modulation of diverse effectors.

CRISPR-Cas9 has yielded a plethora of effectors, including targeted transcriptional activators, base editors and prime editors. Current approaches for inducibly modulating Cas9 activity lack temporal precision and require extensive screening and optimization. We describe a versatile, chemically controlled and rapidly activated single-component DNA-binding Cas9 switch, ciCas9, which we use to confer temporal control over seven Cas9 effectors, including two cytidine base editors, two adenine base editors, a dual base editor, a prime editor and a transcriptional activator. Using these temporally controlled effectors, we analyze base editing kinetics, showing that editing occurs within hours and that rapid early editing of nucleotides predicts eventual editing magnitude. We also reveal that editing at preferred nucleotides within target sites increases the frequency of bystander edits. Thus, the ciCas9 switch offers a simple, versatile approach to generating chemically controlled Cas9 effectors, informing future effector engineering and enabling precise temporal effector control for kinetic studies.

Functional spectrum and specificity of mitochondrial ferredoxins FDX1 and FDX2.

Ferredoxins comprise a large family of iron-sulfur (Fe-S) proteins that shuttle electrons in diverse biological processes. Human mitochondria contain two isoforms of [2Fe-2S] ferredoxins, FDX1 (aka adrenodoxin) and FDX2, with known functions in cytochrome P450-dependent steroid transformations and Fe-S protein biogenesis. Here, we show that only FDX2, but not FDX1, is involved in Fe-S protein maturation. Vice versa, FDX1 is specific not only for steroidogenesis, but also for heme a and lipoyl cofactor biosyntheses. In the latter pathway, FDX1 provides electrons to kickstart the radical chain reaction catalyzed by lipoyl synthase. We also identified lipoylation as a target of the toxic antitumor copper ionophore elesclomol. Finally, the striking target specificity of each ferredoxin was assigned to small conserved sequence motifs. Swapping these motifs changed the target specificity of these electron donors. Together, our findings identify new biochemical tasks of mitochondrial ferredoxins and provide structural insights into their functional specificity.

Carrier-resolved photo-Hall effect.

The fundamental parameters of majority and minority charge carriers-including their type, density and mobility-govern the performance of semiconductor devices yet can be difficult to measure. Although the Hall measurement technique is currently the standard for extracting the properties of majority carriers, those of minority carriers have typically only been accessible through the application of separate techniques. Here we demonstrate an extension to the classic Hall measurement-a carrier-resolved photo-Hall technique-that enables us to simultaneously obtain the mobility and concentration of both majority and minority carriers, as well as the recombination lifetime, diffusion length and recombination coefficient. This is enabled by advances in a.c.-field Hall measurement using a rotating parallel dipole line system and an equation, ΔµH = d(σ2H)/dσ, which relates the hole-electron Hall mobility difference (ΔµH), the conductivity (σ) and the Hall coefficient (H). We apply this technique to various solar absorbers-including high-performance lead-iodide-based perovskites-and demonstrate simultaneous access to majority and minority carrier parameters and map the results against varying light intensities. This information, which is buried within the photo-Hall measurement1,2, had remained inaccessible since the original discovery of the Hall effect in 18793. The simultaneous measurement of majority and minority carriers should have broad applications, including in photovoltaics and other optoelectronic devices.

Plant cell-surface GIPC sphingolipids sense salt to trigger Ca2+ influx.

Salinity is detrimental to plant growth, crop production and food security worldwide. Excess salt triggers increases in cytosolic Ca2+ concentration, which activate Ca2+-binding proteins and upregulate the Na+/H+ antiporter in order to remove Na+. Salt-induced increases in Ca2+ have long been thought to be involved in the detection of salt stress, but the molecular components of the sensing machinery remain unknown. Here, using Ca2+-imaging-based forward genetic screens, we isolated the Arabidopsis thaliana mutant monocation-induced [Ca2+]i increases 1 (moca1), and identified MOCA1 as a glucuronosyltransferase for glycosyl inositol phosphorylceramide (GIPC) sphingolipids in the plasma membrane. MOCA1 is required for salt-induced depolarization of the cell-surface potential, Ca2+ spikes and waves, Na+/H+ antiporter activation, and regulation of growth. Na+ binds to GIPCs to gate Ca2+ influx channels. This salt-sensing mechanism might imply that plasma-membrane lipids are involved in adaption to various environmental salt levels, and could be used to improve salt resistance in crops.

Geometrically programmed self-limited assembly of tubules using DNA origami colloids.

Self-assembly is one of the most promising strategies for making functional materials at the nanoscale, yet new design principles for making self-limiting architectures, rather than spatially unlimited periodic lattice structures, are needed. To address this challenge, we explore the tradeoffs between addressable assembly and self-closing assembly of a specific class of self-limiting structures: cylindrical tubules. We make triangular subunits using DNA origami that have specific, valence-limited interactions and designed binding angles, and we study their assembly into tubules that have a self-limited width that is much larger than the size of an individual subunit. In the simplest case, the tubules are assembled from a single component by geometrically programming the dihedral angles between neighboring subunits. We show that the tubules can reach many micrometers in length and that their average width can be prescribed through the dihedral angles. We find that there is a distribution in the width and the chirality of the tubules, which we rationalize by developing a model that considers the finite bending rigidity of the assembled structure as well as the mechanism of self-closure. Finally, we demonstrate that the distributions of tubules can be further sculpted by increasing the number of subunit species, thereby increasing the assembly complexity, and demonstrate that using two subunit species successfully reduces the number of available end states by half. These results help to shed light on the roles of assembly complexity and geometry in self-limited assembly and could be extended to other self-limiting architectures, such as shells, toroids, or triply periodic frameworks.

FDX1 regulates cellular protein lipoylation through direct binding to LIAS.

Ferredoxins are a family of iron-sulfur (Fe-S) cluster proteins that serve as essential electron donors in numerous cellular processes that are conserved through evolution. The promiscuous nature of ferredoxins as electron donors enables them to participate in many metabolic processes including steroid, heme, vitamin D, and Fe-S cluster biosynthesis in different organisms. However, the unique natural function(s) of each of the two human ferredoxins (FDX1 and FDX2) are still poorly characterized. We recently reported that FDX1 is both a crucial regulator of copper ionophore-induced cell death and serves as an upstream regulator of cellular protein lipoylation, a mitochondrial lipid-based post-translational modification naturally occurring on four mitochondrial enzymes that are crucial for TCA cycle function. Here we show that FDX1 directly regulates protein lipoylation by binding the lipoyl synthase (LIAS) enzyme promoting its functional binding to the lipoyl carrier protein GCSH and not through indirect regulation of cellular Fe-S cluster biosynthesis. Metabolite profiling revealed that the predominant cellular metabolic outcome of FDX1 loss of function is manifested through the regulation of the four lipoylation-dependent enzymes ultimately resulting in loss of cellular respiration and sensitivity to mild glucose starvation. Transcriptional profiling established that FDX1 loss-of-function results in the induction of both compensatory metabolism-related genes and the integrated stress response, consistent with our findings that FDX1 loss-of-function is conditionally lethal. Together, our findings establish that FDX1 directly engages with LIAS, promoting its role in cellular protein lipoylation, a process essential in maintaining cell viability under low glucose conditions.