title: Bacterial microcompartments for isethionate desulfonation in the taurine-degrading human-gut bacterium Bilophila wadsworthia

abstract:
BackgroundBilophila wadsworthia, a strictly anaerobic, sulfite-reducing bacterium and common member of the human gut microbiota, has been associated with diseases such as appendicitis and colitis. It is specialized on organosulfonate respiration for energy conservation, i.e., utilization of dietary and host-derived organosulfonates, such as taurine (2-aminoethansulfonate), as sulfite donors for sulfite respiration, producing hydrogen sulfide (H2S), an important intestinal metabolite that may have beneficial as well as detrimental effects on the colonic environment. Its taurine desulfonation pathway involves a glycyl radical enzyme (GRE), isethionate sulfite-lyase (IslAB), which cleaves isethionate (2-hydroxyethane sulfonate) into acetaldehyde and sulfite.

ResultsWe demonstrate that taurine metabolism in B. wadsworthia 3.1.6 involves bacterial microcompartments (BMCs). First, we confirmed taurine-inducible production of BMCs by proteomic, transcriptomic and ultra-thin sectioning and electron-microscopical analyses. Then, we isolated BMCs from taurine-grown cells by density-gradient ultracentrifugation and analyzed their composition by proteomics as well as by enzyme assays, which suggested that the GRE IslAB and acetaldehyde dehydrogenase are located inside of the BMCs. Finally, we are discussing the recycling of cofactors in the IslAB-BMCs and a potential shuttling of electrons across the BMC shell by a potential ironsulfur (FeS) cluster-containing shell protein identified by sequence analysis.

ConclusionsWe characterized a novel subclass of BMCs and broadened the spectrum of reactions known to take place enclosed in BMCs, which is of biotechnological interest. We also provided more details on the energy metabolism of the opportunistic pathobiont B. wadsworthia and on microbial H2S production in the human gut.


title: Gene-Specific Nonsense-Mediated mRNA Decay Targeting for Cystic Fibrosis Therapy

abstract: Low CFTR mRNA expression due to nonsense-mediated mRNA decay (NMD) is a major hurdle in developing a therapy for cystic fibrosis (CF) caused by the W1282X mutation in the CFTR gene. CFTR-W1282X truncated protein retains partial function, so increasing its levels by inhibiting NMD of its mRNA will likely be beneficial. Because NMD regulates the normal expression of many genes, gene-specific stabilization of CFTR-W1282X mRNA expression is more desirable than general NMD inhibition. Synthetic antisense oligonucleotides (ASOs) designed to prevent binding of exon junction complexes (EJC) downstream of premature termination codons (PTCs) attenuate NMD in a gene-specific manner. We developed a cocktail of three ASOs that specifically increases the expression of CFTR W1282X mRNA and CFTR protein in ASO-transfected human bronchial epithelial cells. This treatment increased the CFTR-mediated chloride current. These results set the stage for clinical development of an allele-specific therapy for CF caused by the W1282X mutation.

category: molecular biology


title: The Bardet-Biedl protein Bbs1 controls photoreceptor outer segment protein and lipid composition

abstract: Primary cilia are key sensory organelles whose dysfunction leads to ciliopathy disorders such as Bardet-Biedl syndrome (BBS). Retinal degeneration is common in ciliopathies, since the outer segments (OSs) of photoreceptors are highly specialized primary cilia. BBS1, encoded by the most commonly mutated BBS-associated gene, is part of the BBSome protein complex. Using a new bbs1 zebrafish mutant, we show that retinal development and photoreceptor differentiation are unaffected by Bbs1-loss, supported by an initially unaffected transcriptome. Quantitative proteomics and lipidomics on isolated OSs show that Bbs1 is required for BBSome-entry into OSs and that Bbs1-loss leads to accumulation of membrane-associated proteins in OSs, with enrichment in proteins involved in lipid homeostasis. Disruption of the tightly regulated OS lipid composition with increased OS cholesterol content are paralleled by early functional visual deficits, which precede progressive OS morphological anomalies. Our findings identify a new role for Bbs1/BBSome in OS lipid homeostasis and suggest a new pathomechanism underlying retinal degeneration in BBS.


title: Host-symbiont population genomics provide insights into partner fidelity, transmission mode and habitat adaptation in deep-sea hydrothermal vent snails

abstract:
Symbiont specificity, both at the phylotype and strain level, can have profound consequences for host ecology and evolution. However, except for insights from a few model symbiosis systems, the degree of partner fidelity and the influence of host versus environmental factors on symbiont composition are still poorly understood. Nutritional symbioses between invertebrate animals and chemosynthetic bacteria at deep-sea hydrothermal vents are examples of relatively selective associations, where hosts affiliate only with particular, environmentally acquired phylotypes of gammaproteobacterial or campylobacterial symbionts. In hydrothermal vent snails of the sister genera Alviniconcha and Ifremeria this phylotype specificity has been shown to play a role in habitat distribution and partitioning among different holobiont species. However, it is currently unknown if fidelity goes beyond species level associations that might influence genetic structuring, connectivity and habitat adaptation of holobiont populations. We used metagenomic analyses to assess sequence variation in hosts and symbionts and identify correlations with geographic and environmental factors. Our analyses indicate that host populations are not differentiated across a ~800 km gradient, while symbiont populations are clearly structured between vent locations due to a combination of neutral and selective processes. Overall, these results suggest that host individuals flexibly associate with locally adapted strains of their specific symbiont phylotypes, which supports a long-standing but untested paradigm of the benefits of horizontal transmission. Strain flexibility in these snails likely enables host populations to exploit a range of habitat conditions, which might favor wide-spread genetic connectivity and ecological resilience unless physical dispersal barriers are present.

Significance StatementSymbiont composition in horizontally transmitted symbioses is influenced by a combination of host genetics, environmental conditions and geographic barriers. Yet the relative importance of these factors and the effects of adaptive versus neutral evolutionary forces on symbiont population structure remain unknown in the majority of marine symbioses. To address these questions, we applied population genomic approaches in four species of deep-sea hydrothermal vent snails that live in obligate association with chemosynthetic bacteria. Our analyses show that host genetics plays a minor role compared to environment for symbiont strain composition despite specificity to symbiont species and corroborate a long-standing hypothesis that vent invertebrates affiliate with locally adapted symbiont strains to cope with the variable habitat conditions characterizing hydrothermal vents.


category: evolutionary biology


title: Co-translational assembly counteracts promiscuous interactions

abstract: During the co-translational assembly of protein complexes, a fully synthesized subunit engages with the nascent chain of a newly synthesized interaction partner. Such events are thought to contribute to productive assembly, but their exact physiological relevance remains underexplored. Here, we examined structural motifs contained in nucleoporins for their potential to facilitate co-translational assembly. We experimentally tested candidate structural motifs and identified several previously unknown co-translational interactions. We demonstrate by selective ribosome profiling that domain invasion motifs of beta-propellers, coiled-coils, and short linear motifs act as co-translational assembly domains. Such motifs are often contained in proteins that are members of multiple complexes (moonlighters) and engage with closely related paralogs. Surprisingly, moonlighters and paralogs assembled co-translationally in only one but not all of the relevant assembly pathways. Our results highlight the regulatory complexity of assembly pathways.


title: Myonuclear maturation dynamics in aged and adult regenerating mouse skeletal muscle

abstract: Skeletal muscle cells are multinucleated syncytial cells arising from cell fusion, yet despite sharing a common cytoplasm individual myonuclei express distinct transcriptional programs. Whether individual myonuclei acquire heterogenous transcriptional states via differences in their progenitors, during differentiation, or once their anatomical position is acquired, is not known. We performed transcriptome and pseudotime analysis of single myogenic nuclei from uninjured and post-injury murine skeletal muscle to assess when myonuclear heterogeneity is acquired. Two distinct progenitors contribute to myonuclei, one a non-myogenic fibroblast subtype, and skeletal muscle stem cells the other. Both progenitors enter a single pseudotime trajectory that bifurcates as myonuclei mature into two branches segregated by myosin isoform expression and metabolic profiles, suggesting transcriptional heterogeneity is acquired as myonuclei mature. In aged skeletal muscle myogenic progenitor expansion is perturbed and nuclei from aged muscle display distinct pseudotemporal kinetics compared to nuclei from young mice. In aged mice, the inferred myogenic differentiation trajectory is delayed, altering the distribution of myogenic nuclei in pseudotime, suggesting that altered transcriptional dynamics in nuclei in aged mice may drive age-associated muscle deficits and bias myonuclei towards acquiring oxidative metabolic profiles.


title: Plexin-A2 Promotes The Proliferation And The Development Of Tumors From Glioblastoma Derived Cells

abstract: The semaphorin guidance factors receptor Plexin-A2 transduces sema6A and sema6B signals and when associated with neuropilins can also transduce sema3C signals. Inhibition of plexin-A2 expression in U87MG glioblastoma cells resulted in strong inhibition of cell proliferation and tumor forming ability. Knock-out of the plexin-A2 gene using CRISPR/Cas9 also inhibited cell proliferation which was rescued following re-expression of the plexin-A2 cDNA or expression of a truncated plexin-A2 lacking the extracellular domain. Inhibition of plexin-A2 expression resulted in cell cycle arrest at the G2/M stage, and was accompanied by changes in cytoskeletal organization, cell flattening, and by the expression of senescence associated {beta}-galactosidase. It was also associated with reduced AKT phosphorylation and enhanced phosphorylation of p38MAPK. We find that the pro-proliferative effects of plexin-A2 are mediated by FARP2 and FYN since mutations in the FARP2 binding domain of plexin-A2 or in the FYN phosphorylation sites of plexin-A2 compromised the rescue of the proliferative activity by the plexin-A2 intracellular domain. Our results suggest that plexin-A2 may represent a novel target for the development of anti-tumorigenic therapeutics.

category: cell biology


title: A Floxed exon (Flexon) approach to Cre-mediated conditional gene expression

