Shorter effective lifespan in laboratory populations of D. melanogaster might reduce sexual selection. The role of sexual selection in mediating levels of sexual conflict has been demonstrated in many experimental evolution studies on Drosophila spp. where competition among males for mating was the target of selection. Sexual selection has also been shown to affect the evolution of life-histories. However, the influence of divergent life-histories on reproductive strategies and, therefore, sexual selection and possibly sexual conflict, has been less well studied. We examined D. melanogaster populations selected for a short development time and early age at reproduction for changes in reproductive behaviour and traits that are proxies of sexual selection. We report a large reduction in reproductive competition experienced by the males of these populations, compared to ancestral populations that are not consciously selected for rapid development or early reproduction, potentially leading to reduced sexual selection. We show that rapidly developing and early reproducing populations have very low levels of mating in their lifetime (females are more or less monandrous), low courtship levels, shorter copulation duration, and longer time from eclosion to first mating, compared to the controls. These results are discussed in the context of the previously demonstrated reduction of inter-locus sexual conflict in these populations. We show that life-history strategies might have a large and significant impact on sexual selection, with each influencing the other and contributing to the complexities of adaptation.

Significance statementSexual conflict, often manifested as an arms-race between males and females trying to enhance their own reproductive success at some cost to the other, is of great evolutionary interest because it can maintain genetic variation in populations, prevent the independent optimization of male and female traits, and also promote speciation. Sexual selection, or variation in mating success, is well known to affect levels of sexual conflict. However, it is not so clear whether, and how, the regular evolution of life-histories also affects sexual selection. Here, we show that life-history evolution in fruit fly populations selected for traits not directly related to sexual conflict might, nevertheless, mediate the possible evolution of altered sexual conflict levels through effects on sexual selection. Populations that evolved to develop to adulthood fast, and reproduce relatively early in life, are shown to potentially experience less sexual selection, which can explain the low sexual conflict levels earlier observed in them.
Phylogenomic analyses of the East Asian endemic Abelia (Caprifoliaceae) shed insights into the temporal and spatial diversification history with widespread hybridizations Background and AimsAbelia (Caprifoliaceae) is a small genus with five species, including one man-made hybrid and several natural hybrids. The genus has a discontinuous distribution in mainland China, Taiwan Island, and the Ryukyu islands, providing a model system to explore mechanisms of species dispersal in the East Asian flora. However, the current phylogenetic relationships within Abelia remain uncertain.

MethodsIn this study, we reconstructed phylogenetic relationships within Abelia using nuclear loci generated by target enrichment and plastomes from genome skimming. Divergence time estimation, ancestral area reconstruction, and ecological niche modelling (ENM) were used to examine the diversification history of Abelia.

Key ResultsWe found extensive cytonuclear discordance across the genus. By integrating lines of evidence from molecular phylogenies, divergence times, and morphology, we propose to merge A. macrotera var. zabelioides into A. uniflora. Network analyses suggested that there have been widespread and multiple hybridization events among Abelia species. These hybridization events may have contributed to the speciation mechanism and resulted in a high observed morphological diversity. The diversification of Abelia began in the early Eocene, followed by A. chinensis var. ionandra colonizing the island of Taiwan in the Middle Miocene. The ENM results suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last Interglacial. Disjunction between the Himalayan-Hengduan Mountain region (HHM) and the island of Taiwan is most likely the consequence of topographic isolation and postglacial contraction.

