Anaerobic conditions unmask antimicrobial resistance in Vibrio cholerae The antibiotic formulary is threatened by high rates of antimicrobial resistance (AMR) among enteropathogens. Enteric bacteria often undergo anaerobic respiration within the gastrointestinal tract, yet little is known about how anaerobic conditions influence AMR. The facultative enteropathogen Vibrio cholerae was chosen as a model to address this knowledge gap. V. cholerae isolates were tested for resistance under aerobic and anaerobic conditions against clinically relevant antibiotics. Using conventional breakpoints established for aerobic conditions, the odds of classifying isolates as resistant increased over 10 times for ciprofloxacin and 100 times for azithromycin under anaerobic conditions compared to aerobic conditions; for doxycycline, nearly all isolates remained sensitive. Genome-wide association studies found associations between genetic elements and AMR phenotypes that varied by oxygen exposure and antibiotic concentrations; AMR phenotypes were more heritable - and more AMR genes were discovered - under anaerobic conditions. These findings challenge the paradigm of testing facultative enteropathogens for AMR under aerobic conditions alone.
Estimating the rate of plasmid transfer with an adapted Luria-Delbrück fluctuation analysis To increase our basic understanding of the ecology and evolution of conjugative plasmids, we need a reliable estimate of their rate of transfer between bacterial cells. However, accurate estimates of plasmid transfer have remained elusive due to biological and experimental complexity. Current methods to measure transfer rate can be confounded by many factors. A notable example involves plasmid transfer between different strains or species where the rate that one type of cell donates the plasmid is not equal to the rate at which the other cell type donates. Asymmetry in these rates has the potential to bias or constrain current transfer estimates, thereby limiting our capabilities for estimating transfer in microbial communities. Inspired by the classic fluctuation analysis of Luria and Delbruck, we develop a novel approach, the Luria-Delbruck method ( LDM), for estimating plasmid transfer rate. Our new approach embraces the stochasticity of conjugation departing from the current deterministic population dynamic methods. In addition, the LDM overcomes obstacles of traditional methods by not being affected by different growth and transfer rates for each population within the assay. Using stochastic simulations and experiments, we show that the LDM has high accuracy and precision for estimation of transfer rates compared to the most widely used methods, which can produce estimates that differ from the LDM estimate by orders of magnitude.

Significance StatementConjugative plasmids play significant roles in the ecology and evolution of microbial communities. Notably, antibiotic resistance genes are often encoded on conjugative plasmids. Thus, conjugation--the transfer of a plasmid copy from one cell to another--is a common way for antibiotic resistance to spread between important clinical pathogens. For both public health modeling and a basic understanding of microbial population biology, accurate estimates of this fundamental rate are of great consequence. We show that widely used methods can lead to biased estimates, deviating from true values by several orders of magnitude. Therefore, we developed a new approach, inspired by the classic fluctuation analysis of Luria and Delbruck, for accurately assessing the rate of plasmid conjugation under a variety of conditions.
Multiple origins of nematode-Wolbachia symbiosis in supergroup F and convergent loss of bacterioferritin in filarial Wolbachia The intracellular endosymbiotic proteobacteria Wolbachia have evolved across the phyla nematoda and arthropoda. In Wolbachia phylogeny, supergroup F is the only clade with members from both arthropod and filarial nematode hosts and therefore can provide unique insights into their evolution and biology. In this study, 4 new supergroup F Wolbachia genomes have been assembled using a metagenomic assembly and binning approach, wMoz and wMpe from the human filarial parasites Mansonella ozzardi and Mansonella perstans, and wOcae and wMoviF from the blue mason bee Osmia caerulescens and the sheep ked Melophagus ovinus respectively. A comprehensive phylogenomic analysis revealed two independent origins of filarial Wolbachia in supergroup F, likely from ancestral arthropod hosts, as well as a convergent pseudogenization and loss of the bacterioferritin gene across all filarial Wolbachia. The new genomes provide a valuable resource for further studies on symbiosis, evolution, and the discovery of new antibiotics to treat mansonellosis.