abstract:
Conditional gene expression allows for genes to be manipulated and lineages to be marked during development. In the established "lox-stop-lox" approach, Cre-mediated tissue-specific gene expression is achieved by excising the stop cassette, a lox-flanked translational stop that is inserted into the 5' untranslated region of a gene to halt its expression. Although lox-stop-lox has been successfully used in many experimental systems, the design of traditional stop cassettes also has common issues and limitations. Here, we describe the Floxed exon (Flexon), a stop cassette within an artificial exon that can be inserted flexibly into the coding region of any gene to cause premature termination of translation and nonsense-mediated decay of the mRNA. We demonstrate its efficacy in C. elegans by showing that, when promoters that cause weak and/or transient cell-specific expression are used to drive Cre in combination with a gfp(flexon) transgene, strong and sustained expression is obtained in specific lineages. We also describe several potential additional applications for using Flexon for developmental studies, including more precise control of gene expression using intersectional methods, tissue-specific protein degradation or RNAi, and generation of genetic mosaics. The Flexon approach should be feasible in any system where any site-specific recombination-based method may be applied.

Summary statementThe Floxed exon (Flexon), a stop cassette that can be inserted flexibly into the coding region of any gene, facilitates Cre-mediated conditional gene expression.

category: genetics


title: Optimized Cas9 expression improves performance of large-scale CRISPR screening

abstract: CRISPR technology is an invaluable tool for large-scale functional genomic screening. Genome editing efficiency and timing are important parameters impacting the performance of pooled CRISPR screens. Here we show that by optimizing Cas9 expression levels, the time necessary for gene editing can be reduced contributing to improved performance of CRISPR based screening.

category: genomics


title: A Plasmid System with Tunable Copy Number

abstract: Plasmids are one of the most commonly used and time-tested molecular biology platforms for genetic engineering and recombinant gene expression in bacteria. Despite their ubiquity, little consideration is given to metabolic effects and fitness costs of plasmid copy numbers on engineered genetic systems. Here, we introduce two systems that allow for the finely-tuned control of plasmid copy number: a plasmid with an anhydrotetracycline-controlled copy number, and a massively parallel assay that is used to generate a continuous spectrum of ColE1-based copy number variants. Using these systems, we investigate the effects of plasmid copy number on cellular growth rates, gene expression, biosynthesis, and genetic circuit performance. We perform single-cell timelapse measurements to characterize plasmid loss, runaway plasmid replication, and quantify the impact of plasmid copy number on the variability of gene expression. Using our massively parallel assay, we find that each plasmid imposes a 0.063% linear metabolic burden on their hosts, hinting at a simple relationship between metabolic burdens and plasmid DNA synthesis. Our plasmid system with tunable copy number should allow for a precise control of gene expression and highlight the importance of tuning plasmid copy number as tool for the optimization of synthetic biological systems.

category: synthetic biology


title: Glutamatergic and GABAergic neurons mediate distinct neurodevelopmental phenotypes of STXBP1 encephalopathy

abstract: Heterozygous pathogenic variants in syntaxin-binding protein 1 (STXBP1, also known as MUNC18-1) cause STXBP1 encephalopathy and are among the most frequent causes of developmental and epileptic encephalopathies and intellectual disabilities. STXBP1 is an essential protein for presynaptic neurotransmitter release, and its haploinsufficiency impairs glutamatergic and GABAergic neurotransmission. However, the mechanism underlying the broad spectrum of neurological phenotypes is poorly understood. Here we show that glutamatergic and GABAergic neurons mediate distinct disease features with few overlaps. Glutamatergic and GABAergic neurons-specific Stxbp1 haploinsufficient mice exhibit different subsets of the cognitive and seizure phenotypes observed in the constitutive Stxbp1 haploinsufficient mice. Developmental delay and most of the motor and psychiatric phenotypes are only recapitulated by GABAergic Stxbp1 haploinsufficiency. Thus, the contrasting roles of excitatory and inhibitory signaling in STXBP1 encephalopathy identify GABAergic dysfunction as a main disease mechanism and reveal the possibility to selectively modulate disease phenotypes by targeting specific neurotransmitter systems.

category: neuroscience


title: Individual believes about temporal continuity explain perceptual biases

abstract: Perception of magnitudes such as duration or distance is often found to be systematically biased. The biases, which result from incorporating prior knowledge in the perceptual process, can vary considerably between individuals. The variations are commonly attributed to differences in sensory precision and reliance on priors. However, another factor not considered so far is the implicit belief about how successive sensory stimuli are generated: independently from each other or with certain temporal continuity. The main types of explanatory models proposed so far - static or iterative - mirror this distinction but cannot adequately explain individual biases. Here we propose a new unifying model that explains individual variation as combination of sensory precision and beliefs about temporal continuity and predicts the experimentally found changes in biases when altering temporal continuity. Thus, according to the model, individual differences in perception depend on beliefs about how stimuli are generated in the world.


title: Individual beliefs about temporal continuity explain variation of perceptual biases


title: Chromosome-level genome assembly of a human fungal pathogen reveals synteny among geographically distinct species

abstract:
Histoplasma capsulatum, a dimorphic fungal pathogen, is the most common cause of fungal respiratory infections in immunocompetent hosts. Histoplasma is endemic in the Ohio and Mississippi River Valleys in the United States and also distributed worldwide. Previous studies revealed at least eight clades, each specific to a geographic location: North American classes 1 and 2 (NAm 1 and NAm 2), Latin American groups A and B (LAm A and LAm B), Eurasian, Netherlands, Australian and African, and an additional distinct lineage (H81) comprised of Panamanian isolates. Previously assembled Histoplasma genomes are highly fragmented, with the highly repetitive G217B (NAm 2) strain, which has been used for most whole genome-scale transcriptome studies, assembled into over 250 contigs. In this study, we set out to fully assemble the repeat regions and characterize the large-scale genome architecture of Histoplasma species. We re-sequenced five Histoplasma strains (WU24 (NAm 1), G217B (NAm 2), H88 (African), G186AR (Panama), and G184AR (Panama)) using Oxford Nanopore Technologies long-read sequencing technology. Here we report chromosomal-level assemblies for all five strains, which exhibit extensive synteny among the geographically distant Histoplasma isolates. The new assemblies revealed that RYP2, a major regulator of morphology and virulence, is duplicated in G186AR. In addition, we mapped previously generated transcriptome datasets onto the newly assembled chromosomes. Our analyses revealed that the expression of transposons and transposon-embedded genes are upregulated in yeast phase compared to mycelial phase in G217B and H88 strains. This study provides an important resource for fungal researchers and further highlights the importance of chromosomal-level assemblies in analyzing high-throughput datasets.

ImportanceHistoplasma species are dimorphic fungi causing significant morbidity and mortality worldwide. These fungi grow as mold in the soil and as budding yeast within the human host. Histoplasma can be isolated from soil in diverse regions, including North America, South America, Africa and Europe. Phylogenetically distinct species of Histoplasma have been isolated and sequenced. However, for the commonly used strains, genome assemblies have been fragmented, leading to underutilization of genome-scale data. This study provides chromosome-level assemblies of the commonly used Histoplasma strains using long-read sequencing technology. Comparative analysis of these genomes shows largely conserved gene order within the chromosomes. Mapping existing transcriptome data on these new assemblies reveals clustering of transcriptionally co-regulated genes. Results of this study highlight the importance of obtaining chromosome-level assemblies in understanding the biology of human fungal pathogens.

category: microbiology


title: Life history traits and dispersal shape neutral genetic diversity in metapopulations

abstract: Genetic diversity at population scale, depends on species life-history traits, population dynamics and local and global environmental factors. We first investigate the effect of life-history traits on the neutral genetic diversity of a single population using a deterministic mathematical model. When the population is stable, we show that semelparous species with precocious maturation and iteroparous species with delayed maturation exhibit higher diversity because their life history traits tend to balance the lifetimes of non reproductive individuals (juveniles) and adults which reproduce. Then, we extend our model to a metapopulation to investigate the additional effect of dispersal on diversity. We show that dispersal may truly modify the local effect of life history on diversity. As a result, the diversity at the global scale of the metapopulation differ from the local diversity which is only described through local life history traits of the populations. In particular, dispersal usually promotes diversity at the global metapopulation scale.

category: ecology


title: Experimental evidence that metapopulation structure can accelerate adaptive evolution

abstract: Whether the spatial arrangement of a population influences adaptive evolution has been a long-standing question in population genetics. In contrast to standard population genetic models, evolutionary graph theory (EGT) predicts certain topologies amplify (increase) the probability that a beneficial mutation will spread in the population relative to a well-mixed population. Here, we test these predictions empirically by tracking the fixation dynamics of an antibiotic resistant mutant under positive selection as it spreads through networks of different topologies both in vitro and in silico. We show that star-like topologies involving bi-directional dispersal between a central hub and peripheral leaves can be amplifiers of selection relative to a well-mixed network, consistent with the predictions of EGT. We further show that the mechanism responsible for amplification is the reduced probability that a rare beneficial mutant will be lost due to drift when it encounters a new patch. Our results provide the first empirical support for the prediction of EGT that spatial structure can amplify the spread of a beneficial mutation and broadens the conditions under which this phenomenon is thought to occur. We also show the importance of considering the migration rate, which is not independently adjustable in most previous models. More generally, our work underscores the potential importance of spatial structure in governing adaptive evolution by showing how the interplay between spatial structure and evolutionary forces determine the fate of a beneficial mutation. It also points the way towards using network topology to amplify the effects of weakly favoured mutations under directed evolution in industrial applications.


title: Climate and mountains shaped human ancestral genetic lineages

abstract: Extensive sequencing of modern and ancient human genomes has revealed that contemporary populations can be explained as the result of recent mixing of a few distinct ancestral genetic lineages1. But the small number of aDNA samples that predate the Last Glacial Maximum means that the origins of these lineages are not well understood. Here, we circumvent the limited sampling by modelling explicitly the effect of climatic changes and terrain on population demography and migrations through time and space, and show that these factors are sufficient to explain the divergence among ancestral lineages. Our reconstructions show that the sharp separation between African and Eurasian lineages is a consequence of only a few limited periods of connectivity through the arid Arabian peninsula, which acted as the gate out of the Arican continent. The subsequent spread across Eurasia was then mostly shaped by mountain ranges, and to a lesser extent deserts, leading to the split of European and Asians, and the further diversification of these two groups. A high tolerance to cold climates allowed the persistence at high latitudes even during the Last Glacial Maximum, maintaining a pocket in Beringia that led to the later, rapid colonisation of the Americas. The advent of food production was associated with an increase in movement2, but mountains and climate have been shown to still play a major role even in this latter period3,4, affecting the mixing of the ancestral lineages that we have shown to be shaped by those two factors in the first place.


