title: No behavioral evidence for rhythmic facilitation of perceptual discrimination

abstract: It has been hypothesized that internal oscillations can synchronize (i.e., entrain) to external environmental rhythms, thereby facilitating perception and behavior. To date, evidence for the link between the phase of neural oscillations and behavior has been scarce and contradictory; moreover, it remains an open question whether the brain can use this tentative mechanism for active temporal prediction. In our present study, we conducted a series of auditory pitch discrimination tasks with 181 healthy participants in an effort to shed light on the proposed behavioral benefits of rhythmic cueing and entrainment. In the three versions of our task, we observed no perceptual benefit of purported entrainment: targets occurring in-phase with a rhythmic cue provided no perceptual benefits in terms of discrimination accuracy or reaction time when compared with targets occurring out-of-phase or targets occurring randomly, nor did we find performance differences for targets preceded by rhythmic vs. random cues. However, we found a surprising effect of cueing frequency on reaction time, in which participants showed faster responses to cue rhythms presented at higher frequencies. We therefore provide no evidence of entrainment, but instead a tentative effect of covert active sensing in which a faster external rhythm leads to a faster communication rate between motor and sensory cortices, allowing for sensory inputs to be sampled earlier in time.


title: Regional structural-functional connectome coupling is heritable and associated with age, sex and cognitive scores in adults

abstract: Large scale white matter brain connections quantified via the structural connectome (SC) act as the backbone for the flow of functional activation, which can be represented via the functional connectome (FC). Many studies have used statistical analysis or computational modeling techniques to relate SC and FC at a global, whole-brain level. However, relatively few studies have investigated the relationship between individual cortical and subcortical regions structural and functional connectivity profiles, here called SC-FC coupling, or how this SC-FC coupling may be heritable or related to age, sex and cognitive abilities. Here, we quantify regional SC-FC coupling in a large group of healthy young adults (22 to 37 years) using diffusion-weighted MRI and resting-state functional MRI data from the Human Connectome Project. We find that while regional SC-FC coupling strengths vary widely across cortical, subcortical and cerebellar regions, they were strongest in highly structurally connected visual and subcortical areas. Additionally, depending on the region, SC-FC coupling varied across sexes and with age and composite cognitive scores. Specifically, SC-FC coupling in the cerebellum tended to decrease with age while coupling in medial fronto-orbital areas tended to increase with age. Males had stronger coupling in many regions, particularly in the right orbito-frontal region and areas in the ventral attention and default mode networks, while females had higher coupling strength in right hippocampus. Furthermore, increased SC-FC coupling in the right insula and decreased coupling in bilateral middle cingulate and supplementary motor areas was associated with higher composite cognitive scores. Finally, we found SC-FC coupling to be highly heritable, particularly in higher order default mode, dorsal/ventral attention and fronto-parietal networks. Taken together, these results suggest regional structure-function coupling in young adults varies with age, is generally stronger in males, is associated with composite cognitive scores and is highly heritable.

category: neuroscience


title: Myogenin controls via AKAP6 non-centrosomal microtubule organizing center formation at the nuclear envelope

abstract: Non-centrosomal microtubule organizing centers (MTOC) are pivotal for the function of multiple cell types, but the processes initiating their formation are unknown. Here, we find that the transcription factor myogenin is required in myoblasts for the localization of MTOC proteins to the nuclear envelope. Moreover, myogenin is sufficient in fibroblasts for nuclear envelope MTOC (NE-MTOC) formation and centrosome attenuation. Bioinformatics combined with loss- and gain-of-function experiments identified induction of AKAP6 expression as one central mechanism for myogenin-mediated NE-MTOC formation. Promoter studies indicate that myogenin preferentially induces the transcription of muscle- and NE-MTOC-specific isoforms of Akap6 and Syne1, which encodes nesprin-1, the NE-MTOC anchor protein in muscle cells. Overexpression of AKAP6{beta} and nesprin-1 was sufficient to recruit endogenous MTOC proteins to the nuclear envelope of myoblasts in the absence of myogenin. Taken together, our results illuminate how mammals transcriptionally control the switch from a centrosomal MTOC to an NE-MTOC and identify AKAP6 as a novel NE-MTOC component in muscle cells.

category: cell biology


title: Non-cell autonomous inhibition of the Shh pathway due to impaired cholesterol biosynthesis requires Ptch1/2

abstract: Congenital birth defects due to mutations of enzymes involved cholesterol synthesis, like Smith-Lemli-Opitz syndrome (SLOS) and Lathosterolosis are characterized by an accumulation of cholesterol precursors. The phenotype of both SLOS and Lathosterolosis have similarities to syndromes associated with abnormal Sonic hedgehog (Shh) signaling, consistent with the observation that cholesterol precursors and their derivative can inhibit Shh signaling. Two types of multipass membrane proteins play central roles in Shh signal transduction, the putative Resistance, Nodulation and Division (RND) antiporters Patched (Ptch)1 and -2, and the G-protein coupled receptor Smoothened (Smo). Sterols have been suggested as cargo for Ptch1/2, while Smo activity can affected both positively and negatively by steroidal molecules. We demonstrate that embryonic stem cells with mutations in the 7-dehydroxycholesterol reductase (7Dhcr) or sterol-C5-desaturase (Sc5d) gene reduce the Hh response in adjacent wildtype cells when grown in mosaic organoids. This non-cell autonomous inhibitory activity of the mutant cells requires the presence of both Ptch1 and Ptch2. These observations support a model in which late cholesterol precursors that accumulate in cells lacking 7Dhcr are the cargo for Ptch1 and Ptch2 efflux activity and mediate the non-cell autonomous inhibition of Smo.

category: developmental biology


title: Improved estimation of phenotypic correlations using summary association statistics

abstract: Estimating the phenotypic correlations between complex traits and diseases based on their genome-wide association summary statistics has been a useful technique in genetic epidemiology and statistical genetics inference. Two state-of-the-art strategies, Z-score correlation across null-effect SNPs and LD score regression intercept, were widely applied to estimate phenotypic correlations. Here, we propose an improved Z-score correlation strategy based on SNPs with low minor allele frequencies (MAFs), and show how this simple strategy can correct the bias generated by the current methods. Comparing to LDSC, the low-MAF estimator improves phenotypic correlation estimation thus is beneficial for methods and applications using phenotypic correlations inferred from summary association statistics.

category: genetics


title: Activity-based urinary biomarkers of response and resistance to checkpoint blockade immunotherapy

abstract:
Immune checkpoint blockade (ICB) therapy has transformed cancer treatment, yet most patients do not derive clinical benefit and responders can acquire resistance to therapy. Noninvasive biomarkers are needed to indicate early on-treatment response and resistance mechanisms. Here we developed ImmuNe Sensors for monItorinG cHeckpoint blockade Therapy (INSIGHT), which comprises a library of mass-barcoded peptide substrates conjugated to PD1 antibodies, as therapeutic sensors. Following systemic administration, INSIGHT carries out the dual role of reinvigorating T cell function and profiling T cell and tumor proteases by the release of cleaved peptides into urine for noninvasive detection by mass spectrometry. We show that an PD1 therapeutic sensor for Granzyme B discriminates early treatment responses before tumor volumes significantly diverge from isotype controls in murine models of colorectal cancer. To differentiate mechanisms of resistance by multivariate analysis, we design a mass-barcoded, 14-plex INSIGHT library to profile proteases differentially expressed by tumors harboring B2m or Jak1 loss-of-function mutations. We find that binary classifiers trained on urine samples indicate response to PD-1 therapy as early as the start of the second dose, and discriminate B2m from Jak1 resistance with high sensitivity and specificity (AUROCs > 0.9). Our data supports the use of activity-based biomarkers for early on-treatment response assessment and classification of refractory tumors based on resistance mechanisms.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=109 SRC="FIGDIR/small/420265v2_ufig1.gif" ALT="Figure 1">
View larger version (26K):
org.highwire.dtl.DTLVardef@12ce799org.highwire.dtl.DTLVardef@118fcc5org.highwire.dtl.DTLVardef@17b227eorg.highwire.dtl.DTLVardef@115a7fc_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: bioengineering


title: Efficacy and accuracy responses of DNA mini-barcodes in species identification under a supervised machine learning approach

abstract: Specific gene regions in DNA, such as cytochrome c oxidase I (COI) in animals, are defined as DNA barcodes and can be used as identifiers to distinguish species. The standard length of a DNA barcode is approximately 650 base pairs (bp). However, because of the challenges associated with sequencing technologies and the unavailability of high-quality genomic DNA, it is not always possible to obtain the full-length barcode sequence of an organism. Recent studies suggest that mini-barcodes, which are shorter (100-300 bp) barcode sequences, can contribute significantly to species identification. Among various methods proposed for the identification task, supervised machine learning methods are effective. However, any prior work indicating the efficacy of mini-barcodes in species identification under a machine learning approach is elusive to find. In this study, we analyzed the effect of different barcode lengths on species identification using supervised machine learning and proposed a general approximation of the required length of the minibarcode. Since Naive Bayes is seen to generally outperform other supervised methods in species identification in other studies, we implemented this classifier and showed the effectiveness of the mini-barcode by demonstrating the accuracy responses obtained after varying the length of the DNA barcode sequences.

category: bioinformatics


title: A pragmatic reevaluation of the efficacy of nonhuman primate optogenetics for psychiatry

abstract: Translational neuroscience is committed to generating discoveries in the laboratory that ultimately can improve human lives. Optogenetics has received considerable attention because of its demonstrated promise in rodent brains to manipulate cells and circuits. In a recent report, Tremblay and colleagues (2020) introduce an open resource detailing optogenetic studies of the nonhuman primate (NHP) brain and make robust claims about the translatability of the technology. We propose that their quantitative (e.g., a 91% success rate) and theoretical claims are questionable because the data were analyzed at a level relevant to the rodent but not NHP brain, injections were clustered within a few monkeys in a few studies in a few brain regions, and their definitions of success was not clearly relevant to human neuropsychiatric disease. A reanalysis of the data with a modified definition of success that included a behavioral and biological effect revealed an 62.5% success rate that was lower when considering only strong outcomes (53.1%). This calls into question the current efficacy of optogenetic techniques in the NHP brain and suggests that we are a long way from being able to leverage them in "the service of patients with neurological or psychiatric conditions" as the Tremblay report claims.


