TMS with fast and accurate electronic control: measuring the orientation sensitivity of corticomotor pathways BackgroundTranscranial magnetic stimulation (TMS) coils allow only a slow, mechanical adjustment of the stimulating electric field (E-field) orientation in the cerebral tissue. Fast E-field control is needed to synchronize the stimulation with the ongoing brain activity. Also, empirical models that fully describe the relationship between evoked responses and the stimulus orientation and intensity are still missing.

ObjectiveWe aimed to (1) develop a TMS transducer for manipulating the E-field orientation electronically with high accuracy at the neuronally meaningful millisecond-level time scale and (2) devise and validate a physiologically based model describing the orientation selectivity of neuronal excitability.

MethodsWe designed and manufactured a two-coil TMS transducer. The coil windings were computed with a minimum-energy optimization procedure, and the transducer was controlled with our custom-made electronics. The electronic E-field control was verified with a TMS characterizer. The motor evoked potential amplitude and latency of a hand muscle were mapped in 3{degrees} steps of the stimulus orientation in 16 healthy subjects for three stimulation intensities. We fitted a logistic model to the motor response amplitude.

ResultsThe two-coil TMS transducer allows one to manipulate the pulse orientation accurately without manual coil movement. The motor response amplitude followed a logistic function of the stimulus orientation; this dependency was strongly affected by the stimulus intensity.

ConclusionThe developed electronic control of the E-field orientation allows exploring new stimulation paradigms and probing neuronal mechanisms. The presented model helps to disentangle the neuronal mechanisms of brain function and guide future non-invasive stimulation protocols.
Dynamic and stationary brain connectivity during movie watching as revealed by functional MRI Spatially remote brain regions show synchronized activity as typically revealed by correlated functional MRI (fMRI) signals. An emerging line of research has focused on the temporal fluctuations of connectivity; however, its relationships with stationary connectivity have not been clearly illustrated. We examined dynamic and stationary connectivity when the participants watched four different movie clips. We calculated point-by-point multiplication between two regional time series to estimate the time-resolved dynamic connectivity, and estimated the inter-individual consistency of the dynamic connectivity time series. Widespread consistent dynamic connectivity was observed for each movie clip, which also showed differences between the clips. For example, a cartoon movie clip, Wall-E, showed more consistent of dynamic connectivity with the posterior cingulate cortex and supramarginal gyrus, while a court drama clip, A Few Good Men, showed more consistent of dynamic connectivity with the auditory cortex and temporoparietal junction, which might suggest the involvement of specific brain processing for different movie contents. In contrast, the stationary connectivity as measured by the correlations between regional time series was highly similar among the movie clips, and showed fewer statistically significant differences. The patterns of consistent dynamic connectivity could be used to classify different movie clips with higher accuracy than the stationary connectivity and regional activity. These results support the functional significance of dynamic connectivity in reflecting functional brain changes, which could provide more functionally relevant information than stationary connectivity.
Outer hair cell electro-mechanical power conversion Outer hair cells are the cellular motors in the mammalian inner ear responsible for sensitive high-frequency hearing. Motor function requires expression of the protein prestin in the OHC lateral membrane and mechano-electrical transduction in the apical hair bundle. In the present report, electrical power consumption and mechanical power output are determined using previously published voltage clamp data from isolated OHCs and membrane patches. Results reveal that power output by the prestin-motor complex peaks at a best frequency much higher than implied by the low-pass character of nonlinear capacitance, and much higher than the whole-cell resistive-capacitive corner frequency. High frequency power output is enabled by a -90{degrees} shift in the phase of prestin-dependent electrical charge displacement, manifested electrically as emergence of imaginary-valued nonlinear capacitance. First detection of herpesvirus and prevalence of mycoplasma infection in free-ranging Hermann s tortoises (Testudo hermanni), and in potential pet vectors AO_SCPLOWBSTRACTC_SCPLOWTwo types of pathogens cause highly contagious upper respiratory tract diseases (URTD) in Chelonians: testudinid herpesviruses (TeHV) and a mycoplasma (Mycoplasma agassizii). In captivity, these infections are frequent and can provoke outbreaks. Pet trade generates international flow of tortoises, often without sanitary checking; individuals intentionally or accidentally released in the wild may spread pathogens. A better understanding of the transmission of infectious agents from captivity to wild tortoises is needed. Many exotic individuals have been introduced in populations of the endangered western Hermanns tortoise (Testudo hermanni hermanni), notably spur-thighed tortoises (Testudo graeca). We assessed the presence of TeHV and mycoplasma in native western Hermanns tortoises and in potential pet vectors in south-eastern France. Using a large sample (N=572 tortoises), this study revealed, by PCR, the worrying presence of herpesvirus in 7 free-ranging individuals (3 sub-populations). Additionally, Mycoplasma agassizii was detected, by PCR, in 15 of the 18 populations sampled with a frequency ranging from 3.4% (1 of 29 tortoises) to 25% (3 of 12 tortoises). Exotic spur-thighed tortoises showed high frequency of mycoplasma infection in captivity (18.2%) and in individuals (50%) found in native Hermanns tortoise sub-populations, suggesting that this species could be a significant vector. The paucity of information of TeHV on European tortoise URTD in natural settings, especially in combination with mycoplasma, prompts for further studies. Indeed, sick tortoises remain concealed and may not be easily detected in the field. Our results indicate a good health for most infected tortoise but it should be screened in the field as well as in captivity. Modeling miRNA-driven post-transcriptional regulatory signals by using exon-intron split analysis (EISA) in pigs The contribution of microRNAs (miRNAs) to mRNA regulation has often been explored by post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes of mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Gottingen minipigs. Exon-intron split analysis (EISA) of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of upregulated miRNAs also in the fed state. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, EISA provides a valuable view about the regulation of gene expression, complementary to canonical differential expression analyses. Distributions of extinction times from fossil ages and tree topologies: the example of mid-Permian synapsid extinctions Given a phylogenetic tree that includes only extinct, or a mix of extinct and extant taxa, where at least some fossil data are available, we present a method to compute the distribution of the extinction time of a given set of taxa under the Fossilized-Birth-Death model.

