Climate warming changes synchrony of plants and pollinators Climate warming changes the phenology of many species. When interacting organisms respond differently, climate change may disrupt their interactions and affect the stability of ecosystems. Here, we used GBIF occurrence records to examine phenology trends in plants and their associated insect pollinators in Germany since the 1980s. We found strong phenological advances in plants, but differences in the extent of shifts among pollinator groups. The temporal trends in plant and insect phenologies were generally associated with interannual temperature variation, and thus likely driven by climate change. When examining the synchrony of species-level plant-pollinator interactions, their temporal trends differed among pollinator groups. Overall, plant-pollinator interactions become more synchronized, mainly since the phenology of plants responded more strongly to climate change than that of the pollinators. However, if the observed trends continue, many interactions may become more asynchronous again in the future. Our study suggests that climate change affects the phenologies of both plants and insects, and that it also influences the synchrony of plant-pollinator interactions.
Rapid global phaseout of animal agriculture has the potential to stabilize greenhouse gas levels for 30 years and offset 68 percent of CO2 emissions Animal agriculture contributes significantly to global warming through ongoing emissions of the potent greenhouse gases methane and nitrous oxide, and displacement of biomass carbon on the land used to support livestock. However, because estimates of the magnitude of the effect of ending animal agriculture often focus on only one factor, the full potential benefit of a more radical change remains underappreciated. Here we quantify the full "climate opportunity cost" of current global livestock production, by modeling the combined, long-term effects of emission reductions and biomass recovery that would be unlocked by a phaseout of animal agriculture. We show that, even in the absence of any other emission reductions, persistent drops in atmospheric methane and nitrous oxide levels, and slower carbon dioxide accumulation, following a phaseout of livestock production would, through the end of the century, have the same cumulative effect on the warming potential of the atmosphere as a 25 gigaton per year reduction in anthropogenic CO2 emissions, providing half of the net emission reductions necessary to limit warming to 2{degrees}C. The magnitude and rapidity of these potential effects should place the reduction or elimination of animal agriculture at the forefront of strategies for averting disastrous climate change.

Significance StatementThe use of animals to produce food has a negative impact on the climate, but the benefits of a global switch to a plant based diet are underappreciated. We show that the global warming impact, through the rest of this century, of eliminating greenhouse gas emissions from livestock and allowing native ecosystems to regrow on the land currently used to house and feed livestock, would be equivalent to a 68% reduction in carbon dioxide emissions. We hope putting clearer numbers on the "climate opportunity cost" of our continued use of animals as food technology will help policymakers and the public properly prioritize dietary change as a climate defense strategy.

