Sexual selection for males with beneficial mutations Sexual selection is the process by which traits providing a mating advantage are favoured. Theoretical treatments of the evolution of sex by sexual selection propose that it operates by reducing the load of deleterious mutations. Here, we postulate instead that sexual selection primarily acts through females preferentially mating with males carrying beneficial mutations. We used simulation and analytical modelling to investigate the evolutionary dynamics of beneficial mutations in the presence of sexual selection. We found that female choice for males with beneficial mutations had a much greater impact on genetic quality than choice for males with low mutational load. We also relaxed the typical assumption of a fixed mutation rate. For deleterious mutations, mutation rate should always be minimized, but when rare beneficial mutations can occur, female choice for males with those rare beneficial mutations could overcome a decline in average fitness and allow an increase in mutation rate. We propose that sexual selection for beneficial mutations could overcome the  two-fold cost of sex much more readily than choice for males with low mutational load and may therefore be a more powerful explanation for the prevalence of sexual reproduction than the existing theory. If sexual selection results in higher fitness at higher mutation rates and if the variability produced by mutation itself promotes sexual selection, then a feedback loop between these two factors could have had a decisive role in driving adaptation.
Hybridization alters the shape of the genotypic fitness landscape, increasing access to novel fitness peaks during adaptive radiation Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 SNPs that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.
Carbon Dot-based Fluorescent Antibody Nanoprobes as Brain Tumour Glioblastoma Diagnostics. The development of efficient and sensitive tools for the detection of brain cancer in patients is of the utmost importance particularly because many of these tumours go undiagnosed until the disease has advanced and when treatment is less effective. Current strategies employ antibodies (Abs) to detect Glial Fibrillary Acid Protein (GFAP) in tissue samples, since GFAP is unique to the brain and not present in normal peripheral blood, and it relies on fluorescent reporters.

Herein we describe a low cost, practical and general method for the labelling of proteins and antibodies with fluorescent carbon dots (CD) to generate diagnostic probes that are robust, photostable and applicable to the clinical setting. The two-step protocol relies on the conjugation of a dibenzocyclooctyne (DBCO)-functionalised CD with azide functionalised proteins by combining amide conjugation and strain promoted alkyne-azide cycloaddition (SPAAC) ligation chemistry. The new class of Ab-CD conjugates developed using this strategy was successfully used for the immunohistochemical staining of human brain tissues of patients with glioblastoma (GBM) validating the approach. Overall, these novel fluorescent probes offer a promising and versatile strategy in terms of costs, photostability and applicability which can be extended to other Abs and protein systems. Dilute phase oligomerization can oppose phase separation and modulate material properties of a ribonucleoprotein condensate Ribonucleoprotein bodies are exemplars of membraneless biomolecular condensates that can form via spontaneous or driven phase transitions. The fungal protein Whi3 forms ribonucleoprotein condensates with different RNA molecules, and these condensates are implicated in key processes such as cell-cycle control and generating cell polarity. Whi3 has a modular architecture that includes a Q-rich intrinsically disordered region (IDR) and a tandem RNA recognition module. Here, we demonstrate that a 21-residue stretch within the Q-rich IDR has a weak intrinsic preference for forming alpha-helical conformations. Through mutagenesis, we find that increased alpha helicity enhances oligomerization in the dilute phase. One consequence of enhanced oligomerization is a dilution of Whi3 in the dense phase. The opposite behavior is observed when helicity within the 21-residue stretch of the Q-rich region is abrogated. Thus, the formation of dilute phase oligomers, driven by a specific sequence motif and potential synergies with the rest of the IDR, opposes incorporation of the Whi3 protein into the dense phase, thereby altering the dense phase stoichiometry of protein to RNA. Our findings, which stand in contrast to other systems where oligomerization has been shown to enhance the drive for phase separation, point to a novel mechanism that might be operative for influencing compositions of condensates. Our work also points to routes for designing synthetic ribonucleoprotein condensates whereby modulation of protein oligomerization via homotypic interactions can impact dense phase concentrations, stoichiometries, and material properties.