title: Divergence in dimerization and activity of primate APOBEC3C

abstract: The APOBEC3 (A3) family of single-stranded DNA cytidine deaminases are host restriction factors that inhibit lentiviruses, such as HIV-1, in the absence of the Vif protein that causes their degradation. Deamination of cytidine in HIV-1 (-)DNA forms uracil that causes inactivating mutations when uracil is used as a template for (+)DNA synthesis. For APOBEC3C (A3C), the chimpanzee and gorilla orthologues are more active than human A3C, and the Old World Monkey A3C from rhesus macaque (rh) is not active against HIV-1. Multiple integrated analyses determined why rhA3C was not active against HIV-1 and how to increase this activity. Biochemical, virological, and coevolutionary analyses combined with molecular dynamics simulations showed that the key amino acids needed to promote rhA3C antiviral activity also promoted dimerization. Although rhA3C shares a similar dimer interface with hominid A3C, the key amino acid contacts were different. Overall, our results determine the basis for why rhA3C is less active than human A3C, establish the amino acid network for dimerization and increased activity, and track the loss and gain of A3C antiviral activity in primates. The coevolutionary analysis of the A3C dimerization interface provides a basis from which to analyze dimerization interfaces of other A3 family members.

category: biochemistry


title: Structural insight into the dual function of LbpB in mediating Neisserial pathogenesis

abstract: Lactoferrin binding protein B (LbpB) is a lipoprotein present on the surface of Neisseria that has been postulated to serve dual functions during pathogenesis in both iron acquisition from lactoferrin, and in providing protection against the cationic antimicrobial peptide lactoferricin. Here, we present the structures of LbpB from N. meningitidis and N. gonorrhoeae in complex with human holo-lactoferrin, forming a 1:1 complex and confirmed by SEC-SAXS. LbpB consists of N- and C-lobes with the N-lobe interacting extensively with the C-lobe of lactoferrin. Our structures provides insight into LbpBs preference towards holo-lactoferrin, and our mutagenesis and binding studies show that lactoferrin and lactoferricin bind independently. Our studies provide the molecular details for how LbpB serves to capture and preserve lactoferrin in an iron-bound state for delivery to the membrane transporter LbpA for iron piracy, and as an antimicrobial peptide sink to evade host immune defenses.

category: biophysics


title: PDF Data Extractor (PDE) - A Free Web Application and R Package Allowing the Extraction of Tables from Portable Document Format (PDF) Files and High-Throughput Keyword Searches of Full-Text Articles

abstract:
The PDF Data Extractor (PDE) R package is designed to perform comprehensive literature reviews for scientists at any stage in a user-friendly way. The PDE_analyzer_i() function permits the user to filter and search thousands of scientific articles using a simple user interface, requiring no bioinformatics skills. In the additional PDE_reader_i() interface, the user can then quickly browse the sentences with detected keywords, open the full-text article, when required, and convert tables conveniently from PDF files to Excel sheets (pdf2table). Specific features of the literature analysis include the adaptability of analysis parameters and the detection of abbreviations of search words in articles.

In this article, we demonstrate and exemplify how the PDE package allows the user-friendly, efficient, and automated extraction of meta-data from full-text articles, which can aid in summarizing the existing literature on any topic of interest. As such, we recommend the use of the PDE package as the first step in conducting an extensive review of the scientific literature. The PDE package is available from the Comprehensive R Archive Network at https://CRAN.R-project.org/package=PDE.

category: scientific communication and education


title: Searching for Consistent Brain Network Topologies Across the Garden of (Shortest) Forking Paths

abstract: The functional interactions between regions of the human brain can be viewed as a network, empowering neuroscientists to leverage tools such as graph theory to obtain insight about brain function. However, obtaining a brain network from functional neuroimaging data inevitably involves multiple steps of data manipulation, which can affect the organisation (topology) of the resulting network and its properties. Test-retest reliability is a gold standard for both basic research and clinical use: a suitable data-processing pipeline for brain networks should recover the same network topology across repeated scan sessions of the same individual. Analyzing resting-state functional Magnetic Resonance Imaging (rs-fMRI) recordings from two test-retest studies across short (45 minutes), medium (2-4 weeks) and long term delays (5-16 months), we investigated the reliability of network topologies constructed by applying 576 unique pipelines to the same fMRI data, obtained from considering combinations of atlas type and size, edge definition and thresholding, and use of global signal regression. We adopted the portrait divergence, an information-theoretic criterion to measure differences in network topology across all scales, enabling us to quantify the influence of different pipelines on the overall organisation of the resulting network. Remarkably, our findings reveal that the choice of pipeline plays a fundamental role in determining how reproducible an individuals brain network topology will be across different scans: there is large and systematic variability across pipelines, such that an inappropriate choice of pipeline can distort the resulting network more than an interval of several months between scans. Across datasets and time-spans, we also identify specific combinations of data-processing steps that consistently yield networks with reproducible topology, enabling us to make recommendations about best practices to ensure high-quality brain networks.


title: Disentangling environmental effects in microbial association networks

abstract:
BackgroundEcological interactions among microorganisms are fundamental for ecosystem function, yet they are mostly unknown or poorly understood. High-throughput-omics can indicate microbial interactions through associations across time and space, which can be represented as association networks. Associations could result from either ecological interactions between microorganisms, or from environmental selection, where the associations are environmentally-driven. Therefore, before downstream analysis and interpretation, we need to distinguish the nature of the association, particularly if it is due to environmental selection or not.

ResultsWe present EnDED (Environmentally-Driven Edge Detection), an implementation of four approaches as well as their combination to predict which links between microorganisms in an association network are environmentally-driven. The four approaches are Sign Pattern, Overlap, Interaction Information, and Data Processing Inequality. We tested EnDED on networks from simulated data of 50 microorganisms. The networks contained on average 50 nodes and 1087 edges, of which 60 were true interactions but 1026 false associations (i.e. environmentally-driven or due to chance). Applying each method individually, we detected a moderate to high number of environmentally-driven edges--87% Sign Pattern and Overlap, 67% Interaction Information, and 44% Data Processing Inequality. Combining these methods in an intersection approach resulted in retaining more interactions, both true and false (32% of environmentally-driven associations). After validation with the simulated datasets, we applied EnDED on a marine microbial network inferred from 10 years of monthly observations of microbial-plankton abundance. The intersection combination predicted that 8.3% of the associations were environmentally-driven, while individual methods predicted 24.8% (Data Processing Inequality), 25.7% (Interaction Information), and up to 84.6% (Sign Pattern as well as Overlap). The fraction of environmentally-driven edges among negative microbial associations in the real network increased rapidly with the number of environmental factors.

ConclusionsTo reach accurate hypotheses about ecological interactions, it is important to determine, quantify, and remove environmentally-driven associations in marine microbial association networks. For that, EnDED offers up to four individual methods as well as their combination. However, especially for the intersection combination, we suggest using EnDED with other strategies to reduce the number of false associations and consequently the number of potential interaction hypotheses.


title: Imaging and single cell sequencing analyses of super-enhancer activation mediated by NF-κB in B cells

abstract: The transcription factor NF-{kappa}B, which plays an important role in cell fate determination, is involved in the activation of super-enhancers (SEs). However, the biological functions of the NF-{kappa}B SEs in gene control are not fully elucidated. We investigated the characteristics of NF-{kappa}B-mediated SE activity using fluorescence live-cell imaging of RelA, single-cell transcriptome, and chromatin accessibility analyses in anti-IgM-stimulated B cells. Cell stimulation induced nuclear foci formation of RelA and gene expression in a switch-like manner. The gained SEs induced a higher fold-change expression and enhanced cell-to-cell variability in transcriptional response. These properties were correlated with the number of gained cis-regulatory interactions, while switch-like gene induction was associated with the number of NF-{kappa}B binding sites in SE. Our study suggests that NF-{kappa}B SEs have unique roles for quantitative control of gene expression through direct binding to accessible DNA and enhanced DNA contacts.


title: Enhanced transcriptional heterogeneity mediated by NF-κB super-enhancers


title: ATP-Independent Nucleosome Unfolding by FACT: Electron Microscopy Analysis

abstract:
FACT is a histone chaperone that unfolds nucleosomes without ATP hydrolysis. We used electron microscopy to study FACT and FACT:nucleosome complexes, and found that both adopt broad ranges of configurations, indicating high flexibility. We found unexpectedly that the DNA binding protein Nhp6 also binds to the C-terminal tails of FACT subunits, inducing more open geometries of FACT even in the absence of nucleosomes. Nhp6 therefore supports nucleosome unfolding by altering both FACT structure and nucleosome properties. Complexes formed with FACT, Nhp6, and nucleosomes also produced a broad range of structures, revealing a large number of potential intermediates along a proposed unfolding pathway. The data suggest that Nhp6 has multiple roles before and during nucleosome unfolding by FACT, and that the process proceeds through a series of energetically similar intermediate structures, ultimately leading to an extensively unfolded form.

One Sentence SummaryElectron microscopy reveals the pathway of ATP-Independent nucleosome unfolding by histone chaperone FACT.


title: Functional peroxisomes are required for heat shock-induced hormesis in Caenorhabditis elegans

abstract: Exact mechanisms of heat shock induced lifespan extension, while documented across species, are still not well understood. Here we put forth evidence that fully functional peroxisomes are required for the activation of the canonical heat shock response and heat-induced hormesis in C. elegans. While during heat shock the HSP-70 chaperone is strongly upregulated in the wild-type (WT) as well as in the absence of peroxisomal catalase ({Delta}ctl-2 mutant), the small heat shock proteins display modestly increased expression in the mutant. Nuclear localization of HSF-1 is reduced in the{Delta} ctl-2 mutant. In addition, heat-induced lifespan extension, observed in the WT, is absent in the{Delta} ctl-2 mutant. Activation of the antioxidant response, the pentose phosphate pathway and increased triglyceride content are the most prominent changes observed during heat shock in the WT worm, but not in the{Delta} ctl-2 mutant. Involvement of peroxisomes in the cell-wide response to transient heat shock reported here gives new insight into the role of organelle communication in the organisms stress response.


