title: Analyses of symbiotic bacterial communities in the plant pest Bemisia tabaci reveal high prevalence of Hemipteriphilus asiaticus on the African continent

abstract: AO_SCPLOWBSTRACTC_SCPLOWMicrobial symbionts are widespread in insects and some of them have been associated to adaptive changes. Primary symbionts (P-symbionts) have a nutritional role that allows their hosts to feed on unbalanced diets (plant sap, wood, blood). Most of them have undergone genome reduction, but their genomes still retain genes involved in pathways that are necessary to synthesize the nutrients that their hosts need. However, in some P-symbionts, essential pathways are incomplete and secondary symbionts (S-symbionts) are required to complete parts of their degenerated functions. The P-symbiont of the phloem sap-feeder Bemisia tabaci, Candidatus Portiera aleyrodidarium, lacks genes involved in the synthesis of vitamins, cofactors, and also of some essential amino-acids. Seven S-symbionts have been detected in the B. tabaci species complex. Phenotypic and genomic analyses have revealed various effects, from reproductive manipulation to fitness benefits, notably some of them have complementary metabolic capabilities to Candidatus Portiera aleyrodidarium, suggesting that their presence may be obligatory. In order to get the full picture of the symbiotic community of this pest, we investigated, through metabarcoding approaches, the symbiont content of individuals from Burkina Faso, a West African country where B. tabaci induces severe crop damage. While no new putative B. tabaci S-symbiont was identified, Candidatus Hemipteriphilus asiaticus, a symbiont only described in B. tabaci populations from Asia, was detected for the first time on this continent. Phylogenetic analyses however reveal that it is a different strain than the reference found in Asia. Specific diagnostic PCRs showed a high prevalence of these S-symbionts and especially of Candidatus Hemipteriphilus asiaticus in different genetic groups. These results suggest that Candidatus Hemipteriphilus asiaticus may affect the biology of B. tabaci and provide fitness advantage in some B. tabaci populations.


title: Analyses of symbiotic bacterial communities in the plant pest Bemisia tabaci reveal high prevalence of Candidatus Hemipteriphilus asiaticus on the African continent

category: evolutionary biology


title: An integrative omics approach reveals posttranscriptional mechanisms underlying circadian temperature compensation

abstract:
A defining property of circadian clocks is temperature compensation, characterized by the resilience of circadian free-running periods against changes in environmental temperature. As an underlying mechanism, the balance or critical reaction hypothesis have been proposed. While the former supposes a temperature-dependent balancing of reactions with opposite effects on circadian period, the latter assumes an insensitivity of certain critical period determining regulations upon temperature changes. Posttranscriptional regulations such as temperature-sensitive alternative splicing or phosphorylation have been described as underlying reactions.

Here, we show that knockdown of cleavage and polyadenylation specificity factor subunit 6 (CPSF6), a key regulator of 3-end cleavage and polyadenylation, abolishes circadian temperature compensation in U-2 OS cells. We apply a combination of 3-End-RNA-seq and mass spectrometry-based proteomics to globally quantify changes in 3 UTR length as well as gene and protein expression between wild type and CPSF6 knock-down cells and their dependency on temperature. Analyzing differential responses upon temperature changes in wild type and CPSF6 knockdown cells reveals candidate genes underlying circadian temperature compensation. We identify that eukaryotic translation initiation factor 2 subunit 1 (EIF2S1) is among these candidates. EIF2S1 is known as a master regulator of cellular stress responses that additionally regulates circadian rhythms. We show that knockdown of EIF2S1 furthermore impairs temperature compensation, suggesting that the role of CPSF6 in temperature compensation may be mediated by its regulation of EIF2S1.

category: genomics


title: Application of machine learning methods in clinical trials for precision medicine

abstract:
ObjectiveA key component for precision medicine is a good prediction algorithm for patients response to treatments. We aim to implement machine learning (ML) algorithms into the response-adaptive randomization (RAR) design and improve the treatment outcomes.

Materials and MethodsWe incorporated nine ML algorithms to model the relationship of patient responses and biomarkers in clinical trial design. Such a model predicted the response rate of each treatment for each new patient and provide guidance for treatment assignment. Realizing that no single method may fit all trials well, we also built an ensemble of these nine methods. We evaluated their performance through quantifying the benefits for trial participants, such as the overall response rate and the percentage of patients who receive their optimal treatments.

ResultsSimulation studies showed that the adoption of ML methods resulted in more personalized optimal treatment assignments and higher overall response rates among trial participants. Compared with each individual ML method, the ensemble approach achieved the highest response rate and assigned the largest percentage of patients to their optimal treatments. For the real-world study, we successfully showed the potential improvements if the proposed design had been implemented in the study.

ConclusionIn summary, the ML-based RAR design is a promising approach for assigning more patients to their personalized effective treatments, which makes the clinical trial more ethical and appealing. These features are especially desirable for late-stage cancer patients who have failed all the FDA-approved treatment options and only can get new treatments through clinical trials.

category: bioinformatics


title: Cell Stiffening Contributes to Complement-mediated Injury of Choroidal Endothelial Cells in Early AMD

abstract: Age-related macular degeneration (AMD) is the leading cause of blindness in the aging population. Yet, no therapies exist for [~]85% of all AMD patients who have the dry form that is marked by degeneration of the retinal pigment epithelium (RPE) and underlying choroidal vasculature. As the choroidal vessels are crucial for RPE development and maintenance, understanding how they degenerate may lead to effective therapies for dry AMD. One likely causative factor for choroidal vascular loss is the cytolytic membrane attack complex (MAC) of the complement pathway that is abundant on choroidal vessels of humans with early dry AMD. To examine this possibility, we studied the effect of complement activation on choroidal endothelial cells (ECs) isolated from a rhesus monkey model of early AMD that, we report, exhibits MAC deposition and choriocapillaris endothelial loss similar to that seen in human early AMD. Treatment of choroidal ECs from AMD eyes with complement-competent normal human serum caused extensive actin cytoskeletal injury that was significantly less pronounced in choroidal ECs from young normal monkey eyes. We further show that ECs from AMD eyes are significantly stiffer than their younger counterparts and exhibit peripheral actin organization that is distinct from the longitudinal stress fibers in young ECs. Finally, these differences in complement susceptibility and mechanostructural properties were found to be regulated by the differential activity of small GTPases Rac and Rho because Rac inhibition in AMD cells led to simultaneous reduction in stiffness and complement susceptibility while Rho inhibition in young cells exacerbated complement injury. Thus, by identifying cell stiffness and cytoskeletal regulators Rac and Rho as important determinants of complement susceptibility, the current findings offer a new mechanistic insight into choroidal vascular loss in early AMD that warrants further investigation for assessment of translational potential.

category: pathology


title: Aberrant Functional Connectivity between Reward and Inhibitory Control Networks in Pre-Adolescent Binge Eating Disorder

abstract:
BackgroundBehavioral features of binge eating disorder (BED) suggest abnormalities in reward and inhibitory control. Studies of adult populations suggest functional abnormalities in reward and inhibitory control networks. Despite behavioral markers often developing in children, the neurobiology of pediatric BED remains unstudied.

Methods58 pre-adolescent children (aged 9-10-years) with BED and 66 age, BMI and developmentally-matched control children were extracted from the 3.0 baseline (Year 0) release of the Adolescent Brain Cognitive Development (ABCD) Study. We investigated group differences in resting-state functional MRI (rs-fMRI) functional connectivity (FC) within and between reward and inhibitory control networks. A seed-based approach was employed to assess nodes in the reward (orbitofrontal cortex, nucleus accumbens, amygdala) and inhibitory control (dorsolateral prefrontal cortex, anterior cingulate cortex) networks via hypothesis-driven seed-to- seed analyses, and secondary seed-to-voxel analyses.

ResultsOur findings revealed reduced FC between the dlPFC and amygdala, and between the anterior cingulate cortex and orbitofrontal cortex in pre-adolescent children with BED, relative to age, gender, BMI and developmentally matched controls. These findings indicating aberrant connectivity between nodes of inhibitory control and reward networks were corroborated by the whole-brain FC analyses.

ConclusionsEarly-onset BED may be characterized by diffuse abnormalities in the functional synergy between reward and cognitive control networks, without perturbations within reward and inhibitory control networks, respectively. The decreased capacity to regulate a reward-driven pursuit of hedonic foods, which is characteristic of BED, may in part, rest on this dysconnectivity between reward and inhibitory control networks.

category: neuroscience


title: Harmoni: a Method for Eliminating Spurious Interactions due to the Harmonic Components in Neuronal Data

abstract: Cross-frequency synchronization (CFS) has been proposed as a mechanism for integrating spatially and spectrally distributed information in the brain. However, investigating CFS in Magneto- and Electroencephalography (MEG/EEG) is hampered by the presence of spurious neuronal interactions due to the non-sinusoidal waveshape of brain oscillations. Such waveshape gives rise to the presence of oscillatory harmonics mimicking genuine neuronal oscillations. Until recently, however, there has been no methodology for removing these harmonics from neuronal data. In order to address this long-standing challenge, we introduce a novel method (called HARMOnic miNImization - Harmoni) that removes the signal components which can be harmonics of a non-sinusoidal signal. Harmonis working principle is based on the presence of CFS between harmonic components and the fundamental component of a non-sinusoidal signal. We extensively tested Harmoni in realistic EEG simulations. The simulated couplings between the source signals represented genuine and spurious CFS and within-frequency phase synchronization. Using diverse evaluation criteria, including ROC analyses, we showed that the within- and cross-frequency spurious interactions are suppressed significantly, while the genuine activities are not affected. Additionally, we applied Harmoni to real resting-state EEG data revealing intricate remote connectivity patterns which are usually masked by the spurious connections. Given the ubiquity of non-sinusoidal neuronal oscillations in electrophysiological recordings, Harmoni is expected to facilitate novel insights into genuine neuronal interactions in various research fields, and can also serve as a steppingstone towards the development of further signal processing methods aiming at refining within- and cross-frequency synchronization in electrophysiological recordings.


title: Beneficial effects of two Ayurvedic formulations, Saraswata Ghrita and Kalyanaka Ghrita on survival and on toxic aggregates in Drosophila models of Huntington's and Alzheimer's disease

