Photosynthesis tunes quantum-mechanical mixing of electronic and vibrational states to steer exciton energy transfer

Photosynthetic species evolved to protect their light-harvesting apparatus from photoxidative damage driven by intracellular redox conditions or environmental conditions. The Fenna–Matthews–Olson (FMO) pigment–protein complex from green sulfur bacteria exhibits redox-dependent quenching behavior partially due to two internal cysteine residues. Here, we show evidence that a photosynthetic complex exploits the quantum mechanics of vibronic mixing to activate an oxidative photoprotective mechanism. We use two-dimensional electronic spectroscopy (2DES) to capture energy transfer dynamics in wild-type and cysteine-deficient FMO mutant proteins under both reducing and oxidizing conditions. Under reducing conditions, we find equal energy transfer through the exciton 4–1 and 4–2-1 pathways because the exciton 4–1 energy gap is vibronically coupled with a bacteriochlorophyll-a vibrational mode. Under oxidizing conditions, however, the resonance of the exciton 4–1 energy gap is detuned from the vibrational mode, causing excitons to preferentially steer through the indirect 4–2-1 pathway to increase the likelihood of exciton quenching. We use a Redfield model to show that the complex achieves this effect by tuning the site III energy via the redox state of its internal cysteine residues. This result shows how pigment–protein complexes exploit the quantum mechanics of vibronic coupling to steer energy transfer.

Photosynthetic organisms convert solar photons into chemical energy by taking advantage of the quantum mechanical nature of their molecular systems and the chemistry of their environment (1–4). Antenna complexes, composed of one or more pigment–protein complexes, facilitate the first steps in the photosynthesis process: They absorb photons and determine which proportion of excitations to move to reaction centers, where charge separation occurs (4). In oxic environments, excitations can generate highly reactive singlet oxygen species. These pigment–protein complexes can quench excess excitations in these environments with molecular moieties such as quinones and cysteine residues (1, 5–7).The Fenna–Matthews–Olson (FMO) complex, a trimer of pigment–protein complexes found in the green sulfur bacterium Chlorobaculum tepidum (8), has emerged as a model system to study the photophysical properties of photosynthetic antenna complexes (9–19). Each subunit in the FMO complex contains eight bacteriochlorophyll-a site molecules (Protein Data Bank, ID code: 3ENI) that are coupled to form a basis of eight partially delocalized excited states called excitons (Fig. 1) (20–23). Previous experiments on FMO have observed the presence of long-lived coherences in nonlinear spectroscopic signals at both cryogenic and physiological temperatures (11, 13). The coherent signals are thought to arise from some combination of electronic (24–26), vibrational (16–18), and vibronic (27) coherences in the system (28–30). One previous study reported that the coherent signals in FMO remain unchanged upon mutagenesis of the protein, suggesting that the signals are ground state vibrational coherences (17). Others discuss the role of vibronic coupling, where electronic and nuclear degrees of freedom become coupled (29). Other dimeric model systems have demonstrated the regimes in which these vibronically coupled states produce coherent or incoherent transport and vibronic coherences (31–33). Recent spectroscopic data has suggested that vibronic coupling plays a role in driving efficient energy transfer through photosynthetic complexes (27, 31, 33, 34), but to date there is no direct experimental evidence suggesting that biological systems use vibronic coupling as part of their biological function.Open in a separate windowFig. 1.(Left) Numbered sites and sidechains of cysteines C353 and C49 in the FMO pigment–protein complex (PDB ID code: 3ENI) (20). (Right) Site densities for excitons 4, 2, and 1 in reducing conditions with the energy transfer branching ratios for the WT oxidized and reduced protein. The saturation of pigments in each exciton denotes the relative contribution number to the exciton. The C353 residue is located near excitons 4 and 2, which have most electron density along one side of the complex, and other redox-active residues such as the Trp/Tyr chain. C353 and C49 surround site III, which contains the majority of exciton 1 density. Excitons 2 and 4 are generally delocalized over sites IV, V, and VII.It has been shown that redox conditions affect excited state properties in pigment-protein complexes, yet little is known about the underlying microscopic mechanisms for these effects (1, 9). Many commonly studied light-harvesting complexes—including the FMO complex (20), light-harvesting complex 2 (LH2) (35), the PC645 phycobiliprotein (36), and the cyanobacterial antenna complex isiA (37)—contain redox-active cysteine residues in close proximity to their chromophores. As the natural low light environment of C. tepidum does not necessitate photoprotective responses to light quantity and quality, its primary photoprotective mechanism concerns its response to oxidative stress. C. tepidum is an obligate anaerobe, but the presence of many active anoxygenic genes such as sodB for superoxide dismutase and roo for rubredoxin oxygen oxidoreductase (38) suggests that it is frequently exposed to molecular oxygen (7, 39). Using time-resolved fluorescence measurements, Orf et al. demonstrated that two cysteine residues in the FMO complex, C49 and C353, quench excitons under oxidizing conditions (1), which could protect the excitation from generating reactive oxygen species (7, 40–42). In two-dimensional electronic spectroscopy (2DES) experiments, Allodi et al. showed that redox conditions in both the wild-type and C49A/C353A double-mutant proteins affect the ultrafast dynamics through the FMO complex (9, 43). The recent discovery that many proteins across the evolutionary landscape possess chains of tryptophan and tyrosine residues provides evidence that these redox-active residues may link the internal protein behavior with the chemistry of the surrounding environment (41, 43).In this paper, we present data showing that pigment–protein complexes tune the vibronic coupling of their chromophores and that the absence of this vibronic coupling activates an oxidative photoprotective mechanism. We use 2DES to show that a pair of cysteine residues in FMO, C49 and C353, can steer excitations toward quenching sites in oxic environments. The measured reaction rate constants demonstrate unusual nonmonotonic behavior. We then use a Redfield model to determine how the exciton energy transfer (EET) time constants arise from changing chlorophyll site energies and their system-bath couplings (44, 45). The analysis reveals that the cysteine residues tune the resonance between exciton 4–1 energy gap and an intramolecular chlorophyll vibration in reducing conditions to induce vibronic coupling and detune the resonance in oxidizing conditions. This redox-dependent modulation of the vibronic coupling steers excitations through different pathways in the complex to change the likelihood that they interact with exciton quenchers.     

From the Cover: Evidence from South Africa for a protracted end-Permian extinction on land

Earth’s largest biotic crisis occurred during the Permo–Triassic Transition (PTT). On land, this event witnessed a turnover from synapsid- to archosauromorph-dominated assemblages and a restructuring of terrestrial ecosystems. However, understanding extinction patterns has been limited by a lack of high-precision fossil occurrence data to resolve events on submillion-year timescales. We analyzed a unique database of 588 fossil tetrapod specimens from South Africa’s Karoo Basin, spanning ∼4 My, and 13 stratigraphic bin intervals averaging 300,000 y each. Using sample-standardized methods, we characterized faunal assemblage dynamics during the PTT. High regional extinction rates occurred through a protracted interval of ∼1 Ma, initially co-occurring with low origination rates. This resulted in declining diversity up to the acme of extinction near the Daptocephalus–Lystrosaurus declivis Assemblage Zone boundary. Regional origination rates increased abruptly above this boundary, co-occurring with high extinction rates to drive rapid turnover and an assemblage of short-lived species symptomatic of ecosystem instability. The “disaster taxon” Lystrosaurus shows a long-term trend of increasing abundance initiated in the latest Permian. Lystrosaurus comprised 54% of all specimens by the onset of mass extinction and 70% in the extinction aftermath. This early Lystrosaurus abundance suggests its expansion was facilitated by environmental changes rather than by ecological opportunity following the extinctions of other species as commonly assumed for disaster taxa. Our findings conservatively place the Karoo extinction interval closer in time, but not coeval with, the more rapid marine event and reveal key differences between the PTT extinctions on land and in the oceans.

Mass extinctions are major perturbations of the biosphere resulting from a wide range of different causes including glaciations and sea level fall (1), large igneous provinces (2), and bolide impacts (3, 4). These events caused permanent changes to Earth’s ecosystems, altering the evolutionary trajectory of life (5). However, links between the broad causal factors of mass extinctions and the biological and ecological disturbances that lead to species extinctions have been difficult to characterize. This is because ecological disturbances unfold on timescales much shorter than the typical resolution of paleontological studies (6), particularly in the terrestrial record (6–8). Coarse-resolution studies have demonstrated key mass extinction phenomena including high extinction rates and lineage turnover (7, 9), changes in species richness (10), ecosystem instability (11), and the occurrence of disaster taxa (12). However, finer time resolutions are central to determining the association and relative timings of these effects, their potential causal factors, and their interrelationships. Achieving these goals represents a key advance in understanding the ecological mechanisms of mass extinctions.The end-Permian mass extinction (ca. 251.9 Ma) was Earth’s largest biotic crisis as measured by taxon last occurrences (13–15). Large outpourings from Siberian Trap volcanism (2) are the likely trigger of calamitous climatic changes, including a runaway greenhouse effect and ocean acidification, which had profound consequences for life on land and in the oceans (16–18). An estimated 81% of marine species (19) and 89% of tetrapod genera became extinct as established Permian ecosystems gave way to those of the Triassic. In the ocean, this included the complete extinction of reef-forming tabulate and rugose corals (20, 21) and significant losses in previously diverse ammonoid, brachiopod, and crinoid families (22). On land, many nonmammalian synapsids became extinct (16), and the glossopterid-dominated floras of Gondwana also disappeared (23). Stratigraphic sequences document a global “coral gap” and “coal gap” (24, 25), suggesting reef and forest ecosystems were rare or absent for up to 5 My after the event (26). Continuous fossil-bearing deposits documenting patterns of turnover across the Permian–Triassic transition (PTT) on land (27) and in the oceans (28) are geographically widespread (29, 30), including marine and continental successions that are known from China (31, 32) and India (33). Continental successions are known from Russia (34), Australia (35), Antarctica (36), and South Africa’s Karoo Basin (Fig. 1 and 37–40), the latter providing arguably the most densely sampled and taxonomically scrutinized (41–43) continental record of the PTT. The main extinction has been proposed to occur at the boundary between two biostratigraphic zones with distinctive faunal assemblages, the Daptocephalus and Lystrosaurus declivis assemblage zones (Fig. 1), which marks the traditional placement of the Permian–Triassic geologic boundary [(37) but see ref. 44]. Considerable research has attempted to understand the anatomy of the PTT in South Africa (38, 39, 45–52) and to place it in the context of biodiversity changes across southern Gondwana (53, 54) and globally (29, 31, 32, 44, 47, 55).Open in a separate windowFig. 1.Map of South Africa depicting the distribution of the four tetrapod fossil assemblage zones (Cistecephalus, Daptocephalus, Lystrosaurus declivis, Cynognathus) and our two study sites where fossils were collected in this study (sites A and B). Regional lithostratigraphy and biostratigraphy within the study interval are shown alongside isotope dilution–thermal ionization mass spectrometry dates retrieved by Rubidge et al., Botha et al., and Gastaldo et al. (37, 44, 80). The traditional (dashed red line) and associated PTB hypotheses for the Karoo Basin (37, 44) are also shown. Although traditionally associated with the PTB, the Daptocephalus–Lystrosaurus declivis Assemblage Zone boundary is defined by first appearances of co-occurring tetrapod assemblages, so its position relative to the three PTB hypotheses is unchanged. The Ripplemead member (*) has yet to be formalized by the South African Committee for Stratigraphy.Decades of research have demonstrated the richness of South Africa’s Karoo Basin fossil record, resulting in hundreds of stratigraphically well-documented tetrapod fossils across the PTT (37, 39, 56). This wealth of data has been used qualitatively to identify three extinction phases and an apparent early postextinction recovery phase (39, 45, 51). Furthermore, studies of Karoo community structure and function have elucidated the potential role of the extinction and subsequent recovery in breaking the incumbency of previously dominant clades, including synapsids (11, 57). Nevertheless, understanding patterns of faunal turnover and recovery during the PTT has been limited by the scarcity of quantitative investigations. Previous quantitative studies used coarsely sampled data (i.e., assemblage zone scale, 2 to 3 Ma time intervals) to identify low species richness immediately after the main extinction, potentially associated with multiple “boom and bust” cycles of primary productivity based on δ¹³C variation during the first 5 My of the Triassic (41, 58). However, many details of faunal dynamics in this interval remain unknown. Here, we investigate the dynamics of this major tetrapod extinction at an unprecedented time resolution (on the order of hundreds of thousands of years), using sample-standardized methods to quantify multiple aspects of regional change across the Cistecephalus, Daptocephalus, and Lystrosaurus declivis assemblage zones.     