ConclusionsWe used genomic data to reconstruct the phylogeny of Abelia and found a clear pattern of reticulate evolution in the group. In addition, our results support shrinkage of postglacial range and the heterogeneity of the terrain have led to the disjunction of the mainland China-Taiwan island. This study provides important new insights into the speciation process and taxonomy of Abelia.
Universal protein misfolding intermediates can bypass the proteostasis network and remain soluble and less-functional Misfolded protein conformations with decreased functionality can bypass the proteostasis machinery and remain soluble in vivo. This is an unexpected phenomenon as several cellular quality control mechanisms have evolved to rid cells of misfolded proteins. Three questions, then, are: how is it structurally possible for long-lived, soluble, misfolded proteins to bypass the proteostasis machinery and processes? How widespread are these soluble, misfolded states across the proteome? And how long do they persist for? Here, we address these questions using coarse-grain molecular dynamics simulations of the synthesis, termination, and post-translational dynamics of a representative set of cytosolic E. coli proteins. We predict that half of all proteins exhibit subpopulations of misfolded conformations that are likely to bypass molecular chaperones, avoid aggregation, and not be rapidly degraded. These misfolded states may persist for months or longer for some proteins. Structurally characterizing these misfolded states, we observe they have a large amount of native structure, but also contain localized misfolded regions from non-native changes in entanglement, in which a protein segment threads through a loop formed by another portion of the protein that is not found in the native state. The surface properties of these misfolded states are native like, suggesting they may bypass the proteostasis machinery and its regulatory processes to remain soluble, while their entanglements make these states long-lived kinetic traps, as disentanglement requires unfolding of already folded portions of the protein. In terms of function, we predict that one-third of proteins have subpopulations that misfold into less-functional states that have structurally perturbed functional sites yet remain soluble. Data from limited-proteolysis mass spectrometry experiments, which interrogate the misfolded conformations populated by proteins upon unfolding and refolding, are consistent with the structural changes seen in the entangled states of glycerol-3-phosphate dehydrogenase upon misfolding. These results provide an explanation for how proteins can misfold into soluble conformations with reduced functionality that can bypass cellular quality controls, and indicate, unexpectedly, this may be a wide-spread phenomenon in proteomes. Such entanglements are observed in many native structures, suggesting the non-native entanglements we observe are plausible. More broadly, these near-native entangled structures suggest a hypothesis for how synonymous mutations can modulate downstream protein structure and function, with these mutations partitioning nascent proteins between these kinetically trapped states. Augmenting Recombinant Antibody Production in HEK293E Cells: Optimising Transfection and Culture Parameters BackgroundOptimising recombinant antibody production is important for cost-effective therapeutics and diagnostics. With impact on commercialisation, higher productivity beyond laboratory scales is highly sought, where efficient production can also accelerate antibody characterisations and investigations.

MethodsInvestigating HEK293E cells for mammalian antibody production, various transfection and culture parameters were systematically analysed for antibody light chain production before evaluating them for whole antibody production. Transfection parameters investigated include seeding cell density, the concentration of the transfection reagent and DNA, complexation time, temperature, and volume, as well as culture parameters such as medium replacement, serum deprivation, use of cell maintenance antibiotic, incubation temperature, medium volume, post-transfection harvest day and common nutrient supplements.

ResultsUsing 2 mL adherent HEK293E cell culture transfections with 25 kDa linear Polyethylenimine in the most optimised parameters, we demonstrated a ~2-fold production increase for light chain alone and for whole antibody production reaching 536 and 49 g respectively in a cost-effective manner. With the addition of peptone, {kappa} light chain increased by ~4-fold to 1032 g while whole antibody increased to a lesser extent by ~2.5-fold to 51 g, with benefits potentially for antibodies limited by their light chains in production.

ConclusionsOur optimised findings show promise for a more efficient and convenient antibody production method through transfection and culture optimisations that can be incorporated to scale up processes and with potential transferability to other mammalian-based recombinant protein production using HEK293E cells.

Statement of SignificanceRecombinant antibody production is crucial for antibody research and development. Systematically investigating transfection and culture parameters such as PEI/DNA concentrations, complexation time, volume, and temperature, supplements, etc., we demonstrated a ~4-fold light chain alone production increase to 1032 g and a 2.5-fold whole antibody production increase to 51 g from 2 mL transfections. Synthetic Cooling Agents in US-marketed E-cigarette Refill Liquids and Popular Disposable E-cigarettes: Chemical Analysis and Risk Assessment. BackgroundMenthol, through its cooling and pleasant sensory effects, facilitates smoking and tobacco product initiation, resulting in the high popularity of mint/menthol-flavored E-cigarettes. More recently, E-cigarette vendors started marketing synthetic cooling agents as additives that impart a cooling effect but lack a characteristic minty odor. Knowledge about content of synthetic coolants in US-marketed E-cigarette products and associated health risks is limited.