SignificanceA large sub-class of biomolecular condensates are linked to RNA regulation and are known as ribonucleoprotein (RNP) bodies. While extensive work has identified driving forces for biomolecular condensate formation, relatively little is known about forces that oppose assembly. Here, using a fungal RNP protein, Whi3, we show that a portion of its intrinsically disordered, glutamine-rich region modulates phase separation by forming transient alpha helical structures that promote the assembly of dilute phase oligomers. These oligomers detour Whi3 proteins from condensates, thereby impacting the driving forces for phase separation, the protein-to-RNA ratio in condensates, and the material properties of condensates. Our findings show how nanoscale conformational and oligomerization equilibria can influence mesoscale phase equilibria. Large protein complex interfaces have evolved to promote cotranslational assembly Assembly pathways of protein complexes should be precise and efficient to minimise misfolding and unwanted interactions with other proteins in the cell. One way to achieve this efficiency is by seeding assembly pathways during translation via the cotranslational assembly of subunits. While recent evidence suggests that such cotranslational assembly is widespread, little is known about the properties of protein complexes associated with the phenomenon. Here, using a combination of proteome-specific protein complex structures and publicly available ribosome profiling data, we show that cotranslational assembly is particularly common between subunits that form large intermolecular interfaces. To test whether large interfaces have evolved to promote cotranslational assembly, as opposed to cotranslational assembly being a non-adaptive consequence of large interfaces, we compared the sizes of first and last translated interfaces of heteromers in bacterial, yeast, and human complexes. We detect a clear preference in N-terminal interfaces, which will be translated first and thus more likely to form cotranslationally, to be larger than C-terminal interfaces, suggesting that large interfaces have evolved as a means to maximise the chance of successful cotranslational subunit binding. A novel camera trap design for studying wildlife in mountain glacier ecosystems yields new insight for glacier biodiversity in the Pacific Northwest, USA ContextThe global recession of glaciers and perennial snowfields is reshaping mountain ecosystems. Beyond physical changes to the landscape and altered downstream hydrology, the implications of glacier decline for biodiversity are poorly known. Before predictions can be made about how climate change will affect wildlife in glacier-associated ecosystems, a more thorough accounting of the role that glaciers play in species life histories is needed. However, typical approaches for documenting wildlife presence and behavior--remote camera traps--are difficult to use in glaciated terrain due to limited options for securing them (e.g., no trees) and dramatic seasonal changes in snowpack.

AimsIn this study, we sought to test a novel camera trap designed for glaciated mountain ecosystems. We also aimed to use this approach to gain insight into wildlife and human usage of a mountain glacier in western North America.

MethodsWe deployed an elevational transect of uniquely designed camera traps along the western margin of the Paradise Glacier, a rapidly receding mountain glacier on the south side of Mount Rainier, WA, USA. Our simple camera trap design consisted of a wildlife camera attached to a camouflaged cylindrical cooler filled with snow and rocks.

Key resultsOur camera design proved ideal for a mountain glacier ecosystem and from June to September 2021, we detected at least 16 vertebrate species (seven birds, nine mammals) over 770 trap nights using glacier-associated habitats. Humans, primarily skiers, were the most common species detected, but we also recorded 99 observations of wildlife (birds and mammals). These included three species of conservation concern in Washington: wolverine (Gulo gulo), Cascade red fox (Vulpes vulpes cascadensis), and white-tailed ptarmigan (Lagopus leucura).

ConclusionsCollectively, our results provide proof-of-concept for a novel camera trap design that is ideal for treeless, perennially snow-covered landscapes and revealed a rich diversity of wildlife using mountain glacier habitat in the Pacific Northwest. We highlight the global need for similar studies to better understand the true scale of biodiversity that will be impacted by glacier recession in mountain ecosystems. Cortical network structure mediates response to stimulation: An optogenetic study in non-human primate cortex As aberrant network-level functional connectivity underlies a variety of neural disorders, the ability to induce targeted functional reorganization would be a profound development towards therapies for neural disorders. Brain stimulation has been shown to alter large-scale network-wide functional connectivity, but the mapping from stimulation to the modification is unclear. Here, we leverage advances in neural interfaces, interpretable machine learning, and graph theory to arrive at a model which accurately predicts stimulation-induced network-wide functional reorganization. The model jointly considers the stimulation protocol and the cortical network structure, departing from the standard approach which only considers the stimulation protocol. We validate our approach in the primary sensorimotor cortex of non-human primates using paired optogenetic stimulation through a large-scale optogenetic interface. We observe that the stimulation protocol only predicts a small portion of the induced functional connectivity changes while the network structure predicts much more, indicating that the underlying network is the primary mediator of the response to stimulation. We extract the relationships linking the stimulation and network characteristics to the functional connectivity changes and observe that the mappings diverge over frequency bands and successive stimulations. Finally, we uncover shared processes governing real-time and longer-term effects of stimulation. Our framework represents a paradigm shift for targeted neural stimulation and can be used to interrogate, improve, and develop stimulation-based interventions for neural disorders.