Cholinergic suppression of hippocampal sharp-wave ripples impairs working memory

Learning and memory are assumed to be supported by mechanisms that involve cholinergic transmission and hippocampal theta. Using G protein–coupled receptor-activation–based acetylcholine sensor (GRAB_ACh3.0) with a fiber-photometric fluorescence readout in mice, we found that cholinergic signaling in the hippocampus increased in parallel with theta/gamma power during walking and REM sleep, while ACh3.0 signal reached a minimum during hippocampal sharp-wave ripples (SPW-R). Unexpectedly, memory performance was impaired in a hippocampus-dependent spontaneous alternation task by selective optogenetic stimulation of medial septal cholinergic neurons when the stimulation was applied in the delay area but not in the central (choice) arm of the maze. Parallel with the decreased performance, optogenetic stimulation decreased the incidence of SPW-Rs. These findings suggest that septo–hippocampal interactions play a task-phase–dependent dual role in the maintenance of memory performance, including not only theta mechanisms but also SPW-Rs.

The neurotransmitter acetylcholine is thought to be critical for hippocampus-dependent declarative memories (1, 2). Reduction in cholinergic neurotransmission, either in Alzheimer’s disease or in experiments with cholinergic antagonists, such as scopolamine, impairs memory function (3–8). Acetylcholine may bring about its beneficial effects on memory encoding by enhancing theta rhythm oscillations, decreasing recurrent excitation, and increasing synaptic plasticity (9–11). Conversely, drugs which activate cholinergic receptors enhance learning and, therefore, are a neuropharmacological target for the treatment of memory deficits in Alzheimer’s disease (5, 12, 13).The contribution of cholinergic mechanisms in the acquisition of long-term memories and the role of the hippocampal–entorhinal–cortical interactions are well supported by experimental data (5, 12, 13). In addition, working memory or “short-term” memory is also supported by the hippocampal–entorhinal–prefrontal cortex (14–16). Working memory in humans is postulated to be a conscious process to “keep things in mind” transiently (16). In rodents, matching to sample task, spontaneous alternation between reward locations, and the radial maze task have been suggested to function as a homolog of working memory [“working memory like” (17)].Cholinergic activity is a critical requirement for working memory (18, 19) and for sustaining theta oscillations (10, 20–22). In support of this contention, theta–gamma coupling and gamma power are significantly higher in the choice arm of the maze, compared with those in the side arms where working memory is no longer needed for correct performance (23–26). It has long been hypothesized that working memory is maintained by persistent firing of neurons, which keep the presented items in a transient store in the prefrontal cortex and hippocampal–entorhinal system (27–31), although the exact mechanisms are debated (32–37). An alternative hypothesis holds that items of working memory are stored in theta-nested gamma cycles (38). Common in these models of working memory is the need for an active, cholinergic system–dependent mechanism (39–41). However, in spontaneous alternation tasks, the animals are not moving continuously during the delay, and theta oscillations are not sustained either. During the immobility epochs, theta is replaced by intermittent sharp-wave ripples (SPW-R), yet memory performance does not deteriorate. On the contrary, artificial blockade of SPW-Rs can impair memory performance (42, 43), and prolongation of SPW-Rs improves performance (44). Under the cholinergic hypothesis of working memory, such a result is unexpected.To address the relationship between cholinergic/theta versus SPW-R mechanism in spontaneous alternation, we used a G protein–coupled receptor-activation–based acetylcholine sensor (GRAB_ACh3.0) (45) to monitor acetylcholine (ACh) activity during memory performance in mice. In addition, we optogenetically enhanced cholinergic tone, which suppresses SPW-Rs by a different mechanism than electrically or optogenetically induced silencing of neurons in the hippocampus (43, 44). We show that cholinergic signaling in the hippocampus increases in parallel with theta power/score during walking and rapid eye movement (REM) sleep and reaches a transient minimum during SPW-Rs. Selective optogenetic stimulation of medial septal cholinergic neurons decreased the incidence of SPW-Rs during non-REM sleep (46–48), as well as during the delay epoch of a working memory task and impaired memory performance. These findings demonstrate that memory performance is supported by complementary theta and SPW-R mechanisms.     

Role of intramolecular hydrogen bonds in promoting electron flow through amino acid and oligopeptide conjugates

Elucidating the factors that control charge transfer rates in relatively flexible conjugates is of importance for understanding energy flows in biology as well as assisting the design and construction of electronic devices. Here, we report ultrafast electron transfer (ET) and hole transfer (HT) between a corrole (Cor) donor linked to a perylene-diimide (PDI) acceptor by a tetrameric alanine (Ala)₄. Selective photoexcitation of the donor and acceptor triggers subpicosecond and picosecond ET and HT. Replacement of the (Ala)₄ linker with either a single alanine or phenylalanine does not substantially affect the ET and HT kinetics. We infer that electronic coupling in these reactions is not mediated by tetrapeptide backbone nor by direct donor–acceptor interactions. Employing a combination of NMR, circular dichroism, and computational studies, we show that intramolecular hydrogen bonding brings the donor and the acceptor into proximity in a “scorpion-shaped” molecular architecture, thereby accounting for the unusually high ET and HT rates. Photoinduced charge transfer relies on a (Cor)NH^…O=C–NH^…O=C(PDI) electronic-coupling pathway involving two pivotal hydrogen bonds and a central amide group as a mediator. Our work provides guidelines for construction of effective donor–acceptor assemblies linked by long flexible bridges as well as insights into structural motifs for mediating ET and HT in proteins.

Many vital biological processes, including photosynthesis (1, 2), enzyme catalysis (3–5), and DNA repair (6–9), rely on electron transfer (ET) through proteins and nucleic acids (10–12). As the redox centers and chromophores mediating charge transfer (CT) are normally separated by relatively long molecular distances (>1 nm), the rates of ET and hole transfer (HT) depend critically on the media between donors and acceptors. It follows that the folding of biopolymer backbones plays a key role in governing CT properties.Owing to the complexities of studying ET processes in biomacromolecules, model systems composed of donor and acceptor units located at the termini of oligomeric bridges have been employed to gain insight into the factors that control the redox behavior of biological systems. The evolution of biomimetic and bioinspired molecular devices also has become a major driving force for such studies (13). Despite the inherent flexibility of natural biomolecular oligomers, rigid bridges are common for the myriad of synthetic CT systems (10, 14–19). The difficulty in accounting for conformational dynamics has plagued interpretations of CT through peptides (20–22) and foldamers (23–25).Long-range ET through covalent bonds in donor–bridge–acceptor (DBA) complexes has been the subject of numerous investigations over the last 40 years (10–12, 26). Of particular interest has been work on rigidly linked donor–acceptor units in Ru-modified proteins that has provided experimentally validated tunneling timetables for interpretation and prediction of ET rates in structurally characterized DBA systems (10). Theoretical work by Beratan and others (27) has shed light on the factors that control these couplings both in small molecular DBA complexes and Ru-proteins.Less well understood are couplings through solvents and other noncovalent molecular assemblies, most especially those in DBA complexes that feature B-group hydrogen-bond interactions (28, 29). In work reported here, we employ amino acids and an oligopeptide as bridges in corrole (Cor)–peptide–perylenediimide (PDI) DBA conjugates in order to elucidate the role of hydrogen bonds in mediating biomimetic ET reactions. Specifically, we focus attention on a tetrapeptide bridging group, (Ala)₄, which if extended would separate D and A and greatly diminish through-bond coupling. For comparison purposes, we also have investigated related DBA systems, one with a single alanine (Ala), another with phenylalanine (Phe) (Fig. 1). Varying the number of bridging residues between one and four does not drastically alter the excited-state dynamics of the DBA conjugates. Here, we show that intramolecular hydrogen bonding opens efficient DBA electronic coupling pathways for these CT reactions.Open in a separate windowFig. 1.Structures of DBA conjugates comprising a corrole derivative (Cor) as electron donor, a perrylenediimide derivative (PDI) as electron acceptor, along with the Cor-Phe and PDI used for controls. The hydrogen-bonding donors are highlighted in blue, and the hydrogen-bond acceptors are in red. The upfield-shifted chemical shifts, obtained from ¹H NMR spectra, indicate sections of the bridge (highlighted) that are above the Cor macrocycle plane.     

Developmental influence on evolutionary rates and the origin of placental mammal tooth complexity

Development has often been viewed as a constraining force on morphological adaptation, but its precise influence, especially on evolutionary rates, is poorly understood. Placental mammals provide a classic example of adaptive radiation, but the debate around rate and drivers of early placental evolution remains contentious. A hallmark of early dental evolution in many placental lineages was a transition from a triangular upper molar to a more complex upper molar with a rectangular cusp pattern better specialized for crushing. To examine how development influenced this transition, we simulated dental evolution on “landscapes” built from different parameters of a computational model of tooth morphogenesis. Among the parameters examined, we find that increases in the number of enamel knots, the developmental precursors of the tooth cusps, were primarily influenced by increased self-regulation of the molecular activator (activation), whereas the pattern of knots resulted from changes in both activation and biases in tooth bud growth. In simulations, increased activation facilitated accelerated evolutionary increases in knot number, creating a lateral knot arrangement that evolved at least ten times on placental upper molars. Relatively small increases in activation, superimposed on an ancestral tritubercular molar growth pattern, could recreate key changes leading to a rectangular upper molar cusp pattern. Tinkering with tooth bud geometry varied the way cusps initiated along the posterolingual molar margin, suggesting that small spatial variations in ancestral molar growth may have influenced how placental lineages acquired a hypocone cusp. We suggest that development could have enabled relatively fast higher-level divergence of the placental molar dentition.