title: A feature-based network analysis and fMRI meta-analysis reveal three distinct types of prosocial decisions

abstract: Tasks that measure correlates of prosocial decision-making share one common feature: agents can make choices that increase the welfare of a beneficiary. However, prosocial decisions vary widely as a function of other task features. The diverse ways that prosociality is defined and the heterogeneity of prosocial decisions have created challenges for interpreting findings across studies and identifying their neural correlates. To overcome these challenges, we aimed to organize the prosocial decision-making task-space of neuroimaging studies. We conducted a systematic search for studies in which participants made decisions to increase the welfare of others during fMRI. We identified shared and distinct features of these tasks and employed an unsupervised graph-based approach to assess how various forms of prosocial decision-making are related in terms of their low-level components (e.g., task features like potential cost to the agent or potential for reciprocity). Analyses uncovered three clusters of prosocial decisions, which we labeled cooperation, equity, and altruism. This feature-based representation of the task structure was supported by results of a neuroimaging meta-analysis that each type of prosocial decisions recruited diverging neural systems. Results clarify some of the existing heterogeneity in how prosociality is conceptualized and generate insight for future research and task paradigm development.


title: Cross-species Transcriptomic Network Analysis Reveals Links Between Ribosomal Protein Mutation and Cancer

abstract: Ribosomal proteins (RPs) are critical to all cellular operations through their key roles in ribosome biogenesis and translation, as well as their extra-ribosomal functions. Leveraging previously identified RP mutants we characterised the RP macro transcriptome and sought to compare it with transcriptomes of pathologies associated with failures of ribosomal function, cancer and Zika virus (ZIKV). Single gene-based analysis revealed highly variable transcriptomes of RP mutations with little overlap in genes that were differentially expressed. In contrast, weighted gene co-expression network analysis revealed a highly conserved transcriptomic network pattern across all RP mutants studied. In addition, when we compared network changes in RP mutants, we observed similarities to transcriptome alterations in human cancer, and thus confirming the oncogenic role of RPs. Finally, it is known that ZIKV infection influences translational machinery, but this study shows infection network changes dissimilar to those of either the RP mutation or cancer.


title: Synchrony of Non-isochronous Signals in an Acoustically Communicating Katydid

abstract: The ability to entrain to auditory stimuli has been a powerful method to investigate the comparative rhythm abilities of different animals. While synchrony to regular (isochronous) rhythms is well documented, synchrony to non-isochronous stimuli, with multiple components at unequal time intervals, is rarer. Several katydid species with isochronous calls have been shown to achieve synchrony as part of their natural calling interactions in multi-individual choruses. Yet no study so far has demonstrated synchrony in any insect with a non-isochronous call. Using natural calling behaviour and playback experiments, we investigated acoustic synchrony and the mechanisms underlying it in the katydid species Mecopoda  Two Part Caller. This species has a complex non-isochronous call consisting of a long trill followed by two or more chirps. We found that individual males synchronized trills and, to a lesser extent, chirps. Further investigation of trill synchrony showed that timing of trills is modified by external trills but not chirps. Chirp synchrony is modified by external chirps but also by trills. We suggest a two-oscillator model underlying synchrony in this species and discuss the implications for the evolution of acoustic synchrony.

category: animal behavior and cognition


title: Concurrent assessment of gait kinematics using marker-based and markerless motion capture

abstract: Kinematic analysis is a useful and widespread tool used in research and clinical biomechanics for the estimation of human pose and the quantification of human movement. Common marker-based optical motion capture systems are expensive, time intensive, and require highly trained operators to obtain kinematic data. Markerless motion capture systems offer an alternative method for the measurement of kinematic data with several practical benefits. This work compared the kinematics of human gait measured using a deep learning algorithm-based markerless motion capture system to those of a common marker-based motion capture system. Thirty healthy adult participants walked on a treadmill while data were simultaneously recorded using eight video cameras (markerless) and seven infrared optical motion capture cameras (marker-based). Video data were processed using markerless motion capture software, marker-based data were processed using marker-based capture software, and both sets of data were compared. The average root mean square distance (RMSD) between corresponding joints was less than 2.5 cm for all joints except the hip, which was 3.6 cm. Lower limb segment angles indicated pose estimates from both systems were very similar, with RMSD of less than 5.5{degrees} for all segment angles except those that represent rotations about the long axis of the segment. Lower limb joint angles captured similar patterns for flexion/extension at all joints, ab/adduction at the knee and hip, and toe-in/toe-out at the ankle. These findings demonstrate markerless motion capture can measure similar 3D kinematics to those from marker-based systems.


title: Assembling Long Accurate Reads Using Multiplex de Bruijn Graphs

abstract: Although most existing genome assemblers are based on the de Bruijn graphs, it remains unclear how to construct these graphs for large genomes and large k-mer sizes. This algorithmic challenge has become particularly important with the emergence of long high-fidelity (HiFi) reads that were recently utilized to generate a semi-manual telomere-to-telomere assembly of the human genome and to get a glimpse into biomedically important regions that evaded all previous attempts to sequence them. To enable automated assemblies of long and accurate reads, we developed a fast LJA algorithm that reduces the error rate in these reads by three orders of magnitude (making them nearly error-free) and constructs the de Bruijn graph for large genomes and large k-mer sizes. Since the de Bruijn graph constructed for a fixed k-mer size is typically either too tangled or too fragmented, LJA uses a new concept of a multiplex de Bruijn graph with varying k-mer sizes. We demonstrate that LJA improves on the state-of-the-art assemblers with respect to both accuracy and contiguity and enables automated telomere-to-telomere assemblies of entire human chromosomes.


title: Functional crosstalk between phosphorylation and disease-causing mutations in the cardiac sodium channel Nav1.5

abstract:
The voltage-gated sodium channel Nav1.5 initiates the cardiac action potential. Alterations of its activation and inactivation properties due to mutations can cause severe, life-threatening arrhythmias. Yet despite intensive research efforts, many functional aspects of this cardiac channel remain poorly understood. For instance, Nav1.5 undergoes extensive post-translational modification in vivo, but the functional significance of these modifications is largely unexplored, especially under pathological conditions. This is because most conventional approaches are unable to insert metabolically stable post-translational modification mimics, thus preventing a precise elucidation of the contribution by these modifications to channel function. Here, we overcome this limitation by using protein semi-synthesis of Nav1.5 in live cells and carry out complementary molecular dynamics simulations. We introduce metabolically stable phosphorylation mimics on both WT and two pathogenic long-QT mutant channel backgrounds and decipher functional and pharmacological effects with unique precision. We elucidate the mechanism by which phosphorylation of Y1495 impairs steady-state inactivation in WT Nav1.5. Surprisingly, we find that while the Q1476R patient mutation does not affect inactivation on its own, it enhances the impairment of steady-state inactivation caused by phosphorylation of Y1495 through enhanced unbinding of the inactivation particle. We also show that both phosphorylation and patient mutations can impact Nav1.5 sensitivity towards the clinically used anti-arrhythmic drugs quinidine and ranolazine, but not flecainide. The data highlight that functional effects of Nav1.5 phosphorylation can be dramatically amplified by patient mutations. Our work is thus likely to have implications for the interpretation of mutational phenotypes and the design of future drug regimens.

Significance statementThe cardiac sodium channel (Nav1.5) is crucial for generating a regular heartbeat. It is thus not surprising that mutations in its sequence have been linked to life-threatening arrhythmias. Interestingly, Nav1.5 activity can also be altered by posttranslational modifications, such as tyrosine phosphorylation. Our combination of protein engineering and molecular modeling studies has revealed that the detrimental effect of a long QT3 patient mutation is only exposed when a proximal tyrosine is phosphorylated. This suggests a dynamic crosstalk between the genetic mutation and a neighboring phosphorylation, a phenomenon that could be important in other classes of proteins. Additionally, we show that phosphorylation can affect the channels sensitivity towards clinically-relevant drugs, a finding that may prove important when devising patient-specific treatment plans.

category: biophysics


title: Integrative genomics identifies lncRNA regulatory networks across 1,044 pediatric leukemias and extra-cranial solid tumors

abstract: Long non-coding RNAs (lncRNAs) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curate RNA sequencing data for 1,044 pediatric leukemia and solid tumors and integrate paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. We report a total of 2,657 robustly expressed lncRNAs across six pediatric cancers, including 1,142 exhibiting histotype-specific expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multi-dimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell-type specific transcription factors further implicated distinct histotype-specific and developmental lncRNAs. We integrated our analyses to prioritize lncRNAs for experimental validation and showed that silencing of TBX2-AS1, our top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming our computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies.

category: cancer biology


title: The Role of Human Platelet-Rich Plasma to Enhance the Differentiation from Adipose derived Mesenchymal Stem Cells into Cardiomyocyte: An Experimental Study

abstract:
BackgroundAdipose derived mesenchymal stem cells (AMSCs) offer great potential to differentiate into cardiomyocyte. However, the optimal method to maximize the proliferation and differentiation is challenging. Platelet rich plasma (PRP) which contains high levels of diverse growth factors that can stimulate stem cell proliferation and differentiation in the context of cardiac tissue regeneration

ObjectiveTo analyze the effect of PRP administration on the AMSCs differentiation into cardiomyocyte and compare to the group without PRP administration.

MethodsThis study is a true experimental randomized post-test design study. AMSCs were isolated from adipose tissues and cultured until 4 passages. The samples were divided into 3 groups, i.e. negative control (-MEM), positive control (differentiation medium), and treatment group (PRP). The assessment of GATA-4 marker expression was conducted using flowcytometry on the fifth day and cTnT was conducted using immunocytochemistry on the tenth day to determine the differentiation to cardiomyocyte. Data analysis was conducted using T-test and One-Way ANOVA on normally distributed data determined through Shapiro Wilk test.