TeaserBrain stimulation rewires the brain, but the pre-existing network structure of the brain controls the rewiring. Intestine-enriched apolipoprotein b orthologs are required for stem cell differentiation and regeneration in planarians Lipid metabolism plays an instructive role in regulating stem cell state and differentiation. However, the roles of lipid mobilization and utilization in stem cell-driven regeneration are unclear. Planarian flatworms readily restore missing tissue due to injury-induced activation of pluripotent somatic stem cells called neoblasts. Here, we identify two intestine-enriched orthologs of apolipoprotein b, apob-1 and apob-2, which mediate transport of neutral lipid stores from the intestine to target tissues including neoblasts, and are required for tissue homeostasis and regeneration. Inhibition of apob function by RNAi causes head regression and lysis in uninjured animals, and delays body axis re-establishment and regeneration of multiple organs in amputated fragments. Furthermore, apob RNAi causes expansion of the population of differentiating neoblast progeny and dysregulates expression of genes enriched in differentiating and mature cells in eight major cell type lineages. We conclude that intestine-derived lipids serve as a source of metabolites required for neoblast differentiation. Hidden dynamics of proteasome autoregulation discovered by cryo-EM data-driven deep learning The cellular functions are executed by macromolecular complexes in nonequilibrium dynamic processes, which exhibit a vast diversity of conformational states. Solving conformational continuum of important complexes at atomic level remains a key challenge in the methodology of computational biophysics. Here we introduce a deep-learning framework, named AlphaCryo4D, which enables atomic-level cryogenic electron microscopy reconstructions of nonequilibrium conformational continuum. AlphaCryo4D integrates 3D deep residual learning with manifold embedding of energy landscapes, which simultaneously improves 3D classification accuracy and reconstruction resolution via an energy-based particle-voting algorithm. Applications of this methodology to six experimental datasets demonstrated its generality in breaking resolution barrier of visualizing dynamic components of functional complexes, in choreographing continuous inter-subunit motions and in exploring their  hidden conformational space. Importantly, AlphaCryo4D identified 64 conformers of the substrate-bound human 26S proteasome and reconstituted single-nucleotide-exchange dynamics of proteasomal AAA-ATPase motor during translocation initiation, which together illuminates a grand hierarchical allostery for proteasome autoregulation. Citation needed?Wikipedia bibliometrics during the first waveof the COVID-19 pandemic. BackgroundWith the COVID-19 pandemics outbreak, millions flocked to Wikipedia for updated information. Amid growing concerns regarding an "infodemic", ensuring the quality of information is a crucial vector of public health. Investigating if and how Wikipedia remained up to date and in line with science is key to formulating strategies to counter misinformation. Using citation analyses, we asked: which sources informed Wikipedias COVID-19-related articles before and during the pandemics first wave (January-May 2020).

ResultsWe found that coronavirus-related articles referenced trusted media sources and high-quality academic research. Moreover, despite a surge in COVID-19 preprints, Wikipedia had a clear preference for open-access studies published in respected journals and made little use of preprints. Building a timeline of English COVID-19 articles from 2001-2020 revealed a nuanced trade-off between quality and timeliness. It further showed how preexisting articles on key topics related to the virus created a framework for integrating new knowledge. Supported by a rigid sourcing policy, this "scientific infrastructure" facilitated contextualization and regulated the influx of new information. Lastly, we constructed a network of DOI-Wikipedia articles, which showed the shifting landscape of pandemic-related knowledge on Wikipedia and how academic citations create a web of shared knowledge supporting topics like COVID-19 vaccine development.