Declaration of Conflict of InterestPatrick Brown is the founder and CEO of Impossible Foods, a company developing alternatives to animals in food-production. Michael Eisen is an advisor to Impossible Foods. Both are shareholders in the company and thus stand to benefit financially from reduction of animal agriculture.
Sexual selection for males with beneficial mutations Sexual selection is the process by which traits providing a mating advantage are favoured. Theoretical treatments of the evolution of sex by sexual selection propose that it operates by reducing the load of deleterious mutations. Here, we postulate instead that sexual selection primarily acts through females preferentially mating with males carrying beneficial mutations. We used simulation and analytical modelling to investigate the evolutionary dynamics of beneficial mutations in the presence of sexual selection. We found that female choice for males with beneficial mutations had a much greater impact on genetic quality than choice for males with low mutational load. We also relaxed the typical assumption of a fixed mutation rate. For deleterious mutations, mutation rate should always be minimized, but when rare beneficial mutations can occur, female choice for males with those rare beneficial mutations could overcome a decline in average fitness and allow an increase in mutation rate. We propose that sexual selection for beneficial mutations could overcome the  two-fold cost of sex much more readily than choice for males with low mutational load and may therefore be a more powerful explanation for the prevalence of sexual reproduction than the existing theory. If sexual selection results in higher fitness at higher mutation rates and if the variability produced by mutation itself promotes sexual selection, then a feedback loop between these two factors could have had a decisive role in driving adaptation. Stimulus manipulations permit activation of fiber subpopulations in the mouse and rat vagus Afferent and efferent vagal fibers mediate bidirectional communication between the brain and visceral organs. Small, unmyelinated C-afferents constitute the majority of vagal fibers, play critical roles in numerous interoceptive circuits and autonomic reflexes in health and disease and may contribute to the efficacy and safety of vagus nerve stimulation (VNS). Selective engagement of C-afferents with electrical stimuli has not been feasible, due to the default fiber recruitment order: larger fibers first, smaller fibers last. Here, we determine and optimize an electrical stimulus that selectively engages vagal C-afferents. Intermittent KHz-frequency electrical stimulation (KES) activates motor and, preferentially, sensory vagal neurons in the brainstem. During KES, asynchronous activity of C-afferents increases, while that of larger fibers remains largely unchanged. In parallel, KES effectively blocks excitability of larger fibers while moderately suppressing excitability of C-afferents. By compiling selectivity indices in individual animals, we find that optimal KES parameters for C-afferents are >5KHz frequency and 7-10 times engagement threshold (xT) intensity in rats, 15-25xT in mice. These effects can be explained in computational models by how sodium channel responses to KES are shaped by axonal size and myelin. Our results indicate that selective engagement of vagal C-afferents is attainable by intermittent KES. Epigenetic scores for the circulating proteome as tools for disease prediction Protein biomarkers have been identified across many age-related morbidities. However, characterising epigenetic influences could further inform disease predictions. Here, we leverage epigenome-wide data to study links between the DNAm signatures of the circulating proteome and incident diseases. Using data from four cohorts, we trained and tested epigenetic scores (EpiScores) for 953 plasma proteins, identifying 109 scores that explained between 1% and 58% of the variance in protein levels after adjusting for known protein quantitative trait loci (pQTL) genetic effects. By projecting these EpiScores into an independent sample, (Generation Scotland; n=9,537) and relating them to incident morbidities over a follow-up of 14 years, we uncovered 137 EpiScore - disease associations. These associations were largely independent of immune cell proportions, common lifestyle and health factors and biological aging. Notably, we found that our diabetes-associated EpiScores highlighted previous top biomarker associations from proteome-wide assessments of diabetes. These EpiScores for protein levels can therefore be a valuable resource for disease prediction and risk stratification. The Hippocampus uses cryptographic principles rather than plasticity to enable animals to secretly cache and retrieve their food For animals, the ability to hide and retrieve valuable information, such as the location of food, can mean the difference between life and death. Here, we propose that to achieve this, their brain uses spatial cells similarly to how we utilize encryption for data security. Some animals are able to cache hundreds of thousands of food items annually by each individual and later retrieve most of what they themselves stashed. Rather than memorizing their cache locations as previously suggested, we propose that they use a single cryptographic-like mechanism during both caching and retrieval. The model we developed is based on hippocampal spatial cells, which respond to an animals positional attention, such as when the animal enters a specific region (place-cells) or gazes at a particular location (spatial-view-cells). We know that the region that activates each spatial cell remains consistent across subsequent visits to the same area but not between areas. This remapping, combined with the uniqueness of cognitive maps, produces a persistent crypto-hash function for both food caching and retrieval. We also show that the model stores temporal information that helps animals in food caching order preference as previously observed. This behavior, which we refer to as crypto-taxis, might also explain consistent differences in decision-making when animals are faced with a large number of alternatives such as in foraging. Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA As the focus for CRISPR edited plants moves from proof-of-concept to real world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one gRNA per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in Populus tremula x P. alba INRA 717-1B4. Unexpected duplications of MYB186 and MYB138 resulted in a total of eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and PCR analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the non-glandular trichomes of poplar.

ONE SENTENCE SUMMARYTargeting conserved sequences with a single gRNA allowed efficient mutagenesis of a multigene family and the recovery of trichomeless and triterpene-free poplar mutants. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. Interpreting ciliopathy-associated missense variants of uncertain significance (VUS) in an animal model PurposeBetter methods are required to interpret the pathogenicity of disease-associated variants of uncertain significance (VUS), which cannot be actioned clinically. In this study, we explore the use of a tractable animal model (Caenorhabditis elegans) for in vivo interpretation of missense VUS alleles of TMEM67, a cilia gene associated with ciliopathies.

