title: The balance of mitochondrial fission and fusion in cortical axons depends on the kinases SadA and SadB

abstract: Neurons are highly polarized cells that display characteristic differences in the organization of their organelles in axons and dendrites. Mitochondria are of particular importance for neuronal homeostasis due to their high metabolic demand. The kinases SadA and SadB (SadA/B) promote the formation of distinct axonal and dendritic extensions during the development of cortical and hippocampal neurons. Here, we show that SadA/B are required for the axon-specific dynamics of mitochondria. The interaction with Ankyrin B (AnkB) stimulates the activity of SadA/B that function as regulators of mitochondrial dynamics through the phosphorylation of Tau. Suppression of SadA/B or AnkB in cortical neurons induces the elongation of mitochondria by disrupting the balance of fission and fusion. The normal dynamics of axonal mitochondria could be restored by mild actin destabilization. Thus, the elongation after a loss of SadA/B results from an excessive stabilization of actin filaments and reduction of Drp1 recruitment to mitochondria.


title: A Reinforcement Meta-Learning Framework of Executive Function and Information Demand

abstract: Gathering information is crucial for maximizing fitness, but requires diverting resources from searching directly for primary rewards to actively exploring the environment. Optimal decision-making thus maximizes information while reducing effort costs, but little is known about the neural implementation of these tradeoffs. We present a Reinforcement Meta-Learning (RML) computational mechanism that solves the trade-offs between the value and costs of gathering information. We implement the RML in a biologically plausible architecture that links catecholaminergic neuromodulators, the medial prefrontal cortex and topographically organized visual maps and show that it accounts for neural and behavioral findings on information demand motivated by instrumental incentives and intrinsic utility. Moreover, the utility function used by the RML, encoded by dopamine, is an approximation of free-energy. Thus, the RML presents a biologically plausible mechanism through which coordinated motivational, executive and sensory systems generate visual information gathering policies that minimize free energy.

category: neuroscience


title: Online tests yield robust thresholds with the right auditory hygiene

abstract: Most human auditory psychophysics research has historically been conducted in carefully controlled environments with calibrated audio-equipment, and over potentially hours of repetitive testing with expert listeners. Here, we operationally define such conditions as having high  auditory hygiene. Hereof, conducting auditory psychophysical paradigms online presents a serious challenge, in that results may hinge on absolute sound presentation level, reliably estimated perceptual thresholds, and sustained motivation and attention. We introduce a set of procedures that address these challenges and facilitate auditory hygiene for online auditory psychophysics. First, we establish a simple means of setting sound presentation levels. Across a set of four level-setting conditions conducted in person, we demonstrate the stability and robustness of this amplitude setting procedure in open air and controlled settings. Second, we test participants tone-in-noise thresholds using widely adopted online experiment platforms and demonstrate that reliable threshold estimates can be derived online in approximately one minute of testing. Third, using these amplitude and threshold setting procedures to establish participant-specific stimulus conditions, we show that an online implementation of the classic probe-signal paradigm can be used to demonstrate frequency-selective attention on an individual participant basis, using a third of the trials used in recent in-lab experiments. Finally, we show how threshold and attentional measures relate to well-validated assays of online participants in-task motivation, fatigue, and confidence. This demonstrates the promise of online auditory psychophysics for asking new auditory neuroscience questions quickly, efficiently, and with more diverse samples. Code for the tests is publicly available through Pavlovia and Gorilla.


title: Robust and efficient online auditory psychophysics with the right auditory hygiene


title: Robust and efficient online auditory psychophysics


title: Pre-trauma memory contextualization as predictor for PTSD-like behavior in male rats

abstract: While many people experience potentially threatening events during their life, only a minority develops posttraumatic stress disorder (PTSD). The identification of individuals at risk among those exposed to trauma is crucial for PTSD prevention in the future. Since re-experiencing trauma elements outside of the original trauma-context is a core feature of PTSD, we investigate if the ability to bind memories to their original encoding context (i.e. memory contextualisation) predicts PTSD vulnerability. We hypothesize that pre-trauma neutral memory contextualization (under stress) negatively relates to PTSD-like behavior, in a prospective design using the cut-off behavioral criteria rat model for PTSD. 72 male Sprague Dawley rats were divided in two experimental groups to assess the predictive value of 1) memory contextualization without acute stress (NS-group) and 2) memory contextualization during the recovery phase of the acute stress-response (S-group) for susceptibility to PTSD-like behavior. A powerful extension to regression analysis-path analysis-was used to test this specific hypothesis, together with secondary research questions. Following traumatic predator scent stress, 19.4% of the rats displayed PTSD-like behavior. Results showed a negative relation between pre-trauma memory contextualization and PTSD-like behavior, but only in the NS-group. Pre-trauma memory contextualization was positively related to fear association in the trauma environment, again only in the NS group. If the predictive value of pre-trauma contextualization of neutral information under non-stressful conditions for PTSD susceptibility is replicated in prospective studies in humans, this factor would supplement already known vulnerability factors for PTSD and improve the identification of individuals at risk among the trauma exposed, especially those at high trauma risk such as soldiers deployed on a mission.


title: Axial spondyloarthritis patients have altered mucosal IgA response to oral and fecal microbiota

abstract:
ObjectiveTo investigate whether axial spondyloarthritis (AxSpA) patients have an altered immunoglobulin A (IgA) response in the gut and oral microbial communities.

MethodsWe performed 16S rRNA gene (16S) sequencing on IgA positive (IgA+) and IgA negative (IgA-) fractions (IgA-SEQ) from feces (n=17 AxSpA; n=14 healthy) and saliva (n=17 AxSpA; n=12 healthy), as well as on IgA-unsorted fecal and salivary samples. PICRUSt2 was used to predict microbial metabolic potential in AxSpA patients and healthy controls (HCs).

ResultsIgA-SEQ revealed enrichment of several microbes in the fecal (Akkermansia, Ruminococcaceae, Lachnospira) and salivary (Prevotellaceae, Actinobacillus) microbiome in AxSpA patients as compared with HCs. Fecal microbiome from AxSpA patients showed a trend towards increased alpha diversity of the IgA+ fraction and decreased diversity in the IgA- fraction in comparison with HCs, while the salivary microbiome exhibits a significant decrease in alpha diversity in both IgA+ and IgA- fractions. Increased IgA coating of Clostridiales Family XIII correlated with disease severity. Inferred metagenomic analysis suggests perturbation of metabolites and metabolic pathways for inflammation (oxidative phosphorylation, glutathione metabolism) and metabolism (propanoate and butanoate metabolism) in AxSpA patients.

ConclusionsAnalyses of fecal and salivary microbes from AxSpA patients reveal distinct populations of immunoreactive microbes using novel IgA-SEQ approach, which were not captured by comparing their relative abundance with HCs. Predictive metagenomic analysis revealed perturbation of metabolites/metabolic pathways in AxSpA patients. Future studies on these immunoreactive microbes may lead to better understanding of the functional role of IgA in maintaining microbial structure and human health.

category: immunology


title: Echidna: integrated simulations of single-cell immune receptor repertoires and transcriptomes

abstract:
Single-cell sequencing now enables the recovery of full-length immune repertoires [B cell receptor (BCR) and T cell receptor (TCR) repertoires], in addition to gene expression information. The feature-rich datasets produced from such experiments require extensive and diverse computational analyses, each of which can significantly influence the downstream immunological interpretations, such as clonal selection and expansion. Simulations produce validated standard datasets, where the underlying generative model can be precisely defined and furthermore perturbed to investigate specific questions of interest. Currently, there is no tool that can be used to simulate a comprehensive ground truth single-cell dataset that incorporates both immune receptor repertoires and gene expression. Therefore, we developed Echidna, an R package that simulates immune receptors and transcriptomes at single-cell resolution. Our simulation tool generates annotated single-cell sequencing data with user-tunable parameters controlling a wide range of features such as clonal expansion, germline gene usage, somatic hypermutation, and transcriptional phenotypes. Echidna can additionally simulate time-resolved B cell evolution, producing mutational networks with complex selection histories incorporating class-switching and B cell subtype information. Finally, we demonstrate the benchmarking potential of Echidna by simulating clonal lineages and comparing the known simulated networks with those inferred from only the BCR sequences as input. Together, Echidna provides a framework that can incorporate experimental data to simulate single-cell immune repertoires to aid software development and bioinformatic benchmarking of clonotyping, phylogenetics, transcriptomics and machine learning strategies.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=60 SRC="FIGDIR/small/452792v1_ufig1.gif" ALT="Figure 1">
View larger version (19K):
org.highwire.dtl.DTLVardef@1c9930eorg.highwire.dtl.DTLVardef@d5a50eorg.highwire.dtl.DTLVardef@47584forg.highwire.dtl.DTLVardef@1f6bb17_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: The age and mouse sperm quality: a flow cytometry investigation

abstract: Postponement of fatherhood is growing worldwide due to socio-economic factors. The choice to conceive the first child above the age of 35 years is often associated with reduced fertility and poor pregnancy outcome. As widely known, several factors (e.g., lifestyle, environment, health problems) can affect spermatogenesis leading to poor reproductive outcome. Currently, the debate on the influence of aging on male gametes and safety/risk of conception at advanced age is still ongoing. Controversial results have been published so far on the changes in semen features of aging men and other mammalian species (mainly rodents). In this study, we aimed to assess how aging affects sperm quality in an inbreed mouse model, without underlying infertility, using a flow cytometry approach. Our data showed that aging is associated with increased sperm chromatin condensation, but not changes in the DNA integrity, metabolic activity or viability. These data suggest a mild effect of aging on sperm quality in a mouse model without underlying infertility.

category: developmental biology


title: Coarse-graining ecological dynamics in the face of unknown microscopic details