Whether developmental processes bias or constrain morphological adaptation is a long-standing question in evolutionary biology (1–4). Many of the distinctive features of a species derive from pattern formation processes that establish the position and number of anatomical structures (5). If developmental processes like pattern formation are biased toward generating only particular kinds of variation, adaptive radiations may often be directed along developmental–genetic “lines of least resistance” (2, 4, 6, 7). Generally, the evolutionary consequences of this developmental bias have been considered largely in terms of how it might influence the pattern of character evolution (e.g., refs. 1, 2, 8–10). But development could also influence evolutionary rates by controlling how much variation is accessible to natural selection in a given generation (11).For mammals, the dentition is often the only morphological system linking living and extinct species (12). Correspondingly, tooth morphology plays a crucial role in elucidating evolutionary relationships, time calibrating phylogenetic trees, and reconstructing adaptive responses to past environmental change (e.g., refs. 13–15). One of the most pervasive features of dental evolution among mammals is an increase in the complexity of the tooth occlusal surface, primarily through the addition of new tooth cusps (16, 17). These increases in tooth complexity are functionally and ecologically significant because they enable more efficient mechanical breakdown of lower-quality foods like plant leaves (18).Placental mammals are the most diverse extant mammalian group, comprising more than 6,000 living species spread across 19 extant orders, and this taxonomic diversity is reflected in their range of tooth shapes and dietary ecologies (12). Many extant placental orders, especially those with omnivorous or herbivorous ecologies (e.g., artiodactyls, proboscideans, rodents, and primates), convergently evolved a rectangular upper molar cusp pattern from a placental ancestor with a more triangular cusp pattern (19–21). This resulted from separate additions in each lineage of a novel posterolingual cusp, the "hypocone'''' [sensu (19)], to the tritubercular upper molar (Fig. 1), either through modification of a posterolingual cingulum (“true” hypocone) or another posterolingual structure, like a metaconule (pseudohypocone) (19). The fossil record suggests that many of the basic steps in the origin of this rectangular cusp pattern occurred during an enigmatic early diversification window associated with the divergence and early radiation of several placental orders (20, 21; Fig. 1). However, there remains debate about the rate and pattern of early placental divergence (22–24). On the one hand, most molecular phylogenies suggest that higher-level placental divergence occurred largely during the Late Cretaceous (25, 26), whereas other molecular phylogenies and paleontological analyses suggest more rapid divergence near the Cretaceous–Paleogene (K–Pg) boundary (21, 24, 27–29). Most studies agree that ecological opportunity created in the aftermath of the K–Pg extinction probably played an important role in ecomorphological diversification within the placental orders (30, 31). But exactly how early placentals acquired the innovations needed to capitalize on ecological opportunity remains unclear. Dental innovations, especially those which facilitated increases in tooth complexity, may have been important because they would have promoted expansion into plant-based dietary ecologies left largely vacant after the K–Pg extinction event (32).Open in a separate windowFig. 1.Placental mammal lineages separately evolved complex upper molar teeth with a rectangular cusp pattern composed of two lateral pairs of cusps from a common ancestor with a simpler, triangular cusp pattern. Many early relatives of the extant placental orders, such as Eritherium, possessed a hypocone cusp and a more rectangular primary cusp pattern. Examples of complex upper molars are the following: Proboscidea, the gomphothere Anancus; Rodentia, the wood mouse Apodemus; and Artiodactyla, the suid Nyanzachoerus.Mammalian tooth cusps form primarily during the “cap” and “bell” stage of dental development, when signaling centers called enamel knots establish the future sites of cusp formation within the inner dental epithelium (33, 34). The enamel knots secrete molecules that promote proliferation and changes in cell–cell adhesion, which facilitates invagination of the dental epithelium into an underlying layer of mesenchymal cells (34, 35). Although a range of genes are involved in tooth cusp patterning (36–38), the basic dynamics can be effectively modeled using reaction–diffusion models with just three diffusible morphogens: an activator, an inhibitor, and a growth factor (39–41). Candidate activator genes in mammalian tooth development include Bmp4, Activin A, Fgf20, and Wnt genes, whereas potential inhibitors include Shh and Sostdc, and Fgf4 and Bmp2 have been hypothesized to act as growth factors (38, 40–43). In computer models of tooth development, activator molecules up-regulated in the underlying mesenchyme stimulate differentiation of overlying epithelium into nondividing enamel knot cells. These in turn secrete molecules that inhibit further differentiation of epithelium into knot cells, while also promoting cell proliferation that creates the topographic relief of the cusp (40). Although many molecular, cellular, and physical processes have the potential to influence cusp formation, and thereby tooth complexity (35, 37), parameters that control the strength and conductance of the activator and inhibitor signals, the core components of the reaction–diffusion cusp patterning mechanism (39, 40) are likely to be especially important.Here, we integrate a previous computer model of tooth morphogenesis called ToothMaker (41), with simulations of trait evolution and data from the fossil record (Fig. 2), to examine the developmental origins of tooth complexity in placental mammals. Specifically, we ask the following: 1) What developmental processes can influence how many cusps form? 2) How might these developmental processes influence the evolution of tooth cusp number, especially rates? And 3) what developmental changes may have been important in the origins of the fourth upper molar cusp, the hypocone, in placental mammal evolution?Open in a separate windowFig. 2.Workflow for simulations of tooth complexity evolution. (A) Tooth shape is varied for five signaling and growth parameters in ToothMaker. (B) From an ancestral state, each parameter is varied in 2.5% increments up to a maximum of ± 50% of the ancestral state. (C) Tooth complexity and enamel knot (EK) pattern were quantified for each parameter combination. Tooth complexity was measured using cusp number/EK number and OPC. ToothMaker and placental upper second molars were classified into categories based on EK/cusp pattern. (D) The parameter space was populated with pattern and tooth complexity datums to build a developmental landscape. (E) Tooth complexity evolution was simulated on each developmental landscape. (F) Resulting diversity and pattern of tooth complexity was compared with placental mammal molar diversity.     

Tropospheric carbonyl sulfide mass balance based on direct measurements of sulfur isotopes

Robust estimates for the rates and trends in terrestrial gross primary production (GPP; plant CO₂ uptake) are needed. Carbonyl sulfide (COS) is the major long-lived sulfur-bearing gas in the atmosphere and a promising proxy for GPP. Large uncertainties in estimating the relative magnitude of the COS sources and sinks limit this approach. Sulfur isotope measurements (³⁴S/³²S; δ³⁴S) have been suggested as a useful tool to constrain COS sources. Yet such measurements are currently scarce for the atmosphere and absent for the marine source and the plant sink, which are two main fluxes. Here we present sulfur isotopes measurements of marine and atmospheric COS, and of plant-uptake fractionation experiments. These measurements resulted in a complete data-based tropospheric COS isotopic mass balance, which allows improved partition of the sources. We found an isotopic (δ³⁴S ± SE) value of 13.9 ± 0.1‰ for the troposphere, with an isotopic seasonal cycle driven by plant uptake. This seasonality agrees with a fractionation of −1.9 ± 0.3‰ which we measured in plant-chamber experiments. Air samples with strong anthropogenic influence indicated an anthropogenic COS isotopic value of 8 ± 1‰. Samples of seawater-equilibrated-air indicate that the marine COS source has an isotopic value of 14.7 ± 1‰. Using our data-based mass balance, we constrained the relative contribution of the two main tropospheric COS sources resulting in 40 ± 17% for the anthropogenic source and 60 ± 20% for the oceanic source. This constraint is important for a better understanding of the global COS budget and its improved use for GPP determination.

The Earth system is going through rapid changes as the climate warms and CO₂ level rises. This rise in CO₂ is mitigated by plant uptake; hence, it is important to estimate global and regional photosynthesis rates and trends (1). Yet, robust tools for investigating these processes at a large scale are scarce (2). Recent studies suggest that carbonyl sulfide (COS) could provide an improved constraint on terrestrial photosynthesis (gross primary production, GPP) (2–12). COS is the major long-lived sulfur-bearing gas in the atmosphere and the main supplier of sulfur to the stratospheric sulfate aerosol layer (13), which exerts a cooling effect on the Earth’s surface and regulates stratospheric ozone chemistry (14).During terrestrial photosynthesis, COS diffuses into leaf stomata and is consumed by photosynthetic enzymes in a similar manner to CO₂ (3–5). Contrary to CO₂, COS undergoes rapid and irreversible hydrolysis mainly by the enzyme carbonic-anhydrase (6, 7). Thus, COS can be used as a proxy for the one-way flux of CO₂ removal from the atmosphere by terrestrial photosynthesis (2, 8–11). However, the large uncertainties in estimating the COS sources weaken this approach (10–12, 15). Tropospheric COS has two main sources: the oceans and anthropogenic emissions, and one main sink–terrestrial plant uptake (8, 10–13). Smaller sources include biomass burning, soil emissions, wetlands, volcanoes, and smaller sinks include OH destruction, stratospheric destruction, and soil uptake (12). The largest source of COS to the atmosphere is the ocean, both as direct COS emission, and as indirect carbon disulfide (CS₂) and dimethylsulfide (DMS) emissions that are rapidly oxidized to COS (10, 16–20). Recent studies suggest oceanic COS emissions are in the range of 200–4,000 GgS/y (19–22). The second major COS source is the anthropogenic source, which is dominated by indirect emissions derived from CS₂ oxidation, mainly from the use of CS₂ as an industrial solvent. Direct emissions of COS are mainly derived from coal and fuel combustion (17, 23, 24). Recent studies suggest that anthropogenic emissions are in the range of 150–585 GgS/y (23, 24). The terrestrial plant uptake is estimated to be in the range of 400–1,360 GgS/y (11). Measurements of sulfur isotope ratios (δ³⁴S) in COS may be used to track COS sources and thus reduce the uncertainties in their flux estimations (15, 25–27). However, the isotopic mass balance approach works best if the COS end members are directly measured and have a significantly different isotopic signature. Previous δ³⁴S measurements of atmospheric COS are scarce and there have been no direct measurements of two important components: the δ³⁴S of oceanic COS emissions, and the isotopic fractionation of COS during plant uptake (15, 25–27). In contrast to previous studies that used assessments for these isotopic values, our aim was to directly measure the isotopic values of these missing components, and to determine the tropospheric COS δ³⁴S variability over space and time.     

Ultrafast X-ray scattering offers a structural view of excited-state charge transfer

Intramolecular charge transfer and the associated changes in molecular structure in N,N′-dimethylpiperazine are tracked using femtosecond gas-phase X-ray scattering. The molecules are optically excited to the 3p state at 200 nm. Following rapid relaxation to the 3s state, distinct charge-localized and charge-delocalized species related by charge transfer are observed. The experiment determines the molecular structure of the two species, with the redistribution of electron density accounted for by a scattering correction factor. The initially dominant charge-localized state has a weakened carbon–carbon bond and reorients one methyl group compared with the ground state. Subsequent charge transfer to the charge-delocalized state elongates the carbon–carbon bond further, creating an extended 1.634 Å bond, and also reorients the second methyl group. At the same time, the bond lengths between the nitrogen and the ring-carbon atoms contract from an average of 1.505 to 1.465 Å. The experiment determines the overall charge transfer time constant for approaching the equilibrium between charge-localized and charge-delocalized species to 3.0 ps.