ResultsFlowcytometry on GATA-4 expression revealed significant improvement on PRP group compared to negative and positive controls (67.04 {+/-} 4.49 vs 58.15 {+/-} 1.23 p < 0.05; 67.04 {+/-} 4.49 vs 52.96 {+/-} 2.02 p < 0.05). This was supported by the results of immunocytochemistry on troponin expression which revealed significant improvement on PRP group compared to negative and positive controls (38.13 {+/-} 5.2 vs 10.73 {+/-} 2.39 p < 0.05; 38.13 {+/-} 5.2 vs 26.00 {+/-} 0.4 p < 0.05). This was concordant to the hypothesis which stated that there was an effect of PRP administration on AMSCs differentiation into cardiomyocyte.

ConclusionPRP administration on AMSCs culture significantly improve the differentiation to cardiomyocyte measured by GATA-4 and cTnT expressions.

category: molecular biology


title: Manual Annotation of Genes within Drosophila Species: the Genomics Education Partnership protocol

abstract: Annotating the genomes of multiple organisms allows us to study their genes as well as the evolution of those genes. While many eukaryotic genome assemblies already include computational gene predictions, these predictions can benefit from review and refinement through manual gene annotation. The Genomics Education Partnership (GEP; thegep.org) has developed an annotation protocol for protein-coding genes that enables undergraduate students and other researchers to create high-quality gene annotations that can be utilized in subsequent scientific investigations. For example, this protocol has been utilized by the GEP faculty to engage undergraduate students in the comparative annotation of genes involved in the insulin signaling pathway in 28 Drosophila species, using D. melanogaster as the informant genome. Students construct gene models using multiple lines of computational and experimental evidence including expression data (e.g., RNA-Seq), sequence similarity (e.g., BLAST, multiple sequence alignments), and computational gene predictions. For quality control, each gene is annotated by at least two students working independently, followed by reconciliation of the submitted gene models by a more experienced student. This article provides an overview of the annotation protocol and describes how discrepancies in student submitted gene models are resolved to produce a final, high-quality gene set suitable for subsequent analyses. This annotation protocol can be adapted to other scientific questions (e.g., expansion of the Drosophila Muller F element) and other species (e.g., parasitoid wasps) to provide additional opportunities for undergraduate students to participate in genomics research. These student annotation efforts can substantially improve the quality of gene annotations in publicly available genomic databases.

category: genomics


title: SOMDE: A scalable method for identifying spatially variable genes with self-organizing map

abstract:
Recent developments of spatial transcriptomic sequencing technologies provide powerful tools for understanding cells in the physical context of tissue micro-environments. A fundamental task in spatial gene expression analysis is to identify genes with spatially variable expression patterns, or spatially variable genes (SVgenes). Several computational methods have been developed for this task. Their high computational complexity limited their scalability to the latest and future large-scale spatial expression data.

We present SOMDE, an efficient method for identifying SVgenes in large-scale spatial expression data. SOMDE uses self-organizing map (SOM) to cluster neighboring cells into nodes, and then uses a Gaussian Process to fit the node-level spatial gene expression to identify SVgenes. Experiments show that SOMDE is about 5-50 times faster than existing methods with comparable results. The adjustable resolution of SOMDE makes it the only method that can give results in [~]5 minutes in large datasets of more than 20,000 sequencing sites. SOMDE is available as a python package on PyPI at https://pypi.org/project/somde free for academic use.


title: Signal processing capacity of the cellular sensory machinery regulates the accuracy of chemotaxis under complex cues

abstract: Chemotaxis is ubiquitous in many biological processes, but it still remains elusive how cells sense and decipher multiple chemical cues. In this study, we postulate a hypothesis that the chemotactic performance of cells under complex cues is regulated by the signal processing capacity of the cellular sensory machinery. The underlying rationale is that cells in vivo should be able to sense and process multiple chemical cues, whose magnitude and compositions are entangled, to determine their migration direction. We experimentally show that the combination of TGF-{beta} and EGF suppresses the chemotactic performance of cancer cells using independent receptors to sense the two cues. Based on this observation, we develop a biophysical framework suggesting that the antagonism is caused by the saturation of the signal processing capacity, but not by the mutual repression. Our framework suggests the significance of the signal processing capacity in the cellular sensory machinery.


title: Stimulus-dependent representational drift in primary visual cortex

abstract: To produce consistent sensory perception, neurons must maintain stable representations of sensory input. However, neurons in many regions exhibit progressive drift across days. Longitudinal studies have found stable responses to artificial stimuli across sessions in primary sensory areas, but it is unclear whether this stability extends to naturalistic stimuli. We performed chronic 2-photon imaging of mouse V1 populations to directly compare the representational stability of artificial versus naturalistic visual stimuli over weeks. Responses to gratings were highly stable across sessions. However, neural responses to naturalistic movies exhibited progressive representational drift across sessions. Differential drift was present across cortical layers, in inhibitory interneurons, and could not be explained by differential response strength or higher order stimulus statistics. However, representational drift was accompanied by similar differential changes in local population correlation structure. These results suggest representational stability in V1 is stimulus-dependent and related to differences in preexisting circuit architecture of co-tuned neurons.


title: Cognitive and GABAergic protection by hAPPSw,Ind and hTauVLW coexpression: a model of Non-Demented with Alzheimer's disease Neuropathology

abstract: Alzheimers disease comprises amyloid-{beta} (A{beta}) and hyperphosphorylated Tau (P-Tau) accumulation, imbalanced neuronal activity, aberrant oscillatory rhythms, and cognitive deficits. Non-Demented with Alzheimers disease Neuropathology (NDAN) defines a novel clinical entity with A{beta} and Tau pathologies, but preserved cognition. The mechanisms underlying such neuroprotection remain undetermined and animal models are currently unavailable for NDAN. We show that J20/VLW mice, accumulating A{beta} and P-Tau, exhibit preserved hippocampal rhythmic activity and cognition, altered in J20 and VLW animals. Furthermore, we show that coexistence with A{beta} renders a particular P-Tau signature in hippocampal interneurons. The GABAergic septohippocampal pathway, responsible for hippocampal rhythmic activity, is preserved in J20/VLW mice, in contrast to single mutants. Our data highlight J20/VLW mice as a suitable animal model to understand the mechanisms driving cognitive preservation in NDAN and suggest that a differential P-Tau pattern in hippocampal interneurons prevents GABAergic septohippocampal innervation loss and alterations in local field potentials, avoiding cognitive deficits.


title: Evaluation of seasonal dynamics of fungal DNA assemblages in a flow-regulated stream in a restored forest using eDNA metabarcoding

abstract: Investigation of seasonal variation in fungal communities is essential for understanding biodiversity and ecosystem functions. However, the conventional sampling method, with substrate removal and high spatial heterogeneity of community composition, makes surveying the seasonality of fungal communities challenging. Recently, water environmental DNA (eDNA) analysis has been explored for its utility in biodiversity surveys. In this study, we assessed whether the seasonality of fungal communities can be detected by monitoring eDNA in a forest stream. We conducted monthly water sampling in a forest stream over two years and used DNA metabarcoding to identify fungal eDNA. The stream water contained DNA from functionally diverse aquatic and terrestrial fungi, such as plant decomposers, parasites, and mutualists. The variation in the fungal assemblage showed a regular annual periodicity, meaning that the assemblages in a given season were similar, irrespective of the year or sampling. Furthermore, the strength of the annual periodicity varied among functional groups. Our results suggest that forest streams may act as a "trap" for terrestrial fungal DNA derived from different habitats, allowing the analysis of fungal DNA in stream water to provide information about the temporal variation in fungal communities in both the aquatic and the surrounding terrestrial ecosystems.

category: ecology


title: Developmental effects on pattern visual evoked potentials characterized by principal component analysis

abstract:
PurposePeak amplitude and peak latency in the pattern reversal visual evoked potential (prVEP) vary with maturation. We considered that principal component analysis (PCA) may be used to describe age-related variation over the entire prVEP time course and provide a means of modeling and removing variation due to developmental age.

MethodsprVEP was recorded from 155 healthy subjects ages 11-19 years at two timepoints. We created a model of the prVEP by identifying principal components (PCs) that explained >95% of the variance in a "training" dataset of 40 subjects. We examined the ability of the PCs to explain variance in an age- and sex-matched "validation" dataset (n=40) and calculated the intra-subject reliability of the PC coefficients between the two timepoints. We explored the effect of subject age and sex upon the PC coefficients.

ResultsSeven PCs accounted for 96.0% of the variability of the training dataset and 90.5% of the variability in the validation dataset with good within-subject reliability across timepoints (R>0.7 for all PCs). The PCA model revealed narrowing and amplitude reduction of the P100 peak with maturation, and a broader and smaller P100 peak in males compared to females.

ConclusionsPCA is a generalizable, reliable, and unbiased method of analyzing prVEP. The PCA model revealed changes across maturation and biological sex not fully described by standard peak analysis.

Translational relevanceWe describe a novel application of PCA to characterize developmental changes of prVEP in youth that can be used to compare healthy and pathologic pediatric cohorts.


title: Temporal Scaling of Human Scalp-Recorded Potentials During Interval Estimation

abstract: Much of human behaviour is governed by common processes that unfold over varying timescales. Standard event-related potential analysis assumes fixed-latency responses relative to experimental events. However, recent single unit recordings in animals have revealed neural activity scales to span different durations during behaviours demanding flexible timing. Here, we employed a general linear modelling approach using a novel combination of fixed-duration and variable-duration regressors to unmix fixed-time and scaled-time components in human magneto/electroencephalography (M/EEG) data. We use this to reveal consistent temporal scaling of human scalp-recorded potentials across four independent EEG datasets, including interval perception, production, prediction and value-based decision making. Between-trial variation in the temporally scaled response predicts between-trial variation in subject reaction times, demonstrating the relevance of this temporally scaled signal for temporal variation in behaviour. Our results provide a general approach for studying flexibly timed behaviour in the human brain.