MethodsCRISPR/Cas9 gene editing was used to generate homozygous knock-in C. elegans worm strains carrying TMEM67 patient variants engineered into the orthologous gene (mks-3). Quantitative phenotypic assays of sensory cilia structure and function measured if the variants affect mks-3 gene function. Results from worms were validated by a genetic complementation assay in a human TMEM67 knock-out hTERT-RPE1 cell line that tests a TMEM67 signaling function.

ResultsAssays in C. elegans accurately distinguished between known benign (Asp359Glu, Thr360Ala) and known pathogenic (Glu361Ter, Gln376Pro) variants. Analysis of eight missense VUS generated evidence that three are benign (Cys173Arg, Thr176Ile, Gly979Arg) and five are pathogenic (Cys170Tyr, His782Arg, Gly786Glu, His790Arg, Ser961Tyr).

ConclusionEfficient genome editing and quantitative functional assays in C. elegans make it a tractable in vivo animal model that allows rapid, cost-effective interpretation of ciliopathy-associated missense VUS alleles. Inhibitors of ERp44, PDIA1, and AGR2 induce disulfide-mediated oligomerization of Death Receptors 4 and 5 and cancer cell death Breast cancer mortality remains unacceptably high, indicating a need for safer and more effective therapeutic agents. Disulfide bond Disrupting Agents (DDAs) were previously identified as a novel class of anticancer compounds that selectively kill cancers that overexpress the Epidermal Growth Factor Receptor (EGFR) or its family member HER2. DDAs kill EGFR+ and HER2+ cancer cells via the parallel downregulation of EGFR, HER2, and HER3 and activation/oligomerization of Death Receptors 4 and 5 (DR4/5). However, the mechanisms by which DDAs mediate these effects are unknown. Affinity purification analyses employing biotinylated-DDAs reveal that the Protein Disulfide Isomerase (PDI) family members AGR2, PDIA1, and ERp44 are DDA target proteins. Further analyses demonstrate that shRNA-mediated knockdown of AGR2 and ERp44, or expression of ERp44 mutants, enhance basal and DDA-induced DR5 oligomerization. DDA treatment of breast cancer cells disrupts PDIA1 and ERp44 mixed disulfide bonds with their client proteins. Together, the results herein reveal DDAs as the first small molecule, active site inhibitors of AGR2 and ERp44, and demonstrate roles for AGR2 and ERp44 in regulating the activity, stability, and localization of DR4 and DR5, and activation of Caspase 8. A novel Boolean network inference strategy to model early hematopoiesis aging Hematopoietic stem cell (HSC) aging is a multifactorial event that leads to changes in HSC properties and function. These changes are intrinsically coordinated and affect the early hematopoiesis, involving hematopoietic stem and progenitor cells (HSPCs). The objective of this work is to better understand the mechanisms and factors controlling these changes. We have therefore developed an original strategy to construct a Boolean network of genes explaining the priming and homeostasis of HSCs (graphical abstract). Based on our previous scRNA-seq data, we performed an exhaustive analysis of the transcriptional network and identified active transcription modules or regulons along the differentiation trajectory of selected HSPC states. This global view of transcriptional regulation led us to focus on 15 components, 13 selected TFs (Tal1, Fli1, Gata2, Gata1, Zfpm1, Egr1, Junb, Ikzf1, Myc, Cebpa, Bclaf1, Klf1, Spi1) and 2 complexes regulating the ability of HSC to cycle (CDK4/6 - Cyclin D and CIP/KIP). We then defined the connections controlling the differentiation dynamics of HSC states and constructed an influence graph between the TFs involved in the dynamics by mixing observations from our scRNA-seq data and knowledge from the literature. Then, using answer set programming (ASP) and in silico perturbation analysis, we obtained a Boolean model which is the solution of a Boolean satisfiability problem. Finally, perturbation of the model based on age-related changes revealed important regulations, such as the overactivation of Egr1 and Junb or the loss of Cebpa activation by Gata2, which were found to be relevant for the myeloid bias of aged HSC. Our work shows the efficiency of the combination of manual and systematic methods to elaborate a Boolean model. The developed strategy led to the proposal of new regulatory mechanisms underlying the differentiation bias of aged HSCs, explaining the decreased transcriptional priming of HSCs to all mature cell types except megakaryocytes. A tool for assessing alignment of biomedical data repositories with open, FAIR, citation and trustworthy principles Increasing attention is being paid to the operation of biomedical data repositories in light of efforts to improve how scientific data is handled and made available for the long term. Multiple groups have produced recommendations for functions that biomedical repositories should support, with many using requirements of the FAIR data principles as guidelines. However, FAIR is but one set of principles that has arisen out of the open science community. They are joined by principles governing open science, data citation and trustworthiness, all of which are important aspects for biomedical data repositories to support. Together, these define a framework for data repositories that we call OFCT: Open, FAIR, Citable and Trustworthy. Here we developed an instrument using the open source PolicyModels toolkit that attempts to operationalize key aspects of OFCT principles and piloted the instrument by evaluating eight biomedical community repositories listed by the NIDDK Information Network (dkNET.org). Repositories included both specialist repositories that focused on a particular data type or domain, in this case diabetes and metabolomics, and generalist repositories that accept all data types and domains. The goal of this work was both to obtain a sense of how much the design of current biomedical data repositories align with these principles and to augment the dkNET listing with additional information that may be important to investigators trying to choose a repository, e.g., does the repository fully support data citation? The evaluation was performed from March to November 2020 through inspection of documentation and interaction with the sites by the authors. Overall, although there was little explicit acknowledgement of any of the OFCT principles in our sample, the majority of repositories provided at least some support for their tenets. Agrobacterium mediated transformation and deciphering SNPs in TaLr67 gene homeologs for gene editing in wheat. Diseases adversely affect grain yield of crop plants. Leaf rust is a major disease of wheat. As race-specific resistance breaks down, introduction of newer sources of resistance often from older accessions is necessary. Linkage drag from the donor accessions adversely impacts grain yield. As a result, CIMMYT breeding effort has shifted to using durable resistance, also referred to as race non-specific host resistance, which allows slow rusting but maintains grain yield. One of the key genes for durable resistance is Lr67 and the resistant form of this gene is absent in CIMMYT elite lines. Hence, we have initiated efforts to convert the susceptible copy of Lr67 into its resistant form directly in these elite lines using gene editing. This would eliminate backcrossing and thus save time as well as eliminate linkage drag that would accompany the resistant copy of the gene if it were to be introgressed from an older accession. As first steps, we have isolated and sequenced the genomic copies from each of the A, B, and D genome of the Lr67 gene from three elite lines and an experimental line. Identification of more than 50 single nucleotide polymorphisms (SNPs) in the open reading frames (ORF) among these lines would be useful in designing the guide RNA molecules with precision. Further, we have streamlined genetic transformation of these elite lines, a prerequisite step for gene editing. How the replication and transcription complex functions in jumping transcription of SARS-CoV-2 BackgroundCoronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although unprecedented efforts are underway to develop therapeutic strategies against this disease, scientists have acquired only a little knowledge regarding the structures and functions of the CoV replication and transcription complex (RTC) and 16 non-structural proteins, named NSP1-16.