abstract: Despite their complexity, microbial ecosystems appear to be at least partially "coarse-grainable" in that some properties of interest can be adequately described by effective models of dimension much smaller than the number of interacting lineages (frequently in the dozens or hundreds). This is especially puzzling since recent studies demonstrate that a surprising amount of functionally relevant diversity is present at all levels of resolution, down to strains differing by 100 nucleotides or fewer. Rigorously defining coarse-grainability and understanding the conditions for its emergence is of critical importance for understanding microbial ecosystems. To begin addressing these questions, we propose a minimal model for investigating hierarchically structured ecosystems within the framework of resource competition. We use our model to operationally define coarse-graining quality based on reproducibility of the outcomes of a specified experiment and show that a coarse-graining can be operationally valid despite grouping together functionally diverse strains. We further show that, at least within our model, a high diversity of strains (while nominally more complex) may in fact facilitate coarse-grainability. However, this only applies if the pool of interacting strains is sampled from the appropriate "native" environment, as we discuss.

category: ecology


title: Defining coarse-grainability in a model of structured microbial ecosystems


title: Genetically-based adaptive trait shifts at an expanding mangrove range margin

abstract:
O_LIMany species are expanding beyond their distributional range margins in response to a warming planet. Due to marginal environmental conditions and novel selection pressures, range margins may foster unique genetic adaptations that can better enable species to thrive under the extreme climatic conditions at and beyond their current distributional limits. Neotropical black mangrove (Avicennia germinans) is expanding poleward into temperate salt marsh along Atlantic Florida, USA, with field evidence of adaptive trait shifts within range-margin A. germinans populations. However, whether these adaptive shifts have a genetic basis remains to be answered.
C_LIO_LIWe monitored twenty A. germinans maternal cohorts from areas in both the Atlantic Florida range core and margin in a greenhouse common garden with annual temperatures analogous to range-margin conditions. We measured variation in a series of phenotypic traits starting at initial planting of field-collected propagules and continuing until two years development.
C_LIO_LIMaternal cohorts from the Atlantic Florida range margin consistently outperformed those from the range core throughout the experiment. Range-margin cohorts survived in greater numbers, established faster, and were less stressed under winter chilling and sub-zero temperatures that are often reached at the Atlantic range margin, but not within the range core. Range-margin cohorts did not grow taller, but instead invested more into lateral growth and biomass accumulation that presumably reflects adaptation to their colder and open canopy environment. Range-margin cohorts also exhibited leaf traits consistent with greater resource acquisition that may compensate for a shorter growing season and reduced light quality at higher latitude.
C_LIO_LISynthesis. We confirmed that there is a genetic basis to adaptive trait shifts towards an expanding mangrove range margin. Our results suggest that genetically-based phenotypic differences better enable these range-margin mangroves to thrive within their stressful environment and may facilitate further poleward expansion in the future. In addition, our documentation of adaptive trait variation among maternal cohorts of an ecologically-important mangrove foundation species, quantitative data that is lacking for mangroves, should help inform mangrove restoration initiatives.
C_LI


title: Gene regulation is commonly selected for high plasticity and low noise

abstract: Bacteria often respond to dynamically changing environments by regulating gene expression. Despite this regulation being critically important for growth and survival, little is known about how selection shapes gene regulation in natural populations. To better understand the role natural selection plays in shaping bacterial gene regulation, here we compare differences in the regulatory behaviour of naturally segregating promoter variants from Escherichia coli (which have been subject to natural selection) to randomly mutated promoter variants (which have never been exposed to natural selection). We quantify gene expression phenotypes (expression level, plasticity, and noise) for hundreds of promoter variants across multiple environments, and show that segregating promoter variants are enriched for mutations with minimal effects on expression level. In many promoters, we infer that there is strong selection to maintain high levels of plasticity, and direct selection to decrease or increase cell-to-cell variability in expression. Finally, taking an integrated view, we show that across all phenotypes combined, segregating promoter variants are far more phenotypically similar than would be expected given their genetic divergence. This is the consequence of both stabilizing and directional selection acting on individual phenotypes to minimize differences among segregating variants. Taken together, these results expand our knowledge of how gene regulation is affected by natural selection and highlight the power of comparing naturally segregating polymorphisms to de novo random mutations to quantify the action of selection.

category: evolutionary biology


title: Discovering cellular programs of intrinsic and extrinsic drivers of metabolic traits using LipocyteProfiler

abstract: A primary obstacle in translating genetics and genomics data into therapeutic strategies is elucidating the cellular programs affected by genetic variants and genes associated with human diseases. Broadly applicable high-throughput, unbiased assays offer a path to rapidly characterize gene and variant function and thus illuminate disease mechanisms. Here, we report LipocyteProfiler, an unbiased high-throughput, high-content microscopy assay that is amenable to large-scale morphological and cellular profiling of lipid-accumulating cell types. We apply LipocyteProfiler to adipocytes and hepatocytes and demonstrate its ability to survey diverse cellular mechanisms by generating rich context-, and process-specific morphological and cellular profiles. We then use LipocyteProfiler to identify known and novel cellular programs altered by polygenic risk of metabolic disease, including insulin resistance, waist-to-hip ratio and the polygenic contribution to lipodystrophy. LipocyteProfiler paves the way for large-scale forward and reverse phenotypic profiling in lipid-storing cells, and provides a framework for the unbiased identification of causal relationships between genetic variants and cellular programs relevant to human disease.

category: genomics


title: Darwin's finches can have taste preferences, but not for human foods

abstract: Urbanization is rapidly changing ecological niches. On the inhabited Galapagos Islands, Darwins finches consume human-introduced foods preferentially; however, it remains unclear why. Here we presented pastry with flavour profiles typical of human foods (oily, salty, sweet) to small ground finches (Geospiza fuliginosa) and medium ground finches (Geospiza fortis) to test if latent taste preferences might drive selection of human foods. If human-food flavours were consumed more than a neutral or bitter control only at sites with human foods, then we predicted tastes were acquired after urbanization; however, if no site-differences were found then this would indicate latent taste preferences. Contrary to both predictions, we found little evidence that human-food flavours were preferred compared to control flavours at any site. Instead, finches showed a weak aversion to oily foods, but only at remote (no human foods present) sites. This was further supported by behavioural responses, with beak wiping occurring more often at remote sites after finches tasted flavours associated with human foods. Our results suggest, therefore, that while Darwins finches may have acquired a tolerance to the flavours of human food, latent taste preferences are unlikely to have played a major role in their dietary response to increased urbanization.

category: animal behavior and cognition


title: Darwin's small and medium ground finches might have taste preferences, but not for human foods


title: Experimental investigations of the human oesophagus: anisotropic properties of the muscular layer in large deformation

abstract: Technological advancements in the field of robotics have led to endoscopic biopsy devices able to extract diseased tissue from between the layers of the gastrointestinal tract. Despite this, the layer-dependent properties of these tissues have yet to be mechanically characterised using human tissue. In this study, the ex vivo mechanical properties of the passive muscularis propia layer of the human oesophagus were extensively investigated. For this, a series of uniaxial tensile tests were conducted. The results displayed hyperelastic behaviour, while the differences between loading the tissue in both the longitudinal and circumferential directions showcased its anisotropy. The anisotropy of the muscular layer was present at different strain rates, with the longitudinal direction being consistently stiffer than the circumferential one. The circumferential direction was found to have little strain-rate dependency, while the longitudinal direction results suggest pronounced strain-rate-dependent behaviour. The repeated trials showed larger variation in terms of stress for a given strain in the longitudinal direction compared to the circumferential direction. The possible causes of variation between trials are discussed, and the experimental findings are linked to the histological analysis which was carried out via various staining methods. Finally, the direction-dependent experimental data was simulated using an anisotropic, hyperelastic model.

category: bioengineering


title: HAMdetector: A Bayesian regression model that integrates information to detect HLA-associated mutations

abstract:
MotivationA key process in anti-viral adaptive immunity is that the Human Leukocyte Antigen system (HLA) presents epitopes as Major Histocompatibility Complex I (MHC I) protein-peptide complexes on cell surfaces and in this way alerts CD8+ cytotoxic T-Lymphocytes (CTLs). This pathway exerts strong selection pressure on viruses, favoring viral mutants that escape recognition by the HLA/CTL system, e.g. by point mutations that decrease binding of viral peptides to MHC I. Naturally, such immune escape mutations often emerge in highly variable viruses, e.g. HIV or HBV, as HLA-associated mutations (HAMs), specific to the host HLA alleles and its MHC I proteins. The reliable identification of HAMs is not only important for understanding viral genomes and their evolution, but it also impacts the development of broadly effective anti-viral treatments and vaccines against variable viruses.

By their very nature HAMs are amenable to detection by statistical methods in paired sequence / HLA data. However, HLA alleles are very polymorphic in the human host population which makes the available data relatively sparse and noisy. Under these circumstances, one way to optimize HAM detection is to integrate all relevant information in a coherent model. Bayesian inference offers a principled approach to achieve this.

ResultsWe present a new regression model for the detection of HAMs. As we choose a Bayesian approach we can include the novel sparsity-inducing priors, and we obtain easily interpretable quantitative information on HAM candidates. The basic model can be extended to include prior information relevant to HAM detection, which we demonstrate by integrating predictions of epitope affinities to MHC I, predictions of epitope peptide processing, and computation of phylogenetic background. This integrative method improves performance in HAM detection considerably over state-of-the-art methods.

AvailabilityThe source code of this software is available at https://github.com/HAMdetector/Escape.jl under a permissive MIT license.

Contactdaniel.habermann@uni-due.de, daniel.hoffmann@uni-due.de

category: bioinformatics


title: Causal inference using deep-learning variable selection identifies and incorporates direct and indirect causalities in complex biological systems

abstract: In complex diseases, causal structure learning across biological variables is critical to identify modifiable triggers or potential therapeutic agents. A limitation of existing causal learning methods is that they cannot identify indirect causal relations, those that would interact through latent mediating variables. We developed the first computational method that identifies both direct and indirect causalities, causal inference using deep-learning variable-selection (causalDeepVASE). To accurately identify indirect causalities and incorporate them with direct causalities, causalDeepVASE develops a deep neural network approach and extends a flexible causal inference method. In simulated and biological data of various contexts, causalDeepVASE outperforms existing methods in identifying expected or validated causal relations. Further, causalDeepVASE facilitates a systematic understanding of complex diseases. For example, causalDeepVASE uniquely identified a possible causal relation between IFN{gamma} and creatinine suggested in a polymicrobial sepsis model. In future biomedical studies, causalDeepVASE can facilitate the identification of driver genes and therapeutic agents.


