No evidence of human genome integration of SARS-CoV-2 found by long-read DNA sequencing A recent study proposed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hijacks the LINE-1 (L1) retrotransposition machinery to integrate into the DNA of infected cells. If confirmed, this finding could have significant clinical implications. Here, we applied deep (>50x) long-read Oxford Nanopore Technologies (ONT) sequencing to HEK293T cells infected with SARS-CoV-2, and did not find the virus integrated into the genome. By examining ONT data from separate HEK293T cultivars, we completely resolved 78 L1 insertions arising in vitro in the absence of L1 overexpression systems. ONT sequencing applied to hepatitis B virus (HBV) positive liver cancer tissues located a single HBV insertion. These experiments demonstrate reliable resolution of retrotransposon and exogenous virus insertions via ONT sequencing. That we found no evidence of SARS-CoV-2 integration suggests such events are, at most, extremely rare in vivo, and therefore are unlikely to drive oncogenesis or explain post-recovery detection of the virus.
Cell type-specific inference of differential expression in spatial transcriptomics Spatial transcriptomics enables spatially resolved gene expression measurements at near single-cell resolution. There is a pressing need for computational tools to enable the detection of genes that are differentially expressed (DE) within specific cell types across tissue context. We show that current approaches cannot learn cell type-specific DE due to changes in cell type composition across space and the fact that measurement units often detect transcripts from more than one cell type. Here, we introduce a statistical method, Cell type-Specific Inference of Differential Expression (C-SIDE), that identifies cell type-specific patterns of differential gene expression while accounting for localization of other cell types. We model spatial transcriptomics gene expression as an additive mixture across cell types of general log-linear cell type-specific expression functions. This approach provides a unified framework for defining and identifying gene expression changes in a wide-range of relevant contexts: changes due to pathology, anatomical regions, physical proximity to specific cell types, and cellular microenvironment. Furthermore, our approach enables statistical inference across multiple samples and replicates when such data is available. We demonstrate, through simulations and validation experiments on Slide-seq and MER-FISH datasets, that our approach accurately identifies cell type-specific differential gene expression and provides valid uncertainty quantification. Lastly, we apply our method to characterize spatially-localized tissue changes in the context of disease. In an Alzheimers mouse model Slide-seq dataset, we identify plaque-dependent patterns of cellular immune activity. We also find a putative interaction between tumor cells and myeloid immune cells in a Slide-seq tumor dataset. We make our C-SIDE method publicly available as part of the open source R package https://github.com/dmcable/spacexr.
The obstacles and potential clues of prime editing applications in tomato, a dicot plant Precision genome editing is highly desired for crop improvement. The recently emerged CRISPR/Cas technology offers great potential applications in precision plant genome engineering. A prime editing (PE) approach combining a reverse transcriptase (RT) with a Cas9 nickase and a "priming" extended guide RNA has shown a high frequency for precise genome modification in mammalian cells and several plant species. Nevertheless, the applications of the PE approach in dicot plants are still limited and inefficient. We designed and tested prime editors for precision editing of a synthetic sequence in a transient assay and for desirable alleles of 10 loci in tomato by stable transformation. Our data obtained by targeted deep sequencing also revealed only low PE efficiencies in both the tobacco and tomato systems. Further assessment of the activities of the PE components uncovered potential reasons for the inefficiency of the PE complexes. Modifying the pegRNA sequences by shortening or introducing mismatches to the primer binding sequences (PBS) in order to reduce their melting temperatures (Tm) did not enhance the PE efficiency at the SlBMP21, SlALC and SlALS1 loci. Our data show challenges of PE approach in tomato, indicating a further improvement of the PE system for the successful applications such as use of improved expression systems. Our work provides an important clue for the successful application of the PE approach in crop improvement. Catchet-MS identifies IKZF1-targeting Thalidomide analogues as novel HIV-1 latency reversal agents A major pharmacological strategy toward HIV cure aims to reverse latency in infected cells as a first step leading to their elimination. While the unbiased identification of molecular targets physically associated with the latent HIV-1 provirus would be highly valuable to unravel the molecular determinants of HIV-1 transcriptional repression and latency reversal, due to technical limitations, this has not been possible. Here we use a dCas9 targeted chromatin and histone enrichment strategy coupled to mass spectrometry (Catchet-MS) to describe the protein composition of the latent and activated HIV-1 5LTR. Catchet-MS identified known and novel latent 5LTR-associated host factors. Among these, IKZF1 is a novel HIV-1 transcriptional repressor, required for Polycomb Repressive Complex 2 recruitment to the LTR. We find the clinically advanced thalidomide analogue iberdomide, and the FDA approved analogues lenalidomide and pomalidomide, to be novel LRAs that, by targeting IKZF1 for degradation, reverse HIV-1 latency in CD4+T-cells isolated from virally suppressed people living with HIV-1.