Our approach differs from the previous ones in that it takes into account (i) the possibility that the taxa or the clade considered may diversify before going extinct and (ii) the whole phylogenetic tree to estimate extinction times, whilst previous methods do not consider the diversification process and deal with each branch independently. Because of this, our method can estimate extinction times of lineages represented by a single fossil, provided that they belong to a clade that includes other fossil occurrences.

We assess and compare our new approach with a standard previous one using simulated data. Results show that our method provides more accurate confidence intervals.

This new approach is applied to the study of the extinction time of three Permo-Carboniferous synapsid taxa (Ophiacodontidae, Edaphosauridae, and Sphenacodontidae) that are thought to have disappeared toward the end of the Cisuralian (early Permian), or possibly shortly thereafter. The timing of extinctions of these three taxa and of their component lineages supports the idea that the biological crisis in the late Kungurian/early Roadian consisted of a progressive decline in biodiversity throughout the Kungurian. 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. Co-infection of Candidatus Piscichlamydia trichopodus (order Chlamydiales) and Henneguya sp (Myxosporea, Myxobolidae) in snakeskin gourami Trichopodus pectoralis (Regan 1910) The present study describes a simultaneous infection of a novel Chlamydia-like organism (CLO) with a Myxozoa parasite, Henneguya sp. in snakeskin gourami Trichopodus pectoralis in Thailand. A new CLO is proposed "Candidatus Piscichlamydia trichopodus" (CPT) based on 16S rRNA phylogenetic analysis. Systemic intracellular CPT infection was confirmed by histological examination, in situ hybridization, PCR assay, and sequencing of 16S rRNA. This novel pathogen belongs to the order Chlamydiales but differs in certain aspects from other species. The histopathological changes associated with CPT infection were different from the typical pathological lesions of epitheliocystis caused by previously known CLO. Unlike other CLO, CPT localized in the connective tissue rather than in the epithelial cells and formed smaller clumps of intracellular bacteria that stained dark blue with hematoxylin. On the other hand, typical myxospores of the genus Henneguya with tails were observed in the gill sections. Infection with Henneguya sp. resulted in extensive destruction of the gill filaments, most likely leading to respiratory distress. Due to the frequency of co-infections and the unavailability of culture methods for CLO and Henneguya sp, it was difficult to determine which pathogens were directly responsible for the associated mortality. However, co-infections may increase the negative impact on the host and the severity of the disease. Given the commercial importance of the snakeskin gourami and its significant aquaculture potential, the findings of this study are important for further studies on disease prevention. Contribution of arm movements to balance recovery after tripping in older adults. Falls are common in daily life, often caused by trips and slips and, particularly in older adults, with serious consequences. Although arm movements play an important role in balance control, there is limited research into the role of arm movements during balance recovery after tripping in older adults. We investigated how older adults use their arms to recover from a trip and the difference in the effects of arm movements between fallers (n=5) and non-fallers (n=11).