ConclusionsUnderstanding how scientific research interacts with the digital knowledge-sphere during the pandemic provides insight into how Wikipedia can facilitate access to science. It also reveals how, aided by what we term its "citizen encyclopedists", it successfully fended off COVID-19 disinformation and how this unique model may be deployed in other contexts. Intermittent Bulk Release of Human Cytomegalovirus Human Cytomegalovirus (HCMV) can infect a variety of cell types by using virions of varying glycoprotein compositions. It is still unclear how this diversity is generated, but spatio-temporally separated envelopment and egress pathways might play a role. So far, one egress pathway has been described in which HCMV particles are individually enveloped into small vesicles and are subsequently exocytosed continuously. However, some studies have also found enveloped virus particles inside multivesicular structures but could not link them to productive egress or degradation pathways.

We used a novel 3D-CLEM workflow allowing us to investigate these structures in HCMV morphogenesis and egress at high spatio-temporal resolution. We found that multiple envelopment events occurred at individual vesicles leading to multiviral bodies (MViBs), which subsequently traversed the cytoplasm to release virions as intermittent bulk pulses at the plasma membrane to form extracellular virus accumulations (EVAs). Our data support the existence of a novel bona fide HCMV egress pathway, which opens the gate to evaluate divergent egress pathways in generating virion diversity. 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. Fusobacterium nucleatum induces pancreatic cancer cell proliferation and migration by regulating host autocrine and paracrine signaling Pancreatic ductal adenocarcinoma (PDAC) harbors a complex tumor microbiome that has been implicated in cancer progression and resistance to chemotherapy. Recent clinical investigations uncovered a correlation between high loads of intratumor Fusobacterium nucleatum and decreased patient survival. Here we show that healthy and cancerous pancreatic cell lines harboring intracellular F. nucleatum secrete increased levels of cancer-associated cytokines including GM-CSF, CXCL1, IL-8, and MIP-3. We report that GM-CSF (granulocyte-macrophage colony stimulating factor) secretion directly increases the proliferation and migration of pancreatic cancer cells via an autocrine mechanism, notably in the absence of immune cell participation. Furthermore, we show that non-cancerous pancreatic epithelial cells do not exhibit increased proliferation or migration in response to these cytokines, but nevertheless, their secreted cytokines stimulate these responses in cancerous cell lines through paracrine signaling. Our results provide evidence that intratumor F. nucleatum in the pancreas elicits an infection-specific cytokine secretion profile from both normal and cancerous cells that adversely contributes to cancer progression through autocrine and paracrine mechanisms. Therefore, these results support the importance of investigating the contributions of both microbiome and host driven processes in pancreatic cancer to guide future therapeutic interventions. Microvesicles Transfer Mitochondria and Increase Mitochondrial Function in Brain Endothelial Cells We have demonstrated, for the first time that microvesicles, a sub-type of extracellular vesicles (EVs) derived from hCMEC/D3: a human brain endothelial cell (BEC) line transfer polarized mitochondria to recipient BECs in culture and to neurons in mice acute brain cortical and hippocampal slices. This mitochondrial transfer increased ATP levels by 100 to 200-fold (relative to untreated cells) in the recipient BECs exposed to oxygen-glucose deprivation, an in vitro model of cerebral ischemia. We have also demonstrated that transfer of microvesicles, the larger EV fraction, but not exosomes resulted in increased mitochondrial function in hypoxic endothelial cultures. Gene ontology and pathway enrichment analysis of EVs revealed a very high association to glycolysis-related processes. In comparison to heterotypic macrophage- derived EVs, BEC-derived EVs demonstrated a greater selectivity to transfer mitochondria and increase endothelial cell survival under ischemic conditions.

HighlightsO_LIMicrovesicles transfer mitochondria to endothelial cells and brain slice neurons
C_LIO_LIMitochondrial transfer increased ATP in ischemic brain endothelial cells (BECs)
C_LIO_LITransfer of microvesicles increased mitochondrial function in hypoxic BECs
C_LIO_LITransfer of exosomes did not affect mitochondrial function in hypoxic BECs
C_LIO_LIHomotypic BEC-derived EVs result in greater ATP levels in the recipient BECs
C_LI An NLR Integrated Domain toolkit to identify plant pathogen effector targets Plant resistance genes (or NLR "Nod-like Receptors") are known to contain atypical domains procuring them with a decoy capacity. Some of these integrated domains (or ID) allow the plant to lure the virulence determinants ("effectors") of pathogens and triggering a specific NLR immune reaction.

