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.     

A rate threshold mechanism regulates MAPK stress signaling and survival

Cells are exposed to changes in extracellular stimulus concentration that vary as a function of rate. However, how cells integrate information conveyed from stimulation rate along with concentration remains poorly understood. Here, we examined how varying the rate of stress application alters budding yeast mitogen-activated protein kinase (MAPK) signaling and cell behavior at the single-cell level. We show that signaling depends on a rate threshold that operates in conjunction with stimulus concentration to determine the timing of MAPK signaling during rate-varying stimulus treatments. We also discovered that the stimulation rate threshold and stimulation rate-dependent cell survival are sensitive to changes in the expression levels of the Ptp2 phosphatase, but not of another phosphatase that similarly regulates osmostress signaling during switch-like treatments. Our results demonstrate that stimulation rate is a regulated determinant of cell behavior and provide a paradigm to guide the dissection of major stimulation rate dependent mechanisms in other systems.

All cells employ signal transduction pathways to respond to physiologically relevant changes in extracellular stressors, nutrient levels, hormones, morphogens, and other stimuli that vary as functions of both concentration and rate in healthy and diseased states (1–7). Switch-like “instantaneous” changes in the concentrations of stimuli in the extracellular environment have been widely used to show that the strength of signaling and overall cellular response are dependent on the stimulus concentration, which in many cases needs to exceed a certain threshold (8, 9). Previous studies have shown that the rate of stimulation can also influence signaling output in a variety of pathways (10–17) and that stimulation profiles of varying rates can be used to probe underlying signaling pathway circuitry (4, 18, 19). However, it is still not clear how cells integrate information conveyed by changes in both the stimulation rate and concentration in determining signaling output. It is also not clear if cells require stimulation gradients to exceed a certain rate in order to commence signaling.Recent investigations have demonstrated that stimulation rate can be a determining factor in signal transduction. In contrast to switch-like perturbations, which trigger a broad set of stress-response pathways, slow stimulation rates activate a specific response to the stress applied in Bacillus subtilis cells (10). Meanwhile, shallow morphogen gradient stimulation fails to activate developmental pathways in mouse myoblast cells in culture, even when concentrations sufficient for activation during pulsed treatment are delivered (12). These observations raise the possibility that stimulation profiles must exceed a set minimum rate or rate threshold to achieve signaling activation. Although such rate thresholds would help cells decide if and how to respond to dynamic changes in stimulus concentration, the possibility of signaling regulation by a rate threshold has never been directly investigated in any system. Further, no study has experimentally examined how stimulation rate requirements impact cell phenotype or how cells molecularly regulate the stimulation rate required for signaling activation. As such, the biological significance of any existing rate threshold regulation of signaling remains unknown.The budding yeast Saccharomyces cerevisiae high osmolarity glycerol (HOG) pathway provides an ideal model system for addressing these issues (Fig. 1A). The evolutionarily conserved mitogen-activated protein kinase (MAPK) Hog1 serves as the central signaling mediator of this pathway (20–22). It is well established that instantaneous increases in osmotic stress concentration induce Hog1 phosphorylation, activation, and translocation to the nucleus (18, 21, 23–30). Activated Hog1 governs the majority of the cellular osmoadaptation response that enables cells to survive (23, 31, 32). Multiple apparently redundant MAPK phosphatases dephosphorylate and inactivate Hog1, which, along with the termination of upstream signaling after adaptation, results in its return to the cytosol (Fig. 1A) (23, 25, 26, 33–39). Because of this behavior, time-lapse analysis of Hog1 nuclear enrichment in single cells has proven an excellent and sensitive way to monitor signaling responses to dynamic stimulation patterns in real time (18, 27–30, 40, 41). Further, such assays have been readily combined with traditional growth and molecular genetic approaches to link observed signaling responses with cell behavior and signaling pathway architecture (27–29).Open in a separate windowFig. 1.Hog1 signaling and cell survival are sensitive to the rate of preconditioning osmotic stress application. (A) Schematic of the budding yeast HOG response. (B) Preconditioning protection assay workflow indicating the first stress treatments to a final concentration of 0.4 M NaCl (Left), high-stress exposure (Middle), and colony formation readout (Right). (C) High-stress survival as a function of each first treatment relative to the untreated first stress condition. Bars and errors are means and SD from three biological replicates. *Statistically significant by Kolmogorov–Smirnov test (P < 0.05). NS = not significant. (D) Treatment concentration over time. (E) Treatment rate over time for quadratic and pulse treatment. The rate for the pulse is briefly infinite (blue vertical line) before it drops to 0. (F) Hog1 nuclear localization during the treatments depicted in D and E. (Inset) Localization pattern in the quadratic-treated sample. Lines represent means and shaded error represents the SD from three to four biological replicates.Here, we use systematically designed osmotic stress treatments imposed at varying rates of increase to show that a rate threshold condition regulates yeast high-stress survival and Hog1 MAPK signaling. We demonstrate that only stimulus profiles that satisfy both this rate threshold condition and a concentration threshold condition result in robust signaling. We go on to show that the protein tyrosine phosphatase Ptp2, but not the related Ptp3 phosphatase, serves as a major rate threshold regulator. By expressing PTP2 under the control of a series of different enhancer–promoter DNA constructs, we demonstrate that changes in the level of Ptp2 expression can alter the stimulation rate required for signaling induction and survival. These findings establish rate thresholds as a critical and regulated component of signaling biology akin to concentration thresholds.     

