title: DomBpred: protein domain boundary predictor using inter-residue distance and domain-residue level clustering

abstract: Domain boundary prediction is one of the most important problems in the study of protein structure and function, especially for large proteins. At present, most domain boundary prediction methods have low accuracy and limitations in dealing with multi-domain proteins. In this study, we develop a sequence-based protein domain boundary predictor, named DomBpred. In DomBpred, the input sequence is firstly classified as either a single-domain protein or a multi-domain protein through a designed effective sequence metric based on a constructed single-domain sequence library. For the multi-domain protein, a domain-residue level clustering algorithm inspired by Ising model is proposed to cluster the spatially close residues according inter-residue distance. The unclassified residues and the residues at the edge of the cluster are then tuned by the secondary structure to form potential cut points. Finally, a domain boundary scoring function is proposed to recursively evaluate the potential cut points to generate the domain boundary. DomBpred is tested on a large-scale test set of FUpred comprising 2549 proteins. Experimental results show that DomBpred better performs than the state-of-the-art methods in classifying whether protein sequences are composed by single or multiple domains, and the Matthews correlation coefficient is 0.882. Moreover, on 849 multi-domain proteins, the domain boundary distance and normalised domain overlap scores of DomBpred are 0.523 and 0.824, respectively, which are 5.0% and 4.2% higher than those of the best comparison method, respectively. Comparison with other methods on the given test set shows that DomBpred outperforms most state-of-the-art sequence-based methods and even achieves better results than the top-level template-based method.


title: Glycosaminoglycans and glycoproteins influence the elastic response of synovial fluid nanofilms on model oxide surfaces

abstract: Synovial fluid (SF) is the natural lubricant found in articulated joints, providing unique cartilage surface protecting films under confinement and relative motion. While it is known that the synergistic interactions of the macromolecular constituents provide its unique load-bearing and tribological performance, it is not fully understood how two of the main constituents, glycosaminoglycans (GAGs) and glycoproteins, regulate the formation and mechanics of robust load-bearing films. Here, we present evidence that the load-bearing capabilities, rather than the tribological performance, of the formed SF films depend strongly on its components integrity. For this purpose, we used a combination of enzymatic treatments, quartz crystal microbalance with dissipation (QCM-D) and the surface forces apparatus (SFA) to characterize the formation and load-bearing capabilities of SF films on model oxide (i.e., silicates) surfaces. We find that, upon cleavage of proteins, the elasticity of the films is reduced and that cleaving GAGs results in irreversible (plastic) molecular re-arrangements of the film constituents when subjected to confinement. Understanding thin film mechanics of SF can provide insight into the progression of diseases, such as arthritis, but may also be applicable to the development of new implant surface treatments or new biomimetic lubricants.

category: biophysics


title: Adaptive ratchets and the evolution of molecular complexity

abstract: Biological systems have evolved to amazingly complex states, yet we do not understand in general how evolution operates to generate increasing genetic and functional complexity. Molecular recognition sites are short genome segments or peptides binding a cognate recognition target of sufficient sequence similarity. Such sites are simple, ubiquitous modules of sequence information, cellular function, and evolution. Here we show that recognition sites, if coupled to a time-dependent target, can rapidly evolve to complex states with larger code length and smaller coding density than sites recognising a static target. The underlying fitness model contains selection for recognition, which depends on the sequence similarity between site and target, and a uniform cost per unit of code length. Site sequences are shown to evolve in a specific adaptive ratchet, which produces selection of different strength for code extensions and compressions. Ratchet evolution increases the adaptive width of evolved sites, accelerating the adaptation to moving targets and facilitating refinement and innovation of recognition functions. We apply these results to the recognition of fast-evolving antigens by the human immune system. Our analysis shows how molecular complexity can evolve as a collateral to selection for function in a dynamic environment.


title: Optimizing reaction coordinate by flux maximization in the transition path ensemble

abstract: Transition path ensemble is a collection of reactive trajectories, all of which largely keep going forward along the transition channel from the reactant state to the product one, and is believed to possess the information necessary for the identification of reaction coordinate. Previously, the full coordinates (both position and momentum) of the snapshots in the transition path ensemble were utilized to obtain the reaction coordinate (J. Chem. Phys. 2016, 144, 114103; J. Chem. Phys. 2018, 148, 084105). Here, with the conformational (or position) coordinates alone, it is demonstrated that the reaction coordinate can be optimized by maximizing the flux of a given coordinate in the transition path ensemble. In the application to alanine dipeptide in vacuum, dihderal angles{phi} and{theta} were identified to be the two best reaction coordinates, which was consistent with the results in existing studies. A linear combination of these two coordinates gave a better reaction coordinate, which is highly correlated with committor. Most importantly, the method obtained a linear combination of pairwise distances between heavy atoms, which was highly correlated with committor as well. The standard deviation of committor at the transition region defined by the optimized reaction coordinate is as small as 0.08. In addition, the effects of practical factors, such as the choice of transition path sub-ensembles and saving interval between frames in transition paths, on reaction coordinate optimization were also considered.


title: Model for ring closure in ER tubular network dynamics

abstract: Tubular networks of endoplasmic reticulum (ER) are dynamic structures whose steady-state conformations are maintained by a dynamic balance between the persistent generation and vanishing of the network elements. While factors producing the ER tubules and inter-tubular junctions have been investigated, the mechanisms behind their elimination remained unknown. Here we addressed the ER ring closure, the process resulting in the tubule and junction removal through constriction of the network unit-cells into junctional knots followed by the knot remodeling into regular junctions. We considered the ring closure to be driven by the tension existing in ER membranes. We modeled, computationally, the structures of the junctional knots containing internal nanopores and analyzed their tension dependence. We predicted an effective interaction between the nanopores facilitating the knot tightening and collapse of additional network unit cells. We analyzed the process of the pore sealing through membrane fission resulting in formation of regular junctions. Considering the hemi-fission as the rate-limiting stage of the fission reaction, we evaluated the membrane tensions guarantying the spontaneous character of the pore sealing. We concluded that feasible membrane tensions explain all stages of the ER ring closure.


title: Simultaneous multicolor DNA-PAINT without sequential fluid exchange using spectral demixing

abstract: Several variants of multicolor single-molecule localization microscopy (SMLM) have been developed to resolve the spatial relationship of nanoscale structures in biological samples. The oligonucleotide-based SMLM approach  DNA-PAINT robustly achieves nanometer localization precision and can be used to count binding sites within nanostructures. However, multicolor DNA-PAINT has primarily been realized by  Exchange-PAINT that requires sequential exchange of the imaging solution and thus leads to extended acquisition times. To alleviate the need for fluid exchange and to speed up the acquisition of current multichannel DNA-PAINT, we here present a novel approach that combines DNA-PAINT with simultaneous multicolor acquisition using spectral demixing (SD). By using newly designed probes and a novel multichannel registration procedure we achieve simultaneous multicolor SD-DNA-PAINT with minimal crosstalk. We demonstrate high localization precision (3 - 6 nm) and multicolor registration of dual and triple-color SD-DNA-PAINT by resolving patterns on DNA origami nanostructures and cellular structures.


title: Cell Painting predicts impact of lung cancer variants

abstract: Most variants in most genes across most organisms have an unknown impact on the function of the corresponding gene. This gap in knowledge is especially acute in cancer, where clinical sequencing of tumors now routinely reveals patient-specific variants whose functional impact on the corresponding gene is unknown, impeding clinical utility. Transcriptional profiling was able to systematically distinguish these variants of unknown significance (VUS) as impactful vs. neutral in an approach called expression-based variant-impact phenotyping (eVIP). We profiled a set of lung adenocarcinoma-associated somatic variants using Cell Painting, a morphological profiling assay that captures features of cells based on microscopy using six stains of cell and organelle components. Using deep-learning-extracted features from each cells image, we found that cell morphological profiling (cmVIP) can predict variants functional impact and, particularly at the single-cell level, reveals biological insights into variants which can be explored in our public online portal. Given its low cost, convenient implementation, and single-cell resolution, cmVIP profiling therefore seems promising as an avenue for using non-gene-specific assays to systematically assess the impact of variants, including disease-associated alleles, on gene function.

category: cancer biology


title: Label-free imaging of collagen fibers in tissue slices using phase imaging with computational specificity

abstract: Evaluating the tissue collagen content in addition to the epithelial morphology has been proven to offer complementary information in histopathology, especially in oncology tumor staging and prediction of survival in cancer patients. One imaging modality widely used for this purpose is second harmonic generation microscopy (SHGM), which reports on the nonlinear susceptibility associated with the collagen fibers. Another method is polarization light microscopy (PLM) combined with picrosirius-red (PSR) tissue staining. However, SHGM requires expensive equipment and provides limited throughput, while PLM and PSR staining are not part of the routine surgical pathology workflow. Here, we utilize phase imaging with computational specificity (PICS) to computationally infer the collagen distribution of unlabeled tissue, with high specificity. PICS utilizes deep learning to translate quantitative phase images (QPI) into corresponding PSR images with high accuracy and inference speed of 200 milisecond per forwardpass through the model once trained. We developed a multimodal imaging instrument that yields both Spatial light Inference Microscopy (SLIM) and polarized light microscopy (PLM) images from the same field of view. Our results indicate that the distributions of collagen fiber orientation, length, and straightness reported by PICS closely match the ones from ground truth as defined by KL-divergence.

category: pathology


title: Proteolysis Targeting Chimeras With Reduced Off-targets

abstract: Proteolysis Targeting Chimeras (PROTACs), a class of heterobifunctional molecules that recruit target proteins to E3 ligases, have gained traction for targeted protein degradation. However, pomalidomide, a widely used E3 ligase recruiter in PROTACs, can independently degrade other targets, such as zinc-finger (ZF) proteins, that hold key functions in normal development and disease progression. This off-target degradation of pomalidomide-based PROTACs raises concerns about their therapeutic applicability and long-term side effects. Therefore, there is a crucial need to develop rules for PROTAC design that minimize off-target degradation. In this study, we developed a high-throughput platform that interrogates the off-target degradation of ZF domains and discovered, using this platform, that PROTACs with the current design paradigm induce significant degradation of several ZF proteins. To identify new rules for PROTAC design, we generated a rationalized library of pomalidomide analogs with distinct exit vector modifications on the C4 and C5 positions of the phthalimide ring and profiled their propensities for ZF protein degradation. We found that modifications on the C5 position with nucleophilic aromatic substitution (SNAr) reduce off-target ZF degradation. We applied our newfound design principles on a previously developed ALK oncoprotein-targeting PROTAC and generated PROTACs with enhanced potency and minimal off-target degradation. We envision the reported off-target profiling platform and pomalidomide analogs will find utility in design of specific PROTACs.