Data depositionRaw Illumina reads for the M. perstans isolate Mpe1 and the M. ozzardi isolate Moz2 have been deposited in the NCBI SRA database under BioProject accession numbers PRJNA666672 and PRJNA666671 respectively. The assembled genome and PGAP annotations for the Wolbachia wMpe and wMoz are available from the NCBI GenBank database with accession numbers JACZHU000000000 and JADAKK000000000 respectively. The annotated genomes for wOcae and wMoviF are available as GenBank accessions JAHDJT000000000 and JAHRBA000000000 respectively, deposited under BioProject accessions PRJNA244005 and PRJNA739984 respectively. Classification of daily crop phenology in PhenoCams using deep learning and hidden markov models Near surface cameras, such as those in the PhenoCam network, are a common source of ground truth data in modelling and remote sensing studies. Despite having locations across numerous agricultural sites, few studies have used near surface cameras to track the unique phenology of croplands. Due to management activities, crops do not have a natural vegetation cycle which many phenological extraction methods are based on. For example, a field may experience abrupt changes due to harvesting and tillage throughout the year. A single camera can also record several different plants due to crop rotations, fallow fields, and cover crops. Current methods to estimate phenology metrics from image time series compress all image information into a relative greenness metric, which discards a large amount of contextual information. This can include the type of crop present, whether snow or water is present on the field, the crop phenology, or whether a field lacking green plants consists of bare soil, fully senesced plants, or plant residue. Here we developed a modelling workflow to create a daily time series of crop type and phenology, while also accounting for other factors such as obstructed images and snow covered fields. We used a mainstream deep learning image classification model, VGG16. Deep learning classification models do not have a temporal component, so to account for temporal correlation among images our workflow incorporates a hidden markov model in the post-processing. The initial image classification model had out of sample F1 scores of 0.83-0.85, which improved to 0.86-0.91 after all post-processing steps. The resulting time series show the progression of crops from emergence to harvest, and can serve as a daily, local scale dataset of field states and phenological stages for agricultural research. Age-associated B cells are long-lasting effectors that restrain reactivation of latent γHV68 Age-associated B cells (ABCs; CD19+CD11c+T-bet+) are a unique population that are increased in an array of viral infections, though their role during latent infection is largely unexplored. Here, we use murine gammaherpesvirus 68 ({gamma}HV68) to demonstrate that ABCs remain elevated long-term during latent infection and express IFN{gamma} and TNF. Using a strain of {gamma}HV68 that is cleared following acute infection, we show that ABCs persist in the absence of latent virus, though their expression of IFN{gamma} and TNF is decreased. With a fluorescent virus we demonstrate that ABCs are infected with {gamma}HV68 at similar rates to other previously activated B cells. We find that mice without ABCs display defects in anti-viral IgG2a/c antibodies and are more susceptible to {gamma}HV68 reactivation when challenged with heterologous infection. Together, these results indicate that ABCs are a persistent effector subset during latent viral infection that restrains {gamma}HV68 reactivation. Preprints in motion: tracking changes between preprint posting and journal publication during a pandemic Amidst the COVID-19 pandemic, preprints in the biomedical sciences are being posted and accessed at unprecedented rates, drawing widespread attention from the general public, press and policymakers for the first time. This phenomenon has sharpened longstanding questions about the reliability of information shared prior to journal peer review. Does the information shared in preprints typically withstand the scrutiny of peer review, or are conclusions likely to change in the version of record? We assessed preprints from bioRxiv and medRxiv that had been posted and subsequently published in a journal through 30th April 2020, representing the initial phase of the pandemic response. We utilised a combination of automatic and manual annotations to quantify how an article changed between the preprinted and published version. We found that the total number of figure panels and tables changed little between preprint and published articles. Moreover, the conclusions of 7.2% of non-COVID-19-related and 17.2% of COVID-19-related abstracts undergo a discrete change by the time of publication, but the majority of these changes do not qualitatively change the conclusions of the paper. 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. Second Heart Field-derived Cells Contribute to Angiotensin II-mediated Ascending Aortopathies BackgroundThe ascending aorta is a common location for aneurysm and dissection. This aortic region is populated by a mosaic of medial and adventitial cells that are embryonically derived from either the second heart field (SHF) or the cardiac neural crest. SHF-derived cells populate areas that coincide with the spatial specificity of thoracic aortopathies. The purpose of this study was to determine whether and how SHF-derived cells contribute to ascending aortopathies.

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

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

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



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


