A high throughput lipidomics method using scheduled multiple reaction monitoring Lipid compositions of cells, tissues and bio-fluids are complex, with varying concentrations and structural diversity, which makes their identification challenging. Newer methods for comprehensive analysis of lipids are thus necessary. Herein, we propose a targeted-mass spectrometry based method for large-scale lipidomics using a combination of variable retention time window and relative dwell time weightage. Using this, we detected more than 1000 lipid species, including structural isomers. The limit of detection varied from femtomolar to nanomolar range and the coefficient of variance <30% for 849 lipid species. We used this method to identify lipids altered due to Vitamin B12 deficiency and found that the levels of lipids with {omega}-3 fatty acid chains decreased while those with {omega}-6 increased. This method enables identification of by far the largest number of lipid species with structural isomers in a single experiment and would significantly advance our understanding of the role of lipids in biological processes.
SUMM4 biochemically complex couples insulator function and DNA replication timingcontrol The mechanisms that establish DNA replication timing programs in eukaryotes remain incompletely understood. Drosophila SNF2-related factor SUUR imparts under-replication (UR) of late-replicating intercalary heterochromatin (IH) in polytene chromosomes. We developed a proteomics technique termed MERCI to isolate a native complex SUMM4 comprising SUUR and chromatin boundary protein Mod(Mdg4)-67.2. Mod(Mdg4) stimulates the ATPase activity of SUUR and is required for its normal spatiotemporal distribution in vivo. Both SUMM4 subunits mediate the activities of gypsy insulator disrupting enhancer-promoter interactions and establishing chromatin barriers. Furthermore, SuUR or mod(mdg4) mutations reverse UR of IH. Our findings uncover a critical role for architectural proteins in attenuating replication fork progression and suggest an alternative mechanism for DNA replication timing that does not depend on an asynchronous firing of replication origins.

One-Sentence SummaryA stable protein complex comprising an insulator factor and a SNF2-like ATPase imparts late replication of heterochromatin.
A Hierarchy of Biomolecular Proportional-Integral-Derivative Feedback Controllers for Robust Perfect Adaptation and Dynamic Performance Proportional-Integral-Derivative (PID) feedback controllers have been the most widely used controllers in industry for almost a century due to their good performance, simplicity, and ease of tuning. Motivated by their success in various engineering disciplines, PID controllers recently found their way into synthetic biology, where the design of feedback molecular control systems has been identified as an important goal. In this paper, we consider the mathematical realization of PID controllers via biomolecular interactions. We propose an array of topologies that offer a compromise between simplicity and high performance. We first demonstrate that different Proportional-Integral (PI) controllers exhibit different capabilities for enhancing the dynamics and reducing variance (cell-to-cell variability). Next, we introduce several derivative controllers that are realized based on incoherent feedforward loops acting in a feedback configuration. Alternatively, we show that differentiators can be realized by placing molecular integrators in a negative feedback loop--an arrangement that can then be augmented by PI components to yield PID feedback controllers. We demonstrate that the derivative component can be exploited for enhancing system stability, dramatically increasing the molecular control systems dynamic performance, and for reducing the noise effect on the output of interest. The PID controller features are established through various deterministic and stochastic analyses as well as numerical simulations. Finally, we provide an experimental demonstration using a recently developed hybrid setup, the cyberloop, where the controller is implemented in silico to control a biological genetic circuit in vivo. The large array of novel biomolecular PID controllers introduced here forms a basis for the design and construction of advanced high-performance biomolecular control systems that robustly regulate the dynamics of living systems. Task-independent metrics of computational hardness predict human cognitive performance The survival of human organisms depends on our ability to solve complex tasks in the face of limited cognitive resources. However, little is known about the factors that drive the complexity of those tasks. Here, building on insights from computational complexity theory, we quantify the computational hardness of cognitive tasks using a set of task-independent metrics related to the computational resource requirements of individual instances of a task. We then examine the relation between those metrics and human behavior and find that they predict both time spent on a task as well as accuracy in three canonical cognitive tasks. Our findings demonstrate that performance in cognitive tasks can be predicted based on generic metrics of their inherent computational hardness.