One Sentence SummarydCas9 targeted chromatin and histone enrichment for mass spectrometry (Catchet-MS) led to the identification of IKZF1-targeting thalidomide analogues as novel HIV-1 latency reversal agents Rec8 cohesin-mediated axis-loop chromatin architecture is required for meiotic recombination During meiotic prophase, cohesin-dependent axial structures are formed in the synaptonemal complex (SC). However, the functional correlation between these structures and cohesion remains elusive. Here, we examined the formation of cohesin-dependent axial structures in the fission yeast Schizosaccharomyces pombe. This organism forms atypical SCs composed of linear elements (LinEs) resembling the lateral elements of SC but lacking the transverse filaments. Hi-C analysis using a highly synchronous population of meiotic S. pombe cells revealed that the axis-loop chromatin structure formed in meiotic prophase was dependent on the Rec8 cohesin complex. In contrast, the Rec8-mediated formation of the axis-loop structure occurred in cells lacking components of LinEs. To dissect the functions of Rec8, we identified a rec8-F204S mutant that lost the ability to assemble the axis-loop structure without losing cohesion of sister chromatids. This mutant showed defects in the formation of the axis-loop structure and LinE assembly and thus exhibited reduced meiotic recombination. Collectively, our results demonstrate that the Rec8-dependent axis-loop structure provides a structural platform essential for LinE assembly, facilitating meiotic recombination of homologous chromosomes, independently of its role in sister chromatid cohesion. Spatial transcriptomic characterization of COVID-19 pneumonitis identifies immune circuits related to tissue injury Severe lung damage in COVID-19 involves complex interactions between diverse populations of immune and stromal cells. In this study, we used a spatial transcriptomics approach to delineate the cells, pathways and genes present across the spectrum of histopathological damage in COVID-19 lung tissue. We applied correlation network-based approaches to deconvolve gene expression data from areas of interest within well preserved post-mortem lung samples from three patients. Despite substantial inter-patient heterogeneity we discovered evidence for a common immune cell signaling circuit in areas of severe tissue that involves crosstalk between cytotoxic lymphocytes and pro-inflammatory macrophages. Expression of IFNG by cytotoxic lymphocytes was associated with induction of chemokines including CXCL9, CXCL10 and CXCL11 which are known to promote the recruitment of CXCR3+ immune cells. The tumour necrosis factor (TNF) superfamily members BAFF (TNFSF13B) and TRAIL (TNFSF10) were found to be consistently upregulated in the areas with severe tissue damage. We used published spatial and single cell SARS-CoV-2 datasets to confirm our findings in the lung tissue from additional cohorts of COVID-19 patients. The resulting model of severe COVID-19 immune-mediated tissue pathology may inform future therapeutic strategies.