title: Morphological integration during postnatal ontogeny: implications for evolutionary biology

abstract: How covariance patterns of phenotypes change during development is fundamental for a broader understanding of evolution. There is compelling evidence that mammalian skull covariance patterns change during ontogeny. However, it is unclear to what extent variation in covariance patterns during ontogeny can impact the response to selection. To tackle this question we explored: i) the extent to which covariance patterns change during postnatal ontogeny; ii) in which ontogenetic stages covariance patterns differ the most, and iii) the extent to which the phenotypic covariance pattern at different ontogenetic stages can be explained by the same processes determining additive genetic covariance. We sampled postnatal ontogenetic series for both marsupials, and placentals. Within each ontogenetic series, we compared covariance matrices (P-matrices) at different ontogenetic stages. Furthermore, we compared these P-matrices to two target matrices [adult P-matrix and an additive genetic covariance matrix (G-matrix)]. Our results show that for all ontogenetic series, covariance patterns from weaning onward are conserved and probably shaped by the same processes determining the G-matrix. We conclude that irrespective of eventual differences in how selection operates during most of postnatal ontogeny, the net response to such pressures will probably not be affected by ontogenetic differences in the covariance pattern.

category: developmental biology


title: Citrus photosynthesis and morphology acclimate to phloem-affecting huanglongbing disease at the leaf and shoot levels

abstract: Huanglongbing (HLB) is a phloem-affecting disease of citrus that reduces growth and has impacted global citrus production. HLB is caused by a phloem-limited bacterium (Candidatus Liberibacter asiaticus; CLas). By inhibiting phloem function, HLB stunts sink growth, including reducing production of new shoots and leaves, and induces hyperaccumulation of foliar starch. HLB induces feedback inhibition of photosynthesis by reducing foliar carbohydrate export. In this work we assessed the relationship of bacterial distribution within the foliage, foliar starch accumulation, and net CO2 assimilation (Anet). Because HLB impacts canopy morphology, we developed a chamber to measure whole-shoot Anet to test the effects of HLB at both leaf and shoot levels. Whole-shoot-level Anet saturated at high irradiance, and green stems had high photosynthetic rates compared to leaves. Starch accumulation was correlated with bacterial population, and starch was negatively correlated with Anet at the leaf level but not at the shoot level. Starch increased initially after infection, then decreased progressively with increasing length of infection. HLB infection reduced Anet at the leaf level, but increased it at the whole shoot level, in association with reduced leaf size and greater relative contribution of stems to photosynthetic surface area. Although HLB-increased photosynthetic efficiency, total carbon fixed per shoot decreased because photosynthetic surface area was reduced. We conclude that the localized effects of infection on photosynthesis are mitigated by whole shoot morphological acclimation over time. Stems contribute important proportions of whole shoot Anet, and these contributions are likely increased by the morphological acclimation induced by HLB.

category: plant biology


title: High-throughput aminoacyl-tRNA synthetase engineering for genetic code expansion in yeast

abstract:
Protein expression with genetically encoded noncanonical amino acids (ncAAs) benefits a broad range of applications, from the discovery of biological therapeutics to fundamental biological studies. A major factor limiting the use of ncAAs is the lack of orthogonal translation systems (OTSs) that support efficient genetic code expansion at repurposed stop codons. Aminoacyl-tRNA synthetases (aaRSs) have been extensively evolved in E. coli but are not always orthogonal in eukaryotes. In this work, we use a yeast display-based ncAA incorporation reporter platform with fluorescence-activated cell sorting to screen libraries of aaRSs in high throughput for 1) incorporation of ncAAs not previously encoded in yeast; 2) improvement of the performance of an existing aaRS; 3) highly selective OTSs capable of discriminating between closely related ncAA analogs; and 4) OTSs exhibiting enhanced polyspecificity to support translation with structurally diverse sets of ncAAs. The number of previously undiscovered aaRS variants we report in this work more than doubles the total number of translationally active aaRSs available for genetic code manipulation in yeast. The success of myriad screening strategies has important implications related to the fundamental properties and evolvability of aaRSs. Furthermore, access to OTSs with diverse activities and specific or polyspecific properties is invaluable for a range of applications within chemical biology, synthetic biology, and protein engineering.

SynopsisA range of flow cytometry-based screens yielded diverse translational machinery for genetic code expansion in yeast, facilitating access to new chemistries and tunable specificity profiles.

O_FIG_DISPLAY_L [Figure 1] M_FIG_DISPLAY C_FIG_DISPLAY


title: Single Serine on TSC2 Exerts Biased Control over mTORC1 Activation by ERK1/2 but Not Akt

abstract: The mammalian target of rapamycin complex 1 (mTORC1) is tightly controlled by tuberous sclerosis complex-2 (TSC2) that is regulated by phosphorylation from kinases responding to environmental cues. Protein kinase G specifically modifies serine-1365 (S1364, human), and its phosphorylation (or phosphomimetic SE mutant) potently blocks mTORC1 co-activation by pathological stress, while a phospho-silenced (SA) mutation does the opposite. Neither alter basal mTORC1 activity. Here we show S1365 exerts biased control over mTORC1 activity (S6K phosphorylation) modifying ERK1/2 but not Akt-dependent stimulation. Whereas mTORC1 activation by endothelin-1 is potently modified by S1365 status, insulin or PDGF stimulation are unaltered. TSC2-S1365 is also phosphorylated upon ET-1 but not insulin stimulation in a PKG-dependent manner, revealing intrinsic bias. Neither energy or nutrient modulation of mTORC1 are impacted by S1365. Consistent with these results, knock-in mice with either TSC2 SA or SE mutations develop identical obesity, glucose intolerance, and fatty liver disease from a high fat diet. Thus, S1365 provides an ERK1/2-selective mTORC1 control mechanism and a genetic means to modify pathological versus physiological mTOR stimuli.


title: Turnover Regulation of the Rho GTPase Cdc42 by Heat Shock Protein Chaperones and the MAPK Pathway Scaffold Bem4

abstract:
All cells maintain an axis of polarity that directs the orientation of growth. Cell polarity can be reorganized during development and in response to extrinsic cues to produce new cell types. Rho GTPases are central regulators of cell polarity and signal-dependent cell differentiation. We show here that one of the best understood Rho GTPases, the highly conserved yeast Cdc42p, is turned over by members of the Heat Shock family of Proteins (HSPs). The Hsp40p chaperone, Ydj1p, was required for turnover of Cdc42p by the NEDD4 E3 ubiquitin ligase, Rsp5p, in the proteosome. Cdc42p turnover was regulated by HSPs at high temperatures, and in aging cells where the protein formed aggregates, implicating HSPs in Rho GTPase quality control. We also show that Cdc42pQ61L, which mimics the active (GTP-bound) conformation of the protein, was turned over at elevated levels by Ydj1p and Rsp5p. A turnover-defective version of Cdc42pQ61L led to multibudding phenotypes, implicating Cdc42 turnover in singularity in cell polarization. Cdc42p turnover also impacted MAP kinase pathway specificity. A pathway-specific scaffold, Bem4p, stabilized Cdc42p levels, which biased Cdc42p function in one MAPK pathway over another. Turnover regulation of Rho GTPases by HSPs and scaffolds provides new dimensions to the regulation of cell polarity and signal-dependent morphogenesis.

Significance StatementRho GTPases are switch-like proteins that govern major decisions in cell polarity and signaling in eukaryotes. We elucidate here a pathway that turns over the yeast Rho GTPase Cdc42p, which is mediated by the heat-shock family of proteins (HSPs) and the NEDD4-type E3 ubiquitin ligase Rsp5p. This finding provides a way for HSPs to exert their widespread effects on morphogenetic responses, phenotypic plasticity, and signaling pathways. We also found that turnover of an active version of Cdc42p is critical for modulating cell polarity. Cdc42p turnover also impacted its function in a pathway specific setting, as stabilization of Cdc42p by Bem4p (SmgGDS-type scaffold) influenced the activity of a specific MAPK pathway. HSPs may regulate Rho GTPase turnover in many systems.


title: Substrate, temperature, and geographical patterns among nearly 2,000 natural yeast isolates

abstract:
Yeasts have broad importance as industrially and clinically relevant microbes and as powerful models for fundamental research, but we are only beginning to understand the roles yeasts play in natural ecosystems. Yeast ecology is often more difficult to study compared to other, more abundant microbes, but growing collections of natural yeast isolates are beginning to shed light on fundamental ecological questions. Here we used environmental sampling and isolation to assemble a dataset of 1,962 isolates collected from throughout the contiguous United States of America (USA) and Alaska, which were then used to uncover geographic patterns, along with substrate and temperature associations among yeast taxa. We found some taxa, including the common yeasts Torulaspora delbrueckii and Saccharomyces paradoxus, to be repeatedly isolated from multiple sampled regions of the US, and we classify these as broadly distributed cosmopolitan yeasts. A number of yeast taxon - substrate associations were identified, some of which were novel and some of which support previously reported associations. Further, we found a strong effect of isolation temperature on the phyla of yeasts recovered, as well as for many species. We speculate that substrate and isolation temperature associations reflect the ecological diversity of and niche partitioning by yeast taxa.

Take AwayO_LIAnalysis of environmental metadata of nearly 2,000 yeast isolates.
C_LIO_LIIndividual yeast taxa associate with specific substrates and plant genera.
C_LIO_LIOptimal yeast isolation temperature differs depending on taxonomic rank.
C_LIO_LISubstrate type and isolation temperatures affect isolated yeast diversity.
C_LI


title: Long-read isoform sequencing reveals tissue-specific isoform expression between active and hibernating brown bears (Ursus arctos)