category: pharmacology and toxicology


title: Glyphosate and its formulations Roundup Bioflow and RangerPro alter bacterial and fungal community composition in the rat caecum microbiome

abstract: The potential health consequences of glyphosate-induced gut microbiome alterations have become a matter of intense debate. As part of a multifaceted study investigating toxicity, carcinogenicity and multigenerational effects of glyphosate and its commercial herbicide formulations, we assessed changes in bacterial and fungal populations in the caecum microbiota of rats exposed prenatally until adulthood (13 weeks after weaning) to three doses of glyphosate (0.5, 5, 50 mg/kg body weight/day), or to the formulated herbicide products Roundup Bioflow and RangerPro at the same glyphosate-equivalent doses. Caecum bacterial microbiota were evaluated by 16S rRNA sequencing whilst the fungal population was determined by ITS2 amplicon sequencing. Results showed that both fungal and bacterial diversity were affected by the Roundup formulations in a dose-dependent manner, whilst glyphosate alone significantly altered only bacterial diversity. At taxa level, a reduction in Bacteroidota abundance, marked by alterations in the levels of Alloprevotella, Prevotella and Prevotellaceae UCG-003, was concomitant to increased levels of Firmicutes (e.g., Romboutsia, Dubosiella, Eubacterium brachy group or Christensenellaceae) and Actinobacteria (e.g., Enterorhabdus, Adlercreutzia, or Asaccharobacter). Treponema and Mycoplasma also had their levels reduced by the pesticide treatments. Analysis of fungal composition indicated that the abundance of the rat gut commensal Ascomycota Kazachstania was reduced while the abundance of Gibberella, Penicillium, Claviceps, Cornuvesica, Candida, Trichoderma and Sarocladium were increased by exposure to the Roundup formulations, but not to glyphosate. Altogether, our data suggest that glyphosate and its Roundup RangerPro and Bioflow caused profound changes in caecum microbiome composition by affecting the fitness of major commensals, which in turn reduced competition and allowed opportunistic fungi to grow in the gut, in particular in animals exposed to the herbicide formulations. This further indicates that changes in gut microbiome composition might influence the long-term toxicity, carcinogenicity and multigenerational effects of glyphosate-based herbicides.


title: Algorithmic Parameter Estimation and Uncertainty Quantification for Hodgkin-Huxley Neuron Models

abstract: Experimental data-based parameter search for Hodgkin-Huxley-style (HH) neuron models is a major challenge for neuroscientists and neuroengineers. Current search strategies are often computationally expensive, are slow to converge, have difficulty handling nonlinearities or multimodalities in the objective function, or require good initial parameter guesses. Most important, many existing approaches lack quantification of uncertainties in parameter estimates even though such uncertainties are of immense biological significance. We propose a novel method for parameter inference and uncertainty quantification in a Bayesian framework using the Markov chain Monte Carlo (MCMC) approach. This approach incorporates prior knowledge about model parameters (as probability distributions) and aims to map the prior to a posterior distribution of parameters informed by both the model and the data. Furthermore, using the adaptive parallel tempering strategy for MCMC, we tackle the highly nonlinear, noisy, and multimodal loss function, which depends on the HH neuron model. We tested the robustness of our approach using the voltage trace data generated from a 9-parameter HH model using five levels of injected currents (0.0, 0.1, 0.2, 0.3, and 0.4 nA). Each test consisted of running the ground truth with its respective currents to estimate the model parameters. To simulate the condition for fitting a frequency-current (F-I) curve, we also introduced an aggregate objective that runs MCMC against all five levels simultaneously. We found that MCMC was able to produce many solutions with acceptable loss values (e.g., for 0.0 nA, 889 solutions were within 0.5% of the best solution and 1,595 solutions within 1% of the best solution). Thus, an adaptive parallel tempering MCMC search provides a "landscape" of the possible parameter sets with acceptable loss values in a tractable manner. Our approach is able to obtain an intelligently sampled global view of the solution distributions within a search range in a single computation. Additionally, the advantage of uncertainty quantification allows for exploration of further solution spaces, which can serve to better inform future experiments.

category: neuroscience


title: Antigen presentation by type 3 innate lymphoid cells instructs the differentiation of gut microbiota-specific regulatory T cells

abstract: The mutualistic relationship of gut-resident microbiota and cells of the host immune system promotes homeostasis that ensures maintenance of the microbial community and of a poised, but largely non-aggressive, immune cell compartment1, 2. Consequences of disturbing this balance, by environmental or genetic factors, include proximal inflammatory conditions, like Crohns disease, and systemic illnesses, both metabolic and autoimmune. One of the means by which this equilibrium is achieved is through induction of both effector and suppressor or regulatory arms of the adaptive immune system. In mice, Helicobacter species induce regulatory (iTreg) and follicular helper (Tfh) T cells in the colon-draining mesenteric lymph nodes under homeostatic conditions, but can instead induce inflammatory Th17 cells when iTreg cells are compromised3, 4. How Helicobacter hepaticus and other gut bacteria direct T cells to adopt distinct functions remains poorly understood. Here, we investigated which cells and molecular components are required to convey the microbial instruction for the iTreg differentiation program. We found that antigen presentation by cells expressing ROR{gamma}t, rather than by classical dendritic cells, was both required and sufficient for iTreg induction. These ROR{gamma}t+ cells, likely to be type 3 innate lymphoid cells (ILC3) and/or a recently-described population of Aire+ cells termed Janus cells5, require the MHC class II antigen presentation machinery, the chemokine receptor CCR7, and v integrin, which activates TGF-{beta}, for iTreg cell differentiation. In the absence of any of these, instead of iTreg cells there was expansion of microbiota-specific pathogenic Th17 cells, which were induced by other antigen presenting cells (APCs) that did not require CCR7. Thus, intestinal commensal microbes and their products target multiple APCs with pre-determined features suited to directing appropriate T cell differentiation programs, rather than a common APC that they endow with appropriate functions. Our results illustrate the ability of microbiota to exploit specialized functions of distinct innate immune system cells, targeting them to achieve the desired composition of equipoised T cells, thus maintaining tolerance.

category: immunology


title: Microbiota-instructed regulatory T cell differentiation is mediated by a distinct RORγt+ antigen presenting cell subset


title: Insights into the trafficking of human leukocytes to colostrum evidences a modulation of the B lymphocyte compartment in obesity

abstract: Breastmilk is a dynamic fluid which initial goal is to provide the most adapted nutrition to the neonate. Additional functions have been recently attributed to breastmilk, with the evidence of a specific microbiota and the presence of a variety of components of the immune system, such as cytokines and leukocytes. The composition of breastmilk varies through time, according to the health status of mother and child, and altogether contributes to future health of the infant. Obesity is a rising condition worldwide, that creates a state of systemic, chronic inflammation including leukocytosis. Here, we asked whether colostrum, the milk produced within the first 48 h post-partum, would contain a distinct leukocyte composition depending on the body mass index (BMI) of the mother. We applied a panel of 6 antibodies plus viability marker to the peripheral blood and colostrum obtained from obese (BMI > 30) and lean (BMI < 25) mothers to characterize 10 major leukocyte subpopulations using flow cytometry. While lymphoid cells were otherwise unaffected by their tissue of origin, the phenotypes of granulocyte and monocyte populations significantly contrasted between blood and colostrum, including variations in morphology and surface expression of CD45 and CD16. These differences recapitulated across groups, which suggests a generalized cell-specific phenotype alteration caused by trafficking to colostrum. The B lymphocyte compartment was significantly reduced in obese colostrum and these cells did not exhibit enhanced CD16 shedding in this tissue, unlike B lymphocytes from lean mothers colostrum. This is the first exhaustive characterization of major leukocyte subsets in obese mothers colostrum, and the first report of leukocyte subpopulations from Latin-American womens colostrum. This pioneering study is a steppingstone to further investigate active immunity in human breastmilk.


title: T cell-intrinsic vitamin A metabolism and its signaling are targets for memory T cell-based cancer immunotherapy

abstract: Memory T cells play an essential role in infectious and tumor immunity. Vitamin A metabolites such as retinoic acid are immune modulators, but the role of vitamin A metabolism in memory T- cell differentiation is unclear. In this study, we identified retinol dehydrogenase 10 (Rdh10), which metabolizes vitamin A to retinal (RAL), as a key molecule for regulating T cell differentiation. T cell-specific Rdh10 deficiency enhanced memory T-cell formation through blocking RAL production in infection model. Epigenetic profiling revealed that retinoic acid receptor (RAR) signaling activated by vitamin A metabolites induced comprehensive epigenetic repression of memory T cell-associated genes, including TCF7, thereby promoting effector T-cell differentiation. Importantly, memory T cells generated by Rdh10 deficiency and blocking RAR signaling elicited potent anti-tumor responses in adoptive T-cell transfer setting. Thus, T cell differentiation is regulated by vitamin A metabolism and its signaling, which should be novel targets for memory T cell-based cancer immunotherapy.


title: Loss of self-tolerance leads to altered gene expression and IMD pathway activation in Drosophila melanogaster

abstract: Immune self-tolerance is the ability of a hosts immune system to recognize and avoid triggering immune responses against self-tissue. This allows the host to avoid self-directed immune damage while still responding appropriately to pathogen infection. A breakdown of self-tolerance can lead to an autoimmune state in which immune cells target healthy self-tissue, leading to inflammation and tissue damage. In order to better understand the basic biology of autoimmunity and the role of the innate immune system in maintaining self-tolerance, we have recently characterized the Drosophila melanogaster tuSz autoimmune mutant. This mutant strain can serve as a model of innate immune mediated self-tolerance, and here we identify transcripts that are deregulated in flies experiencing a loss of self-tolerance. We found that these changes include the ectopic activation of pro-inflammatory signaling through the Relish/NF{kappa}B transcription factor, alterations in transcripts encoding proteins predicted to mediate organismal metabolism, and a downregulation of transcripts linked to developmental processes. This study can provide insight into the transcriptional and physiological changes underlying self-tolerance and autoimmunity.