title: Temporal Scaling of Human Scalp-Recorded Potentials


title: Identification of a family of Vibrio type III secretion system effectors that contain a conserved serine/threonine kinase domain

abstract:
Vibrio parahaemolyticus is a marine Gram-negative bacterium that is a leading cause of seafood-borne gastroenteritis. Pandemic strains of V. parahaemolyticus rely on a specialized protein secretion machinery known as the type III secretion system 2 (T3SS2) to cause disease. The T3SS2 mediates the delivery of effector proteins into the cytosol of infected cells, where they subvert multiple cellular pathways. Here, we identify a new T3SS2 effector protein encoded by VPA1328 (VP_RS21530) in V. parahaemolyticus RIMD2210633. Bioinformatic analysis revealed that VPA1328 is part of a larger family of uncharacterized T3SS effector proteins with homology to the VopG effector protein in V. cholerae AM-19226. These VopG-like proteins are found in many but not all T3SS2 gene clusters and are distributed among diverse Vibrio species including V. parahaemolyticus, V. cholerae, V. mimicus, and V. diabolicus and also in Shewanella baltica. Structure-based prediction analyses uncovered the presence of a conserved C-terminal kinase domain in VopG orthologs, similar to the serine/threonine kinase domain found in the NleH family of T3SS effector proteins. However, in contrast to NleH effector proteins, in tissue culture-based infections, VopG did not impede host cell death or suppress IL-8 secretion, suggesting a yet undefined role for VopG during V. parahaemolyticus infection. Collectively, our work reveals that VopG effector proteins, a new family of likely serine/threonine kinases, is widely distributed in the T3SS2 effector armamentarium among marine bacteria.

IMPORTANCEVibrio parahaemolyticus is the leading bacterial cause of seafood-borne gastroenteritis worldwide. The pathogen relies on a type III secretion system to deliver a variety of effector proteins into the cytosol of infected cells to subvert cellular function. In this study, we identified a novel Vibrio parahaemolyticus effector protein that is similar to the VopG effector of Vibrio cholerae. VopG-like effectors were found in diverse Vibrio species and contain a conserved serine/threonine kinase domain that bears similarity to the kinase domain in the EHEC and Shigella NleH effectors that manipulate host cell survival pathways and host immune responses. Together our findings identify a new family of Vibrio effector proteins and highlight the role of horizontal gene transfer events among marine bacteria in shaping T3SS gene clusters.

category: microbiology


title: A Synthetic Transcription Platform for Programmable Gene Expression in Mammalian Cells

abstract: Precise, scalable, and sustainable control of genetic and cellular activities in mammalian cells is key to developing precision therapeutics and smart biomanufacturing. We created a highly tunable, modular, versatile CRISPR-based synthetic transcription system for the programmable control of gene expression and cellular phenotypes in mammalian cells. Genetic circuits consisting of well-characterized libraries of guide RNAs, binding motifs of synthetic operators, transcriptional activators, and additional genetic regulatory elements expressed mammalian genes in a highly predictable and tunable manner. We demonstrated the programmable control of reporter genes episomally and chromosomally, with up to 25-fold more activity than seen with the EF1 promoter, in multiple cell types. We used these circuits to program the secretion of human monoclonal antibodies and to control T-cell effector function marked by interferon-{gamma} production. Antibody titers and interferon-{gamma} concentrations significantly correlated with synthetic promoter strengths, providing a platform for programming gene expression and cellular function in diverse applications.

category: synthetic biology


title: Hippocampal cells multiplex positive and negative engrams

abstract:
The hippocampus is involved in processing a variety of mnemonic computations specifically the spatiotemporal components and emotional dimensions of contextual memory.1-3 Recent studies have demonstrated vast structural and functional heterogeneity along the dorsal-ventral axis1, 5 of the hippocampus, and while much is known about how the dorsal hippocampus processes spatial- temporal content, much less is known about whether or not the ventral hippocampus (vHPC) partitions positive and negative experiences into distinct sets of cells.4-9 Here, we combine transgenic and all-virus based activity-dependent tagging strategies to visualize multiple valence- specific engrams in the vHPC and demonstrate two partially segregated cell populations and projections that respond to positive and negative experiences. Next, using an RNA sequencing approach, we find that vHPC positive and negative engram cells display distinct transcriptional programs compared to a neutral engram population. Additionally, while optogenetic manipulation of tagged cell bodies in vHPC is not sufficient to drive preference or avoidance, stimulation of tagged vHPC terminals projecting to the amygdala and nucleus accumbens (NAc), but not the prefrontal cortex (PFC), drives preference and avoidance. These terminals can also undergo a "switch" or "reset" in their capacity to drive either, thereby demonstrating their flexible contributions to behavior. We conclude that the vHPC contains genetically, cellularly, and behaviorally distinct populations of cells processing positive and negative memory engrams. Together, our findings provide a novel means by which to visualize multiple engrams within the same brain and point to their unique genetic signatures as reference maps for the future development of new therapeutic strategies.

One sentence summaryThe hippocampus contains neurons that correspond to positive and negative engrams, which are segregated by their molecular, cellular, and projection-specific features


title: Target2DeNovoDrugPropMax : a novel programmatic tool incorporating deep learning and in silico methods for automated de novo drug design for any target of interest

abstract:
The past decade has seen a surge in the range of application data science, machine learning, deep learning, and AI methods to drug discovery. The presented work involves an assemblage of a variety of AI methods for drug discovery along with the incorporation of in silico techniques to provide a holistic tool for automated drug discovery. When drug candidates are required to be identified for a particular drug target of interest, the user is required to provide the tool target signatures in the form of an amino acid sequence or its corresponding nucleotide sequence. The tool collects data registered on PubChem required to perform an automated QSAR and with the validated QSAR model, prediction and drug lead generation are carried out. This protocol we call Target2Drug. This is followed by a protocol we call Target2DeNovoDrug wherein novel molecules with likely activity against the target are generated de novo using a generative LSTM model. It is often required in drug discovery that the generated molecules possess certain properties like drug-likeness, and therefore to optimize the generated de novo molecules toward the required drug-like property we use a deep learning model called DeepFMPO, and this protocol we call Target2DeNovoDrugPropMax. This is followed by the fast automated AutoDock-Vina based in silico modeling and profiling of the interaction of optimized drug leads and the drug target. This is followed by an automated execution of the Molecular Dynamics protocol that is also carried out for the complex identified with the best protein-ligand interaction from the AutoDock-Vina based virtual screening. The results are stored in the working folder of the user. The code is maintained, supported, and provide for use in the following GitHub repository

https://github.com/bengeof/Target2DeNovoDrugPropMax


title: Optogenetic inhibition and activation of Rac and Rap1 using a modified iLID system

abstract: The small GTPases Rac1 and Rap1 can fulfill multiple cellular functions because their activation kinetics and localization are precisely controlled. To probe the role of their spatio-temporal dynamics, we generated optogenetic tools that activate or inhibit endogenous Rac and Rap1 in living cells. An improved version of the light-induced dimerization (iLID) system was used to control plasma membrane localization of protein domains that specifically activate or inactivate Rap1 and Rac (Tiam1 and Chimerin for Rac, RasGRP2 and Rap1GAP for Rap1). Irradiation yielded a 50-230% increase in the concentration of these domains at the membrane, leading to effects on cell morphodynamics consistent with the known roles of Rac1 and Rap1.


title: A viral histone-like protein exploits antagonism between linker histones and HMGB proteins to obstruct the cell cycle

abstract: Virus infection necessarily requires redirecting cellular resources towards viral progeny production. Adenovirus encodes the histone-like protein VII that causes catastrophic global reorganization of host chromatin to promote virus infection. Protein VII recruits the family of high mobility group box (HMGB) proteins to chromatin along with the histone chaperone SET. As a consequence of this recruitment, we find that protein VII causes chromatin-depletion of several linker histone H1 isoforms. The relationship between linker histone H1 and the functionally opposite HMGB proteins is critical for higher order chromatin structure. However, the physiological consequences of perturbing this relationship are largely unknown. Here, we employ complementary systems in Saccharomyces cerevisiae and human cells to demonstrate that adenovirus protein VII disrupts the H1-HMGB balance to obstruct the cell cycle. We find that protein VII causes an accumulation of G2/M cells both in yeast and human systems, underscoring the high conservation of this chromatin vulnerability. In contrast, adenovirus E1A and E1B proteins are well-established to override cell cycle regulation and promote transformation of human cells. Strikingly, we find that protein VII obstructs the cell cycle even in the presence of E1A and E1B, suggesting that protein VII-directed cell cycle disruption ensures host resources are directed towards viral proliferation during infection. Together, our results demonstrate that protein VII targets H1-HMGB1 antagonism to obstruct cell cycle progression, revealing an unexpected chromatin vulnerability exploited for viral benefit.


title: Evolution of the recombination regulator PRDM9 in minke whales

abstract: We explored the structure and variability of the Prdm9 gene, which codes for the PRDM9 protein, in samples of the minke whales from the Atlantic, Pacific and Southern Oceans. The PRDM9 protein controls the reshuffling of parental genomes in most metazoans and we show that minke whale possess all the features characteristic of PRDM9-directed recombination initiation, including complete KRAB, SSXRD and SET domains and a rapidly evolving array of C2H2-type-Zincfingers (ZnF). We uncovered eighteen novel PRDM9 variants and evidence of rapid evolution, particularly at DNA-recognizing positions that evolve under positive selection. At different geographical scales, we observed extensive Prdm9 allelic diversity in Antarctic minke whales (Balaenoptera bonarensis) that, conversely, lack observable population differentiation in mitochondrial DNA and microsatellites. In contrast, a single PRDM9 variant is shared between all Common Minke whales and even across subspecies boundaries of North Atlantic (B. a. acutorostrata) and North Pacific (B. a. scammoni) minke whale, which do show clear population differentiation. PRDM9 variation of whales predicts distinct recombination initiation landscapes genome-wide, which has possible consequences for speciation.


title: Measuring Compound Eye Optics with Microscope and MicroCT Images