MethodsE-liquid vendor sites were searched with the terms "koolada", "kool/cool", "ice", or WS-3/WS-23, denoting individual cooling agents, and relevant refill E-liquids were purchased. "Ice" flavor varieties of Puffbar, the most popular disposable E-cigarette brand, were compared with non-"Ice" varieties. E-liquids were characterized, and synthetic coolants quantified using GC/MS. Margin of exposure (MOE), a risk assessment parameter, was calculated to assess the risk associated with synthetic coolant exposure from E-cigarette use.

ResultsWS-3 was detected in 24/25 refill E-liquids analyzed. All Puffbar flavor varieties contained either WS-23 (13/14) or WS-3 (5/14), in both "Ice"- and non-"Ice" flavors. Modeling consumption of WS-3 from vaped E-liquids, resulted in MOEs below the safe margin of 100 for most daily use scenarios. MOEs for WS-23 were <100 for 10/13 Puffbar flavors in all use scenarios. Puffbar power specifications are identical to Juul devices.

ConclusionsSynthetic cooling agents (WS-3/WS-23) were present in US-marketed E-cigarettes, at levels that may result in consumer exposures exceeding safety thresholds set by regulatory agencies. Synthetic coolants are not only found in mint-or menthol-flavored products, but also in fruit- and candy-flavored products, including popular disposable E-cigarette products such as Puffbar.

ImplicationsSynthetic cooling agents are widely used in "kool/cool"- and "ice"-flavored E-liquids and in E-liquids without these labels, both as a potential replacement for menthol or to add cooling  notes to non-menthol flavors. These agents may be used to bypass current and future regulatory limits on menthol content in tobacco products, and not just E-cigarettes. Since synthetic cooling agents are odorless, they may not fall under the category of "characterizing flavor", potentially circumventing regulatory measures based on this concept. Regulators need to consider the additional health risks associated with exposure to synthetic cooling agents. Long-term culture of fetal monocyte precursors in vitro allowing the generation of bona fide alveolar macrophages in vivo Tissue-resident macrophage-based immune therapies have been proposed for various diseases. However, generation of sufficient numbers that possess tissue-specific functions remains a major handicap. Here, we show that fetal liver monocytes (FLiMo) cultured with GM-CSF (also known as CSF2) rapidly differentiate into a long-lived, homogeneous alveolar macrophage (AM)-like population in vitro. CSF2-cultured FLiMo remain the capacity to develop into bona fide AM upon transfer into Csf2ra-/- neonates and prevent development of alveolar proteinosis and efferocytosis of apoptotic cells for at least 1 year in vivo. Compared to transplantation of AM-like cells derived from bone marrow macrophages (BMM), CSF2-cFliMo more efficiently engraft empty AM niches in the lung and protect mice from respiratory viral infection. Harnessing the potential of this approach for gene therapy, we restored a disrupted Csf2ra gene in FLiMo and their capacity to develop into AM in vivo. Together, we provide a novel platform for generation of immature AM-like precursors amenable for genetic manipulation, which will be useful to study to dissect AM development and function and pulmonary transplantation therapy. Spatial regulation of the anaphase-promoting complex/cyclosome by a microtubule adaptor The anaphase-promoting complex/cyclosome (APC/C) coordinates advancement through mitosis via temporally controlled polyubiquitination of effector proteins. Despite the long-appreciated spatial organization of key events in mitosis mediated largely by cytoskeletal networks, the spatial regulation of APC/C, the major mitotic E3 ligase, is poorly understood. Here, we describe a microtubule-resident protein, PLEKHA5, as an interactor of APC/C and spatial regulator of its activity in mitosis. PLEKHA5 knockdown delayed mitotic progression, causing accumulation of APC/C substrates dependent upon the PLEKHA5-APC/C interaction. To assess the spatial importance of this interaction, we developed microtubule-localized proximity biotinylation tools, which revealed that depletion of PLEKHA5 decreased the extent of APC/C association with the microtubule cytoskeleton. This decreased APC/C microtubule-localization in turn prevented efficient loading of APC/C with its co-activator CDC20, leading to defects in E3 ligase catalytic activity. We propose that PLEKHA5 functions as an adaptor of APC/C that promotes its subcellular localization to microtubules and facilitates its activation by CDC20, thus ensuring the timely turnover of key mitotic APC/C substrates and proper progression through mitosis. Methylome dynamic upon proteasome inhibition by the Pseudomonas syringae virulence factor Syringolin A DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, where modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effectors molecule, several of which act as proteasome inhibitors. Here we investigated the effect of proteasome inhibition by the bacterial virulence factor Syringolin A on genome-wide DNA methylation. We show that Syringolin A treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH DMRs that are enriched in the proximity of transcriptional start sites. Syringolin A treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway: AGO3, AGO9 and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of a epi-genomic arms race against pathogens. Remarked Suppression of Aβ42 Protomer-protomer Dissociation Reaction via Pentamer Dimerization Multimeric protein complexes are molecular apparatuses to regulate biological systems and often determine their fate. Among proteins forming such molecular assemblies, amyloid proteins have drawn attention over a half-century since amyloid fibril formation of these proteins is supposed to be a common pathogenic cause for neurodegenerative diseases. This process is triggered by the accumulation of fibril-like aggregates, while the microscopic mechanisms are mostly elusive due to technical limitation of experimental methodologies in individually observing each of diverse aggregate species in the aqueous solution. We then addressed this problem by employing atomistic molecular dynamics simulations for the paradigmatic amyloid protein, amyloid-{beta} (1-42) (A{beta}42). Seven different dimeric forms of oligomeric A{beta}42 fibril-like aggregate in aqueous solution, ranging from tetramer to decamer, were considered. We found additive effects of the size of these fibril-like aggregates on their thermodynamic stability and have clarified kinetic suppression of protomer-protomer dissociation reactions at and beyond the point of pentamer dimer formation. This observation was obtained from the specific combination of the A{beta}42 protomer structure and the physicochemical condition that we here examined, while it is worthwhile to recall that several amyloid fibrils take dimeric forms of their protomers. We could thus conclude that the stable formation of fibril-like protomer dimer should be involved in a turning point where rapid growth of amyloid fibrils is triggered. RapidoPGS: A rapid polygenic score calculator for summary GWAS data without a test dataset MotivationPolygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RapidoPGS, a flexible and fast method to compute PGS requiring summary-level GWAS datasets only, with little computational requirements and no test data required for parameter tuning.