From the Cover: Datura quids at Pinwheel Cave,California, provide unambiguous confirmation of the ingestion of hallucinogens at a rock art site

While debates have raged over the relationship between trance and rock art, unambiguous evidence of the consumption of hallucinogens has not been reported from any rock art site in the world. A painting possibly representing the flowers of Datura on the ceiling of a Californian rock art site called Pinwheel Cave was discovered alongside fibrous quids in the same ceiling. Even though Native Californians are historically documented to have used Datura to enter trance states, little evidence exists to associate it with rock art. A multianalytical approach to the rock art, the quids, and the archaeological context of this site was undertaken. Liquid chromatography−mass spectrometry (LC-MS) results found hallucinogenic alkaloids scopolamine and atropine in the quids, while scanning electron microscope analysis confirms most to be Datura wrightii. Three-dimensional (3D) analyses of the quids indicate the quids were likely masticated and thus consumed in the cave under the paintings. Archaeological evidence and chronological dating shows the site was well utilized as a temporary residence for a range of activities from Late Prehistory through Colonial Periods. This indicates that Datura was ingested in the cave and that the rock painting represents the plant itself, serving to codify communal rituals involving this powerful entheogen. These results confirm the use of hallucinogens at a rock art site while calling into question previous assumptions concerning trance and rock art imagery.

Since the late 1980s, the role that altered states of consciousness (or ASC) played in the making of rock art has been one of the most contentious questions confronted by rock art researchers across the globe (1–7). The ASC model purports that humans universally experience three distinct visual phases during trance, which are replicated in rock art imagery (1). The ASC model can be induced in a number of ways including the use of hallucinogenic substances (1). However, there remains no clear evidence for the preparation and consumption of hallucinogenic substances directly associated with any rock art site in the world. Indeed, fierce debate has occurred over the last 30 y, with many researchers questioning the validity of the ACS model and the idea of shamanism as a viable explanation for the creation of rock art (2–6). California has been central within this debate (8, 9). Whitley (9) has argued that the many south-central Californian rock paintings were shamanic self-portraits depicting a shaman’s experience during ASC while rock art sites were owned by individual shamans, and avoided by the local populace. In this view, trance, shamanism, and rock art are inextricably linked in their separation from normal activity of the wider populace. However, evidence from systematic archaeological work in south-central California has clearly shown that the majority of rock art sites were integrated into habitation sites, and are not separated from public view (10, 11). Recent analyses also suggest that the pictographs were probably not self-depictions of shamans in trance but, instead, stock iconographic images drawing upon mythology and the personifying of insects, animals, plant, and astronomical elements such as the sun (12, 13).Even so, ethnographic documentation details how hallucinogens played a pivotal role in Native California, especially Datura wrightii (14). A member of the Solanacae family, Datura is distinctive by large white “trumpet” flowers that uncoil in a five-pointed pinwheeling fashion. Datura as a genus can be found across multiple continents, including the Americas, Asia, Europe, and South Africa (15). Its wide availability and hallucinogenic properties, due to the presence of the tropane alkaloids atropine and scopolamine, are behind its use across different cultures (16, 17). The most noted usage of Datura in Native California is in youth initiations where the root was processed into a drink or “tea” known historically as toloache (14, 18–21). Initiates would often be instructed in cultural rules of entering adulthood and how to interpret the visions themselves (14, 19, 21). For some, these ceremonies where highly codified, such as the Chinigchinich religion of Southern California, which ended with the making of a sand painting in which boys learned religious principles (19, 21). The sand paintings did not depict the visions induced by Datura, but, instead, were cosmological maps detailing the ontological principles of the Southern Californian Native societies making them (22). After undergoing the puberty ceremony, Datura could be taken throughout one’s lifetime for a variety of reasons, including to gain supernatural power for doctoring, to counteract negative supernatural events, to ward off ghosts, and to see the future or find lost objects, but, most especially, as a mendicant for a variety of ailments (14, 18, 23). Datura consumption could occur prior to hunting to increase stamina and power (24). Importantly, Datura could be consumed in a variety of ways, including drinking toloache, but also by roasting the roots, eating the flowers or seeds, applying poultices on wounds, or often simply chewing the roots or other parts of the plant (14). For the Tübatulabal, Datura originated as a man named Mo mo ht who subsequently turned into the plant in its present form (25), while, in Chumash mythology, the plant was a prominent supernatural grandmother called Momoy (26). Since Datura, with its psychoactive substances, was used within spiritual, ritual, and mythological contexts, it should be considered as an entheogen (27).Worldwide, different hallucinogens have been suggested as inspiring rock art making, such as mushrooms, Peyote, Datura, San Pedro cactus, Brunsvigia, and others (2, 28–35). Datura is of particular focus in the North American West. Malotki (33) argues that Datura influenced archaic Basketmaker rock art in Arizona, while Boyd’s (32) extensive analyses have shown that Lower Pecos rock art iconography likely relates to mythological narratives concerning Peyote and Datura. Images include hornworms, the larval stages of hawkmoths, Datura’s primary pollinator. Mimbres pottery and Kiva murals also include representations of Datura, plus anthropomorphized versions of the hawkmoth (36, 37). Evidence of Datura alkaloids have been found in ceramics (38), while ancient peyote buttons and Datura seeds have been found in Lower Pecos archaeological deposits, but none have been reported specifically at rock art sites (39).Datura has also been suggested as relating to the making of Chumash rock art in California (40). Found in the Chumash borderlands of interior south-central California, we detail here our investigations at a rock art site called Pinwheel Cave (CA-KER-5836) (Fig. 1) and its associated food processing bedrock mortar (BRM) complex (CA-KER-5837). The name originates from a large red pinwheel motif, which we hypothesized may represent the opening Datura flower (Fig. 1, Bottom Left and Bottom Right). Dozens of fibrous clumps known as quids are also located within crevices in the cave ceiling (41). Quids, usually found in archaeological deposits, are typically made of yucca, agave, tule, or tobacco and are thought to have been chewed to extract nutrients or stimulants (41–44). With the painting likely representing the opening of the Datura flower, and with the very unusual insertion of quids in the ceiling, we investigated the possibility that the quids could contain Datura. We present the results of a multianalytical investigation of the quids and the archaeological context to investigate the potential use of Datura in association with rock art iconography.Open in a separate windowFig. 1.Pinwheel Cave, California. (Top) Interior of cave during laser scanning. (Bottom Left) Pinwheel painting within cave. Image credit: Rick Bury (photographer). (Bottom Right) Unfurling flower of D. wrightii from plant near cave site. Image credit: Melissa Dabulamanzi (photographer).     

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.     

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.     

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.     