abstract: Understanding hibernation in brown bears (Ursus arctos) can provide insight into many human diseases. During hibernation, brown bears experience states of insulin resistance, physical inactivity, extreme bradycardia, obesity, and the absence of urine production. These states closely mimic human diseases such as type 2 diabetes, muscle atrophy, renal and heart failure, cachexia, and obesity. The reversibility of these states from hibernation to active season allows for the identification of novel mediators with possible therapeutic value for humans. Recent studies have identified genes and pathways that are differentially expressed between active and hibernation seasons. However, little is known about the role of differential expression of gene isoforms on hibernation physiology. To identify both distinct and novel mRNA isoforms, we performed full-length RNA-sequencing (Iso-Seq) on three tissue types from three individuals sampled during both active and hibernation seasons. We combined the long-read data with the reference annotation for an improved transcriptome and mapped RNA-seq data from six individuals to the improved transcriptome to quantify differential isoform usage between tissues and seasons. We identified differentially expressed isoforms in all study tissues and showed that adipose has a high level of differential isoform usage with isoform switching, regardless of whether the genes were differentially expressed. Our analyses provide a comprehensive evaluation of isoform usage between active and hibernation states, revealing that differential isoform usage, even in the absence of differential gene expression, is an important mechanism for modulating genes during hibernation. These findings demonstrate the value of isoform expression studies and will serve as the basis for deeper exploration into hibernation biology.


title: Whole animal modelling reveals neuronal mechanisms of decision-making and reproduces unpredictable swimming in frog tadpoles

abstract: Animal behaviour is based on interaction between nervous, musculoskeletal and environmental systems. How does an animal process sensory stimuli, use it to decide whether and how to respond, and initiate the locomotor behaviour? We build the whole body computer models of a simple vertebrate with a complete chain of neural circuits and body units for sensory information processing, decision-making, generation of spiking activities, muscle innervation, body flexion, body-water interaction, and movement. Our Central Nervous System (CNS) model generates biologically-realistic spiking and reveals that sensory memory populations on two hindbrain sides compete for swimming initiation and first body flexion. Biomechanical 3-dimensional "Virtual Tadpole" (VT) model is constructed to evaluate if motor outputs of CNS model can produce swimming-like movements in a volume of "water". We find that whole animal modelling generates reliable and realistic swimming. The combination of CNS and VT models opens a new perspective for experiments with immobilised tadpoles.

category: bioinformatics


title: Statistical Transfer Learning with Generative Encoding for Spatial Transcriptomics

abstract: AO_SCPLOWBSTRACTC_SCPLOWA flexible model is introduced which shares ideas with the Autoencoder, Canonical Correlation Analysis, Singular Value Decomposition, and Procrustes Analysis. It is proposed to find relevant maps to transform multiple datasets of various types from one modality to another. Here, the Generative Encoder is used to transform spatial gene expression from breast tissue, to the images of histology tissue measured with Spatial Transcriptomics. The model is directly interpretable given all parameters are linked to the data space. It is scalable on Big Data, training reasonably on several thousand RGB images of 100 by 100 pixels in under an hour, which equates to 30,000 pixel features per sample image.


title: Scalable parametric encoding of multiple modalities


title: Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier

abstract: It is increasingly recognized that sex chromosomes are not only the "battlegrounds" between sexes, but also the "Great Walls" fencing-off introgression between diverging lineages. Here we describe conflicting roles of nascent sex chromosomes on patterns of introgression in an experimental hybrid swarm. Drosophila nasuta and D. albomicans are recently diverged, fully fertile sister species that have different sex chromosome systems. The fusion between an autosome (Muller CD) with the ancestral X and Y gave rise to neo-sex chromosomes in D. albomicans, while Muller CDs remains unfused in D. nasuta. We found that a large block containing overlapping inversions on the neo-sex chromosome stood out as the strongest barrier to introgression. Intriguingly, the neo-sex chromosome introgression barrier is asymmetrical in a sex-dependent manner. Female hybrids showed significant D. albomicans biased introgression on Muller CD (neo-X excess), while males showed heterosis with excessive (neo-X, D. nasuta Muller CD) genotypes. While the neo-Y is a more compatible pairing partner of the neo-X, it also shows moderate levels of degeneration and may thus be selectively disfavored, and sex ratio assay revealed heterospecific meiotic drive. We used a population genetic model to dissect the interplay of sex chromosome drive, heterospecific pairing incompatibility between the neo-sex chromosomes and unfused Muller CD, neo-Y disadvantage, and neo-X advantage in generating the observed neo-X excess in females and heterozygous (neo-X, D. nasuta Muller CD) genotypes in males. We show that moderate neo-Y disadvantage and D. albomicans specific meiotic drive are required to counteract the effect of heterospecific meiotic drive observed in our cross, in concert with pairing incompatibility and neo-X advantage to explain observed genotype frequencies. Together, this hybrid swarm between a young species pair shed light onto the dual roles of neo-sex chromosome evolution in creating a sex-dependent asymmetrical introgression barrier at species boundary.


title: Long-read transcriptome sequencing analysis with IsoTools

abstract: Long-read transcriptome sequencing (LRTS) holds the promise to boost our understanding of alternative splicing. Recent advances in accuracy and throughput have diminished the major limitations and enabled the direct quantification of isoforms. Considering the complexity of the data and the broad range of potential applications, it is clear that highly flexible, accurate analysis tools are crucial. Here, we present IsoTools, a comprehensive Python-based analysis package, for the improvement of alternative and differential splicing analysis. Iso-Tools provides a comprehensive data structure that integrates genomic information from LRTS transcripts together with the reference annotation, and enables broad functionality to quality control, visualize and analyze the data. Additionally, we implemented a graph-based method for the identification of alternative splicing events and a statistical approach based on the beta binomial distribution for the detection of differential events. To demonstrate our methods, we generated PacBio Iso-Seq data of human hepatocytes treated with the HDAC inhibitor valproic acid, a compound known to induce widespread transcriptional changes. Contrasted with short read RNA-Seq of the same samples, this analysis shows that LRTS provides valuable additional insights for a better understanding of alternative splicing, in particular with respect to complex novel and differential splicing events. IsoTools is made available for the community along with extensive documentation at https://github.com/MatthiasLienhard/isotools.


title: Quantifying the size and healing of volumetric muscle wounds using 3D Slicer on CT scans in a canine model

abstract: Volumetric soft tissue and muscle wounds can arise from trauma or necrotizing soft tissue infection. Quantifying the size of these wounds can be challenging, as they often have irregular borders and contours and invariably involve skin loss. 3-dimensional Computed Tomography (3dCT) has been used to characterize the volume of numerous tissue structures, but these use cases invariably involve structures for which clear anatomic borders exist. This is not the case for volumetric soft tissue or muscle wounds, where the volume of the wound being assessed, which is actually a void representing the absence of tissue, does not contain an explicit border at the superficial surface. We present a method that allows quantification of the void size of volumetric muscle wounds using CT scans processed with the software package 3D Slicer. This quantification allows us to chart the progression of healing in such wounds with sequential scans. The development of a means to quantify wound size and healing rate is a necessary capability in order to assess the efficacy of potential therapeutic interventions aimed at enhancing healing of such wounds.

category: physiology


title: Parallel processing by distinct classes of principal neurons in the olfactory cortex

abstract: Understanding the specific roles that different neuron types play within a neural circuit is essential for understanding what that circuit does and how it does it. Primary olfactory (piriform, PCx) cortex contains two main types of principal neurons: semilunar (SL) and pyramidal (PYR). SLs and PYRs have different morphologies, connectivity, biophysical properties, and downstream projections, predicting distinct roles in cortical odor processing. The prevailing model is that odor processing in PCx occurs in two stages, where SLs are the primary recipients of olfactory bulb (OB) input, and PYRs receive and transform information from SLs. Here we recorded from opto-genetically-identified SLs and PYRs in awake, head-fixed mice. We found differences in SL and PYR odor-evoked activity that reflect their different connectivity profiles. But SL responses did not precede PYR responses and suppressing SL activity had little effect on PYR odor responses. These results suggest that SLs and PYRs form parallel odor processing channels.


title: The ternary complex CEP90, FOPNL and OFD1 specifies the future location of centriolar distal appendages, and promotes their assembly

abstract:
Cilia assembly starts with centriole to basal body maturation, migration to the cell surface and docking to the plasma membrane. The basal body docking process involves the interaction of both the distal end of the basal body and the transition fibers (or mature distal appendages), with the plasma membrane. During this process, the transition zone assembles and forms the structural junction between the basal body and the nascent cilium. Mutations in numerous genes involved in basal body docking and transition zone assembly are associated with the most severe ciliopathies, highlighting the importance of these events in cilium biogenesis.

The conservation of this sequence of events across phyla is paralleled by a high conservation of the proteins involved. We identified CEP90 by BioID using FOPNL as a bait. Ultrastructure expansion microscopy showed that CEP90, FOPNL and OFD1 localized at the distal end of both centrioles/basal bodies in Paramecium and mammalian cells. These proteins are recruited early after duplication on the procentriole. Finally, functional analysis performed both in Paramecium and mammalian cells demonstrate the requirement of this complex for distal appendage assembly and basal body docking. Altogether, we propose that this ternary complex is required to determine the future position of distal appendages.


title: IL7R is required for hematopoietic stem cell reconstitution of tissue-resident lymphoid cells

abstract: Traditional, adult-derived lymphocytes that circulate provide adaptive immunity to infection and pathogens. However, subsets of lymphoid cells are also found in non-lymphoid tissues and are called tissue-resident lymphoid cells (TLCs). TLCs encompass a wide array of cell types that span the spectrum of innate-to-adaptive immune function. Unlike traditional lymphocytes that are continuously generated from hematopoietic stem cells (HSCs), many TLCs are of fetal origin and poorly generated from adult HSCs. Here, we sought to understand the development of murine TLCs across multiple tissues and therefore probed the roles of Flk2 and IL7R, two cytokine receptors with known roles in traditional lymphopoiesis. Using Flk2- and Il7r-Cre lineage tracing models, we found that peritoneal B1a cells, splenic marginal zone B (MZB) cells, lung ILC2s and regulatory T cells (Tregs) were highly labeled in both models. Despite this high labeling, highly quantitative, in vivo functional approaches showed that the loss of Flk2 minimally affected the generation of these cells in situ. In contrast, the loss of IL7R, or combined deletion of Flk2 and IL7R, dramatically reduced the cell numbers of B1a cells, MZBs, ILC2s, and Tregs both in situ and upon transplantation, indicating an intrinsic and more essential role for IL7R. Surprisingly, reciprocal transplants of WT HSCs showed that an IL7R-/- environment selectively impaired reconstitution of TLCs when compared to TLC numbers in situ. Taken together, our data revealed functional roles of Flk2 and IL7R in the establishment of tissue-resident lymphoid cells.