abstract: In order to understand the health promotive, rejuvenative and disease preventive approach of the Ayurvedic system of medicine in the light of current principles, we examined two Rasayana formulations, viz., Kalayanaka Ghrita (KG) and Saraswata Ghrita (SG) for their effects in Alzheimers (AD) and Huntingtons (HD) neurodegenerative disease models of Drosophila. Initial experiments involving feeding of wild type flies on food supplemented with 0.05%, 0.25% and 0.5% (w/v) KG or SG revealed 0.05% to be without any adverse effect while higher concentrations caused dose-dependent reduction in pupation frequency and adult life span in wild type flies. Rearing GMR-GAL4>127Q (HD model) and ey-GAL4>A{beta}42 (AD model) larvae and adults on 0.05% or 0.25% SG or KG supplemented food enhanced the otherwise significantly reduced larval lethality and enhanced their median life span, with the 0.25% SG or KG concentrations being less effective than the 0.05%. In parallel with the better larval survival and enhanced adult life span, feeding the HD and AD model larvae on either of the Ghrita supplemented food (0.05% and 0.25%) substantially reduced the polyQ aggregates or amyloid plaques, respectively, in the larval eye discs. The present first in vivo organismic model study results have clinical implications for the increasing burden of age-associated dementia and neurodegenerative diseases like AD and HD in human populations.


title: Somnotate: A robust automated sleep stage classifier that exceeds human performance and identifies ambiguous states in mice

abstract: Manual sleep stage annotation is a time-consuming but often essential step in the analysis of sleep data. To address this bottleneck we need automated approaches that exhibit high levels of performance, are robust under different experimental conditions, are accessible, and meet the specific needs of sleep scientists. Here we develop an unbiased framework for assessing automated performance against a consensus derived from multiple experienced researchers. We then construct a new sleep stage classifier that combines automated feature extraction using linear discriminant analysis, with inference based on vigilance state-dependent contextual information using a hidden Markov model. This produces annotation accuracies that exceed expert performance on rodent electrophysiological data. We demonstrate that the classifier is robust to errors in the training data, compatible with different recording configurations, and maintains high performance during experimental interventions including sleep deprivation and optogenetic manipulations. Finally, the classifier quantifies and reports its certainty, which can be leveraged to identify ambiguous epochs for further examination. Ambiguous epochs are shown to include unsuccessful transitions between vigilance states, which may offer new insight into the mechanisms underlying sleep-wake dynamics. We call our classifier  Somnotate and make an implementation available to the neuroscience community.


title: Somnotate: An accurate, robust, and flexible sleep stage classifier for the experimentalist


title: Prediction and preview strongly affect reading times but not skipping during natural reading

abstract: In a typical text, readers look much longer at some words than at others and fixate some words multiple times, while skipping others altogether. Historically, researchers explained this variation via low-level visual or oculomotor factors, but today it is primarily explained in terms of cognitive factors, such as how well word identity can be predicted from context or discerned from parafoveal preview. While the existence of these effects has been well established in experiments, the relative importance of prediction, preview and low-level factors for eye movement variation in natural reading is unclear. Here, we address this question in three large datasets (n=104, 1.5 million words), using a deep neural network and Bayesian ideal observer to model linguistic prediction and parafoveal preview from moment to moment in natural reading. Strikingly, neither prediction nor preview was important for explaining word skipping - the vast majority of skipping was explained by a simple oculomotor model. For reading times, by contrast, we found strong but independent contributions of both prediction and preview, with effect sizes matching those from controlled experiments. Together, these results challenge dominant models of eye movements in reading by showing that linguistic prediction and parafoveal preview are not important determinants of word skipping.


title: Measurements of EEG Alpha Peak Frequencies Over the Lifespan: Validating Target Ranges on an In-Clinic Platform

abstract:
BackgroundPeak individual alpha frequencies (IAF) extracted from an EEG exam can provide novel sources of information regarding brain function. This information can help measure and track changes in cognition arising from conditions such as concussion or unhealthy aging.

Objective 1To validate a method for combining eyes-closed EEG with eyes-closed audio P300 ERP in order to streamline testing times involved in IAF extraction.

Objective 2To validate age-stratified target ranges of IAF collected in clinic against published research

Objective 3To validate the stability of IAF for data collected in routine clinical settings.

ParticipantsTwo thousand twenty-five subjects aged 13-90.

MethodsEEG with audio P300 was collected as part of a health screening exam for studies through Colorado University, Childrens Hospital Colorado, Boone Heart Institute, WAVi Co., and various clinics alongside other clinical evaluations.

Results(1) No differences were seen between IAF extracted during an eyes-closed resting and the P300 protocol. (2) The age-related AF trends measured in clinic match the age-related trends from previous research. (3) IAF remained stable over the course of 0-2 years in a test-retest dataset.

ConclusionIn-clinic measures of peak EEG frequency corroborate the age-related trends of published research taken over the last several decades and IAF. These results also confirm that IAF is a stable trait, making it useful for within-person longitudinal tracking. By following changes in IAF over time, deviations from normal CNS functioning, such as onset or progression of disease, can be monitored.


title: Diversity-induced plant history and soil history effects modulate plant responses to global change

abstract: Global change has dramatic impacts on grassland diversity. However, little is known about how fast species can adapt to these changes and how this affects their responses to global change. To close this gap, we performed a common garden experiment testing whether plant responses to global change are influenced by the selection history of the plants and the conditioning history of soil at different levels of plant diversity. Therefore, we collected seeds and took soil samples from 14-year old plant communities of a biodiversity experiment. Offspring of plants from low- and high-diversity communities were either grown in their own soil or in soil of a different community, and were either exposed to drought, increased nitrogen input, or a combination of both. Results show that, under nitrogen addition, offspring of plants selected at high diversity produced more biomass than those selected at low diversity, while drought neutralized differences in biomass production. Moreover, under the influence of global change drivers, mainly soil, and to a lesser extent plant history, influenced the expression of plant traits. Our results show that plant diversity modulates plant-soil interactions and growth strategies of plants, which feedback on the eco-evolutionary pathways of the plants and thus their responses to global change.

category: ecology


title: Kilometer-scale larval dispersal processes predict metapopulation connectivity pathways for Paramuricea biscaya in the northern Gulf of Mexico

abstract: Fine-scale larval dispersal and connectivity processes are key to species survival, growth, recovery and adaptation under rapidly changing disturbances. Quantifying both are required to develop any effective management strategy. In the present work, we examine the dispersal pattern and potential connectivity of a common deep-water coral, Paramuricea biscaya, found in the northern Gulf of Mexico by evaluating predictions of physical models with estimates of genetic connectivity. While genetic approaches provide estimates of realized connectivity, they do not provide information on the dispersal process. Physical circulation models can now achieve kilometer-scale resolution sufficient to provide detailed insight into the pathways and scales of larval dispersal. A high-resolution regional ocean circulation model is integrated for 2015 and its advective pathways are compared with the outcome of the genetic connectivity estimates of corals collected at six locations over the continental slope at depths comprised between 1000 and 3000 meters. Furthermore, the likely interannual variability is extrapolated using ocean hindcasts available for this basin. The general connectivity pattern exhibits a dispersal trend from east to west following 1000 to 2000-meter isobaths, corresponding to the overall westward near-bottom circulation. The connectivity networks predicted by our model were mostly congruent with the estimated genetic connectivity patterns. Our results show that although dispersal distances of 100 km or less are common, depth differences between tens to a few hundred meters can effectively limit larval dispersal. A probabilistic graphic model suggests that stepping-stone dispersal mediated by intermediate sites provides a likely mechanism for long-distance connectivity between the populations separated by distances of 300 km or greater, such as those found in the DeSoto and Keathley canyons.


title: Range contractions of the Broad-winged Hawk in the Northeast United States

abstract: The Broad-winged Hawk (BWHA, Buteo platypterus) is a small, secretive hawk with distinguishing broad black tail bands that breeds in northeastern North America. The hawk nests in deciduous or mixed forest, often near water, and close to clearings or forest edges. Land conversion and fragmentation alters the landscape and reduces the area of contiguous forest used by BWHA. This study seeks to determine the landscape characteristics influencing the apparent breeding range declines of the BWHA at the landscape scale. Landscape characteristics and BWHA presence data from 18,684 Breeding Bird Atlas blocks (each about 25km2) from Ohio, West Virginia, Maryland, Pennsylvania, and New York for two atlas period (1st Atlas: 1980s, 2nd Atlas: 2000s) were analyzed. Bayesian latent Gaussian models were fitted using INLA to determine best fit model for predicting the landscape characteristics associated with BWHA presence. The best models included landscape changes in land cover, including forest, water, urban, barren, farmland, and wetland and fragmentation of the landscape. Trends in loss were especially prevalent around the regions largest cities: New York, Philadelphia, Baltimore and Washington DC. Loss of BWHA at the block-level was associated with areas with less forest in the 2000s, a decline in size of largest forest patches, lower elevations and lower latitudes. We suggest that both habitat loss and climate change may be contributing to the range contraction of the Broad-winged Hawk in the northeast United States.


title: Identification of sex chromosomes using genomic and cytogenetic methods in a range-expanding spider, Argiope bruennichi (Araneae: Araneidae)

abstract: Differences between sexes in growth, ecology and behavior strongly shape species biology. In some animal groups, such as spiders, it is difficult or impossible to identify the sex of juveniles based on external morphology. This information would be useful for field surveys, behavioral experiments, and ecological studies on e.g. sex ratios and dispersal. In species with sex chromosomes, sex can be determined based on the specific sex chromosome complement. Additionally, information on the sequence of sex chromosomes provides the basis for studying sex chromosome evolution. We combined cytogenetic and genomic data to identify the sex chromosomes in the sexually dimorphic spider Argiope bruennichi, and designed RT-qPCR sex markers. We found that genome size and GC content of this spider falls into the range reported for the majority of araneids. The male karyotype is formed by 24 acrocentric chromosomes with an X1X20 sex chromosome system, with little similarity between X chromosomes, suggesting origin of these chromosomes by X chromosome fission or early duplication of an X chromosome and subsequent independent differentiation of the copies. Our data suggest similarly sized X chromosomes in A. bruennichi. They are smaller chromosomes of the complement. Our findings open the door to new directions in spider evolutionary and ecological research.

category: genetics