In this work, our goal was to generate a library of known IDs that could be screened with plant pathogen effectors in order to identify putative new effector virulence targets and NLR-effector pairs.

We curated the IDs contained in NLRs from seven model and crop plant species. We cloned 52 IDs representing 31 distinct Pfam domains. This library was screened for interaction by yeast-two-hybrid with a set of 31 conserved Ralstonia solanacearum type III effectors. This screening and the further in planta interaction assay allowed us to identify three interactions, involving different IDs (kinase, DUF3542, WRKY) and two type III effectors (RipAE and PopP2).

PopP2 was found to physically interact with ID#85, an atypical WRKY domain integrated in the GmNLR-ID85 NLR protein from Soybean. Using a imaging method in living plant cells, we showed that PopP2 associates with ID#85 in the nucleus. But unlike the known WRKY-containing Arabidopsis RRS1-R NLR receptor, this newly identified soybean WRKY domain could not be acetylated by PopP2 and its atypical sequence (WRKYGKR) also probably renders it inefficient in plant immunity triggering.

This ID toolkit is available for screening with other plant pathogen effectors and should prove useful to discover new effectors targets and potentially engineer new plant resistance genes. Mitochondrial genomes of two parasitic Cuscuta species lack clear evidence of horizontal gene transfer and retain unusually fragmented ccmFC genes BackgroundThe intimate association between parasitic plants and their hosts favours the exchange of genetic material, potentially leading to horizontal gene transfer (HGT) between plants. With the recent publication of several parasitic plant nuclear genomes, there has been considerable focus on such non-sexual exchange of genes. To enhance the picture on HGT events in a widely distributed parasitic genus, Cuscuta (dodders), we assembled and analyzed the organellar genomes of two recently sequenced species, C. australis and C. campestris, making this the first account of complete mitochondrial genomes (mitogenomes) for this genus.

ResultsThe mitogenomes are 265,696 and 275,898 bp in length and contain a typical set of mitochondrial genes, with ten missing or pseudogenized genes often lost from angiosperm mitogenomes. Each mitogenome also possesses a structurally unusual ccmFC gene, which exhibits splitting of one exon and a shift to trans-splicing of its intron. Based on phylogenetic analysis of mitochondrial genes from across angiosperms and similarity-based searches, there is little to no indication of HGT into the Cuscuta mitogenomes. A few candidate regions for plastome-to-mitogenome transfer were identified, with one suggestive of possible HGT.