title: The roles of NADPH oxidases during adult zebrafish fin regeneration

abstract: Sustained elevated levels of reactive oxygen species (ROS) have been shown to be essential for whole body, appendage and organ regeneration in various organisms, including planarians, Hydra, zebrafish, axolotl, Xenopus, geckos and mice. In the majority of cases these roles have been shown via the use of NADPH oxidase pharmacological inhibitors, which generally target all NAPDH oxidases (NOXes). To identify the specific NOX or NOXes essential for ROS production during adult fin regeneration in zebrafish, we generated nox mutants for duox, nox5 and cyba (a key subunit of NOXes 1-4). We also crossed these mutant lines to a transgenic line ubiquitously expressing HyPer, which permits the measurement of ROS levels in adult zebrafish fins. Using this approach, we found that homozygous duox mutants have significantly attenuated ROS levels following fin amputation, and this correlated with a significantly diminished rate of fin regeneration. While the other nox homozygous mutants (nox5 and cyba) showed less of an effect on ROS levels or adult fin regeneration, duox/cyba double mutants showed a more diminished rate of fin regeneration than duox mutants alone, suggesting that Nox1-4 do play a role during regeneration, but one that is secondary to that of Duox. This work also serendipitously found that ROS levels in amputated adult zebrafish fins oscillate during the day with a circadian rhythm.


title: Pulling on multiple threads of evidence to reveal tree species coexistence

abstract: Understanding drivers of species coexistence is a central challenge in ecology (1, 2). Coexistence cannot be observed directly, and while species co-occurrence in time and space is necessary for coexistence, it is not sufficient to prove coexistence (3, 4). Species exclusion from a region is potentially observable, but can take decades to occur, and still might just occur stochastically (5, 6). Thus, ecologists generally use theory to identify indirect observations that are indicative of mechanisms driving coexistence or exclusion (7-9). Various methods have been developed to indirectly infer coexistence, each of which requires different data, none of which are usually conclusive on their own, and are not typically combined (e.g. (9-12)). Here, we combine three different approaches which are normally used in isolation to build a case about a lack of coexistence of multiple hardwood species. Importantly, though each method has its flaws, we find that all three point to the same conclusion which strengthens the case for a lack of coexistence. First, in an experimental planting of three mature tree species we found no relationship between productivity and species diversity, which has been hypothesized to occur due to a lack of niche differences among species. Second, we used modern coexistence theory to calculate niche and fitness differences for each pair of species, which confirmed the high niche overlap among species and competitive dominance of one species. Third, we used the United States Department of Agriculture Forest Inventory and Analysis data to examine co-occurrence patterns of our species across thousands of natural forest stands and found that, these three species were distributed randomly throughout the USA, and do not significantly co-occur as would be expected if coexistence occurred naturally. Given that these independent methods agree; we take this as strong evidence about a lack of coexistence.


title: Multiple lines of evidence combine to understand tree species coexistence


title: Statistics of the number of equilibria in random social dilemma evolutionary games with mutation

abstract: In this paper, we study analytically the statistics of the number of equilibria in pairwise social dilemma evolutionary games with mutation where a games payoff entries are random variables. Using the replicator-mutator equations, we provide explicit formulas for the probability distributions of the number of equilibria as well as other statistical quantities. This analysis is highly relevant assuming that one might know the nature of a social dilemma game at hand (e.g., cooperation vs coordination vs anti-coordination), but measuring the exact values of its payoff entries is difficult. Our delicate analysis shows clearly the influence of the mutation probability on these probability distributions, providing insights into how varying this important factor impacts the overall behavioural or biological diversity of the underlying evolutionary systems.


title: When Does Mutualism Offer a Competitive Advantage? A Game-Theoretic Analysis of Host-Host Competition in Mutualism

abstract: Plants due to their non-motile nature rely heavily on mutualistic interactions to obtain resources and carry out services. One key mutualism is the plant-microbial mutualism in which a plant trades away carbon to a microbial partner for nutrients like nitrogen and phosphorous. Plants show much variation in the use of this partnership from the individual level to entire lineages depending upon ecological, evolutionary, and environmental context. We sought to determine how this context dependency could result in the promotion, exclusion, or coexistence of the microbial mutualism by seeing if and when the partnership provided a competitive advantage to the plant. To that end, we created a simple 2 x 2 evolutionary game in which plants could either be a mutualist and pair with a microbe or a non-mutualist and forgo the partnership. This model included nutrients freely available to the plant, nutrients obtained only through mutualism with microbes, the cost of producing roots, the cost of trade with microbes, and the size of the local competitive neighborhood. Not surprisingly, we found that mutualism could offer a competitive advantage if its net benefit was positive. Coexistence between strategies is possible though due to competition between mutualists over the microbially obtained nutrient. In addition, the greater the size of the local competitive neighborhood, the greater the region of coexistence but only at the expense of mutualist fixation (non-mutualist fixation was unaffected). Our model, though simple, shows that plants can gain a competitive advantage from using a mutualism depending upon the context and points to basic experiments that can be done to verify the results.


title: Benchmarking ultra-high molecular weight DNA preservation methods for long-read and long-range sequencing

abstract: Studies in vertebrate genomics require sampling from a broad range of tissue types, taxa, and localities. Recent advancements in long-read and long-range genome sequencing have made it possible to produce high-quality chromosome-level genome assemblies for almost any organism. However, adequate tissue preservation for the requisite ultra-high molecular weight DNA (uHMW DNA) remains a major challenge. Here we present a comparative study of preservation methods for field and laboratory tissue sampling, across vertebrate classes and different tissue types. We find that no single method is best for all cases. Instead, the optimal storage and extraction methods vary by taxa, by tissue, and by down-stream application. Therefore, we provide sample preservation guidelines that ensure sufficient DNA integrity and amount required for use with long-read and long-range sequencing technologies across vertebrates. Our best practices generated the uHMW DNA needed for the high-quality reference genomes for Phase 1 of the Vertebrate Genomes Project (VGP), whose ultimate mission is to generate chromosome-level reference genome assemblies of all [~]70,000 extant vertebrate species.


title: PAC: Highly accurate quantification of allelic gene expression for population and disease genetics

abstract: Analysis of allele-specific gene expression (ASE) is a powerful approach for studying gene regulation. However, detection of ASE events relies on accurate alignment of RNA-sequencing reads, where challenges still remain. We have developed PAC, a method that combines multiple steps to improve the quantification of allelic reads, including personalised (i.e. diploid) read alignment with improved allocation of multi-mapping reads. We show that PAC outperforms standard alignment approaches for ASE detection in both accuracy and in the number of sites it can reliably quantify.


title: Robust Acquisition of High Resolution Spatial Transcriptomes from Preserved Tissues with Immunofluorescence Based Laser Capture Microdissection

abstract: The functioning of tissues is fundamentally dependent upon not only the phenotypes of the constituent cells but also their spatial organization in the tissue. However, acquisition of comprehensive transcriptomes of spatially- and phenotypically-defined cells in situ remains challenging. Here we present a general and robust method based on immunofluorescence-guided laser capture microdissection (immuno-LCM-RNAseq) to acquire finely resolved spatial transcriptomes including isoforms with as few as tens of cells from snap-frozen or RNAlater-treated clinical tissues, circumventing the problem of significant RNA degradation during this time-consuming process. The efficacy of this approach is exemplified by the characterization of differences at the transcript isoform level between the mouse small intestine lacteal cells at the tip versus the main capillary body. With the extensive repertoire of specific antibodies that are presently available, our method provides a powerful means by which spatially resolved cellular states can now be delineated in situ with preserved tissues. Moreover, such high quality spatial transcriptomes defined by immunomarkers can be used to compare with the phenotypes derived from single-cell RNAseq of dissociated cells as well as applied to bead-based spatial transcriptomic approaches that require such information a priori for cell identification.


title: Active enhancers strengthen TAD insulation by bRNA mediated CTCF enrichment at the TAD boundaries

abstract: Vertebrate genomes are partitioned into Topologically Associating Domains (TADs), which are typically bound by head-to-head pairs of CTCF binding sites. Transcription at TAD boundaries correlates with better insulation, however, it is not known whether the boundary transcripts themselves contribute to boundary function. Here we characterize boundary-associated RNAs genome-wide, focusing on the disease-relevant INK4a/ARF TAD. Using CTCF site deletions and boundary-associated RNA knockdowns, we observe that boundary-associated RNAs facilitate recruitment and clustering of CTCF at TAD borders. The resulting CTCF enrichment enhances TAD insulation, enhancer:promoter interactions and TAD gene expression. Importantly, knockdown of boundary-associated RNAs results in loss of boundary insulation function. Using enhancer deletions and CRISPRi of promoters we show that active TAD enhancers but not promoters induce boundary-associated RNA transcription, thus defining a novel class of regulatory enhancer RNAs.


title: Active enhancers strengthen insulation by RNA-mediated CTCF binding at TAD boundaries


title: Equatorial to Polar genomic variability of the microalgae Bathycoccus prasinos

abstract: Phytoplankton plays a fundamental role in the ecology of ocean systems and is the key player in the global carbon cycle. At a time of global warming, understanding the mechanisms of its adaptation to temperature is therefore of paramount importance. Cosmopolitan planktonic species abundant in different marine environments provide both a unique opportunity and an efficient methodological tool to study the genomic bases of their adaptation. This is the case for the eukaryotic picoalga Bathycoccus prasinos, whose genomic variability we chose to study in temperate and polar oceanic waters. Using multiple metagenomic datasets, we found that ~5% of B. prasinos genomic positions are variable, with an overwhelming majority of biallelic motifs. Cold and temperate waters are clearly associated with changes in variant frequencies, whereas in transitional waters we found more balanced polymorphism at most of these positions. Mesophilic and psychrophilic gene variants are distinguished by only a few amino acid changes located at positions critical for physical and functional protein properties. These results provide new information on the genomic diversity of a cosmopolitan eukaryotic planktonic specie and reveal "minimal mutational strategies" which finely tune the properties of specific proteins at different temperatures.