title: Fecal microbiota and their association with heat stress in Bos taurus

abstract: Humans have been influencing climate changes by burning fossil fuels, farming livestock, and cutting down rainforests, which has led to global temperature rise. This problem of global warming affects animals by causing heat stress, which negatively affects their health, biological functions, and reproduction. On the molecular level, it has been proved that heat stress changes the expression level of genes and therefore causes changes in proteome and metabolome. The importance of a microbiome in many studies showed that it is considered as individuals "second genome". Physiological changes caused by heat stress may impact the microbiome composition. In this study, we identified fecal microbiota associated with heat stress that was quantified by three metrics - rectal temperature, drooling, and respiratory scores and represented by their Estimated Breeding Values. For this purpose, the 16S rRNA sequencing technique was used. We analyzed the microbiota from 136 fecal samples of Chinese Holstein cows through a 16S rRNA gene sequencing approach. Sequence data were processed using a pipeline involving QIIME2 software together with SILVA database. Statistical modeling was performed using a negative binomial regression. The analysis revealed the total number of 24 genera and 12 phyla associated with heat stress metrics. Rhizobium and Pseudobutyrivibrio turned out to be the most significant genera, while Acidobacteria and Gemmatimonadetes were the most significant phyla. Phylogenetic analysis revealed that that three heat stress indicators quantify different metabolic ways of animals reaction to heat stress. Other studies already identified that those genera had significantly increased abundance in mice exposed to stressor-induced changes. Moreover, identified microbiota significantly associated with heat stress measures were mostly mesophilic, so their association seems to be due to heat stress-induced secondary, metabolic changes, and not directly by temperature. Moreover, high fold changes of many genera suggest that they may be used as biomarkers for monitoring the level of heat stress in cattle.


title: The long-term effects of repeated heroin vapor inhalation during adolescence on measures of nociception and anxiety-like behavior in adult Wistar rats

abstract: The United States continues to experience a public health crisis related to the nonmedical use of opioid drugs. Adolescents represent a vulnerable group due to increased experimentation with illicit substances that is often associated with the adolescent period, and because adolescent drug use can result in long-term effects that differ from those caused by drug use initiated during adulthood. The present study examined the effects of repeated heroin inhalation, a common route of administering opioids, during adolescence, on measures of nociception and anxiety-like behavior in adulthood. Rats were exposed twice daily to 30-minutes of heroin vapor from post-natal day (PND) 36 to PND 45. At 12 weeks of age, baseline thermal nociception was assessed across a range of temperatures with a warm-water tail-withdrawal assay. Anxiety-like behavior was assessed in an elevated plus-maze (EPM) and activity was measured in an open field arena. Starting at 23 weeks of age, baseline thermal nociception was re-assessed, nociception was determined after acute heroin or naloxone injection, and anxiety-like behavior was redetermined in the EPM. Adolescent heroin inhalation altered baseline thermal nociception in female rats at 12 weeks of age and in both female and male rats at [~]23 weeks. Heroin-treated animals exhibited an increase in anxiety-like behavior when tested in the elevated plus-maze, showed blunted heroin-induced analgesia, but exhibited no effect on naloxone-induced hyperalgesia. The present study demonstrates for the first time that heroin vapor inhalation during adolescence produces behavioral and physiological consequences in rats that persist well into adulthood.

category: pharmacology and toxicology


title: Software choice and depth of sequence coverage can impact plastid genome assembly - A case study in the narrow endemic Calligonum bakuense

abstract: Most plastid genome sequences are assembled from short-read whole-genome sequencing data, yet the impact that sequence coverage and the choice of assembly software can have on the accuracy of the resulting assemblies is poorly understood. In this study, we test the impact of both factors on plastid genome assembly in the threatened and rare endemic shrub Calligonum bakuense, which forms a distinct lineage in the genus Calligonum. We aim to characterize the differences across plastid genome assemblies generated by different assembly software tools and levels of sequence coverage and to determine if these differences are large enough to affect the phylogenetic position inferred for C. bakuense. Four assembly software tools (FastPlast, GetOrganelle, IOGA, and NOVOPlasty) and three levels of sequence coverage (original depth, 2,000x, and 500x) are compared in our analyses. The resulting assemblies are evaluated with regard to reproducibility, contig number, gene complement, inverted repeat length, and computation time; the impact of sequence differences on phylogenetic tree inference is also assessed. Our results show that software choice can have a considerable impact on the accuracy and reproducibility of plastid genome assembly and that GetOrganelle produced the most consistent assemblies for C. bakuense. Moreover, we found that a cap in sequence coverage can reduce both the sequence variability across assembly contigs and computation time. While no evidence was found that the sequence variability across assemblies was large enough to affect the phylogenetic position inferred for C. bakuense, differences among the assemblies may influence genotype recognition at the population level.

category: plant biology


title: Temporal phenomic predictions from unoccupied aerial systems can outperform genomic predictions

abstract: A major challenge of genetic improvement and selection is to accurately predict individuals with the highest fitness in a population without direct measurement. Over the last decade genomic predictions (GP) based on genome-wide markers have become reliable and routine. Now phenotyping technologies, including unoccupied aerial systems (UAS also known as drones), can characterize individuals with a data depth comparable to genomics when used throughout growth. This study, for the first time, demonstrated that the prediction power of temporal UAS phenomic data can achieve or exceed that of genomic data. UAS data containing red-green-blue (RGB) bands over fifteen growth time points and multispectral (RGB, red-edge and near infrared) bands over twelve time points were compared across 280 unique maize hybrids. Through cross validation of untested genotypes in tested environments (CV2), temporal phenomic prediction (TPP) outperformed GP (0.80 vs 0.71); TPP and GP performed similarly in three other cross validation scenarios. Genome wide association mapping using area under temporal curves of vegetation indices (VIs) revealed 24.5 percent of a total of 241 discovered loci (59 loci) had associations with multiple VIs, explaining up to 51 percent of grain yield variation, less than GP and TPP predicted. This suggests TPP, like GP, integrates small effect loci well improving plant fitness predictions. More importantly, temporal phenomic prediction appeared to work successfully on unrelated individuals unlike genomic prediction.


title: In silico analysis for potential proteins and microRNAs in Glioblastoma and Parkinsonism.

abstract: In todays world, neurodegenerative diseases such as Alzheimers disease, Parkinsons Disease, Huntingtons Disease as well as brain cancers such as astrocytomas, ependymomas, glioblastomas have become a great threat to us. In this study, we are trying to find a probable molecular connection associated with two very much different diseases, Glioblastoma, also known as Glioblastoma Multiforme (cancers of microglial cells of our brain) and Parkinsons disease. We at first downloaded the microarray datasets of these two diseases from Gene Expression Omnibus (GEO) and then analyzed them by the GEO2R tool. After analysis, we found 249 common upregulated differential expressed genes and 135 common downregulated differential expressed genes of these two diseases. Therefore the common differentially expressed genes, both upregulated and downregulated, were imported into STRING online tool to find out the protein-protein interactions. Now, this whole network was subjected to Cytoscape and the top ten hub genes were found by Cyto-Hubba plug-in. The top then hub genes are EGFR, CCNB1, CDK1, CCNA2, CHEK1, RAD51, MAD2L1, KIF20A, BUB1, and CCNB2. These all genes are upregulated in both diseases. To find out the biological processes, molecular functions, cellular components, and pathways associated with these hub genes Enrichr online software was used. We used miRNet software to determine the interactions of hub genes with microRNAs. This study will be useful in the future for drug targets discovery for these diseases.


title: Novelty detection in early olfactory processing of the honey bee, Apis mellifera

abstract: Animals are constantly bombarded with stimuli, which presents a fundamental problem of sorting among pervasive uninformative stimuli and novel, possibly meaningful stimuli. We evaluated novelty detection behaviorally in honey bees as they position their antennae differentially in an air stream carrying familiar or novel odors. We then characterized neuronal responses to familiar and novel odors in the first synaptic integration center in the brain - the antennal lobes. We found that the neurons that exhibited stronger initial responses to the odor that was to be familiarized are the same units that later distinguish familiar and novel odors, independently of chemical identities. These units, including both projection neurons and local neurons, showed a decreased response to the familiar odor but an increased response to the novel odor. Our results suggest that the antennal lobe may assign a category of familiarity or novelty to an odor stimulus in addition to its chemical identity code. Therefore, the mechanisms for novelty detection may be present in early sensory processing, either as a result of local synaptic interaction or via feedback from higher brain centers.


title: Seascape genomics reveals metapopulation connectivity network of Paramuricea biscaya in the northern Gulf of Mexico

abstract: The degree of connectivity among populations influences their ability to respond to natural and anthropogenic stressors. In marine systems, determining the scale, rate, and directionality of larval dispersal is therefore central to understanding how coral metapopulations are interconnected and the degree of resiliency in the event of a localized disturbance. Understanding these source-sink dynamics is essential to guide restoration efforts and for the study of ecology and evolution in the ocean. The patterns and mechanisms of connectivity in the deep-sea (> 200 meters deep) are largely understudied. In this study, we investigated the spatial diversity patterns and metapopulation connectivity of the octocoral Paramuricea biscaya throughout the northern Gulf of Mexico (GoM). Paramuricea biscaya is one of the most abundant corals on the lower continental slope (between 1200 and 2500 m) in the GoM. The 2010 Deepwater Horizon oil spill (DWH) directly impacted populations of this species and thus are considered primary targets for restoration. We used a combination of seascape genomic analyses, high-resolution ocean circulation modeling, and larval dispersal simulations to quantify the degree of population structuring and connectivity among P. biscaya populations. Evidence supports the hypotheses that the genetic diversity of P. biscaya is predominantly structured by depth, and that larval dispersal among connected populations is asymmetric due to dominant ocean circulation patterns. Our results suggest that there are intermediate unsampled populations in the central GoM that serve as stepping stones for dispersal. The data suggest that the DeSoto Canyon area, and possibly the West Florida Escarpment, critically act as sources of larvae for areas impacted by the DWH oil spill in the Mississippi Canyon. This work illustrates that the management of deep-sea marine protected areas should incorporate knowledge of connectivity networks and depth-dependent processes throughout the water column.


title: Predicting selection-response gradients of heat tolerance in a wide-ranging reef-building coral