ResultsWe show that RapidoPGS performs slightly less well than two out of three other widely-used PGS methods (LDpred2, PRScs, and SBayesR) for case-control datasets, with median r2 difference: -0.0092, -0.0042, and 0.0064, respectively, but up to 17,000-fold faster with reduced computational requirements. RapidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog.

Availability and implementationOur method is available with a GPL license as an R package from GitHub. The Resting-State Causal Human Connectome is Characterized by Hub Connectivity of Executive and Attentional Networks We demonstrate a data-driven approach for calculating a "causal connectome" of directed connectivity from resting-state fMRI data using a greedy adjacency search and pairwise non-Gaussian edge orientations. We used this approach to construct n=442 causal connectomes. These connectomes were very sparse in comparison to typical Pearson correlation-based graphs (roughly 2.25% edge density) yet were fully connected in nearly all cases. Prominent highly connected hubs of the causal connectome were situated in attentional (dorsal attention) and executive (frontoparietal and cingulo-opercular) networks. These hub networks had distinctly different connectivity profiles: attentional networks shared incoming connections with sensory regions and outgoing connections with higher cognitive networks, while executive networks primarily connected to other higher cognitive networks and had a high degree of bidirected connectivity. Virtual lesion analyses accentuated these findings, demonstrating that attentional and executive hub networks are points of critical vulnerability in the human causal connectome. These data highlight the central role of attention and executive control networks in the human cortical connectome and set the stage for future applications of data-driven causal connectivity analysis in psychiatry. Foraging with friends, moving with anyone: effects of the social environment on movement-integrated habitat selection O_LIMovement links the distribution of habitats with the social environment of animals using those habitats; yet integrating movement, habitat selection, and socioecology remains an opportunity for further study.
C_LIO_LIHere, our objective was to disentangle the roles of habitat selection and social association as drivers of collective movement in a gregarious ungulate. To accomplish this objective, we (1) assessed whether socially familiar individuals form discrete social communities and whether social communities have high spatial, but not necessarily temporal, overlap; and (2) we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection analysis.
C_LIO_LIWe used social network analysis to assign individuals to social communities and determine short and long-term social preference among individuals. Using integrated step selection functions (iSSF), we then modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.
C_LIO_LIBased on assignment of individuals to social communities and home range overlap analyses, individuals assorted into discrete social communities, and these communities had high spatial overlap. By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals foraged with more familiar individuals, but moved collectively with any between foraging patches.
C_LIO_LIOur study demonstrates that social behaviour and space use are inter-related based on spatial overlap of social communities and movement-dependent habitat selection. Movement, habitat selection, and social behaviour are linked in theory. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.
C_LI A Eukaryotic Circuit for Secrete-and-Sense Autonomy Cancers represent complex autonomous robust systems, displaying interconnectivity with feedback control. Autonomous growth is fueled by a cancer cells ability to  secrete-and-sense growth factors: a poorly understood phenomenon. Using an integrated systems and experimental approach, here we dissect the impact of a feedback-coupled GTPase circuit within the secretory pathway that imparts secretion-coupled autonomy. The circuit is assembled when the Ras-superfamily monomeric GTPase Arf1, and the heterotrimeric GTPase Giabg and their corresponding GAPs and GEFs are coupled by GIV/Girdin, a protein that is known to fuel aggressive traits in diverse cancers. One forward and two key negative feedback loops within the circuit create closed-loop control (CLC), allow the two GTPases to coregulate each other, and convert the expected switch-like behavior of Arf1-dependent secretion into an unexpected dose response alignment behavior of sensing and secretion. Such behavior translates into cell survival that is self-sustained by stimuli-proportionate secretion. Findings highlight how enhanced coupling of two biological switches in cancer cells is critical for multiscale feedback control to achieve secretion-coupled autonomy.