One Sentence SummarySpatial analysis identifies IFN{gamma} response signatures as focal to severe alveolar damage in COVID-19 pneumonitis. Improved AlphaFold modeling with implicit experimental information Machine learning prediction algorithms such as AlphaFold1 and RoseTTAFold2 can create remarkably accurate protein models, but these models usually have some regions that are predicted with low confidence or poor accuracy3-6. We hypothesized that by implicitly including experimental information, a greater portion of a model could be predicted accurately, and that this might synergistically improve parts of the model that were not fully addressed by either machine learning or experiment alone. An iterative procedure was developed in which AlphaFold models are automatically rebuilt based on experimental density maps and the rebuilt models are used as templates in new AlphaFold predictions. We find that including experimental information improves prediction beyond the improvement obtained with simple rebuilding guided by the experimental data. This procedure for AlphaFold modeling with density has been incorporated into an automated procedure for crystallographic and electron cryo-microscopy map interpretation. Rho Kinase regulates neutrophil NET formation that is involved in UVB-induced skin inflammation ObjectiveUltraviolet B (UVB) is an important trigger of skin inflammation and lupus with leukocyte recruitment to inflamed skin. We recently reported the involvement of neutrophil NETosis in UVB-induced skin inflammation, and that NETotic nuclear envelope rupture is driven by PKC-mediated nuclear lamin B disassembly. To address the role of Actin cytoskeleton in NETosis, we investigated the effects of Rho kinase (ROCK) and its downstream actomyosin cytoskeletal networks on PKC nuclear translocation and NET formation, as well as their involvement in UVB-induced skin inflammation.

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

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

ConclusionsROCK regulated NETosis through modulation of PKC nuclear translocation via actomyosin cytoskeletal networks in neutrophils. ROCK1 deficiency ameliorated UVB- induced skin inflammation by attenuation of NETosis and NET-associated cytokines. Microbial colonization of tannin-rich tropical plants: interplay between degradability, methane production and tannin disappearance in the rumen AO_SCPLOWBSTRACTC_SCPLOWCondensed tannins in plants are found free and attached to protein and fibre but it is not known whether these fractions influence rumen degradation and microbial colonization. This study explored the rumen degradation of tropical tannins-rich plants and the relationship between condensed tannins fractions disappearance and microbial communities colonising plant particles. Leaves from Calliandra calothyrsus, Gliricidia sepium, and Leucaena leucocephala, pods from Acacia nilotica and the leaves of two agricultural by-products: Manihot esculenta and Musa spp. were incubated in situ in the rumen of dairy cows. An in vitro approach was also used to assess the effects of these plants on rumen fermentation parameters. All plants contained more than 100 g/kg of condensed tannins with a large proportion (32 to 61%) bound to proteins. Calliandra calothyrsus had the highest concentration of condensed tannins at 361 g/kg, whereas Acacia nilotica was particularly rich in hydrolysable tannins (350 g/kg). Hydrolysable and free condensed tannins from all plants completely disappeared after 24 h incubation in the rumen. Disappearance of protein-bound condensed tannins was variable with values ranging from 93% for Gliricidia sepium to 21% for Acacia nilotica. In contrast, fibre-bound condensed tannins disappearance averaged ~82% and did not vary between plants. Disappearance of bound fractions of condensed tannins was not associated with degradability of plant fractions. The presence of tannins interfered with the microbial colonisation of plants. Each plant had distinct bacterial and archaeal communities after 3 and 12 h of incubation in the rumen and distinct protozoal communities at 3 h. Adherent communities in tannin-rich plants had a lower relative abundance of fibrolytic microbes, notably Fibrobacter spp. Whereas, archaea diversity was reduced in high tannin-containing Calliandra calothyrsus and Acacia nilotica at 12 h of incubation. Concurrently, in vitro methane production was lower for Calliandra calothyrsus, Acacia nilotica and Leucaena leucocephala although for the latter total volatile fatty acids production was not affected and was similar to control. Here we show that the total amount of hydrolysable and condensed tannins contained in a plant govern the interaction with rumen microbes affecting degradability and fermentation. The effect of protein- and fibre-bound condensed tannins on degradability is less important. Multiple spillovers and onward transmission of SARS-Cov-2 in free-living and captive White-tailed deer Many animal species are susceptible to SARS-CoV-2 and could potentially act as reservoirs, yet transmission of the virus in non-human free-living animals has not been documented. White-tailed deer (Odocoileus virginianus), the predominant cervid in North America, are susceptible to SARS-CoV-2 infection, and experimentally infected fawns can transmit the virus. To test the hypothesis that SARS-CoV-2 may be circulating in deer, we tested 283 retropharyngeal lymph node (RPLN) samples collected from 151 free-living and 132 captive deer in Iowa from April 2020 through December of 2020 for the presence of SARS-CoV-2 RNA. Ninety-four of the 283 deer (33.2%; 95% CI: 28, 38.9) samples were positive for SARS-CoV-2 RNA as assessed by RT-PCR. Notably, between November 23, 2020 and January 10, 2021, 80 of 97 (82.5%; 95% CI 73.7, 88.8) RPLN samples had detectable SARS-CoV-2 RNA by RT-PCR. Whole genome sequencing of the 94 positive RPLN samples identified 12 SARS-CoV-2 lineages, with B.1.2 (n = 51; 54.5%), and B.1.311 (n = 19; 20%) accounting for ~75% of all samples. The geographic distribution and nesting of clusters of deer and human lineages strongly suggest multiple zooanthroponotic spillover events and deer-to-deer transmission. The discovery of sylvatic and enzootic SARS-CoV-2 transmission in deer has important implications for the ecology and long-term persistence, as well as the potential for spillover to other animals and spillback into humans. These findings highlight an urgent need for a robust and proactive "One Health" approach to obtaining a better understanding of the ecology and evolution of SARS-CoV-2.