abstract:
Ocean temperatures continue to rise due to climate change but it is unclear if heat tolerance of marine organisms will keep pace with warming. Understanding how tolerance scales from individuals to species and quantifying adaptive potentials is essential to forecasting responses to warming. We reproductively crossed corals from a globally distributed species (Acropora tenuis) on the Great Barrier Reef (Australia) from three thermally distinct reefs to create 85 offspring lineages. Individuals were experimentally exposed to temperatures (27.5, 31, and 35.5 {degrees}C) in adult and two critical early life stages (larval and settlement) to assess acquired heat tolerance via outcrossing of offspring phenotypes by comparing five physiological responses (photosynthetic yields, bleaching, necrosis, settlement, and survival). Adaptive potentials and physiological reaction norms were calculated across three stages to integrate heat tolerance at different biological scales. Selective breeding improved larval survival to heat by 1.5 - 2.5x but did not result in substantial enhancement of settlement, although population crosses were significantly different. At heat, adults were less variable compared to larval responses in warmer reefs compared to the cooler reef. Adults and offspring also differed in their mean population responses, likely underpinned by heat stress imposing strong divergent selection on adults. These results have implications for downstream selection during reproduction, evidenced by variability in a conserved heat tolerance response across offspring lineages. These results inform our ability to forecast the impacts of climate change on wild populations of corals and will aid in developing novel conservation tools like the assisted evolution of at-risk species.

SUMMARY STATEMENTHeat stress exerts disruptive selection on adult corals. This likely underpins variability in offspring survival and results in differences in offspring responses to selection.


title: Phylogenetic analyses for orthogroup-based classification of GDSL-type esterase/lipase (GELP) family in angiosperms

abstract:
GDSL-type esterase/lipase (GELP) enzymes have multiple functions in plants, spanning from developmental processes to the response to biotic and abiotic stresses. Genes encoding GELP belong to a large gene family with several tens to more than hundred members per angiosperm species.

Here, we applied iterative phylogenetic analyses to identify 10 main clusters subdivided into 44 expert-curated reference orthogroups (OGs) using three monocot and five dicot genomes. Our results show that some GELP OGs expanded while others were maintained as single copy genes.

This semi-automatic approach proves to be effective to characterize large gene families and provides a solid classification framework for the GELP members in angiosperms. The orthogroup-based reference will be useful to perform comparative studies, infer gene functions and better understand the evolutionary history of this gene family.


title: Phylogenetic analyses for orthogroup-based classification of GDSL-type esterase/lipase (GELP) family in angiosperm representative species


title: Multi-omics analysis of multiple glucose-sensing receptor systems in yeast

abstract: The yeast Saccharomyces cerevisiae has long been used to produce alcohol from glucose and other sugars. While much is known about glucose metabolism, relatively little is known about the receptors and signaling pathways that indicate glucose availability. Here we compare the two glucose receptor systems in S. cerevisiae. The first is a heterodimer of transporter-like proteins (transceptors), while the second is a seven-transmembrane receptor coupled to a large G protein (Gpa2) and two small G proteins (Ras1 and Ras2). Through comprehensive measurements of glucose-dependent transcription and metabolism, we demonstrate that the two receptor systems have distinct roles in glucose signaling: the G protein-coupled receptor directs carbohydrate and energy metabolism, while the transceptors regulate ancillary processes such as ribosome, amino acids, cofactor and vitamin metabolism. The large G protein transmits the signal from its cognate receptor, while the small G protein Ras2 (but not Ras1) integrates responses from both receptor pathways. Collectively, our analysis reveals the molecular basis for glucose detection and the earliest events of glucose-dependent signal transduction in yeast.


title: Exo-enzymatic addition of diazirine-modified sialic acid to cell surfaces enables photocrosslinking of glycoproteins

abstract: Glycan binding often mediates extracellular macromolecular recognition events. Accurate characterization of these binding interactions can be difficult because of dissociation and scrambling that occur during purification and analysis steps. Use of photocrosslinking methods has been pursued to covalently capture glycan-dependent interactions in situ however use of metabolic glycan engineering methods to incorporate photocrosslinking sugar analogs is limited to certain cell types. Here we report an exo-enzymatic labeling method to add a diazirine-modified sialic acid (SiaDAz) to cell surface glycoconjugates. The method involves chemoenzymatic synthesis of diazirine-modified CMP-sialic acid (CMP-SiaDAz), followed by sialyltransferase-catalyzed addition of SiaDAz to desialylated cell surfaces. Cell surface SiaDAz-ylation is compatible with multiple cell types and is facilitated by endogenous extracellular sialyltransferase activity present in Daudi B cells. This method for extracellular addition of 2-6-linked SiaDAz enables UV-induced crosslinking of CD22, demonstrating the utility for covalent capture of glycan-mediated binding interactions.

category: biochemistry


title: Impairment of the cellulose degradation machinery enhances fungal virulence but limits reproductive fitness

abstract:
Fungal endophytes grow in the apoplastic space, in constant contact with the plant cell wall (CW) that hinders microbe progression, while representing a source of nutrients. Although numerous fungal CW modifying proteins have been identified, their role during host colonization remains underexplored. Here we show that the root-infecting plant pathogen Fusarium oxysporum (Fo) does not require its complete arsenal of cellulases to infect the host plant. Quite the opposite, Fo mutants impaired in cellulose degradation become hypervirulent by enhancing the secretion of virulence factors. On the other hand, the reduction on cellulase activity had a severe negative effect on saprophytic growth and microconidia production during the final stages of the Fo infection cycle. These findings enhance our understanding on the function of plant CW degradation on the outcome of host-microbe interactions and reveal an unexpected role of cellulose degradation in a pathogens evolutionary success.

TeaserUnexpectedly, fungi compromised in their capacity to degrade plant cellulose are hypervirulent but impaired in sporulation.

category: microbiology


title: A Systematic Interrogation of MHC Class I Antigen Presentation Identifies Constitutive and Compensatory Protein Degradation Pathways

abstract: Protein degradation products are constitutively presented as peptide antigens by MHC Class I. While hypervariability of Class I genes is known to tremendously impact antigen presentation, whether differential function of protein degradation pathways (comprising >1000 genes) could alter antigen generation remains poorly understood apart from a few model substrates. In this study, we introduce normalization methods for quantitative antigen mass spectrometry and confirm that most Class I antigens are dependent on ubiquitination and proteasomal degradation. Remarkably, many antigens derived from mitochondrial inner membrane proteins are not. Additionally, we find that atypical antigens can arise from compensatory protein degradation pathways, such as an increase in mitochondrial and membrane protein antigen presentation upon proteasome inhibition. Notably, incomplete inhibition of protein degradation pathways may have clinical utility in cancer immunotherapy, as evidenced by appearance of novel antigens upon partial proteasome inhibition.

category: molecular biology


title: Analysis of DYRK1B, PPARG, and CEBPB Expression Patterns in Adipose-Derived Stem Cells from Patients carrying DYRK1B R102C and Healthy Individuals During Adipogenesis

abstract:
BackgroundMetabolic syndrome (MetS) is a group of signs and symptoms that are associated with higher risk of Type 2 Diabetes Mellitus (T2DM) and Cardiovascular Diseases (CVDs). The major risk factor for developing MetS is abdominal obesity that is caused by increase in adipocyte size or number. Adipocyte number multiplication is caused by differentiation of mesenchymal stem cells into adipose tissue. Numerous studies have evaluated the expression of key transcription factors including PPARG and CEBPB during adipocyte differentiation in murine cells such as 3T3-L1 cell line. In order to comprehend the expression changes during the process of fat accumulation in adipose tissue derived stem cells (ASCs), we compared the expression of DYRK1B, PPARG, and CEBPB in undifferentiated and differentiated ASCs into mature adipocytes between the patient (harboring DYRK1b R102C) and control (healthy individuals) groups.

MethodsGene expression was evaluated on eighth days pre-induction and day 1, 5 and 15 post-induction. The pluripotent capacity of ASCs and the potential for differentiation into adipocyte were confirmed by flow cytometry analysis of surface markers (CD34, CD44, CD105 and CD90), and Oil red O staining, respectively. Expression of DYRK1B, PPARG, and CEBPB were assessed by RT-PCR in patients and normal individuals samples.

ResultsThe expression of DYRK1B kinase and transcription factors (CEBPB and PPARG) are significantly higher in adipose derived stem cells harboring DYRK1b R102C compared to non-carriers on day 5 and 15 during adipocyte differentiation. These proteins may be suitable targets for therapeutic strategies in obesity and obesity related disorders like metabolic syndrome. Furthermore, AZ191 exhibited a potent and selectively inhibitory activity toward DYRK1B and CEBPB.

ConclusionCEBPB, PPARA and DYRK1B contribute to adipogenesis and the development of metabolic syndrome; thus, they can be harnessed in developing therapeutic agents against metabolic syndrome.

category: cell biology


title: Cytoplasmic dynein-1 cargo diversity is mediated by the combinatorial assembly of FTS-Hook-FHIP complexes

abstract: In eukaryotic cells, intracellular components are organized by the microtubule motors cytoplasmic dynein-1 (dynein) and kinesins, which are linked to cargos via adaptor proteins. While [~]40 kinesins transport cargo toward the plus end of microtubules, a single dynein moves cargo in the opposite direction. How dynein transports a wide variety of cargos remains an open question. The FTS-Hook-FHIP ("FHF") cargo adaptor complex links dynein to cargo in mammals and fungi. As human cells have three Hooks and four FHIP proteins, we hypothesized that the combinatorial assembly of different Hook and FHIP proteins could underlie dynein cargo diversity. Using proteomic approaches, we determine the protein  interactome of each FHIP protein. Live-cell imaging and biochemical approaches show that different FHF complexes associate with distinct motile cargos. These complexes also move with dynein and its cofactor dynactin in single-molecule in vitro reconstitution assays. Complexes composed of FTS, FHIP1B, and Hook1/Hook3 co-localize with Rab5-tagged early endosomes via a direct interaction between FHIP1B and GTP-bound Rab5. In contrast, complexes composed of FTS, FHIP2A and Hook2 colocalize with Rab1A-tagged ER-to-Golgi cargos and FHIP2A is involved in the motility of Rab1A tubules. Our findings suggest that combinatorial assembly of different FTS-Hook-FHIP complexes is one mechanism dynein uses to achieve cargo specificity.