title: Open access in silico tools to predict the ADMET profiling and PASS (Prediction of Activity Spectra for Substances of Bioactive compounds of Garlic (Allium sativum L.)

abstract:
BackgroundGarlic (Allium sativum L.) is a common spice with many health benefits, mainly due to its diverse bioactive compounds, (see below) such as organic sulphides, saponins, phenolic compounds, and polysaccharides. Several studies have demonstrated its functions such as anti-inflammatory, antibacterial, and antiviral, antioxidant, cardiovascular protective and anticancer property. In this work we have investigated the main bioactive components of garlic through a bioinformatics approach. Indeed, we are in an era of bioinformatics where we can predict data in the fields of medicine. Approaches with open access in silico tools have revolutionized disease management due to early prediction of the absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiles of the chemically designed and eco-friendly next-generation drugs.

MethodsThis paper encompasses the fundamental functions of open access in silico prediction tools, as PASS database (Prediction of Activity Spectra for Substances) that it estimates the probable biological activity profiles for compounds. This paper also aims to help support new researchers in the field of drug design and to investigate best bioactive compounds in garlic.

Resultsscreening through each of pharmacokinetic criteria resulted in identification of Garlic compounds that adhere to all the ADMET properties.

ConclusionsIt was established an open-access database (PASS database, available bioinformatics tool SwissADME, PreADMET pkCSM database) servers were employed to determine the ADMET (metabolism, distribution, excretion, absorption, and toxicity) attributes of garlic molecules and to enable identification of promising molecules that follow ADMET properties.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=113 SRC="FIGDIR/small/452815v1_ufig1.gif" ALT="Figure 1">
View larger version (53K):
org.highwire.dtl.DTLVardef@1344e43org.highwire.dtl.DTLVardef@fe1674org.highwire.dtl.DTLVardef@1742577org.highwire.dtl.DTLVardef@47a966_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: Sex-Specific Computational Models of Kidney Function in Human Patients with Diabetes

abstract: The kidney plays an essential role in homeostasis, accomplished through the regulation of pH, electrolytes and fluids, by the building blocks of the kidney, the nephrons. One of the important markers of the proper functioning of a kidney is the glomerular filtration rate. Diabetes is characterised by an enlargement of the glomerular and tubular size of the kidney, affecting the afferent and efferent arteriole resistance and hemodynamics, ultimately leading to chronic kidney disease. We postulate that the diabetes-induced changes in kidney may exhibit significant sex differences as the distribution of renal transporters along the nephron may be markedly different between women and men, as recently shown in rodents. The goals of this study are to (i) analyze how kidney function is altered in male and female patients with diabetes, and (ii) assess the renal effects, in women and men, of an anti-hyperglycemic therapy that inhibits the sodium-glucose cotransporter 2 (SGLT2) in the proximal convoluted tubules. To accomplish these goals, we have developed computational models of kidney function, separate for male and female patients with diabetes. The simulation results indicate that diabetes enhances Na+ transport, especially along the proximal tubules and thick ascending limbs, to similar extents in male and female patients, which can be explained by the diabetes-induced increase in glomerular filtration rate. Additionally, we conducted simulations to study the effects of diabetes and SGLT2 inhibition on solute and water transport along the nephrons. Model simulations also suggest that SGLT2 inhibition, which constricts the afferent arteriole to attenuate glomerular hyperfiltration, can then limit Na+-glucose transport, consequently raising luminal [Cl-] at the macula densa and finally restoring the tubuloglomerular feedback signal. By inducing osmotic diuresis in the proximal tubules, SGLT2 inhibition reduces paracellular transport, eventually leading to diuresis and natriuresis. Those effects on urinary excretion are blunted in women, in part due to their higher distal transport capacity.

category: physiology


title: Lane-maze for preference testing in flies

abstract: Drosophila melanogaster is a candidate species to replace rodents in some neurobiological studies, encouraging attempts to develop behavioural tests for these flies. This study aimed to develop a behavioural test to simultaneously evaluate ethological (categorical) aspects of the motor and fluid intake activities, which may be used to assess sucrose preference in flies. For that, a lane-maze was 3D-printed to accommodate up to 14 individual flies in a single trial. Each lane had a capillary filled with 5% sucrose solution attached to one of the extremities. To validate a 5-min lane-maze test, male and female flies (adults, 5-6 days of age) underwent 0, 2, 8 or 20 h of food deprivation (FD, n=9-11/group) before testing. Duration of locomotion, immobility and grooming in the lane or capillary were scored from the video-recorded trials using EthoWatcher software. Minor effects of sex or FD were observed in the behaviours of flies. Independent of sex or FD, flies spent proportionally longer on the capillary than on the lane. Flies exhibited a significantly higher preference than expected for the capillary zone when food-deprived for 2h (males) or 20 h (females). Data suggest that short lane-maze test is a feasibly high throughput assessment of sucrose preference in flies, which may be sexually dimorphic as in other species studied so far.


title: Towards on-chip real-time classification of extra-cellular neural recordings

abstract: On-line classification of neural recordings can be extremely useful in brain-machine interface, prosthetic applications or therapeutic intervention. In this work we present a feasibility study for developing compact low-power VLSI systems able to classify neural recordings in real-time, using spike-based neuromorphic circuits. We developed a framework for classifying extra-cellular recordings made in rat auditory cortex in response to different auditory stimuli and porting the classification algorithm onto a spiking multi-neuron VLSI chip with programmable synaptic weights. We present recording methods and software classification algorithms; we demonstrate real-time classification in hardware and quantify the system performance; finally, we identify the potential sources of problems in developing such types of systems and propose strategies for overcoming them.


title: Traces of semantization - from episodic to semantic memory in a spiking cortical network model

abstract: Episodic memory is the recollection of past personal experiences associated with particular times and places. This kind of memory is commonly subject to loss of contextual information or "semantization", which gradually decouples the encoded memory items from their associated contexts while transforming them into semantic or gist-like representations. Novel extensions to the classical Remember/Know behavioral paradigm attribute the loss of episodicity to multiple exposures of an item in different contexts. Despite recent advancements explaining semantization at a behavioral level, the underlying neural mechanisms remain poorly understood. In this study, we suggest and evaluate a novel hypothesis proposing that Bayesian-Hebbian synaptic plasticity mechanisms might cause semantization of episodic memory. We implement a cortical spiking neural network model with a Bayesian-Hebbian learning rule called Bayesian Confidence Propagation Neural Network (BCPNN), which captures the semantization phenomenon and offers a mechanistic explanation for it. Encoding items across multiple contexts leads to item-context decoupling akin to semantization. We compare BCPNN plasticity with the more commonly used spike-timing dependent plasticity (STDP) learning rule in the same episodic memory task. Unlike BCPNN, STDP does not explain the decontextualization process. We also examine how selective plasticity modulation of isolated salient events may enhance preferential retention and resistance to semantization. Our model reproduces important features of episodicity on behavioral timescales under various biological constraints whilst also offering a novel neural and synaptic explanation for semantization, thereby casting new light on the interplay between episodic and semantic memory processes.


title: Decoding molecular and cellular heterogeneity of nucleus accumbens with high-throughput scRNA-seq and MERFISH

abstract: The nucleus accumbens (NAc) plays an important role in regulating multiple behaviors and its dysfunction has been linked to many neural disorders. However, the molecular, cellular and anatomic heterogeneity underlying its functional diversity remains incompletely understood. Here, we generate a cell census of the mouse NAc using high-throughput single cell RNA sequencing and multiplexed error-robust FISH, revealing a high level of cell heterogeneity in this brain region. We show that the transcriptional and spatial diversity of neuron subtypes underlie NAcs anatomic and functional heterogeneity, and possibly contribute to the pathogenesis of different neurological disorders. These findings explain how the seemingly simple neuronal composition of the NAc achieves its highly heterogenous structure and diverse functions. Collectively, our study generates a spatially resolved cell taxonomy for understanding the NAc structure and function, which demonstrates the importance of combining molecular and spatial information in revealing the fundamental features of the nervous system.


title: Isotropic Reconstruction of Electron Tomograms with Deep Learning

abstract: Cryogenic electron tomography (cryoET) allows visualization of cellular structures in situ. However, anisotropic resolution arising from the intrinsic "missing-wedge" problem has presented major challenges in visualization and interpretation of tomograms. Here, we have developed IsoNet, a deep learning-based software package that iteratively reconstructs the missing-wedge information and increases signal-to-noise ratio, using the knowledge learned from raw tomograms. Without the need for sub-tomogram averaging, Isonet generates tomograms with significantly reduced resolution anisotropy. Applications of IsoNet to three representative types of cryoET data demonstrate greatly improved structural interpretability: resolving lattice defects in immature HIV particles, establishing architecture of the paraflagellar rod in Eukaryotic flagella, and identifying heptagon-containing clathrin cages inside a neuronal synapse of cultured cells. Therefore, by overcoming two fundamental limitations of cryoET, IsoNet enables functional interpretation of cellular tomograms without sub-tomogram averaging. Its application to high-resolution cellular tomograms should also help identify differently oriented complexes of the same kind for near-atomic resolution sub-tomogram averaging.

category: cell biology


title: Dosage compensation in Bombyx mori is achieved by partial repression of both Z chromosomes in males

abstract:
Interphase chromatin is organized precisely to facilitate accurate gene expression. The structure-function relationship of chromatin is epitomized in sex chromosome dosage compensation (DC), where sex-linked gene expression is balanced between males and females via sex-specific alterations to 3D chromosome structure. Studies in ZW-bearing species suggest that DC is absent or incomplete in most lineages except butterflies and moths, where male (ZZ) chZ expression is reduced by half to equal females (ZW). However, whether one chZ is inactivated (as in mammals) or both are partially repressed (as in C. elegans) is unknown. Using Oligopaints in the silkworm, Bombyx mori, we visualize autosome and chZ organization in somatic cells from both sexes for the first time. We find that B. mori interphase chromosomes are highly compact relative to Drosophila chromosomes. Importantly, we show that in B. mori males, both chZs are similar in size and shape and are more compact than autosomes or the female chZ after DC establishment, suggesting that both male chZs are partially and equally downregulated. We also find that in the early stages of DC, the female chZ repositions toward the nuclear center concomitant with increased Z-linked gene expression, revealing the first non-sequencing-based support for Ohnos hypothesis. These studies represent the first visualization of interphase genome organization and chZ structure in Lepidoptera. We uncover striking similarities between DC in B. mori and C. elegans, despite these lineages harboring evolutionarily distinct sex chromosomes (ZW/XY), suggesting convergent evolution of DC mechanisms and a possible role for holocentricity in DC evolution.

SignificanceGenes on sex chromosomes (for example, the X chromosome in humans) are regulated by distinct processes that do not affect non-sex chromosomes (autosomes). The expression of genes of sex chromosomes needs to be balanced between males and females due to differences in sex chromosome dosage (XY versus XX). This sex-specific gene regulation is called dosage compensation (DC). In most species, DC is achieved by altering the shape and compaction of sex chromosomes specifically in one sex. In this study, we use a chromosome painting approach called Oligopaints to examine how DC is achieved in silkworms for the first time. This is the first study to visualize this phenomenon in a species with ZW sex chromosomes, which evolved completely independently of XY. We find evidence supporting a long-standing model for how DC mechanisms evolved across all organisms, and we show high similarity between DC in silkworms and nematodes, suggesting these mechanisms have emerged multiple independent times throughout evolution.

category: genetics


title: Structural Characterization of the PawL-Derived Peptide Family, an Ancient Subfamily of Orbitides

abstract: Plants are an excellent source of bioactive peptides, often with disulfide bonds and/or a cyclic backbone. While focus has predominantly been directed at disulfide-rich peptides, a large family of small, cyclic but non-disulfide bonded plant peptides, known as orbitides, has been relatively ignored. A recently discovered subfamily of orbitides is the PawL-Derived Peptides (PLPs), produced during the maturation of precursors for seed storage albumins. Although their evolutionary origins have been dated, in-depth exploration of the familys structural characteristics and potential bioactivities remains to be conducted. Here we present an extensive and systematic characterization of the PLP family. Nine PLPs were chosen and prepared by solid phase peptide synthesis. Their structural features were studied using solution NMR spectroscopy and seven were found to possess regions of backbone order. Ordered regions consist of {beta}-turns, with some PLPs adopting two well-defined {beta}-turns within sequences as short as seven residues, which are largely the result of side chain interactions. Our data highlight that the sequence diversity within this family results in equally diverse molecular scaffolds. None of these nine PLPs showed antibacterial or antifungal activity.

category: biochemistry


title: PLIGHT: A tool to assess privacy risk by inferring identifying characteristics from sparse, noisy genotypes

abstract: The leakage of identifying information in genetic and omics data has been established in many studies, with single nucleotide polymorphisms (SNPs) shown to carry a strong risk of reidentification for individuals and their genetic relatives. While the ability of thousands or hundreds of thousands of SNPs (especially rare ones) to identify individuals has been demonstrated, here we sought to measure the informativeness of even a sparse set of tens of noisy, common SNPs from an individual, by putting the genotype-based privacy leakage from an individual on quantitative footing. We present a computational tool, PLIGHT ("Privacy Leakage by Inference across Genotypic HMM Trajectories"), that employs a population-genetics-based Hidden Markov Model of recombination and mutation to find piecewise matches of a sparse query set of SNPs to a reference genotype panel. Given the ready availability of auxiliary sources of noisy genotype data - such as acquiring small samples of environmental DNA or learning about someones Mendelian diseases and physical characteristics - inference on sparse data becomes a genuine concern. We explore cases where query individuals are either known to be in databases or not, and consider both simulated "mosaics" of genotypes (i.e. genotypes stitched together from diploid segments sampled from two or more source individuals) and actual genotype data obtained from swabs of coffee cups used by a known individual. Our findings are as follows: (1) Even 10 common SNPs (minor allele frequency > 0.05) often are sufficient to identify individuals in conventional genomic databases. (2) We are able to identify first-order relatives (parents, children and siblings) of query individuals with 20-30 common SNPs. (3) We find some potential for leakage of phenotypic information, based on a simulated attack by combining polygenic risk scores (PRSs) of the piecewise genotypic matches. We also found, for simulated mosaics of two individuals, that 20 common SNPs were often sufficient to find the correct identities of both component individuals. Finally, applying PLIGHT to coffee-cup-derived SNPs, we find that our tool is able identify the individual (when present in the reference database) using as little as 30 SNPs; alternatively, when the individual is not present in the reference database, we reconstruct possible genomes for the individual based on just 30-90 query SNPs by piecewise matching to the reference haplotype database. In this way, we are able to perform a small degree of imputation of unobserved query SNPs. Overall, the tool could be used to determine the value of selectively masking released SNPs, in a way that is agnostic to any explicit assumptions about underlying population membership or allele frequencies.


title: Docking Molecular analysis of potential Aldosterone antagonists

abstract: BackgroundAldosterone antagonists (spironolactone, eplerenone) inhibit the action of aldosterone in the collecting duct; as such, these agents cause modest diuresis but inhibit potassium and hydrogen ion secretion. We report first time Potential Aldosterone antagonists by in Silico approach, using AutoDock Vina and AutoDock 4 (or MGL Tool), estimated with Pyrx and AM Dock Software, calculating three different important parameters: Binding Affinity (kcal/mol), estimated Ki (in nM units) and Ligand Efficiency (L.E. in kcal/mol). After a selective analysis of over 1000 drugs, processed with Pyrx (a Virtual Screening software for Computational Drug Discovery) in the Ligand Binding site pocket of the protein (ID PDB 2OAX Chain A:), we noticed high values of Binding Energy, about -13.55 kcal/mol estimated by AutoDock 4 with AM Dock Software, concluding that it could be an excellent candidate drug, compared to everyone else Aldosterone antagonists. Indeed, from the results of AutoDock Vina and AutoDock 4 (or AutoDock 4.2), implemented with Lamarckian genetic algorithm, LGA, trough AMDock Software, our results of Binding Energy are very similar to the crystallized Spironolactone in PDB 2OAX Chain A protein.


title: Application of personalized differential expression analysis in human cancer proteome

abstract: Owing to the recent technological advances, liquid chromatography-mass spectrometry (LC-MS)-based quantitative proteomics can measure expression of thousands of proteins from biological specimens. Currently, several studies have used the LC-MS-based proteomics to measure protein expression levels in human cancer. Identifying differentially expressed proteins (DEPs) between tumors and normal controls is a common way to investigate carcinogenesis mechanisms. However, most statistical methods used for DEPs analysis can only identify deregulated proteins at the population-level and ignore the heterogeneous differential expression of proteins in individual patients. Thus, to identify patient-specific molecular defects for personalized medicine, it is necessary to perform personalized differential analysis at the scale of a single sample. To date, there is a scarcity of systematic and easy-to-handle tool that could be used to evaluate the performance of individualized difference expression analysis algorithms in human cancer proteome. Herein, we developed a user-friendly tool kit, IDEP, to enable implementation and evaluation of personalized difference expression analysis algorithms. IDEP evaluates five rank-based tools (RankComp v1/v2, PENDA, Peng and Quantile) through classic computational and functional criteria in lung, gastric and liver cancer proteome. The results show that the within-sample relative expression orderings (REOs) of protein pairs in normal tissues were highly stable, which provided the basis for individual level DEPs analysis. Moreover, these individualized difference analysis tools could reach much higher efficiency in detecting sample-specific deregulated proteins than the group-based methods. Pathway enrichment and survival analysis results were dataset and analysis method dependent. In summary, IDEP has integrated necessary toolkits for individualized identification of DEPs and supported flexible methods evaluation analysis and visualization modules. It could provide a robust and scalable framework to extract personalized deregulation patterns and could also be used for the discovery of prognostic biomarkers for personalized medicine.


title: CoLoRd: Compressing long reads

abstract: The costs of maintaining exabytes of data produced by sequencing experiments every year has become a major issue in todays genomics. In spite of the increasing popularity of the third generation sequencing, the existing algorithms for compressing long reads exhibit minor advantage over general purpose gzip. We present CoLoRd, an algorithm able to reduce 3rd generation sequencing data by an order of magnitude without affecting the accuracy of downstream analyzes.


title: Prior Signal Acquisition Software Versions for Orbitrap Underestimate Low Isobaric Mass Tag Intensities, Without Detriment to Differential Abundance Experiments

abstract: Tandem mass tags (TMT) enable simple and accurate quantitative proteomics for multiplexed samples by relative quantification of tag reporter ions. Orbitrap quantification of reporter ions has been associated with a characteristic notch region in intensity distribution, within which few reporter intensities are recorded. This has been resolved in version 3 of the instrument acquisition software, Tune. However, 53 % of Orbitrap Fusion, Lumos or Eclipse submissions to PRIDE were generated using prior software versions. To quantify the impact of the notch on existing quantitative proteomics data, we generated a mixed species benchmark and acquired quantitative data using Tune versions 2 and 3. Sub-notch intensities are systemically underestimated with Tune version 2, leading to over-estimation of the true differences in intensities between samples. However, when summarising reporter ion intensities to higher level features, such as peptides and proteins, few features are significantly affected. Targeted removal of spectra with reporter ion intensities below the notch is not beneficial for differential peptide or protein testing. Overall, we find the systematic quantification bias associated with the notch is not detrimental for a typical proteomics experiment.


title: Pleiotropic Variability Score: A Genome Interpretation Metric to Quantify Phenomic Associations of Genomic Variants

abstract: A more complete understanding of phenomic space is critical for elucidating genome-phenome relationships and for assessing disease risk from genome sequencing. To incorporate knowledge of how related a variants associations are, we developed a new genome interpretation metric called Pleiotropic Variability Score (PVS). PVS uses semantic reasoning to score the relatedness of a genetic variants associated phenotypes based on those phenotypes relationships in the human phenotype ontology (HPO) and disease ontology (DO). We tested 78 unique semantic similarity methods and integrated six robust metrics to define the pleiotropy score of SNPs. We computed PVS for 12,541 SNPs which were mapped to 382 HPO and 317 DO unique phenotype terms in a genotype-phenotype catalog (10,021 SNPs mapped to DO phenotypes and 8,569 SNPs mapped to HPO phenotypes). We validated the utility of PVS by computing pleiotropy using an electronic health record linked genomic database (BioME, n=11,210). Further we demonstrate the application of PVS in personalized medicine using "personalized pleiotropy score" reports for individuals with genomic data that could potentially aid in variant interpretation. We further developed a software framework to incorporate PVS into VCF files and to consolidate pleiotropy assessment as part of genome interpretation pipelines. As the genome-phenome catalogs are growing, PVS will be a useful metric to assess genetic variation to find SNPs with highly pleiotropic effects. Additionally, variants with varying degree of pleiotropy can be prioritized for explorative studies to understand specific roles of SNPs and pleiotropic hubs in mediating novel phenotypes and drug development.


title: Common molecular mechanisms underlie the transfer of alpha-synuclein, Tau and huntingtin and modulate spontaneous activity in neuronal cells

abstract:
The misfolding and accumulation of disease-related proteins are common hallmarks among several neurodegenerative diseases. Alpha-synuclein (aSyn), Tau and huntingtin (wild-type and mutant, 25QHtt and 103QHtt, respectively) were recently shown to be transferred from cell-to-cell through different cellular pathways, thereby contributing to disease progression and neurodegeneration. However, the relative contribution of each of these mechanisms towards the spreading of these different proteins and the overall effect on neuronal function is still unclear.

To address this, we exploited different cell-based systems to conduct a systematic comparison of the mechanisms of release of aSyn, Tau and Htt, and evaluated the effects of each protein upon internalization in microglial, astrocytic, and neuronal cells. In the models used, we demonstrate that 25QHtt, aSyn and Tau are released to the extracellular space at higher levels than 103QHtt, and their release can be further augmented with the co-expression of USP19. Furthermore, cortical neurons treated with recombinant monomeric 43QHtt exhibited alterations in neuronal activity that correlated with the toxicity of the polyglutamine expansion. Tau internalization resulted in an increase in neuronal activity, in contrast to slight effects observed with aSyn. Interestingly, all these disease-associated proteins were present at higher levels in ectosomes than in exosomes. The internalization of both types of extracellular vesicles (EVs) by microglial or astrocytic cells elicited the production of pro-inflammatory cytokines and promoted an increase in autophagy markers. Additionally, the uptake of the EVs modulated neuronal activity in cortical neurons.

Overall, our systematic study demonstrates the release of neurodegenerative disease-associated proteins through similar cellular pathways. Furthermore, it emphasizes that protein release, both in a free form or in EVs, might contribute to a variety of detrimental effects in receiving cells and to progression of pathology, suggesting they may be exploited as valid targets for therapeutic intervention in different neurodegenerative diseases.

Graphical abstract

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=153 SRC="FIGDIR/small/452825v1_ufig1.gif" ALT="Figure 1">
View larger version (31K):
org.highwire.dtl.DTLVardef@1400e8org.highwire.dtl.DTLVardef@1478fccorg.highwire.dtl.DTLVardef@1319abcorg.highwire.dtl.DTLVardef@1a5ad40_HPS_FORMAT_FIGEXP  M_FIG C_FIG


title: FMRP regulates mRNAs encoding distinct functions in the cell body and dendrites of CA1 pyramidal neurons

abstract: Neurons are believed to rely on dendritic localization and translation of mRNAs in order to generate activity-dependent changes in the synaptic plasticity. Here, we develop a strategy combining compartment-specific CLIP and TRAP in conditionally tagged mice to precisely define the ribosome-bound dendritic transcriptome of CA1 pyramidal neurons. This revealed transcripts that have differentially localized alternative 3UTR and splicing isoforms. FMRP targets are overrepresented among dendritic mRNAs, and compartment-specific FMRP-CLIP defined 383 dendritic FMRP targets, and also allowed for segregation of whole-cell FMRP targets into functional modules that are locally regulated by FMRP. In the absence of FMRP, dendritic FMRP targets show increased ribosome association, consistent with reported roles for FMRP in translational repression. Together, the data support a model in which distinct patterns of FMRP localization allow it to differentially regulate the expression of nuclear proteins and synaptic proteins within different compartments of a single neuronal cell type.


title: EEG Biomarkers of reduced inhibition in human cortical microcircuits in depression

abstract: Reduced cortical inhibition by somatostatin-expressing (SST) interneurons has been strongly associated with treatment-resistant depression. However, whether the effects of reduced SST interneuron inhibition on microcircuit activity have signatures detectible in electroencephalography (EEG) signals remains unknown. We simulated resting-state activity and EEG using detailed models of human cortical microcircuits with normal (healthy) or reduced SST interneuron inhibition (depression). Healthy microcircuit models showed emergent key features of resting-state EEG, and depression microcircuits exhibited increased theta, alpha and low beta power (4 - 15 Hz). The changes in depression involved a combination of an aperiodic broadband, and periodic theta and low beta components. We then demonstrated the specificity of the EEG signatures of reduced SST interneuron inhibition by showing they were distinct from those corresponding to reduced parvalbumin-expressing (PV) interneuron inhibition. Our study thus links SST interneuron inhibition level to distinct features in EEG simulated from detailed human microcircuits, which can serve to better identify mechanistic subtypes of depression using EEG, and non-invasively monitor modulation of cortical inhibition.


title: EEG biomarkers of reduced inhibition in human cortical microcircuits in depression


title: Altered cortical processing of sensory input in Huntington disease mouse models

abstract: Huntington disease (HD), a hereditary neurodegenerative disorder, manifests as progressively impaired movement and cognition. Although early abnormalities of neuronal activity in striatum are well established in HD models, there are fewer in vivo studies of the cortex. Here, we record local field potentials (LFPs) in YAC128 HD model mice versus wild-type mice. In multiple cortical areas, limb sensory stimulation evokes a greater change in LFP power in YAC128 mice. Mesoscopic imaging using voltage-sensitive dyes reveal more extensive spread of evoked sensory signals across the cortical surface in YAC128 mice. YAC128 layer 2/3 sensory cortical neurons ex vivo show increased excitatory events, which could contribute to enhanced sensory responses in vivo. Cortical LFP responses to limb stimulation, visual and auditory input are also significantly increased in zQ175 HD mice. Results presented here extend knowledge of HD beyond ex vivo studies of individual neurons to the intact cortical network.


title: Senescence of stromal cells contributes to endometrium dysfunction and embryo implantation failure

abstract: Successful implantation requires two-sided interaction between genetically normal embryo and properly prepared endometrium. At the cellular level the latter means hormone-regulated differentiation of endometrial stromal cells (ESCs) into decidual ones that create appropriate microenvironment for invading embryo. Impaired decidualization is proved to mediate implantation failures. Here we elicited ESCs senescence as the cause for disturbed decidualization of endometrial stroma and impaired blastocyst implantation. Ability to decidualize and to accept modeled blastocysts inversely correlated with senescence in patients ESCs lines. Reduced hormonal responsiveness of senescent ESCs led to inappropriate decidualization dynamics resulting in altered receptivity, disturbed ligand-receptor interaction with trophoblasts and modified architecture of extracellular matrix what hindered blastocysts invasion. Furthermore, senescent ESCs caused  bystander quenching of decidual reaction in adjacent cells reinforcing dysfunction of stromal compartment. Implementation of senomorphics reducing senescence phenotype diminished adverse effects of senescent ESCs on decidualization and implantation using both in vitro models and patients lines, what suggests a promising strategy to increase in vitro fertilization efficacy.


title: CLARA: A web portal for interactive exploration of the cardiovascular cellular landscape in health and disease

abstract: Mammalian cardiovascular tissues are comprised of complex and diverse collections of cells. Recent advances in single-cell profiling technologies have accelerated our understanding of tissue cellularity and the molecular networks that orchestrate cardiovascular development, maintain homeostasis, and are disrupted in pathological states. Despite the rapid development and application of these technologies, many cardiac single-cell functional genomics datasets remain inaccessible for most cardiovascular biologists. Access to custom visual representations of the data, including querying changes in cellular phenotypes and interactions in diverse contexts, remains unavailable in publicly accessible data portals. Visualizing data is also challenging for scientists without expertise in processing single-cell genomic data. Here we present CLARA--CardiovascuLAR Atlas--a web portal facilitating exploration of the cardiovascular cellular landscape. Using mouse and human single-cell transcriptomic datasets, CLARA enables scientists unfamiliar with single-cell-omic data analysis approaches to examine gene expression patterns and the cell population dynamics of cardiac cells in a range of contexts. The web-application also enables investigation of intercellular interactions that form the cardiac cellular niche. CLARA is designed for ease-of-use and we anticipate that the portal will aid deeper exploration of cardiovascular cellular landscapes in the context of development, homeostasis and disease. CLARA is freely available at https://clara.baker.edu.au.


title: Facultative variation across a shallow to deep torpor spectrum in hummingbirds

abstract: Many small endotherms use torpor, saving energy by a controlled reduction of their body temperature and metabolic rate. Some species (e.g. arctic ground squirrels, hummingbirds) enter deep torpor, dropping their body temperatures by 23-37 {degrees}C, while others can only enter shallow torpor (e.g., pigeons, 3-10 {degrees}C reductions). However, deep torpor in mammals can increase predation risk (unless animals are in burrows or caves), inhibit immune function, and result in sleep deprivation, so even for species that can enter deep torpor, facultative shallow torpor might help balance energy savings with these potential costs. Deep torpor occurs in three avian orders. Although the literature hints that some bird species can use both shallow and deep torpor, little empirical evidence of such an avian torpor spectrum exists. We infrared imaged three hummingbird species that are known to use deep torpor, under natural temperature and light cycles, to test if they were also capable of shallow torpor. All three species used both deep and shallow torpor, often on the same night. Depending on the species, they used shallow torpor for 5-35% of the night. The presence of a bird torpor spectrum indicates a capacity for fine-scale physiological and genetic regulation of avian torpid metabolism.


title: Evolution of Plant Niche Construction Traits in Biogeomorphic Landscapes

abstract: By virtue of niche construction traits, plants play a significant role in shaping landscapes. The resultant outcome is a change in the selective environment, which influences the evolution of these same plants. So far almost all biogeomorphic models have assumed that niche construction traits are invariant in time. On the other hand, niche construction studies have assumed that independent abiotic changes are either nonexistent or are simply linear. Here, I considered the concomitant evolution of plant niche construction traits during landscape development. I constructed a geo-evolutionary model that couples a population genetic module with a landscape development module. Allowing plants to evolve always results in landforms different from those that appear when evolution is not accounted for. The topographic difference between cases with and without evolution ranges from a small difference in the steady-state topography, to drastic differences in landforms. The amount of difference is contingent upon forms of landscape development and the strength of geo-evolutionary coupling. Allowing the landscape to develop while evolution occurs changes evolutionary trajectories for niche construction traits. The landscape can even develop spatial structures that suppress selection. Model results clearly support the need to integrate niche construction theory and biogeomorphology to better understand both.


title: Empirical prediction of variant-associated cryptic-donors with 87% sensitivity and 95% specificity

abstract: Predicting which cryptic-donors may be activated by a genetic variant is notoriously difficult. Through analysis of 5,145 cryptic-donors activated by 4,811 variants (versus 86,963 decoy-donors not used; any GT or GC), we define an empirical method predicting cryptic-donor activation with 87% sensitivity and 95% specificity. Strength (according to four algorithms) and proximity to the authentic-donor appear important determinants of cryptic-donor activation. However, other factors such as auxiliary splicing elements, which are difficult to identify, play an important role and are likely responsible for current prediction inaccuracies. We find that the most frequent mis-splicing events at each exon-intron junction, mined from 40,233 RNA-sequencing samples, predict with remarkable accuracy which cryptic-donor will be activated in rare disease. Aggregate RNA-Sequencing splice-junction data provides an accurate, evidence-based method to predict variant-activated cryptic-donors in genetic disorders, assisting pathology consideration of possible consequences of a variant for the encoded protein and RNA diagnostic testing strategies.


title: Single cell analysis reveals X chromosome upregulation is not global and primarily belongs to ancestral genes in pre-gastrulation embryos

abstract:
Evolution of sex chromosome dosage compensation in mammals remains poorly understood. Ohnos hypothesis state that evolution of dosage compensation occurred through two steps: first, to compensate the dosage imbalance created due to the degradation of Y chromosome in male, upregulation of X-chromosome happened. Subsequently, transmission of X-chromosome upregulation (XCU) in female led to the evolution of X-chromosome inactivation (XCI) to counteract extra dosage of X-linked genes in female cells. Here, we have profiled gene-wise dynamics of XCU in pre-gastrulation mouse embryos at single cell level and found that XCU is dynamically linked with XCI, however, XCU is not global or chromosome-wide like XCI. Therefore, our result raises question whether XCU driven the evolution of XCI. If so, then why XCI is chromosome wide while XCU is not. We propose that XCI might have evolved independent of XCU and therefore refining the current model is necessary. Separately, we show that higher occupancy of different activating factors at upregulated X-linked genes leads to enhanced transcriptional burst frequency and thereby leads to upregulation. On the other hand, our analysis indicates that extent of upregulation, enrichment of different activating marks differs between ancestral and newly acquired X-linked genes. Altogether, our study provides significant insight into the dynamics and mechanistic basis of evolution of sex chromosome dosage compensation.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY


title: X upregulation is not global and extent of upregulation differs between ancestral and acquired X-linked genes


title: ACVR1 antibodies exacerbate heterotopic ossification in fibrodysplasia ossificans progressiva (FOP) by activating FOP-mutant ACVR1

abstract: Fibrodysplasia ossificans progressiva (FOP) is a rare genetic disorder whose most debilitating pathology is progressive and cumulative heterotopic ossification (HO) of skeletal muscles, ligaments, tendons, and fascia. FOP is caused by mutations in the type I BMP receptor gene ACVR1, which enable ACVR1 to utilize its natural antagonist, Activin A, as an agonistic ligand. The physiological relevance of this property is underscored by the fact that HO in FOP is exquisitely dependent on activation of FOP-mutant ACVR1 by Activin A, an effect countered by inhibition of Activin A via monoclonal antibody treatment. Hence, we surmised that ACVR1 antibodies that block activation of ACVR1 by ligand should also inhibit HO in FOP and provide an additional therapeutic option for this condition. Therefore, we generated ACVR1 monoclonal antibodies that block ACVR1s activation by its ligands. Surprisingly, in vivo, these ACVR1 antibodies stimulate HO and activate signaling of FOP-mutant ACVR1. This property is restricted to FOP-mutant ACVR1 and results from ACVR1 antibody-mediated dimerization of ACVR1. Conversely, wild type ACVR1 is inhibited by ACVR1 antibodies. These results uncover an additional novel property of FOP-mutant ACVR1 and indicate that ACVR1 antibodies should not be considered as therapeutics for FOP.


title: Fast and efficient CRISPR-mediated genome editing in Aureobasidium pullulans using Cas9 ribonucleoproteins

abstract: Aureobasidium pullulans is a ubiquitous, polyextremotolerant, "yeast-like" ascomycete used for the industrial production of pullulan and other products and as biocontrol agent in the agriculture. Its application potential and its wide-spread occurrence make A. pullulans an interesting study object. The availability of a fast and efficient genome editing method is an obvious advantage for future basic and applied research on A. pullulans. In this study, we describe the development of a CRISPR/Cas9-based genome editing method using ribonucleoproteins (RNPs). We demonstrate that this method can be used for single and multiplex genome editing using only RNPs by targeting ura3 (encoding for orotidine-5'-phosphate decarboxylase), praics (encoding for phosphoribosyl aminoimidazole-succinocarboxamide synthase) and asl (encoding for arginine succinate lyase). We demonstrate the applicability of Trichoderma reesei pyr4 and Aspergillus fumigatus pyrG to complement the ura3 deficiency. Further, we show that the usage of RNPs can boost the homologous recombination rate up to nearly 100%, even when using only 20bp long homologous flanks. Therefore, the repair cassettes can be constructed by a single PCR, abolishing the need for laborious and time-consuming cloning. The here presented method allows fast and efficient genome editing for gene deletions, modifications, and insertions in A. pullulans.


title: gmxapi: a Gromacs-native Python interface for molecular dynamics with ensemble and plugin support

abstract:
Gmxapi provides an integrated, native Python API for both standard and advanced molecular dynamics simulations in GROMACS. The Python interface permits multiple levels of integration with the core GROMACS libraries, and legacy support is provided via an interface that mimics the command-line syntax, so that all GROMACS commands are fully available. Gmxapi has been officially supported since the GROMACS 2019 release and is enabled by default in current versions of the software. Here we describe gmxapi 0.3 and later. Beyond simply wrapping GROMACS library operations, the API permits several advanced operations that are not feasible using the prior command-line interface. First, the API allows custom user plugin code within the molecular dynamics force calculations, so users can execute custom algorithms without modifying the GROMACS source. Second, the Python interface allows tasks to be dynamically defined, so high-level algorithms for molecular dynamics simulation and analysis can be coordinated with loop and conditional operations. Gmxapi makes GROMACS more accessible to custom Python scripting while also providing support for high-level data-flow simulation algorithms that were previously feasible only in external packages.

Author SummaryThe gmxapi software provides a Python interface for molecular dynamics simulations in GROMACS. In addition to simply wrapping GROMACS commands, it supports custom user plugin code, ensemble simulation, and data-flow chaining of commands. As such, gmxapi enables the writing and execution of high-level simulation algorithms. The software ships with GROMACS and is freely available under an LGPL2 license.


title: Auxin and pectin remodeling interplay during rootlet emergence in white lupin

abstract: Secondary root emergence is a crucial trait that shapes the plants underground system. Virtually every developmental step of root primordium morphogenesis is controlled by auxin. However, how the hormone controls cell separation in primordium-overlaying tissues through wall loosening is poorly understood. Here, we took advantage of white lupin and its spectacular cluster root development to assess the contribution of auxin to this process. We show that auxins positive role on rootlet emergence is associated with an upregulation of cell wall pectin modifying and degrading genes. Downregulation of a pectinolytic enzyme gene expressed in cells surrounding the primordium resulted in delayed emergence. Pectins were demethylesterified in the emergence zone and auxin treatment further enhanced this effect. Additionally, we report specific rhamnogalacturonan-I modifications during cortical cell separation. In conclusion, we propose a model in which auxin has a dual role during rootlet emergence: Firstly, through active pectin demethylesterification and secondly by regulating the expression of cell wall remodeling enzymes.

category: plant biology


title: A unified framework for measuring selection on cellular lineages and traits

abstract: Intracellular states probed by gene expression profiles and metabolic activities are intrinsically noisy, causing phenotypic variations among cellular lineages. Understanding the adaptive and evolutionary roles of such variations requires clarifying their linkage to population growth rates. Extending a cell lineage statistics framework, here we show that a populations growth rate can be expanded by fitness cumulants of any cell lineage trait. The expansion enables quantifying the contribution of each fitness cumulant, such as variance and skewness, to population growth. We introduce a function that contains all the essential information of cell lineage statistics, including mean lineage fitness and selection strength. We reveal a relation between fitness heterogeneity and population growth rate response to perturbation. We apply the framework to experimental cell lineage data from bacteria to mammalian cells, revealing that third or higher-order cumulants contributions are negligible under constant growth conditions but could be significant in regrowing processes from growth-arrested conditions. Furthermore, we identify cellular populations in which selection leads to an increase of fitness variance among lineages. The framework assumes no particular growth models or environmental conditions, and is thus applicable to various biological phenomena for which phenotypic heterogeneity and cellular proliferation are important.

category: systems biology


title: Ensemble encoding of conditioned fear by prefrontal somatostatin interneurons

abstract: Neurons preferentially activated by learning have been ascribed the unique potential to encode memory. However, it remains unclear which genetically-defined cell types are recruited as part of such an ensemble, or what role discrete subpopulations play in behavior. Here we show that fear conditioning activates a heterogeneous neural ensemble in the medial prefrontal cortex (mPFC), comprised to a large degree of GABAergic interneurons immunoreactive for somatostatin (SST-INs). Using an intersectional genetic approach, we demonstrate that fear learning-activated SST-INs exhibit distinct circuit properties, are preferentially reactivated during memory retrieval, and mediate the expression of defensive freezing. We further show that a rewarding experience, morphine treatment, activates an orthogonal SST-IN population that exerts opposing control over fear. These results outline an important role for discrete GABAergic ensembles in fear memory encoding, and point to an unappreciated capacity for functional specialization among SST-INs.


title: The critical role of spiny stellate cells in seizure onset based on dynamic analysis of a neural mass model

abstract: Growing evidence suggests that excitatory neurons in the brain play a significant role in seizure generation. Nonetheless, spiny stellate cells are cortical excitatory non-pyramidal neurons in the brain, which their basic role in seizure occurrence is not well understood. In the present research, we study the critical role of spiny stellate cells or the excitatory interneurons (EI), for the first time, in epileptic seizure generation using an extended neural mass model inspired by Liu and Wang model in 2017. Applying bifurcation analysis on this modified model, we investigated the rich dynamics corresponding to the epileptic seizure onset and transition between interictal and ictal states caused by EI connectivity to other cell types. Our results indicate that the transition between interictal and ictal states (preictal signal) corresponds to a supercritical Hopf bifurcation, and thus, the extended model suggests that before seizure onset, the amplitude and frequency of neural activities gradually increase. Moreover, we showed that 1) the altered function of GABAergic and glutamatergic receptors of EI can cause seizure, and 2) the pathway between the thalamic relay nucleus and EI facilitates the transition from interictal to the ictal activity by decreasing the preictal period. Thereafter, we considered both sensory and cortical periodic inputs to study model responses to various harmonic stimulations. Bifurcation analysis of the model, in this case, suggests that the initial state of the model might be the main cause for the transition between interictal and ictal states as the stimulus frequency changes. The extended thalamocortical model shows also that the amplitude jump phenomenon and nonlinear resonance behavior result from the preictal state of the modified model. These results can be considered as a step forward to a deeper understanding of the mechanisms underlying the transition from normal activities to epileptic activities.