ConclusionsThe lack of HGT is surprising given examples from the nuclear genomes, and may be due in part to the relatively small size of the Cuscuta mitogenomes, limiting the capacity to integrate foreign sequences. Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain. An epigenomic shift in amygdala marks the transition to maternal behaviors in alloparenting virgin female mice In many species, adults care for young offspring that are not their own, a phenomenon called alloparenting. However, most nonparental adults must be sensitized by repeated or extended exposures to newborns before they will robustly display parental-like behaviors. To capture neurogenomic events underlying the transition to active parental caring behaviors, we analyzed brain gene expression and chromatin profiles of virgin female mice co-housed with mothers during pregnancy and after birth. After an initial display of antagonistic behaviors and a surge of defense-related gene expression, we observed a dramatic shift in the chromatin landscape specifically in amygdala of the pup-exposed virgin females, accompanied by a dampening of anxiety-related gene expression. This epigenetic shift coincided with hypothalamic expression of the oxytocin gene and the emergence of behaviors and gene expression patterns classically associated with maternal care. The results outline a neurogenomic program associated with dramatic behavioral changes and suggest molecular networks relevant to human postpartum mental health. Trade-offs in biosensor optimization for dynamic pathway engineering Recent progress in synthetic biology allows the construction of dynamic control circuits for metabolic engineering. This technology promises to overcome many challenges encountered in traditional pathway engineering, thanks to their ability to self-regulate gene expression in response to bioreactor perturbations. The central components in these control circuits are metabolite biosensors that read out pathway signals and actuate enzyme expression. However, the construction of metabolite biosensors is a major bottleneck for strain design, and a key challenge is to understand the relation between biosensor dose-response curves and pathway performance. Here we employ multiobjective optimization to quantify performance trade-offs that arise in the design and calibration of metabolite biosensors. Our approach reveals strategies for tuning dose-response curves along an optimal trade-off between production flux and the cost of an increased expression burden on the host. We explore properties of control architectures built in the literature, and identify their advantages and caveats in terms of performance and robustness to growth conditions and leaky promoters. We demonstrate the optimality of a control circuit for glucaric acid production in Escherichia coli, which has been shown to increase titer by 2.5-fold as compared to static designs. Our results lay the groundwork for the automated design of control circuits for pathway engineering, with applications in the food, energy and pharmaceutical sectors. Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ-9-tetrahydrocannabinol In the light of legislative changes and the widespread use of cannabis as a recreational and medicinal drug, delayed effects of cannabis upon brief exposure during embryonic development are of high interest as early pregnancies often go undetected. Here, zebrafish embryos were exposed to cannabidiol (CBD) and {Delta}-9-tetrahydrocannabinol l (THC) until the end of gastrulation (1-10 hours post-fertilization) and analyzed later in development (4-5 days post-fertilization). In order to measure neural activity, we implemented CaMPARI (Calcium-Modulated Photoactivatable Ratiometric Integrator) and optimized the protocol for a 96-well format complemented by locomotor analysis. Our results revealed that neural activity was decreased by CBD more than THC. At higher doses, both cannabinoids could dramatically reduce neural activity and locomotor activity. Interestingly, the decrease was more pronounced when CBD and THC were combined. At the receptor level, CBD-mediated reduction of locomotor activity was partially prevented using cannabinoid type 1 and 2 receptor inhibitors. Overall, we report that CBD toxicity occurs via two cannabinoid receptors and is synergistically enhanced by THC exposure to negatively impact neural activity late in larval development. Future studies are warranted to reveal other cannabinoids and receptors involved in this pathway to understand the subsequent health implications of cannabis consumption on fetal development. A mixture of innate cryoprotectants is key for freeze tolerance and cryopreservation of a drosophilid fly larva. Insects that naturally tolerate internal freezing produce complex mixtures of multiple cryoprotectants (CPs). Better knowledge on composition of these mixtures, and on mechanisms of how the individual CPs interact, could inspire development of laboratory CP formulations optimized for cryopreservation of cells and other biological material. Here we identify and quantify (using high resolution mass spectrometry) a range of putative CPs in larval tissues of a subarctic fly, Chymomyza costata that survives long-term cryopreservation in liquid nitrogen. The CPs (proline, trehalose, glutamine, asparagine, glycine betaine, glycerophosphoethanolamine, glycerophosphocholine, and sarcosine) accumulate in hemolymph in a ratio of 313:108:55:26:6:4:3:0.5 mmol.L-1. Using calorimetry, we show that the artificial mixtures, mimicking the concentrations of major CPs in hemolymph of freeze-tolerant larvae, suppress the melting point of water and significantly reduce the ice fraction. We demonstrate in a bioassay that mixtures of CPs administered through the diet act synergistically rather than additively to enable cryopreservation of otherwise freeze-sensitive larvae. Using MALDI-MSI, we show that during slow extracellular freezing trehalose becomes concentrated in partially dehydrated hemolymph where it stimulates transition to the amorphous glass phase. In contrast, proline moves to the boundary between extracellular ice and dehydrated hemolymph and tissues where it likely forms a layer of dense viscoelastic liquid. We propose that amorphous glass and viscoelastic liquids may protect macromolecules and cells from thermomechanical shocks associated with freezing and transfer into and out of liquid nitrogen.

Summary statementThe composition of natural cryoprotectant mixture of the extremely freeze-tolerant insect is revealed. Components of the mixture work in synergy and behave differently during organismal freezing and cryopreservation. Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin BackgroundEvolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation at many enhancers, distal regulatory elements that control spatio-temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific, code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

ResultsWe present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that is able to predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Then, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity and liver enhancers associated with liver regeneration have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

