Potent but transient immunosuppression of T-cells is a general feature of CD71+ erythroid cells CD71+ erythroid cells (CECs) have been recently recognized in both neonates and cancer patients as potent immunoregulatory cells. Here, we show that in mice early-stage CECs expand in anemia, have high levels of arginase 2 (ARG2) and reactive oxygen species (ROS). In the spleens of anemic mice, CECs expansion-induced L-arginine depletion suppresses T-cell responses. In humans with anemia, CECs expand and express ARG1 and ARG2 that suppress T-cells IFN-{gamma} production. Moreover, bone marrow CECs from healthy human donors suppress T-cells proliferation. CECs differentiated from peripheral blood mononuclear cells potently suppress T-cell activation, proliferation, and IFN-{gamma} production in an ARG- and ROS-dependent manner. These effects are the most prominent for early-stage CECs (CD71highCD235adim cells). The suppressive properties disappear during erythroid differentiation as more differentiated CECs and mature erythrocytes lack significant immunoregulatory properties. Our studies provide a novel insight into the role of CECs in the immune response regulation.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=179 SRC="FIGDIR/small/427109v2_ufig1.gif" ALT="Figure 1">
View larger version (26K):
org.highwire.dtl.DTLVardef@136fbd8org.highwire.dtl.DTLVardef@d05d00org.highwire.dtl.DTLVardef@189baceorg.highwire.dtl.DTLVardef@1472b2b_HPS_FORMAT_FIGEXP  M_FIG C_FIG
Multispecific DARPin(R) therapeutics demonstrate very high potency against SARS-CoV-2 variants in vitro SARS-CoV-2 has infected millions of people globally and continues to undergo evolution. Emerging variants can be partially resistant to vaccine induced immunity and therapeutic antibodies, emphasizing the urgent need for accessible, broad-spectrum therapeutics. Here, we report a comprehensive study of ensovibep, the first trispecific clinical DARPin candidate, that can simultaneously engage all three units of the spike protein trimer to potently inhibit ACE2 interaction, as revealed by structural analyses. The cooperative binding of the individual modules enables ensovibep to retain inhibitory potency against all frequent SARS-CoV-2 variants, including Omicron BA.1 and BA.2, as of February 2022. Moreover, viral passaging experiments show that ensovibep, when used as a single agent, can prevent development of escape mutations comparably to a cocktail of monoclonal antibodies (mAb). Finally, we demonstrate that the very high in vitro antiviral potency also translates into significant therapeutic protection and reduction of pathogenesis in Roborovski dwarf hamsters infected with either the SARS-CoV-2 wild-type or the Alpha variant. In this model, ensovibep prevents fatality and provides substantial protection equivalent to the standard of care mAb cocktail. These results support further clinical evaluation and indicate that ensovibep could be a valuable alternative to mAb cocktails and other treatments for COVID-19.
HYBRID REACTION-DIFFUSION AND CLOCK-AND-WAVEFRONT MODEL FOR THE ARREST OF OSCILLATIONS IN THE SOMITOGENESIS SEGMENTATION CLOCK AO_SCPLOWBSTRACTC_SCPLOWThe clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model, and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Somitogenesis, the process by which somites are formed, is an essential developmental stage in many vertebrates. This process occurs with a strikingly regular periodicity, that is preserved among embryos of a single species. The clock and wavefront paradigm is arguably the most widely accepted model for explaining somitogenesis. However, it is incapable of explaining some experimental facts, like the appearance of somites in the absence of an external wavefront (i.e. a receding morphogen gradient). Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation, and called it a progressive oscillatory reaction-diffusion (PORD) model. The Meinhardt-PORD model tackles some of the deficiencies of the clock and wave-front models, but it has some serious issues of its own. In the present work, we introduce an equivalent Meinhardt-PORD model, and then amend it to couple it with a receding morphogen gradient. By doing so, we get a hybrid model that incorporates characteristics of the Meinhardt-PORD and clock-and-wavefront models. We show that this hybrid model under-goes a bifurcation, from a stable to an unstable limit cycle, as the value of the parameter accounting for a background regulatory input (associated to the receding morphogen gradient) decreases. This bifurcation allows the model to explain why somites can form in the absence of an external wavefront, reassesses the role of the receding morphogen gradient as a conductor for somitogenesis, and makes the model behavior robust to random fluctuations, as well as independent from specific initial conditions (the latter, are two of the weak points of the Meinhardt-PORD model). We argue that this findings provide convincing evidence that reaction-diffusion and positional information (receding morphogen gradient) mechanisms could work together in somitogenesis. Herbivore life histories are altered by drought stress in their hosts plants Drought associated with climate change can stress plants, altering their interactions with phytophagous arthropods. Dryness not only impacts cultivated plants but also their parasites, which in some cases are favoured by drought. Herbivorous arthropods feeding on drought-stressed plants typically produce bigger offspring and develop faster. However, how much responses to abiotic factors differ among populations of a species remains poorly documented. Here, we evaluate variability among populations of a major agricultural pest, the two spotted spider mite, Tetranychus urticae, in response to drought-stress. We compare key life-history parameters of twelve populations that originate from climates ranging from wet and cool Atlantic locations to medium to dry hot Mediterranean locations. We evaluated how plant drought stress affects four life history traits: development time, fecundity, sex-ratio and emigration rate in an experiment comparing well-watered and drought-stressed bean plants. Mites feeding on drought-stressed plants developed faster and attempted to leave leaves less often, and young females were more fecund. The mites from wet temperate climates exhibited greater plasticity between the two water regimes than mites originating from dryer and hot climates. Thus, climate in the area of origin influences mite response to drought. GCN5L1 promotes diastolic dysfunction by inhibiting cardiac pyruvate oxidation Left ventricular diastolic dysfunction is a structural and functional condition that precedes the development of heart failure with preserved ejection fraction (HFpEF). The etiology of diastolic dysfunction includes alterations in fuel substrate metabolism that negatively impact cardiac bioenergetics, and may precipitate the eventual transition to heart failure. To date, the molecular mechanisms that regulate early changes in fuel metabolism leading to diastolic dysfunction remain unclear. In this report, we use a diet-induced obesity model and quantitative acetylproteomics in aged mice to show that inhibitory lysine acetylation of the pyruvate dehydrogenase (PDH) complex promotes energetic deficits and diastolic dysfunction in mouse hearts. Cardiomyocyte-specific deletion of the mitochondrial lysine acetylation regulatory protein GCN5L1 prevented hyperacetylation of the PDH complex subunit PDHA1, allowing aged obese mice to continue using pyruvate as a bioenergetic substrate in the heart. Our findings suggest that changes in mitochondrial protein lysine acetylation represent a key metabolic component of diastolic dysfunction that precedes the development of heart failure. Cis-Regulatory Element Hijacking by Structural Variants Overshadows Higher-Order Topological Changes in Primary Prostate Cancer Prostate cancer is a heterogeneous disease whose progression is linked to genome instability. However the impact of this instability on the three-dimensional chromatin organization and how this drives progression is unclear. Using primary benign and tumour tissue, we find a high concordance in the higher-order three-dimensional genome organization across normal and prostate cancer cells. This concordance argues for constraints to the topology of prostate tumour genomes. Nonetheless, we identify changes to focal chromatin interactions and show how structural variants can induce these changes to guide cis-regulatory element hijacking. Such events result in opposing differential expression on genes found at antipodes of rearrangements. Collectively, our results argue that cis-regulatory element hijacking from structural variant-induced altered focal chromatin interactions overshadows higher-order topological changes in the development of primary prostate cancer. Experimental sexual selection affects the evolution of physiological and life history traits Sexual selection and sexual conflict are expected to affect all aspects of the phenotype, not only traits that are directly involved in reproduction. Here, we show coordinated evolution of multiple physiological and life history traits in response to long-term experimental manipulation of the mating system in populations of Drosophila pseudoobscura. Development time was extended under polyandry relative to monogamy in both sexes, potentially due to higher investment in traits linked to sexual selection and sexual conflict. Individuals (especially males) evolving under polyandry had higher metabolic rates and locomotor activity than those evolving under monogamy. Polyandry individuals also invested more in metabolites associated with increased endurance capacity and efficient energy metabolism and regulation, namely lipid and glycogen. Finally, polyandry males were less desiccation- and starvation-resistant than monogamy males, suggesting trade-offs between resistance and sexually selected traits. Our results provide experimental evidence that mating systems can impose selection that influences the evolution of non-sexual phenotypes such as development, activity, metabolism, and nutrient homeostasis. Decreased Interfacial Dynamics Caused by the N501Y Mutation in the SARS-CoV-2 S1 Spike:ACE2 Complex Corona Virus Disease of 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome Corona Virus 2 (SARS-CoV-2) has caused a massive health crisis across the globe, with some genetic variants gaining enhanced infectivity and competitive fitness, and thus significantly aggravating the global health concern. In this regard, the recent SARS-CoV-2 alpha variant, B.1.1.7 lineage, reported from the United Kingdom (UK), is of great significance in that it contains several mutations that increase its infection and transmission rates as evident from clinical reports. Specifically, the N501Y mutation in the SARS-CoV-2 spike S1 receptor binding domain (S1-RBD) has been shown to possess an increased affinity for ACE2, although the basis for this is not entirely clear yet. Here, we dissect the mechanism underlying the increased affinity using molecular dynamics (MD) simulations of the available ACE2-S1-RBD complex structure (6M0J) and show a prolonged and stable interfacial interaction of the N501Y mutant S1-RBD with ACE2 compared to the wild type S1-RBD. Additionally, we find that the N501Y mutant S1-RBD displays altered dynamics that likely aids in its enhanced interaction with ACE2. By elucidating a mechanistic basis for the increased affinity of the N501Y mutant S1-RBD for ACE2, we believe that the results presented here will aid in developing therapeutic strategies against SARS-CoV-2 including designing drugs targeting the ACE2-S1-RBD interaction.