title: Activation of the essential kinase PDK1 by phosphoinositide-driven autophosphorylation

abstract: 3-phosphoinositide-dependent kinase 1 (PDK1) is an essential serine/threonine protein kinase, which plays a crucial role in cell growth and proliferation. It is often referred to as a  master kinase due to its ability to activate at least 23 downstream protein kinases implicated in various signaling pathways. In this study, we have elucidated the mechanism of phosphoinositide-driven PDK1 auto-activation. We show that PDK1 trans- autophosphorylation is mediated by a PIP3-mediated face-to-face dimer. We report regulatory motifs in the kinase-PH interdomain linker that allosterically activate PDK1 autophosphorylation via a linker-swapped dimer mechanism. Finally, we show that PDK1 is autoinhibited by its PH domain and that positive cooperativity of PIP3 binding drives switch- like activation of PDK1. Our work implies that the PDK1-mediated activation of effector kinases, including Akt, PKC, Sgk, S6K and RSK, many of whom are not directly regulated by phosphoinositides, is also likely to be dependent on PIP3 or PI(3, 4)P2.


title: Membrane fission during bacterial spore development requires DNA-driven cellular inflation

abstract:
Bacteria require membrane fission for cell division and endospore formation. FisB catalyzes membrane fission during sporulation, but the molecular basis is unclear as it cannot remodel membranes by itself. Sporulation initiates with an asymmetric division that generates a large mother cell and a smaller forespore that contains only 1/4 of its complete genome. As the mother cell membranes engulf the forespore, a DNA translocase pumps the rest of the chromosome into the small forespore compartment, inflating it due to increased turgor. When the engulfing membranes undergo fission, the forespore is released into the mother cell cytoplasm. Here we show that forespore inflation and FisB accumulation are both required for efficient membrane fission. We suggest that high membrane tension in the engulfment membrane caused by forespore inflation drives FisB-catalyzed membrane fission. Collectively our data indicate that DNA-translocation has a previously unappreciated second function in energizing FisB-mediated membrane fission under energy-limited conditions.

HIGHLIGHTS- Membrane fission during endospore formation requires rapid forespore inflation by ATP-driven DNA translocation.
- Forespore inflation is fast enough to increase the tension of the engulfment and forespore membranes to near lysis tensions.
- FisB catalyzes membrane fission by impeding the flux of lipids that partially support forespore and engulfment membrane growth.
- Membrane fission utilizes chemical energy transduced to mechanical energy during DNA-packing into the forespore.

category: biophysics


title: Membrane fission during bacterial spore development requires DNA-pumping driven cellular inflation


title: RNA length has a non-trivial effect in the stability of biomolecular condensates formed by RNA-binding proteins

abstract: Biomolecular condensates formed via liquid-liquid phase separation (LLPS) play a crucial role in the spatiotemporal organization of the cell material. Nucleic acids can act as critical modulators in the stability of these protein condensates. Here, we present a multiscale computational strategy, exploiting the advantages of both a sequence-dependent coarse-grained representation of proteins and a minimal coarse-grained model that represent proteins as patchy colloids, to unveil the role of RNA length in regulating the stability of RNA-binding protein (RBP) condensates. We find that for a constant RNA/protein ratio in which phase separation is enhanced, the protein fused in sarcoma (FUS), which can phase separate on its own--i.e., via homotypic interactions--only exhibits a mild dependency on the RNA strand length, whereas, the 25-repeat proline-arginine peptide (PR25), which does not undergo LLPS on its own at physiological conditions but instead exhibits complex coacervation with RNA--i.e., via heterotypic interactions--shows a strong dependence on the length of the added RNA chains. Our minimal patchy particle simulations, where we recapitulate the modulation of homotypic protein LLPS and complex coacervation by RNA length, suggest that the strikingly different effect of RNA length on homotypic LLPS versus complex coacervation is general. Phase separation is RNA-length dependent as long as the relative contribution of heterotypic interactions sustaining LLPS is comparable or higher than that committed by protein homotypic interactions. Taken together, our results contribute to illuminate the intricate physicochemical mechanisms that influence the stability of RBP condensates through RNA inclusion.


title: Single-cell analysis of basal cell carcinoma reveals heat shock proteins promote tumor growth in response to WNT5A-mediated inflammatory signals

abstract: How basal cell carcinoma (BCC) interacts with its tumor microenvironment to promote growth is unclear. Here we use singe-cell RNA sequencing to define the human BCC ecosystem and discriminate between normal and malignant epithelial cells. We identify spatial biomarkers of both tumors and their surrounding stroma that reinforce the heterogeneity of each tissue type. Combining pseudotime, RNA velocity, cellular entropy, and regulon analysis in stromal cells reveal a cancer-specific rewiring of fibroblasts where STAT1, TGF-{beta}, and inflammatory signals induce a non-canonical WNT5A program that maintains the stromal inflammatory state. Cell-cell communication modeling suggests that tumors respond to the sudden burst of fibroblast-specific inflammatory signaling pathways by producing heat shock proteins, which we validated in situ. Finally, dose-dependent treatment with an HSP70 inhibitor suppresses in vitro BCC cell growth and Hedgehog signaling and in vivo tumor growth in a BCC mouse model, validating HSP70s essential role in tumor growth and reinforcing the critical nature of tumor microenvironment crosstalk in BCC progression.

category: cancer biology


title: CRISPR-Cas9-mediated knockout of CYP79D1 and CYP79D2 in cassava attenuates toxic cyanogen production

abstract: Cassava (Manihot esculenta Crantz) is a starchy root crop that supports over a billion people in tropical and subtropical regions of the world. This staple, however, produces toxic cyanogenic compounds and requires processing for safe consumption. Excessive consumption of insufficiently processed cassava, in combination with protein-poor diets, can have neurodegenerative impacts. Reducing the cyanogen content by conventional breeding is problematic due to the heterozygous nature of the crop; recombination will generally disrupt a clonally propagated cultivars suite of desirable traits. To reduce cyanide levels in cassava, we used CRISPR-mediated mutagenesis to disrupt the cytochrome P450 genes CYP79D1 and CYP79D2 whose protein products catalyze the first step in cyanogenic glucoside biosynthesis. Knockout of both genes eliminated cyanide in leaves and storage roots of cassava accession 60444 and the West African, farmer-preferred cultivar TME 419. Although knockout of CYP79D2 alone resulted in significant reduction of cyanide, mutagenesis of CYP79D1 did not, indicating these paralogs have diverged in their function. Our work demonstrates cassava genome editing for food safety, reduced processing requirements, and environmental benefits that could be readily extended to other farmer-preferred cultivars.


title: Infection-induced miR-126 suppresses tsc1- and cxcl12a-dependent permissive macrophages during mycobacterial infection

abstract: Regulation of host microRNA (miRNA) expression is a contested node that controls the host immune response to mycobacterial infection. The host must overcome concerted subversive efforts of pathogenic mycobacteria to launch and maintain a protective immune response. Here we examine the role of miR-126 in the zebrafish model of Mycobacterium marinum infection and identify a protective role for this infection-induced miRNA through multiple effector pathways. Specifically, we analyse the impact of the miR-126 knockdown-induced tsc1a and cxcl12a/ccl2/ccr2 signalling axes during early host-M. marinum interactions. We find a strong detrimental effect of tsc1a upregulation that renders zebrafish embryos susceptible to higher bacterial burden and increased cell death despite dramatically higher recruitment of macrophages to the site of infection. We demonstrate that infection-induced miR-126 suppresses tsc1 and cxcl12a expression thus improving macrophage function early in infection, partially through activation of mTOR signalling and strongly through preventing the recruitment of Ccr2+ permissive macrophages, resulting in the recruitment of protective tnfa-expressing macrophages. Together our results demonstrate an important role for infection-induced miR-126 in shaping an effective immune response to M. marinum infection in zebrafish embryos.

category: immunology


title: Suppression of B-Cell Activation by Human Cord Blood-Derived Stem Cells (CB-SC) through the Galectin-9-Dependent Cell Contact Mechanism

abstract:
BackgroundWe developed the Stem Cell Educator therapy among multiple clinical trials based on the immune modulations of multipotent cord blood-derived stem cells (CB-SC) on different compartments of immune cells such as T cells and monocytes/macrophages in diabetes and other autoimmune diseases. However, the effects of CB-SC on the B cells remained unclear. To better understand the molecular mechanisms underlying the immune education of CB-SC, we explored the modulations of CB-SC on human B cells.

MethodsCB-SC were isolated from human cord blood units and confirmed by flow cytometry with different markers for their purity. B cells were purified by using anti-CD19 immunomagnetic beads from human peripheral blood mononuclear cells (PBMC). Next, the activated B cells were treated in the presence or absence of coculture with CB-SC for 7 days before undergoing flow cytometry analysis of phenotypic change with different markers. RT-PCR was utilized to evaluate the levels of galectin expressions with or without treatment of activated B cells in order to find the key galectin contributing to the B-cell modulation.

ResultsFlow cytometry demonstrated that the proliferation of activated B cells was markedly suppressed in the presence of CB-SC, leading to the down-regulation of immunoglobulin productions from the activated B cells. Phenotypic analysis revealed that treatment with CB-SC increased the percentage of IgD+CD27- naive B cells, but decreased the percentage of IgD-CD27+ switched B cells. Transwell assay showed that the immune suppression of CB-SC on B cells was dependent on the manner of cell-cell contact via Gal-9 molecule, as confirmed by the blocking experiment with the anti-Gal-9 monoclonal antibody. Mechanistic studies demonstrated that both calcium levels of cytoplasm and mitochondria were down-regulated after the treatment with CB-SC, causing the decline of mitochondrial membrane potential in the activated B cells. Western blot exhibited that the levels of phosphorylated Akt and Erk1/2 signaling proteins in the activated B cells were also markedly reduced in the presence of CB-SC.

ConclusionsCB-SC displayed multiple immune modulations on B cells through the Gal-9-mediated cell-cell contact mechanism and calcium flux/Akt/Erk1/2 signaling pathways. The data advances current understanding about the molecular mechanisms underlying the Stem Cell Educator therapy to treat autoimmune diseases in clinics.