Understanding the dynamic process of photoinduced charge transfer is expected to lead to many practical applications, including efficient photovoltaic systems, the development of photocatalysts, and better materials for energy storage (1–3). Charge transfer redistributes the electrons in a molecule and is typically associated with changes in the molecular geometry (2). On the fastest time scale, electrons move so rapidly that the nuclei appear frozen, a phenomenon known as charge migration (4–10). When time scales approach the typical vibrational motions of molecules, that is, tens of femtoseconds (10⁻¹⁴ s), the nuclei can adjust their positions, often resulting in localization of electronic charge and permanent changes in molecular geometry (11). While exhibiting a rich phenomenology on different time scales, it is evident that electron charge transfer and nuclear dynamics are intrinsically coupled (12–17). An accurate determination of the changes in molecular structure during charge transfer is therefore of great interest from both applied and fundamental perspectives.New scientific technologies, in particular X-ray free-electron lasers (XFEL) (18, 19) and MeV ultrafast electron diffraction (20), have made it possible to study structural dynamics in the ultrafast regime. Recent femtosecond gas-phase scattering experiments have successfully tracked structural changes during chemical reactions (21–26) and probed specific signatures of excited electronic states (27, 28). Given the emerging ability of ultrafast gas-phase scattering to record both nuclear and electronic structure (25, 27, 28), we use time-resolved gas-phase X-ray scattering to study the photoinduced intramolecular charge transfer in an organic molecule, N,N′-dimethylpiperazine (DMP, C₆H₁₄N₂), shown in Fig. 1. In its ground electronic state, DMP has C_2h symmetry with two equivalent ionization centers, one on each nitrogen atom. Valence ionization, or in the present case excitation to an electronic Rydberg state (29), induces charge transfer between the two nitrogen atoms, making DMP a prototype for exploring electron lone-pair interactions and charge transfer (30–32).Open in a separate windowFig. 1.A schematic illustration of the experimental setup. The ground-state molecules (DMP) were excited by 200 nm ultraviolet pump pulses, and the transient structures were probed by 9.5 keV X-ray pulses at variable time delays. The scattering signals were recorded on a CSPAD detector. The inset shows the calculated spin density, which gives the difference in density of electrons with spin up and spin down, of the charge-localized DMP (3sL) and charge-delocalized DMP (3sD) in the 3s Rydberg states at isovalues of 0.1 electron/ų.Previously, energy relaxation pathways and charge transfer in electronically excited DMP were explored using Rydberg fingerprint spectroscopy (33), a form of photoelectron spectroscopy. As depicted in Scheme 1, the investigations by Deb et al. (33) found that optical excitation at 207 nm prepares the molecule in a 3p Rydberg state, creating a state with a localized charge in the molecular core (3pL). Internal conversion to 3s then leads to charge-localized (3sL) and charge-delocalized (3sD) conformers with 230 and 480 fs time constants, respectively. The charge transfer proceeds as the molecules explore the 3s potential energy surface. An equilibrium between 3sL and 3sD structures is eventually established with an overall time constant of 2.65 ps, with the forward and backward first-order kinetic time constants for the transformation 3.4 and 12.0 ps, respectively (33).Open in a separate windowScheme 1.Reaction pathway for Rydberg-excited DMP as determined previously (33).A limitation in the prior spectroscopic work is that the photoelectron peaks are assigned by comparing measured binding energies with computational results and that the molecular structures that underlie the calculations cannot not be independently determined in the experiments. This is compounded by the fact that theoretical calculation of charge-localized and charge-delocalized excited states is challenging. Widely used computational methods sometimes give unsatisfactory results (34), and even results from high-level computations can be controversial (35–37). This has led to interesting discussions about whether a stable charge-localized structure exists in the DMP cation at all and about the validity of different exchange functionals within density functional theory to study the charge localization in DMP molecules (34). Considering the general experimental and theoretical interest in this system, direct structural measurements would be invaluable. The current study uses ultrafast time-resolved X-ray scattering to observe the structural relaxation dynamics in DMP. This provides a wealth of information, making it possible to test key assumptions of the photoelectron study. Most importantly, we determine the molecular structures of the charge-localized and -delocalized excited states.     

Identifying hydrophobic protein patches to inform protein interaction interfaces

Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces.

Protein–protein interactions play a crucial role in numerous biological processes, ranging from signal transduction and immune response to protein aggregation and phase behavior (1–3). Consequently, being able to understand, predict, and modulate protein interactions has important implications for understanding cellular processes and mitigating the progression of disease (4, 5). A necessary first step toward this ambitious goal is uncovering the interfaces through which proteins interact (6–8). In principle, identifying hydrophobic protein regions, which interact weakly with water, should be a promising strategy for uncovering protein interaction interfaces (9, 10). Indeed, the release of weakly interacting hydration waters from hydrophobic regions can drive protein interactions, as well as other aqueous assemblies (11–13). However, even when the structure of a protein is available at atomistic resolution, it is challenging to identify its hydrophobic patches because they are not uniformly nonpolar, but display variations in polarity and charge at the nanoscale. Moreover, the emergent hydrophobicity of a protein patch stems from the collective response of protein hydration waters to the nanoscale chemical and topographical patterns displayed by the patch (14–20) and cannot be captured by simply counting the number of nonpolar groups in the patch, or even through more involved additive approaches, such as hydropathy scales or surface-area models (21–28).To address this challenge, we build upon seminal theoretical advances and molecular simulation studies, which have shown that near a hydrophobic surface, it is easier to disrupt surface–water interactions and form interfacial cavities (29–34). To uncover protein regions that have the weakest interactions with water, here, we employ specialized molecular simulations, wherein protein–water interactions are disrupted by systematically displacing water molecules from the protein hydration shell (35–37). By identifying the protein patches that nucleate cavities most readily in our simulations, we are then able to uncover the most hydrophobic protein regions. Interestingly, we find that both hydrophobic and hydrophilic protein patches are highly heterogeneous and contain comparable numbers of nonpolar and polar atoms. Our results thus highlight the nontrivial relationship between the chemical composition of protein patches and their emergent hydrophobicity (24–26), and further emphasize the importance of accounting for the collective solvent response in characterizing protein hydrophobicity (16). We then interrogate whether the most hydrophobic protein patches, which nucleate cavities readily, are also likely to mediate protein interactions. Across five proteins that participate in either homodimer or heterodimer formation, we find that roughly 60 to 70% of interfacial contacts and only about 10 to 20% of noncontacts nucleate cavities. Our work thus provides a versatile computational framework for characterizing hydrophobicity and uncovering interaction interfaces of not just proteins, but also of other complex amphiphilic solutes, such as cavitands, dendrimers, and patchy nanoparticles (38–41).     

Seasonal drought events in tropical East Asia over the last 60,000 y

The cause of seasonal hydrologic changes in tropical East Asia during interstadial/stadial oscillations of the last glaciation remains controversial. Here, we show seven seasonal drought events that occurred during the relatively warm interstadials by phytolith and pollen records. These events are significantly manifested as high percentages of bilobate phytoliths and are consistent with the large zonal sea-surface temperature (SST) gradient from the western to eastern tropical Pacific, suggesting that the reduction in seasonal precipitation could be interpreted by westward shifts of the western Pacific subtropical high triggered by changes of zonal SST gradient over the tropical Pacific and Hadley circulation in the Northern Hemisphere. Our findings highlight that both zonal and meridional ocean–atmosphere circulations, rather than solely the Intertropical Convergence Zone or El Niño-Southern Oscillation, controlled the hydrologic changes in tropical East Asia during the last glaciation.

The tropics are home to regions of central importance to global hydrologic cycles. As an important hydrological parameter, seasonal precipitation underwent profound changes in the tropical oceans and their adjacent continents on different timescales during the last glaciation (1). However, the features and mechanisms that caused seasonal precipitation changes during interstadial/stadial oscillations of the last glacial–interglacial cycle in tropical regions remain controversial (2). Two contrasting hypotheses were proposed to explain these changes. Both of these hypotheses involve interactions between the meridional Hadley and zonal Walker circulations as well as changes in sea-surface temperature (SST) (3).The first hypothesis invokes tropical SST changes in response to a slowdown of the Atlantic Meridional Overturning Circulation (AMOC) triggered by large iceberg discharges in the North Atlantic during the stadials (4, 5) and the resulting shift of the intertropical convergence zone (ITCZ). This AMOC-forced southward movement of the ITCZ caused antiphase changes in seasonal rainfall between north and south of the equator in the tropical Pacific and its adjacent continents on the millennium scale (6–9). The second hypothesis attributes the changes in the tropical atmosphere–ocean dynamics stimulated by tropical SST changes (e.g., El Niño-Southern Oscillation, ENSO) (10, 11) as the cause of seasonal heavy/light rainfall in the eastern/western tropical Pacific during the Heinrich stadials (12).However, most of the evidence supporting these two hypotheses comes from the Atlantic (13) and the tropical Pacific (1, 12, 14). Limited evidence is from tropical continents, and particularly little evidence is available from terrestrial tropical East Asia, where the mechanisms responsible for seasonal precipitation changes are still under debate. Although the oxygen isotopes of stalagmites on land have been viewed as one of the most robust East Asian summer monsoon records, the interpretation of speleothem δ¹⁸O in China remains controversial (15–17). Therefore, to demonstrate these hypotheses, both unique research sites, which could simultaneously receive continental and tropical ocean signals, and unambiguous proxies, which could reflect annual and seasonal precipitation, are urgently needed.The core studied in this paper is located at Huguangyan Maar Lake (HML) in Guangdong, southern coast of China (110°17′ E, 21°9′ N, 23 m above sea level) (Fig. 1). This region is closely influenced by Hadley circulation (HC) over the continent and Walker circulation (WC) over the tropical Pacific (18, 19), as well as the related changes in the western Pacific subtropical high (WPSH) (19, 20). Here, we present two sets of proxy datasets. One proxy is the phytolith record, a reliable indicator for the evaluation of seasonal–annual Poaceae in ecosystems (21, 22); and Poaceae taxa are more sensitive than arboreal taxa in revealing seasonal hydrological changes (21, 22). The other proxy is the pollen record, which is an indicator of annual precipitation and reflects the general paleoclimate features at low latitudes in East Asia (23–29). Therefore, based on the two contrasting proxies, the successive 60,000-y phytolith and pollen sequences in this study could substantially reveal the mechanism of regional seasonal precipitation changes during stadials/interstadials since the last glacial period in northern tropical East Asia.Open in a separate windowFig. 1.Geographical and climatological settings of study site and region. (A) Topographic map from ENVI 5.1 showing the location of study site (red star, Huguangyan Maar Lake, HML) in south China and seasonal positions of ITCZ (30). (B) Schematic map showing the settings of atmospheric circulations and SSTs in the tropical Pacific. Black solid circles indicate cores mentioned in this paper: 1 = MD98-2181, 2 = ODP806B, 3 = MD02-2529, 4 = ME0005A-24JC, 5 = ODP846, 6 = TR163-22, and 7 = TR163-19. Colors are annual mean SSTs from the period 1981–2010 obtained from https://psl.noaa.gov.     

Replisome bypass of a protein-based R-loop block by Pif1