What are the clinical implications?- Heterogeneity of the embryonic origins of SMCs and fibroblasts contributes to complex mechanisms of vasculopathy formation, which should be considered when investigating the pathogenesis of thoracic aortopathies. Population differences in the length and early-life dynamics of telomeres among European pied flycatchers Telomere length and shortening rate are increasingly used as biomarkers for long-term costs in ecological and evolutionary studies because of their relationships with survival and fitness. Both early-life conditions and growth, and later-life stressors can create variation in telomere shortening rate. Studies on between- population telomere length and dynamics are scarce, despite the expectation that populations exposed to varying environmental constraints would present divergent telomere length patterns. The pied flycatcher (Ficedula hypoleuca) is a passerine bird breeding across Eurasia (from Spain to western Siberia) and migrating through the Iberian Peninsula to spend the non-breeding period in sub-Saharan Africa. Thus, different populations show marked differences in migration distance. We studied the large-scale variation of telomere length and early-life dynamics in the pied flycatcher by comparing six European populations across a north- south gradient (Finland, Estonia, England, and Spain) predicting negative effect of migration distance on adult telomere length, and of nestling growth on nestling telomere dynamics. There were clear population differences in telomere length, with English birds from mid-latitudes having the longest telomeres. Telomere length did not thus show consistent latitudinal variation and was not linearly linked to differences in migration distance. Early-life telomere shortening rate tended to vary between populations. Fast growth was associated with shorter telomeres in the early life, but faster nestling growth affected telomeres more negatively in northern than southern populations. While the sources of between-population differences in telomere-related biology remain to be more intensively studied, our study illustrates the need to expand telomere studies at the between-population level. Intermittent compression increases flow through a gel: a mechanical model perhaps relevant to flow of cerebrospinal fluid Exchange of molecules between cerebrospinal fluid (CSF) and brain cells contributes to brain function and protection from dementia. Despite widespread acceptance of the  glymphatic theory, the route by which CSF is brought close enough to the neural tissue for solutes to be exchanged by extracellular diffusion is not entirely clear. Exogenous molecules injected into CSF are observed to reach the basement lamina that surrounds the dense capillary network. Transport of solutes by diffusion along the basement lamina, a gel of macromolcules about 100 nm thick, would be too slow; bulk flow in a static geometry would require unphysiologically high pressures. However, it is known that the pulsation of blood aids transport of CSF, and we hypothesized that this is because the pulsation intermittently squeezes the pericapillary lamina. In a primitive mimicry, we have tested whether intermittent squeezing increases flow through an agar gel. In all but one of 216 tests, pulsation caused a reversible increase, sometimes by a factor of 100 or more. The enhancement was greatest for frequencies 5-11 Hz and, over the tested range of pressure heads (20 - 50 cmH2O), was greatest for the lowest pressure. The results suggest a reason why exercise slows the aging of the brain. Universality of evolutionary trajectories under arbitrary competition dynamics The assumption of constant population size is central in population genetics. It led to a large body of results, that are robust to modeling choices and that have proven successful to understand evolutionary dynamics. In reality, allele frequencies and population size are both determined by the interaction between a population and the environment. Relaxing the constant-population assumption have two big drawbacks. It increases the technical difficulty of the analysis, and it requires specifying a mechanism for the saturation of the population size, possibly making the results contingent on model details. Here, we develop a framework that encompasses a great variety of systems with an arbitrary mechanism for population growth limitation. By using techniques based on scale separation for stochastic processes, we are able to calculate analytically properties of evolutionary trajectories, such as the fixation probability. Remarkably, these properties assume a universal form with respect to our framework, which depends on only three parameters related to the inter-generation timescale, the invasion fitness, and the carrying capacity of the strains. In other words, different systems, such as Lotka-Volterra or a chemostat model (contained in our framework), share the same evolutionary outcomes after a proper re-mapping of their parameters. An important and surprising consequence of our results is that the direction of selection can be inverted, with a population evolving to reach lower values of invasion fitness. Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable and highly reproducible in terms of microbial composition, function and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation. 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. Dynamic modulation of enhancer responsiveness by core promoter elements in living Drosophila embryo Regulatory interactions between enhancers and core promoters are fundamental for the temporal and spatial specificity of gene expression in development. The central role of core promoters is to initiate productive transcription in response to enhancers activation cues. However, it has not been systematically assessed how individual core promoter elements affect the induction of transcriptional bursting by enhancers. Here, we provide evidence that each core promoter element dynamically alters "enhancer responsiveness" by changing functional parameters of transcriptional bursting in developing Drosophila embryos. Quantitative live imaging analysis revealed that the timing and the continuity of burst induction are common regulatory steps on which core promoter elements impact. We also show that core promoters use the upstream TATA box as an amplifier of transcriptional bursting upon interaction with enhancers. Genome editing analysis of the pair-rule gene fushi tarazu further revealed that the endogenous TATA box and the downstream core promoter element are both essential for the correct expression and the function of key developmental genes in early embryos. We suggest that core promoter elements serve as a key regulatory module in converting enhancer activity into transcription dynamics during animal development. Cast away on Mindoro island: could density dependence threaten tamaraw (Bubalus mindorensis) on the long term? Endangered species, despite often living at low population densities, may undergo unexpected density-dependent feedbacks in the case of successful recovery or marked reduction in range. Because density-dependence dynamics can increase risk of extinction, these effects can hamper conservation efforts. In this study, we analyse the dynamics of the largest population of the tamaraw (Bubalus mindorensis), a critically endangered species endemic to the Mindoro island, Philippines. The population is located within a < 3,000 ha area in Mounts Iglit-Baco Natural Park, with limited expansion possibilities. We took advantage of a 21 year time series of tamaraw counts to estimate annual population growth rate and possible density-dependence, accounting for sampling errors in the counts. The tamaraw population has been increasing in size at an average rate of +5% per year, as would be expected given its protected status by law. Population growth showed strong spatial structuring, with a population growth close to +10% in the core area of protection, and a reduction of abundance of - 5% at the periphery of its range, inside the protected area. This range constriction is concerning because our best population dynamic model suggests significant negative density-dependence (Bayes factor = 0.9). The contraction of tamaraw range is likely caused by anthropogenic pressures forcing herbivores to live at relatively high densities in the core zone where protection is most effective. Our study highlights the fact that, despite the continuous population growth over the last two decades, the long-term viability of the Mounts Iglit-Baco Natural Park tamaraw population remains uncertain. Bayesian Additive Regression Trees for Genotype by Environment Interaction Models We propose a new class of models for the estimation of genotype by environment (GxE) interactions in plant-based genetics. Our approach, named AMBARTI, uses semi-parametric Bayesian additive regression trees to accurately capture marginal genotypic and environment effects along with their interaction in a fully Bayesian model. We demonstrate that our approach is competitive or superior to similar models widely used in the literature via both simulation and a real world data set. Furthermore, we introduce new types of visualisation to properly assess both the marginal and interactive predictions from the model. An R package that implements our approach is available at https://github.com/ebprado/ambarti. A rat model of pregnancy in the male parabiont Male pregnancy is a unique phenomenon in syngnathidae which refers to embryo or fetus incubation by males. However, whether male mammalian animals have the potential to conceive and maintain pregnancy remains unclear. Here, we constructed a rat model of pregnancy in the male parabiont using a four-step strategy. First, a heterosexual parabiotic pair was produced by surgically joining a castrated male rat and a female rat. After 8 weeks, uterus transplantation was performed on the male parabiont. After recovery, blastocyst-stage embryos were transplanted to the grafted uterus of the male parabiont and the native uterus of the female parabiont. Caesarean section was performed at embryonic day 21.5. The success rate of the model was 3.68%, and pups were only delivered from the male parabiont under pregnant blood exposure. Our results reveal the potential for rat embryonic development in male parabionts, which may have a profound impact on reproductive biology.