title: In silico analyses of penicillin binding proteins in Burkholderia pseudomallei uncovers SNPs with utility for phylogeography, species differentiation, and sequence typing

abstract:
BackgroundBurkholderia pseudomallei causes melioidosis. Sequence typing this pathogen can reveal geographical origin and uncover epidemiological associations. Here, we describe B. pseudomallei genes encoding putative penicillin binding proteins (PBPs) and investigate their utility for determining phylogeography and differentiating closely related species.

Methodology & Principal FindingsWe performed in silico analysis to characterize 10 PBP homologs in B. pseudomallei 1026b. As PBP active site mutations can confer {beta}-lactam resistance in Gram-negative bacteria, PBP sequences in two resistant B. pseudomallei strains were examined for similar alterations. Sequence alignments revealed single amino acid polymorphisms (SAAPs) unique to the multidrug resistant strain Bp1651 in the transpeptidase domains of two PBPs, but not directly within the active sites. Using BLASTn analyses of complete assembled genomes in the NCBI database, we determined genes encoding PBPs were conserved among B. pseudomallei (n=101) and Burkholderia mallei (n=26) strains. Within these genes, single nucleotide polymorphisms (SNPs) useful for predicting geographic origin of B. pseudomallei were uncovered. SNPs unique to B. mallei were also identified. Based on 11 SNPs identified in two genes encoding predicted PBP-3s, a dual-locus sequence typing (DLST) scheme was developed. The robustness of this typing scheme was assessed using 1,523 RefSeq genomes from B. pseudomallei (n=1,442) and B. mallei (n=81) strains, resulting in 32 sequence types (STs). Compared to multi-locus sequence typing (MLST), the DLST scheme demonstrated less resolution to support the continental separation of Australian B. pseudomallei strains. However, several STs were unique to strains originating from a specific country or region. The phylogeography of Western Hemisphere B. pseudomallei strains was more highly resolved by DLST compared to internal transcribed spacer (ITS) typing, and all B. mallei strains formed a single ST.

SignificanceConserved genes encoding PBPs in B. pseudomallei are useful for strain typing, can enhance predictions of geographic origin, and differentiate strains of closely related Burkholderia species.

Author SummaryBurkholderia pseudomallei causes the life-threatening disease melioidosis and is considered a biological threat and select agent by the United States government. This soil-dwelling bacterium is commonly found in regions of southeast Asia and northern Australia, but it is also detected in other tropical and sub-tropical areas around the world. With a predicted global burden of 165,000 annual cases and mortality rate that can exceed 40% without prompt and appropriate antibiotic treatment, understanding the epidemiology of melioidosis and mechanisms of antibiotic resistance in B. pseudomallei can benefit public health and safety. Recently, we identified ten conserved genes encoding putative penicillin binding proteins (PBPs) in B. pseudomallei. Here, we examined B. pseudomallei PBP sequences for amino acid mutations that may contribute to {beta}-lactam resistance. We also uncovered nucleotide mutations with utility to predict the geographical origin of B. pseudomallei strains and to differentiate closely related Burkholderia species. Based on 11 informative single nucleotide polymorphisms in two genes each encoding a PBP-3, we developed a simple, targeted dual-locus typing approach.


title: A Snapshot of the Global Drinking Water Virome: Diversity and Metabolic Potential Vary with Residual Disinfectant Use

abstract: Viruses are important drivers of microbial community ecology and evolution, influencing microbial mortality, metabolism, and horizontal gene transfer. However, the effects of viruses remain largely unknown in many environments, including in drinking water systems. Drinking water metagenomic studies have offered a whole community perspective of bacterial impacts on water quality, but have not yet considered the influences of viruses. In this study, we address this gap by mining viral DNA sequences from publicly available drinking water metagenomes from distribution systems in six countries around the world. These datasets provide a snapshot of the taxonomic diversity and metabolic potential of the global drinking water virome; and provide an opportunity to investigate the effects of geography, climate, and drinking water treatment practices on viral diversity. Both environmental conditions and differences in sample processing were found to influence the viral composition. Using free chlorine as the residual disinfectant was associated with clear differences in viral taxonomic diversity and metabolic potential, with significantly fewer viral populations and less even viral community structures than observed in distribution systems without residual disinfectant. Additionally, drinking water viruses carry antibiotic resistance genes (ARGs), as well as genes to survive oxidative stress and nitrogen limitation. Through this study, we have demonstrated that viral communities are diverse across drinking water systems and vary with the use of residual disinfectant. Our findings offer directions for future research to develop a more robust understanding of how virus-bacteria interactions in drinking water distribution systems affect water quality.


title: Vpr shapes the proviral landscape and polyclonal HIV-1 reactivation patterns in cultured cells

abstract:
Cell culture models suggest that the HIV-1 viral protein R (Vpr) is dispensable for latency establishment. However, whether Vpr affects the persistent proviral landscape and responsiveness to latency reversing agents (LRAs) is unclear. Here, integration site landscape, clonal dynamics, and latency reversal effects of Vpr were studied by comparing barcoded vpr+ and vpr- populations arising after infection of Jurkat cells in vitro. The results showed that individual integrant clones differed in fractions of LTR-active daughter cells: some clones gave rise to few to no LTR-active cells while for others almost all daughter cells were LTR-active. Integrant clones with at least 60% LTR-active cells (high LTR-active clones) contained proviruses positioned closer to preexisting enhancers (H3K27ac) and promoters (H3K4me3) than clones with <30% LTR-active cells (low LTR-active clones). Comparing vpr+ and vpr- populations revealed that the vpr+ population was depleted of high LTR-active clones. Complementing vpr-defective proviruses by transduction with vpr 16 days after infection led to rapid loss of high LTR-active clones, indicating that the effect of Vpr on proviral populations occurs post-integration. Comparing vpr+ and vpr- integration sites revealed that predominant vpr+ proviruses were farther from enhancers and promoters. Correspondingly, distances to these marks among previously reported intact HIV proviruses in ART-suppressed patients were more similar to those in the vpr+ pool than to vpr- integrants. To compare latency reactivation agent (LRA) responsiveness, the LRAs prostratin and JQ1 were applied separately or in combination. vpr+ and vpr- population-wide trends were similar, but combination treatment reduced virion release in a subset of vpr- clones relative to when JQ1 was applied separately, an effect not observed in vpr+ pools. Together, these observations highlight the importance of Vpr to proviral population dynamics, integration site landscapes, and responsiveness to latency reversing agents.

One Sentence Summary

Expression properties and responsiveness to latency reactivation agents of individual HIV-1 proviral clones within polyclonal populations are masked by dominant clones and influenced by proviral proximity to certain epigenetic marks and by Vpr, a viral factor not previously known to affect latency and reactivation.


title: 3D imaging for the quantification of spatial patterns in microbiota of the intestinal mucosa

abstract:
Improving our understanding of host-microbe relationships in the gut requires the ability to both visualize and quantify the spatial organization of microbial communities in their native orientation with the host tissue. We developed a systematic procedure to quantify the 3D spatial structure of the native mucosal microbiota in any part of the intestines with taxonomic and high spatial resolution. We performed a 3D biogeographical analysis of the microbiota of mouse cecal crypts at different stages of antibiotic exposure. By tracking eubacteria and four dominant bacterial taxa, we found that the colonization of crypts by native bacteria is a dynamic and spatially organized process. Ciprofloxacin treatment drastically reduced bacterial loads and eliminated Muribaculaceae (or all Bacteroidetes entirely) even 10 days after recovery when overall bacterial loads returned to pre-antibiotic levels. Our 3D quantitative imaging approach revealed that the bacterial colonization of crypts is organized in a spatial pattern that consists of clusters of adjacent colonized crypts that are surrounded by unoccupied crypts, and that this spatial pattern was resistant to the elimination of Muribaculaceae or of all Bacteroidetes by ciprofloxacin. Our approach also revealed that the composition of cecal crypt communities is diverse and that bacterial taxa are distributed differently within crypts, with Lactobacilli laying closer to the lumen than Bacteroidetes, Ruminococcaceae, and Lachnospiraceae. Finally, we found that crypts communities with similar taxonomic composition were physically closer to each other than communities that were taxonomically different.

Significance StatementMany human diseases are causally linked to the gut microbiota, yet the field still lacks mechanistic understanding of the underlying complex interactions because existing tools cannot simultaneously quantify microbial communities and their intact native context. In this work, we provide a new approach to tissue clearing and preservation that enables visualization, in 3D and at scales ranging from centimeters to micrometers, of the complete biogeography of the host-microbiota interface. We combine this new tool with sequencing and multiplexed labelling of the microbiota to provide the field with a platform on which to discover patterns in the spatial distribution of microbes. We validated this platform by quantifying the distribution of bacteria in the cecal mucosa at different stages of antibiotic exposure. This approach will enable researchers to formulate and test new hypotheses about host-microbe and microbe-microbe interactions.


title: Growth and development of trabecular structure in the calcaneus of Japanese macaques (Macaca fuscata) reflects locomotor behavior, life-history, and neuromuscular development.

abstract: We aim to broaden the analysis of bone structure by suggesting a new way to incorporate the interactions between behavior, neuromuscular development, and life-history. We examine the associations between these variables and age-related variation in trabecular structure in the calcaneus of Japanese macaques (Macaca fuscata). If skeletal markers linking these variables can be established, our inferences of the biology and behavior of fossil species would be significantly improved. We {micro}CT scanned the calcaneus in a cross-sectional sample of 36 juveniles aged between 0 and 7 years old and 5 adults at the Primate Research Institute, Japan. We calculated whole bone averages of standard trabecular properties and generated whole-bone morphometric maps of bone volume fraction and Youngs modulus. Trabecular structure is increasingly heterogeneous in older individuals. BV/TV decreases during the first month of life and increases afterwards, coinciding with the onset of independent locomotion. At birth, primary Youngs modulus is oriented orthogonal to the ossification center, but after locomotor onset bone structure becomes stiffest in the direction of joint surfaces and muscle attachments. Age-related variation in bone volume fraction is best predicted by an interaction between neuromaturation, body mass, and locomotor independence. Results support the common assumption that trabecular structure dynamically adapts to novel joint loading conditions during ontogeny. The timing of independent locomotion, body size, and neuromuscular development, are all correlated to age-related variation in the trabecular structure of the macaque calcaneus. The causal mechanisms behind the observed patterns cannot be directly inferred from our cross-sectional study. If the model presented in this paper holds up under longitudinal experimental conditions, trabecular structure can be used both to infer behavior from fossil morphology and to serve as a valuable proxy for neuromuscular maturation and life history events like locomotor onset and the achievement of an adult-like gait.