Heading perception depends on time-varying evolution of optic flow

There is considerable support for the hypothesis that perception of heading in the presence of rotation is mediated by instantaneous optic flow. This hypothesis, however, has never been tested. We introduce a method, termed “nonvarying phase motion,” for generating a stimulus that conveys a single instantaneous optic flow field, even though the stimulus is presented for an extended period of time. In this experiment, observers viewed stimulus videos and performed a forced-choice heading discrimination task. For nonvarying phase motion, observers made large errors in heading judgments. This suggests that instantaneous optic flow is insufficient for heading perception in the presence of rotation. These errors were mostly eliminated when the velocity of phase motion was varied over time to convey the evolving sequence of optic flow fields corresponding to a particular heading. This demonstrates that heading perception in the presence of rotation relies on the time-varying evolution of optic flow. We hypothesize that the visual system accurately computes heading, despite rotation, based on optic acceleration, the temporal derivative of optic flow.

James Gibson first remarked that the instantaneous motion of points on the retina (Fig. 1A) can be formally described as a two-dimensional (2D) field of velocity vectors called the “optic flow field” (or “optic flow”) (1). Such optic flow, caused by an observer’s movement relative to the environment, conveys information about self-motion and the structure of the visual scene (1–15). When an observer translates in a given direction along a straight path, the optic flow field radiates from a point in the image with zero velocity, or singularity, called the focus of expansion (Fig. 1B). It is well known that under such conditions, one can accurately estimate one’s “heading” (i.e., instantaneous direction of translation in retinocentric coordinates) by simply locating the focus of expansion (SI Appendix). However, if there is angular rotation in addition to translation (by moving along a curved path or by a head or eye movement), the singularity in the optic flow field will be displaced such that it no longer corresponds to the true heading (Fig. 1 C and D). In this case, if one estimates heading by locating the singularity, the estimate will be biased away from the true heading. This is known as the rotation problem (14).Open in a separate windowFig. 1.Projective geometry, the rotation problem, time-varying optic flow, and the optic acceleration hypothesis. (A) Viewer-centered coordinate frame and perspective projection. Because of motion between the viewpoint and the scene, a 3D surface point traverses a path in 3D space. Under perspective projection, the 3D path of this point projects onto a 2D path in the image plane (retina), the temporal derivative of which is called image velocity. The 2D velocities associated with all visible points define a dense 2D vector field called the optic flow field. (B–D) Illustration of the rotation problem. (B) Optic flow for pure translation (1.5-m/s translation speed, 0° heading, i.e., heading in the direction of gaze). Optic flow singularity (red circle) corresponds to heading (purple circle). (C) Pure rotation, for illustrative purposes only and not corresponding to any experimental condition (2°/s rightward rotation). (D) Translation + rotation (1.5 m/s translation speed, 0° heading, 2°/s rightward rotation). Optic flow singularity (red circle) is displaced away from heading (purple circle). (E) Three frames from a video depicting movement along a circular path with the line-of-sight initially perpendicular to a single fronto-parallel plane composed of black dots. (F) Time-varying evolution of optic flow. The first optic flow field reflects image motion between the first and second frames of the video. The second optic flow field reflects image motion between the second and third frames of the video. For this special case (circular path), the optic flow field evolves (and the optic flow singularity drifts) only due to the changing depth of the environment relative to the viewpoint. (G) Illustration of the optic acceleration hypothesis. Optic acceleration is the derivative of optic flow over time (here, approximated as the difference between the second and first optic flow fields). The singularity of the optic acceleration field corresponds to the heading direction. Acceleration vectors autoscaled for visibility.Computer vision researchers and vision scientists have developed a variety of algorithms that accurately and precisely extract observer translation and rotation from optic flow, thereby solving the rotation problem. Nearly all of these rely on instantaneous optic flow (i.e., a single optic flow field) (4, 9, 16–25) with few exceptions (26–29). However, it is unknown whether these algorithms are commensurate with the neural computations underlying heading perception.The consensus of opinion in the experimental literature is that human observers can estimate heading (30, 31) from instantaneous optic flow, in the absence of additional information (5, 10, 15, 32–34). Even so, there are reports of systematic biases in heading perception (11); the visual consequences of rotation (eye, head, and body) can bias heading judgments (10, 15, 35–37), with the amount of bias typically proportional to the magnitude of rotation. Other visual factors, such as stereo cues (38, 39), depth structure (8, 10, 40–43), and field of view (FOV) (33, 42–44) can modulate the strength of these biases. Errors in heading judgments have been reported to be greater when eye (35–37, 45, 46) or head movements (37) are simulated versus when they are real, which has been taken to mean that observers require extraretinal information, although there is also evidence to the contrary (10, 15, 33, 40, 41, 44, 47–50). Regardless, to date no one has tested whether heading perception (even with these biases) is based on instantaneous optic flow or on the information available in how the optic flow field evolves over time. Some have suggested that heading estimates rely on information accumulated over time (32, 44, 51), but no one has investigated the role of time-varying optic flow without confounding it with stimulus duration (i.e., the duration of evidence accumulation).In this study, we employed an application of an image processing technique that ensured that only a single optic flow field was available to observers, even though the stimulus was presented for an extended period of time. We called this condition “nonvarying phase motion” or “nonvarying”: The phases of two component gratings comprising each stationary stimulus patch shifted over time at a constant rate, causing a percept of motion in the absence of veridical movement (52). Phase motion also eliminated other cues that may otherwise have been used for heading judgments, including image point trajectories (15, 32) and their spatial compositions (i.e., looming) (53, 54). For nonvarying phase motion, observers exhibited large biases in heading judgments in the presence of rotation. A second condition, “time-varying phase motion,” or “time-varying,” included acceleration by varying the velocity of phase motion over time to match the evolution of a sequence of optic flow fields. Doing so allowed observers to compensate for the confounding effect of rotation on optic flow, making heading perception nearly veridical. This demonstrates that heading perception in the presence of rotation relies on the time-varying evolution of optic flow.     