Sixteen older males and females (69.7{+/-}2.3 years) walked along a walkway and were occasionally tripped over suddenly appearing obstacles. We analysed the first trip using a biomechanical model based on full-body kinematics and force-plate data to calculate whole body orientation during the trip and recovery phase. With this model, we simulated the effects of arm movements at foot-obstacle impact and during trip recovery on body orientation.

Apart from an increase in sagittal plane forward body rotation at touchdown in fallers, we found no significant differences between fallers and non-fallers in the effects of arm movements on trip recovery. Like earlier studies in young adults, we found that arm movements during the recovery phase had most favourable effects in the transverse plane: by delaying the transfer of angular momentum of the arms to the body, older adults rotated the tripped side more forward thereby allowing for a larger recovery step. Older adults that are prone to falling might improve their balance recovery after tripping by learning to prolong ongoing arm movements. 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. Multivariate genome-wide association study on tissue-sensitive diffusion metrics identifies key molecular pathways for axonal growth, synaptogenesis, and astrocyte-mediated neuroinflammation The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discovered 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% were validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development were found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders. Innovative GenExpA software for selecting suitable reference genes for reliable normalization of gene expression in melanoma The algorithms commonly used to select the best stable reference gene in RT-qPCR data analysis have their limitations. We showed that simple selection of the reference gene or pair of genes with the lowest stability value from the pool of potential reference genes - a commonly used approach - is not sufficient to accurately and reliably normalize the target gene transcript and can lead to biologically incorrect conclusions. For reliable assessment of changes in a target gene expression level, we propose our innovative GenExpA software, which works in a manner independent of the experimental model and the normalizer used. GenExpA software selects the best reference by combining the NormFinder algorithm with progressive removal of the least stable gene from the candidate genes in a given experimental model and in the set of daughter models assigned to it. The reliability of references is validated based on the consistency of the statistical analyses of normalized target gene expression levels through all models, described by the coherence score (CS). The use of the CS value imparts a new quality to qPCR analysis because it clarifies how low the stability value of reference must be in order for biologically correct conclusions to be drawn. We tested our method on qPCR data for the B4GALT genes family in melanoma, which is characterized by a high mutation rate, and in melanocytes. GenExpA is available at https://github.com/DorotaHojaLukowicz/GenExpA or https://www.sciencemarket.pl/baza-programow-open-source#oferty.

HighlightsGenExpA - next-generation software for normalizer selection and validation

The GenExpA tool defines how low the stability value of the reference should be

The GenExpA tool determines the level of gene expression analysis robustness