GRAPHIC ABSTRACT

O_FIG O_LINKSMALLFIG WIDTH=163 HEIGHT=200 SRC="FIGDIR/small/436048v3_ufig1.gif" ALT="Figure 1">
View larger version (51K):
org.highwire.dtl.DTLVardef@183ff5corg.highwire.dtl.DTLVardef@f21de7org.highwire.dtl.DTLVardef@1893683org.highwire.dtl.DTLVardef@159d966_HPS_FORMAT_FIGEXP  M_FIG C_FIG In BriefThis work delivers an experimentally validated dynamical systems model of the cooperativity between two distinct classes of biological switches in eukaryotic cells and reveals the basis of secretion-coupled cellular autonomy in cancer.

HIGHLIGHTSO_LIModeling and experimental approaches were used to dissect a coupled GTPase circuit
C_LIO_LICoupling enables closed loop feedback and mutual control of GTPases
C_LIO_LICoupling generates dose response alignment behavior of sensing and secretion.
C_LIO_LICoupling is critical for multiscale feedback control to achieve secretion-coupled autonomy.
C_LI The C. elegans gonadal sheath Sh1 cells extend asymmetrically over a differentiating germ cell population in the proliferative zone The C. elegans adult hermaphrodite germ line is surrounded by a thin tube formed by somatic sheath cells that support germ cells as they mature from the stem-like mitotic state through meiosis, gametogenesis and ovulation. Recently, we discovered that the distal-most Sh1 sheath cells associate with mitotic germ cells as they exit the niche. Here we report that these distal sheath-associated germ cells differentiate first in animals with temperature-sensitive mutations affecting germ cell state, and stem-like germ cells are maintained distal to the Sh1 boundary. We analyze several markers of the distal sheath, which is best visualized with endogenously tagged membrane proteins, as overexpressed fluorescent proteins fail to localize to distal membrane processes and can cause gonad morphology defects. However, such reagents with highly variable expression can be used to determine the relative positions of the two Sh1 cells, one of which often extends further distal than the other. pkgndep: a tool for heaviness analysis on R package dependencies SummaryNumerous R packages have been developed for bioinformatics analysis in the last decade and dependencies among packages have become critical issues to consider. In this work, we proposed a new metric named dependency heaviness that measures the number of unique dependencies that a parent brings to a package, and we proposed possible solutions for reducing the complexity of dependencies by optimizing the use of heavy parents. We implemented the metric in a new R package pkgndep which provides an intuitive way for dependency heaviness analysis. Based on pkgndep, we additionally performed a global analysis of dependency heaviness on CRAN and Bioconductor ecosystems and we revealed top packages that have significant contributions of high dependency heaviness to their child packages.