TeaserThe ability of humans to solve cognitive tasks is affected by generic mathematical properties of problems related to their computational complexity. High school science fair: Experiences of two groups of undergraduate bioscience students Science fairs offer potential opportunities for students to learn first-hand about the practices of science. Over the past six years we have been carrying out voluntary and anonymous surveys with regional and national groups of high school and post high school students to learn about their high school science fair experiences regarding help received, obstacles encountered, and opinions about the value and impact of science fair. Understanding what students think about science fairs will help educators make science fairs more effective learning opportunities. In this paper, we focus on the findings with two national groups of post high school students - undergraduate research fellows (SURF students) who did research at UT Southwestern Medical Center during 2014-2019 and undergraduates biology students attending the 2019 Howard Hughes Medical Institute Science Education Alliance (SEA) summer symposium. About 25% of the students who completed surveys indicated that they had participated in high school science fair, but more than half attended high schools where science fairs were unavailable. Effectively, 6 out of every 10 students participated in science fair if available. Students who could have participated in high school science fair but chose not to do so identified not enough time and coming up with their project idea as major reasons why not. About half the SURF students favored requiring non-competitive science fair regardless whether they themselves had participated in science fair. On the other hand, less than 1 in 5 thought that competitive science fair should be required. Introduction to the scientific process and general learning were mentioned most frequently as the reasons to require non-competitive science fair; these reasons were mentioned rarely in connection with competitive science fair. Unlike the national cohort of high school students we surveyed previously, who mostly did science fair in 9th and 10th grades, SURF students participated in science fair throughout high school and were twice as likely as high school students to have carried out science fair more than once. In conclusion, our findings suggest that participation of the undergraduate bioscience majors in high school science fairs occurs far more frequently than recognized previously. The findings emphasize further the importance of incentivizing rather than requiring science fair participation, especially in 9th and 10th grades, and the potential value of developing non-competitive science fairs. Spontaneous differentiation leads to emergence of hybrid cell states relate to poor prognosis in oral cancer PurposeCellular dynamics between phenotypically heterogeneous subpopulations of cancer cells within individual tumor is shown to be responsible for drug tolerance and overall poor prognosis; however, evidences were largely missing in oral cancer. Therefore, this study was undertaken to describe the dynamic phenotypic states among oral cancer cells, its influence on transcriptomic heterogeneity as well as its clinical significance.

Experimental DesignWe multiplexed phenotypic markers of putative oral-stem-like cancer cells (SLCCs) and characterized diversity among CD44-positive oral cancer cell subpopulations with respect to distinct expression of CD24 and aldehyde dehydrogenase (ALDH)-activity in multiple cell lines. Population trajectories were characterized by Markov model and cell states were defined based on the population specific RNA sequencing (RNAseq). ssGSEA based gene expression signatures were explore for prognostic significance.

ResultsOral cancer cells followed two distinct patterns of spontaneous repopulation dynamics with stochastic inter-conversions on  ALDH-axis, however a strict non-interconvertible transition on  CD24-axis. Interestingly, plastic  ALDH-axis was harnessed to enrich ALDHHigh subpopulations in response to Cisplatin treatment, to adapt a drug tolerant state. Phenotype-specific RNAseq results suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate states of stemness within the differentiating oral cancer cells. Further, survival analysis with each subpopulation-specific gene signature strongly suggested that the cell-state dynamics may act as possible mechanism to drive ITH, resulting in poor prognosis in patient.

ConclusionsOur results emphasized the prognostic power of the population dynamics in oral cancer. Importantly, we have described the phenotypic-composition of heterogeneous subpopulations critical for global tumor behaviour in oral cancer; which is a prerequisite knowledge important for precision treatment, however largely lacking for most solid tumors.

Graphical AbstractWe have characterized diversity among CD44-positive oral cancer cells lines with respect to distinct expression of CD24 and ALDH-activity. Subpopulations showed stochastic inter-conversions on ALDH-axis but a strict non-interconvertible transition of CD24Low to CD24High phenotype, even in response to chemotherapy-induced stress. RNAseq study suggested the organization of subpopulations into hierarchical structure with possible maintenance of intermediate alternate states of stemness within the differentiating oral cancer cells. The described population dynamics demonstrtaed influence tumor behaviour possibly by increasing intratumoral heterogeneity in aggressive oral tumors.