title: The critical modulatory role of spiny stellate cells in seizure onset based on dynamic analysis of a neural mass model


title: Macrophages govern antiviral responses in human lung tissues protected from SARS-CoV-2 infection

abstract:
The majority of SARS-CoV-2 infections among healthy individuals result in asymptomatic to mild disease. However, the immunological mechanisms defining effective lung tissue protection from SARS-CoV-2 infection remain elusive. Unlike mice solely engrafted with human fetal lung xenograft (fLX), mice co-engrafted with fLX and a myeloid-enhanced human immune system (HNFL mice) are protected against SARS-CoV-2 infection, severe inflammation, and histopathology. Effective control of viral infection in HNFL mice associated with significant macrophage infiltration, and the induction of a potent macrophage-mediated interferon response. The pronounced upregulation of the USP18-ISG15 axis (a negative regulator of IFN responses), by macrophages was unique to HNFL mice and represented a prominent correlate of reduced inflammation and histopathology. Altogether, our work shed light on unique cellular and molecular correlates of lung tissue protection during SARS-CoV-2 infection, and underscores macrophage IFN responses as prime targets for developing immunotherapies against coronavirus respiratory diseases.

HIGHLIGHTSO_LIMice engrafted with human fetal lung xenografts (fLX-mice) are highly susceptible to SARS-CoV-2.
C_LIO_LICo-engraftment with a human myeloid-enriched immune system protected fLX-mice against infection.
C_LIO_LITissue protection was defined by a potent and well-balanced antiviral response mediated by infiltrating macrophages.
C_LIO_LIProtective IFN response was dominated by the upregulation of the USP18-ISG15 axis.
C_LI