The GenExpA tool increases the speed of analysis and reduces cost A General Biphasic-Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate and Drug Clearance from Birth to Adulthood Understanding the maturation process of human physiology and metabolism has broad medical and pharmaceutical implications. Age and bodyweight are frequently considered as separate variables in modeling the dynamical changes of human organ functions and of drug clearance from birth to adulthood. The objective of this study is to propose a unified, continuous and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR) and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from age-incorporated maturation-and-growth models with reported parameter values. The model used in the simulation is independent of the proposed model. A statistical approach was used to simulate the model generated CL and GFR data for a hypothetical population with 26 age groups (ranging from 0 to 20 years). Besides, individual CL data for one drug, and sparse PBPK-modeled CL values for two drugs were also included for analysis. A 4-parameter, mixed-allometry equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. Compared with simple allometry, the biphasic model fits satisfactorily to all datasets (based on Akaikes Information Criterion and residual plots). The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (slope >1; ca. 1.5 - 3.5) below the characteristic bodyweight at the phase transition; and above which, a sublinear scaling phase (slope <1; ca. 0.5 - 0.7) is evident. The phase-transition bodyweight is ranging from 5 to 20 kg (corresponding to 0.5 - 9 years) and the mean value is around 10 kg ([~]2 years) for all data sets. The dimensionless analysis generalizes, and offers quantitative realization of, the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition occurring in the human maturation process of BMR, GFR and drug CL. Single-dose AAV-based vaccine induces a high level of neutralizing antibodies against SARS-CoV-2 in rhesus macaques Coronavirus disease 2019 (COVID-19), which is triggered by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, continues to threaten global public health. Developing a vaccine that only requires single immunization but provides long-term protection for the prevention and control of COVID-19 is important. Here, we developed an adeno-associated virus (AAV)-based vaccine expressing a stable receptor-binding domain (SRBD) protein. The vaccine requires only a single shot but provides effective neutralizing antibodies (NAbs) over 598 days in rhesus macaques (Macaca mulatta). Importantly, our results showed that the NAbs were kept in high level and long lasting against authentic wild-type SARS-CoV-2, Beta, Delta and Omicron variants using plaque reduction neutralization test. Of note, although we detected pre-existing AAV2/9 antibodies before immunization, the vaccine still induced high and effective NAbs against COVID-19 in rhesus macaques. AAV-SRBD immune serum also efficiently inhibited the binding of ACE2 with RBD in the SARS-CoV-2 B.1.1.7 (Alpha), B.1.351 (Beta), P.1/P.2 (Gamma), B.1.617.2 (Delta), B.1.617.1/3(Kappa), and C.37 (Lambda) variants. Thus, these data suggest that the vaccine has great potential to prevent the spread of SARS-CoV-2. Suppression of sost/sclerostin and dickkopf-1 promote intervertebral disc structure in mice Intervertebral disc (IVD) degeneration is a leading cause of low back pain and characterized by accelerated extracellular matrix breakdown and IVD height loss but there is no approved pharmacological therapeutic. Deletion of Wnt signaling receptor Lrp5 induces IVD degeneration and suggests that Wnt signaling in the IVD may be responsive to inhibition of Wnt signaling inhibitors sost(gene)/sclerostin(protein) or dickkopf-1 (dkk1). Anti-sclerostin antibody (Scl-Ab) is an FDA-approved bone therapeutic that activates Wnt signaling. We (1) determined if pharmacological neutralization of sclerostin, dkk1 or their combination stimulate Wnt signaling and promote IVD structure and (2) determined the extent of the response of the IVD to global, persistent deletion of sost. Nine-week-old C57Bl/6J female mice (n=6-7/grp) were subcutaneously injected 2x/wk for 5.5 wk with scl-Ab (25 mg/kg), dkk1-Ab (25 mg/kg), 3:1 scl-Ab/dkk1-Ab (18.75:6.25 mg/kg) or vehicle (Veh). Separately, IVD of sost KO and WT (wildtype) mice (n=8, grp) were harvested at 16 weeks of age. First, compared to vehicle, scl-Ab, dkk1-Ab and 3:1 scl-Ab/dkk1-Ab similarly increased lumbar IVD height and {beta}-catenin gene expression. Despite these similarities, scl-Ab decreased cellular stress-related heat shock protein gene expressions while neither dkk1-Ab nor scl-Ab/dkk1-Ab altered the same. Genetically and compared to WT, sost KO increased MRI-determined hydration and proteoglycan staining in the IVD. Notably, persistent deletion of sost was compensated by upregulation of dkk1, which consequently reduced the cell nuclear fraction for Wnt signaling transcription factor {beta}-catenin in whole IVD. Lastly, RNA-sequencing pathway analysis confirmed the compensatory suppression of Wnt signaling and determined a reduction of cellular stress pathways. Together, suppression of sost/sclerostin or dkk1 each promote IVD structure by stimulating Wnt signaling, but sclerostin and dkk1 may differentially regulate cellular stress pathways. Ultimately, postmenopausal women prescribed scl-Ab injections to prevent vertebral fracture may also benefit from a restoration of IVD height and health.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=72 SRC="FIGDIR/small/449486v2_ufig1.gif" ALT="Figure 1">
View larger version (20K):
org.highwire.dtl.DTLVardef@db2e17org.highwire.dtl.DTLVardef@1af1bbborg.highwire.dtl.DTLVardef@145cc4corg.highwire.dtl.DTLVardef@167a19e_HPS_FORMAT_FIGEXP  M_FIG C_FIG Suppression of Wnt signaling inhibitors by genetic or pharmacological approaches promotes intervertebral disc structure and hydration by Wnt signaling. However, persistent activation of Wnt signaling induces a compensatory reduction of Wnt signaling that shifts IVD cells toward a chondrocyte-like (CLC) phenotype. AF: annulus fibrosus, NC: notochordal cell, NP: nucleus pulposus, PG: proteoglycan Hypoxia reduces cell attachment of SARS-CoV-2 spike protein by modulating the expression of ACE2, neuropilin-1, syndecan-1 and cellular heparan sulfate A main clinical parameter of COVID-19 pathophysiology is hypoxia. Here we show that hypoxia decreases the attachment of the receptor binding domain (RBD) and the S1 subunit (S1) of the spike protein of SARS-CoV-2 to epithelial cells. In Vero E6 cells, hypoxia reduces the protein levels of ACE2 and neuropilin-1 (NRP1), which might in part explain the observed reduction of the infection rate. In addition, hypoxia inhibits the binding of the spike to NCI-H460 human lung epithelial cells by decreasing the cell surface levels of heparan sulfate (HS), a known attachment receptor of SARS-CoV-2. This interaction is also reduced by lactoferrin, a glycoprotein that blocks HS moieties on the cell surface. The expression of syndecan-1, an HS-containing proteoglycan expressed in lung, is inhibited by hypoxia on a HIF-1-dependent manner. Hypoxia or deletion of syndecan-1 results in reduced binding of the RBD to host cells. Our study indicates that hypoxia acts to prevent SARS-CoV-2 infection, suggesting that the hypoxia signaling pathway might offer therapeutic opportunities for the treatment of COVID-19. Macroscopic control of synchronous electrical signaling with chemically-excited gene expression Excitable cells can convert electrical signals into chemical outputs to facilitate the active transport of information across larger distances. This electrical-to-chemical conversion requires a tightly regulated expression of ion channels. Alterations of ion channel expression provide landmarks of numerous pathological diseases, such as cardiac arrhythmia, epilepsy, or cancer. Although the activity of ion channels can be locally regulated by external light or chemical stimulus, it remains challenging to coordinate the expression of ion channels on extended spatial-temporal scales in a non-invasive manner. Here, we have engineered yeast S. cerevisiae to read and convert local chemical concentrations into a dynamic electrical field distributed across cell populations. The core mechanism encodes a chemically-excitable dual-feedback gene circuit that precisely tunes the expression domain of potassium channels, globally coordinating cyclic firing of the plasma membrane potential (PMP). We demonstrate that this mechanism leverages an engineered constitutively open bacterial potassium channel KcsA to directly couple chemical stimuli with ion flux through gene expression and it can interface with the host ion channels through the pulsatile production of toxins. Our study provides a robust synthetic transcriptional toolbox underlying the conversion of local chemical environments into spatiotemporally organized electrical impulses for various cellular engineering, synthetic biology, and potential therapeutic applications. Deep learning-based quantification of arbuscular mycorrhizal fungi in plant roots Soil fungi establish mutualistic interactions with the roots of most vascular land plants. Arbuscular mycorrhizal (AM) fungi are among the most extensively characterised mycobionts to date. Current approaches to quantifying the extent of root colonisation and the abundance of hyphal structures in mutant roots rely on staining and human scoring involving simple, yet repetitive tasks prone to variations between experimenters.