ResultsIn the present study, we proposed a two-route model to answer how the RTC functions in the jumping transcription of CoVs. The key step leading to this model was that the motif AAACH for METTL3 recognition flanking the transcription regulatory sequence (TRS) motif was discovered to determine the m6A methylation of SARS-CoV-2 RNAs, by reanalyzing public Nanopore RNA-seq data. As the most important finding, TRS hairpins were reported for the first time to interpret NSP15 cleavage, RNA methylation of CoVs and their association at the molecular level. In addition, we reported canonical TRS motifs of all CoVs to prove the importance of our findings.

ConclusionsThe main conclusions are: (1) TRS hairpins can be used to identify recombination regions in CoV genomes; (2) RNA methylation of CoVs participates in the determination of the RNA secondary structures by affecting the formation of base pairing; and (3) The eventual determination of the CoV RTC global structure needs to consider METTL3 in the experimental design. Our findings enrich fundamental knowledge in the field of gene expression and its regulation, providing a crucial basis for future studies. Classification of Pacemaker Dynamics in the Mouse Intestine by Field Potential Microimaging The flexible and sophisticated movement of the gastrointestinal (GI) tract implies the involvement of mechanisms other than enteric neural circuits, to coordinate excitation in microregions. We thus performed microimaging of pacemaker dynamics in the small intestine of mice since it contains typical network-forming pacemaker cells. A dialysis membrane-reinforced low-impedance microelectrode array (MEA) enabled field potentials over a wide frequency range to be stably measured in microregions. The pacemaker dynamics were classified into basic patterns despite large variations. In the developmental process, pacemaker activity was categorized as either an  expanding or a  migrating pattern that was initiated in or propagated to the MEA sensing area, respectively. The intercellular current of the volume conductor complicated the waveform of both activities. The existence of  expanding and  migrating patterns was attributable to duplicated pacemaker systems such as intracellular Ca2+ oscillation-activated and voltage-gated mechanisms. Additionally, from the spatio-temporal feature during the period of pacemaker events, the  bumpy/aberrant pattern was defined by aberrant, incoherent propagation, and associated with local impairment of excitability, while the  colliding/converging pattern involved the interaction of multiple activities in the MEA area. Interconversion between the four micro-coordination patterns occurred in the same microregion. 5-Hydroxytryptamine (5-HT) promoted  migrating activity, implying an improvement or restoration of spatial conductivity. These results agree well with the action of 5-HT to change GI movement toward propulsion. In conclusion, our MEA method of microimaging classification enables the quantitative assessment of spatio-temporal electric coordination underlying GI motility, suggesting its application to small model animals. Age differences in the functional architecture of the human brain The intrinsic functional organization of the brain changes into older adulthood. Age differences are observed at multiple spatial scales, from global reductions in modularity and segregation of distributed brain systems, to network-specific patterns of dedifferentiation. Whether dedifferentiation reflects an inevitable, global shift in brain function with age, circumscribed, experience dependent changes, or both, is uncertain. We employed a multi-method strategy to interrogate dedifferentiation at multiple spatial scales. Multi-echo (ME) resting-state fMRI was collected in younger (n=181) and older (n=120) healthy adults. Cortical parcellation sensitive to individual variation was implemented for precision functional mapping of each participant, while preserving group-level parcel and network labels. ME-fMRI processing and gradient mapping identified global and macroscale network differences. Multivariate functional connectivity methods tested for microscale, edge-level differences. Older adults had lower BOLD signal dimensionality, consistent with global network dedifferentiation. Gradients were largely age-invariant. Edge-level analyses revealed discrete, network-specific dedifferentiation patterns in older adults. Visual and somatosensory regions were more integrated within the functional connectome; default and frontoparietal control network regions showed greater connectivity; and the dorsal attention network was more integrated with heteromodal regions. These findings highlight the importance of multi-scale, multi-method approaches to characterize the architecture of functional brain aging. Rescue of stalled clathrin-mediated endocytosis by asymmetric Arp2/3-mediated actin assembly Actin assembly facilitates vesicle formation in several trafficking pathways. Clathrin-mediated endocytosis (CME) shows elevated actin assembly dependence under high membrane tension. Why actin assembly at CME sites occurs heterogeneously even within the same cell, and how assembly forces are harnessed, are not fully understood. Here, endocytic dynamics, actin presence, and geometry of CME proteins from three different functional modules, were analyzed using three-dimensional (3D) super-resolution microscopy, live-cell imaging, and machine-learning-based computation. When hundreds of CME events were compared, sites with actin assembly showed a distinct signature, a delay between completion of coat expansion and vesicle scission, indicating that actin assembly occurs preferentially at stalled CME sites. N-WASP is recruited to one side of CME sites where it is positioned to stimulate asymmetric actin assembly. We propose that asymmetric actin assembly rescues stalled CME sites by pulling vesicles into the cell much like a bottle opener pulls off a bottle cap. Cytoplasmic Colocalization of Granulins and TDP-43 Prion-like Domain Involves Electrostatically Driven Coacervation Tuned by the Redox State of Cysteines. Cytoplasmic inclusions containing aberrant proteolytic fragments of TDP-43 are associated with frontotemporal lobar degeneration (FTLD) and other related pathologies. In FTLD, TDP-43 is translocated into the cytoplasm and proteolytically cleaved to generate a prion-like domain (PrLD) containing C-terminal fragments (C25 and C35) that form toxic inclusions. Under stress, TDP-43 partitions into membraneless organelles called stress granules (SGs) by coacervating with RNA and other proteins. We were interested in understanding if and how cysteine-rich granulins (GRNs 1-7), which are the proteolytic products of a genetic risk factor in FTLD called progranulin, interact with TDP-43. We show that extracellular GRNs internalize and colocalize with PrLD as puncta in the cytoplasm of neuroblastoma cells but show no presence in SGs. In addition, we show GRNs and PrLD coacervate to undergo liquid-liquid phase separation (LLPS) or form gel- or solid-like aggregates. Identification of the sequence determinants within GRNs for the observed phase transitions reveal the negative charges to be the drivers of LLPS modulated by the positive charges and the redox state of cysteines. Furthermore, RNA and GRNs compete and expunge one another from PrLD condensates, providing a basis for GRNs absence in SGs. Together, the results illustrate the potential mechanisms by which extracellular GRNs, formed during chronic inflammatory conditions, could internalize, and modulate cytoplasmic TDP-43 inclusions in proteinopathies. Tobramycin-induced secretion of P. aeruginosa 5' tRNA-fMet halves suppresses lung inflammation via AGO2 gene silencing Although inhaled tobramycin increases lung function in people with cystic fibrosis (pwCF), the density of P. aeruginosa in the lungs is only modestly reduced by tobramycin; hence, the mechanism whereby tobramycin improves lung function is unclear. Here, we demonstrate that tobramycin increases the abundance of two 5' tRNA-fMet halves in outer membrane vesicles (OMVs) secreted by P. aeruginosa and that the 5' tRNA-fMet halves reduce IL-8 secretion by CF bronchial epithelial cells (CF-HBECs). In mouse lung, the 5' tRNA-fMet halves attenuate KC secretion and neutrophil recruitment. We also report that the 5' tRNA-fMet halves suppress pro-inflammatory network gene expression by an Argonaut 2 (AGO2)-mediated gene silencing mechanism, thereby reducing IL-8 secretion in CF-HBECs. Moreover, tobramycin reduces the IL-8 concentration and neutrophil content in bronchoalveolar lavage fluid of pwCF. Thus, we conclude that tobramycin improves lung function in part by reducing chronic inflammation and neutrophil-mediated lung damage in pwCF. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. SARS-CoV-2 Can Infect Human Embryos The spread of SARS-CoV-2 has led to a devastating pandemic, with infections resulting in a range of symptoms collectively known as COVID-19. The full repertoire of human tissues and organs susceptible to infection is an area of active investigation, and some studies have implicated the reproductive system. The effects of COVID-19 on human reproduction remain poorly understood, and particularly the impact on early embryogenesis and establishment of a pregnancy are not known. In this work, we explore the susceptibility of early human embryos to SARS-CoV-2 infection. We note that ACE2 and TMPRSS2, two canonical cell entry factors for SARS-CoV-2, are co-expressed in cells of the trophectoderm in blastocyst-stage preimplantation embryos. Using fluorescent reporter virions pseudotyped with Spike (S) glycoprotein from SARS-CoV-2, we observe robust infection of trophectoderm cells, and this permissiveness could be attenuated with blocking antibodies targeting S or ACE2. When exposing human blastocysts to the live, fully infectious SARS-CoV-2, we detected cases of infection that compromised embryo health. Therefore, we identify a new human target tissue for SARS-CoV-2 with potential medical implications for reproductive health during the COVID-19 pandemic and its aftermath. Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer The Epithelial-to-Mesenchymal Transition (EMT) is a developmental cellular program frequently coopted by cancer cells1 and is a key contributor to both heterogeneity in solid tumors2-4 and later stage chemo-resistance and metastasis5,6. Rather than being a switch from an epithelial to a mesenchymal state, increasing evidence points to the existence of intermediate EMT states, wherein cells co-express both epithelial and mesenchymal traits7-13. Multiple stable intermediate EMT states possessing unique characteristics exist across the EMT spectrum7,8,14,15, contributing to the complex heterogeneity of tumors to promote metastasis16. While much work has been carried out identifying and characterizing EMT-inducing transcription factors17-20, the transcriptional and epigenetic networks responsible for the stability and maintenance of the midpoints along the EMT spectrum are poorly defined. In addition, there are currently no approaches to identifying and quantifying intermediate EMT subpopulations within patient tumors to evaluate their prognostic significance. Using clonally isolated derivatives of the SUM149PT breast cancer cell line, we systematically interrogate how each EMT state independently contributes to heterogeneity and influences metastatic progression, uncovering the role of RUNX2 in stabilizing certain intermediate states. Using SUM149PT-derived tumors as a training set, we develop an entropy-based model to quantify phenotypic heterogeneity and EMT status. Remarkably, the majority of cell states captured in the SUM149PT model are represented in a selection of patient tumors, laying the foundation for quantification of epithelial-mesenchymal heterogeneity and understanding the role of the intermediate EMT state in tumor progression. A multiscale lipid and cellular atlas of the human kidney The human kidney is composed of over 26 cell types that actively coordinate with each other to form higher-order structures, such as the nephron. It is not yet understood how these structures vary throughout a single organ or amongst the same organs within the human population. We have developed an extensive lipid and cellular atlas of the human kidney consisting of over 3 million cells comprising 75,000 functional tissue units (i.e., glomeruli, proximal tubules, distal tubules, and collecting ducts) from 13 human subjects. This atlas was developed using spatially registered and integrated technologies consisting of imaging mass spectrometry, multiplexed immunofluorescence, stained microscopy, and autofluorescence microscopy, to comprehensively probe large (i.e., centimeter-sized) areas of tissue. The cellular organization and lipid profiles of glomeruli, proximal tubules, distal tubules, and collecting ducts were discovered through these multimodal imaging data as well as their intra- and inter-subject variance. Relating the lipid profiles obtained from imaging mass spectrometry to distinct cell types obtained from immunofluorescence allowed us to hypothesize the functional role of specific phospholipids that have not previously been described. These hypotheses include subject characteristics, such as BMI and sex. The integrated data from the aforementioned datasets provide a valuable reference for kidney researchers, are publicly available through the NIH Human Biomolecular Atlas Program (https://portal.hubmapconsortium.org/), and discussed below. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