category: microbiology


title: Vaccinia Virus Arrests and Shifts the Cell Cycle

abstract: Modulation of the host cell cycle is a common strategy used by viruses to create a pro-replicative environment. To facilitate viral genome replication, vaccinia virus (VACV) has been reported to alter cell cycle regulation and trigger the host cell DNA damage response. However, the cellular factors and viral effectors that mediate these changes remain unknown. Here, we set out to investigate the effect of VACV infection on cell proliferation and host cell cycle progression. Using a subset of VACV mutants we characterize the stage of infection required for inhibition of cell proliferation and define the viral effectors required to dysregulate the host cell cycle. Consistent with previous studies, we show that VACV inhibits, and subsequently shifts the host cell cycle. We demonstrate that these two phenomena are independent of one another, with viral early genes being responsible for cell cycle inhibition, and post-replicative viral gene(s) responsible for the cell cycle shift. Extending previous findings, we show that the viral kinase F10 is required to activate the DNA damage checkpoint and that the viral B1/B12 (pseudo) kinases mediate degradation of checkpoint effectors p53 and p21 during infection. We conclude that VACV modulates host cell proliferation and host cell cycle progression through temporal expression of multiple VACV effector proteins.


title: Rhodanese Rdl2 produces reactive sulfur species to scavenge hydroxyl radical and protect mitochondria