abstract: The arthropod compound eye is the most prevalent eye type in the animal kingdom, with an impressive range of shapes and sizes. Studying its natural range of morphologies provides insight into visual ecology, development, and evolution. In contrast to the camera-type eyes we possess, external structures of compound eyes often reveal resolution, sensitivity, and field of view if the eye is spherical. Non-spherical eyes, however, require measuring internal structures using imaging technology like MicroCT (CT). CT is a burgeoning 3D X-Ray imaging technique that has been used to image arthropod muscles, brains, ocelli, and eyes. Thus far, there is no efficient tool to automate characterizing compound eye optics in either 2D or 3D data. We present two open-source programs: (1) the ommatidia detecting algorithm (ODA), which automatically measures ommatidia count and diameter in 2D images, and (2) a CT pipeline (ODA-3D), which calculates anatomical acuity, sensitivity, and field of view across the eye by applying the ODA to 3D data. We validate these algorithms on images, images of replicas, and CT scans from eyes of ants, fruit flies, moths, and a bee.


title: Bridges, not walls: the hippocampus builds narrative memories across distant events

abstract: Lifes events are scattered throughout time, yet we often recall different events in the context of an integrated narrative. Prior research suggests that the hippocampus, which supports memory for past events, can sometimes support integration of overlapping associations or separate events, but the conditions which lead to hippocampus-dependent memory integration are unclear. We used functional brain imaging to test whether the ability to form a larger narrative (narrative coherence) drives hippocampal memory integration. During encoding of fictional stories, the hippocampus supported patterns of activity, including activity at boundaries between events, which were more similar between distant events that formed one coherent narrative, compared with events taken from unrelated narratives. Twenty-four hours later, the hippocampus preferentially supported detailed recall of coherent narrative events, through reinstatement of hippocampal activity patterns from encoding. These findings reveal a key function of the hippocampus: the dynamic integration of events into a narrative structure for memory.


title: Building bridges: the hippocampus constructs narrative memories across distant events


title: A new method shows that a minority of liver-localized CD8 T cells display hard-to-detect attraction to Plasmodium-infected hepatocytes

abstract: Malaria is a disease caused by Plasmodium parasites, resulting in over 200 million infections and 400,000 deaths every year. A critical step of malaria infection is when sporozoites, injected by mosquitoes, travel to the liver and form liver stages. Malaria vaccine candidates which induce large numbers of malaria-specific CD8 T cells in mice are able to eliminate all liver stages, preventing fulminant malaria. However, how CD8 T cells find all parasites in 48 hours of the liver stage lifespan is not well understood. Using intravital microscopy of murine livers, we generated unique data on T cell search for malaria liver stages within a few hours after infection. To detect attraction of T cells to an infection site, we used the von Mises-Fisher distribution in 3D, similar to the 2D von Mises distribution previously used in ecology. Our results suggest that the vast majority (70-95%) of malaria-specific and non-specific liver-localized CD8 T cells did not display attraction towards the infection site, suggesting that the search for malaria liver stages occurs randomly. However, a small fraction (15-20%) displayed weak but detectable attraction towards parasites which already had been surrounded by several T cells. We found that speeds and turning angles correlated with attraction, suggesting that understanding mechanisms that determine the speed of T cell movement in the liver may improve the efficacy of future T cell-based vaccines. Stochastic simulations suggest that a small movement bias towards the parasite dramatically reduces the number of CD8 T cells needed to eliminate all malaria liver stages, but to detect such attraction by individual cells requires data from long imaging experiments which are not currently feasible. Importantly, as far as we know this is the first demonstration of how activated/memory CD8 T cells might search for the pathogen in nonlymphoid tissues a few hours after infection. We have also established a framework for how attraction of individual T cells towards a location in 3D can be rigorously evaluated.

category: immunology


title: A new method based on the von Mises-Fisher distribution shows that a minority of liver-localized CD8 T cells display hard-to-detect attraction to Plasmodium-infected hepatocytes


title: Uncovering complementary sets of variants for predicting quantitative phenotypes

abstract:
MotivationGenome-wide association studies show that variants in individual genomic loci alone are not sufficient to explain the heritability of complex, quantitative phenotypes. Many computational methods have been developed to address this issue by considering subsets of loci that can collectively predict the phenotype. This problem can be considered a challenging instance of feature selection in which the number of dimensions (loci that are screened) is much larger than the number of samples. While currently available methods can achieve decent phenotype prediction performance, they either do not scale to large datasets or have parameters that require extensive tuning.

ResultsWe propose a fast and simple algorithm, Macarons, to select a small, complementary subset of variants by avoiding redundant pairs that are in linkage disequilibrium. Our method features two interpretable parameters that control the time/performance trade-off without requiring parameter tuning. In our computational experiments, we show that Macarons consistently achieves similar or better prediction performance than state-of-the-art selection methods while having a simpler premise and being at least 2 orders of magnitude faster. Overall, Macarons can seamlessly scale to the human genome with ~107 variants in a matter of minutes while taking the dependencies between the variants into account.

ConclusionMacarons can offer a reasonable trade-off between phenotype predictivity, runtime and the complementarity of the selected subsets. The framework we present can be generalized to other high-dimensional feature selection problems within and beyond biomedical applications.

AvailabilityMacarons is implemented in Matlab and the source code is available at: https://github.com/serhan-yilmaz/macarons


title: Physiological and transcriptomic response of Medicago truncatula to high and low benefit mycorrhizal fungi

abstract: Arbuscular mycorrhizal (AM) fungi form a root endosymbiosis with many agronomically important crop species and both enhance the ability of their host to obtain nutrients from the soil and increase host tolerance to biotic and abiotic stressors. However, AM fungal species differ in the benefits they provide to their host plants. Here, we examined the putative molecular mechanisms involved in the regulation of the physiological response of Medicago truncatula to either Rhizophagus irregularis or Glomus aggregatum, a high or a low benefit AM fungus, respectively. Colonization with R. irregularis led to higher growth and nutrient uptake benefits than the colonization with G. aggregatum. These benefits were linked to an elevated expression in the roots of genes involved in strigolactone biosynthesis (NSP1, NSP2, CCD7, and MAX1a), mycorrhiza-induced phosphate (PT8), ammonium (AMT2;3), and nitrate (NPF4.12) transporters and the putative ammonium transporter NIP1;5. R. irregularis also stimulated the expression of photosynthesis-related genes in the shoot and the upregulation of the sugar transporters SWEET1.2, SWEET3.3 and SWEET 12 and the lipid biosynthesis gene RAM2 in the roots. In contrast, G. aggregatum induced the expression of biotic stress defense response genes in the shoots and several genes associated with abiotic stress in the roots. This suggests that either the host perceives colonization by G. aggregatum as a pathogen attack or that G. aggregatum can prime host defense responses. Our findings reveal novel insights into the molecular mechanisms that control the host plant response to colonization with high- and low-benefit arbuscular mycorrhizal fungal symbionts.

category: plant biology


title: Pioneer activity of an oncogenic fusion transcription factor at inaccessible chromatin

abstract: Recent characterizations of pioneer transcription factors have led to new insights into their structures and patterns of chromatin recognition that are instructive for understanding their role in cell fate commitment and transformation. Intersecting with these basic science concepts, the identification of pioneer factors (PFs) fused together as driver translocations in childhood cancers raises questions of whether these fusions retain the fundamental ability to invade repressed chromatin, consistent with their monomeric PF constituents. In this study, we define the cellular and chromatin localization of the translocation, PAX3-FOXO1, an oncogenic driver of childhood rhabdomyosarcoma (RMS), derived from a genetic fusion of PFs. To quantitatively define its chromatin-targeting functions and capacity to drive epigenetic reprogramming, we developed a new method for ChIP-seq with per-cell normalization (pc-ChIP-seq). Our quantitative localization studies address structural variation in RMS genomes and reveal novel insights into heterochromatin localization of PAX3-FOXO1. From these studies, we report novel pioneer function for the major driver oncogene in RMS, with repressed chromatin binding and nucleosome-motif targeting in human cells.


title: Toward a universal model for spatially structured populations

abstract: A key question in evolution is how likely a mutant is to take over. This depends on natural selection and on stochastic fluctuations. Population spatial structure can impact mutant fixation probabilities. We introduce a model for structured populations on graphs that generalizes previous ones by making migrations independent of birth and death. We demonstrate that by tuning migration asymmetry, the star graph transitions from amplifying to suppressing natural selection. Our results are universal in the sense that they do not hinge on a modeling choice of microscopic dynamics or update rules. Instead, they depend on migration asymmetry, which can be experimentally tuned and measured.

category: evolutionary biology


title: Gene family evolution underlies cell type diversification in the hypothalamus of teleosts

abstract: Hundreds of cell types form the vertebrate brain, but it is largely unknown how similar these cellular repertoires are between or within species, or how cell type diversity evolves. To examine cell type diversity across and within species, we performed single-cell RNA sequencing of [~]130,000 hypothalamic cells from zebrafish (Danio rerio) and surface- and cave-morphs of Mexican tetra (Astyanax mexicanus). We found that over 75% of cell types were shared between zebrafish and Mexican tetra, which last shared a common ancestor over 150 million years ago. Orthologous cell types displayed differential paralogue expression that was generated by sub-functionalization after genome duplication. Expression of terminal effector genes, such as neuropeptides, was more conserved than the expression of their associated transcriptional regulators. Species-specific cell types were enriched for the expression of species-specific genes, and characterized by the neo-functionalization of members of recently expanded or contracted gene families. Within species comparisons revealed differences in immune repertoires and transcriptional changes in neuropeptidergic cell types associated with genomic differences between surface- and cave-morphs. The single-cell atlases presented here are a powerful resource to explore hypothalamic cell types, and reveal how gene family evolution and the neo- and sub-functionalization of paralogs contribute to cellular diversity.


title: Identification of a Golgi-localized peptide encoded by the CENP-R transcript