MethodsAscending aortic pathologies were examined in patients with sporadic thoracic aortopathies and angiotensin II (AngII)-infused mice. Ascending aortas without overt pathology from AngII-infused mice were subjected to mass spectrometry assisted proteomics, and molecular features of SHF-derived cells were determined by single cell transcriptomic analyses. Genetic deletion of either low-density lipoprotein receptor-related protein 1 (Lrp1) or transforming growth factor-{beta} receptor 2 (Tgfbr2) in SHF- derived cells was conducted to examine the impact of SHF-derived cells on vascular integrity.

ResultsPathologies in human ascending aortic aneurysmal tissues were predominant in outer medial layers and adventitia. This gradient was mimicked in mouse aortas following AngII infusion that was coincident with the distribution of SHF-derived cells. Proteomics indicated that brief AngII infusion, prior to overt pathology, evoked downregulation of SMC proteins and differential expression of extracellular matrix proteins, including several LRP1 ligands. LRP1 deletion in SHF-derived cells augmented AngII-induced ascending aortic aneurysm and rupture. Single cell transcriptomic analysis revealed that brief AngII infusion decreased Lrp1 and Tgfbr2 mRNA abundance in SHF-derived cells and induced a unique fibroblast population with low abundance of Tgfbr2 mRNA. SHF-specific Tgfbr2 deletion led to embryonic lethality at E12.5 with dilatation of the outflow tract and retroperitoneal hemorrhage. Integration of proteomic and single cell transcriptomics results identified plasminogen activator inhibitor 1 (PAI1) as the most increased protein in SHF-derived SMCs and fibroblasts during AngII infusion. Immunostaining revealed a transmural gradient of PAI1 in both ascending aortas of AngII-infused mice and human ascending aneurysmal aortas that mimicked the gradient of medial and adventitial pathologies.