SignificanceThe emergence of the new SARS-CoV-2 lineage in the UK in December 2020 has further aggravated the COVID-19 pandemic due to an increased ability of the variant to infect human hosts, likely due to mutations in the viral S1 spike protein including the N501Y S1-RBD mutation that is located at the interface of S1-RBD and ACE2, the host cell receptor for SARS-CoV-2. Given its location at the interface, N501Y S1-RBD mutation can therefore potentially alter the interfacial interaction. Multiple, all-atom, explicit solvent MD simulations of the ACE2-S1-RBD complex carried here indicated a more stable interaction between the N501Y mutant S1-RBD and ACE2 through stabilizing interfacial interactions of residues at one end of the interface that are either sequentially or physically near the mutation site. These mechanistic details will aid in better understanding the mechanism by which the alpha variant has increased infectivity as well as in designing better therapeutics including ACE2-S1 spike protein inhibitors that will, in turn, help thwarting the current and future pandemic.

HighlightsO_LIN501 in the wild type SARS-CoV-2 S1-RBD forms unsustained hydrogen bonds with residues in the ACE2, namely Y41 and K353
C_LIO_LIY501 in the N501Y mutant SARS-CoV-2 S1-RBD is not capable of forming substantial hydrogen bonds with ACE2 within the time span of the current simulation
C_LIO_LIEvidence from analyzing the simulation results suggests that Y501 of S1-RBD could form other types of non-covalent interactions with ACE2, such as van der Waals interactions
C_LIO_LIN501Y S1-RBD mutation stabilizes the position of interfacial residues neighboring to the mutation site, as well as other non-interfacial residues that are distant from the mutation site
C_LIO_LIThese altered dynamics results in more stable interaction of S1-RBD with ACE2 which could be the main reason underlying the reported enhanced affinity of S1-RBD in the SARS-CoV-2 alpha variant (UK B.1.1.7 lineage) to ACE2
C_LI Comparison of hybrid learning and remote education in the implementation of the 'Adopt a Microorganism' methodology The Internet has changed the way teachers and students access information and build knowledge. The recent COVID-19 pandemic has created challenges for both teachers and students and a demand for new methodologies of remote learning. In the life sciences, mixing online content with practical activities represents an even greater challenge. In microbiology, the implementation of an active teaching methodology, the #Adopt project, based on the social network Facebook(R), represents an excellent option for connecting remote education with classroom activities. In 2020, the version applied in high school, "Adopt a Microorganism", was adapted to meet the demands of emergency remote education owing to the suppression of face-to-face activities caused by the pandemic. In the present study, we assessed how the change in methodology impacted the discourse richness of students from high school integrated with technical education in the Business Administration program of the Federal Institute of Sao Paulo, Sorocaba Campus. Three questionnaires related to the groups of microorganisms (Archaea, Bacteria, Virus, Fungi, and Protozoan) were applied. The students responses in the 2019 and 2020 classes were compared concerning content richness and multiplicity of concepts through the application of the Shannon diversity index, an approach that is generally used to assess biodiversity in different environments. The observed results suggest that remote learning provided students with a conceptual basis and richness of content equivalent to that achieved by students subjected to the hybrid teaching model. In conclusion, this study suggests that the #Adopt project methodology increases students discourse richness in microbiology even without face-to-face traditional classes. Behavioral effects of ethanol in the Red Swamp Crayfish (Procambarus clarkii) Alcohol abuse and dependence remains one of the primary preventable sources of human mortality in the United States. Model laboratory species can be used to evaluate behavioral, neuropharmacological and other biological changes associated with alcohol exposure and to identify novel treatment modalities. This report describes methods for evaluating the behavioral effects of ethyl alcohol (ethanol; EtOH) in a crayfish model. Crayfish (Procambarus clarkii) were immersed in ethanol baths with concentrations ranging from 0.1 to 1.0 molar, and for durations of 10-30 minutes. Studies evaluated hemolymph alcohol concentration, locomotor behavior in an open field and anxiety-like behavior using a Light/Dark transfer approach. EtOH immersion produced dose dependent increases in hemolymph EtOH concentration (up to 249 mg/dL) and reductions in open field locomotor behavior that depended on EtOH concentration or exposure duration. Under baseline conditions, crayfish exhibit avoidance of the open parts of the locomotor arena and a preference for the covered portion, when available. Acute EtOH immersion decreased time spent in the covered portion of the Light/Dark arena, consistent with a decrease in anxiety-like behavior. EtOH immersion once per day for 5 days did not alter the acute locomotor behavioral response, however increased activity was observed 3 days after the repeated EtOH regimen. Overall, this study shows that this inexpensive, easily maintained species can be used for behavioral pharmacological experiments designed to assess the acute and repeated effects of EtOH. Tourmaline: a containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake BackgroundAmplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