One-Sentence SummarySARS-CoV-2 was detected in one-third of sampled white-tailed deer in Iowa between September 2020 and January of 2021 that likely resulted from multiple human-to-deer spillover and deer-to-deer transmission events. 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. Multi-layered regulation of neuroectoderm differentiation by retinoic acid in a primitive streak-like context The formation of the primitive streak (PS) and the subsequent induction of neuroectoderm are hallmarks of gastrulation. Combining an in vitro reconstitution of this process based on mouse embryonic stem cells (mESCs) with a collection of knockouts in reporter mESC lines, we reassessed the contribution of retinoic acid (RA) signaling at early stages of neural commitment and its cross-talk with TGF{beta} and Wnt signaling inhibition. Single-cell RNA sequencing analysis captured the temporal unfolding of cell type diversification from epiblast and primitive streak-like cells up to the emergence of anterior and posterior neural fates. In conditions thought to lack RA synthesis, we discovered a hitherto unidentified residual RA production via a sensitive RA reporter. Genetic perturbations proved that the RA-degrading enzyme Cyp26a1 safeguard the developmental capabilities of the PS-like cells, limiting neural differentiation in wild type and in Chrd-/-Nog-/- or Dkk1-/- cells. Finally, the knockout of the three RAR receptors highlighted their function as negative regulators of loci critical for neural induction. Overall, we identified two mechanisms whereby components of the RA pathway can control the formation of neural progenitors in our PS-like context: a RA-dependent neural induction gated by RARs, and a receptor-mediated repression in the absence of ligand. Linked machine learning classifiers improve species classification of fungi when using error-prone long-reads on extended metabarcodes BackgroundThe increased usage of error-prone long-read sequencing for metabarcoding of fungi has not been matched with adequate public databases and concomitant analysis approaches. We address this gap and present a proof-of-concept study for classifying fungal taxa using linked machine learning classifiers. We demonstrate the capability of linked machine learning classifiers to accurately classify species and strains using real-world and simulated fungal ribosomal DNA datasets, including plant and human pathogens. We benchmark our new approach in comparison to current alignment and k-mer based methods based on synthetic mock communities. We also assess real world applications of species identification in complex unlabelled datasets.

ResultsOur machine learning approach assigned individual nanopore long-read amplicon sequences to fungal species with high recall rates and low false positive rates. Importantly, our approach successfully distinguished between closely-related species and strains when individual read errors were higher than the genetic distance between individual taxa, which the alignment and k-mer methods could not do. The machine learning approach showed an ability to identify key species with high recall rates, even in complex samples of unknown species composition.