Efficient and faithful replication of the genome is essential to maintain genome stability. Replication is carried out by a multiprotein complex called the replisome, which encounters numerous obstacles to its progression. Failure to bypass these obstacles results in genome instability and may facilitate errors leading to disease. Cells use accessory helicases that help the replisome bypass difficult barriers. All eukaryotes contain the accessory helicase Pif1, which tracks in a 5′–3′ direction on single-stranded DNA and plays a role in genome maintenance processes. Here, we reveal a previously unknown role for Pif1 in replication barrier bypass. We use an in vitro reconstituted Saccharomyces cerevisiae replisome to demonstrate that Pif1 enables the replisome to bypass an inactive (i.e., dead) Cas9 (dCas9) R-loop barrier. Interestingly, dCas9 R-loops targeted to either strand are bypassed with similar efficiency. Furthermore, we employed a single-molecule fluorescence visualization technique to show that Pif1 facilitates this bypass by enabling the simultaneous removal of the dCas9 protein and the R-loop. We propose that Pif1 is a general displacement helicase for replication bypass of both R-loops and protein blocks.

Efficient and faithful replication of the genome is essential to maintain genome stability and is carried out by a multiprotein complex called the replisome (1–4). There are numerous obstacles to progression of the replisome during the process of chromosome duplication. These obstacles include RNA-DNA hybrids (R-loops), DNA secondary structures, transcribing RNA polymerases, and other tightly bound proteins (5–9). Failure to bypass these barriers may result in genome instability, which can lead to cellular abnormalities and genetic disease. Cells contain various accessory helicases that help the replisome bypass these difficult barriers (10–20). A subset of these helicases act on the opposite strand of the replicative helicase (1, 2, 14, 19).All eukaryotes contain an accessory helicase, Pif1, which tracks in a 5′–3′ direction on single-stranded DNA (ssDNA) (11–16). Pif1 is important in pathways such as Okazaki-fragment processing and break-induced repair that require the removal of DNA-binding proteins as well as potential displacement of R-loops (11–13, 21, 15–18, 22–25). Genetic studies and immunoprecipitation pull-down assays indicate that Pif1 interacts with PCNA (the DNA sliding clamp), Pol ε (the leading-strand polymerase), the MCMs (the motor subunits of the replicative helicase CMG), and RPA (the single-stranded DNA-binding protein) (15, 26, 27). Pif1 activity in break-induced repair strongly depends on its interaction with PCNA (26). These interactions with replisomal components suggest that Pif1 could interact with the replisome during replication. In Escherichia coli, the replicative helicase is the DnaB homohexamer that encircles the lagging strand and moves in a 5′–3′ direction (20). E. coli accessory helicases include the monomeric UvrD (helicase II) and Rep, which move in the 3′–5′ direction and operate on the opposite strand from the DnaB hexamer. It is known that these monomeric helicases promote the bypass of barriers during replication such as stalled RNA polymerases (5). The eukaryotic replicative helicase is the 11-subunit CMG (Cdc45, Mcm2–7, GINS) and tracks in the 3′–5′ direction, opposite to the direction of Pif1 (25, 28). Once activated by Mcm10, the MCM motor domains of CMG encircle the leading strand (29–32). We hypothesized that, similar to UvrD and Rep in E. coli, Pif1 interacts with the replisome tracking in the opposite direction to enable bypass of replication obstacles.In this report, we use an in vitro reconstituted Saccharomyces cerevisiae replisome to study the role of Pif1 in bypass of a “dead” Cas9 (dCas9), which is a Cas9 protein that is deactivated in DNA cleavage but otherwise fully functional in DNA binding. As with Cas9, dCas9 is a single-turnover enzyme that can be programmed with a guide RNA (gRNA) to target either strand. The dCas9–gRNA complex forms a roadblock consisting of an R-loop and a tightly bound protein (dCas9), a construct that is similar to a stalled RNA polymerase. This roadblock (hereafter dCas9 R-loop) arrests replisomes independent of whether the dCas9 R-loop is targeted to the leading or lagging strand (30). Besides its utility due to its programmable nature (33), the use of the dCas9 R-loop allows us to answer several mechanistic questions. For example, the ability to program the dCas9 R-loop block to any specific sequence enables us to observe whether block removal is different depending on whether the block is on the leading or lagging strand. Furthermore, the inner diameter of CMG can accommodate double-stranded DNA (dsDNA) and possibly an R-loop, but not a dCas9 protein. Using the dCas9 R-loop block allows us to determine the fate of each of its components.Here, we report that Pif1 enables the bypass of the dCas9 R-loop by the replisome. Interestingly, dCas9 R-loops targeted to either the leading or lagging strand are bypassed with similar efficiency. In addition, the PCNA clamp is not required for bypass of the block, indicating that Pif1 does not need to interact with PCNA during bypass of the block. We used a single-molecule fluorescence imaging to show that both the dCas9 and the R-loop are displaced as an intact nucleoprotein complex. We propose that Pif1 is a general displacement helicase for replication bypass of both R-loops and protein blocks.     

Chemokine CCL5 promotes robust optic nerve regeneration and mediates many of the effects of CNTF gene therapy

Ciliary neurotrophic factor (CNTF) is a leading therapeutic candidate for several ocular diseases and induces optic nerve regeneration in animal models. Paradoxically, however, although CNTF gene therapy promotes extensive regeneration, recombinant CNTF (rCNTF) has little effect. Because intraocular viral vectors induce inflammation, and because CNTF is an immune modulator, we investigated whether CNTF gene therapy acts indirectly through other immune mediators. The beneficial effects of CNTF gene therapy remained unchanged after deleting CNTF receptor alpha (CNTFRα) in retinal ganglion cells (RGCs), the projection neurons of the retina, but were diminished by depleting neutrophils or by genetically suppressing monocyte infiltration. CNTF gene therapy increased expression of C-C motif chemokine ligand 5 (CCL5) in immune cells and retinal glia, and recombinant CCL5 induced extensive axon regeneration. Conversely, CRISPR-mediated knockdown of the cognate receptor (CCR5) in RGCs or treating wild-type mice with a CCR5 antagonist repressed the effects of CNTF gene therapy. Thus, CCL5 is a previously unrecognized, potent activator of optic nerve regeneration and mediates many of the effects of CNTF gene therapy.

Like most pathways in the mature central nervous system (CNS), the optic nerve cannot regenerate once damaged due in part to cell-extrinsic suppressors of axon growth (1, 2) and the low intrinsic growth capacity of adult retinal ganglion cells (RGCs), the projection neurons of the eye (3–5). Consequently, traumatic or ischemic optic nerve injury or degenerative diseases such as glaucoma lead to irreversible visual losses. Experimentally, some degree of regeneration can be induced by intraocular inflammation or growth factors expressed by inflammatory cells (6–10), altering the cell-intrinsic growth potential of RGCs (3–5), enhancing physiological activity (11, 12), chelating free zinc (13, 14), and other manipulations (15–19). However, the extent of regeneration achieved to date remains modest, underlining the need for more effective therapies.Ciliary neurotrophic factor (CNTF) is a leading therapeutic candidate for glaucoma and other ocular diseases (20–23). Activation of the downstream signal transduction cascade requires CNTF to bind to CNTF receptor-α (CNTFRα) (24), which leads to recruitment of glycoprotein 130 (gp130) and leukemia inhibitory factor receptor-β (LIFRβ) to form a tripartite receptor complex (25). CNTFRα anchors to the plasma membrane through a glycosylphosphatidylinositol linkage (26) and can be released and become soluble through phospholipase C-mediated cleavage (27). CNTF has been reported to activate STAT3 phosphorylation in retinal neurons, including RGCs, and to promote survival, but it is unknown whether these effects are mediated by direct action of CNTF on RGCs via CNTFRα (28). Our previous studies showed that CNTF promotes axon outgrowth from neonate RGCs in culture (29) but fails to do so in cultured mature RGCs (8) or in vivo (6). Although some studies report that recombinant CNTF (rCNTF) can promote optic nerve regeneration (20, 30, 31), others find little or no effect unless SOCS3 (suppressor of cytokine signaling-3), an inhibitor of the Jak-STAT pathway, is deleted in RGCs (5, 6, 32). In contrast, multiple studies show that adeno-associated virus (AAV)-mediated expression of CNTF in RGCs induces strong regeneration (33–40). The basis for the discrepant effects of rCNTF and CNTF gene therapy is unknown but is of considerable interest in view of the many promising clinical and preclinical outcomes obtained with CNTF to date.Because intravitreal virus injections induce inflammation (41), we investigated the possibility that CNTF, a known immune modulator (42–44), might act by elevating expression of other immune-derived factors. We report here that the beneficial effects of CNTF gene therapy in fact require immune system activation and elevation of C-C motif chemokine ligand 5 (CCL5). Depletion of neutrophils, global knockout (KO) or RGC-selective deletion of the CCL5 receptor CCR5, or a CCR5 antagonist all suppress the effects of CNTF gene therapy, whereas recombinant CCL5 (rCCL5) promotes axon regeneration and increases RGC survival. These studies point to CCL5 as a potent monotherapy for optic nerve regeneration and to the possibility that other applications of CNTF and other forms of gene therapy might similarly act indirectly through other factors.     

A solution to a sex ratio puzzle in Melittobia wasps

The puzzling sex ratio behavior of Melittobia wasps has long posed one of the greatest questions in the field of sex allocation. Laboratory experiments have found that, in contrast to the predictions of theory and the behavior of numerous other organisms, Melittobia females do not produce fewer female-biased offspring sex ratios when more females lay eggs on a patch. We solve this puzzle by showing that, in nature, females of Melittobia australica have a sophisticated sex ratio behavior, in which their strategy also depends on whether they have dispersed from the patch where they emerged. When females have not dispersed, they lay eggs with close relatives, which keeps local mate competition high even with multiple females, and therefore, they are selected to produce consistently female-biased sex ratios. Laboratory experiments mimic these conditions. In contrast, when females disperse, they interact with nonrelatives, and thus adjust their sex ratio depending on the number of females laying eggs. Consequently, females appear to use dispersal status as an indirect cue of relatedness and whether they should adjust their sex ratio in response to the number of females laying eggs on the patch.