abstract: Prior work has identified signal sequences and motifs that are necessary and sufficient to target proteins to specific subcellular regions and organelles such as the plasma membrane, nucleus, Endoplasmic Reticulum, and mitochondria. In contrast, minimal sequence motifs that are sufficient for Golgi localization remain largely elusive. In this work, we identified a 37 amino acid alternative open reading frame (altORF) within the mRNA of the centromere protein, CENP-R. This altORF peptide localizes specifically to the cytoplasmic surface of the Golgi apparatus. Through mutational analysis, we identify a minimal 10 amino acid sequence and a critical cysteine residue that are necessary and sufficient for Golgi localization. Pharmacological perturbations suggest that this peptide is modified through palmitoylation to promote its localization. Together, our work defines a minimal sequence that is sufficient for Golgi targeting and provide a valuable Golgi marker for live cell imaging.


title: Identification of a Golgi-localized altORF peptide reveals a minimal Golgi targeting motif


title: Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

abstract: Mutational signature analysis is commonly performed in genomic studies surveying cancer and normal somatic tissues. Here we present SigProfilerExtractor, an automated tool for accurate de novo extraction of mutational signatures for all types of somatic mutations. Benchmarking with a total of 34 distinct scenarios encompassing 2,500 simulated signatures operative in more than 60,000 unique synthetic genomes and 20,000 synthetic exomes demonstrates that SigProfilerExtractor outperforms thirteen other tools across all datasets with and without noise. For genome simulations with 5% noise, reflecting high-quality genomic datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true positive signatures while yielding more than 5-fold less false positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome sequenced and 19,184 whole-exome sequenced cancers reveals four previously missed mutational signatures. Two of the signatures are confirmed in independent cohorts with one of these signatures associating with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting mutational signatures, and several novel mutational signatures including a signature putatively attributed to direct tobacco smoking mutagenesis in bladder cancer and in normal bladder epithelium.


title: 5-HT2C receptors regulate the balance between instrumental vigour and restraint

abstract: The serotonin (5-HT) system, particularly the 5-HT2C receptor, has consistently been implicated in behavioural control. However, while some studies have focused on the role 5-HT2C receptors plays in regulating motivation to work for reward, others have highlighted its importance in response restraint. To date, it is unclear how 5-HT transmission at this receptor regulates the balance of response invigoration and restraint in anticipation of future reward. In addition, it remains to be established how 5-HT2C receptors gate the influence of internal versus cue-driven processes over reward-guided actions. To elucidate these issues, we investigated the effects of administering the 5-HT2C receptor antagonist SB242084, both systemically and directly into the nucleus accumbens core (NAcC), in rats performing a Go/No-Go task for small or large rewards. The results were compared to administration of d-amphetamine into the NAcC, which has previously been shown to promote behavioural activation. Systemic perturbation of 5-HT2C receptors - but crucially not intra-NAcC infusions - consistently boosted rats performance and instrumental vigour on Go trials when they were required to act. Concomitantly, systemic administration also reduced their ability to withhold responding for rewards on No-Go trials, particularly late in the holding period. Notably, these effects were often apparent only when the reward on offer was small. By contrast, inducing a hyperdopaminergic state in the NAcC with d-amphetamine strongly impaired response restraint on No-Go trials both early and late in the holding period, as well as speeding action initiation. Together, these findings suggest that 5-HT2C receptor transmission, outside the NAcC, shapes the vigour of ongoing goal-directed action as well as the likelihood of responding as a function of expected reward.

category: neuroscience


title: 5-HT2C receptor perturbation has bidirectional influence over instrumental vigour and restraint


title: Translational control as a novel regulator of gradient sensing and chemotropism in yeast

abstract:
The yeast mating pathway regulates haploid cell fusion in response to pheromone signaling via a mitogen-activated protein kinase (MAPK) cascade that controls directional growth (chemotropism). However, the regulators of chemotropic morphogenesis are ill-defined. By using a non-biased genome-wide screen, we identified hundreds of genes that affect mating. An additional screens identified and validated >20 novel positive and negative regulators of pheromone gradient sensing, chemotropism, shmoo development, and mating. Aside from known regulators of exocytosis and endocytosis, genes involved in translational control downstream of the G-protein-regulated pheromone and filamentous growth MAPK pathways were identified. These include the Scp160 RNA-binding protein and the Asc1, Rpl12b, and Rpl19b ribosomal proteins (RPs). Importantly, we demonstrate that pheromone treatment and G (Gpa1) activation stimulate Scp160 binding to (and inhibition of) Asc1, which acts downstream of glucose-activated G (Gpa2) on the filamentous growth pathway. Moreover, we identify both Rpl12b and Rpl19b as RP paralog-specific positive regulators of translation of mating components, including Scp160. Thus, opposing MAPK pathways may converge at the level of translational control to regulate signaling output.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=165 SRC="FIGDIR/small/422562v2_ufig1.gif" ALT="Figure 1">
View larger version (25K):
org.highwire.dtl.DTLVardef@3944e1org.highwire.dtl.DTLVardef@1d9204eorg.highwire.dtl.DTLVardef@1d7e902org.highwire.dtl.DTLVardef@2bfc6_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: cell biology


title: Plasticity of fibroblast transcriptional response to physical and biochemical cues revealed by dynamic network analysis

abstract: Data on human skin fibroblast transcriptional responses to external cues were used to reconstruct dynamic gene regulatory networks. The goal of the reconstruction was to determine dynamic network interactions (quantitative predictive relationships of mutual regulatory influences of and on transcription factor expression) from time course data on 56 transcript expression levels obtained following different external cues. The inherently under-determined nature of this problem was addressed in part by excluding putative regulatory motifs that did not appear to be functional in multiple independent experiments from different independent external perturbations. Data were obtained from a previously published experiment in which the 56 transcripts were assayed by bioluminescence in live cells cultured on substrates of varying levels of stiffness and exposed to different levels of arginylglycylaspartic acid (RGD) peptide. The inferred dynamical networks were validated via comparison of predictions to a priori known interactions from gene databases. We discovered that exposures to different substrate stiffnesses and to RGD stimulate responses that are mediated through GATA4, SMAD3/4, ETS-1, and STAT5 and other genes, which can initiate hypertrophic, fibrotic, and inflammatory responses. The developed dynamical system identification method for discovering new mechanotransduction pathways is applicable to the identification of gene regulatory networks in numerous emerging applications where time-series data on multiple state variables and from multiple external perturbations are available.

category: systems biology


title: Dietary folic acid deficiency impacts hippocampal morphology and cortical acetylcholine metabolism in adult male and female mice

abstract:
ObjectiveOne-carbon metabolism is a metabolic network that integrates nutritional signals with biosynthesis, redox homeostasis, and epigenetics. There are sex differences in hepatic one-carbon metabolism. However, it is unclear whether there are sex differences in dietary deficiencies of one-carbon metabolism in the brain. The aim of this study was to investigate if sex modulates the effects of dietary folic acid deficiency in brain tissue using a mouse model.

MethodsMale and female C57Bl/6J mice were placed on a folic acid deficient (FD) or control diet (CD) at six weeks of age. Mice were maintained on these diets for six months, after which animals were euthanized and brain tissue and serum were collected for analysis. Serum folate levels were measured. In brain tissue, hippocampal volume and morphology including Cornu Ammonis 1 and 3 (CA1; CA3), and dentate gyrus thickness were measured. Apoptosis within the hippocampus was assessed using active caspase-3 immunofluorescence staining. Additionally, cortical acetylcholine metabolism was measured in brain tissue using immunofluorescence staining of acetylcholinesterase (AChE), or choline acetyltransferase (ChAT), and neuronal nuclei (NeuN).

ResultsMale and female FD mice had reduced serum levels of folate. Both males and females maintained on a FD showed a decrease in the thickness of the hippocampal CA1-CA3 region. Interestingly, there was a sex difference in the levels of active caspase-3 within the CA3 region of the hippocampus. In cortical tissue, there were increased levels of neuronal ChAT and reduced levels of AChE in FD females and male mice.

ConclusionsThe results indicated that FD impacts hippocampal morphology and cortical neuronal acetylcholine metabolism. The data from our study indicate that there was only one sex difference and that was in hippocampal apoptosis. Our study provides little evidence that sex modulates the effects of dietary folate deficiency on hippocampal morphology and cortical neuronal acetylcholine metabolism.


title: Analysing economic costs of invasive alien species with the invacost R package

abstract:
O_LIThe reported costs of invasive alien species from the global database InvaCost are heterogenous and cover different spatio-temporal scales. A standard procedure for aggregating invasive species cost estimates is necessary to ensure the repeatability and comparativeness of studies.
C_LIO_LIWe introduce here the invacost R package, an open-source software designed to query and analyse the InvaCost database. We illustrate this package and its framework with cost data associated with invasive alien invertebrates.
C_LIO_LIFirst, the invacost package provides updates of this dynamic database directly in the analytical environment R. Second, it helps understand the heteregoneous nature of monetary cost data for invasive species, processes to harmonize the data, and the inherent biases associated with such data. Third, it readily provides complementary methods to investigate the costs of invasive species at different scales, all the while accounting for econometric statistical issues.
C_LIO_LIThis tool will be useful for scientists working on invasive alien species, by (i) facilitating access to and use of this multi-disciplinary data resource and (ii) providing a standard procedure which will facilitate reproducibility and comparability among studies, one of the major critics of this topic until now. It should facilitate further interdisciplinary works including economists and invasion ecology researchers.
C_LI


title: Large-scale geographic survey provides insights into the colonization history of a major aphid pest on its cultivated apple host in Europe, North America and North Africa