abstract: During aerobic respiration, mitochondria generate superoxide anion (O2*-), hydrogen peroxide (H2O2), and hydroxyl radical (HO*), and these reactive oxygen species (ROS) are detrimental to mitochondria. Mitochondrial damage is linked to a broad spectrum of pathologies such as Alzheimers disease, hemochromatosis, and diabetes. Mitochondria contain several enzymes for rapidly removing superoxide anion and hydrogen peroxide, but how they antagonize HO* is elusive, representing a loophole in the anti-ROS system. Herein, we discovered that Rhodanese 2 (Rdl2) is critical for maintaining the functionality and integrity of mitochondria under sub-lethal ROS stress in Saccharomyces cerevisiae. Rdl2 converts stable sulfur species (thiosulfate and dialkyl polysulfide) to reactive sulfane sulfur including persulfide that protects mitochondrial DNA via scavenging HO*. Surprisingly, hydrogen sulfide (H2S) promotes HO* production through stimulating the Fenton reaction, leading to increased DNA damage. Our study may reveal an ex-ante mean for antagonizing HO*, patching the loophole of the anti-ROS system in mitochondria.


title: Parkinson' disease medication alters rat small intestinal motility and microbiota composition

abstract:
Parkinsons disease (PD) is known to be associated with altered gastrointestinal function and microbiota composition. Altered gastrointestinal function is key in the development of small intestinal bacterial overgrowth, which is a comorbidity often observed in PD patients. Although PD medication could be an important confounder in the reported alterations, its effect per se on the microbiota composition or the gastrointestinal function at the site of drug absorption has not been studied. To this end, healthy (i.e., not PD model) wild-type Groningen rats were employed and treated with dopamine, pramipexole (in combination with levodopa/carbidopa), or ropinirole (in combination with levodopa/carbidopa) for 14 sequential days. Rats treated with dopamine agonists showed a significant reduction in the small intestinal motility and an increase in bacterial overgrowth in the distal small intestine. Importantly, significant alterations in microbial taxa were observed between the treated and vehicle groups, analogous to the changes previously reported in human PD vs HC microbiota studies. These microbial changes included an increase in Lactobacillus and Bifidobacterium, and decrease in Lachnospiraceae and Prevotellaceae. Importantly, certain Lactobacillus species correlated negatively with the plasma levels of levodopa. Overall, the results highlight the significant effect of PD medication per se on the gut microbiota, and the disease-associated comorbidities, including gastrointestinal dysfunction and small intestinal bacterial overgrowth.