Sex allocation has produced many of the greatest success stories in the study of social behaviors (1–4). Time and time again, relatively simple theory has explained variation in how individuals allocate resources to male and female reproduction. Hamilton’s local mate competition (LMC) theory predicts that when n diploid females lay eggs on a patch and the offspring mate before the females disperse, the evolutionary stable proportion of male offspring (sex ratio) is (n − 1)/2n (Fig. 1) (5). A female-biased sex ratio is favored to reduce competition between sons (brothers) for mates and to provide more mates (daughters) for those sons (6–8). Consistent with this prediction, females of >40 species produce female-biased sex ratios and reduce this female bias when multiple females lay eggs on the same patch (higher n; Fig. 1) (9). The fit of data to theory is so good that the sex ratio under LMC has been exploited as a “model trait” to study the factors that can constrain “perfect adaptation” (4, 10–13).Open in a separate windowFig. 1.LMC. The sex ratio (proportion of sons) is plotted versus the number of females laying eggs on a patch. The bright green dashed line shows the LMC theory prediction for the haplodiploid species (5, 39). A more female-biased sex ratio is favored in haplodiploids because inbreeding increases the relative relatedness of mothers to their daughters (7, 32). Females of many species adjust their offspring sex ratio as predicted by theory, such as the parasitoid Nasonia vitripennis (green diamonds) (82). In contrast, the females of several Melittobia species, such as M. australica, continue to produce extremely female-biased sex ratios, irrespective of the number of females laying eggs on a patch (blue squares) (15).In stark contrast, the sex ratio behavior of Melittobia wasps has long been seen as one of the greatest problems for the field of sex allocation (3, 4, 14–21). The life cycle of Melittobia wasps matches the assumptions of Hamilton’s LMC theory (5, 15, 19, 21). Females lay eggs in the larvae or pupae of solitary wasps and bees, and then after emergence, female offspring mate with the short-winged males, who do not disperse. However, laboratory experiments on four Melittobia species have found that females lay extremely female-biased sex ratios (1 to 5% males) and that these extremely female-biased sex ratios change little with increasing number of females laying eggs on a patch (higher n; Fig. 1) (15, 17–20, 22). A number of hypotheses to explain this lack of sex ratio adjustment have been investigated and rejected, including sex ratio distorters, sex differential mortality, asymmetrical male competition, and reciprocal cooperation (15–18, 20, 22–26).We tested whether Melittobia’s unusual sex ratio behavior can be explained by females being related to the other females laying eggs on the same patch. After mating, some females disperse to find new patches, while some may stay at the natal patch to lay eggs on previously unexploited hosts (Fig. 2). If females do not disperse, they can be related to the other females laying eggs on the same host (27–31). If females laying eggs on a host are related, this increases the extent to which relatives are competing for mates and so can favor an even more female-biased sex ratio (28, 32–35). Although most parasitoid species appear unable to directly assess relatedness, dispersal behavior could provide an indirect cue of whether females are with close relatives (36–38). Consequently, we predict that when females do not disperse and so are more likely to be with closer relatives, they should maintain extremely female-biased sex ratios, even when multiple females lay eggs on a patch (28, 35).Open in a separate windowFig. 2.Host nest and dispersal manners of Melittobia. (A) Photograph of the prepupae of the leaf-cutter bee C. sculpturalis nested in a bamboo cane and (B) a diagram showing two ways that Melittobia females find new hosts. The mothers of C. sculpturalis build nursing nests with pine resin consisting of individual cells in which their offspring develop. If Melittobia wasps parasitize a host in a cell, female offspring that mate with males inside the cell find a different host on the same patch (bamboo cane) or disperse by flying to other patches.We tested whether the sex ratio of Melittobia australica can be explained by dispersal status in a natural population. We examined how the sex ratio produced by females varies with the number of females laying eggs on a patch and whether or not they have dispersed before laying eggs. To match our data to the predictions of theory, we developed a mathematical model tailored to the unique population structure of Melittobia, where dispersal can be a cue of relatedness. We then conducted a laboratory experiment to test whether Melittobia females are able to directly access the relatedness to other females and adjust their sex ratio behavior accordingly. Our results suggest that females are adjusting their sex ratio in response to both the number of females laying eggs on a patch and their relatedness to the other females. However, relatedness is assessed indirectly by whether or not they have dispersed. Consequently, the solution to the puzzling behavior reflects a more-refined sex ratio strategy.     

Superelastic oxide micropillars enabled by surface tension–modulated 90° domain switching with excellent fatigue resistance

Superelastic materials capable of recovering large nonlinear strains are ideal for a variety of applications in morphing structures, reconfigurable systems, and robots. However, making oxide materials superelastic has been a long-standing challenge due to their intrinsic brittleness. Here, we fabricate ferroelectric BaTiO₃ (BTO) micropillars that not only are superelastic but also possess excellent fatigue resistance, lasting over 1 million cycles without accumulating residual strains or noticeable variation in stress–strain curves. Phase field simulations reveal that the large recoverable strains of BTO micropillars arise from surface tension–modulated 90° domain switching and thus are size dependent, while the small energy barrier and ultralow energy dissipation are responsible for their unprecedented cyclic stability among superelastic materials. This work demonstrates a general strategy to realize superelastic and fatigue-resistant domain switching in ferroelectric oxides for many potential applications.

Superelastic materials are capable of recovering large amount of nonlinear “plastic” strains, way beyond their linear elastic regimes (1–4). They are ideal for a variety of applications from morphing structures, reconfigurable systems, to robots (5–8). The effects have traditionally been associated with macroscopically compliant/ductile rubbers (2) or microscopically phase-transforming shape memory alloys (SMAs) (7–11). The only macroscopically brittle oxide recently discovered to be superelastic is ZrO₂-based micropillars or particles (12–20), which is realized via austenite-martensite phase transformation similar to SMAs. Although ultimate strengths approaching the theoretical limit have been demonstrated in nanoscale samples (21, 22), long fatigue life is elusive, which is arguably more important for most applications. As a matter of fact, poor fatigue life has been a long-standing challenge for oxide ceramics in general (23, 24). Even for ductile SMAs that enjoy excellent fatigue life, irrecoverable residual strains gradually accumulate over cycling, leading to substantial variations in stress–strain curves at different cycles (9, 10, 25). We overcome these difficulties by reporting superelastic barium titanate (BaTiO₃ [BTO]) micropillars enabled by surface tension–modulated 90° domain switching, which exhibit excellent fatigue resistance, while bulk BTO crystals or ceramics are rather brittle. The demonstration of over one million cycles of loading and unloading without accumulating residual strains or noticeable variation in stress–strain curves is unprecedented among superelastic materials.BTO is a ferroelectric oxide exhibiting modest piezoelectric strains around 0.1 to 0.2% (26) and fracture toughness of ∼1 MPa ⋅ m^1/2, and thus it is quite brittle (27). Considerable research efforts have been devoted to enhancing its electric field–induced strain via 90° ferroelectric domain switching (28–30). However, the process is often irreversible, and external mechanisms such as restoring force (28, 29) and internal mechanisms such as defect pinning (30) have to be invoked to make the electrostrain recoverable. Nevertheless, it hints at the possibility of BTO being made superelastic by taking advantage of the stress-induced 90° domain switching (6). Earlier works suggest that surface tension induces an in-plane compressive stress that favors the axial polarization in one-dimensional ferroelectrics at small size (31, 32), which may provide the necessary restoring mechanism for the stress-switched domains. Thus, if a compressive axial force is applied, reversible domain switching may occur during unloading, leading to superelasticity. To verify this hypothesis, we fabricated single-crystalline BTO micropillars from [001]-oriented bulk crystals (SI Appendix, Fig. S1A) via focused ion beam (FIB), as detailed in Materials and Methods and SI Appendix, Fig. S1B. The diameters (Φ) of the micropillars range from 0.5 μm to 5 μm, with their height to diameter ratio fixed at 3. No visible defects can be seen from the scanning electron microscopy (SEM) images of these micropillars shown in Fig. 1 A–D, and their surfaces appear to be quite smooth, suggesting that no apparent damages are induced by FIB.Open in a separate windowFig. 1.Superelastic BTO micropillars below a critical size. (A–D) SEM images of the micropillars with Φ = 5, 3, 2, and 0.5 μm. (E–G) The first and second cycles of stress–strain curves for BTO micropillars with Φ = 5, 2, and 0.5 μm. (H) S_r/S_max and ΔW/W_max during the first cycle for BTO micropillars of different diameters. Here, S_r and S_max denote the residual strain and the maximum strain (SI Appendix, Fig. S6A), while ΔW and W_max are energy dissipated and stored in the first cycle, respectively (SI Appendix, Fig. S6F).     

Cognitive aging is associated with redistribution of synaptic weights in the hippocampus

Behaviors that rely on the hippocampus are particularly susceptible to chronological aging, with many aged animals (including humans) maintaining cognition at a young adult-like level, but many others the same age showing marked impairments. It is unclear whether the ability to maintain cognition over time is attributable to brain maintenance, sufficient cognitive reserve, compensatory changes in network function, or some combination thereof. While network dysfunction within the hippocampal circuit of aged, learning-impaired animals is well-documented, its neurobiological substrates remain elusive. Here we show that the synaptic architecture of hippocampal regions CA1 and CA3 is maintained in a young adult-like state in aged rats that performed comparably to their young adult counterparts in both trace eyeblink conditioning and Morris water maze learning. In contrast, among learning-impaired, but equally aged rats, we found that a redistribution of synaptic weights amplifies the influence of autoassociational connections among CA3 pyramidal neurons, yet reduces the synaptic input onto these same neurons from the dentate gyrus. Notably, synapses within hippocampal region CA1 showed no group differences regardless of cognitive ability. Taking the data together, we find the imbalanced synaptic weights within hippocampal CA3 provide a substrate that can explain the abnormal firing characteristics of both CA3 and CA1 pyramidal neurons in aged, learning-impaired rats. Furthermore, our work provides some clarity with regard to how some animals cognitively age successfully, while others’ lifespans outlast their “mindspans.”