One-sentence summaryA rat model of pregnancy was constructed in the male parabiont using a four-step strategy.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=152 SRC="FIGDIR/small/447686v2_ufig1.gif" ALT="Figure 1">
View larger version (46K):
org.highwire.dtl.DTLVardef@1f0694eorg.highwire.dtl.DTLVardef@91997forg.highwire.dtl.DTLVardef@f47086org.highwire.dtl.DTLVardef@7a33f_HPS_FORMAT_FIGEXP  M_FIG C_FIG Distinct cellular origins and differentiation process account for distinct oncogenic and clinical behaviors of leiomyosarcomas Differentiated phenotypes in solid tumors mostly lead to a better outcome than in their more immature counterparts. In leiomyosarcoma (LMS), a very aggressive tumor for which neither chemotherapy nor targeted or immune therapies have demonstrated efficacy, a relatively transcriptionally uniform subgroup of well differentiated tumors has been described and is associated with poor survival. In this disease, differentiation does not systematically hinder tumor progression, so it is a good model for studying the relation between these two mechanisms. By applying ad hoc bioinformatic methods, we conducted a comprehensive analysis in order to detect subgroups of leiomyosarcoma and finely investigate their transcriptomic and genomic profiles. We first used unsupervised hierarchical clustering to identify the most transcriptionally homogeneous LMS subgroup, which we named  hLMS. We derived a transcriptional signature that could consistently stratify patients with similar clinical and genomic characteristics in TCGA and ICGC cohorts. Functional analysis of the signature showed over-expression of SRF and E2Fs targets harboring histone modifications that were specifically identified in smooth muscles in hLMS. By integrating multi-omics data, we identified punctual and structural variations as well as the potential cellular origins of each group of LMS using GTEX normal tissue gene expression. Altogether, our results suggest that hLMS, but not oLMS (other LMS), may originate from vascular smooth muscle cells in which oncogenesis is driven by the MYOCD/SRF axis, which is essential for hLMS survival. We propose that the phenotypic plasticity of vascular smooth muscle cells coupled with MYOCD/SRF pathway overactivation could explain how hLMS balance this uncommon interplay between differentiation and aggressiveness. F1F0 ATP Hydrolysis is a Determinant of Metabolic Rate, a Correlate of Lifespan, and a Weakness of Cancer Dramatic rectal temperature decrease, in mice administered with a drug that specifically inhibits F1F0 ATP hydrolysis, suggests that F1F0 ATP hydrolysis is a major determinant of metabolic rate in vivo. Across twelve investigated species, less F1F0 ATP hydrolysis correlates with greater maximal lifespan. Specific drug inhibition of F1F0 ATP hydrolysis exerts potent anti-cancer activity in vitro. Natural Selection as the Sum Over All Histories. As evolution can be connected to the principle of least action, and if it is depicted in evolution-space versus time then it corresponds to the direction of ultimate causation. As an organism evolves and follows a path of proximate causation, if the vector is closely parallel to that of the Ultimate Causation then the changes will confer desirable attributes which will lead to further development. If, however, the variations do not occur in a direction close to that of the ultimate causation vector the evolved organism will quickly die out. Therefore Natural Selection may be viewed as similar to Feynmans "sum over all histories". This approach is compatible with both Neutral Theory and Selection, as it includes both positive and negative mutations and selection. Therefore, the principle of least action gives a direction, but not a purpose, to evolution. Exploring the Interactome of Cytochrome P450 2E1 in Human Liver Microsomes with Chemical Crosslinking Mass Spectrometry Aiming to elucidate the grounds of system-wide effects of the alcohol-induced increase in the content of cytochrome P450 2E1 (CYP2E1) on drug metabolism, we explored the array of its protein-protein interactions (proteome) in human liver microsomes (HLM) with chemical crosslinking mass spectrometry (CXMS). Our strategy employs membrane incorporation of purified CYP2E1 modified with photoreactive crosslinkers benzophenone-4-maleimide (BPM) and 4-(N-succinimidylcarboxy)benzophenone (BPS). Exposure of bait-incorporated HLM samples to light was followed by isolating the His-tagged bait protein and its crosslinked aggregates on Ni-NTA agarose. Analyzing the individual bands of SDS-PAGE slabs of thereby isolated protein with the toolset of untargeted proteomics, we detected the crosslinked dimeric and trimeric complexes of CYP2E1 with other drug-metabolizing enzymes. Among the most extensively crosslinked partners of CYP2E1 are the cytochromes P450 2A6, 2C8, 3A4, 4A11, and 4F2. Of particular interest are the interactions of CYP2E1 with the latter two enzymes, which are involved in the synthesis and disposal of vasoactive and pro-inflammatory eicosanoids. We also detected the conjugates of CYP2E1 with UDP-glucuronosyltransferases (UGTs) 1A and 2B, fatty aldehyde dehydrogenase (ALDH3A2), epoxide hydrolase 1 (EPHX1), disulfide oxidase 1 (ERO1L), and ribophorin II (RPN2). These results demonstrate the exploratory power of the proposed CXMS strategy and corroborate the concept of tight functional integration in the human drug-metabolizing ensemble through protein-protein interactions of the constituting enzymes. Shedding Light on the Circadian Clock of the Threespine Stickleback The circadian clock is an internal timekeeping system shared by most organisms, and knowledge about its functional importance and evolution in natural environments is still needed. Here, we investigated the circadian clock of wild-caught threespine sticklebacks (Gasterosteus aculeatus) at the behavioural and molecular levels. While their behaviour, ecology, and evolution are well studied, information on their circadian rhythms are scarce. We quantified the daily locomotor activity rhythm under a light-dark cycle (LD) and under constant darkness (DD). Under LD, all fish exhibited significant daily rhythmicity, while under DD, only 18% of individuals remained rhythmic. This interindividual variation suggests that the circadian clock controls activity only in certain individuals. Moreover, under LD, some fish were almost exclusively nocturnal, while others were active around the clock. Furthermore, the most nocturnal fish were also the least active. These results suggest that light masks activity (i.e. suppresses activity without entraining the internal clock) more strongly in some individuals than others. Finally, we quantified the expression of five clock genes in the brain of sticklebacks under DD using qPCR. We did not detect circadian rhythmicity, which could either indicate that the clock molecular oscillator is highly light-dependent, or that there was an oscillation but that we were unable to detect it. Overall, our study suggests that a strong circadian control on behavioural rhythms may not necessarily be advantageous in a natural population of sticklebacks and that the daily phase of activity varies greatly between individuals because of a differential masking effect of light.