title: Polyethylene glycol (PEG) methods are superior to acidification for secondary concentration of Adenovirus and MS2 in water

abstract:
Enteric viruses are a leading cause of waterborne illness worldwide and surveillance studies lack standardization in method selection. The most common and cost-effective approach to concentrating viruses from water samples involves virus adsorption and elution (VIRADEL) procedures, followed by secondary concentration. There is a lack of consistency in how secondary concentration methods are practiced and some methods may have better recovery for particular groups of viruses. Secondary concentration methods typically involve precipitation and the most common methods employ organic flocculation (OF) by acidification at a pH of 3.5, or precipitation by polyethylene glycol (PEG) in combination with the addition of NaCl. In this study, the recovery of coliphage MS2 using the plaque assay and human adenovirus strain 41 (HAdV41) using cell-culture and qPCR assays were evaluated by OF and PEG secondary concentration of spiked samples of wastewater, surface water, and groundwater. The recovery of MS2 and HAdV41 by PEG precipitation was significantly higher than that by OF (p<0.0001) when viruses were detected by culture based methods and marginally better when HAdV41 was enumerated by qPCR (p<0.019). The recovery of HAdV41 by qPCR ranged from 75.3% to 94.4% (n=36). The mean recovery of MS2 by OF was 4.4% (0.9%-7.7%; n=14) and ranged from 57.1% to 87.9% (n=28) for the PEG methods. The poor recovery of MS2 by OF was attributed to inactivation or poor stability at acidic conditions as MS2 were not recovered in the supernatant following OF and centrifugation. The inconsistency and lack of justification for method selection in many studies calls for a systematic study to inform guidance and standardization with respect to the application of concentration methods for various water types and viral pathogens.

IMPORTANCEMS2 should not be used as a process control for methods involving acidification and culture-based detection. The dense floc produced by the PEG method may have contributed to higher recoveries as the pellet was more compact and stable than the loose pellet formed by OF. Standard methods for the detection of enteric viruses and surrogates that involve acidification could be modified with PEG precipitation to uphold virus recovery and minimize inactivation.

category: microbiology


title: The global spread and invasion capacities of alien ants

abstract:
The ecological and economic impacts of biological invasions are usually highly conspicuous, but these are the outcome of a global, multistage process that is obscured from view. For most taxa, we lack a large-scale picture of the movements of alien species, the biases and barriers that promote or inhibit their spread at each stage, and blind spots in our ability to detect species during their spread. For instance, countries rely heavily on customs interceptions to prevent new species introductions, but their efficacy for detecting invaders remains unclear. To address these gaps, we synthesize and provide data at unprecedented geographic resolution on the global diversity and distribution of alien ants, a pervasive group strongly impacting humans and ecosystems. From >145,000 records spanning 602 regions, we identify 522 ant species exhibiting human-mediated spread, doubling recent estimates of their diversity. We show that movements of alien ant species across regions globally are non-random and, moreover, that these flows differ by the extents to which species invade--ranging from arrival to indoor establishment, naturalization, and harmful status. Importantly, we find that almost two-thirds of the 309 species that naturalize globally--most of which are ground- and litter-dwelling--are absent from customs interceptions, which record disproportionately high numbers of arboreal species. Our results illustrate the vast, yet uneven extent of ant invasions globally, and suggest that most alien species bypass biosecurity controls while spreading successfully worldwide. This raises doubts on the efficacy of current customs interceptions procedures and highlights a need for radically new approaches.

Significance statementBiological invasions impact humans and ecosystems tremendously. Invasions are difficult to address because little is known about large-scale patterns of spread, species capacities to invade ecosystems, and the efficacy of existing biosecurity measures. This paper tackles these issues for alien ants, one of the most damaging groups of invasive animals. An unprecedented dataset reveals that 522 ant species exhibit human-mediated spread, 270 more than previously known. Species are grouped under four levels of invasion capacity corresponding to different invasion barriers. Geographic sources and sinks for the groups differ globally. Two-thirds of species with naturalization capacity have never been recorded at customs interceptions, which fail to detect many litter-dwelling species. Novel detection and control strategies for inconspicuous alien organisms are needed.

category: ecology


title: A Time-calibrated Firefly (Coleoptera: Lampyridae) Phylogeny: Using Genomic Data for Divergence Time Estimation

abstract: Fireflies (Coleoptera: Lampyridae) consist of over 2,000 described extant species. A well-resolved phylogeny of fireflies is important for the study of their population genetics, bioluminescence, evolution, and conservation. We used a recently published anchored hybrid enrichment dataset (AHE; 436 loci for 88 Lampyridae species and 10 outgroup species) and state-of-the-art statistical methods (the fossilized birth-death-range process implemented in a Bayesian framework) to estimate a time-calibrated phylogeny of Lampyridae. Unfortunately, estimating calibrated phylogenies using AHE and the latest and most robust time-calibration strategies is not possible because of computational constraints. As a solution, we subset the full dataset by applying three different strategies: (i) using the most complete loci, (ii) using the most homogeneous loci, and (iii) using the loci with the highest accuracy to infer the well established Photinus clade. The estimated topology using the three data subsets agreed on almost all major clades and only showed minor discordance within less supported nodes. The estimated divergence times overlapped for all nodes that are shared between the topologies. Thus, divergence time estimation is robust as long as the topology inference is robust and any well selected data subset suffices. Additionally, we observed an un-expected amount of gene tree discordance between the 436 AHE loci. Our assessment of model adequacy showed that standard phylogenetic substitution models are not adequate for any of the 436 AHE loci which is likely to bias phylogenetic inferences. We performed a simulation study to explore the impact of (a) incomplete lineage sorting, (b) uniformly distributed and systematic missing data, and (c) systematic bias in the position of highly variable and conserved sites. For our simulated data, we observed less gene tree variation which shows that the empirically observed amount of gene tree discordance for the AHE dataset is unexpected and needs further investigation.

category: evolutionary biology


title: A genome assembly of the Atlantic chub mackerel (Scomber colias): a valuable teleost fishing resource

abstract: The Atlantic chub mackerel, Scomber colias Gmelin, 1789, is a medium-size pelagic fish with substantial importance in the fisheries of the Atlantic Ocean and the Mediterranean Sea. Over the past decade, this species has gained special relevance being one of the main targets of pelagic fisheries in the NE Atlantic. Here, we sequenced and annotated the first high-quality draft genome assembly of S. colias, produced with Pacbio HiFi long reads and Illumina Paired-End short reads. The estimated genome size is 814 Mb distributed into 2,028 scaffolds and 2,093 contigs with an N50 length of 4,19 and 3,34 Mb, respectively. We annotated 27,675 protein-coding genes and the BUSCO analyses indicated high completeness, with 97.3 % of the single-copy orthologs in the Actinopterygii library profile. The present genome assembly represents a valuable resource to address the biology and management of this relevant fishery. Finally, this is the fourth high-quality genome assembly within the Order Scombriformes and the first in the genus Scomber.

category: genomics


title: Weighted Shoes in the Wild: Initial Insights into the Relationship Between the Effort of Walking and the Amount of Walking Performed

abstract:
PurposeWalking comprises a large portion of active energy expenditure in humans. Interventions that increase or decrease the energy expended on each step, such as ankle weights or energy-saving orthoses, may therefore strongly impact fitness. The overall effect of increasing or decreasing per-step energy use is unclear, however, because people may choose to walk less or more, respectively, in response.

MethodsIn this study, healthy college students with normal body mass index wore weighted and unweighted shoes embedded with an inertial measurement unit for one week each. Community-based walking data were analyzed for number of steps, distance traveled and walking speed. Oxygen consumption using each set of shoes at a range of speeds were measured in a laboratory setting and used to estimate metabolic energy expended during community-based walking. A survey measured subjective response to each pair of shoes.

ResultsThe weighted shoes increased per-step energy cost by about 26%. Subjects strongly disliked the weighted shoes (P = 0.001) and found them tiring (P = 0.003). Despite this dislike, subjects did not significantly reduce distance walked (P = 0.6), number of steps (P = 0.7), or average speed (P = 0.9) compared to normal shoes. This led to a small but not statistically significant increase in energy expended during walking over a five-day period (12.3 {+/-} 9.6% increase, P = 0.2). On the final collection day this trend appeared to reverse, with fewer steps taken and lower metabolic energy expended with the weighted shoes. Twenty-four subjects were recruited but only ten completed the protocol, with dislike of the weighted shoe condition being the primary reason for dropout.

ConclusionsIncreasing the energy cost of each step led to greater energy expended through walking. However, there are indications that behavioral changes would be greater with a longer intervention or increased retention. For example, the large dropout rate suggests that some subjects avoided walking with the weighted shoes entirely, simply by leaving the study. Follow-on studies among patient groups may reveal a fitness benefit to either increasing or reducing the energy cost of walking.

category: animal behavior and cognition


title: Validity of freehand three-dimensional ultrasound system in measurement of three-dimensional surface shape of shoulder muscles

abstract: Freehand three-dimensional ultrasound (3DUS) system is a promising technique for accurately assessing muscle morphology. However, its accuracy has been validated mainly in terms of volume by examining lower limb muscles. This study aimed to validate 3DUS in the measurements of 3D surface shape and volume by comparing them with MRI measurements while ensuring the reproducibility of participant posture by focusing on the shoulder muscles. The supraspinatus, infraspinatus, and posterior deltoid muscles of 10 healthy males were scanned using 3DUS and MRI while secured by an immobilization support customized for each participant. A 3D surface model of each muscle was created from the 3DUS and MRI methods, and the agreement between them was assessed. For the muscle volume, the mean difference between the two models was within -0.51 cm3. For the 3D surface shape, the distances between the closest points of the two models and Dice similarity coefficient were calculated. The results showed that the median surface distance was less than 1.21 mm and Dice similarity coefficient was larger than 0.85. These results suggest that, given the above error is permitted, 3DUS can be used as an alternative to MRI in measuring volume and surface shape, even for the shoulder muscles.