Research in contextO_ST_ABSEvidence before this studyC_ST_ABSIn the last years many studies focused on the faecal microbiota profiles of Parkinsons disease (PD) patients and compared them to healthy control (HC) subjects and reported a difference between the microbiota profiles. Although some bacterial taxa have been reported to be differential abundant across studies there is a low consensus between the overall changes observed. It remains to be elucidated whether and how the disease status itself, the PD medication or the gastrointestinal dysfunction, an often reported comorbidity, plays a role in the altered microbiota profiles of Parkinsons disease patients.

Added value of this studyTwo major confounding factors potentially affecting the microbiota profiles differentiating PD patients from HC subjects are the PD medication and gastrointestinal dysfunction. We showed that PD medication elicits an decreasing effect on the small intestinal gut motility in a healthy (i.e. non Parkinson) rat model which was a contributing factor to the altered microbiota profile observed. We found that taxa reported previously to be differentially abundant were mirrored in our healthy rat model showing the impact of PD medication on microbiota profiles.

Implications of all the available evidenceThis study showed that PD medication and gastrointestinal motility are important factors in the microbiota profiles and could explain some of the differential abundant taxa reported in the cross-sectional PD microbiota studies. Future studies should take the potential side effects of medication that could elicit alterations of the microbiota composition into account when seeking for microbial biomarkers of disease status.


title: Identification and characterization of hPSC-derived FOXA2+ progenitor cells with ventricular cardiac differentiation potential