title: Prolactin and prolactin receptor expression in the HPG axis and crop sac during parental care in both sexes of a biparental bird (Columba livia)

abstract: During breeding, multiple circulating hormones, including prolactin, facilitate reproductive transitions in species that exhibit parental care. Prolactin underlies parental behaviors and related physiological changes across many vertebrates, including birds and mammals. While circulating prolactin levels often fluctuate across breeding, less is known about how relevant target tissues vary in their prolactin responsiveness via prolactin receptor (PRLR) expression. Recent studies have also investigated prolactin (PRL) gene expression outside of the pituitary (i.e., extra-pituitary PRL), but how PRL gene expression varies during parental care in non-pituitary tissue (e.g., hypothalamus, gonads) remains largely unknown. Further, it is unclear if and how tissue-specific PRL and PRLR vary between the sexes during biparental care. To address this, we measured PRL and PRLR gene expression in tissues relevant to parental care, the endocrine reproductive hypothalamic-pituitary-gonadal (HPG) axis and the crop (a tissue with a similar function as the mammalian mammary gland), across various reproductive stages in both sexes of a biparental bird, the rock dove (Columba livia). We also assessed how these genes responded to changes in offspring presence by adding chicks mid-incubation, simulating an early hatch when prolactin levels were still moderately low. We found that pituitary PRL expression showed similar increases as plasma prolactin levels, and detected extra-pituitary PRL in the hypothalamus, gonads and crop. Hypothalamic and gonadal PRLR expression also changed as birds began incubation. Crop PRLR expression correlated with plasma prolactin, peaking when chicks hatched. In response to replacing eggs with a novel chick mid-incubation, hypothalamic and gonadal PRL and PRLR gene expression differed significantly compared to mid-incubation controls, even when plasma prolactin levels did not differ. We also found sex differences in PRL and PRLR that suggest gene expression may allow males to compensate for lower levels in prolactin by upregulating PRLR in all tissues. Overall, this study advances our understanding of how tissue-specific changes in responsiveness to parental hormones may differ across key reproductive transitions, in response to offspring cues, and between the sexes.

category: animal behavior and cognition


title: Aspiration-assisted freeform bioprinting of mesenchymal stem cell spheroids within alginate microgels

abstract: Aspiration-assisted freeform bioprinting (AAfB) has emerged as a promising technique for precise placement of tissue spheroids in three-dimensional (3D) space for fabrication of tissues. For successful embedded bioprinting using AAfB, an ideal support bath should possess shear-thinning behavior and yield-stress to obtain tightly fused assembly of bioprinted spheroids. Several studies have demonstrated support baths for embedded bioprinting, but these materials pose major challenges due to their low biocompatibility, opaqueness, complex and prolonged preparation procedures, and limited spheroid fusion efficacy. In this study, to circumvent the aforementioned limitations, we present the feasibility of AAfB of human mesenchymal stem cell (hMSC) spheroids in alginate microgels as a support bath. First, alginate microgels were prepared with different particle sizes modulated by blending time and concentration, followed by determination of the optimal bioprinting conditions by the assessment of rheological properties, bioprintability, and spheroid fusion efficiency. The bioprinted and consequently self-assembled tissue structures made of hMSC spheroids were osteogenically induced for bone tissue formation. Alongside, we investigated the effects of peripheral blood monocyte-derived osteoclast incorporation into the hMSC spheroids in heterotypic bone tissue formation. We demonstrated that alginate microgels enabled unprecedented positional accuracy (~5%), transparency for visualization, and improved fusion efficiency (~97%) of bioprinted hMSC spheroids for bone fabrication. This study demonstrates the feasibility of using alginate microgels as a support bath for many different applications including but not limited to freeform bioprinting of spheroids, cell-laden hydrogels, and fugitive inks to form viable tissue constructs.

category: bioengineering


title: A Structure-based B-cell Epitope Prediction Model Through Combing Local and Global Features

abstract: B-cell epitopes (BCEs) are a set of specific sites on the surface of an antigen that binds to an antibody produced by B-cell. The recognition of BCEs is a major challenge for drug design and vaccines development. Compared with experimental methods, computational approaches have strong potential for BCEs prediction at much lower cost. Moreover, most of the currently methods focus on using local information around target residue without taking the global information of the whole antigen sequence into consideration. We propose a novel deep leaning method through combing local features and global features for BCEs prediction. In our model, two parallel modules are built to extract local and global features from the antigen separately. For local features, we use Graph Convolutional Networks(GCNs) to capture information of spatial neighbors of a target residue. For global features, Attention-Based Bidirectional Long Short-Term Memory(Att-BLSTM) networks are applied to extract information from the whole antigen sequence. Then the local and global features are combined to predict BCEs. The experiments show that the proposed method achieves superior performance over the state-of-the-art BCEs prediction methods on benchmark datasets. Also, we compare the performance differences between data with or without global features. The experimental results show that global features play an important role in BCEs prediction. Our detailed case study on the BCEs prediction for SARS-Cov-2 receptor binding domain confirms that our method is effective for predicting and clustering true BCEs.


title: MetaFun: Unveiling sex differences in multiple omics studies through comprehensive functional meta-analysis

abstract:
SummarySex and gender differences in different health scenarios has been thoroughly acknowledged in the literature, and yet, very scarcely analyzed. To fill the gap, here we present MetaFun, which allows to meta-analyze multiple omics datasets with a sex-based perspective, and to combine different datasets to gain major statistical power and to assist the researcher in understanding these sex differences in the diseases under study. Metafun is freely available at bioinfo.cipf.es/metafun Availability and implementation: MetaFun is available under http://bioinfo.cipf.es/metafun. The backend has been implemented in R and Java and the frontend has been developed using Angular.

Supplementary informationR code available at https://gitlab.com/ubb-cipf/metafunr


title: Binding Free energy Decomposition and Multiple Unbinding Paths of Buried Ligands in a PreQ1 Riboswitch

abstract: Riboswitches are naturally occurring RNA elements that control bacterial gene expression by binding to specific small molecules. They serve as important models for RNA-small molecule recognition and have also become a novel class of targets for developing antibiotics. Here, we carried out conventional and enhanced-sampling molecular dynamics (MD) simulations, totaling 141.5 s, to characterize the determinants of binding free energies and unbinding paths for the cognate and synthetic ligands of a PreQ1 riboswitch. Binding free energy analysis showed that two triplets of nucleotides U6-C15-A29 and G5-G11-C16, contribute the most to the binding of the cognate ligands, by hydrogen bonding and by base stacking, respectively. Mg2+ ions are essential in stabilizing the binding pocket. For the synthetic ligands, the hydrogen-bonding contributions of the U6-C15-A29 triplet are significantly compromised, and the bound state resembles the apo state in several respects, including the disengagement of the C15-A14-A13 and A32-G33 base stacks. The bulkier synthetic ligands lead to significantly loosening of the binding pocket, including extrusion of the C15 nucleobase and a widening of the C15-C30 groove. Enhanced-sampling simulations further revealed that the cognate and synthetic ligands unbind in almost opposite directions. Our work offers new insight for designing riboswitch ligands.


title: Infinite re-reading of single proteins at single-amino-acid resolution using nanopore sequencing

abstract:
As identifying proteins is of paramount importance for cell biology and applications, it is of interest to develop a protein sequencer with the ultimate sensitivity of decoding individual proteins. Here, we demonstrate a nanopore-based single-molecule sequencing approach capable of reliably detecting single amino-acid substitutions within individual peptides. A peptide is linked to a DNA molecule that is pulled through the biological nanopore MspA by a DNA helicase in single amino-acid steps. The peptide sequence yields clear stepping ion current signals which allows to discriminate single-amino-acid substitutions in single reads. Molecular dynamics simulations show these signals to result from size exclusion and pore binding. Notably, we demonstrate the capability to  rewind peptide reads, obtaining indefinitely many independent reads of the same individual molecule, yielding virtually 100% read accuracy in variant identification, with an error rate less than 10-6. These proof-of-concept experiments constitute a promising basis for developing a single-molecule protein sequencer.

One-sentence summaryThis paper presents proof-of-concept experiments and simulations of a nanopore-based approach to sequencing individual proteins.


title: Capsaicin alters human NaV1.5 mechanosensitivity

abstract: SCN5A-encoded NaV1.5 is a voltage-gated Na+ channel expressed in cardiac myocytes and human gastrointestinal (GI) smooth muscle cells (SMCs). NaV1.5 contributes to electrical excitability in the heart and slow waves in the gut. NaV1.5 is also mechanosensitive, and mechanical force modulates several modes of NaV1.5s voltage-dependent function. NaV1.5 mutations in patients with cardiac arrhythmias and gastrointestinal diseases lead to abnormal mechano- and voltage-sensitivity. Membrane permeable amphipathic drugs that target NaV1.5 in the heart and GI tract alter NaV1.5 mechanosensitivity (MS), suggesting that amphipaths may be a viable therapeutic option for modulating NaV1.5 function. We therefore searched for membrane-permeable amphipathic agents that would modulate NaV1.5 MS with minimal effect on NaV1.5 voltage-gating intact to more selectively target mechanosensitivity. We used two methods to assess NaV1.5 MS: (1) membrane suction in cell-attached macroscopic patches and (2) fluid shear stress on whole cells. We tested the effect of capsaicin on NaV1.5 MS by examining macropatch and whole-cell Na+ current parameters with and without force. The pressure- and shear-mediated peak current increase and acceleration were effectively abolished by capsaicin. Capsaicin abolished the mechanosensitive shifts in the voltage-dependence of activation (shear) and inactivation (pressure and shear). Exploring the recovery from inactivation and use-dependent entry into inactivation, we found divergent stimulus-dependent effects that could potentiate or mitigate the effect of capsaicin, suggesting that mechanical stimuli may differentially modulate NaV1.5 MS. We conclude that selective modulation of MS makes capsaicin is a novel modulator of NaV1.5 MS and a promising therapeutic candidate.