category: bioengineering


title: Substrate stiffness controls the cell cycle of human mesenchymal stem cells via cellular traction

abstract:
The microenvironment of human mesenchymal stem cells (hMSCs) regulates their self-renewal and differentiation properties. Previously it was shown that hMSCs remained quiescent on soft (0.25 kPa) polyacrylamide (PA) gels but re-entered into cell cycle on a stiff (7.5 kPa) gel. However, how cells behave on intermediate stiffness and what intracellular factors transmit mechanical changes to cell interior thereby regulating cell cycle remained unknown. In this work we demonstrated that PA gels between 1 and 5 kPa act as a mechanical switch in regulating cell cycle of hMSCs. By experiments on cell-cycle exit and re-entry, we found that hMSCs demonstrated a sharp transition from quiescence to proliferation between 1 and 5 kPa. Further studies with ROCK inhibitor Y-27632 revealed that contractile proteins, but not cell spread area, accounts for the sensitivity of hMSCs towards substrate stiffness and hence correlates with their changes in cell cycle. These observations therefore suggest that substrate stiffness regulates hMSC proliferation through contractile forces as generated by cellular contractile proteins in a unique pattern which is distinct from other cell types as studied.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=182 SRC="FIGDIR/small/469207v1_ufig1.gif" ALT="Figure 1">
View larger version (83K):
org.highwire.dtl.DTLVardef@a4fff3org.highwire.dtl.DTLVardef@9fb415org.highwire.dtl.DTLVardef@e7d633org.highwire.dtl.DTLVardef@9fba04_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Robust Identification of Temporal Biomarkers in Longitudinal Omics Studies

abstract: Longitudinal studies increasingly collect rich  omics data sampled frequently over time and across large cohorts to capture dynamic health fluctuations and disease transitions. However, the generation of longitudinal omics data has preceded the development of analysis tools that can efficiently extract insights from such data. In particular, there is a need for statistical frameworks that can identify not only which omics features are differentially regulated between groups but also over what time intervals. Additionally, longitudinal omics data may have inconsistencies, including nonuniform sampling intervals, missing data points, subject dropout, and differing numbers of samples per subject. In this work, we developed a statistical method that provides robust identification of time intervals of temporal omics biomarkers. The proposed method is based on a semi-parametric approach, in which we use smoothing splines to model longitudinal data and infer significant time intervals of omics features based on an empirical distribution constructed through a permutation procedure. We benchmarked the proposed method on five simulated datasets with diverse temporal patterns, and the method showed specificity greater than 0.99 and sensitivity greater than 0.72. Applying the proposed method to the Integrative Personal Omics Profiling (iPOP) cohort revealed temporal patterns of amino acids, lipids, and hormone metabolites that are differentially regulated in male versus female subjects following a respiratory infection. In addition, we applied the longitudinal multi-omics dataset of pregnant women with and without preeclampsia, and the method identified potential lipid markers that are temporally significantly different between the two groups. We provide an open-source R package, OmicsLonDA (Omics Longitudinal Differential Analysis): https://bioconductor.org/packages/OmicsLonDA to enable widespread use.

category: bioinformatics


title: DETIRE: A Hybrid Deep Learning Model for identifying Viral Sequences from Metagenomes

abstract: A metagenome contains all DNA sequences from an environmental sample, including viruses, bacteria, fungi, actinomycetes and so on. Since viruses are of huge abundance and have caused vast mortality and morbidity to human society in history as a kind of major pathogens, detecting viruses from metagenomes plays a crucial role in analysing the viral component of samples and is the very first step for clinical diagnosis. However, detecting viral fragments directly from the metagenomes is still a tough issue because of the existence of huge number of short sequences. In this paper, a hybrid Deep lEarning model for idenTifying vIral sequences fRom mEtagenomes (DETIRE), is proposed to solve the problem. Firstly, the graph-based nucleotide sequence embedding strategy is utilized to enrich the expression of DNA sequences by training an embedding matrix. Then the spatial and sequential features are extracted by trained CNN and BiLSTM networks respectively to improve the feature expression of short sequences. Finally, the two set of features are weighted combined for the final decision. Trained by 220,000 sequences of 500bp subsampled from the Virus and Host RefSeq genomes, DETIRE identifies more short viral sequences (<1,000bp) than three latest methods, DeepVirFinder, PPR-Meta and CHEER. DETIRE is freely available at https://github.com/crazyinter/DETIRE.


title: TCR-BERT: learning the grammar of T-cell receptors for flexible antigen-binding analyses

abstract: The T-cell receptor (TCR) allows T-cells to recognize and respond to antigens presented by infected and diseased cells. However, due to TCRs staggering diversity and the complex binding dynamics underlying TCR antigen recognition, it is challenging to predict which antigens a given TCR may bind to. Here, we present TCR-BERT, a deep learning model that applies self-supervised transfer learning to this problem. TCR-BERT leverages unlabeled TCR sequences to learn a general, versatile representation of TCR sequences, enabling numerous downstream applications. We demonstrate that TCR-BERT can be used to build state-of-the-art TCR-antigen binding predictors with improved generalizability compared to prior methods. TCR-BERT simultaneously facilitates clustering sequences likely to share antigen specificities. It also facilitates computational approaches to challenging, unsolved problems such as designing novel TCR sequences with engineered binding affinities. Importantly, TCR-BERT enables all these advances by focusing on residues with known biological significance. TCR-BERT can be a useful tool for T-cell scientists, enabling greater understanding and more diverse applications, and provides a conceptual framework for leveraging unlabeled data to improve machine learning on biological sequences.


title: nucGEMs probe the biophysical properties of the nucleoplasm

abstract: The cell interior is highly crowded and far from thermodynamic equilibrium. This environment can dramatically impact molecular motion and assembly, and therefore influence subcellular organization and biochemical reaction rates. These effects depend strongly on length-scale, with the least information available at the important mesoscale (10-100 nanometers), which corresponds to the size of crucial regulatory molecules such as RNA polymerase II. It has been challenging to study the mesoscale physical properties of the nucleoplasm because previous methods were labor-intensive and perturbative. Here, we report nuclear Genetically Encoded Multimeric nanoparticles (nucGEMs). Introduction of a single gene leads to continuous production and assembly of protein-based bright fluorescent nanoparticles of 40 nm diameter. We implemented nucGEMs in budding and fission yeast and in mammalian cell lines. We found differences in particle motility between the nucleus and the cytosol at the mesoscale, that mitotic chromosome condensation ejects nucGEMs from the nucleus, and that nucGEMs are excluded from heterochromatin and the nucleolus. nucGEMs enable hundreds of nuclear rheology experiments per hour, and allow evolutionary comparison of the physical properties of the cytosol and nucleoplasm.


title: Heterogeneous run-and-tumble motion accounts for transient non-Gaussian super-diffusion in haematopoietic multi-potent progenitor cells

abstract: Multi-potent progenitor (MPP) cells act as a key intermediary step between haematopoietic stem cells and the entirety of the mature blood cell system. Their eventual fate determination is thought to be achieved through migration in and out of spatially distinct niches. Here we first analyze statistically MPP cell trajectory data obtained from a series of long time-course 3D in-vivo imaging experiments on irradiated mouse calvaria, and report that MPPs display transient super-diffusion with apparent non-Gaussian displacement distributions. Second, we explain these experimental findings using a run-and-tumble model of cell motion which incorporates the observed dynamical heterogeneity of the MPPs. Third, we use our model to extrapolate the dynamics to time-periods currently inaccessible experimentally, which enables us to quantitatively estimate the time and length scales at which super-diffusion transitions to Fickian diffusion. Our work sheds light on the potential importance of motility in early haematopoietic progenitor function.


title: Restoring single-molecule localizations with wavefront sensing adaptive optics for deep-tissue super-resolution imaging

abstract: Specimen-induced aberration has been one of the major factors limiting the imaging depth in single-molecule localization microscopy (SMLM). In this study, we measured the wavefront of intrinsic reflectance signal at the fluorescence emission wavelength to construct a time-gated reflection matrix and find complex tissue aberration without resorting to fluorescence detection. Physically correcting the identified aberration via wavefront shaping with a liquid-crystal spatial light modulator (SLM)enables super-resolution imaging even when the aberration is too severe for initiating localization processes. We demonstrate the correction of complex tissue aberration, the root-mean-square (RMS)wavefront distortion of which is more than twice the 1 rad limit presented in previous studies; this leads to the recovery of single molecules by 77 times increased localization number. We visualised dendritic spines in mouse brain tissues and early myelination processes in a whole zebrafish at up to 102 m depth with 28-39 nm localization precision. The proposed approach can expand the application range of SMLM to thick samples that cause the loss of localization points owing to severe aberration.


title: Study on the effect of a triple cancer treatment of propolis, thermal cycling-hyperthermia, and low-intensity ultrasound on PANC-1 cells

abstract:
BackgroundPancreatic cancer is a deadly cancer around the world. To reduce side effects and enhance treatment efficacy, study on combination therapy for pancreatic cancer has gained much attention in recent years.

MethodsIn this paper, we propose a novel triple treatment combining propolis and two physical stimuli-thermal cycling-hyperthermia (TC-HT) and low-intensity ultrasound (US) on a human pancreatic cancer cell line PANC-1. MTT assay was used to determine the viability of PANC-1 cells. Flow cytometry was used to detect apoptosis, mitochondrial membrane potential (MMP) loss, and intracellular reactive oxygen species (ROS) levels. Western blot analysis was further performed to measure protein expression and phosphorylation.

ResultsThe experiments found that, after the triple treatment, the cell viability of the PANC-1 cells decreased to a level 80% less than the control, without affecting the normal pancreatic cells. Another result was excessive accumulation of ROS after the triple treatment, leading to the amplification of apoptotic pathway through the mitogen-activated protein kinase (MAPK) family and mitochondrial dysfunction. Moreover, the combination of TC-HT and US also promotes the anticancer effect of the heat-sensitive chemotherapy drug cisplatin on PANC-1 cells.