abstract: While much progress has been made in understanding early cardiac development, the precise mechanisms that specify the different cardiomyocyte subtypes remain poorly understood. Recent data from our lab have shown that transient Foxa2 expression identifies a progenitor population with exclusive ventricular differentiation potential in the mouse heart. Here we have translated this concept to the human pluripotent stem cell (hPSC) system. Using a FOXA2-GFP reporter cell line we characterized expression of FOXA2 during hPSC cardiac differentiation and found that a subset of cardiac mesoderm precursors transiently expresses FOXA2. Gene expression analysis of FOXA2+ and FOXA2- cardiac mesoderm revealed that both populations similarly express early cardiac specification markers such as PDGFRA, TBX5, and ISL1, while other key candidates including TBX20 and GATA4 are significantly upregulated in the FOXA2+ population. Isolation and subsequent differentiation of FOXA2+ and FOXA2- populations demonstrates their comparable differentiation potential to both cardiomyocytes and epicardial cells. However, cardiomyocytes derived from FOXA2+ precursors showed enhanced differentiation efficiency toward ventricular cardiomyocytes compared to cardiomyocytes derived from FOXA2- precursors. To identify new mechanisms that regulate ventricular specification, we performed small molecule screening and found that inhibition of the EGFR pathway strongly increased the cardiac mesoderm population in general, and the FOXA2+ precursors in particular. Finally, we have identified a combination of cell surface markers to specifically isolate FOXA2+ cardiac precursors. In summary, our results suggest that FOXA2+ cardiac mesoderm harbors ventricular-specific differentiation potential and isolation of these cells permits the generation of cultures enriched for ventricular cardiomyocytes. Generating such enriched cardiac populations will be relevant for regenerative medicine approaches, as well as for disease modeling from induced pluripotent stem cells.


title: Severe introduced predator impacts despite attempted functional eradication

abstract: Established non-native species can have significant impacts on native biodiversity without any possibility of complete eradication. In such cases, one management approach is functional eradication, the reduction of introduced species density below levels that cause unacceptable effects on the native community. Functional eradication may be particularly effective for species with limited dispersal ability, which may limit rates of reinvasion from distant populations. Here, we evaluate the potential for functional eradication of introduced predatory oyster drills (Urosalpinx cinerea) using a community science approach in San Francisco Bay. We combined observational surveys, targeted removals, and a caging experiment to evaluate the effectiveness of this approach in mitigating the mortality of prey Olympia oysters (Ostrea lurida), a conservation and restoration priority species. Despite the efforts of over 300 volunteers that removed over 30,000 oyster drills, we report limited success and discuss several possible mechanisms for this result with broad relevance to management for this and other introduced species. We also found a strong negative relationship between oyster drills and oysters, showing virtually no coexistence across eight sites. At two removal sites, there was no effect of oyster drill removal on oyster survival, which was only observed by caging treatment (0 and 1.6% survival in open and partial cage treatments, as compared to 89.1% in predator exclusion treatments). We conclude that functional eradication of this species requires significantly greater effort and may not be a viable management strategy. Oyster restoration efforts should not be undertaken where Urosalpinx is established or is likely to invade.


title: ACE2 binding is an ancestral and evolvable trait of sarbecoviruses

abstract: Two different sarbecoviruses have caused major human outbreaks in the last two decades1,2. Both these sarbecoviruses, SARS-CoV-1 and SARS-CoV-2, engage ACE2 via the spike receptor-binding domain (RBD)2-6. However, binding to ACE2 orthologs from humans, bats, and other species has been observed only sporadically among the broader diversity of bat sarbecoviruses7-11. Here, we use high-throughput assays12 to trace the evolutionary history of ACE2 binding across a diverse range of sarbecoviruses and ACE2 orthologs. We find that ACE2 binding is an ancestral trait of sarbecovirus RBDs that has subsequently been lost in some clades. Furthermore, we demonstrate for the first time that bat sarbecoviruses from outside Asia can bind ACE2. In addition, ACE2 binding is highly evolvable: for many sarbecovirus RBDs there are single amino-acid mutations that enable binding to new ACE2 orthologs. However, the effects of individual mutations can differ markedly between viruses, as illustrated by the N501Y mutation which enhances human ACE2 binding affinity within several SARS-CoV-2 variants of concern12 but severely dampens it for SARS-CoV-1. Our results point to the deep ancestral origin and evolutionary plasticity of ACE2 binding, broadening consideration of the range of sarbecoviruses with spillover potential.


title: Molecular evolution and structural analyses of the spike glycoprotein from Brazilian SARS-CoV-2 genomes: the impact of the fixation of selected mutations

abstract: The COVID-19 pandemic caused by Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has reached by July 2021 almost 200 million cases and more than 4 million deaths worldwide since its beginning in late 2019, leading to enhanced concern in the scientific community and the general population. One of the most important pieces of this host-pathogen interaction is the spike protein, which binds to the human Angiotensin-converting enzyme 2 (hACE2) cell receptor, mediates the membrane fusion and is the major target of neutralizing antibodies against SARS-CoV-2. The multiple amino acid substitutions observed in this region, specially in the Receptor Binding Domain (RBD), mainly after almost one year of its emergence (late 2020), have enhanced the hACE2 binding affinity and led to several modifications in the mechanisms of SARS-CoV-2 pathogenesis, improving the viral fitness and/or promoting immune evasion, with potential impact in the vaccine development. In this way, the present work aimed to evaluate the effect of positively selected mutations fixed in the Brazilian SARS-CoV-2 lineages and to check for mutational evidence of coevolution. Additionally, we evaluated the impact of selected mutations identified in some of the VOC and VOI lineages (C.37, B.1.1.7, P.1, and P.2) of Brazilian samples on the structural stability of the spike protein, as well as their possible association with more aggressive infection profiles by estimating the binding affinity in the RBD-hACE2 complex. We identified 48 sites under selective pressure in Brazilian spike sequences, 17 of them with the strongest evidence by the HyPhy tests, including VOC related mutation sites 138, 142, 222, 262, 484, 681, and 845, among others. The coevolutionary analysis identified a number of 28 coevolving sites that were found not to be conditionally independent, such as the couple E484K - N501Y from P.1 and B.1.351 lineages. Finally, the molecular dynamics and free energy estimates showed the structural stabilizing effect and the higher impact of E484K for the improvement of the binding affinity between the spike RBD and the hACE2 in P.1 and P.2 lineages, as well as the stabilizing and destabilizing effects for the positively selected sites.


title: Huntington's disease age at motor onset is modified by the tandem hexamer repeat in TCERG1

abstract:
BackgroundHuntingtons disease is caused by an expanded CAG tract in HTT. The length of the CAG tract accounts for over half the variance in age at onset of disease, and is influenced by other genetic factors, mostly implicating the DNA maintenance machinery. We examined a single nucleotide variant, rs79727797, on chromosome 5 in the TCERG1 gene, previously reported to be associated with Huntingtons disease and a quasi-tandem repeat (QTR) hexamer in exon 4 of TCERG1 with a central pure repeat.

MethodsWe developed a novel method for calling perfect and imperfect repeats from exome sequencing data, and tested association between the QTR in TCERG1 and residual age at motor onset (after correcting for the effects of CAG length in the HTT gene) in 610 individuals with Huntingtons disease via regression analysis.

ResultsWe found a significant association between age at onset and the sum of the repeat lengths from both alleles of the QTR (p = 2.1x10-9), with each added repeat hexamer reducing age at onset by one year (95% confidence interval [0.7, 1.4]). This association explained that previously observed with rs79727797.

ConclusionsThe association with age at onset in the genome-wide association study is due to a QTR hexamer in TCERG1, translated to a glutamine/alanine tract in the protein. We could not distinguish whether this was due to cis-effects of the hexamer repeat on gene expression or of the encoded glutamine/alanine tract in the protein. These results motivate further study of the mechanisms by which TCERG1 modifies onset of HD.


title: The genomic formation of Tanka people, an isolated 'Gypsies in water' in the coastal region of Southeast China

abstract:
ObjectivesThree different hypotheses proposed via the controversial evidence from cultural, anthropological and uniparental genetic analysis respectively stated that Tanka people probably originated from Han Chinese, ancient Baiyue tribe, or the admixture of them. Therefore, the genetic origin and admixture history of the Tanka people, an isolated "Gypsies in water" in the coastal region of Southeast China, are needed to be genetically clarified.

Materials and methodsTo elucidate the genetic origin of the Southeast Tanka people and explore their genetic relationship with surrounding indigenous Tai-Kadai and Austronesian people and Neolithic-to-historic ancients from the Yellow River Basin (YRB) and Fujian, we conducted a large-scale population genomic study among 1498 modern and ancient Eurasians, in which 73 Tanka and 4 Han people were first reported here. Both allele-shared and haplotype-based statistical methods were used here, including PCA, ADMIXTURE, f-statistics, ALDER, qpGraph-/TreeMix and qpAdm/qpWave, ChromoPainter and FineSTRUCTURE

ResultsWe found a specific genetic cline in PCA plots and detected the Tanka-specific homogeneous ancestry in model-based ADMIXTURE results, suggesting differentiated demographic history between Tanka and surrounding Hans. Formal tests based on sharing allele patterns showed a close relationship between Tanka people and Han Chinese, but the Tanka population harbored more southern East Asian ancestry related to Austronesian and Tai-Kadai people compared with southern Hans. Besides, the reconstructed differentiated demographic history revealed that southern Xinshizhou Tanka harbored more ancestry related to the Tai-Kadai or coastal ancient Neolithic to Bronze Age East Asians compared with northern Shacheng Tanka. The qpGraph-/TreeMix-based phylogenetic framework, qpAdm/qpWave-based admixture modeling and FineSTRUCTURE-based dendrogram among ancient northern and southern East Asians further demonstrated that the primary ancestry of modern Tanka derived from ancient millet farmers in the YRB with additional admixture from multiple southern East Asian sources.

DiscussionSharing ancestry estimated from the f-statistics and sharing haplotypic landscape inferred from the ChromoPainter and FineSTRUCTURE showed that Southeast Tanka people not only had a close genetic relationship with both Northern Hans and YRB millet farmers but also possessed more southern East Asian ancestry related to Austronesian, Tai-Kadai and Hmong-Mien people. Our genomic data and fitted admixture models supported modern Tanka originated from ancient North China and obtained additional gene flow from ancient southern East Asians in the processes of southward migrations.