Structure,self-assembly,and properties of a truncated reflectin variant

Naturally occurring and recombinant protein-based materials are frequently employed for the study of fundamental biological processes and are often leveraged for applications in areas as diverse as electronics, optics, bioengineering, medicine, and even fashion. Within this context, unique structural proteins known as reflectins have recently attracted substantial attention due to their key roles in the fascinating color-changing capabilities of cephalopods and their technological potential as biophotonic and bioelectronic materials. However, progress toward understanding reflectins has been hindered by their atypical aromatic and charged residue-enriched sequences, extreme sensitivities to subtle changes in environmental conditions, and well-known propensities for aggregation. Herein, we elucidate the structure of a reflectin variant at the molecular level, demonstrate a straightforward mechanical agitation-based methodology for controlling this variant’s hierarchical assembly, and establish a direct correlation between the protein’s structural characteristics and intrinsic optical properties. Altogether, our findings address multiple challenges associated with the development of reflectins as materials, furnish molecular-level insight into the mechanistic underpinnings of cephalopod skin cells’ color-changing functionalities, and may inform new research directions across biochemistry, cellular biology, bioengineering, and optics.

Materials from naturally occurring and recombinant proteins are frequently employed for the study of fundamental biological processes and leveraged for applications in fields as diverse as electronics, optics, bioengineering, medicine, and fashion (1–13). Such broad utility is enabled by the numerous advantageous characteristics of protein-based materials, which include sequence modularity, controllable self-assembly, stimuli-responsiveness, straightforward processability, inherent biological compatibility, and customizable functionality (1–13). Within this context, unique structural proteins known as reflectins have recently attracted substantial attention because of their key roles in the fascinating color-changing capabilities of cephalopods, such as the squid shown in Fig. 1A, and have furthermore demonstrated their utility for unconventional biophotonic and bioelectronic technologies (11–40). For example, in vivo, Bragg stack-like ultrastructures from reflectin-based high refractive index lamellae (membrane-enclosed platelets) are responsible for the angle-dependent narrowband reflectance (iridescence) of squid iridophores, as shown in Fig. 1B (15–20). Analogously, folded membranes containing distributed reflectin-based particle arrangements within sheath cells lead to the mechanically actuated iridescence of squid chromatophore organs, as shown in Fig. 1C (15, 16, 21, 22). Moreover, in vitro, films processed from squid reflectins not only exhibit proton conductivities on par with some state-of-the-art artificial materials (23–27) but also support the growth of murine and human neural stem cells (28, 29). Additionally, morphologically variable coatings assembled from different reflectin isoforms can enable the functionality of chemically and electrically actuated color-changing devices, dynamic near-infrared camouflage platforms, and stimuli-responsive photonic architectures (27, 30–34). When considered together, these discoveries and demonstrations constitute compelling motivation for the continued exploration of reflectins as model biomaterials.Open in a separate windowFig. 1.(A) A camera image of a D. pealeii squid for which the skin contains light-reflecting cells called iridophores (bright spots) and pigmented organs called chromatophores (colored spots). Image credit: Roger T. Hanlon (photographer). (B) An illustration of an iridophore (Left), which shows internal Bragg stack-like ultrastructures from reflectin-based lamellae (i.e., membrane-enclosed platelets) (Inset). (C) An illustration of a chromatophore organ (Left), which shows arrangements of reflectin-based particles within the sheath cells (Inset). (D) The logo of the 28-residue-long N-terminal motif (RMN), which depicts the constituent amino acids (Upper) and their predicted secondary structures (Lower). (E) The logo of the 28-residue-long internal motif (RMI), which depicts the constituent amino acids (Upper) and their predicted secondary structures (Lower). (F) The logo of the 21-residue-long C-terminal motif (RMC), which depicts the constituent amino acids (Upper) and their predicted secondary structures (Lower). (G) The amino acid sequence of full-length D. pealeii reflectin A1, which contains a single RMN motif (gray oval) and five RMI motifs (orange ovals). (H) An illustration of the selection of the prototypical truncated reflectin variant (denoted as RfA1TV) from full-length D. pealeii reflectin A1.Given reflectins’ demonstrated significance from both fundamental biology and applications perspectives, some research effort has been devoted to resolving their three-dimensional (3D) structures (30, 31, 35–39). For example, fibers drawn from full-length Euprymna scolopes reflectin 1a and films processed from truncated E. scolopes reflectin 1a were shown to possess secondary structural elements (i.e., α-helices or β-sheets) (30, 31). In addition, precipitated nanoparticles and drop-cast films from full-length Doryteuthis pealeii reflectin A1 have exhibited β-character, which was seemingly associated with their conserved motifs (35, 36). Moreover, nanoparticles assembled from both full-length and truncated Sepia officinalis reflectin 2 variants have demonstrated signatures consistent with β-sheet or α-helical secondary structure, albeit in the presence of surfactants (38). However, such studies were made exceedingly challenging by reflectins’ atypical primary sequences enriched in aromatic and charged residues, documented extreme sensitivities to subtle changes in environmental conditions, and well-known propensities for poorly controlled aggregation (12, 14, 15, 30–32, 34–39). Consequently, the reported efforts have all suffered from multiple drawbacks, including the need for organic solvents or denaturants, the evaluation of only polydisperse or aggregated (rather than monomeric) proteins, a lack of consensus among different experimental techniques, inadequate resolution that precluded molecular-level insight, imperfect agreement between computational predictions and experimental observations, and/or the absence of conclusive correlations between structure and optical functionality. As such, there has emerged an exciting opportunity for investigating reflectins’ molecular structures, which remain poorly understood and the subject of some debate.Herein, we elucidate the structure of a reflectin variant at the molecular level, demonstrate a robust methodology for controlling this variant’s hierarchical assembly, and establish a direct correlation between its structural characteristics and optical properties. We first rationally select a prototypical reflectin variant expected to recapitulate the behavior of its parent protein by using a bioinformatics-guided approach. We next map the conformational and energetic landscape accessible to our selected protein by means of all-atom molecular dynamics (MD) simulations. We in turn produce our truncated reflectin variant with and without isotopic labeling, develop solution conditions that maintain the protein in a monomeric state, and characterize the variant’s size and shape with small-angle X-ray scattering (SAXS). We subsequently resolve our protein’s dynamic secondary and tertiary structures and evaluate its backbone conformational fluctuations with NMR spectroscopy. Finally, we demonstrate a straightforward mechanical agitation-based approach to controlling our truncated reflectin variant’s secondary structure, hierarchical self-assembly, and bulk refractive index distribution. Overall, our findings address multiple challenges associated with the development of reflectins as materials, furnish molecular-level insight into the mechanistic underpinnings of cephalopod skin cells’ color-changing functionalities, and appear poised to inform new directions across biochemistry, cellular biology, bioengineering, and optics.     