category: developmental biology


title: Can DNA help trace the local trade of pangolins? A genetic assessment of white-bellied pangolins from the Dahomey Gap (West Africa)

abstract: We conducted in the Dahomey Gap (DG) a pioneer study on the genetic tracing of the African pangolin trade. We sequenced and genotyped 189 white-bellied pangolins from 18 forests and 12 wildlife markets using one mitochondrial fragment and 20 microsatellites loci. Tree-based assignment procedure showed the  endemicity of the pangolin trade, as strictly fed by the lineage endemic to the DG (DGL). DGL populations were characterized by low levels of genetic diversity, an overall absence of equilibrium, inbreeding depression and lack of geographic structure. We identified a 92-98% decline in DGL effective population size 200-500 ya -concomitant with major political transformations along the  Slave Coast- leading to contemporaneous estimates inferior to minimum viable population size. Genetic tracing suggested that wildlife markets from the DG sourced through the entire DGL range. Our loci provided the necessary power to distinguish among all the genotyped pangolins, tracing the dispatch of same individuals on the markets and within local communities. We developed an approach combining rarefaction analysis of private allele frequencies and cross-validation with observed data that could trace five traded pangolins to their forest origin, c. 200-300 km away from the markets. Although the genetic toolkit that we designed from traditional markers can prove helpful to trace the pangolin trade, our tracing ability was limited by the lack of population structure within DGL. Given the deleterious combination of genetic, demographic and trade-related factors affecting DGL populations, the conservation status of white-bellied pangolins in the DG should be urgently re-evaluated.


title: Polyunsaturated fatty acids inhibit a pentameric ligand-gated ion channel through one of two binding sites

abstract: Polyunsaturated fatty acids (PUFAs) inhibit pentameric ligand-gated ion channels (pLGICs) but the mechanism of inhibition is not well understood. The PUFA, docosahexaenoic acid (DHA), inhibits agonist responses of the pLGIC, ELIC, more effectively than palmitic acid, similar to the effects observed in the GABAA receptor and nicotinic acetylcholine receptor. Using photo-affinity labeling and coarse-grained molecular dynamics simulations, we identified two fatty acid binding sites in the outer transmembrane domain (TMD) of ELIC. Fatty acid binding to the photolabeled sites is selective for DHA over palmitic acid, and specific for an agonist-bound state. Hexadecyl-methanethiosulfonate modification of one of the two fatty acid binding sites in the outer TMD recapitulates the inhibitory effect of PUFAs in ELIC. The results demonstrate that DHA selectively binds to multiple sites in the outer TMD of ELIC, but that state-dependent binding to a single intrasubunit site mediates DHA inhibition of ELIC.


title: Automated Optimization of Residual Reduction Algorithm Parameters in OpenSim

abstract: The residual reduction algorithm (RRA) in OpenSim is designed to improve dynamic consistency of kinematics and ground reaction forces in movement simulations of musculoskeletal models. RRA requires the user to select numerous tracking weights for the joint kinematics to reduce residual errors. Selection is often performed manually, which can be time-consuming and is unlikely to yield optimal tracking weights. A multi-heuristic optimization algorithm was used to expedite tracking weight decision making to reduce residual errors. This method produced more rigorous results than manual iterations and although the total computation time was not significantly reduced, this method does not require the user to monitor the algorithms progress to find a solution, thereby reducing manual tuning. Supporting documentation and code to implement this optimization is freely provided to assist the community with developing movement simulations.

category: bioengineering


title: Pseudofinder: detection of pseudogenes in prokaryotic genomes

abstract: Prokaryotic genomes are generally gene dense and encode relatively few pseudogenes, or nonfunctional/inactivated remnants of genes. However, in certain contexts, such as recent ecological shifts or extreme population bottlenecks (such as those experienced by symbionts and pathogens), pseudogenes can quickly accumulate and form a substantial fraction of the genome. Identification of pseudogenes is, thus, a critical step for understanding the evolutionary forces acting upon, and the functional potential encoded within, prokaryotic genomes. Here, we present Pseudofinder, an open-source software dedicated to pseudogene identification and analysis. With Pseudofinders multi-pronged, reference-based approach, we demonstrate its capacity to detect a wide variety of pseudogenes, including those that are highly degraded and typically missed by gene-calling pipelines, as well newly formed pseudogenes, which can have only one or a few inactivating mutations. Additionally, Pseudofinder can detect intact genes undergoing relaxed selection, which may indicate incipient pseudogene formation. Implementation of Pseudofinder in annotation pipelines will not only clarify the functional potential of sequenced microbes, but will also generate novel insights and hypotheses regarding the evolutionary dynamics of bacterial and archaeal genomes.


title: Kinetic interplay between droplet maturation and coalescence modulates shape of aged protein condensates

abstract: Biomolecular condensates formed by the process of liquid-liquid phase separation (LLPS) play diverse roles inside cells, from spatiotemporal compartmentalisation to speeding up chemical reactions. Upon maturation, the liquid-like properties of condensates, which underpin their functions, are gradually lost, eventually giving rise to solid-like states with potential pathological implications. Enhancement of inter-protein interactions is one of the main mechanisms suggested to trigger the formation of solid-like condensates. To gain a molecular-level understanding of how the accumulation of stronger interactions among proteins inside condensates affect the kinetic and thermodynamic properties of biomolecular condensates, and their shapes over time, we develop a tailored coarse-grained model of proteins that transition from establishing weak to stronger inter-protein interactions inside condensates. Our simulations reveal that the fast accumulation of strongly binding proteins during the nucleation and growth stages of condensate formation results in aspherical solid-like condensates. In contrast, when strong inter-protein interactions appear only after the equilibrium condensate has been formed, or when they accumulate slowly over time, with respect to the time needed for droplets to fuse and grow, spherical solid-like droplets emerge. By conducting atomistic potential-of-mean-force simulations of NUP-98 peptides--prone to forming inter-protein {beta}-sheets--we observe that formation of inter-peptide {beta}-sheets increases the strength of the interactions consistently with the loss of liquid-like condensate properties we observe at the coarse-grained level. Overall, our work aids in elucidating fundamental molecular, kinetic, and thermodynamic mechanisms linking the rate of change in protein interaction strength to condensate shape and maturation during ageing.


title: Effects of cryo-EM cooling on structural ensembles

abstract: Structure determination by cryo electron microscopy (cryo-EM) provides information on structural heterogeneity and ensembles at atomic resolution. To obtain cryo-EM images of macromolecules, the samples are first rapidly cooled down to cryogenic temperatures. To what extent the structural ensemble is perturbed during cooling is currently unknown. Here, to quantify the effects of cooling, we combined continuum model calculations of the temperature drop, molecular dynamics simulations of a ribosome complex before and during cooling with kinetic models. Our results suggest that three effects markedly contribute to the narrowing of the structural ensembles: thermal contraction, reduced thermal motion within local potential wells, and the equilibration into lower free-energy conformations by overcoming separating free-energy barriers. During cooling, barrier heights below 10 kJ/mol were found to be over-come, which is expected to reduce B-factors in ensembles imaged by cryo-EM. Our approach now enables the quantification of the heterogeneity of room-temperature ensembles from cryo-EM structures.


title: A mechanobiological model for tumor spheroids evolution: application to glioblastoma

abstract:
Spheroids are in vitro spherical structures of cell aggregates, eventually cultured within a hydrogel matrix, that are used, among other applications, as a technological platform to investigate tumor formation and evolution. Several interesting features can be replicated using this methodology, such as cell communication mechanisms, the effect of gradients of nutrients, or the creation of realistic 3D biological structures. In this paper, we propose a continuum mechanobiological model which accounts for the most relevant phenomena that take place in tumor spheroids evolution under in vitro suspension, namely, nutrients diffusion in the spheroid, kinetics of cellular growth and death, and mechanical interactions among the cells. The model is qualitatively validated, after calibration of the model parameters, versus in vitro experiments of spheroids of different glioblastoma cell lines. This preliminary validation allowed us to conclude that glioblastoma tumor spheroids evolution is mainly driven by mechanical interactions of the cell aggregate and the dynamical evolution of the cell population. In particular, it is concluded that our model is able to explain quite different setups, such as spheroids growth (up to six times the initial configuration for U-87 MG cell line) or shrinking (almost half of the initial configuration for U-251 MG cell line); as the result of the mechanical interplay of cells driven by cellular evolution. Indeed, the main contribution of this work is to link the spheroid evolution with the mechanical activity of cells, coupled with nutrient consumption and the subsequent cell dynamics. All this information can be used to further investigate mechanistic effects in the evolution of tumors and their role in cancer disease.

Author summarySpheroids structures of cell aggregates are an available experimental platform to analyze the evolution and drug response of solid tumors. In particular, the dynamics of different glioblastoma cell lines have been studied in this work using spheroids. Interestingly, very different behaviors were observed, from a half of the initial configuration shrinking for U-251 MG cell line to six times the initial configuration growth for U-87 MG cell line. These results were replicated by means of a coupled mathematical model which accounts for nutrients diffusion in the spheroid, kinetics of cellular growth and death, and mechanical interactions among the cells. Tumor growth or shrinkage can be explained from a continuum mechanics view driven by cell activity and nutrients availability. This modeling put the focus on mechanistic effects and is aligned with novel experimental techniques to quantify the mechanical microenvironment in tumors. These techniques may be combined with the approach presented in this work to further investigate the role of mechanics in cancer disease.


title: High Expression of Glycolytic Genes in Clinical Glioblastoma Patients Correlates with Lower Survival