ConclusionsA proof of concept machine learning approach using a tree-descent approach on a decision tree of classifiers can identify known taxa with high accuracy, and precisely detect known target species from complex samples with high recall rates. We propose this approach is suitable for detecting the known knowns of pathogens or invasive species in any environment of mostly unknown composition, including agriculture and wild ecosystems. Anaplerotic flux into the Calvin-Benson cycle. Hydrogen isotope evidence for in vivo occurrence in C3 metabolism - As the central carbon uptake pathway in photosynthetic cells, the Calvin-Benson cycle is among the most important biochemical cycles for life on Earth. A carbon flux of anaplerotic origin (i.e., through the chloroplast-localised oxidative branch of the pentose phosphate pathway) into the Calvin-Benson cycle was proposed recently.
- Here, we measured intramolecular deuterium abundances in leaf starch of Helianthus annuus grown at varying ambient CO2 concentrations, Ca. Additionally, we modelled deuterium fractionations expected for the anaplerotic pathway and compared modelled with measured fractionations.
- We report deuterium fractionation signals at H1 and H2 of starch glucose. Below a Ca change point, these signals increase with decreasing Ca consistent with modelled fractionations by anaplerotic flux. Under standard conditions (Ca=450 ppm corresponding to intercellular CO2 concentrations, Ci, of 328 ppm), we estimate negligible anaplerotic flux. At Ca=180 ppm (Ci=140 ppm), more than 10% of the glucose 6-phosphate entering the starch biosynthesis pathway is diverted into the anaplerotic pathway.
- In conclusion, we report evidence consistent with anaplerotic carbon flux into the Calvin-Benson cycle in vivo. We propose the flux may help to (i) maintain high levels of ribulose 1,5-bisphosphate under source-limited growth conditions to facilitate photorespiratory nitrogen assimilation required to build-up source strength and (ii) counteract oxidative stress. A theoretical model of mitochondrial ATP synthase deficiencies. The role of mitochondrial carriers. The m.8993T>G mutation of the mitochondrial MT-ATP6 gene is associated with NARP syndrome (Neuropathy, Ataxia and Retinitis Pigmentosa). The equivalent point mutation introduced in yeast Saccharomyces cerevisiae mitochondrial DNA considerably reduced the activity of ATP synthase and of cytochrome-c-oxidase preventing yeast growth on oxidative substrates. The over-expression of the mitochondrial oxodicarboxylate carrier (Odc1p) is able to rescue the growth on oxidative substrate by increasing the substrate-level phosphorylation of ADP coupled to conversion of -ketoglutarate (AKG) into succinate with an increase in Complex IV activity. Previous studies showed that equivalent point mutations in ATP synthase behave similarly and can be rescued by Odc1p overexpression and/or uncoupler of OXPHOS from ATP synthesis. In order to better understand the mechanism of ATP synthase mutation bypass, we developed a core model of mitochondrial metabolism based on AKG as respiratory substrate. We describe the different possible metabolite output and the ATP/O ratio values as a function of ATP synthase inhibition. On Rank Deficiency in Phenotypic Covariance This paper is concerned with rank deficiency in phenotypic covariance matrices: first to establish it is a problem by measuring it, and then proposing methods to treat for it. Significant rank deficiency can mislead current measures of whole-shape phenotypic integration, because they rely on eigenvalues of the covariance matrix, and highly rank deficient matrices will have a large percentage of meaningless eigenvalues. This paper has three goals. The first is to examine a typical geometric morphometric data set and establish that its covariance matrix is rank deficient. We employ the concept of information, or Shannon, entropy to demonstrate that a sample of dire wolf jaws is highly rank deficient. The different sources of rank deficiency are identified, and include the Generalized Procrustes analysis itself, use of the correlation matrix, insufficient sample size, and phenotypic covariance. Only the last of these is of biological interest.

Our second goal is to examine a test case where a change in integration is known, allowing us to document how rank deficiency affects two measures of whole shape integration (eigenvalue standard deviation and standardized generalized variance). This test case utilizes the dire wolf data set from Part 1, and introduces another population that is 5000 years older. Modularity models are generated and tested for both populations, showing that one population is more integrated than the other. We demonstrate that eigenvalue variance characterizes the integration change incorrectly, while the standardized generalized variance lacks sensitivity. Both metrics are impacted by the inclusion of many small eigenvalues arising from rank deficiency of the covariance matrix. We propose a modification of the standardized generalized variance, again based on information entropy, that considers only the eigenvalues carrying non-redundant information. We demonstrate that this metric is successful in identifying the integration change in the test case.