ConclusionThis study, to the best of our knowledge, is the first attempt to combine TC-HT, US and a nature compound in cancer treatment. We demonstrate that physical stimuli could augment the therapeutical effect of anticancer agents. It is expected that optimized parameters for different agents and different types of cancer will expand the methodology on oncological therapy in a safe manner.


title: A chemerin peptide analog stimulates tumor growth in two xenograft mouse models of human colorectal carcinoma

abstract:
BackgroundChemerin plasma concentration has been reported to be positively correlated with the risk of colorectal cancer. However, the potential regulation of CRC tumorigenesis and progression has not yet been investigated in an experimental setting. This study addresses this hypothesis by investigating proliferation, colony formation and migration of CRC cell lines in vitro as well as in animal models.

MethodsIn vitro, microscopic assays to study proliferation as well as a scratch-wound assay for migration monitoring were applied using the human CRC cell lines HCT116, HT29 and SW620 under the influence of the chemerin analog CG34. The animal study investigated HCT116-luc and HT29-luc subcutaneous tumor size and bioluminescence during treatment with CG34 versus control, followed by ex-vivo analysis of vessel density and mitotic activity.

ResultsWhile proliferation of the three CRC cell lines in monolayers was not clearly stimulated by CG34, the chemerin analog promoted colony formation in three-dimensional aggregates. An effect on cell migration was not observed. In the treatment study, CG34 significantly stimulated both growth and bioluminescence signal of HCT116-luc and HT29-luc xenografts.

ConclusionsThe results of this study represent the first indication of a tumor growth-stimulating effect of chemerin signaling in CRC.


title: Preclinical and Toxicology Studies of BRD5529, a Selective Inhibitor of CARD9

abstract:
BackgroundExuberant inflammation during Pneumocystis pneumonia leads to lung injury. CARD9 is a central mediator of inflammatory signaling mediated by C-type lectin receptors. CARD9 inhibitor BRD5529 has been shown to be an effective in vitro inhibitor of Pneumocystis {beta}-glucan-induced proinflammatory signaling and downstream TNF-alpha production, suggesting its viability as a candidate for preliminary drug testing as an anti- inflammatory agent in the rodent Pneumocystis pneumonia model (PCP).

MethodsTo assess for potential toxicity, mice were injected intraperitoneally (IP) daily either with vehicle or BRD5529 at 0.1 mg/kg or 1.0 mg/kg for two weeks. Mouse weights were taken daily. At day 14, mice were euthanized, weighed, and analyzed by flexiVent for lung stiffness. Lungs, liver, and kidney were then harvested for H&E staining and pathology scoring. Lung samples were further analyzed for proinflammatory cytokines via ELISA and extracellular matrix generation via quantitative PCR (q-PCR). Blood collection postmortem was performed for blood chemistry analysis.

ResultsBRD5529 at both doses of IP administration resulted in no significant changes in daily or final weight gain. Analysis of lung stiffness by flexiVent showed no significant differences between the control or treated groups. Furthermore, ELISA results for proinflammatory IL-1 Beta, IL-6, and TNF-alpha showed no major differences in the respective groups. qPCR analysis of extracellular matrix transcripts collagen type I, alpha 1 (Col1a1) and fibronectin (Fn) were statically similar as well in the treated and control groups. Examination and pathology scoring of H&E slides from lung, liver, and kidney from the each of the mice in all groups and subsequent pathology scoring showed no significant change. Blood chemistry analysis revealed similar, non-significant patterns.

ConclusionsBRD5529 in our initial general safety and toxicology assessments displayed no inherent safety concerns in the analyzed parameters. These data support broader in vivo testing of the inhibitor as a timed adjunct therapy to the deleterious proinflammatory host immune response often associated with anti-Pneumocystis therapy.


title: Partial prevention of glucocorticoid-induced osteocyte deterioration with osteocrin gene therapy

abstract: Glucocorticoid (GC)-induced osteoporosis and subsequent bone fragility are preceded by death and dysfunction at the cellular level. In particular, short-term glucocorticoid excess suppresses osteocyte remodeling of the surrounding bone mineral, causes apoptosis of osteoblasts and osteocytes, and disrupts homeostatic bone remodeling. Preventing apoptosis and preserving osteocyte morphology and function could be effective means of preventing bone loss during glucocorticoid excess. We hypothesized that osteocrin, which preserves osteocyte viability and morphology in other models where osteocyte defects exist, could prevent osteocyte death and dysfunction in a GC excess model. We used a liver-targeted adeno-associated virus (AAV8) to induce osteocrin overexpression in mice one week prior to implantation with prednisolone or placebo pellets. After 28 days, tissues were collected for micro-CT and histological analysis. GC excess caused the expected reduction in cortical bone thickness and osteocyte canalicular length in control AAV8-treated mice, and these effects were blunted in mice overexpressing osteocrin. However, GC-induced changes in cortical porosity, trabecular bone mass, and gene expression were not prevented by osteocrin. While the mechanism of osteocrins effects on osteocyte morphology warrants further investigation, this study does not support a role for this model of osteocrin supplementation to combat the full skeletal effects of GC excess.

category: physiology


title: Metaeffector interactions modulate the type III effector-triggered immunity load of Pseudomonas syringae

abstract: The bacterial plant pathogen Pseudomonas syringae requires type III secreted effectors (T3SEs) for pathogenesis. However, a major facet of plant immunity entails the recognition of a subset of P. syringaes T3SEs by intracellular host receptors in a process called Effector- Triggered Immunity (ETI). Prior work has shown that ETI-eliciting T3SEs are pervasive in the P. syringae species complex raising the question of how P. syringae mitigates its ETI load to become a successful pathogen. While pathogens can evade ETI by T3SE mutation, recombination, or loss, there is increasing evidence that effector-effector (a.k.a., metaeffector) interactions can suppress ETI. To study the ETI-suppression potential of P. syringae T3SE repertoires, we compared the ETI-elicitation profiles of two genetically divergent strains: P. syringae pv. tomato DC3000 (PtoDC3000) and P. syringae pv. maculicola ES4326 (PmaES4326), which are both virulent on Arabidopsis thaliana but harbour largely distinct effector repertoires. Of the 529 T3SE alleles screened on A. thaliana Col-0 from the P. syringae T3SE compendium (PsyTEC) [1], 69 alleles from 21 T3SE families elicited ETI in at least one of the two strain backgrounds, while 50 elicited ETI in both backgrounds, resulting in 19 differential ETI responses including two novel ETI-eliciting families: AvrPto1 and HopT1. Although most of these differences were quantitative, three ETI responses were completely absent in one of the pathogenic backgrounds. We performed ETI suppression screens to test if metaeffector interactions contributed to these ETI differences, and found that HopQ1a suppressed AvrPto1m-mediated ETI, while HopG1c and HopF1g suppressed HopT1b-mediated ETI. Overall, these results show that P. syringae strains leverage metaeffector interactions and ETI suppression to overcome the ETI load associated with their native T3SE repertoires.

category: plant biology


title: Smart soils to observe hidden rhizosphere processes

abstract: Agriculture must reduce green-house gas emission and pollution, produce safer and healthier food, closer to home, reducing waste whilst delivering more diverse diets to a growing world population. Soils could enable this transformation, but unfortunately, they have a hugely complex and opaque structure and studies of its myriad of mechanisms are difficult. Here, the fabrication of smart soils for the screening of below-ground bio-processes is demonstrated. Particles were generated from fluoropolymer waste with functionalisation converting them into sensors able to report on key chemical dynamics. Tailored functionalization was obtained by radical terpolymerisation to improve growth conditions and sensing capabilities. The study demonstrates the potential for the development of accelerated genetic or agrochemical screens and could pave the way for improved models for rhizosphere dynamics.


title: A human kinase yeast array for the identification of kinases modulating phosphorylation-dependent protein-protein interactions

abstract: Protein kinases play an important role in cellular signaling pathways and their dysregulation leads to multiple diseases, making kinases prime drug targets. While more than 500 human protein kinases are known to collectively mediate phosphorylation of over 290,000 S/T/Y sites, the activities have been characterized only for a minor, intensively studied subset. To systematically address this discrepancy, we developed a human kinase array in Saccharomyces cerevisiae as a simple readout tool to systematically assess kinase activities. For this array, we expressed 266 human kinases in four different Saccharomyces cerevisiae strains and profiled ectopic growth as a proxy for kinase activity across 33 conditions. More than half of the kinases showed an activity-dependent phenotype across many conditions and in more than one strain. We then employed the kinase array to identify the kinase(s) that can modulate protein-protein-interactions (PPIs). Two characterized, phosphorylation-dependent PPIs with unknown kinase-substrate relationships were analyzed in a phospho-yeast two-hybrid assay. CK21 and SGK2 kinases can abrogate the interaction between the spliceosomal proteins AAR2 and PRPF8 and NEK6 kinase was found to mediate the estrogen receptor (ER) interaction with 14-3-3 proteins. The human kinase yeast array can thus be used for a variety of kinase activity-dependent readouts.

category: systems biology


title: Passive Exposure Sparsifies Neural Activity in the Primary Visual Cortex

abstract: A notable feature of neural activity is sparseness - namely, that only a small fraction of neurons in a local circuit have high activity at any moment. Not only is sparse neural activity observed experimentally in most areas of the brain, but sparseness has been proposed as an optimization or design principle for neural circuits. Sparseness can increase the energy efficiency of the neural code as well as allow for beneficial computations to be carried out. But how does the brain achieve sparse-ness? Here, we found that when neurons in the primary visual cortex were passively exposed to a set of images over several days, neural responses became more sparse. Sparsification was driven by a decrease in the response of neurons with low or moderate activity, while highly active neurons retained similar responses. We also observed a net decorrelation of neural activity. These changes sculpt neural activity for greater coding efficiency.


title: A Biophysical Model of Nonquantal Transmission at the Vestibular Hair Cell-Calyx Synapse: KLV currents Modulate Fast Electrical and Slow K+ potentials in the Synaptic Cleft