O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/457509v3_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@1dc6a7org.highwire.dtl.DTLVardef@dd0948org.highwire.dtl.DTLVardef@18c9603org.highwire.dtl.DTLVardef@cca8fb_HPS_FORMAT_FIGEXP  M_FIG C_FIG Translational relevanceIntratumoral heterogeneity (ITH) has been the clinically important factor, impacting aggressive cancer behaviour, drug tolerance and overall poor prognosis. Recent high-throughput studies have provided better cellular and molecular resolution of ITH; however, the prerequisite knowledge which defines the composition of subpopulations critical for global tumor behaviour is majorly lacking for most of the solid tumors. By combining phenotypic markers, we have defined four subpopulations of oral cancer cells. These subpopulations showed stochastic inter-conversions as well as a strict non-interconvertible transition among them to acheive heterogeneity. Importantly, transcriptional states of each subpopulations indicated a clinically relevant signatures for patient prognosis. Also, we observed interconversions of these subpopulations in response to Cisplatin to accumulate drug-tolerant cell state, as rapid and reversible strategy to respond to chemotherapy induced stress. Thus, the characteristics of described phenotypic subgroups may be translated to the clinic for estimating the extent of intratumoral heterogeneity in oral cancer patients. Genomic insights into longan evolution from a chromosome-level genome assembly and population genomics of longan accessions Longan (Dimocarpus longan) is a subtropical fruit best known for its nutritious fruit and regarded as a precious tonic and traditional medicine since ancient times. High-quality chromosome-scale genome assembly is valuable for functional genomic study and genetic improvement of longan. Here, we report a chromosome-level reference genome sequence for longan cultivar JDB with an assembled genome of 455.5 Mb in size anchored to fifteen chromosomes, representing a significant improvement of contiguity (contig N50=12.1 Mb, scaffold N50= 29.5 Mb) over a previous draft assembly. A total of 40,420 protein-coding genes were predicted in D. longan genome. Synteny analysis suggests longan shares the widespread gamma event with core eudicots, but has no other whole genome duplications. Comparative genomics showed that D. longan genome experienced significant expansions of gene families related to phenylpropanoid biosynthesis and UDP-glucosyltransferase. Deep genome sequencing analysis of longan cultivars identified longan biogeography as a major contributing factor for genetic diversity, and revealed a clear population admixture and introgression among cultivars of different geographic origins, postulating a likely migration trajectory of longan overall confirmed by existing historical records. Finally, genome-wide association studies (GWAS) of longan cultivars identified quantitative trait loci (QTL) for six different fruit quality traits and revealed a shared QTL containing three genes for total soluble solid and seed weight. The chromosome-level reference genome assembly, annotation and population genetic resource for D. longan will facilitate the molecular studies and breeding of desirable longan cultivars in the future. Streptococcus agalactiae MprF synthesizes a novel cationic glycolipid and contributes to brain entry and meningitis Bacterial membrane lipids are critical for membrane bilayer formation, cell division, protein localization, stress responses, and pathogenesis. Despite their critical roles, membrane lipids have not been fully elucidated for many pathogens. Here, we report the discovery of a novel cationic glycolipid, Lysyl-Glucosyl-Diacylglycerol (Lys-Glc-DAG) that is synthesized in high abundance by the bacterium Streptococcus agalactiae (Group B Streptococcus, GBS). To our knowledge, Lys-Glc-DAG is more positively charged than any other known lipids. Lys-Glc-DAG carries two positive net charges per molecule, distinct from the widely described lysylated phospholipid Lysyl-phosphatidylglycerol (Lys-PG) which carries one positive net charge due to the presence of a negatively charged phosphate moiety. We use normal phase liquid chromatography (NPLC) coupled with electrospray ionization (ESI) high-resolution tandem mass spectrometry (HRMS/MS) and genetic approaches to determine that Lys-Glc-DAG is synthesized by the enzyme MprF in GBS, which covalently modifies the neutral glycolipid Glc-DAG with the cationic amino acid lysine. GBS is a leading cause of neonatal meningitis, which requires traversal of the endothelial blood-brain barrier (BBB). We demonstrate that GBS strains lacking mprF exhibit a significant decrease in the ability to invade BBB endothelial cells. Further, mice challenged with a GBS{Delta}mprF mutant developed bacteremia comparably to Wild-Type infected mice yet had less recovered bacteria from brain tissue and a lower incidence of meningitis. Thus, our data suggest that Lys-Glc-DAG may contribute to bacterial uptake into host cells and disease progression. Importantly, our discovery provides a platform for further study of cationic lipids at the host-pathogen interface. Direct Observation of Heterogeneous Formation of Amyloid Spherulites in Real-time by Super-resolution Microscopy The misfolding of proteins and their aggregation in the form of fibrils or amyloid-like spherulites are involved in a spectrum of pathological abnormalities. Our current understanding of protein amyloid aggregation mechanisms has primarily relied on the use of spectrometric methods to determine the average growth rates and diffraction limited microscopes with low temporal resolution to observe the large-scale morphologies of intermediates. We developed a REal-time kinetics via binding and Photobleaching LOcalisation Microscopy (REPLOM) super-resolution method to directly observe and quantify the existence and abundance of diverse aggregate morphologies below the diffraction limit and extract their heterogeneous growth kinetics. Our results revealed that even the growth of a microscopically identical aggregates, e.g. amyloid spherulites, may follow distinct pathways. Specifically, spherulites do not exclusively grow isotropically but, surprisingly, may also grow anisotropically, following similar pathways as reported for minerals and polymers. Combining our technique with machine learning approaches, we associated growth rates to specific morphological transitions and provided energy barriers and the energy landscape at the level of single aggregate morphology. Our unifying framework for the detection and analysis of spherulite growth can be extended to other self-assembled systems characterized by a high degree of heterogeneity, disentangling the broad spectrum of diverse morphologies at the single-molecule level. A parametric logistic equation with light flux and medium concentration for cultivation planning of microalgae Microalgae are considered to be promising producers of bioactive chemicals, feeds, and fuels from carbon dioxide by photosynthesis. Thus, the prediction of microalgal growth profiles is important for the planning of cost-effective and sustainable cultivation-harvest cycles. This paper proposes a mathematical model capable of predicting the effect of light flux into culture and medium concentration on the growth profiles of microalgae by incorporating these growth-limiting factors into the logistic equation. The specific form of the equation is derived based on the experimentally measured growth profiles of Monoraphidium sp., a microalgal strain isolated by the authors, under 16 conditions consisting of the combinations of incident light fluxes into culture and initial medium concentrations. Using a cross-validation method, it is shown that the proposed model has the ability to predict necessary incident light flux into culture and initial medium concentration for harvesting target biomass at target time. Finally, the model-guided cultivation planning is performed and is evaluated by comparing the result with the experimental data. Algal photosystem I dimer and high resolution model of PSI:plastocyanin complex Photosystem I (PSI) enables photo-electron transfer and regulates photosynthesis in the bioenergetic membranes of cyanobacteria and chloroplasts. Being a multi-subunit complex, its macromolecular organization affects the dynamics of photosynthetic membranes. Here, we reveal a chloroplast PSI from the green alga Chlamydomonas reinhardtii that is organized as a homodimer, comprising 40 protein subunits with 118 transmembrane helices that provide scaffold for 568 pigments. Our cryo-EM structure identifies that the absence of PsaH and Lhca2 gives rise to a head-to-head relative orientation of the PSI-LHCI monomers in a way that is essentially different from the oligomer formation in cyanobacteria. The light-harvesting protein Lhca9 is the key element for mediating this dimerization. The interface between the monomers is lacking PsaH, and thus partially overlaps with the surface area that would bind one of the LHCII complexes in state transitions. We also define the most accurate available PSI-LHCI model at 2.3 [A] resolution, including a flexibly bound electron donor plastocyanin, and assign correct identities and orientations of all the pigments, as well as 621 water molecules that affect energy transfer pathways. A bacterial extracellular vesicle-based intranasal vaccine against SARS-CoV-2 protects against disease and elicits neutralizing antibodies to wild-type virus and Delta variant Several vaccines have been introduced to combat the coronavirus infectious disease-2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Current SARS-CoV-2 vaccines include mRNA-containing lipid nanoparticles or adenoviral vectors that encode the SARS-CoV-2 Spike (S) protein of SARS-CoV-2, inactivated virus, or protein subunits. Despite growing success in worldwide vaccination efforts, additional capabilities may be needed in the future to address issues such as stability and storage requirements, need for vaccine boosters, desirability of different routes of administration, and emergence of SARS-CoV-2 variants such as the Delta variant. Here, we present a novel, well-characterized SARS-CoV-2 vaccine candidate based on extracellular vesicles (EVs) of Salmonella typhimurium that are decorated with the mammalian cell culture-derived Spike receptor-binding domain (RBD). RBD-conjugated outer membrane vesicles (RBD-OMVs) were used to immunize the golden Syrian hamster (Mesocricetus auratus) model of COVID-19. Intranasal immunization resulted in high titers of blood anti-RBD IgG as well as detectable mucosal responses. Neutralizing antibody activity against wild-type and Delta variants was evident in all vaccinated subjects. Upon challenge with live virus, hamsters immunized with RBD-OMV, but not animals immunized with unconjugated OMVs or a vehicle control, avoided body mass loss, had lower virus titers in bronchoalveolar lavage fluid, and experienced less severe lung pathology. Our results emphasize the value and versatility of OMV-based vaccine approaches. Extracellular Matrix Hydrogels Promote Expression of Muscle-Tendon Junction Proteins Muscle and tendon injuries are prevalent and range from minor sprains and strains to traumatic, debilitating injuries. However, the interactions between these tissues during injury and recovery remain unclear. Three-dimensional tissue models that incorporate both tissues and a physiologically relevant junction between muscle and tendon may help understand how the two tissues interact. Here, we use tissue specific extracellular matrix (ECM) derived from muscle and tendon to determine how cells of each tissue interact with the microenvironment of the opposite tissue resulting in junction specific features. ECM materials were derived from the Achilles tendon and gastrocnemius muscle, decellularized, and processed to form tissue specific pre-hydrogel digests. ECM materials were unique in respect to protein composition and included many types of ECM proteins, not just collagens. After digestion and gelation, ECM hydrogels had similar complex viscosities which were less than type I collagen hydrogels at the same concentration. C2C12 myoblasts and tendon fibroblasts were cultured in tissuespecific ECM conditioned media or encapsulated in tissue-specific ECM hydrogels to determine cell-matrix interactions and the effects on a muscle-tendon junction marker, paxillin. ECM conditioned media had only a minor effect on upregulation of paxillin in cells cultured in monolayer. However, cells cultured within ECM hydrogels had 50-70% higher paxillin expression than cells cultured in type I collagen hydrogels. Contraction of the ECM hydrogels varied by the type of ECM used. Subsequent experiments with varying density of type I collagen (and thus contraction) showed no correlation between paxillin expression and the amount of gel contraction, suggesting that a constituent of the ECM was the driver of paxillin expression in the ECM hydrogels. In addition, the extracellular matrix protein type XXII collagen had similar expression patterns as paxillin, with smaller effect sizes. Using tissue specific ECM allowed for the de-construction of the cell-matrix interactions similar to muscletendon junctions to study the expression of MTJ specific proteins.