The third goal of this paper is to generalize the new metric to the case of arbitrary sample size. This is done by normalizing the new metric to the amount of information present in a permuted covariance matrix. We term the resulting metric the  relative dispersion, and it is sample size corrected. As a proof of concept we us the new metric to compare the dire wolf data set from the first part of this paper to a third data set comprising jaws of Smilodon fatalis. We demonstrate that the Smilodon jaw is much more integrated than the dire wolf jaw. Finally, this information entropy-based measures of integration allows comparison of whole shape integration in dense semilandmark environments, allowing characterization of the information content of any given shape, a quantity we term  latent dispersion. Experimental evolution in TP53 deficient human gastric organoids recapitulates tumorigenesis The earliest events during human tumor initiation are poorly characterized but may hold clues as to how to detect and prevent malignancy. Here we model occult pre-neoplasia by engineering TP53 deficiency, a common and early event in gastric cancer, into human gastric organoids. By performing experimental evolution in multiple clonally derived cultures over two years we define causal relationships between this common initiating genetic lesion and resulting phenotypes. TP53 loss elicited progressive aneuploidy, including copy number alterations and complex structural variants that are common in gastric cancers and which follow preferred temporal orders. Longitudinal single cell sequencing of TP53 deficient gastric organoids similarly indicates progression towards malignant transcriptional programs. Moreover, lineage tracing with expressed cellular barcodes demonstrates reproducible dynamics whereby initially rare subclones with shared transcriptional programs repeatedly attain clonal dominance. This powerful platform for experimental evolution exposes stringent selection, clonal interference and a striking degree of phenotypic convergence in pre-malignant epithelial organoids, implying that the earliest stages of tumorigenesis may be predictable while illuminating evolutionary constraints and barriers to malignant transformation.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/487529v2_ufig1.gif" ALT="Figure 1">
View larger version (63K):
org.highwire.dtl.DTLVardef@122debborg.highwire.dtl.DTLVardef@4cd2daorg.highwire.dtl.DTLVardef@1a720d6org.highwire.dtl.DTLVardef@15ccaa7_HPS_FORMAT_FIGEXP  M_FIG C_FIG 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. Quantifying the effect of caloric and non-caloric sweeteners in the brain response using EEG and convolutional neural network Sweetener type can influence sensory properties and consumers acceptance and preference for low-calorie products. An ideal sweetener does not exist, and each sweetener must be used in situations to which it is best suited. Aspartame and sucralose can be good substitutes for sucrose in passion fruit juice. Despite the interest in artificial sweeteners, little is known about how artificial sweeteners are processed in the human brain. Here, we evaluated brain signals of 11 healthy subjects when they tasted passion fruit juice equivalently sweetened with sucrose (9.4 g/100 g), sucralose (0.01593 g/100 g), or aspartame (0.05477 g/100 g). Electroencephalograms were recorded for two sites in the gustatory cortex (i.e., C3 and C4). Data with artifacts were disregarded, and the artifact-free data were used to feed a CNN. Our results indicated that the brain responses distinguish juice sweetened with different sweeteners with an average accuracy of 0.823.