abstract:
Vestibular hair cells transmit information about head position and motion across synapses to primary afferent neurons. At some of these synapses, the afferent neuron envelopes the hair cell, forming an enlarged synaptic terminal called a calyx. The vestibular hair cell-calyx synapse supports a mysterious form of electrical transmission that does not involve gap junctions termed nonquantal transmission (NQT). The NQT mechanism is thought to involve the flow of ions from the pre-synaptic hair cell to the post-synaptic calyx through low-voltage-activated channels driven by changes in cleft [K+] as K+ exits the hair cell. However, this hypothesis has not been tested with a quantitative model and the possible role of an electrical potential in the cleft has remained speculative. Here we present a computational model that captures salient experimental observations of NQT and identifies overlooked features that corroborate the existence of an electrical potential ({phi}) in the synaptic cleft. We show that changes in cleft{phi} reduce transmission latency and illustrate the relative contributions of both cleft [K+] and{phi} to the gain and phase of NQT. We further demonstrate that the magnitude and speed of NQT depend on calyx morphology and that increasing calyx height reduces action potential latency in the calyx afferent. These predictions are consistent with the idea that the calyx evolved to enhance NQT and speed up vestibular signals that drive neural circuits controlling gaze, balance, and orientation.

Significance StatementThe ability of the vestibular system to drive the fastest reflexes in the nervous system depends on rapid transmission of mechanosensory signals at vestibular hair cell synapses. In mammals and other amniotes, afferent neurons form unusually large calyx terminals on certain hair cells, and communication at these synapses includes nonquantal transmission (NQT), which avoids the synaptic delay of quantal transmission. We present a quantitative model that shows how NQT depends on the extent of the calyx covering the hair cell and attributes the short latency of NQT to changes in synaptic cleft electrical potential caused by current flowing through open potassium channels in the hair cell. This previously undescribed mechanism may act at other synapses.


title: Nonquantal Transmission at the Vestibular Hair Cell-Calyx Synapse: KLV currents Modulate Fast Electrical and Slow K+ potentials in the Synaptic Cleft


title: WWOX-Mediated Degradation of AMOTp130 Negatively Affects Egress of Filovirus VP40 VLPs

abstract:
Ebola (EBOV) and Marburg (MARV) viruses continue to emerge and cause severe hemorrhagic disease in humans. A comprehensive understanding of the filovirus-host interplay will be crucial for identifying and developing antiviral strategies. The filoviral VP40 matrix protein drives virion assembly and egress, in part by recruiting specific WW-domain-containing host interactors via its conserved PPxY Late (L) domain motif to positively regulate virus egress and spread. In contrast to these positive regulators of virus budding, a growing list of WW-domain-containing interactors that negatively regulate virus egress and spread have been identified, including BAG3, YAP/TAZ and WWOX. In addition to host WW-domain regulators of virus budding, host PPxY-containing proteins also contribute to regulating this late stage of filovirus replication. For example, angiomotin (AMOT) is a multi-PPxY-containing host protein that functionally interacts with many of the same WW-domain-containing proteins that regulate virus egress and spread. In this report, we demonstrate that host WWOX, which negatively regulates egress of VP40 VLPs and recombinant VSV-M40 virus, interacts with and suppresses the expression of AMOT. We found that WWOX disrupts AMOTs scaffold-like tubular distribution and reduces AMOT localization at the plasma membrane via lysosomal degradation. In sum, our findings reveal an indirect and novel mechanism by which modular PPxY/WW-domain interactions between AMOT and WWOX regulate PPxY-mediated egress of filovirus VP40 VLPs. A better understanding of this modular network and competitive nature of protein-protein interactions will help to identify new antiviral targets and therapeutic strategies.

IMPORTANCEFiloviruses (Ebola [EBOV] and Marburg [MARV]) are zoonotic, emerging pathogens that cause outbreaks of severe hemorrhagic fever in humans. A fundamental understanding of the virus-host interface is critical for understanding the biology of these viruses and for developing future strategies for therapeutic intervention. Here, we reveal a novel mechanism by which host proteins WWOX and AMOTp130 interact with each other and with the EBOV matrix protein VP40 to regulate VP40-mediated egress of virus like particles (VLPs). Our results highlight the biological impact of competitive interplay of modular virus-host interactions on both the virus lifecycle and the host cell.


title: The developmentally-timed decay of an essential microRNA family is seed sequence-dependent

abstract: MicroRNA (miRNA) abundance is tightly controlled by regulation of biogenesis and decay. Here we show that the mir-35 miRNA family undergoes regulated decay at the transition from embryonic to larval development in C. elegans. The seed sequence of the miRNA is necessary and sufficient for this regulation. Sequences outside the seed (3 end) regulate mir-35 abundance in the embryo but are not necessary for sharp decay at the transition to larval development. Enzymatic modifications of the miRNA 3 end are neither prevalent nor correlated with changes in decay, suggesting that miRNA 3 end display is not a core feature of this mechanism and further supporting a seed-driven decay model. Our findings demonstrate that seed sequence-specific decay can selectively and coherently regulate all redundant members of a miRNA seed family, a class of mechanism that has great biological and therapeutic potential for dynamic regulation of a miRNA familys target repertoire.

category: molecular biology


title: Connexinplexity: The spatial and temporal expression of connexin genes during vertebrate organogenesis

abstract: Animal development requires coordinated communication between cells. The Connexin family of proteins is a major contributor to intercellular communication in vertebrates by forming gap junction channels that facilitate the movement of ions, small molecules, and metabolites between cells. Additionally, individual hemichannels can provide a conduit to the extracellular space for paracrine and autocrine signaling. Connexin-mediated communication is well appreciated in epithelial, neural, and vascular development and homeostasis, and most tissues likely use this form of communication. In fact, Connexin disruptions are of major clinical significance contributing to disorders developing from all major germ layers. Despite the fact that Connexins serve as an essential mode of cellular communication, the temporal and cell-type specific expression patterns of connexin genes remain unknown in vertebrates. A major challenge is the large and complex connexin gene family. To overcome this barrier, we probed the expression of all connexins in zebrafish using single-cell RNA-sequencing of entire animals across several stages of organogenesis. Our analysis of expression patterns has revealed that few connexins are broadly expressed, but rather, most are expressed in tissue- or cell-type-specific patterns. Additionally, most tissues possess a unique combinatorial signature of connexin expression with dynamic temporal changes across the organism, tissue, and cell. Our analysis has identified new patterns for well-known connexins and assigned spatial and temporal expression to genes with no-existing information. We provide a field guide relating zebrafish and human connexin genes as a critical step towards understanding how Connexins contribute to cellular communication and development throughout vertebrate organogenesis.

category: developmental biology


title: CLASS-II KNOX genes coordinate spatial and temporal patterns of the tomato ripening

abstract:
Ripening is a complex developmental change of a mature organ, the fruit. In plants like a tomato, it involves softening, pigmentation, and biosynthesis of metabolites beneficial for the human diet. Examination of the transcriptional changes towards ripening suggests that redundant uncharacterized factors may be involved in the coordination of the ripening switch. Previous studies have demonstrated that Arabidopsis CLASS-II KNOX genes play a significant role in controlling the maturation of siliques and their transition to senescence. Here we examined the combined role of all four tomato CLASS-II KNOX genes in the maturation and ripening of fleshy fruits using an artificial microRNA targeting them simultaneously. As expected, the knockdown plants (35S::amiR-TKN-CL-II) exhibited leaves with increased complexity, reminiscent of the leaf phenotype of plants overexpressing CLASS-I KNOX, which antagonize CLASS-II KNOX gene functions. The fruits of 35S::amiR-TKN-CL-II plants were notably smaller than the control. While their internal gel/placenta tissue softened and accumulated the typical pigmentation, the pericarp color break took place ten days later than control, and eventually, it turned yellow instead of red.

Additionally, the pericarp of 35S::amiR-TKN-CL-II fruits remained significantly firmer than control even after three weeks of shelf storage. Strikingly, the 35S::amiR-TKN-CL-II fruits showed early ethylene release and respiration peak, but these were correlated only with liquefaction and pigmentation of the internal tissues. Our findings suggest that CLASS-II KNOX genes are required to coordinate the spatial and temporal patterns of tomato fruit ripening.

One sentence summaryTomato CLASS-II KNOX genes play antagonistic roles in the regulation of ripening at the internal fruit domains and pericarp.


title: Interrogating Genomic Data in the Phylogenetic Placement of Treeshrews Reveals Potential Sources of Conflict

abstract: The position of some taxa on the Tree of Life remains controversial despite the increase in genomic data used to infer phylogenies. While analyzing large datasets alleviates stochastic errors, it does not prevent systematic errors in inference, caused by both biological (e.g., incomplete lineage sorting, hybridization) and methodological (e.g., incorrect modeling, erroneous orthology assessments) factors. In our study, we systematically investigated factors that could result in these controversies, using the treeshrew (Scandentia, Mammalia) as a study case. Recent studies have narrowed the phylogenetic position of treeshrews to three competing hypotheses: sister to primates and flying lemurs (Primatomorpha), sister to rodents and lagomorphs (Glires), or sister to a clade comprising all of these. We sampled 50 mammal species including three treeshrews, a selection of taxa from the potential sister groups, and outgroups. Using a large diverse set of loci, we assessed support for the alternative phylogenetic position of treeshrews. A plurality of loci support treeshrews as sister to rodents and lagomorphs; however, only a few loci exhibit strong support for any hypothesis. Surprisingly, we found that a subset of loci that strongly support the monophyly of Primates, support treeshrews as sister to primates and flying lemurs. The overall small magnitude of differences in phylogenetic signal among the alternative hypotheses suggests that these three groups diversified nearly simultaneously. However, with our large dataset and approach to examining support, we provide evidence for the hypothesis of treeshrews as sister to rodents and lagomorphs, while demonstrating why support for alternate hypotheses has been seen in prior work. We also suggest that locus selection can unwittingly bias results.


title: Indirect effects of parental conflict on conspecific offspring development

abstract: Hybrid Seed Inviability (HSI) is a common barrier in angiosperms. Recent work suggests that the rapid evolution of HSI may, in part, be due to conflict between maternal and paternal optima for resource allocation to developing offspring (i.e. parental conflict). However, parental conflict requires that paternally-derived resource acquiring alleles impose a maternal cost. I test this requirement using three closely related species in the Mimulus guttatus species complex that exhibit significant HSI and differ in their inferred histories of parental conflict. I show that the presence of hybrid seeds significantly affects conspecific seed size for almost all crosses, such that conspecific seeds are smaller after developing with hybrids from fathers with a stronger history of conflict, and larger after developing with hybrids from fathers with a weaker history of conflict. This work demonstrates a cost of paternally-derived alleles, and also has implications for species fitness in secondary contact.