Aging is the biggest risk factor for Alzheimer’s disease, but many aged individuals nevertheless retain the ability to perform cognitive tasks with young adult (YA)-like competency, and are thus resilient to age-related cognitive decline and dementias (1, 2). The mechanisms of such resilience are unknown, but are thought to involve neural or cognitive reserve, brain network adaptations, or simply the ability to maintain cognitive brain circuits in a YA-like state (3–5). Much of the cellular and functional insight into the concept or risk of/resilience against age-related cognitive impairments has come from animal models of normal/nonpathological aging (6–10). Many of these studies have shown that circuit function abnormalities are associated with behavioral impairments. The cellular and structural bases for such functional aberrations, however, remain largely unknown.Two of the most well-studied cognitive domains that show susceptibility to chronological aging in both rodents and nonhuman primates are working memory and spatial/temporal memory (6–10). Importantly, these cognitive domains engage anatomically distinct neurocognitive systems, with the former relying on prefrontal/orbitofrontal cortical circuits and the latter relying on hippocampal circuitry. Interestingly, although behavioral deficits in these two domains (in the case of rat models of cognitive aging) begin to emerge, worsen, and become increasingly prevalent between 12 and 18 mo of age in most strains (reviewed in refs. 9 and 11), cognitive aging within hippocampus-dependent forms of learning and memory are relatively independent of those that engage the prefrontal/orbitofrontal cortical neural systems (6–9, 12–15).Neither the mechanisms underlying the conservation of memory function across chronological aging nor those contributing to the age-related emergence and exacerbation of memory impairments are clearly understood for either neurocognitive system. It is clear, however, that neither frank neuronal loss (16, 17) nor overall synapse loss (18) contributes to cognitive aging within the medial temporal lobe/hippocampal memory system. Rather, the intriguing idea that has emerged from work in both the hippocampal and the prefrontal/orbitofrontal cortical memory systems is that there are functional alterations in the synaptic connections in individual microcircuits embedded within these larger neuroanatomical systems (6–10, 19–31).Axospinous synapses (including those in hippocampal and cortical circuits) are characterized on the basis of the three-dimensional morphology of their postsynaptic densities (PSDs) (20, 32–34). The most-abundant axospinous synaptic subtype has a continuous, macular, disk-shaped PSD, as compared to the less-abundant perforated synaptic subtype, which has at least one discontinuity in its PSD (34). In addition to differing substantially with regard to relative frequency, perforated and nonperforated synapses also harbor major differences in size and synaptic AMPA-type and NMDA-type receptor expression levels (AMPAR and NMDAR, respectively) (34–38). There is also evidence that perforated and nonperforated synapses are differentially involved in synaptic plasticity (39–44) and in preservation of—or reductions in—memory function during chronological aging (6, 20, 45). Layered onto these general distinctions between perforated and nonperforated synapses are more specific differences in their characteristics when considered within neural circuits. For example, perforated synapses have a stronger and more consistent influence on neuronal computation within hippocampal region CA1 than their nonperforated counterparts, which nevertheless outnumber the former by a roughly 9-to-1 ratio (34, 46, 47).These and other circuit-specific differences necessitate a circuit-based approach to understanding the synaptic bases underlying the retention or loss of YA-like memory function in the aging brain. In many ways, the hippocampal system is particularly convenient for such circuit-based approaches (48, 49). Information about the internal and external world is funneled to the parahippocampal system and then relayed via the entorhinal cortex to the dentate gyrus, the first component of the so-called trisynaptic circuit in the hippocampus proper. Granule cells in the dentate gyrus then transmit their computations to hippocampal region CA3 via the mossy fibers, which form very large and anatomically distinct synapses called mossy fiber bouton–thorny excrescence synaptic complexes in the stratum lucidum (SL). CA3 pyramidal neurons then integrate information from their autoassociational connections in the stratum radiatum (SR) and stratum oriens (SO), with both direct entorhinal inputs in stratum lacunosum-moleculare (SLM) and the dentate gyrus inputs in the SL, and convey this information to hippocampal CA1 pyramidal neurons. Neurons in hippocampal CA1 then integrate this information in their basal and apical SR dendrites with direct entorhinal cortical inputs in their most distal, tufted dendrites in the SLM, and represent the first and largest extrahippocampal output from the hippocampus proper. Thus, the computations performed both within individual hippocampal subregions and between them as an interconnected neurocognitive system are complex, and involve a combination of intrinsic (i.e., membrane-bound ion channels that regulate membrane excitability) and synaptic (i.e., ligand-gated ion channels expressed at both excitatory and inhibitory synapses) influences. Additionally, age-related changes at any level of these complex circuits will have downstream consequences on the accuracy/reliability of the information being relayed to extrahippocampal regions via CA1 pyramidal neurons.Given the amount of evidence supporting a possible synaptic explanation for age-related learning and memory impairments in hippocampus-dependent forms of cognition (6–10), we combined patch-clamp physiology, serial section conventional and immunogold electron microscopy (EM), quantitative Western blot analyses, and behavioral characterization using two hippocampus-dependent forms of learning in YA (6- to 8-mo old) and aged rats (28- to 29-mo old) to examine two interconnected hippocampal regions implicated in cognitive aging: Regions CA1 and CA3. We focused on CA1 and CA3 because of their central location in the hippocampal circuit (48–50), their similar laminar dendritic structure (48–50), and their well-documented age-related changes in place field specificity and reliability (51–56).We find that the synaptic architecture and balance of synaptic weights in YA and aged, learning-unimpaired (AU) rats is remarkably similar, but that both are different in aged, learning-impaired (AI) rats. Moreover, this restructuring among “unsuccessful” cognitive agers has an intriguing specificity: It involves only AMPARs, only perforated axospinous synapses, and only hippocampal region CA3, which together shift the balance of synaptic weights that drive action potential output in CA3 pyramidal neurons maladaptively toward an overemphasis of the autoassociational synapses that interconnect CA3 pyramidal neurons.     

Circular swimming motility and disordered hyperuniform state in an algae system

Active matter comprises individually driven units that convert locally stored energy into mechanical motion. Interactions between driven units lead to a variety of nonequilibrium collective phenomena in active matter. One of such phenomena is anomalously large density fluctuations, which have been observed in both experiments and theories. Here we show that, on the contrary, density fluctuations in active matter can also be greatly suppressed. Our experiments are carried out with marine algae (Effrenium voratum), which swim in circles at the air–liquid interfaces with two different eukaryotic flagella. Cell swimming generates fluid flow that leads to effective repulsions between cells in the far field. The long-range nature of such repulsive interactions suppresses density fluctuations and generates disordered hyperuniform states under a wide range of density conditions. Emergence of hyperuniformity and associated scaling exponent are quantitatively reproduced in a numerical model whose main ingredients are effective hydrodynamic interactions and uncorrelated random cell motion. Our results demonstrate the existence of disordered hyperuniform states in active matter and suggest the possibility of using hydrodynamic flow for self-assembly in active matter.

Active matter exists over a wide range of spatial and temporal scales (1–6) from animal groups (7, 8) to robot swarms (9–11), to cell colonies and tissues (12–16), to cytoskeletal extracts (17–20), to man-made microswimmers (21–25). Constituent particles in active matter systems are driven out of thermal equilibrium at the individual level; they interact to develop a wealth of intriguing collective phenomena, including clustering (13, 22, 24), flocking (11, 26), swarming (12, 13), spontaneous flow (14, 20), and giant density fluctuations (10, 11). Many of these observed phenomena have been successfully described by particle-based or continuum models (1–6), which highlight the important roles of both individual motility and interparticle interactions in determining system dynamics.Current active matter research focuses primarily on linearly swimming particles which have a symmetric body and self-propel along one of the symmetry axes. However, a perfect alignment between the propulsion direction and body axis is rarely found in reality. Deviation from such a perfect alignment leads to a persistent curvature in the microswimmer trajectories; examples of such circle microswimmers include anisotropic artificial micromotors (27, 28), self-propelled nematic droplets (29, 30), magnetotactic bacteria and Janus particles in rotating external fields (31, 32), Janus particle in viscoelastic medium (33), and sperm and bacteria near interfaces (34, 35). Chiral motility of circle microswimmers, as predicted by theoretical and numerical investigations, can lead to a range of interesting collective phenomena in circular microswimmers, including vortex structures (36, 37), localization in traps (38), enhanced flocking (39), and hyperuniform states (40). However, experimental verifications of these predictions are limited (32, 35), a situation mainly due to the scarcity of suitable experimental systems.Here we address this challenge by investigating marine algae Effrenium voratum (41, 42). At air–liquid interfaces, E. voratum cells swim in circles via two eukaryotic flagella: a transverse flagellum encircling the cellular anteroposterior axis and a longitudinal one running posteriorly. Over a wide range of densities, circling E. voratum cells self-organize into disordered hyperuniform states with suppressed density fluctuations at large length scales. Hyperuniformity (43, 44) has been considered as a new form of material order which leads to novel functionalities (45–49); it has been observed in many systems, including avian photoreceptor patterns (50), amorphous ices (51), amorphous silica (52), ultracold atoms (53), soft matter systems (54–61), and stochastic models (62–64). Our work demonstrates the existence of hyperuniformity in active matter and shows that hydrodynamic interactions can be used to construct hyperuniform states.     

Modulation of immune cell reactivity with cis-binding Siglec agonists

Inflammatory pathologies caused by phagocytes lead to numerous debilitating conditions, including chronic pain and blindness due to age-related macular degeneration. Many members of the sialic acid-binding immunoglobulin-like lectin (Siglec) family are immunoinhibitory receptors whose agonism is an attractive approach for antiinflammatory therapy. Here, we show that synthetic lipid-conjugated glycopolypeptides can insert into cell membranes and engage Siglec receptors in cis, leading to inhibitory signaling. Specifically, we construct a cis-binding agonist of Siglec-9 and show that it modulates mitogen-activated protein kinase (MAPK) signaling in reporter cell lines, immortalized macrophage and microglial cell lines, and primary human macrophages. Thus, these cis-binding agonists of Siglecs present a method for therapeutic suppression of immune cell reactivity.

Sialic acid-binding immunoglobulin (IgG)-like lectins (Siglecs) are a family of immune checkpoint receptors that are on all classes of immune cells (1–5). Siglecs bind various sialoglycan ligands and deliver signals to the immune cells that report on whether the target is healthy or damaged, “self” or “nonself.” Of the 14 human Siglecs, 9 contain cytosolic inhibitory signaling domains. Accordingly, engagement of these inhibitory Siglecs by sialoglycans suppresses the activity of the immune cell, leading to an antiinflammatory effect. In this regard, inhibitory Siglecs have functional parallels with the T cell checkpoint receptors CTLA-4 and PD-1 (6–9). As with these clinically established targets for cancer immune therapy, there has been a recent surge of interest in antagonizing Siglecs to potentiate immune cell reactivity toward cancer (10). Conversely, engagement of Siglecs with agonist antibodies can suppress immune cell reactivity in the context of antiinflammatory therapy. This approach has been explored to achieve B cell suppression in lupus patients by agonism of CD22 (Siglec-2) (11, 12), and to deplete eosinophils for treatment of eosinophilic gastroenteritis by agonism of Siglec-8 (13). Similarly, a CD24 fusion protein has been investigated clinically as a Siglec-10 agonist for both graft-versus-host disease and viral infection (14, 15).Traditionally, Siglec ligands have been studied as functioning in trans, that is, on an adjacent cell (16–18), or as soluble clustering agents (9, 19). In contrast to these mechanisms of action, a growing body of work suggests that cis ligands for Siglecs (i.e., sialoglycans that reside on the same cell membrane) cluster these receptors and maintain a basal level of inhibitory signaling that increases the threshold for immune cell activation. Both Bassik and coworkers (20) and Wyss-Coray and coworkers (21) have linked the depletion of cis Siglec ligands with increased activity of macrophages and microglia, and other studies have shown that a metabolic blockade of sialic acid renders phagocytes more prone to activation (22).Synthetic ligands are a promising class of Siglec agonists (17, 23, 24). Many examples rely on clustering architectures (e.g., sialopolymers, nanoparticles, liposomes) to induce their effect (19, 23–26). Indeed, we have previously used glycopolymers to study the effects of Siglec engagement in trans on natural killer (NK) cell activity (16). We and other researchers have employed glycopolymers (16, 23), glycan-remodeling enzymes (27, 28), chemical inhibitors of glycan biosynthesis (22), and mucin overexpression constructs (29, 30) to modulate the cell-surface levels of Siglec ligands. However, current approaches lack specificity for a given Siglec.We hypothesized that Siglec-specific cis-binding sialoglycans displayed on immune cell surfaces could dampen immune cell activity with potential therapeutic applications. Here we test this notion with the synthesis of membrane-tethered cis-binding agonists of Siglec-9 (Fig. 1). Macrophages and microglia widely express Siglec-9 and are responsible for numerous pathologies including age-related inflammation (31), macular degeneration (32), neural inflammation (33), and chronic obstructive pulmonary disease (34). We designed and developed a lipid-linked glycopolypeptide scaffold bearing glycans that are selective Siglec-9 ligands (pS9L-lipid). We show that pS9L-lipid inserts into macrophage membranes, binds Siglec-9 specifically and in cis, and induces Siglec-9 signaling to suppress macrophage activity. By contrast, a lipid-free soluble analog (pS9L-sol) binds Siglec-9 but does not agonize Siglec-9 or modulate macrophage activity. Membrane-tethered glycopolypeptides are thus a potential therapeutic modality for inhibiting phagocyte activity.Open in a separate windowFig. 1.Lipid-tethered glycopolypeptides cluster and agonize Siglecs in cis on effector cells. (A) Immune cells express activating receptors that stimulate inflammatory signaling. (B) Clustering of Siglec-9 by cis-binding agonists stimulates inhibitory signaling that quenches activation.     