ConclusionSHF-derived cells exert a critical role in maintaining vascular integrity through LRP1 and TGF-{beta} signaling associated with increases of aortic PAI1.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=89 SRC="FIGDIR/small/930917v6_ufig1.gif" ALT="Figure 1">
View larger version (34K):
org.highwire.dtl.DTLVardef@185e3d9org.highwire.dtl.DTLVardef@131d702org.highwire.dtl.DTLVardef@11afa42org.highwire.dtl.DTLVardef@242786_HPS_FORMAT_FIGEXP  M_FIG C_FIG Clinical PerspectiveO_ST_ABSWhat is new?C_ST_ABS- SHF-derived SMCs and fibroblasts associate with AngII-induced aortic pathologies.
- AngII induces a distinct fibroblast sub-cluster that is less abundant for mRNAs related to major extracellular components and TGF{beta} ligands and receptors, but more abundant for proliferative genes.
- TGFBR2 deletion in SHF-derived cells are embryonic lethal with significant dilatation of the outflow tract in mice.
- SHF-specific deletion of LRP1 leads to aortic pathologies in mice, supporting the importance of SHF-derived cells in maintaining ascending aortic wall integrity.


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Recovery of high-quality assembled genomes via metagenome binning guided with single-cell amplified genomes BackgroundObtaining high-quality (HQ) reference genomes from microbial communities is crucial for understanding the phylogeny and function of uncultured microbes in complex microbial ecosystems. Despite the improved bioinformatic approaches to generate curated metagenome-assembled genomes (MAGs), existing metagenomic binners often fail to obtain reliable MAGs, and thus, they are nowhere comparable to genomes sequenced from isolates in terms of strain level resolution. Here, we present a single-cell genome-guided metagenome binning (MetaSAG) to reconstruct the strain-resolved genomes from microbial communities at once.