Availability and implementationThe package pkgndep and documentations are freely available from the Comprehensive R Archive Network (CRAN) https://cran.r-project.org/package=pkgndep. The dependency heaviness analysis for all 21,741 CRAN and Bioconductor packages retrieved on 2021-10-28 are available at https://pkgndep.github.io/. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Design, Mutate, Screen: High-throughput creation of genetic clocks with different period-amplitude characteristics Directed evolution has become an invaluable tool in protein engineering and has also greatly influenced the construction of synthetic gene circuits. The ability to generate diversity at precise targets for directed evolution approaches has improved vastly, allowing researchers to create large, specific mutant libraries with relative ease. Screening approaches for large mutant libraries have similarly come a long way, especially when the desired behavior can easily be tested for with static, single time-point assays. For more complex gene circuits with dynamic phenotypes that change over time, directed evolution approaches to controlling and tuning circuit behavior have been hindered by the lack of sufficiently high-throughput screening methods to isolate variants with desired characteristics. Here we utilize directed mutagenesis and multiplexed microfluidics to develop a workflow for creating, screening and tuning dynamic gene circuits that operate at the population level. Specifically, we create a mutant library of an existing oscillator, the synchronized lysis circuit, and tune its dynamics while uncovering principles regarding its behavior. Lastly, we utilize this directed evolution workflow to construct a new synchronized genetic oscillator that exhibits robust dynamics over long time scales. Primary brain cell infection by Toxoplasma gondii reveals the extent and dynamics of parasite differentiation and its impact on neuron biology. Toxoplasma gondii is a eukaryotic parasite that forms latent cysts in the brain of immunocompetent individuals. The latent parasite infection of the immune-privileged central nervous system is linked to most complications. With no drug currently available to eliminate the latent cysts in the brain of infected hosts, the consequences of neurons long-term infection are unknown. It has long been known that T. gondii specifically differentiates into a latent form (bradyzoite) in neurons, but how the infected neuron responds to the infection remains to be elucidated. We have established a new in vitro model resulting in the production of mature bradyzoite cysts in brain cells. Using dual, host and parasite RNA-seq, we characterized the dynamics of differentiation of the parasite, revealing the involvement of key pathways in this process. Moreover, we identified how the infected brain cells responded to the parasite infection revealing the drastic changes that take place. We showed that neuronal-specific pathways are strongly affected, with synapse signaling being particularly affected, especially glutamatergic synapse signaling. The establishment of this new in vitro model allows investigating both the dynamics of parasite differentiation and the specific response of neurons to the long-term infection by this parasite. OncoLoop: A network-based precision cancer medicine framework At present, prioritizing cancer treatments at the individual patient level remains challenging, and performing co-clinical studies using patient-derived models in real-time is often not feasible. To circumvent these challenges, we introduce OncoLoop, a precision medicine framework to predict and validate drug sensitivity in human tumors and their pre-existing high-fidelity (cognate) model(s) by leveraging perturbational profiles of clinically-relevant oncology drugs. As proof-of-concept, we applied OncoLoop to prostate cancer (PCa) using a series of genetically-engineered mouse models (GEMMs) that recapitulate a broad spectrum of disease states, including castration-resistant, metastatic, and neuroendocrine prostate cancer. Interrogation of published cohorts using Master Regulator (MR) conservation analysis revealed that most patients were represented by at least one cognate GEMM-derived tumor (GEMM-DT). Drugs predicted to invert MR activity in patients and their cognate GEMM-DTs were successfully validated, including in two cognate allografts and one cognate patient-derived xenograft (PDX). OncoLoop is a highly generalizable framework that can be extended to other cancers and potentially other diseases.