title: Genetic mapping of genotype-by-ploidy effects in Arabidopsis thaliana

abstract: Plants can express different phenotypic responses following polyploidization, but ploidy-dependent phenotypic variation has so far not been assigned to specific genetic factors. To map such effects, segregating populations at different ploidy levels are required. The availability of an efficient haploid-inducer line in Arabidopsis thaliana allows for the rapid development of large populations of segregating haploid offspring. Because Arabidopsis haploids can be self-fertilised to give rise to homozygous doubled haploids, the same genotypes can be phenotyped at both the haploid and diploid ploidy level. Here, we compared the phenotypes of recombinant haploid and diploid offspring derived from a cross between two late flowering accessions to map genotype x ploidy (GxP) interactions. Ploidy-specific quantitative trait loci (QTLs) were detected at both ploidy levels. This implies that mapping power will increase when phenotypic measurements of monoploids are included in QTL analyses. A multi-trait analysis further revealed pleiotropic effects for a number of the ploidy specific QTLs as well as opposite effects at different ploidy levels for general QTLs. Taken together, we provide evidence of genetic variation between different Arabidopsis accessions being causal for dissimilarities in phenotypic responses to altered ploidy levels, revealing a GxP effect. Additionally, by investigating a population derived from late flowering accessions we revealed a major vernalisation specific QTL for variation in flowering time, countering the historical bias of research in early flowering accessions.

category: genetics


title: Evidence for the role of transcription factors in the co-transcriptional regulation of intron retention

abstract:
Alternative splicing is a widespread regulatory phenomenon that enables a single gene to produce multiple transcripts. Among the different types of alternative splicing, intron retention is one of the least explored despite its high prevalence in both plants and animals. The recent discovery that the majority of splicing is co-transcriptional has led to the finding that chromatin state affects alternative splicing. Therefore it is plausible that transcription factors can regulate splicing outcomes. We provide evidence for this hypothesis by studying regions of open chromatin in retained and excised introns. Using deep learning models designed to distinguish between regions of open chromatin in retained introns and non-retained introns, we identified motifs enriched in IR events with significant hits to known human transcription factors. Our model predicts that the majority of transcription factors that affect intron retention come from the zinc finger family. We demonstrate the validity of these predictions using ChIP-seq data for multiple zinc finger transcription factors and find strong over-representation for their peaks in intron retention events.

AvailabilitySource code available at https://github.com/fahadahaf/chromir


title: On the hindlimb biomechanics of the avian take-off leap

abstract: Although extant land birds take to the air by leaping, generating the initial take-off velocity primarily from the hindlimbs, the detailed musculoskeletal mechanics remain largely unknown. We therefore simulated in silico the take-off leap of the zebra finch, Taeniopygia guttata, a model species of passerine, a class of bird which includes over half of all extant bird species. A 3D computational musculoskeletal model of the zebra finch hindlimb, comprising of 43 musculotendon units was developed and driven with previously published take-off ground reaction forces and kinematics. Using inverse dynamics, the external moments at the ankle, knee, and hip joints were calculated and contrasted to the cumulative muscle capability to balance these moments. Mean peak external flexion moments at the hip and ankle were 0.55 bodyweight times leg length (BWL) each whilst peak knee extension moments were about half that value (0.29 BWL). Muscles had the capacity to generate 146%, 230%, and 212 % of the mean peak external moments at the hip, knee, and ankle, respectively. Similarities in hindlimb morphology and external loading across passerine species suggest that the effective take-off strategy employed by the zebra finch may be shared across the passerine clade and therefore half of all birds.


title: Characterization and simulation of metagenomic nanopore sequencing data with Meta-NanoSim

abstract: Nanopore sequencing is crucial to metagenomic studies as its kilobase-long reads can contribute to resolving genomic structural differences among microbes. However, platform-specific challenges, including high base-call error rate, non-uniform read lengths, and the presence of chimeric artifacts, necessitate specifically designed analytical tools. Here, we present Meta-NanoSim, a fast and versatile utility that characterizes and simulates the unique properties of nanopore metagenomic reads. Further, Meta-NanoSim improves upon state-of-the-art methods on microbial abundance estimation through a base-level quantification algorithm. We demonstrate that Meta-NanoSim simulated data can facilitate the development of metagenomic algorithms and guide experimental design through a metagenomic assembly benchmarking task.


title: Disorder in Ca2+ Release Unit Locations Confers Robustness but Cuts Flexibility of Heart Pacemaking

abstract:
Excitation-contraction coupling kinetics are dictated by the rate and rhythm of the excitations generated by sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which accelerates diastolic depolarization and determines the rate and rhythm of the excitations. The LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While the spatial CRU distribution in pacemaker cells exhibits substantial heterogeneity, it remains unknown if it has any functional importance. Using numerical modeling, here we showed that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad lower range of realistic ICaL densities. Allowing each CRU to deviate from its original lattice position fundamentally changes the model behavior: during diastolic depolarization sparks propagate, forming LCRs observed experimentally. As disorder in the CRU positions increases, the CRU distribution exhibits larger empty spaces but simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing AP firing rate and reviving pacemaker function within lower ICaL densities at which cells with lattice CRU distribution were dormant/non-firing. However, cells with fully disordered CRU positions cannot reach low firing rates and their {beta}-adrenergic receptor stimulation effect was substantially decreased. Thus, order/disorder in CRU locations regulates Ca release propagation and could be harnessed by pacemaker cells to regulate their function. Excessive disorder is expected to limit heart rate range that may contribute to heart rate range decline with age and in disease.

SummaryThe present numerical modeling study shows that disorder in locations of Ca release units in cardiac pacemaker cells has substantial functional impact by creating release clusters, similar to Poisson clumping, and opportunity of Ca release to propagate within the clusters.

category: biophysics


title: Abcc6 null mice a model for mineralization disorder PXE show vertebral osteopenia without enhanced intervertebral disc calcification with aging

abstract:
Chronic low back pain is a highly prevalent health condition intricately linked to intervertebral disc degeneration. One of the prominent features of disc degeneration that is commonly observed with aging is dystrophic calcification. ATP-binding cassette sub-family C member 6 (ABCC6), a presumed ATP efflux transporter, is a key regulator of systemic levels of the mineralization inhibitor pyrophosphate (PPi). Mutations in ABCC6 result in pseudoxanthoma elasticum (PXE), a progressive human metabolic disorder characterized by mineralization of the skin and elastic tissues. The implications of ABCC6 loss-of-function on pathological mineralization of structures in the spine, however, are unknown. Using the ABCC6-/- mouse model of PXE, we investigated age-dependent changes in the vertebral bone and intervertebral disc. ABCC6-/- mice exhibited diminished trabecular bone quality parameters at 7-months which remained significantly lower than the wild-type mice at 18 months-of-age. ABCC6-/- vertebrae showed increased TRAP staining along with decreased TNAP staining, suggesting an enhanced bone resorption as well as decreased bone formation. Surprisingly, however, loss of ABCC6 resulted only in a mild, aging disc phenotype without evidence of dystrophic mineralization. Finally, we tested the utility of oral K3Citrate to treat the vertebral phenotype since it is shown to regulate hydroxyapatite mechanical behavior. The treatment resulted in inhibition of osteoclastic response and an early improvement in mechanical properties of the bone underscoring the promise of potassium citrate as a therapeutic agent. Our data suggest that although ectopic mineralization is tightly regulated in the disc, loss of ABCC6 compromises vertebral bone quality and dysregulates osteoblast-osteoclast coupling.

Author SummaryInherited mutations in the ABCC6 transporter gene results in mineralization, often in the form of hydroxyapatite, of connective tissues throughout the body, predominantly affecting the skin, eyes, and blood vessels. Functional loss of ABCC6 causes reduced levels of the potent mineralization inhibitor pyrophosphate (PPi) in blood resulting in these pathologies. Pathological mineralization is also a prominent feature of intervertebral disc degeneration, but the role of ABCC6 and systemic PPi levels and its correlation to disc mineralization and vertebral bone health has remained unexplored. In this study, we show for the first time that loss of ABCC6 in mice results in significant decline in vertebral bone quality and mild age-related disc degeneration without increased incidence of abnormal mineralization. Importantly, treatment of ABCC6 deficient mice with K3Citrate resulted in restoration of early cellular changes which drive bone loss and mechanical function of the vertebrae. In summary, our data reveal that ABCC6 is dispensable for mineralization prevention in the intervertebral disc. Unexpectedly, we found that vertebral bone quality and bone cell activities are linked to ABCC6 function.

category: cell biology


title: A tripartite structure, the complex nuclear receptor element (cNRE), is a cis-regulatory module of viral origin required for atrial chamber preferential gene expression

abstract: Optimal cardiac function requires appropriate contractile proteins in each heart chamber. Atria require slow myosins to act as variable reservoirs, while ventricles demand fast myosin for swift pumping functions. Hence, myosin is under chamber-biased cis-regulatory control to achieve this functional distribution. Failure in proper regulation of myosin genes can lead to severe congenital heart dysfunction. The precise regulatory input leading to cardiac chamber-biased expression remains uncharted. To address this, we computationally and molecularly dissected the quail Slow Myosin Heavy Chain III (SMyHC III) promoter that drives specific gene expression to the atria to uncover the regulatory information leading to chamber expression and understand their evolutionary origins. We show that SMyHC III gene states are autonomously orchestrated by a complex nuclear receptor cis-regulatory element (cNRE), a 32- bp sequence with hexanucleotide binding repeats. Using in vivo transgenic assays in zebrafish and mouse models, we demonstrate that preferential atrial expression is achieved by the combinatorial regulatory input composed of atrial activation motifs and ventricular repression motifs. Through comparative genomics, we provide evidence that the cNRE emerged from an endogenous viral element, most likely through infection of an ancestral host germline. Our study reveals an evolutionary pathway to cardiac chamber-specific expression.


title: Allele segregation analysis of F1 hybrids between independent Brassica  allohexaploid lineages