FindingsWe developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files--a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences--it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm, assigns taxonomy to the resulting representative sequences, performs analyses of taxonomic, alpha, and beta diversity, and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S rRNA amplicon dataset from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

ConclusionsAutomated workflows like Tourmaline enable rapid analysis of environmental and biomedical amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline. Contribution of arm movements to balance recovery after tripping in older adults. Falls are common in daily life, often caused by trips and slips and, particularly in older adults, with serious consequences. Although arm movements play an important role in balance control, there is limited research into the role of arm movements during balance recovery after tripping in older adults. We investigated how older adults use their arms to recover from a trip and the difference in the effects of arm movements between fallers (n=5) and non-fallers (n=11).

Sixteen older males and females (69.7{+/-}2.3 years) walked along a walkway and were occasionally tripped over suddenly appearing obstacles. We analysed the first trip using a biomechanical model based on full-body kinematics and force-plate data to calculate whole body orientation during the trip and recovery phase. With this model, we simulated the effects of arm movements at foot-obstacle impact and during trip recovery on body orientation.

Apart from an increase in sagittal plane forward body rotation at touchdown in fallers, we found no significant differences between fallers and non-fallers in the effects of arm movements on trip recovery. Like earlier studies in young adults, we found that arm movements during the recovery phase had most favourable effects in the transverse plane: by delaying the transfer of angular momentum of the arms to the body, older adults rotated the tripped side more forward thereby allowing for a larger recovery step. Older adults that are prone to falling might improve their balance recovery after tripping by learning to prolong ongoing arm movements. A chromosome-level assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and improves time of divergence estimates for the Cannabaceae family O_LIWe present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly resulted in fragmented assemblies. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur.
C_LIO_LIThe chromosome-level assembly of Cascade was developed by scaffolding the previously-reported Cascade assembly generated with PacBio long-read sequencing, and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations, and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assess molecular evolution in gene sequences, gene family expansion and contraction, and time divergence using Bayesian inference.
C_LIO_LIWe identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families.
C_LIO_LIThe Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.
C_LI Extended Graphical Lasso for Multiple Interaction Networks for High Dimensional Omics Data There has been a spate of interest in association networks in biological and medical research, for example, genetic interaction networks. In this paper, we propose a novel method, the extended joint hub graphical lasso (EDOHA), to estimate multiple related interaction networks for high dimensional omics data across multiple distinct classes. To be specific, we construct a convex penalized log likelihood optimization problem and solve it with an alternating direction method of multipliers (ADMM) algorithm. The proposed method can also be adapted to estimate interaction networks for high dimensional compositional data such as microbial interaction networks. The performance of the proposed method in the simulated studies shows that EDOHA has remarkable advantages in recognizing class-specific hubs than the existing comparable methods. We also present three applications of real datasets. Biological interpretations of our results confirm those of previous studies and offer a more comprehensive understanding of the underlying mechanism in disease.