From the Cover: Neurophysiological coordination of duet singing

Coordination of behavior for cooperative performances often relies on linkages mediated by sensory cues exchanged between participants. How neurophysiological responses to sensory information affect motor programs to coordinate behavior between individuals is not known. We investigated how plain-tailed wrens (Pheugopedius euophrys) use acoustic feedback to coordinate extraordinary duet performances in which females and males rapidly take turns singing. We made simultaneous neurophysiological recordings in a song control area “HVC” in pairs of singing wrens at a field site in Ecuador. HVC is a premotor area that integrates auditory feedback and is necessary for song production. We found that spiking activity of HVC neurons in each sex increased for production of its own syllables. In contrast, hearing sensory feedback produced by the bird’s partner decreased HVC activity during duet singing, potentially coordinating HVC premotor activity in each bird through inhibition. When birds sang alone, HVC neurons in females but not males were inhibited by hearing the partner bird. When birds were anesthetized with urethane, which antagonizes GABAergic (γ-aminobutyric acid) transmission, HVC neurons were excited rather than inhibited, suggesting a role for GABA in the coordination of duet singing. These data suggest that HVC integrates information across partners during duets and that rapid turn taking may be mediated, in part, by inhibition.

Animals routinely rely on sensory feedback for the control of their own behavior. In cooperative performances, such sensory feedback can include cues produced by other participants (1–8). For example, in interactive vocal communication, including human speech, individuals take turns vocalizing. This “turn taking” is a consequence of each participant responding to auditory cues from a partner (4–6, 9, 10). The role of such “heterogenous” (other-generated) feedback in the control of vocal turn taking and other cooperative performances is largely unknown.Plain-tailed wrens (Pheugopedius euophrys) are neotropical songbirds that cooperate to produce extraordinary duet performances but also sing by themselves (Fig. 1A) (4, 10, 11). Singing in plain-tailed wrens is performed by both females and males and used for territorial defense and other functions, including mate guarding and attraction (1, 11–16). During duets, female and male plain-tailed wrens take turns, alternating syllables at a rate of between 2 and 5 Hz (Fig. 1A) (4, 11).Open in a separate windowFig. 1.Neural control of solo and duet singing in plain-tailed wrens. (A) Spectrogram of a singing bout that included male solo syllables (blue line, top) followed by a duet. Solo syllables for both sexes (only male solo syllables are shown here) are sung at lower amplitudes than syllables produced in duets. Note that the smeared appearance of wren syllables in spectrograms reflects the acoustic structure of plain-tailed wren singing. (B and C) Each bird has a motor system that is used to produce song and sensory systems that mediate feedback. (B) During solo singing, the bird hears its own song, which is known as autogenous feedback (orange). (C) During duet singing, each bird hears both its own singing and the singing of its partner, known as heterogenous feedback (green). The key difference between solo and duet singing is heterogenous feedback that couples the neural systems of the two birds. This coupling results in changes in syllable amplitude and timing in both birds.There is a categorical difference between solo and duet singing. In solo singing, the singing bird receives only autogenous (hearing its own vocalization) feedback (Fig. 1B). The partner may hear the solo song if it is nearby, a heterogenous (other-generated) cue. In duet singing, birds receive both heterogenous and autogenous feedback as they alternate syllable production (Fig. 1C). Participants use heterogenous feedback during duet singing for precise timing of syllable production (4, 11). For example, when a male temporarily stops participating in a duet, the duration of intersyllable intervals between female syllables increases (4), showing an effect of heterogenous feedback on the timing of syllable production.How does the brain of each wren integrate heterogenous acoustic cues to coordinate the precise timing of syllable production between individuals during duet performances? To address this question, we examined neurophysiological activity in HVC, a nucleus in the nidopallium [an analogue of mammalian cortex (17, 18)]. HVC is necessary for song learning, production, and timing in species of songbirds that do not perform duets (19–24). Neurons in HVC are active during singing and respond to playback of the bird’s own learned song (25–27). In addition, recent work has shown that HVC is also involved in vocal turn taking (19).To examine the role of heterogenous feedback in the control of duet performances, we compared neurophysiological activity in HVC when female or male wrens sang solo syllables with syllables sung during duets. Neurophysiological recordings were made in awake and anesthetized pairs of wrens at the Yanayacu Biological Station and Center for Creative Studies on the slopes of the Antisana volcano in Ecuador. We found that heterogenous cues inhibited HVC activity during duet performances in both females and males, but inhibition was only observed in females during solo singing.     

Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic

The mode and extent of rapid evolution and genomic change in response to human harvesting are key conservation issues. Although experiments and models have shown a high potential for both genetic and phenotypic change in response to fishing, empirical examples of genetic responses in wild populations are rare. Here, we compare whole-genome sequence data of Atlantic cod (Gadus morhua) that were collected before (early 20th century) and after (early 21st century) periods of intensive exploitation and rapid decline in the age of maturation from two geographically distinct populations in Newfoundland, Canada, and the northeast Arctic, Norway. Our temporal, genome-wide analyses of 346,290 loci show no substantial loss of genetic diversity and high effective population sizes. Moreover, we do not find distinct signals of strong selective sweeps anywhere in the genome, although we cannot rule out the possibility of highly polygenic evolution. Our observations suggest that phenotypic change in these populations is not constrained by irreversible loss of genomic variation and thus imply that former traits could be reestablished with demographic recovery.

As anthropogenic activities rapidly transform the environment, a fundamental question is whether wild populations have the capacity to adapt and evolve fast enough in response (1–3). Phenotypic change can result from phenotypic plasticity, but emerging examples of genomic change over only a few generations have made clear that rapid evolution is also possible (4–6). In the literature, one of the most dramatic and widely cited cases involves the declining age and size at maturation of Atlantic cod (Gadus morhua) following several generations of high fishing pressure (3, 7–10). Fisheries produce some of the fastest rates of phenotypic change ever observed in wild populations (2, 11), but the extent to which fisheries-induced evolution has occurred in the wild and the degree to which it is reversible remain strongly debated (12).The hypothesis that evolution underlies these phenotypic changes is supported by a range of observations. For example, theory on the selective nature of many fisheries reveals that higher rates of harvesting will—with only a few exceptions—favor earlier sexual maturation, greater investment in reproduction, and slower growth (13). In addition, experiments in the laboratory that selectively remove large or small individuals from a population reveal rapid evolution of body size and maturation time in only a few generations, as well as substantial impacts on fishery yields (14–16). Fisheries-induced evolution experiments in the laboratory also reveal selective sweeps through dramatic shifts in allele frequencies, loss of genetic diversity, and increases in linkage disequilibrium at specific locations in the genome (15, 17, 18).However, translating these findings to wild populations has been substantially more difficult. One concern is that phenotypic plasticity, gene flow, or spatial shifts in populations can also explain the substantial phenotypic and limited genotypic changes reported from the wild to date (10, 13, 19–23). The magnitude and rate of fisheries-induced evolution may also be quite small in the wild (19). While theory provides strong evidence that fishing can be a potent driver of evolutionary changes, a clear empirical demonstration of fisheries-induced evolution would require evidence that the observed change is genetic (13). Whether and to what extent the widespread genomic reorganization observed in experiments also occurs in wild-harvested populations therefore remain unknown.Genomic analyses of temporal samples before and after selective events have provided key opportunities to test for rapid adaptive evolution from standing genetic variation in wild populations by identifying unusually strong shifts in allele frequencies over time (4, 5). In addition, the history of genomic research with Atlantic cod (24, 25) provides a unique opportunity to test for genomic signatures of fisheries-induced evolution in particular. Archival samples collected by fisheries scientists decades or even centuries ago represent a valuable source of historical genomic material that can provide rare insight into the genetic patterns of the past (26). Here, we obtained whole-genome sequence data from well-preserved archives of Atlantic cod scales and otoliths (ear bones) that were originally collected from two populations on either side of the Atlantic Ocean: the northeast Arctic population sampled near Lofoten, Norway in 1907 and the Canadian northern cod population sampled near Twillingate, Newfoundland in 1940 (Fig. 1A and SI Appendix, Table S1). The Canadian northern population collapsed from overfishing in the early 1990s, while the northeast Arctic population experienced high fishing rates but smaller declines in biomass (10, 27, 28). Both populations have shown marked reductions in age at maturation, though with slight increases in maturation age in northeast Arctic cod after 2005 (Fig. 1B). We compared these historical genomes with modern data from the same locations (Fig. 1A and SI Appendix, Table S2). In total, we analyzed 113 individual genomes (Methods) from these two unique populations that had independently experienced intensive fishing during the last century (7, 10). We found a marked lack of large genomic changes or selective sweeps through time, suggesting instead that phenotypic plasticity or, potentially, highly polygenic evolution can explain the observed changes in phenotype.Open in a separate windowFig. 1.Spatiotemporal population structure based on genome-wide data in Atlantic cod. (A) In total, 113 modern and historical specimens were analyzed from northern cod collected in Newfoundland, Canada (1940, yellow; 2013, dark yellow) and from northeast Arctic cod collected in the Lofoten archipelago, Norway (1907, orange; modern: 2011, red; 2014, dark red). (B) Age at 50% maturity over time in each population. (C) PCA as implemented in PCAngsd. Velicier’s minimum average partial (MAP) test identified a single significant PC and only one PC is shown. Individuals are colored according to A. (D) Model-based ADMIXTURE ancestry components for historical (1907, 1940) and modern (2013, 2014) populations (k = 2; NGSadmix). Each individual is represented by a column colored to show the proportion of each ancestry component for Canada (dark yellow) and Norway (orange). Population differentiation based on pairwise weighted F_ST is also shown. (E) The correlation between the allele frequencies in historical and modern populations. Colors reflect the relative density of points, from darker (more density) to lighter (less density). R², coefficient of correlation.     