abstract: AO_SCPLOWBSTRACTC_SCPLOWWith frequent host shifts involving the colonization of new hosts across large geographical ranges, crop pests are good models for examining the mechanisms of rapid colonization. The microbial partners of pest insects may also be involved in or affected by colonization processes, which has been little studied so far. We investigated the demographic history of the rosy apple aphid, Dysaphis plantaginea, a major pest of the cultivated apple (Malus domestica) in Europe, North Africa and North America, as well as the diversity of its microbiota. We genotyped a comprehensive sample of 714 colonies from Europe, Morocco and the US using mitochondrial (CytB and CO1), bacterial (16s rRNA and TrnpB), and 30 microsatellite markers. We detected five populations spread across the US, Morocco, Western and Eastern Europe and Spain. Populations showed weak genetic differentiation and high genetic diversity, except the ones from Morocco and North America that are likely the result of recent colonization events. Coalescent-based inferences revealed high levels of gene flow among populations during the colonization but did not allow determining the sequence of colonization of Europe, North America and Morroco by D. plantaginea, likely because of the weak genetic differentiation and the occurrence of gene flow among populations. We found that D. plantaginea rarely hosts other endosymbiotic bacteria than its obligate nutritional symbiont Buchnera aphidicola. This suggests that secondary endosymbionts did not play an important role in the rapid spread of the rosy apple aphid. These findings have fundamental importance for understanding pest colonization processes and implications for sustainable pest control programs.


title: Identification of bis-benzylisoquinoline alkaloids as SARS-CoV-2 entry inhibitors from a library of natural products in vitro

abstract: Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major public health issue. To screen for antiviral drugs for COVID-19 treatment, we constructed a SARS-CoV-2 spike (S) pseudovirus system using an HIV-1-based lentiviral vector with a luciferase reporter gene to screen 188 small potential antiviral compounds. Using this system, we identified nine compounds, specifically, bis-benzylisoquinoline alkaloids, that potently inhibited SARS-CoV-2 pseudovirus entry, with EC50 values of 0.1-10 M. Mechanistic studies showed that these compounds, reported as calcium channel blockers (CCBs), inhibited Ca2+-mediated membrane fusion and consequently suppressed coronavirus entry. These candidate drugs showed broad-spectrum efficacy against the entry of several coronavirus pseudotypes (SARS-CoV, MERS-CoV, SARS-CoV-2 [S-D614, S-G614, N501Y.V1 and N501Y.V2]) in different cell lines (293T, Calu-3, and A549). Antiviral tests using native SARS-CoV-2 in Vero E6 cells confirmed that four of the drugs (SC9/cepharanthine, SC161/hernandezine, SC171, and SC185/neferine) reduced cytopathic effect and supernatant viral RNA load. Among them, cepharanthine showed the strongest anti-SARS-CoV-2 activity. Collectively, this study offers new lead compounds for coronavirus antiviral drug discovery.


title: Transcription of macromolecular N degrading genes provides a window into organic nitrogen decomposition in soil

abstract:
BackgroundNitrogen (N) is commonly limiting in soil in part because most N is present in macromolecular organic compounds and not directly available to plants. Using extracellular enzymes, the soil microbial community present near roots (rhizosphere) is largely responsible for transforming organic substrates, which yields nitrogen in forms that are plant-available. Microbial genetic potential (genes) and the cycling of inorganic nitrogen pools have been extensively documented and widely explored in rhizospere soil. But there is scarce information on microbial gene expression for macromolecular nitrogen decomposition and how it relates in space and time to life stages of plant roots. Here we use a suite of time-resolved metatranscriptomes from rhizosphere and bulk soil to follow bacterial and fungal extracellular protease and chitinase expression during rhizosphere aging. We also explore the effect of adding plant root litter as a distinct source of macromolecular carbon and nitrogen.

ResultsExpression of extracellular proteases in rhizosphere soil increased with plant growth in the absence of litter, whereas chitinase (chit1) was highly upregulated in the detritusphere. Structural groups of proteases were dominated by serineproteases, despite the higher abundance of metalloproteases previously observed in soil and aquatic systems. Extracellular proteases of Betaprotebacteria were more highly expressed in the rhizosphere, whereas those of Deltaroteobacteria and Fungi responded strongly to the presence of litter. We identified distinct functional guilds of taxa specializing in decomposition of proteins in the rhizosphere, detritusphere and in the vicinity of aging roots. We also identify a sub-guild that appears to specialize in protein decomposition in the presence of growing roots as well as litter that increases its activity in aging rhizosphere soil, which may be particularly N limited. Taxonomic membership of these guilds partially overlapped with functional guilds identified by carbon decomposing enzymes, making up multifunctional guilds.

ConclusionsWe found organic N decomposers have distinct spatiotemporal preferences that are selected during plant growth. Rhizosphere microbes which could potentially benefit plants by degrading macromolecular nitrogen were not comparatively active in bulk soil, even in the presence of root litter, implying that they are enhanced near plant roots and less competitive in bulk soils.


title: Pathogenic and Indigenous Denitrifying Bacteria are Transcriptionally Active and Key Multi-Antibiotic Resistant Players in Wastewater Treatment Plants

abstract: The global rise and spread of antibiotic resistance greatly challenge the treatment of bacterial infections. Wastewater treatment plants (WWTPs) harbor and discharge antibiotic resistance genes (ARGs) as environmental contaminants. However, the knowledge gap on the host identity, activity and functionality of ARGs limits transmission and health risk assessment of WWTPs resistome. Hereby, a genome-centric quantitative metatranscriptomic approach was exploited to realize high-resolution qualitative and quantitative analyses of bacterial hosts of ARGs (i.e., multi-resistance, pathogenicity, activity and niches) throughout 12 urban WWTPs. We found that [~]45% of 248 recovered genomes expressed ARGs against multiple classes of antibiotics, among which bacitracin and aminoglycoside resistance genes in Proteobacteria was the most prevalent scenario. Both potential pathogens and indigenous denitrifying bacteria were transcriptionally active hosts of ARGs. The almost unchanged relative expression levels of ARGs in the most resistant populations (66.9%) and the surviving ARG hosts including globally emerging pathogens (e.g., Aliarcobacter cryaerophilus) in treated WWTP effluent prioritizes future examination on the health risks related with resistance propagation and human exposure in the receiving environment.


title: Chromatin remodeler Arid1a regulates subplate neuron identity and wiring of cortical connectivity

abstract:
Loss-of-function mutations in chromatin remodeler gene ARID1A are a cause of Coffin-Siris syndrome, a developmental disorder characterized by dysgenesis of corpus callosum. Here, we characterize Arid1a function during cortical development and find unexpectedly selective roles for Arid1a in subplate neurons. Subplate neurons (SPNs), strategically positioned at the interface of cortical grey and white matter, orchestrate multiple developmental processes indispensable for neural circuit wiring. We find that pan-cortical deletion of Arid1a leads to extensive mistargeting of intracortical axons and agenesis of corpus callosum. Sparse Arid1a deletion, however, does not autonomously misroute callosal axons, implicating non-cell autonomous Arid1a functions in axon guidance. Supporting this possibility, the ascending axons of thalamocortical neurons, which are not autonomously affected by cortical Arid1a deletion, are also disrupted in their pathfinding into cortex and innervation of whisker barrels. Coincident with these miswiring phenotypes, which are reminiscent of subplate ablation, we unbiasedly find a selective loss of SPN gene expression following Arid1a deletion. In addition, multiple characteristics of SPNs crucial to their wiring functions, including subplate organization, subplate-thalamocortical axon co-fasciculation ("handshake"), and extracellular matrix, are severely disrupted. To empirically test Arid1a sufficiency in subplate, we generate a cortical plate deletion of Arid1a that spares SPNs. In this model, subplate Arid1a expression is sufficient for subplate-thalamocortical axon co-fasciculation and extracellular matrix assembly. Consistent with these wiring functions, subplate Arid1a sufficiently enables normal callosum formation, thalamocortical axon targeting, and whisker barrel development. Thus, Arid1a is a multifunctional regulator of subplate-dependent guidance mechanisms essential to cortical circuit wiring.

SignificanceThe cognitive, perceptive, and motor capabilities of the mammalian cerebral cortex depend on assembly of circuit connectivity during development. Subplate neurons, strategically located at the junction of grey and white matter, orchestrate the wiring of cortical circuits. Using a new approach to study gene necessity and sufficiency in subplate neurons, we uncover an essential role for chromatin remodeler Arid1a in subplate neuron gene expression and axon guidance functions. Cortical deletion of Arid1a disrupts subplate-dependent formation of corpus callosum, targeting of thalamocortical axons, and development of sensory maps. Together, our study identifies Arid1a as a central regulator of subplate-dependent axon pathfinding, establishes subplate function as essential to callosum development, and highlights non-cell autonomous mechanisms in neural circuit formation and disorders thereof.


title: Public good exploitation in natural bacterioplankton communities

abstract: Microorganisms such as bacteria often interact with their environment through extracellular molecules that increase access to limiting resources. These secretions can act as public goods, creating incentives for exploiters, a.k.a cheaters, to invade and  steal public goods away from producers. This phenomenon has been studied extensively in microbiology due to its implications for the evolution of cooperation, but little is known about the occurrence and impact of public good exploiters in the environment. Here, we develop a new genomic approach to systematically identify bacteria that can exploit public goods produced during the degradation of polysaccharides. Focusing on chitin - the second most abundant biopolymer on the planet, we show that public good exploiters are active in natural marine microbial communities that assemble on chitin particles, invading during early stages of colonization and potentially hindering degradation. Unlike in classical studies of social evolution, exploiters and polysaccharide degraders are not isogenic and instead belong to distant lineages, facilitating their coexistence. Our approach opens novel avenues to use the wealth of genomic data available to infer ecological roles and interactions among microbes.


title: Recurrent emergence and transmission of a SARS-CoV-2 Spike deletion H69/V70