Author summaryReconstruction of multiple association networks from high dimensional omics data is an important topic, especially in biology. Previous studies focused on estimating different networks and detecting common hubs among all classes. Integration of information over different classes of data while allowing difference in the hub nodes is also biologically plausible. Therefore, we propose a method, EDOHA, to jointly construct multiple interaction networks with capacity in finding different hub networks for each class of data. Simulation studies show the better performance over conventional methods. The method has been demonstrated in three real world data. The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. 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. The quantitative genetics of the prevalence of infectious diseases: hidden genetic variation due to Indirect Genetic Effects dominates heritable variation and response to selection Infectious diseases have profound effects on life, both in nature and agriculture. Despite the availability of well-established epidemiological theory, however, a quantitative genetic theory of the host population for the endemic prevalence of infectious diseases is almost entirely lacking. While several studies have demonstrated the relevance of the transmission dynamics of infectious diseases for heritable variation and response to selection of the host population, our current theoretical framework of quantitative genetics excludes these dynamics. As a consequence, we do not know which genetic effects of the host population determine the prevalence of an infection, and have no concepts of breeding value and heritable variation for endemic prevalence.

Here we integrate quantitative genetics and epidemiology, and propose a quantitative genetic theory for R0 and for the endemic prevalence of an infectious disease. We first identify the genetic factors that determine the prevalence of an infection, using an approach founded in epidemiological theory. Subsequently we investigate the population level consequences of individual genetic variation, both for R0 and for the endemic prevalence. Next, we present expressions for the breeding value and heritable variation, for both prevalence and individual binary disease status, and show that these parameters depend strongly on the level of the prevalence. Results show that heritable variation for endemic prevalence is substantially greater than currently believed, and increases when prevalence approaches zero, while heritability of individual disease status goes to zero. As a consequence, response of prevalence to selection accelerates considerably when prevalence goes down, in contrast to predictions from classical theory. Finally, we show that most of the heritable variation for the endemic prevalence of an infection is hidden due to indirect genetic effects, suggesting a key role for kin-group selection both in the evolutionary history of current populations and for genetic improvement strategies in animals and plants. Cluster-specific gene markers enhance Shigella and Enteroinvasive Escherichia coli in silico serotyping Shigella and enteroinvasive Escherichia coli (EIEC) cause human bacillary dysentery with similar invasion mechanisms and share similar physiological, biochemical and genetic characteristics. The ability to differentiate Shigella and EIEC from each other is important for clinical diagnostic and epidemiologic investigations. The existing genetic signatures may not discriminate between Shigella and EIEC. However, phylogenetically, Shigella and EIEC strains are composed of multiple clusters and are different forms of E. coli. In this study, we identified 10 Shigella clusters, 7 EIEC clusters and 53 sporadic types of EIEC by examining over 17,000 publicly available Shigella/EIEC genomes. We compared Shigella and EIEC accessory genomes to identify the cluster-specific gene markers or marker sets for the 17 clusters and 53 sporadic types. The gene markers showed 99.63% accuracy and more than 97.02% specificity. In addition, we developed a freely available in silico serotyping pipeline named Shigella EIEC Cluster Enhanced Serotype Finder (ShigEiFinder) by incorporating the cluster-specific gene markers and established Shigella/EIEC serotype specific O antigen genes and modification genes into typing. ShigEiFinder can process either paired end Illumina sequencing reads or assembled genomes and almost perfectly differentiated Shigella from EIEC with 99.70% and 99.81% cluster assignment accuracy for the assembled genomes and mapped reads respectively. ShigEiFinder was able to serotype over 59 Shigella serotypes and 22 EIEC serotypes and provided a high specificity with 99.40% for assembled genomes and 99.38% for mapped reads for serotyping. The cluster markers and our new serotyping tool, ShigEiFinder (https://github.com/LanLab/ShigEiFinder), will be useful for epidemiologic and diagnostic investigations.

Impact statementThe differentiation of Shigella strains from enteroinvasive E. coli (EIEC) is important for clinical diagnosis and public health epidemiologic investigations. The similarities between Shigella and EIEC strains make this differentiation very difficult as both share common ancestries within E. coli. However, Shigella and EIEC are phylogenetically separated into multiple clusters, making high resolution separation using cluster specific genomic markers possible. In this study, we identified 17 Shigella or EIEC clusters including five that were newly identified through examination of over 17,000 publicly available Shigella and EIEC genomes. We further identified an individual or a set of cluster-specific gene markers for each cluster using comparative genomic analysis. These markers can then be used to classify isolates into clusters and were used to develop an in silico pipeline, ShigEiFinder (https://github.com/LanLab/ShigEiFinder) for accurate differentiation, cluster typing and serotyping of Shigella and EIEC from Illumina sequencing reads or assembled genomes. This study will have broad application from understanding the evolution of Shigella/EIEC to diagnosis and epidemiology.

Data summarySequencing data have been deposited at the National Center for Biotechnology Information under BioProject number PRJNA692536.

RepositoriesRaw sequence data are available from NCBI under the BioProject number PRJNA692536. Supervised Learning Model Predicts Protein Adsorption to Carbon Nanotubes Engineered nanoparticles are advantageous for numerous biotechnology applications, including biomolecular sensing and delivery. However, testing the compatibility and function of nanotechnologies in biological systems requires a heuristic approach, where unpredictable biofouling via protein corona formation often prevents effective implementation. Moreover, rational design of biomolecule-nanoparticle conjugates requires prior knowledge of such interactions or extensive experimental testing. Toward better applying engineered nanoparticles in biological systems, herein, we develop a random forest classifier (RFC) trained with proteomic mass spectrometry data that identifies proteins that adsorb to nanoparticles, based solely on the proteins amino acid sequence. We model proteins that populate the corona of a single-walled carbon nanotube (SWCNT)-based optical nanosensor and study whether there is a relationship between the proteins amino acid-based properties and the proteins adsorption to SWCNTs. We optimize the classifier and characterize the classifier performance against other models. To evaluate the predictive power of our model, we apply the classifier to rapidly identify proteins with high binding affinity to SWCNTs, followed by experimental validation. We further determine protein features associated with increased likelihood of SWCNT binding: high content of solvent-exposed glycine residues and non-secondary structure-associated amino acids. Conversely, proteins with high content of leucine residues and beta-sheet-associated amino acids are less likely to form the SWCNT protein corona. The classifier presented herein provides a step toward undertaking the otherwise intractable problem of predicting protein-nanoparticle interactions, which is needed for more rapid and effective translation of nanobiotechnologies from in vitro synthesis to in vivo use.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=66 SRC="FIGDIR/small/449132v2_ufig1.gif" ALT="Figure 1">
View larger version (22K):
org.highwire.dtl.DTLVardef@bbc175org.highwire.dtl.DTLVardef@9a0b94org.highwire.dtl.DTLVardef@16e22e1org.highwire.dtl.DTLVardef@1b40a84_HPS_FORMAT_FIGEXP  M_FIG C_FIG Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

MethodsWe studied the dynamic changes of ROCK and actomyosin cytoskeletal networks during NETosis induction and their involvement in PKC nuclear translocation. Using mice with hematopoietic-specific ROCK1 deficiency, we investigated the effects of ROCK1 deficiency on NETosis, and its involvement in UVB-induced skin inflammation.

ResultsOur time course studies demonstrated the dynamic changes of actin polymerization and ROCK activation, support the role of actin cytoskeleton in nuclear translocation of cytosolic PKC in early stage of NETosis induction. Inhibition of actin polymerization or key molecules of the ROCK/MLCK/myosin pathway decreased PKC nuclear translocation and NET formation. Genetic deficiency of ROCK1, inhibited NETosis ex vivo and in vivo, decreased extracellular display of NET-associated IL-17A, TNF, IFN{gamma}, and IFN in inflamed skin, which were correlated with the ameliorated skin inflammation in UVB-irradiated mice with hematopoietic-specific ROCK1 deficiency.

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Temperature drives variation in flying insect biomass across a German malaise trap network O_LIAmong the many concerns for biodiversity in the Anthropocene, recent reports of flying insect loss are particularly alarming, given their importance as pollinators and as a food source for many predators. Few insect monitoring programs cover large spatial scales required to provide more generalizable estimates of insect responses to global change drivers.
C_LIO_LIWe ask how climate and surrounding habitat affect flying insect biomass and day of peak biomass using data from the first year of a new standardized distributed monitoring network at 84 locations across Germany comprising spatial gradient of land-cover types from protected to urban areas.
C_LIO_LIFlying insect biomass increased linearly with monthly temperature across Germany. However, the effect of temperature on flying insect biomass flipped to negative in the hot months of June and July when local temperatures most exceeded long-term averages.
C_LIO_LILand-cover explained little variation in insect biomass, but biomass was lowest in forested sites. Grasslands, pastures and orchards harbored the highest insect biomass. The date of peak biomass was primarily driven by surrounding land-cover type, with grasslands especially having earlier insect biomass phenologies.
C_LIO_LIStandardized, large-scale monitoring is pivotal to uncover underlying processes of insect decline and to develop climate-adapted strategies to promote insect diversity. In a temperate climate region, we find that the benefits of temperature on flying insect biomass diminish in a German summer at locations where temperatures most exceeded long-term averages. These results highlighting the importance of local adaptation in climate change-driven impacts on insect communities.
C_LI Region-specific microRNA alterations in marmosets carrying SLC6A4 polymorphisms are associated with anxiety-like behavior. Psychiatric diseases such as depression and anxiety are multifactorial conditions, highly prevalent in western societies. Human studies have identified a number of high-risk genetic variants for these diseases. Among them, polymorphisms in the promoter region of the serotonin transporter gene (SLC6A4) have attracted much attention. However, due to the paucity of experimental models, molecular alterations induced by these genetic variants and how they correlate to behavioral deficits have not been examined. Marmosets have emerged as a powerful model in translational neuroscience to investigate molecular underpinnings of complex behaviors. Here, we took advantage of naturally occurring genetic polymorphisms in marmoset SLC6A4 gene that have been linked to anxiety-like behaviors. Using FACS-sorted cells from different brain regions, we revealed that marmosets bearing different SLC6A4 variants exhibit distinct microRNAs signatures in a region of the prefrontal cortex whose activity has been consistently altered in patients with depression/anxiety. We also identified DCC, a gene previously linked to these diseases, as a downstream target of the dysregulated microRNAs. Significantly, we showed that levels of both microRNAs and DCC in this region were highly correlated to anxiety-like behaviors as well as to the response to citalopram, a selective serotonin re-uptake inhibitor and widely prescribed anti-depressant. Our findings establish links between genetic variants, molecular modifications in specific cortical regions and complex behavioral/pharmacological responses, providing new insights into gene-behavior relationships underlying human psychopathology. Pseudomonas aeruginosa C-terminal processing protease CtpA assembles into a hexameric structure that requires activation by a spiral-shaped lipoprotein binding partner Pseudomonas aeruginosa CtpA is a carboxyl terminal-processing protease that partners with the outer membrane lipoprotein LbcA to degrade at least five cell wall-associated proteins, four of which are cell wall hydrolases. This activity plays an important role in supporting P. aeruginosa virulence in a mouse model of acute pneumonia. However, almost nothing is known about the molecular mechanisms underlying CtpA and LbcA function. Here, we used structural analysis to show that CtpA alone assembles into an inactive hexamer comprising a trimer of dimers, which limits its substrate access and prevents nonspecific degradation. The adaptor protein LbcA is a right-handed open spiral with 11 tetratricopeptide repeats, which might wrap around a substrate to deliver it to CtpA for degradation. By structure-guided mutagenesis and functional assays, we also showed that the interfaces of the CtpA trimer-of-dimers, and an N-terminal helix of LbcA, are important for LbcA-mediated substrate degradation by CtpA both in vitro and in vivo. This work improves our understanding of the molecular mechanism of a CTP within the C-terminal processing peptidase-3 group.

IMPORTANCECarboxyl-terminal processing proteases (CTPs) are found in all three domains of life. In bacteria, some CTPs have been associated with virulence, raising the possibility that they could be theraputic targets. However, relatively little is known about their molecular mechanisms of action. In Pseudomonas aeruginosa, CtpA supports virulence by working in complex with the outer membrane lipoprotein LbcA to degrade cell wall hydrolases. Here, we report structure-function analyses of CtpA and LbcA, which reveals that CtpA assembles into an inactive hexamer comprising a trimer of dimers. LbcA is monomeric, with an N-terminal region important for binding to and activating CtpA, followed by a spiral structure composed of 11 tetratricopetide repeats, which could wrap around a substrate for delivery to CtpA. This work provides the first structure of a CTP-3 group member, revealing a unique mutimeric arrangement and insight into how this important proteolytic system functions. New insights into tomato CLE peptide repertoire and perception mechanisms Precision in sensing the environmental cues and adjusting the growth and the physiology of the root system are necessary for plant robustness. Plants achieve their phenotypic plasticity by tightly controlling and buffering developmental decisions. In addition to the classical plant hormones that mediate plant development and stress responses, the CLE peptides constitute an additional crucial level of regulation. While the CLV3-CLV1 module appears to be highly conserved to control the proliferation of the shoot apical meristem stem cells, we do not yet fully understand the function of the additional CLE genes and whether they act in a similar way across the plant species, including tomato. Due to the small gene size and high sequence variability, it is extremely difficult to precisely annotate CLE genes in plant genomes. Here we present our analysis of the CLE family in tomato, based on a combination of iterative tBLASTn and Hidden-Markov-Model (HMM), which allowed us to identify thirty-seven new SlCLEs in addition to the fifteen reported previously. We could confirm the biological activities of selected SlCLEs in suppressing root meristematic cell divisions. We show that root response is mediated by SlCLAVATA2, indicating the conservation of CLE perception mechanism.

One-sentence summaryUsing a combination of iterative tBLASTn and Hidden-Markov-Model approaches, we uncovered 37 new tomato CLE genes predominantly expressed in roots, and we showed a conserved effect on root meristem arrest, that was SlCLAVATA2-dependent.

HighlightsO_LIWe applied a combined approach of iterative tBLASTn and Hidden-Markov-Model to identify fifty-two tomato SlCLE genes, including thirty-seven new genes
C_LIO_LIAll identified genes encode for pre-propeptides with a single CLE-domain containing conserved residues similar to Arabidopsis
C_LIO_LIAnalyzing the publicly available RNAseq datasets, we could confirm the expression of SlCLE genes that was often associated specifically with root or shoot, a certain developmental stage of the fruit, or with drought stress conditions
C_LIO_LIRemarkably, the majority of SlCLE genes are predominantly expressed in the root tissues
C_LIO_LIWe showed the conserved inhibitory effect on the root meristem and columella cells division for the selected SlCLE peptides that were SlCLAVATA2-dependent.
C_LI Meiotic DSB-independent role of protein phosphatase 4 in Hop1 assembly to promote meiotic chromosome axis formation in budding yeast Dynamic changes in chromosomal structure that occur during meiotic prophase play an important role in the progression of meiosis. Among these structures, meiosis-specific chromosomal axis-loop structures are important as a scaffold for integrated control between the meiotic recombination reaction and the associated checkpoint system to ensure accurate chromosome segregation. However, the molecular mechanism of the initial step of chromosome axis-loop construction is not well understood. Here, we showed that in budding yeast, protein phosphatase 4 (PP4) that primarily counteracts Mec1/Tel1 phosphorylation is required to promote the assembly of chromosomal axis components Hop1 and Red1 onto meiotic chromatin via interaction with Hop1. PP4, in contrast, did not affect Rec8 assembly. Notably, unlike the previously known function of PP4, this novel function of PP4 was independent of meiotic DSB-dependent Tel1/Mec1 kinase activities. The defect in Hop1/Red1 assembly in the absence of PP4 function was not suppressed by dysfunction of Pch2, which removes Hop1 protein from the chromosome axis. This suggests that PP4 is required for the initial step of chromatin loading of Hop1 rather than stabilization of Hop1 on axes. These results indicate that phosphorylation/dephosphorylation-mediated regulation of Hop1 recruitment onto chromatin occurs during chromosome axis construction before meiotic double-strand break formation. ATP synthase is essential for high conductance but not low conductance Mitochondrial Permeability Transition A sudden increase in permeability of the mitochondrial inner membrane, known as mitochondrial permeability transition (PT), is the central event responsible for cell death and tissue damage in conditions such as stroke and heart attack. PT is caused by the opening of the Cyclosporin A (CSA) dependent calcium-induced large pore, the Permeability Transition Pore (PTP). The molecular details of PTP opening remain incompletely understood. We utilized holographic microscopy imaging to directly assess PTP induction in HAP1 living cells. We found that calcium stress caused CSA-sensitive depolarization of the mitochondrial inner membrane. PTP did not occur until after nearly complete mitochondrial membrane depolarization. In cells lacking c and {delta} subunits of the ATP synthase, we observed CSA-sensitive membrane depolarization but not high-conductance PTP. We propose that PTP activation occurs as a two-phase process, with the initial membrane depolarization followed by the second phase of large pore opening that requires ATP synthase.