Significance StatementOncoLoop is a transcriptomic-based experimental and computational framework that can support rapid-turnaround co-clinical studies to identify and validate drugs for individual patients, which can then be readily adapted to clinical practice. This framework should be applicable in many cancer contexts for which appropriate models and drug perturbation data are available. Modeling miRNA-driven post-transcriptional regulatory signals by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. 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. Spatio-temporal, optogenetic control of gene expression in organoids Organoids derived from stem cells become increasingly important to study human development and to model disease. However, methods are needed to control and study spatio-temporal patterns of gene expression in organoids. To this aim, we combined optogenetics and gene perturbation technologies to activate or knock-down RNA of target genes, at single-cell resolution and in programmable spatio-temporal patterns. To illustrate the usefulness of our approach, we locally activated Sonic Hedgehog (SHH) signaling in an organoid model for human neurodevelopment. High-resolution spatial transcriptomic and single-cell analyses showed that this local induction was sufficient to generate stereotypically patterned organoids in three dimensions and revealed new insights into SHHs contribution to gene regulation in neurodevelopment.

With this study, we propose optogenetic perturbations in combination with spatial transcriptomics as a powerful technology to reprogram and study cell fates and tissue patterning in organoids. Confronting population models with experimental microcosm data: from trajectory matching to state-space models O_LIPopulation and community ecology traditionally has a very strong theoretical foundation with well-known models, such as the logistic and its many variations, and many modification of the classical Lotka-Volterra predator-prey and interspecific competition models. More and more, these classical models are confronted to data via fitting to empirical time-series for purposes of projections or for estimating model parameters of interest. However, the interface between mathematical population or community models and data, provided by a statistical model, is far from trivial.
C_LIO_LIIn order to help empiricists, especially researchers working with experimental laboratory populations in micro- and mesocosms, make informed decisions, we here ask which error structure one should use when fitting classical deterministic ODE models to empirical data, from single species to community dynamics and trophic interactions. We use both realistically simulated data and empirical data from microcosms to answer this question in a Bayesian framework.
C_LIO_LIWe find that many model parameters can be estimated precisely with an appropriate choice of error structure using pure observation error or state-space models, if observation errors are not too high. However, Allee effect models are typically hard to identify and state-space models should be preferred with when model complexity increases.
C_LIO_LIOur work shows that, at least in the context of experimental laboratory populations, deterministic models can be used to describe stochastic population dynamics that include process variability and observation error. Also, many scenarios do not require a complex state-space model formulation and simpler trajectory matching is sufficient for accurate parameter estimates. Finally, we provide a comprehensive tutorial for fitting these models in R.
C_LI 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. Reconsider phase reconstruction in signals with dynamic periodicity from the modern signal processing perspective Phase is the most fundamental physical quantity when we study an oscillatory time series. There are many tools aiming to estimate phase, most of them are developed based on the analytic function model. Unfortunately, this approach might not be suitable for modern signals with intrinsic nonstartionary structure, including multiple oscillatory components, each with time-varying frequency, amplitude, and non-sinusoidal oscillation, e.g., biomedical signals. Specifically, due to the lack of consensus of model and algorithm, phases estimated from signals simultaneously recorded from different sensors for the same physiological system from the same subject might be different. This fact might challenge reproducibility, communication, and scientific interpretation and thus we need a standardized approach with theoretical support over a unified model. In this paper, after summarizing existing models for phase and discussing the main challenge caused by the above-mentioned intrinsic nonstartionary structure, we introduce the adaptive non-harmonic model (ANHM), provide a definition of phase called fundamental phase, which is a vector-valued function describing the dynamics of all oscillatory components in the signal, and suggest a time-varying bandpass filter (tvBPF) scheme based on time-frequency analysis tools to estimate the fundamental phase. The proposed approach is validated with a simulated database and a real-world database with experts labels, and it is applied to two real-world databases, each of which has biomedical signals recorded from different sensors, to show how to standardize the definition of phase in the real-world experimental environment. Specifically, we report that the phase describing a physiological system, if properly modeled and extracted, is immune to the selected sensor for that system, while other approaches might fail. In conclusion, the proposed approach resolves the above-mentioned scientific challenge. We expect its scientific impact on a broad range of applications.