abstract: In the Brassica genus we find both diploid species (one genome) and allotetraploid species (two different genomes) but no naturally occurring hexaploid species (three different genomes, AABBCC). Although hexaploids can be produced via human intervention, these neo- polyploids have quite unstable genomes and usually suffer from severe genome reshuffling. Whether these genome rearrangements continue in later generations and whether genomic arrangements follow similar, reproducible patterns between different lines is still unknown. We crossed Brassica hexaploids resulting from different species combinations to produce five F1 hybrids, and analyzed the karyotypes of the parents and the F1 hybrids, as well as allele segregation in a resulting test-cross population via molecular karyotyping using SNP array genotyping. Although some genomic regions were found to be more likely to be duplicated, deleted or rearranged, a consensus pattern was not shared between genotypes. Brassica hexaploids had a high tolerance for fixed structural rearrangements, but which rearrangements occur and become fixed over many generations does not seem to show either strong reproducibility or to indicate selection for stability. On average, we observed 10 de novo chromosome rearrangements contributed almost equally from both parents to the F1 hybrids. At the same time, the F1 hybrid meiosis produced on average 8.6 new rearrangements. Hence, the increased heterozygosity in the F1 hybrid did not significantly improve genome stability in our hexaploid hybrids, and might have had the opposite effect. However, hybridization between lineages was readily achieved and may be exploited for future genetics and breeding purposes.


title: Set1-mediated histone H3K4 methylation is required for azole induction of the ergosterol biosynthesis genes and antifungal drug resistance in Candida glabrata.

abstract:
Candida glabrata is an opportunistic pathogen that has developed the ability to adapt and thrive under azole treated conditions. The common mechanisms that can result in Candida drug resistance are due to mutations or overexpression of the drug efflux pump or the target of azole drugs, Cdr1 and Erg11, respectively. However, the role of epigenetic histone modifications in azole-induced gene expression and drug resistance are poorly understood in C. glabrata. In this study, we show for the first time that Set1 mediates histone H3K4 mono-, di-, and trimethylation in C. glabrata. In addition, loss of SET1 and histone H3K4 methylation results in increased susceptibility to azole drugs in both C. glabrata and S. cerevisiae. Intriguingly, this increase in susceptibility to azole drugs in strains lacking Set1-mediated histone H3K4 methylation is not due to altered transcript levels of CDR1, PDR1 or Cdr1s ability to efflux drugs. Genome-wide transcript analysis revealed that Set1 is necessary for azole-induced expression of 12 genes involved in the late biosynthesis of ergosterol including ERG11 and ERG3. Importantly, chromatin immunoprecipitation analysis showed that histone H3K4 trimethylation was detected on chromatin of actively transcribed ERG genes. Furthermore, H3K4 trimethylation increased upon azole-induced gene expression which was also found to be dependent on the catalytic activity of Set1. Altogether, our findings show that Set1-mediated histone H3K4 methylation governs the intrinsic drug resistant status in C. glabrata via epigenetic control of azole-induced ERG gene expression.

IMPORTANCEC. glabrata is the second most commonly isolated species from Candida infections, coming in second to C. albicans. Treatment of C. glabrata infections are difficult due to their natural resistance to antifungal azole drugs and their ability to adapt and become multidrug resistant. In this study, we investigated the contributing cellular factors for controlling drug resistance. We have determined that an epigenetic mechanism governs the expression of genes involved in the late ergosterol biosynthesis pathway, an essential pathway that antifungal drugs target. This epigenetic mechanism involves histone H3K4 methylation catalyzed by the Set1 methyltransferase complex (COMPASS). We also show that Set1-mediated histone H3K4 methylation is needed for expression of specific azole induced genes and azole drug resistance in C. glabrata. Identifying epigenetic mechanisms contributing to drug resistance and pathogenesis could provide alternative targets for treating patients with fungal infections.


title: Mosaic loss of chromosome Y in aged human microglia.

abstract: Mosaic loss of chromosome Y (LOY) is a particularly common acquired structural mutation in the leukocytes of aging men and it has been shown to correlate with several age-related diseases including Alzheimers disease (AD). To derive the molecular basis of LOY in brain cells, we create an integrated resource by aggregating data from 21 single-cell and single-nuclei RNA brain studies, yielding 763,410 cells to investigate the presence and cell-type specific burden of LOY. We created robust quantification metrics for assessing LOY, which were validated using a multi-modal dataset. Using this new resource and LOY-quantification approach, we found that LOY frequencies differed widely between CNS cell-types and individual donors. Among five common neural cell types, microglia were most affected by LOY (7.79%, n=41,949), while LOY in neurons was rare (0.48%, n=220,010). Differential gene expression analysis in microglia found 188 autosomal genes, 6 X-linked genes, and 11 pseudoautosomal genes, pointing to broad dysregulation in lipoprotein metabolism, inflammatory response, and antigen processing that coincides with loss of Y. To our knowledge, we provide the first evidence of LOY in the microglia, and highlight its potential roles in aging and the pathogenesis of neurodegenerative disorders such as AD.


title: Shotgun metagenomics of soil invertebrate communities reflects taxonomy, biomass and reference genome properties

abstract:
O_LIMetagenomics - shotgun sequencing of all DNA fragments from a community DNA extract - is routinely used to describe the composition, structure and function of microorganism communities. Advances in DNA sequencing and the availability of genome databases increasingly allow the use of shotgun metagenomics on eukaryotic communities. Metagenomics offers major advances in the recovery of biomass relationships in a sample, in comparison to taxonomic marker gene based approaches (metabarcoding). However, little is known about the factors which influence metagenomics data from eukaryotic communities, such as differences among organism groups, the properties of reference genomes and genome assemblies.
C_LIO_LIWe evaluated how shotgun metagenomics records composition and biomass in artificial soil invertebrate communities. We generated mock communities of controlled biomass ratios from 28 species from all major soil mesofauna groups: mites, springtails, nematodes, tardigrades and potworms. We shotgun-sequenced these communities and taxonomically assigned them with a database of over 270 soil invertebrate genomes.
C_LIO_LIWe recovered 90% of the species, and observed relatively high false positive detection rates. We found strong differences in reads assigned to different taxa, with some groups (e.g. springtails) consistently attracting more hits than others (e.g. enchytraeids). Original biomass could be predicted from read counts after considering these taxon-specific differences. Species with larger genomes, and with more complete assemblies consistently attracted more reads than species with smaller genomes. The GC content of the genome assemblies had no effect on the biomass-read relationships.
C_LIO_LIThe results show considerable differences in taxon recovery and taxon specificity of biomass recovery from metagenomic sequence data. The properties of reference genomes and genome assemblies also influence biomass recovery, and they should be considered in metagenomic studies of eukaryotes. We provide a roadmap for investigating factors which influence metagenomics-based eukaryotic community reconstructions. Understanding these factors is timely as accessibility of DNA sequencing, and momentum for reference genomes projects show a future where the taxonomic assignment of DNA from any community sample becomes a reality.
C_LI


title: Sequence characterization of eccDNA content in glyphosate sensitive and resistant Palmer amaranth from geographically distant populations

abstract: The discovery of non-chromosomal circular DNA offers new directions in linking genome structure with function in plant biology. Glyphosate resistance through EPSPS gene copy amplification in Palmer amaranth was due to an autonomously replicating extra-chromosomal circular DNA mechanism (eccDNA). CIDER-Seq analysis of geographically distant glyphosate sensitive (GS) and resistant (GR) Palmer Amaranth (Amaranthus palmeri) revealed the presence of numerous small extra-chromosomal circular DNAs varying in size and with degrees of repetitive content, coding sequence, and motifs associated with autonomous replication. In GS biotypes, only a small portion of these aligned to the 399 kb eccDNA replicon, the vehicle underlying gene amplification and genetic resistance to the herbicide glyphosate. The aligned eccDNAs from GS were separated from one another by large gaps in sequence. In GR biotypes, the eccDNAs were present in both abundance and diversity to assemble into a nearly complete eccDNA replicon. Mean sizes of eccDNAs were similar in both biotypes and were around 5kb with larger eccDNAs near 25kb. Gene content for eccDNAs ranged from 0 to 3 with functions that include ribosomal proteins, transport, metabolism, and general stress response genetic elements. Repeat content among smaller eccDNAs indicate a potential for recombination into larger structures. Genomic hotspots were also identified in the Palmer amaranth genome with a disposition for gene focal amplifications as eccDNA. The presence of eccDNA may serve as a reservoir of genetic heterogeneity in this species and may be functionally important for survival.


title: The natural compound TMYX reduces SAN cells rate by antagonizing the cAMP modulation of f-channels

abstract:
Tongmai Yangxin (TMYX), is a complex compound of a Traditional Chinese Medicine (TCM) used to treat several cardiac rhythm disorders; however, no information regarding its mechanism of action is available. In this study we provide a detailed characterization of the effects of TMYX on the electrical activity of pacemaker cells and unravel its mechanism of action.

Single-cell electrophysiology revealed that TMYX elicits a reversible and dose-dependent (2/6 mg/ml) slowing of spontaneous action potentials rate (-20.8/-50.2%) by a selective reduction of the diastolic phase (-50.1/-76.0%). This action is mediated by a negative shift of the If activation curve (-6.7/-11.9 mV) and is caused by a reduction of the cAMP-induced stimulation of pacemaker channels. We provide evidence that TMYX acts by directly antagonizes the cAMP-induced allosteric modulation of the pacemaker channels. Noticeably, this mechanism functionally resembles the pharmacological actions of muscarinic stimulation or {beta}-blockers, but it does not require generalized changes in cytoplasmic cAMP levels thus ensuring a selective action on rate.

In agreement with a competitive inhibition mechanism, TMYX exerts its maximal antagonistic action at submaximal cAMP concentrations and then progressively becomes less effective thus ensuring a full contribution of If to pacemaker rate during high metabolic demand and sympathetic stimulation.

Funding sourcesThis work was supported by grants from Tianjin Zhongxin Pharmaceutical Group Co., Ltd. Le Ren Tang Pharmaceutical Factory P.R. China. The financial supporter played no role in the study design, data collection and analysis, decision to publish, or the preparation of the manuscript.