Impact StatementThe muscle-tendon junction is an important feature of muscle-tendon units; however, despite crosstalk between the two tissue types, it is overlooked in current research. Deconstructing the cell-matrix interactions will provide the opportunity to study significant junction specific features and markers that should be included in tissue models of the muscletendon unit, while gaining a deeper understanding of the natural junction. This research aims to inform future methods to engineer a more relevant multi-tissue platform to study the muscletendon unit. ISL1 is necessary for auditory neuron development and tonotopic organization A cardinal feature of the auditory pathway is frequency selectivity, represented in a tonotopic map from the cochlea to the cortex. The molecular determinants of the auditory frequency map are unknown. Here, we discovered that the transcription factor ISL1 regulates the molecular and cellular features of auditory neurons, including the formation of the spiral ganglion and peripheral and central processes that shape the tonotopic representation of the auditory map. We selectively knocked out Isl1 in auditory neurons using Neurod1Cre strategies. In the absence of Isl1, spiral ganglion neurons migrate into the central cochlea and beyond, and the cochlear wiring is profoundly reduced and disrupted. The central axons of Isl1 mutants lose their topographic projections and segregation at the cochlear nucleus. Transcriptome analysis of spiral ganglion neurons shows that Isl1 regulates neurogenesis, axonogenesis, migration, neurotransmission-related machinery, and synaptic communication patterns. We show that peripheral disorganization in the cochlea affects the physiological properties of hearing in the midbrain and auditory behavior. Surprisingly, auditory processing features are preserved despite the significant hearing impairment, revealing central auditory pathway resilience and plasticity in Isl1 mutant mice. Mutant mice have a reduced acoustic startle reflex, altered prepulse inhibition, and characteristics of compensatory neural hyperactivity centrally. Our findings show that ISL1 is one of the obligatory factors required to sculpt auditory structural and functional tonotopic maps. Still, upon Isl1 deletion, the ensuing central compensatory plasticity of the auditory pathway does not suffice to overcome developmentally induced peripheral dysfunction of the cochlea. Snrpb, the gene mutated in CCMS, is haploinsufficient and required in neural crest cells for proper splicing of developmentally important genes essential for craniofacial morphogenesis Heterozygous mutations in SNRPB, an essential core component of the five small ribonucleoprotein particles of the spliceosome, are responsible for Cerebrocostomandibular Syndrome (CCMS). We show that Snrpb heterozygous embryos arrest shortly after implantation. Additionally, heterozygous deletion of Snrpb in the developing brain and neural crest cells models craniofacial malformations found in CCMS, and results in death shortly after birth. RNAseq analysis of mutant heads prior to morphological defects revealed increased exon-skipping and intron-retention in association with increased 5 splice site strength. We found increased exon-skipping in negative regulators of the P53-pathway, increased levels of nuclear P53, P53-target genes. However, removing TrP53 in Snrpb heterozygous mutant neural crest cells did not completely rescue craniofacial development. We also found a small but significant increase in exon-skipping of several transcripts required for head and midface development, including Smad2 and Rere. Furthermore, mutant embryos exhibited ectopic or missing expression of Fgf8 and Shh, which are required to coordinate face and brain development. Thus, we propose that mis-splicing of transcripts that regulate P53-activity and craniofacial-specific genes both contribute to craniofacial malformations. 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. Glucose inhibits haemostasis and accelerates diet-induced hyperlipidaemia in zebrafish larvae Hyperglycaemia damages the microvasculature in part through the reduced recruitment of immune cells and interference with platelet signalling, leading to poor wound healing and accelerated lipid deposition in mammals. We investigated the utility of zebrafish larvae to model the effect of glucose on neutrophil and macrophage recruitment to a tail wound, wound-induced haemostasis, and chicken egg yolk feed challenge-induced hyperlipidaemia by supplementing larvae with exogenous glucose by immersion or injection. Neither method of glucose supplementation affected the recruitment of neutrophils and macrophages following tail transection. Glucose injection reduced thrombocyte retention and fibrin plug formation while only thrombocyte retention was reduced by glucose immersion following tail transection. We observed accelerated lipid accumulation in glucose-injected larvae challenged with high fat chicken egg yolk feeding. Our study identifies conserved and divergent effects of high glucose on inflammation, haemostasis, and hyperlipidaemia in zebrafish larvae compared to mammals. AFFORDABLE, REPLICABLE ROBOTIC MINIRHIZOTRON SAMPLING FOR ROOT DYNAMICS STUDIES O_LIMinirhizotrons (paired camera systems and buried observatories) are the best current method to make repeatable measurements of fine root dynamics in the field. Automating data collection is the only way to assess fine-scale temporal changes in fine roots in natural systems. When paired with above-ground proximal remote sensing, automatic minirhizotron data could reveal coordination between above and below ground phenology, supporting both plant allocation and carbon exchange understanding. However, there are many challenges to automating minirhizotrons for field experiments. First, instruments must be cheap enough to replicate given shallow field of view and variable natural environments. Second, automated analysis must be robust to changeable soil and root conditions because image properties extracted must have biological meaning.
C_LIO_LIHere, we demonstrate an automatic minirhizotron scheme, built with off-the-shelf parts and capable to sample at sub-daily resolution paired with a neural network-based image analysis method. As a proof of concept, we tested the instruments and the methodology in a mesocosm study and in two short real-world field trials alongside datasets from other instruments (e.g. proximal remote sensed leaf dynamics).
C_LIO_LIWe show that the method produces robust daily time-series of root dynamics in the mesocosm and the field. We applied the same model across multiple instruments, demonstrating good reproducibility of measurements. Temporal variations of root cover were a strong driver of soil CO2 efflux in mesocosm, and vegetation dynamics from above ground proximal remote sensing and minirhizotron data differed both in mesocosm and field. We found sensitivity of root cover extracted to soil moisture and sometimes time of day (potentially relating to soil moisture), which may only affect high frequency imagery and not commonly reported as encountered when training neural networks on traditional, time-distinct minirhizotron studies.
C_LIO_LIWe discuss avenues for dealing with such issues in future field applications of minirhizotrons. If such issues are dealt with to a satisfactory manner in the field, timeseries of root biomass and traits from replicated instruments can add a new dimension to phenology studies at the ecosystem level by understanding the dynamics of root properties and traits.
C_LI LazySampling and LinearSampling: Fast Stochastic Sampling of RNA Secondary Structure with Applications to SARS-CoV-2 Many RNAs fold into multiple structures at equilibrium. The classical stochastic sampling algorithm can sample secondary structures according to their probabilities in the Boltzmann ensemble, and is widely used. However, this algorithm, consisting of a bottom-up partition function phase followed by a top-down sampling phase, suffers from three limitations: (a) the formulation and implementation of the sampling phase are unnecessarily complicated; (b) the sampling phase repeatedly recalculates many redundant recursions already done during the partition function phase; (c) the partition function runtime scales cubically with the sequence length. These issues prevent stochastic sampling from being used for very long RNAs such as the full genomes of SARS-CoV-2. To address these problems, we first adopt a hypergraph framework under which the sampling algorithm can be greatly simplified. We then present three sampling algorithms under this framework, among which the LazySampling algorithm is the fastest by eliminating redundant work in the sampling phase via on-demand caching. Based on LazySampling, we further replace the cubic-time partition function by a linear-time approximate one, and derive LinearSampling, an end-to-end linear-time sampling algorithm that is orders of magnitude faster than the standard one. For instance, LinearSampling is 176x faster (38.9s vs. 1.9h) than Vienna RNAsubopt on the full genome of Ebola virus (18,959 nt). More importantly, LinearSampling is the first RNA structure sampling algorithm to scale up to the full-genome of SARS-CoV-2 without local window constraints, taking only 69.2 seconds on its reference sequence (29,903 nt). The resulting sample correlates well with the experimentally-guided structures. On the SARS-CoV-2 genome, LinearSampling finds 23 regions of 15 nt with high accessibilities, which are potential targets for COVID-19 diagnostics and drug design.