ConclusionThe framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements. Frequency modulated timer regulates mammalian hibernation Mammalian hibernators decrease basal metabolism and body temperature (Tb) to minimize energy expenditure in harsh seasons. During hibernation, Tb drops to low temperature (<10 {degrees}C) and remains constant for days, known as deep torpor in small mammalian hibernators. Spontaneous interbout arousals interrupt torpor bouts, when Tb recovers to euthermic state ~37 {degrees}C. Torpor-interbout arousal event repeats during several months of hibernation. However, little is known about mechanisms governing Tb fluctuation across torpor-interbout arousal cycles during hibernation. Recent improvement in data-logging techniques enables us to monitor Tb for more than hundred days with high precision, opening up new avenues for quantitative analysis to address the principle governing Tb fluctuation. Here, we analyzed Tb fluctuation across torpor-interbout arousal cycle of Syrian hamster, which can hibernate in response to chronic cold and short photoperiod under a laboratory condition, using generalized harmonic analysis and discovered a model with frequency modulation quantitatively reproducing Tb fluctuation. This analysis also identified an unexpectedly longer period of 120-430 days as the period that modulates another period of several days, generating Tb fluctuation for Syrian hamster. We propose that concerted action of two endogenous periods governs torpor-interbout arousal cycles during hibernation. Tatajuba -- Exploring the distribution of homopolymer tracts Length variation of homopolymeric tracts, which induces phase variation, is known to regulate gene expression leading to phenotypic variation in a wide range of bacterial species. There is no specialised bioinformatics software which can, at scale, exhaustively explore and describe these features from sequencing data. Identifying these is non-trivial as sequencing and bioinformatics methods are prone to introducing artefacts when presented with homopolymeric tracts due to the decreased base diversity. We present tatajuba, which can automatically identify potential homopolymeric tracts and their putative phenotypic impact, allowing for rapid investigation. We use it to detect all tracts in two separate datasets, one of Campylobacter jejuni and one of three Bordetella species, and to highlight those tracts that are polymorphic across samples. With this we confirm homopolymer tract variation with phenotypic impact found in previous studies and additionally find many more with potential variability. The software is written in C and is available under the open source license GNU GPL version 3 from https://github.com/quadram-institute-bioscience/tatajuba. Variable susceptibility of intestinal organoid-derived monolayers to SARS-CoV-2 infection Gastrointestinal effects associated with COVID-19 are highly variable for reasons that are not understood. In this study, we used intestinal organoid-derived cultures differentiated from primary human specimens as a model to examine inter-individual variability. Infection of intestinal organoids derived from different donors with SARS-CoV-2 resulted in orders of magnitude differences in virus replication in small intestinal and colonic organoid-derived monolayers. Susceptibility to infection correlated with ACE2 expression level and was independent of donor demographic or clinical features. ACE2 transcript levels in cell culture matched the amount of ACE2 in primary tissue indicating this feature of the intestinal epithelium is retained in the organoids. Longitudinal transcriptomics of organoid-derived monolayers identified a delayed yet robust interferon signature, the magnitude of which corresponded to the degree of SARS-CoV-2 infection. Interestingly, virus with the Omicron variant spike protein infected the organoids with the highest infectivity, suggesting increased tropism of the virus for intestinal tissue. These results suggest that heterogeneity in SARS-CoV-2 replication in intestinal tissues results from differences in ACE2 levels, which may underlie variable patient outcomes. Quantifying the Morphology and Mechanisms of Cancer Progression in 3D in-vitro environments: Integrating Experiments and Multiscale Models Mathematical models of cancer growth have become increasingly more accurate both in the space and time domains. However, the limited amount of data typically available has resulted in a larger number of qualitative rather than quantitative studies. In the present study, we provide an integrated experimental-computational framework for the quantification of the morphological characteristics and the mechanistic modelling of cancer progression in 3D environments. The proposed framework allows for the calibration of multiscale, spatiotemporal models of cancer growth using state-of-the-art 3D cell culture data, and their validation based on the resulting experimental morphological patterns. Its implementation enabled us to pursue two goals; first, the quantitative description of the morphology of cancer progression in 3D cultures, and second, the relation of tumour morphology with underlying biophysical mechanisms that govern cancer growth and migration. We applied this framework to the study of the spatiotemporal progression of Triple Negative Breast Cancer cells cultured in Matrigel scaffolds, and validated the hypothesis of chemotactic migration using a multiscale Keller-Segel model. The results revealed transient, non-random spatial distributions of cancer cells that consist of clustered, and dispersion patterns. The proposed model was able to describe the general characteristics of the experimental observations and suggests that cancer cells exhibited chemotactic migration and accumulation, as well as random motion during the examined time period of development. To our knowledge, this is the first time that a multiscale model is used to quantify the relationship between the spatial patterns and the underlying mechanisms of cancer growth in 3D environments.