abstract: Glioblastoma (GBM), the most aggressive brain tumor, is associated with a median survival at diagnosis of 16-20 months and limited treatment options. The key hallmark of GBM is altered tumor metabolism and marked increase in the rate of glycolysis. Aerobic glycolysis along with elevated glucose consumption and lactate production supports rapid cell proliferation and GBM growth. In this study, we examined the gene expression profile of metabolic targets in GBM samples from patients with low grade glioma (LGG) and GBM. We found that gene expression of glycolytic enzymes is up-regulated in GBM samples and significantly associated with an elevated risk for developing GBM. Our findings of clinical outcomes showed that GBM patients with high expression of HK2 and PKM2 in the glycolysis related genes and low expression of genes involved in mitochondrial metabolism-SDHB and COX5A related to tricarboxylic acid (TCA) cycle and oxidative phosphorylation (OXPHOS), respectively, was associated with poor patient overall survival. Surprisingly, expression levels of genes involved in mitochondrial oxidative metabolism are markedly increased in GBM compared to LGG but was lower compared to normal brain. The fact that in GBM the expression levels of TCA cycle and OXPHOS-related genes are higher than those in LGG patients suggests the metabolic shift in GBM cells when progressing from LGG to GBM. These results are an important step forward in our understanding of the role of metabolic reprogramming in glioma as drivers of the tumor and could be potential prognostic targets in GBM therapies.


title: Single-cell RNA transcriptome analysis reveals CXCL16/CXCR6 as maintenance factors for CNS tissue resident T cells that drive synapse elimination during viral recovery

abstract:
BackgroundEmerging RNA viruses that target the central nervous system (CNS) lead to cognitive sequelae in survivors. Studies in humans and mice infected with West Nile virus (WNV), a re-emerging RNA virus associated with learning and memory deficits, revealed microglial-mediated synapse elimination within the hippocampus. Moreover, CNS resident memory T (TRM) cells activate microglia, limiting synapse recovery and inducing spatial learning defects in WNV-recovered mice. The signals involved in T cell-microglia interactions are unknown.

MethodsHere, we examined the murine WNV-recovered forebrain using single-cell RNA sequencing to identify putative ligand-receptor pairs involved in intercellular communication between T cells and microglia. Clustering and differential gene analyses were followed by protein validation, genetic and antibody-based approaches utilizing an established murine model of WNV recovery in which microglia and complement promote ongoing hippocampal synaptic loss.

ResultsProfiling of host transcriptome at 25 days post-infection revealed a shift in forebrain homeostatic microglia to activated subpopulations with transcriptional signatures that have previously been observed in studies of neurodegenerative diseases. Importantly, CXCL16/CXCR6, a chemokine signaling pathway involved in TRM cell biology, was identified as critically regulating CXCR6 expressing CD8+ TRM cell numbers within the WNV-recovered forebrain. We demonstrate that CXCL16 is highly expressed by all myeloid cells, and its unique receptor, CXCR6, is highly expressed on all CD8+ T cells. Using genetic and pharmacological approaches, we demonstrate that CXCL16/CXCR6 is required not only for the maintenance of WNV-specific, CD8 TRM cells in the post-infectious CNS, but also contributes to their expression of TRM cell markers. Moreover, CXCR6+CD8+ T cells are required for glial activation and ongoing synapse elimination.

ConclusionsWe provide a comprehensive assessment of the role of CXCL16/CXCR6 as an interaction link between microglia and CD8+ T cells that maintains forebrain TRM cells, microglial and astrocyte activation, and ongoing synapse elimination in virally recovered animals. We also show that therapeutic targeting of CXCL16 during recovery may reduce CNS CD8+ TRM cells.


title: High-throughput mutagenesis identifies mutations and RNA binding proteins controlling CD19 splicing and CART-19 therapy resistance

abstract:
During CART-19 immunotherapy for B-cell acute lymphoblastic leukaemia (B-ALL), many patients relapse due to loss of the cognate CD19 epitope. Since epitope loss can be caused by aberrant CD19 exon 2 processing, we herein investigate the regulatory code that controls CD19 splicing. We combine high-throughput mutagenesis with mathematical modelling to quantitatively disentangle the effects of all mutations in the region comprising CD19 exons 1-3. Thereupon, we identify ~200 single point mutations that alter CD19 splicing and thus could predispose B-ALL patients to CART-19 resistance. Furthermore, we report almost 100 previously unknown splice isoforms that emerge from cryptic splice sites and likely encode non-functional CD19 proteins. We further identify cis-regulatory elements and trans-acting RNA-binding proteins that control CD19 splicing (e.g., PTBP1 and SF3B4) and validate that loss of these factors leads to enhanced CD19 mis-splicing. Our dataset represents a comprehensive resource for potential prognostic factors predicting success of CART-19 therapy.

HighlightsO_LIMutations in relapsed CART-19 patients lead to CD19 mis-splicing
C_LIO_LIHigh-throughput mutagenesis uncovers ~200 single point mutations with a potential role in CART-19 therapy resistance
C_LIO_LIMany mutations generate non-functional CD19 proteins by activating cryptic splice sites
C_LIO_LIRNA-binding proteins such as PTBP1 are key to the expression of properly spliced, CART-19 immunotherapy-sensitive isoforms
C_LI


title: WUSCHEL-Related Homeobox Genes Cooperate with Cytokinin Signalling to Promote Bulbil Formation in Lilium lancifolium

abstract: The bulbil is an important vegetative reproductive organ in triploid Lilium lancifolium. Based on our previously obtained transcriptome data, we screened two WUSCHCEL-related homeobox (WOX) genes closely related to bulbil formation, LlWOX9 and LlWOX11. However, the biological functions and regulatory mechanisms of LlWOX9 and LlWOX11 are unclear. In this study, we cloned the full-length coding sequences of LlWOX9 and LlWOX11. Transgenic Arabidopsis showed increased branch numbers, and the overexpression of LlWOX9 and LlWOX11 in stem segments promoted bulbil formation, while the silencing of LlWOX9 and LlWOX11 inhibited bulbil formation, indicating that LlWOX9 and LlWOX11 are positive regulators of bulbil formation. Cytokinins acting through type-B response regulators (type-B RRs) could bind to the promoters of LlWOX9 and LlWOX11 and promote their transcription. LlWOX11 could enhance cytokinin pathway signalling by inhibiting the transcription of type-A LlRR9. Our study enriches the understanding of the regulation of plant development by the WOX gene family and lays a foundation for further research on the molecular mechanism of bulbil formation in lily.


title: Investigation of vitamin B12 concentrations and tissue distributions in larval and adult Pacific oysters and related bivalves

abstract: Vitamin B12 (B12) is an essential micronutrient for all animals, but is not present in plants and is produced de novo only by bacteria or archaea. Accordingly, humans must derive required B12 from eating animal products or vitamin supplements, as deficiencies can lead to severe health issues including neuropathy. An often overlooked source in the human diet of B12 is shellfish, in particular bivalves, which have significantly higher levels of B12 than other animal sources, including all vertebrate meats. Origins and key metabolic processes involving B12 in bivalves remain largely unknown, despite the exceptionally high levels. In this study, we examined in several Australian bivalve species, hypotheses concerning B12 utilisation and uptake through diet or microorganism symbiosis. Vitamin B12 is not distributed evenly across different tissues types of the Pacific oyster, the commercial scallop and Goolwa cockle (pipi), with higher accumulation in the oyster adductor muscle and gill, and mantle and syphons of the Goolwa cockle. Oyster larvae before first feeding already contained high amount of B12; however, a significant decrease in B12 concentration post metamorphosis indicates a higher utilisation of B12 during this life event. We demonstrated that microalgal feed can be supplemented with B12, resulting in an enriched feed, but this did not result in an increase in larval B12 concentrations when oyster larvae were fed with this diet relative to controls, thus supporting the theory that a B12 producing microbiome within bivalves was the potential source of B12 rather than feed. However, B12 concentrations in the digestive tract of adult oysters were low compared to other tissue types, which might challenge this theory, at least in adults. Our findings provide insight into B12 uptake and function in bivalve species, which will aid the promotion of bivalves as suitable B12 source for humans as well as provide crucial information to the aquaculture industry in relation to optimisation of vitamin supplementation in bivalve hatchery production.


title: Coevolution of interacting proteins through non-contacting and non-specific mutations

abstract: Proteins often accumulate neutral mutations that do not affect current functions1 but can profoundly influence future mutational possibilities and functions2-4. Understanding such hidden potential has major implications for protein design and evolutionary forecasting5-7, but has been limited by a lack of systematic efforts to identify potentiating mutations8,9. Here, through the comprehensive analysis of a bacterial toxin-antitoxin system, we identified all possible single substitutions in the toxin that enable it to tolerate otherwise interface-disrupting mutations in its antitoxin. Strikingly, the majority of enabling mutations in the toxin do not contact, and promote tolerance non-specifically to, many different antitoxin mutations, despite covariation in homologs occurring primarily between specific pairs of contacting residues across the interface. In addition, the enabling mutations we identified expand future mutational paths that both maintain old toxin-antitoxin interactions and form new ones. These non-specific mutations are missed by widely used covariation and machine learning methods10,11. Identifying such enabling mutations will be critical for ensuring continued binding of therapeutically relevant proteins, such as antibodies, aimed at evolving targets12-14.

category: evolutionary biology


title: VIP interneurons in mouse whisker S1 exhibit sensory and action-related signals during goal-directed behavior

abstract: Vasoactive intestinal peptide-expressing (VIP) interneurons, which constitute 10-15% of the cortical inhibitory neuron population1,2, have emerged as an important cell type for regulating excitatory cell activity based on behavioral state. VIP cells in sensory cortex are potently engaged by neuromodulatory and motor inputs during active exploratory behaviors like locomotion and whisking, which in turn promote pyramidal cell firing via disinhibition3-5. Such state-dependent modulation of activity by VIP cells in sensory cortex has been studied widely in recent years. However, the function of VIP cells during goal-directed behavior is less well understood. It is not clear how task-related events like sensory stimuli, motor actions, or reward activate VIP cells in sensory cortex since there is often temporal overlap in the occurrence of these events. We developed a Go/NoGo whisker touch detection task which incorporates a post-stimulus delay period to separate sensory-driven activity from action- or reward-related activity during behavior. We used 2-photon calcium imaging to measure task-related signals of L2/3 VIP neurons in S1 of behaving mice. We report for the first time that VIP cells in mouse whisker S1 are activated by both whisker stimuli and goal-directed licking. Whisker- and lick-related signals were spatially organized in relation to anatomical columns in S1. Sensory responses of VIP cells were tuned to specific whiskers, whether or not they also displayed lick-related activity.