ResultsMetaSAG employs single-cell amplified genomes (SAGs) generated with microfluidic technology as binning guides to recover improved draft genomes with the metagenomic data. To assess the performance of reconstructing genomes from various microbial communities, we compared MetaSAG with four conventional metagenomic binners using a cell mock community, human gut microbiota, and skin microbiota samples. MetaSAG showed precise contig binning and higher recovery rates (>97%) of rRNA and plasmids compared to conventional binners in genome reconstruction from the cell mock community. In human microbiota samples, MetaSAG recovered the largest number of genomes with a total of 103 gut microbial genomes (21 HQ and 65 showed >90% completeness) and 45 skin microbial genomes (10 HQ and 40 showed >90% completeness), respectively. Conventional binners recovered one Staphylococcus hominis genome, whereas MetaSAG recovered two S. hominis genomes from the identical skin microbiota sample. Single-cell sequencing indicated that these S. hominis genomes clearly derived from two distinct strains harboring specifically different plasmids. We found that all conventional S. hominis MAGs had substantial lack or excess of the genome sequences and contamination of other Staphylococcus bacteria (S. epidermidis).

ConclusionsMetaSAG enabled us to obtain the strain-resolved genomes in the mock community and human microbiota samples by assigning metagenomic sequences correctly and covering both highly conserved genes such as rRNA genes and unique extrachromosomal elements, including plasmids. MetaSAG will provide HQ genomes that are difficult to obtain with metagenomic analyses alone and will facilitate the understanding of microbial ecosystems by elucidating detailed metabolic pathways and horizontal gene transfer networks through microbial genomes. MetaSAG is available at https://github.com/kojiari/metasag.