See code: https://github.com/LinearFold/LinearSampling Rapid expansion of SARS-CoV-2 variants of concern is a result of adaptive epistasis The SARS-CoV-2 pandemic recently entered an alarming new phase with the emergence of the variants of concern (VOC) and understanding their biology is paramount to predicting future ones. Current efforts mainly focus on mutations in the spike glycoprotein (S), but changes in other regions of the viral proteome are likely key. We analyzed more than 900,000 SARS-CoV-2 genomes with a computational systems biology approach including a haplotype network and protein structural analyses to reveal lineage-defining mutations and their critical functional attributes. Our results indicate that increased transmission is promoted by epistasis, i.e., combinations of mutations in S and other viral proteins. Mutations in the non-S proteins involve immune-antagonism and replication performance, suggesting convergent evolution. Furthermore, adaptive mutations appear in geographically disparate locations, suggesting that either independent, repeat mutation events or recombination among different strains are generating VOC. We demonstrate that recombination is a stronger hypothesis, and may be accelerating the emergence of VOC by bringing together cooperative mutations. This emphasizes the importance of a global response to stop the COVID-19 pandemic. A nontoxic tetrodotoxin analog, 5,6,11-trideoxytetrodotoxin, as an olfactory chemoattractant for the grass puffer Toxic puffers contain the potent neurotoxin, tetrodotoxin (TTX). Although TTX is considered to serve as a defense substance, previous behavioral studies have demonstrated that TTX acts as an attractive pheromone for some toxic puffers. To elucidate the physiological mechanism of putative pheromonal action of TTX, we examined whether grass puffers Takifugu alboplumbeus can detect TTX. Electroolfactogram (EOG) results suggest that the olfactory epithelium (OE) of grass puffers responded to a type of TTX analog (5,6,11-trideoxyTTX), although it did not respond to TTX. We also examined the attractive action of 5,6,11-trideoxyTTX on grass puffers by recording their swimming behavior under dark conditions. Grass puffers preferred to stay on the side of the aquarium where 5,6,11-trideoxyTTX was administered, and their swimming speed decreased. Additionally, odorant-induced labeling of olfactory sensory neurons by immunohistochemistry against neural activity marker (phosphorylated extracellular signal regulated kinase; pERK) revealed that labeled olfactory sensory neurons were localized in the region surrounding "islets" where there was considered as nonsensory epithelium. 5,6,11-trideoxyTTX has been known to accumulate in grass puffers, but its toxicity is much lower (almost nontoxic) than TTX. Our results suggest that toxic puffers may positively use this unknown function TTX analog as an olfactory chemoattractant. Cell Cycle Delay, Pro-metaphase Arrest and C-metaphase Inducing Effects of Petroleum Ether Fraction of Leaf Aqueous Extract of Clerodendrum viscosum Vent. Clerodendrum viscosum is a traditionally used medicinal plant. The present study aimed to analyze cell cycle delay, pro-metaphase arrest, and c-metaphase inducing effects of the petroleum ether fraction (AQPEF) of leaf aqueous extract of C. Viscosum (LAECV). The LAECV was fractionated with petroleum ether and its metaphase arrest, cell cycle delay, and c-metaphase inducing activities were tested on A. cepa root tip cells. The AQPEF induced cell cycle delay, and colchicine like metaphase, c-metaphase, in A. cepa root tip cells. Thus, the present study explores AQPEF as an active fraction of LAECV having metaphase arresting, cell cycle delay, and c-metaphase inducing potentials. The BAF A12T mutation associated with premature aging impedes lamin A/C recruitment to sites of nuclear rupture, contributing to nuclear envelope fragility. The premature aging disorder Nestor Guillermo Progeria Syndrome (NGPS) is caused by a homozygous Alanine to Threonine mutation at position 12 (A12T) in Barrier-to- Autointegration Factor (BAF). BAF is a small essential protein that binds to DNA and nuclear envelope proteins. It contributes to important cellular processes including transcription regulation and nuclear envelope reformation after mitosis. More recently, BAF was identified as an important factor for nuclear envelope repair upon rupture in interphase. However, the mechanism by which the BAF A12T mutation causes NGPS has remained unclear. To investigate the effects of this mutation on nuclear envelope integrity, we used NGPS-derived patient cells and engineered an isogenic cell line by reversing the BAF A12T homozygous mutation using CRISPR/Cas9. Using a combination of cellular models, structural data and in vitro assays, we identified that the A12T mutation reduces the affinity of BAF for lamin A/C by tenfold. As a result, BAF A12T is unable to recruit lamin A/C to sites of nuclear envelope rupture. This leads to persistent lamin A/C gaps at sites of ruptures, and contributes to nuclear fragility in NGPS patient cells, which show increased frequency of nuclear envelope re- rupturing. Overexpression of wild-type BAF in a NGPS context rescues lamin A/C recruitment to sites of nuclear rupture, which could explain why the heterozygous A12T mutation does not cause premature aging. Changes to the mtDNA copy number during yeast culture growth We show that the mtDNA copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest ([~]10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains, but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype. FACT regulates pluripotency through distal regulation of gene expression in murine embryonic stem cells The FACT complex is a conserved histone chaperone with essential roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear. To determine this mechanism, we used inducible depletion of FACT subunit SPT16 in murine embryonic stem cells paired with genomic factor localization, nascent transcription, and chromatin accessibility analyses. Over a timecourse of SPT16 depletion, nucleosomes invade loci bound by master pluripotency factors and gene-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (OCT4), Sox2, Nanog, and enhancer RNAs produced at the genes' associated enhancers are downregulated, suggesting that FACT regulates expression of the pluripotency factors themselves. We find that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency. Modulations in motor unit discharge are related to changes in fascicle length during isometric contractions The integration of electromyography (EMG) and ultrasound imaging has provided important information about the mechanisms of muscle activation and contraction. Unfortunately, EMG does not allow an accurate assessment of the interplay between the neural drive received by muscles, changes in fascicle length (FL) and the force/torque produced. We aimed to assess the relationship between modulations in tibialis anterior (TA) motor unit (MU) firing rate, FL and dorsiflexion torque (DT) using ultrasound-transparent high-density EMG electrodes. EMG and ultrasound images were recorded simultaneously from TA, using a 32-electrode silicon matrix, while performing isometric dorsiflexion, at diverse ankle joint positions (0{degrees} and 30{degrees} plantar flexion) and torques (20% and 40% of maximum). EMG signals were decomposed into individual MUs and changes in FL were assessed with a fascicle-tracking algorithm. MU firings were converted into a cumulative spike train (CST) that was cross-correlated with DT (CST-DT) and FL (CST-FL). High cross-correlations were found for CST-FL, 0.60 (range: 0.31-0.85) and CST-DT 0.71 (range: 0.31-0.88). Cross-correlation lags revealed that the delay between CST-FL (~75ms) was significantly smaller than CST-DT (~150ms, p<0.001). These delays affected the interpretation of MU recruitment/de-recruitment thresholds, with FL showing similar lengths for both recruitment and de-recruitment. This study is the first to demonstrate that changes in TA FL are closely related to both modulations in MU firing frequency and DT. These relationships allow assessment of the interplay between neural drive, muscle contraction and resultant torque, thereby providing a better understanding of the mechanisms responsible for the generation of muscle force.

NEW AND NOTEWORTHYBy employing ultrasound-transparent high-density surface EMG electrodes, we show that modulations in tibialis anterior motor unit discharge rate were closely related to both changes in its fascicle length and resultant torque. These relationships allowed quantifying delays between neural drive and muscle shortening as well as muscle shortening and torque during submaximal isometric contractions, providing an accurate estimation of the time required to generate muscle force and subsequent production of torque via the tendon.