We developed AMFinder which allows for automatic computer vision-based identification and quantification of AM fungal colonisation and intraradical hyphal structures on ink-stained root images using convolutional neural networks.

AMFinder delivered high-confidence predictions on image datasets of roots of multiple plant hosts (Nicotiana benthamiana, Medicago truncatula, Lotus japonicus, Oryza sativa) and captured the altered colonisation in ram1-1, str and smax1 mutants. A streamlined protocol for sample preparation and imaging allowed us to quantify mycobionts from the genera Rhizophagus, Claroideoglomus, Rhizoglomus and Funneliformis via flatbed scanning or digital microscopy including dynamic increases in colonisation in whole root systems over time.

AMFinder adapts to a wide array of experimental conditions. It enables accurate, reproducible analyses of plant root systems and will support better documentation of AM fungal colonisation analyses. AMFinder can be accessed here: https://github.com/SchornacklabSLCU/amfinder.git 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. Predicting the re-distribution of antibiotic molecules caused byinter-species interactions in microbial communities Microbes associate in nature forming complex communities, but they are often studied in purified form. Here I show that differences in antibiotic sensitivity between co-existing bacterial species can modify the diffusion of antibiotic molecules, and alter the drug efficacy measured in conventional in vitro assays. I developed a model, validated experimentally using an artificial two-species community, showing antibiotic molecules re-distribute evenly across all species if they have similar drug sensitivity. The result is lower drug efficacy with respect to conventional in vitro assays. However, when species have different sensitivities, the re-distribution of antibiotic is uneven and most sensitive species accumulate more drug. The result is increased drug efficacy against these species. Despite this discrepancy with conventional laboratory assays, my data suggests the changes in drug efficacy are predictable and provides insight into how communities can affect the diffusion of other molecules in the environment. Natural variants in SARS-CoV-2 S protein pinpoint structural and functional hotspots; implications for prophylaxis and therapeutics strategies In December 2019, a novel coronavirus, termed severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was identified as the cause of pneumonia with severe respiratory distress and outbreaks in Wuhan, China. The rapid and global spread of SARS-CoV-2 resulted in the coronavirus 2019 (COVID-19) pandemic. Earlier during the pandemic, there were limited genetic viral variations. As millions of people became infected, multiple single amino acid substitutions emerged. Many of these substitutions have no consequences. However, some of the new variants show a greater infection rate, more severe disease, and reduced sensitivity to current prophylaxes and treatments. Of particular importance in SARS-CoV-2 transmission are mutations that occur in the Spike (S) protein, the protein on the viral outer envelope that binds to the human angiotensin-converting enzyme receptor (hACE2). Here, we conducted a comprehensive analysis of 441,168 individual virus sequences isolated from humans throughout the world. From the individual sequences, we identified 3,540 unique amino acid substitutions in the S protein. Analysis of these different variants in the S protein pinpointed important functional and structural sites in the protein. This information may guide the development of effective vaccines and therapeutics to help arrest the spread of the COVID-19 pandemic. Affordable Caenorhabditis elegans tracking system for classroom use For decades the nematode C. elegans has served as an outstanding research organism owing to its unsurpassed experimental amenability. This advantage has also made this tiny worm an attractive vehicle for science instruction across higher learning institutions. However, the prohibitive cost associated with the automated behavioral assessment of these animals remains an obstacle preventing their full adoption in undergraduate and high school settings. To improve this situation, we developed an inexpensive worm tracking system for use by high school interns and undergraduate students. Over the past two years this tracker has been successfully used by undergraduate students in our introductory Cell and Molecular lab (BSC220) at Illinois State University. Here we describe and demonstrate the use of our inexpensive worm tracking system. GPI-anchored FGF directs cytoneme-mediated bidirectional signaling to self-regulate tissue-specific dispersion How signaling proteins generate a multitude of information to organize tissue patterns is critical to understanding morphogenesis. In Drosophila, FGF produced in wing-disc cells regulates the development of the disc-associated air-sac-primordium (ASP). Here, we show that FGF is Glycosylphosphatidylinositol-anchored to the producing cell surface and that this modification both inhibits free FGF secretion and activates target-specific bidirectional FGF-FGFR signaling through cytonemes. FGF-source and recipient ASP cells extend cytonemes that present FGF and FGFR on their surfaces and reciprocally recognize each other over distance by contacting through CAM-like FGF-FGFR binding. Contact-mediated FGF-FGFR binding induces bidirectional signaling, which, in turn, promotes ASP and source cells to polarize cytonemes toward each other and reinforce signaling contacts. Subsequent un-anchoring of FGFR-bound-FGF from the source cell membrane dissociates cytoneme contacts and delivers FGF target-specifically to ASP cytonemes for paracrine functions. Thus, GPI-anchored FGF organizes both source and recipient cells and self-regulates its cytoneme-mediated tissue-specific dispersion and signaling. Modelling porcine reproductive and respiratory syndrome virus dissemination dynamics to quantifythe contribution of multiple modes of transmission: between-farm pig and vehicle movements, farmto-farm proximity, feed ingredients, and re-break Accounting for multiple modes of livestock disease dissemination in epidemiological models remains a challenge. We developed and calibrated a mathematical model for transmission of porcine reproductive and respiratory syndrome virus (PRRSV), tailored to fit nine modes of between-farm transmission pathways including: farm-to-farm proximity (local transmission), contact network of batches of pigs transferred between farms (pig movements), re-break probabilities for farms with previous PRRSV outbreaks, with the addition of four different contact networks of transportation vehicles (vehicles to transport pigs to farms, pigs to markets, feed and crew) and the amount of animal by-products within feed ingredients (e.g. animal fat or meat and bone meal). The model was calibrated on weekly PRRSV outbreaks data. We assessed the role of each transmission pathway considering the dynamics of specific types of production (i.e., sow, nursery). Although our results estimated that the networks formed by transportation vehicles were more densely connected than the network of pigs transported between-farms, pig movements and farm proximity were the main PRRSV transmission routes regardless of farm types. Among the four vehicle networks, vehicles transporting pigs to farms explained a large proportion of infections, sow = 20.9%; nursery = 15%; and finisher = 20.6%. The animal by-products showed a limited association with PRRSV outbreaks through descriptive analysis, and our model results showed that the contribution of animal fat contributed only 2.5% and meat and bone meal only 0.03% of the infected sow farms. Our work demonstrated the contribution of multiple routes of PRRSV dissemination, which has not been deeply explored before. It also provides strong evidence to support the need for cautious, measured PRRSV control strategies for transportation vehicles and further research for feed by-products modeling. Finally, this study provides valuable information and opportunities for the swine industry to focus effort on the most relevant modes of PRRSV between-farm transmission. Directed evolution of angiotensin-converting enzyme 2 (ACE2) peptidase activity profiles for therapeutic applications Angiotensin Converting-Enzyme 2 (ACE2) is currently being investigated for its ability to beneficially modulate the Angiotensin receptor (ATR) therapeutic axis to treat multiple human diseases, but its broad substrate scope and diverse physiological roles limit its potential as a therapeutic agent. In this work we engineer ACE2 variants with enhanced hydrolytic activity and specificity toward Angiotensin-II (Ang-II) that may lead to new therapeutic candidates with improved efficacy and reduced off-target side effects. We established a yeast display-based liquid chromatography screen that enabled use of directed evolution to discover ACE2 variants with improved Ang-II activity and specificity relative to the off-target peptide substrates Apelin-13 and Angiotensin-I (Ang-I). We screened ACE2 active site libraries to reveal three substitution-tolerant positions (M360, T371 and Y510) that can be mutated to enhance ACE2s activity profile and followed up on these hits with focused double mutant libraries to further improve the enzyme. Relative to wildtype ACE2, our top variant (T371L/Y510Ile) displayed a sevenfold increase in Ang-II turnover number (kcat), a sixfold diminished catalytic efficiency (kcat/Km) on Apelin-13, and an overall decreased activity on other ACE2 substrates that were not directly assayed in directed evolution screen. At physiologically relevant substrate concentrations, T371L/Y510Ile hydrolyzes more Ang-II than wildtype ACE2 with concomitant Ang-II:Apelin-13 specificity improvements reaching 30-fold. Our efforts have delivered ATR axis-acting therapeutic candidates with relevance to both established and unexplored ACE2 therapeutic applications and provide a foundation for further ACE2 engineering efforts. The immunopeptidome from a genomic perspective: Establishing immune-relevant regions for cancer vaccine design A longstanding disconnect between the growing number of MHC Class I immunopeptidomic studies and genomic medicine hinders cancer vaccine design. We develop COD-dipp to genomically map the full spectrum of detected canonical and non-canonical (non-exonic) MHC Class I antigens from 26 cancer studies. We demonstrate that patient mutations in regions overlapping physically identified antigens better predict immunotherapy response when compared to neoantigen predictions. We suggest a vaccine design approach using 140,966 highly immune-visible regions of the genome annotated by their expression and haplotype frequency in the human population. These regions tend to be highly conserved, mutated in cancer and harbor 7.8 times more immunogenicity. Intersecting pan-cancer mutations with these immune surveilled regions revealed a potential to create off-the-shelf multi-epitope vaccines against public neoantigens. Here we release COD-dipp, a cancer vaccine toolkit as a web-application (https://www.proteogenomics.ca/COD-dipp) and open-source high-throughput resource.