Cingulo-opercular control network and disused motor circuits joined in standby mode

Whole-brain resting-state functional MRI (rs-fMRI) during 2 wk of upper-limb casting revealed that disused motor regions became more strongly connected to the cingulo-opercular network (CON), an executive control network that includes regions of the dorsal anterior cingulate cortex (dACC) and insula. Disuse-driven increases in functional connectivity (FC) were specific to the CON and somatomotor networks and did not involve any other networks, such as the salience, frontoparietal, or default mode networks. Censoring and modeling analyses showed that FC increases during casting were mediated by large, spontaneous activity pulses that appeared in the disused motor regions and CON control regions. During limb constraint, disused motor circuits appear to enter a standby mode characterized by spontaneous activity pulses and strengthened connectivity to CON executive control regions.

Disuse is a powerful paradigm for inducing plasticity that has uncovered key organizing principles of the human brain (1–4). Monocular deprivation—prolonged covering of one eye—revealed that multiple afferent inputs can compete for representational territory in the primary visual cortex (1). Similar competition between afferents also shapes the somatomotor system. Manipulations such as peripheral nerve deafferentation, whisker trimming, and limb constraint all drive plasticity in the primary somatosensory and motor cortex (2–4). Most plasticity studies to date have used focal techniques, such as microelectrode recordings, to study local changes in brain function. As a result, little is known about how behavior and experience shape the brain-wide functional networks that support complex cognitive operations (5).The brain is composed of networks of regions that cooperate to perform specific cognitive functions (5–8). These functional networks show synchronized spontaneous activity while the brain is at rest, a phenomenon known as resting-state functional connectivity (FC) (9–11). FC can be measured noninvasively in humans using resting-state functional MRI (rs-fMRI) and has been used to parse the brain into canonical functional networks (12, 13), including visual, auditory, and somatomotor networks (14, 15); ventral and dorsal attention networks (8, 16); a default mode network with roles in internally directed cognition and episodic memory (7, 11); a salience network thought to assess the homeostatic relevance of external stimuli (17); a frontoparietal control network supporting error processing and moment-to-moment adjustments in behavior (18–20); and a cingulo-opercular control network (CON), which maintains executive control during goal-directed behavior (18, 19, 21). Each functional network likely carries out a variety of additional functions.A more recent advance in human neuroscience has been the recognition of individual variability in network organization (22–25). Most early rs-fMRI studies examined central tendencies in network organization using group-averaged FC measurements (10, 12, 13). Recent work has demonstrated that functional networks can be identified in an individual-specific manner if sufficient rs-fMRI data are acquired, an approach termed precision functional mapping (PFM) (22, 23, 26–30). PFM respects the unique functional anatomy of each person and avoids averaging together functionally distinct brain regions across individuals.We recently demonstrated that PFM can be used to follow the time course of disuse-driven plasticity in the human brain (31). Three adult participants (Nico, Ashley, and Omar) were scanned at the same time of day for 42 to 64 consecutive days (30 min of rs-fMRI per day) before, during, and after 2 wk of dominant upper-extremity casting (Fig. 1 A and B). Casting caused persistent disuse of the dominant upper extremity during daily behaviors and led to a marked loss of strength and fine motor skill in all participants. During casting, the upper-extremity regions of the left primary somatomotor cortex (L-SM1_ue) and right cerebellum (R-Cblm_ue) functionally disconnected from the remainder of the somatomotor network. Disused motor circuits also exhibited large, spontaneous pulses of activity (Fig. 1C). Disuse pulses did not occur prior to casting, started to occur frequently within 1 to 2 d of casting, and quickly waned after cast removal.Open in a separate windowFig. 1.Experimental design and spontaneous activity pulses. (A) Three participants (Nico, Ashley, and Omar) wore casts covering the entire dominant upper extremity for 2 wk. (B) Participants were scanned every day for 42 to 64 consecutive days before, during, and after casting. All scans included 30 min of resting-state functional MRI. (C) During the Cast period, disused somatomotor circuits exhibited large pulses of spontaneous activity. (C, Left) Whole-brain ANOVA showing which brain regions contained disuse-driven pulses. (C, Right) Time courses of all pulses recorded from the disused primary somatomotor cortex.Somatomotor circuits do not function in isolation. Action selection and motor control are thought to be governed by complex interactions between the somatomotor network and control networks, including the CON (18). Prior studies of disuse-driven plasticity, including our own, have focused solely on somatomotor circuits. Here, we leveraged the whole-brain coverage of rs-fMRI and the statistical power of PFM to examine disuse-driven plasticity throughout the human brain.