Summary statementWe found that in wild-caught threespine sticklebacks, the circadian clock does not control locomotor activity in most, but not all, individuals. Sticklebacks are mostly nocturnal, although interindividual variation exists. OrganoID: a versatile deep learning platform for tracking and analysis of single-organoid dynamics Organoids have immense potential as ex vivo disease models for drug discovery and personalized drug screening. Dynamic changes in individual organoid morphology, number, and size can indicate important drug responses, however these metrics are difficult and labor-intensive to obtain for high-throughput image datasets. Here, we present OrganoID, a robust image analysis platform that automatically recognizes, labels, and tracks single organoids, pixel-by-pixel, in brightfield and phase-contrast microscopy experiments. The platform was trained on images of pancreatic cancer organoids and validated on separate images of pancreatic, lung, colon, and adenoid cystic carcinoma organoids, which showed excellent agreement with manual measurements of organoid count (96%) and size (95%) without any parameter adjustments. Single-organoid tracking accuracy remained above 89% over a four-day time-lapse microscopy study. Automated single-organoid morphology analysis of a chemotherapy dose-response experiment identified decreased organoid circularity as an important morphological feature reflecting drug response. OrganoID enables straightforward, detailed, and accurate image analysis to accelerate the use of organoids in high-throughput, data-intensive biomedical applications. 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. Endometriotic Organoids: A Novel In Vitro Model of Endometriotic Lesion Development. The development and progression of endometriotic lesions are poorly understood, but immune cell dysfunction and inflammation are closely associated with the pathophysiology of endometriosis. A lack of suitable 3D in vitro models permitting the study of interactions between cell types and the microenvironment is a contributing factor. To address this limitation, we developed endometriotic organoids (EO) to explore the role of epithelial-stromal interactions and model peritoneal cell invasion associated with lesion development. Using a non-adherent microwell culture system, spherical organoids were generated with endometriotic epithelial cells (12Z) combined with immortalized endometriotic stromal cells (iEc-ESC) or immortalized uterine stromal cells (iHUF). Organoids self-organized with stromal cells occupying the center and epithelial cells on the periphery of the organoid. Endometriotic organoids (EO), containing iEc-ESC, resulted in the development of stratified 12Z epithelial cells compared to those with iHUF where the 12Z cells developed as a single layered epithelium. Transcriptomic analysis found 4,522 differentially expressed genes (DEG) between EO and 12Z/iHUF organoids, and the top DEG included increased expression of interleukins and prostaglandin synthase enzymes. An overlap of the EO DEG with baboon endometriotic lesions was highly significant. Finally, to mimic invasion of endometrial tissue into the peritoneum, a model was developed using EO and extracellular matrix containing human peritoneal mesothelial cells (LP9). Invasion of EO into the extracellular matrix-LP9 layer was increased in presence of estrogen or THP1-derived proinflammatory macrophages. Taken together, our results strongly support the concept that EO are an appropriate model for dissecting mechanisms that contribute to endometriotic lesion development.