Practical ApplicationsFinding sweeteners that best fit consumer preferences evolves understanding how the gustatory cortex processes sweeteners. Ideal equivalence will occur when the brain is no longer able to distinguish stimuli that are consciously perceived. This study presents a method of signal acquisition using a single channel and an open-source processing environment. This would allow, for example, to disregard the use of a commercial electroencephalograph and expand the studies in this area and offering to food industry additional tools in the development of products sweetened with non-caloric sweeteners. 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. Classification of Pacemaker Dynamics in the Mouse Intestine by Field Potential Microimaging The flexible and sophisticated movement of the gastrointestinal (GI) tract implies the involvement of mechanisms other than enteric neural circuits, to coordinate excitation in microregions. We thus performed microimaging of pacemaker dynamics in the small intestine of mice since it contains typical network-forming pacemaker cells. A dialysis membrane-reinforced low-impedance microelectrode array (MEA) enabled field potentials over a wide frequency range to be stably measured in microregions. The pacemaker dynamics were classified into basic patterns despite large variations. In the developmental process, pacemaker activity was categorized as either an  expanding or a  migrating pattern that was initiated in or propagated to the MEA sensing area, respectively. The intercellular current of the volume conductor complicated the waveform of both activities. The existence of  expanding and  migrating patterns was attributable to duplicated pacemaker systems such as intracellular Ca2+ oscillation-activated and voltage-gated mechanisms. Additionally, from the spatio-temporal feature during the period of pacemaker events, the  bumpy/aberrant pattern was defined by aberrant, incoherent propagation, and associated with local impairment of excitability, while the  colliding/converging pattern involved the interaction of multiple activities in the MEA area. Interconversion between the four micro-coordination patterns occurred in the same microregion. 5-Hydroxytryptamine (5-HT) promoted  migrating activity, implying an improvement or restoration of spatial conductivity. These results agree well with the action of 5-HT to change GI movement toward propulsion. In conclusion, our MEA method of microimaging classification enables the quantitative assessment of spatio-temporal electric coordination underlying GI motility, suggesting its application to small model animals. Genetic architecture of natural variation of cardiac performance in flies Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP)1. Genome Wide Associations Studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict Transcription Factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans. A general theory for temperature-dependence in biology At present, there is no simple, complete, and first principles-based model for quantitatively describing the full range of observed biological temperature responses. Here, we derive a theory exhibiting these features based on the Eyring-Evans-Polanyi theory governing chemical reaction rates, and which is applicable across all scales from the micro to the macro. Assuming only that the conformational entropy of molecules changes with temperature, we derive a theory for the temperature dependence which takes the form of an exponential function modified by a power-law. Our framework leads to six deductions applicable to any biological trait that depends on temperature, and elucidates novel aspects of universal temperature responses across the tree of life, from quantum to classical scales. All predictions are well supported by data for a wide variety of biological rates and steady states, from molecular to ecological scales and across multiple taxonomic groups. In addition, we provide novel explanations of several empirical relationships including optimal values in temperature response curves.

One-Sentence SummaryWe derive a simple and universal formulae to characterize temperature responses of biological processes across the tree of life. Temporal Scaling of Human Scalp-Recorded Potentials Much of human behaviour is governed by common processes that unfold over varying timescales. Standard event-related potential analysis assumes fixed-latency responses relative to experimental events. However, recent single unit recordings in animals have revealed neural activity scales to span different durations during behaviours demanding flexible timing. Here, we employed a general linear modelling approach using a novel combination of fixed-duration and variable-duration regressors to unmix fixed-time and scaled-time components in human magneto/electroencephalography (M/EEG) data. We use this to reveal consistent temporal scaling of human scalp-recorded potentials across four independent EEG datasets, including interval perception, production, prediction and value-based decision making. Between-trial variation in the temporally scaled response predicts between-trial variation in subject reaction times, demonstrating the relevance of this temporally scaled signal for temporal variation in behaviour. Our results provide a general approach for studying flexibly timed behaviour in the human brain. Adaptive Wild Introgression Mediate Temporal Wheat Root Water Influxes under Drought Stress through Plant-Atmospheric Interrelations Drought intensity as experienced by plants depends upon soil moisture status and atmospheric variables such as temperature, radiation, and air vapour pressure deficit (VPD). Although the role of shoot architecture with these edaphic and atmospheric factors is well-characterized, the extent to which shoot and root dynamic interactions as a continuum are controlled by genotypic variation is less known. Here, we targeted these interactions using a wild emmer introgression line (IL20) with a distinct drought-induced shift in the shoot-to-root ratio and its drought-sensitive recurrent parent Svevo. Using a gravimetric platform, we show that IL20 maintained higher root water influx and gas exchange under drought stress, which supported a greater growth. Interestingly, the advantage of IL20 in root water influx and transpiration was expressed earlier during the daily diurnal cycle under lower VPD and therefore supported higher transpiration efficiency. Application of structural equation model indicates that under drought, VPD and radiation are antagonistic to transpiration rate, whereas the root water influx operates as feedback for the higher atmospheric responsiveness of leaves. Collectively, our results suggest that a drought-induced shift in root-to-shoot ratio can improve plant water uptake potential in a short preferable time window determined by both water and atmospheric parameters. 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.