Plant species richness at archaeological sites suggests ecological legacy of Indigenous subsistence on the Colorado Plateau

Humans have both intentional and unintentional impacts on their environment, yet identifying the enduring ecological legacies of past small-scale societies remains difficult, and as such, evidence is sparse. The present study found evidence of an ecological legacy that persists today within an semiarid ecosystem of western North America. Specifically, the richness of ethnographically important plant species is strongly associated with archaeological complexity and ecological diversity at Puebloan sites in a region known as Bears Ears on the Colorado Plateau. A multivariate model including both environmental and archaeological predictors explains 88% of the variation in ethnographic species richness (ESR), with growing degree days and archaeological site complexity having the strongest effects. At least 31 plant species important to five tribal groups (Navajo, Hopi, Zuni, Ute Mountain Ute, and Apache), including the Four Corners potato (Solanum jamesii), goosefoot (Chenopodium sp.), wolfberry (Lycium pallidum), and sumac (Rhus trilobata), occurred at archaeological sites, despite being uncommon across the wider landscape. Our results reveal a clear ecological legacy of past human behavior: even when holding environmental variables constant, ESR increases significantly as a function of past investment in habitation and subsistence. Consequently, we suggest that propagules of some species were transported and cultivated, intentionally or not, establishing populations that persist to this day. Ensuring persistence will require tribal input for conserving and restoring archaeo-ecosystems containing “high-priority” plant species, especially those held sacred as lifeway medicines. This transdisciplinary approach has important implications for resource management planning, especially in areas such as Bears Ears that will experience greater visitation and associated impacts in the near future.

Local resource abundance is important for determining where in a given landscape humans decide to live. Nearby water, game, soil, and plants provide readily available wild resources for foraging and conditions that allow for cultivation (1–5). However, humans also modify their surrounding environments in order to increase the abundance and diversity of local plant (6–11) and animal (12–15) resources. Such “human niche construction” is a hallmark of ancient and modern societies (16, 17), having positive and negative impacts on global biodiversity while possibly creating enduring ecological legacies (18–21). This may be especially true for more sedentary and dense populations (22, 4) that are more likely to find investment worthwhile (23) and to produce unintentional impacts. Thus, variation in contemporary ecological diversity may in part reflect past land use dynamics and, therefore, be revealed through coupled archaeological and ecological research (24–33).Coupled ecological and archaeological research has led to the discovery of altered patterns of succession resulting from 1) forest clearing and changes in canopy light regime (34, 35), 2) alterations of soil especially linked to food refuse (36, 37), 3) changes in fire regimes (38, 39), and, more rarely, 4) the importation of plant propagules from distant sites of collection (40, 41). Identifying such long-lost dynamics between humans and landscapes can inform conservation aimed at restoring site-specific artifacts, features, and the associated resource base past and present, here termed “archaeo-ecosystems” (42, 43). This would greatly facilitate cross-cultural management of public lands (44) in ways that promote Indigenous health, cultural reclamation, and sovereignty (7, 45). The linkages, however, between ecological legacies, archaeo-ecosystem restoration and cross-cultural management have yet to be systematically tested or practically applied.Here, we offer a formal evaluation of this archaeo-ecosystem approach by using paired archaeological and ecological survey data focused on Puebloan occupation of a region known as Bears Ears on the Colorado Plateau in southeastern Utah (Fig. 1). Puebloan populations modified their environment by constructing terraces and check dams, developing blinds and wing traps, importing exogenous species, and setting fires (4, 22, 46), but investments were not uniform across the region. We test the hypothesis that locations with greater investment indicated by larger and more complex archaeological sites should today have higher richness of culturally significant plant species, here termed ethnographic species richness (ESR), as an enduring legacy of past investment. Our study expands previous work on ecological legacies by using field surveys to develop an explanatory model applied to 265 sites across one million acres of semiarid public lands. It documents the occurrence of uncommon and ethnographically significant plant species associated with those sites and infuses traditional ecological knowledge into proposed management actions for conserving these archaeo-ecosystems. Controlling for underlying environmental variation, our results indicate that past human habitation increases the diversity of plant species important for Indigenous subsistence.Open in a separate windowFig. 1.Location of Bears Ears National Monument in southeastern Utah. The predicted ESR at 265 known archaeological sites across the original and reduced monument boundaries and surrounding region are shown.     

ICAM-1 induced rearrangements of capsid and genome prime rhinovirus 14 for activation and uncoating

Most rhinoviruses, which are the leading cause of the common cold, utilize intercellular adhesion molecule-1 (ICAM-1) as a receptor to infect cells. To release their genomes, rhinoviruses convert to activated particles that contain pores in the capsid, lack minor capsid protein VP4, and have an altered genome organization. The binding of rhinoviruses to ICAM-1 promotes virus activation; however, the molecular details of the process remain unknown. Here, we present the structures of virion of rhinovirus 14 and its complex with ICAM-1 determined to resolutions of 2.6 and 2.4 Å, respectively. The cryo-electron microscopy reconstruction of rhinovirus 14 virions contains the resolved density of octanucleotide segments from the RNA genome that interact with VP2 subunits. We show that the binding of ICAM-1 to rhinovirus 14 is required to prime the virus for activation and genome release at acidic pH. Formation of the rhinovirus 14–ICAM-1 complex induces conformational changes to the rhinovirus 14 capsid, including translocation of the C termini of VP4 subunits, which become poised for release through pores that open in the capsids of activated particles. VP4 subunits with altered conformation block the RNA–VP2 interactions and expose patches of positively charged residues. The conformational changes to the capsid induce the redistribution of the virus genome by altering the capsid–RNA interactions. The restructuring of the rhinovirus 14 capsid and genome prepares the virions for conversion to activated particles. The high-resolution structure of rhinovirus 14 in complex with ICAM-1 explains how the binding of uncoating receptors enables enterovirus genome release.

Human rhinoviruses are the cause of more than half of common colds (1). Medical visits and missed days of school and work cost tens of billions of US dollars annually (2, 3). There is currently no cure for rhinovirus infections, and the available treatments are only symptomatic. Rhinoviruses belong to the family Picornaviridae, genus Enterovirus, and are classified into species A, B, and C (4). Rhinoviruses A and B can belong to either “major” or “minor” groups, based on their utilization of intercellular adhesion molecule-1 (ICAM-1) or low-density lipoprotein receptor for cell entry (5–7). Type C rhinoviruses use CDHR3 as a receptor (8). Rhinovirus 14 belongs to the species rhinovirus B and uses ICAM-1 as a receptor. Receptors recognized by rhinoviruses and other enteroviruses can be divided into two groups based on their function in the infection process (9). Attachment receptors such as DAF, PSGL1, KREMEN1, CDHR3, and sialic acid enable the binding and endocytosis of virus particles into cells (10–13). In contrast, uncoating receptors including ICAM-1, CD155, CAR, and SCARB2 enable virus cell entry but also promote genome release from virus particles (5, 14–16).Virions of rhinoviruses are nonenveloped and have icosahedral capsids (17). Genomes of rhinoviruses are 7,000 to 9,000 nucleotide-long single-stranded positive-sense RNA molecules (1, 17). The rhinovirus genome encodes a single polyprotein that is co- and posttranslationally cleaved into functional protein subunits. Capsid proteins VP1, VP3, and VP0, originating from one polyprotein, form a protomer, 60 of which assemble into a pseudo-T = 3 icosahedral capsid. To render the virions mature and infectious, VP0 subunits are cleaved into VP2 and VP4 (18, 19). VP1 subunits form pentamers around fivefold symmetry axes, whereas subunits VP2 and VP3 form heterohexamers centered on threefold symmetry axes. The major capsid proteins VP1 through 3 have a jelly roll β-sandwich fold formed by two β-sheets, each containing four antiparallel β-strands, which are conventionally named B to I (20–22). The two β-sheets contain the strands BIDG and CHEF, respectively. The C termini of the capsid proteins are located at the virion surface, whereas the N termini mediate interactions between the capsid proteins and the RNA genome on the inner surface of the capsid. VP4 subunits are attached to the inner face of the capsid formed by the major capsid proteins. The surfaces of rhinovirus virions are characterized by circular depressions called canyons, which are centered around fivefold symmetry axes of the capsids (21).The VP1 subunits of most rhinoviruses, but not those of rhinovirus 14, contain hydrophobic pockets, which are filled by molecules called pocket factors (17, 21, 23, 24). It has been speculated that pocket factors are fatty acids or lipids (25). The pockets are positioned immediately below the canyons. The exposure of rhinoviruses to acidic pH induces expulsion of the pocket factors, which leads to the formation of activated particles and genome release (17, 26–32). The activated particles are characterized by capsid expansion, a reduction in interpentamer contacts, the release of VP4 subunits, externalization of N termini of VP1 subunits, and changes in the distribution of RNA genomes (17, 26–29, 33, 34). Artificial hydrophobic compounds that bind to VP1 pockets with high affinity inhibit infection by rhinoviruses (35, 36).ICAM-1 is an endothelial- and leukocyte-associated protein that stabilizes cell–cell interactions and facilitates the movement of leukocytes through endothelia (37). ICAM-1 can be divided into an extracellular amino-terminal part composed of five immunoglobulin domains, a single transmembrane helix, and a 29-residue–long carboxyl-terminal cytoplasmic domain. The immunoglobulin domains are characterized by a specific fold that consists of seven to eight β-strands, which form two antiparallel β-sheets in a sandwich arrangement (38–40). The immunoglobulin domains of ICAM-1 are stabilized by disulfide bonds and glycosylation (38–41). The connections between the immunoglobulin domains are formed by flexible linkers that enable bending of the extracellular part of ICAM-1. For example, the angle between domains 1 and 2 differs by 8° between molecules in distinct crystal forms (38, 42). As a virus receptor, ICAM-1 enables the virus particles to sequester at the cell surface and mediates their endocytosis (5). The structures of complexes of rhinoviruses 3, 14, and 16, and CVA21 with ICAM-1 have been determined to resolutions of 9 to 28 Å (42–46). It was shown that ICAM-1 molecules bind into the canyons at the rhinovirus surface, approximately between fivefold and twofold symmetry axes (42–46). ICAM-1 interacts with residues from all three major capsid proteins. It has been speculated that the binding of ICAM-1 triggers the transition of virions of rhinovirus 14 to activated particles and initiates genome release (45, 47). However, the limited resolution of the previous studies prevented characterization of the corresponding molecular mechanism.Here, we present the cryo-electron microscopy (cryo-EM) reconstruction of the rhinovirus 14 virion, which contains resolved density of octanucleotide segments of the RNA genome that interact with VP2 subunits. Furthermore, we show that the binding of ICAM-1 to rhinovirus 14 induces changes in its capsid and genome, which are required for subsequent virus activation and genome release at acidic pH.