title: Stochastic model of T Cell repolarization during target elimination (II)

abstract: Cytotoxic T lymphocytes (T cells) and natural killer cells form a tight contact, the immunological synapse (IS), with target cells, where they release their lytic granules containing perforin/granzyme and cytokine containing vesicles. During this process the cell repolarizes and moves the microtubule organizing center (MTOC) towards the IS. In the first part of our work we developed a computational model for the molecular-motor-driven motion of the MT cytoskeleton during T cell polarization and analyzed effects of cortical sliding and capture-shrinkage mechanisms. Here we use this model to analyze the dynamics of the MTOC repositioning in situations in which a) the IS is in an arbitrary position with respect to the initial position of the MTOC and b) the T cell has two IS at two arbitrary positions. In the case of one IS, we found that the initial position determines which mechanism is dominant and that the time of repositioning does not rise monotonously with the MTOC-IS distance. In the case of two IS, we observe several scenarios that have also been reported experimentally: the MTOC alternates stochastically (but with a well defined average transition time) between the two IS; it wiggles in between the two IS without transiting to one of the two; or it is at some point pulled to one of the two IS and stays there. Our model allows to predict which scenario emerges in dependency of the mechanisms in action and the number of dyneins present. We report that the presence of capture-shrinkage mechanism at least in one IS is necessary to assure the transitions in every cell configuration. Moreover, the frequency of transitions does not decrease with the distance between the two IS and is the highest when both mechanisms are present in both IS.


title: Dual effects of presynaptic membrane mimetics on α-synuclein amyloid aggregation

abstract: Aggregation of intrinsically disordered -synuclein (SN) under various conditions is closely related to synucleinopathies. Although various biological membranes have shown to alter the structure and aggregation propensity of SN, a thorough understanding of the molecular and mechanical mechanism of amyloidogenesis in membranes remains unanswered. Herein, we examined the structural changes, binding properties, and amyloidogenicity of three variations of SN mutants under two types of liposomes, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and presynaptic vesicle mimetic (Mimic) membranes. While neutrally charged DOPC membranes elicited marginal changes in the structure and amyloid fibrillation of SNs, negatively charged Mimic membranes induced dramatic helical folding and biphasic amyloid generation. At low concentration of Mimic membranes, the amyloid fibrillation of SNs was promoted in a dose-dependent manner. However, further increases in the concentration constrained the fibrillation process. These results suggest the dual effect of Mimic membranes on regulating the amyloidogenesis of SN, which is rationalized by the amyloidogenic structure of SN and condensation-dilution of local SN concentration. Finally, we propose physicochemical properties of SN and membrane surfaces, and their propensity to drive electrostatic interactions as decisive factors of amyloidogenesis.


title: Aberrant enteric neuromuscular system and dysbiosis in amyotrophic lateral sclerosis

abstract:
BackgroundEmerging evidence has demonstrated that microbiota directly affects the enteric neuron system (ENS) and smooth muscle cell functions via metabolic products or endogenous bacterial components. Amyotrophic Lateral Sclerosis is a neuromuscular disease characterized by the progressive death of motor neurons and muscle atrophy. The GI symptoms in patients were largely ignored or underestimated, especially before the diagnosis of ALS. The relationship between enteric neuromuscular system and microbiome in ALS progression is unknown.

MethodsWe performed longitudinal studies on the ENS and microbiome in the ALS human-SOD1G93A transgenic G93A mice. We treated age-matched wild-type and ALS mice with bacterial product butyrate or antibiotics to investigate microbiome and neuromuscular functions. Intestinal motility, microbiome, an ENS marker GFAP, a smooth muscle marker (SMMHC), and human colonoids have been examined. The distribution of human-G93A-SOD1 (Superoxide Dismutase 1) protein was tested as an indicator of ALS progression.

ResultsAt 2-month-old before ALS onset, G93A mice had significant lower intestinal motility, decreased grip strength, and reduced time in the rotarod. We observed increased GFAP and decreased SMMHC expression. These changes correlated with consistent increased aggregation of mutated SOD1G93A in the colon, small intestine, and spinal cord. Butyrate and antibiotic treatment showed a significantly longer latency to fall in the rotarod test, reduced SOD1G93A aggregation, and enhanced ENS and muscle function. Feces from 2-month-old SOD1G93A mice significantly enhanced SOD1G93A aggregation in human colonoids transfected with a SOD1G93A-GFP plasmid. Longitudinal studies of microbiome data further showed the altered bacterial community related with autoimmunity (e.g., Clostridium sp. ASF502, Lachnospiraceae bacterium A4), inflammation (e.g., Enterohabdus Muris,), and metabolism (e.g., Desulfovibrio fairfieldensis) at 1- and 2-month-old SOD1G93A mice, suggesting the early microbial contribution to the pathological changes.

ConclusionsWe have demonstrated a novel link between microbiome, hSOD1G93A aggregation, and intestinal mobility. Dysbiosis occurred at the early stage of the ALS mice before observed mutated-SOD1 aggregation, slow intestinal motility, and dysfunction of ENS. Manipulating the microbiome improves the muscle performance of SOD1G93A mice. Our study provides insights into fundamentals of intestinal neuromuscular structure/function and microbiome in ALS.


title: Engineering sorghum for higher 4-hydroxybenzoic acid content

abstract: Lignocellulosic biomass represents a renewable source of sugars for the manufacturing of bioproducts such as biofuels. The high cost associated with deconstruction of plant biomass to simple sugars remains one of the challenges preventing the deployment of economically sustainable advanced bioproducts. The accumulation in-planta of value-added coproducts such as platform chemicals can improve the economics of biofuels. Among other crops, sorghum is an ideal bioenergy feedstock due to its low input requirements, efficient nitrogen recycling, and high water use efficiency and biomass yields. In this work, we engineered sorghum to overproduce the valuable chemical 4-hydroxybenzoic acid (4-HBA) by co-expressing plastid-targeted versions of Escherichia coli chorismate pyruvate-lyase (UbiC) and feedback-resistant 3-deoxy-D-arabino-heptulonate-7-phosphate synthase (AroG*). Two independent lines containing the dual aroG*-ubiC construct were selected for characterization in the T2 generation. Using liquid chromatography-mass spectrometry, analysis of methanolic extracts obtained from biomass samples revealed an accumulation of 4-HBA in the two transgenic lines (corresponding to 1.56% and 1.72% dry weight, respectively), with 4-HBA glucose conjugates representing major forms. Measurements of biomass composition and several agronomic traits showed no significant difference between transgenic and wild-type control plants grown under controlled environment. This work demonstrates the transferability of the ubiC engineering approach to sorghum; generated lines will be useful to assess the agronomic performance of modified 4-HBA-rich sorghum under natural conditions.


title: Image-based assessment of plant disease progression identifies new genetic loci for resistance

abstract: A major challenge in global crop production is mitigating yield loss due to plant diseases. One of the best means of disease control is plant resistance, but the identification of genes that promote resistance has been limited by the subjective quantification of disease, which is typically scored by the human eye. We hypothesized that image-based, non-destructive quantification of disease phenotypes would enable the rapid identification of new disease resistance loci. We tested this using the interaction between tomato and Ralstonia solanacearum, a soilborne pathogen that causes bacterial wilt disease. We acquired over 40,000 time-series images of disease progression in a tomato recombinant inbred line population, and developed an image analysis pipeline providing a suite of ten traits to quantify wilt disease based on plant shape and size. Quantitative trait loci (QTL) analyses using image-based phenotyping identified QTL that were both unique and shared compared with those identified by human assessment of wilting. When shared loci were identified, image-based phenotyping could detect some QTL several days earlier than human assessment. Thus, expanding the phenotypic space of disease with image-based, non- destructive phenotyping allowed both earlier detection and identified new genetic components of resistance.


title: Inferring viral occurrence patterns through a synthetic data simulation

abstract: Viruses infecting humans are manifold and several of them provoke significant morbidity and mortality. Simulations creating large synthetic datasets from observed multiple viral strain infections in a limited population sample can be a powerful tool to infer significant pathogen occurrence and interaction patterns, particularly if limited number of observed data units is available. Here, to demonstrate diverse human papillomavirus (HPV) strain occurrence patterns, we used log-linear models combined with Bayesian framework for graphical independence network (GIN) analysis. That is, to simulate datasets based on modeling the probabilistic associations between observed viral data points, i.e different viral strain infections in a set of population samples. Our GIN analysis outperformed in precision all oversampling methods tested for simulating large synthetic viral strain-level prevalence dataset from observed set of HPVs data. Altogether, we demonstrate that network modeling is a potent tool for creating synthetic viral datasets for comprehensive pathogen occurrence and interaction pattern estimations.


title: gpuZoo: Cost-effective estimation of gene regulatory networks using the Graphics Processing Unit

abstract:
Gene regulatory network inference allows for the study of transcriptional control to identify the alteration of cellular processes in human diseases. Our group has developed several tools to model a variety of regulatory processes, including transcriptional (PANDA, SPIDER) and post-transcriptional (PUMA) gene regulation, and gene regulation in individual samples (LIONESS). These methods work by performing repeated operations on data matrices in order to integrate information across multiple lines of biological evidence. This limits their use for large-scale genomic studies due to the associated high computational burden. To address this limitation, we developed gpuZoo, which includes GPU-accelerated implementations of these algorithms. The runtime of the gpuZoo implementation in MATLAB and Python is up to 61 times faster and 28 times less expensive than the multi-core CPU implementation of the same methods. gpuZoo takes advantage of the modern multi-GPU device architecture to build a population of sample-specific gene regulatory networks with similar runtime and cost improvements by combining GPU acceleration with an efficient on-line derivation. Taken together, gpuZoo allows parallel and on-line gene regulatory network inference in large-scale genomic studies with cost-effective performance.

gpuZoo is available in MATLAB through the netZooM package https://github.com/netZoo/netZooM and in Python through the netZooPy package https://github.com/netZoo/netZooPy.

category: systems biology