abstract:
SARS-CoV-2 amino acid replacements in the receptor binding domain (RBD) occur relatively frequently and some have a consequence for immune recognition. Here we report recurrent emergence and significant onward transmission of a six-nucleotide out of frame deletion in the S gene, which results in loss of two amino acids: H69 and V70. We report that in human infections {Delta}H69/V70 often co-occurs with the receptor binding motif amino acid replacements N501Y, N439K and Y453F, and in the latter two cases has followed the RBD mutation. One of the {Delta}H69/V70+ N501Y lineages, now known as B.1.1.7, has undergone rapid expansion and includes eight S gene mutations: RBD (N501Y and A570D), S1 ({Delta}H69/V70 and {Delta}144) and S2 (P681H, T716I, S982A and D1118H). In vitro, we show that {Delta}H69/V70 does not reduce serum neutralisation across multiple convalescent sera. However, {Delta}H69/V70 increases infectivity and is associated with increased incorporation of cleaved spike into virions. {Delta}H69/V70 is able to compensate for small infectivity defects induced by RBD mutations N501Y, N439K and Y453F. In addition, replacement of H69 and V70 residues in the B.1.1.7 spike reduces its infectivity and spike mediated cell-cell fusion. Based on our data {Delta}H69/V70 likely acts as a permissive mutation that allows acquisition of otherwise deleterious immune escape mutations. Enhanced surveillance for the {Delta}H69/V70 deletion with and without RBD mutations should be considered as a global priority not only as a marker for the B.1.1.7 variant, but potentially also for other emerging variants of concern. Vaccines designed to target the deleted spike protein could mitigate against its emergence as increased selective forces from immunity and vaccines increase globally.

HighlightsO_LI{Delta}H69/V70 is present in at least 28 SARS-CoV-2 lineages
C_LIO_LI{Delta}H69/V70 does not confer escape from convalescent sera
C_LIO_LI{Delta}H69/V70 increases spike infectivity and compensates for RBD mutations
C_LIO_LI{Delta}H69/V70 is associated with greater spike cleavage
C_LIO_LIB.1.1.7 requires {Delta}H69/V70 for optimal spike cleavage and infectivity
C_LI


title: Massively parallel transposon mutagenesis identifies temporally essential genes for biofilm formation in Escherichia coli

abstract: Biofilms complete a life cycle where cells aggregate, grow and produce a structured community before dispersing to seed biofilms in new environments. Progression through this life cycle requires temporally controlled gene expression to maximise fitness at each stage. Previous studies have largely focused on the essential genome for the formation of a mature biofilm, but here we present an insight into the genes involved at different stages of biofilm formation. We used TraDIS-Xpress; a massively parallel transposon mutagenesis approach using transposon-located promoters to assay the impact of disruption or altered expression of all genes in the genome on biofilm formation. We determined temporal differences in the importance of genes in E. coli growing as a biofilm on glass beads after 12, 24 and 48 hours. A selection of genes identified as important were then validated independently by assaying biofilm biomass, aggregation, curli production and adhesion ability of defined mutants. We identified 48 genes that affected biofilm fitness including genes with known roles and those not previously implicated in biofilm formation. Regulation of type 1 fimbriae and motility were important at all time points. Adhesion and motility were important for the early biofilm, whereas matrix production and purine biosynthesis were only important as the biofilm matured. We found strong temporal contributions to biofilm fitness for some genes including some where expression changed between being beneficial or detrimental depending on the stage at which they are expressed, including dksA and dsbA. Novel genes implicated in biofilm formation included zapE and truA involved in cell division, maoP in DNA housekeeping and yigZ and ykgJ of unknown function. This work provides new insights into the requirements for successful biofilm formation through the biofilm life cycle and demonstrates the importance of understanding expression and fitness through time.


title: Decomposing the sources of SARS-CoV-2 fitness variation in the United States

abstract: The fitness of a pathogen is a composite phenotype determined by many different factors influencing growth rates both within and between hosts. Determining what factors shape fitness at the host population-level is especially challenging because both intrinsic factors like pathogen genetics and extrinsic factors such as host behaviour influence between-host transmission potential. These challenges have been highlighted by controversy surrounding the population-level fitness effects of mutations in the SARS-CoV-2 genome and their relative importance when compared against non-genetic factors shaping transmission dynamics. Building upon phylodynamic birth-death models, we develop a new framework to learn how hundreds of genetic and non-genetic factors have shaped the fitness of SARS-CoV-2. We estimate the fitness effects of all amino acid variants and several structural variants that have circulated in the United States between February 2020 and March 2021 from viral phylogenies. We also estimate how much fitness variation among pathogen lineages is attributable to genetic versus non-genetic factors such as spatial heterogeneity in transmission rates. Before September 2020, most fitness variation between lineages can be explained by background spatial heterogeneity in transmission rates across geographic regions. Starting in late 2020, genetic variation in fitness increased dramatically with the emergence of several new lineages including B.1.1.7, B.1.427, B.1.429 and B.1.526. Our analysis also indicates that genetic variants in less well-explored genomic regions outside of Spike may be contributing significantly to overall fitness variation in the viral population.


title: Genomic surveillance framework and global population structure for Klebsiella pneumoniae

abstract: K. pneumoniae is a leading cause of antimicrobial-resistant (AMR) healthcare-associated infections, neonatal sepsis and community-acquired liver abscess, and is associated with chronic intestinal diseases. Its diversity and complex population structure pose challenges for analysis and interpretation of K. pneumoniae genome data. Here we introduce Kleborate, a tool for analysing genomes of K. pneumoniae and its associated species complex, which consolidates interrogation of key features of proven clinical importance. Kleborate provides a framework to support genomic surveillance and epidemiology in research, clinical and public health settings. To demonstrate its utility we apply Kleborate to analyse publicly available Klebsiella genomes, including clinical isolates from a pan-European study of carbapenemase-producing Klebsiella, highlighting global trends in AMR and virulence as examples of what could be achieved by applying this genomic framework within more systematic genomic surveillance efforts. We also demonstrate the application of Kleborate to detect and type K. pneumoniae from gut metagenomes.


title: The SEQC2 Epigenomics Quality Control (EpiQC) Study: Comprehensive Characterization of Epigenetic Methods, Reproducibility, and Quantification

abstract: Cytosine modifications in DNA such as 5-methylcytosine (5mC) underlie a broad range of developmental processes, maintain cellular lineage specification, and can define or stratify cancer and other diseases. However, the wide variety of approaches available to interrogate these modifications has created a need for harmonized materials, methods, and rigorous benchmarking to improve genome-wide methylome sequencing applications in clinical and basic research. Here, we present a multi-platform assessment and a global resource for epigenetics research from the FDAs Epigenomics Quality Control (EpiQC) Group. The study design leverages seven human cell lines that are designated as reference materials and publicly available from the National Institute of Standards and Technology (NIST) and Genome in a Bottle (GIAB) consortium. These samples were subject to a variety of genome-wide methylation interrogation approaches across six independent laboratories, with a primary focus was on 5-methylcytosine modifications. Each sample was processed in two or more technical replicates by three whole-genome bisulfite sequencing (WGBS) protocols (TruSeq DNA methylation, Accel-NGS MethylSeq, and SPLAT), oxidative bisulfite sequencing (TrueMethyl), one enzymatic deamination method (EMseq), targeted methylation sequencing (Illumina Methyl Capture EPIC), and single-molecule long-read nanopore sequencing from Oxford Nanopore Technologies. After rigorous quality assessment and comparison to Illumina EPIC methylation microarrays and testing on a range of algorithms (Bismark, BitmapperBS, BWAMeth, and GemBS), we found overall high concordance between assays (R=0.87-R0.93), differences in efficency of read mapping and CpG capture and coverage, and platform performance. The data provided herein can guide continued used of these reference materials in epigenomics assays, as well as provide best practices for epigenomics research and experimental design in future studies.


title: Multiplex genotyping method to validate the multiallelic genome editing outcomes using machine learning-assisted long-read sequencing

abstract: Genome editing can introduce designed mutations into a target genomic site. Recent research has revealed that it can also induce various unintended events such as structural variations, small indels, and substitutions at, and in some cases, away from the target site. These rearrangements may result in confounding phenotypes in biomedical research samples and cause a concern in clinical or agricultural applications. However, current genotyping methods do not allow a comprehensive analysis of diverse mutations for phasing and mosaic variant detection. Here, we developed a genotyping method with an on-target site analysis software named Determine Allele mutations and Judge Intended genotype by Nanopore sequencer (DAJIN) that can automatically identify and classify both intended and unintended diverse mutations, including point mutations, deletions, inversions, and cis double knock-in at single-nucleotide resolution. Our approach with DAJIN can handle approximately 100 samples under different editing conditions in a single run. With its high versatility, scalability, and convenience, DAJIN-assisted multiplex genotyping may become a new standard for validating genome editing outcomes.


title: Fetal estrogens are not involved in sex determination but critical for early ovarian differentiation.

abstract: AROMATASE, encoded by the CYP19A1 gene, is the cytochrome enzyme responsible for the synthesis of estrogens in vertebrates. In most mammals a peak of expression of the CYP19A1 gene occurs in the fetal XX gonad when sexual differentiation starts up. To elucidate the role of this peak, we produced 3 lines of TALEN genetically edited CYP19A1 KO rabbits, that were void of any production of estradiol. All KO XX rabbits developed as females, with aberrantly small sized ovaries at adulthood, an almost empty reserve of primordial follicles and very few large antrum follicles. Ovulation never occurred. Our histological, immunohistological and transcriptomic analyses showed that the surge of estradiol in the XX fetal rabbit gonad is dispensable for its determination as an ovary, or for meiosis. However, it is mandatory for the high proliferation and differentiation of both somatic and germ cells, and consequently for the establishment of the ovarian reserve.

category: developmental biology


title: Paracrine behaviors arbitrate parasite-like interactions between tumor subclones

abstract: Explaining the emergence and maintenance of intratumor heterogeneity is an important question in cancer biology. Tumor cells can generate considerable subclonal diversity, which influences tumor growth rate, treatment resistance, and metastasis, yet we know remarkably little about how cells from different subclones interact. Here, we confronted two murine mammary cancer cell lines to determine both the nature and mechanisms of subclonal cellular interactions in vitro. Surprisingly, we found that, compared to monoculture, growth of the  winner was enhanced by the presence of the  loser cell line, whereas growth of the latter was reduced. Mathematical modeling and laboratory assays indicated that these interactions are mediated by the production of paracrine metabolites resulting in the winner subclone effectively  farming the loser. Our findings add a new level of complexity to the mechanisms underlying subclonal growth dynamics.