One Sentence SummaryEndometriotic organoids are an appropriate model to study epithelial-stromal interactions and model cell invasion associated with lesion development. Dynamic trajectories of connectome state transitions are heritable The brains functional connectome is dynamic, constantly reconfiguring in an individual-specific manner. However, which characteristics of such reconfigurations are subject to genetic effects, and to what extent, is largely unknown. Here, we identified heritable dynamic features, quantified their heritability, and determined their association with cognitive phenotypes. In resting-state fMRI, we obtained multivariate features, each describing a temporal or spatial characteristic of connectome dynamics jointly over a set of connectome states. We found strong evidence for heritability of temporal features, particularly fractional occupancy (FO) and transition probability (TP), describing the trajectory of state transitions. Genetic effects explained a substantial proportion of phenotypic variance of these features (h2=.39, 95% CI= [.24.54] for FO; h2=.43, 95% CI=[.29.57] for TP). Moreover, these temporal phenotypes were associated with cognitive performance. Contrarily, we found no robust evidence for heritability of spatial features of the dynamic states (i.e., states Modularity and connectivity pattern). Genetic effects may therefore primarily contribute to how the connectome transitions across states, rather than the precise spatial instantiation of the states in individuals. In sum, genetic effects impact the duration spent in each connectivity configuration and the frequency of shifting between configurations, and such temporal features may act as endophenotypes for cognitive abilities.