title: Conserved patterns across ion channels correlate with variant pathogenicity and clinical phenotypes

abstract:
Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. Consequently, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated.

We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion channel families. We collected and curated 3,049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12,546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures.

We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our novel 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5[A] distance from the pore axis center and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and fucntional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate (NMDA) receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1,422 neurodevelopmental disorder pathogenic patient variants, and 679 electrophysiological experiments that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain vs. loss-of-channel function.

In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional read-outs and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future.


title: An integrative analysis of the biological clock hypothesis in human gut microbiome

abstract:
BackgroundWhile previous studies have explored the relationship between aging and the gut microbiome, it remains unclear how consistent and reproducible this association is across different cultures and groups. We performed an integrative analysis with 11 independent datasets from nine different countries to test the idea that the aging gut microbiome can be viewed as a biological clock, in which microbial changes associated with age are consistent and measurable across distinct datasets.

ResultsAs expected, our Principal Coordinate Analysis found strong batch effects with study ID by far the strongest signal across datasets. Despite this large batch effect, we found a consistent signal across studies that was largely driven by sample size with only our larger cohorts showing taxa in common associated with age. Likewise, Shannon diversity and richness were not consistently associated with age across the 11 datasets, but some positive correlations with richness and host age were observed among the four largest cohorts. The taxon with the most potential as a biomarker for the aging human gut microbiome is genus Bifidobacterium, with significantly negative associations with host age in three out of the four datasets that had more than 200 samples.

ConclusionThe driving force behind low reproducibility of association of age with the microbiome in previous studies appears to be inadequate sample size rather than structural differences in the microbial community based on cohort characteristics. Results from a power analysis suggest that future studies on the aging human gut microbiome consider on the order of 100-300 samples to consistently observe an age signal. With these larger sample sizes, parametric and non-parametric model yield broadly similar power.

category: bioinformatics


title: Seqrutinator: Non-Functional Homologue Sequence Scrutiny for the Generation of large Datatsets for Protein Superfamily Analysis

abstract:
BackgroundIn recent years protein protein bioinformatics has resulted in much improved algorithms for multiple sequence alignment (MSA) and phylogeny. Few attention has been paid to sequence selection whereas particularly recently published complete proteomes often have many sequences that are partial or derive from pseudogenes. Not only do these sequences add noise to the MSA, phylogeny and other downstream computational analyses, they actually also instigate many errors in the processing of the MSAs and therewith all other downstream analyses, including the phylogeny.

ObjectiveHere we aim to provide and test an objective, automated and flexible pipeline for the scrutiny of sequence sets from large, complex protein superfamilies, resulting in the classification of sequences into either functional or non-functional. The pipeline should identify few Swissprot sequences and result in an improved MSA.

ResultsSeqrutinator is a pipeline that consists of five modules written in Python that identify and remove sequences that are likely Non-Functional Homologues. Here we tested the pipeline in its default order using default settings using three complex, related plant superfamilies (BAHD, CYP and UGT) that act in specialized metabolism, using the complete proteome of 16 plant species as input and using Swissprot as a control. Only 3.4% of Swissprot sequences were identified as NFH. Most NFH sequences are partial but, interestingly, their removal results in much shorter MSAs with many more reliable columns. Few but significant sequences that instigate gaps are found. The three superfamilies furthermore show similar scrutiny patterns.

Conclusion and perspectiveSeqrutinator forms a consistent pipeline for sequence scrutiny that does result in sequence sets that generate high fidelity MSAs. Additional tests will be performed in order to further demonstrate Seqrutinators performance.


title: The traits of skull anatomy in green frogs of the genus Pelophylax (Anura, Ranidae)

abstract: During dissecting skulls of green frogs of the genus Pelophylax it has found that the description of a number bones and cartilaginous structures requires addition and clarification. In addition, previously undescribed details of anatomical structure of the skull of green frogs has revealed. The traits of anatomical structure of skull of green frogs, characteristic for ancestral forms of Anura found too.

category: zoology


title: The Rise and Fall of Slow Wave Tides: Vacillations of Slow Wave/Spindle Coupling Shift the Composition of Slow Wave Activity Through Sleep Cycles in Accordance with Depth of Sleep

abstract: Slow wave activity (SWA) during sleep is associated with synaptic regulation and memory processing functions. Each cycle of non-rapid-eye-movement (NREM) sleep demonstrates a waxing and waning amount of SWA during the transitions between stages N2 and N3 sleep, and the deeper N3 sleep is associated with an increased density of SWA. Further, SWA is an amalgam of different types of slow waves, each identifiable by their temporal coupling to spindle subtypes with distinct physiological features. The objectives of this study were to better understand the neurobiological properties that distinguish different slow wave and spindle subtypes, and to examine the composition of SWA across cycles of NREM sleep. We further sought to explore changes in the composition of NREM cycles that occur among aging adults. To address these goals, we analyzed subsets of data from two well-characterized cohorts of healthy adults: 1) The DREAMS Subjects Database (n=20), and 2) The Cleveland Family Study (n=60). Our analyses indicate that slow wave/spindle coupled events can be characterized as frontal versus central in their relative distribution between electroencephalography (EEG) channels. The frontal predominant slow waves are identifiable by their coupling to late-fast spindles and occur more frequently during stage N3 sleep. Conversely, the central-associated slow waves are identified by coupling to early-fast spindles and favor occurrence during stage N2 sleep. Together, both types of slow wave/spindle coupled events form the composite of SWA, and their relative contribution to the SWA rises and falls across cycles of NREM sleep in accordance with depth of sleep. Exploratory analyses indicated that older adults produce a different composition of SWA, with a shift toward the N3, frontal subtype, which becomes increasingly predominant during cycles of NREM sleep. Overall, these data demonstrate that subtypes of slow wave/spindle events have distinct cortical propagation patterns and differ in their distribution across lighter versus deeper NREM sleep. Future efforts to understand how slow wave sleep and slow wave/spindle coupling impact memory performance and neurological disease may benefit from examining the composition of SWA to avoid potential confounds that may occur when comparing dissimilar neurophysiological events.

category: neuroscience


title: Microglial Piezo1 senses Aβ fibrils stiffness to restrict Alzheimer's disease

abstract: The pathology of Alzheimers disease (AD) is associated with the extracellular amyloid-{beta} (A{beta}) plaques that perturb the mechanical properties of brain tissue. Microglia sense and integrate biochemical cues in their local microenvironment, intimate linking with AD progress. However, neither the microglial mechanosensing pathway nor its impact on AD pathogenesis is well studied. Here, we showed that the mechanosensitive ion channel Piezo1 is increased in microglia upon stiffness stimuli of A{beta} fibrils. The upregulation of Piezo1 in A{beta} plaque-associated microglia was observed in AD mouse models and human patients. Microglia lacking Piezo1 disturbed microglial clustering, phagocytosis, and compaction of A{beta} plaques, resulting in the exacerbation of A{beta} and neurodegenerative pathologies in AD. Conversely, pharmacological activation of Piezo1 ameliorated brain A{beta} burden and cognitive impairment in the 5xFAD mouse. Together, our results reveal Piezo1, as a mechanosensor of A{beta} fibrils stiffness in microglia, could represent a promising therapeutic target for AD.


title: Connexin36 expression in the cochlear nucleus complex of the echolocating bat, Eptesicus fuscus

abstract: Gap junctions and electrical synapses in the central nervous system are associated with rapid temporal processing and coincidence detection. Using histology, immunohistochemistry, and in situ hybridization, we investigated the distribution of Connexin36 (Cx36), a protein that comprises neuronal gap junctions, throughout the cochlear nucleus complex of the echolocating big brown bat, Eptesicus fuscus, a species exhibiting extreme behavioral sensitivity to minute temporal changes in ultrasonic echoes. For comparison, we visualized Cx36 expression in the cochlear nucleus of transgenic Cx36 reporter mice, species that hear ultrasound but do not echolocate. We observed Cx36 expression in the anteroventral and dorsal cochlear nucleus, with more limited expression in the posteroventral cochlear nucleus, of both species. Several different morphological cell types were labeled, including globular and spherical bushy, octopus, stellate, and fusiform cells. Labeled Cx36 puncta were also observed. Cx36 expression in the bat was spread throughout a relatively smaller area of the cochlear nucleus than in the mouse, even though the bat cochlear nucleus is hypertrophied. In the bat, the anteroventral cochlear nucleus showed higher percent area label than the dorsal cochlear nucleus, with a trend towards the opposite result in the mouse. The presence of gap junctions appears to be a conserved feature of the mammalian cochlear nucleus and thus not uniquely tied to the temporal hyperacuity of echolocation.


title: Marked regional glial heterogeneity in the human white matter of the central nervous system

abstract: The myelinated white matter tracts of the central nervous system (CNS) are essential for fast transmission of electrical impulses and are commonly affected in neurodegenerative diseases. However, these often uniquely human diseases differentially affect white matter regions, at various ages and between males and females, and we hypothesised that this is secondary to physiological variation in white matter glia with region, age and sex. Using single nucleus RNA sequencing of healthy human post-mortem samples, we find marked glial heterogeneity with tissue region (primary motor cortex, cerebellum, cervical spinal cord), with tissue-specific cell populations of oligodendrocyte precursor cells and astrocytes, and a spinal cord-enriched oligodendrocyte type that appears human-specific. Spinal cord microglia but not astrocytes show a more activated phenotype compared to brain. These regional effects, with additional differentially expressed genes with age and sex in all glial lineages, help explain pathological patterns of disease - essential knowledge for therapeutic strategies.


title: CD4 binding-site antibodies induced by a subtype B HIV-1 envelope trimer

abstract: In the last decade considerable advances have been made towards the design of HIV-1 vaccines that induce neutralizing antibodies (NAbs). Despite the progress, no vaccine is able to consistently elicited broadly neutralizing antibodies (bNAbs). Here we present a case study of a rabbit that was immunized with a subtype B native like envelope glycoprotein (Env) trimer, AMC016 SOSIP.v4.2, with a dense and intact glycan shield, followed by a trivalent combination of subtype B trimers. After the priming phase serum from this animal neutralized several heterologous subtype B neutralization resistant (tier 2) viruses. Subsequent immunization with the trivalent combination of subtype B trimers further increased the breadth and potency of the NAb response. EM based polyclonal epitope mapping revealed that a cross reactive CD4 binding-site (CD4bs) antibody response, that was present after priming with the monovalent trimer and boosting with the trivalent combination, was most likely responsible for the broad neutralization. While anecdotal, this study provides proof-of-concept that native-like Env trimers are capable of inducing CD4bs-directed bNAb responses and should guide efforts to improve the consistency with which such responses are generated.

category: immunology


title: Effects of mucin and DNA concentrations in airway mucus on Pseudomonas aeruginosa biofilm recalcitrance

abstract: The pathological properties of airway mucus in cystic fibrosis (CF) are dictated by mucus concentration and composition, with mucins and DNA being responsible for mucus viscoelastic properties. As CF pulmonary disease progresses, the concentrations of mucins and DNA increase and are associated with increased mucus viscoelasticity and decreased transport. Similarly, the biophysical properties of bacterial biofilms are heavily influenced by the composition of their extracellular polymeric substances (EPS). While the roles of polymer concentration and composition in mucus and biofilm mechanical properties have been evaluated independently, the relationship between mucus concentration and composition and the biophysical properties of biofilms grown therein remains unknown. Pseudomonas aeruginosa biofilms were grown in airway mucus as a function of overall concentration and DNA concentration to mimic healthy and CF physiology and biophysical properties were evaluated with macro- and microrheology. Biofilms were also characterized after exposure to DNase or DTT to examine the effects of DNA and mucin degradation, respectively. Identifying critical targets in biofilms for disrupting mechanical stability in highly concentrated mucus may lead to the development of efficacious biofilm therapies and ultimately improve CF patient outcomes. Overall mucus concentration was the predominant contributor to biofilm viscoelasticity and both DNA degradation and mucin reduction resulted in compromised biofilm mechanical strength.

category: microbiology


title: Efficient rescue of a newly classified Ebinur lake orthobunyavirus with GFP reporter and its application in rapid antiviral screening

abstract: Orthobunyaviruses have been reported to cause severe diseases in humans or animals, posing a threat to human health and social economy. Ebinur lake virus (EBIV) is a newly classified orthobunyavirus, which needs further intensive study and therapies to cope with its potential infection risk to human and animals. Here, through the reverse genetics system, the recombinant EBIV of wild type (rEBIV/WT) and NP-conjugated-eGFP (rEBIV/eGFP/S) were rescued for the application of the rapid antiviral drug screening. The eGFP fluorescence signal of the rEBIV/eGFP/S was stable in the process of successive passage in BHK-21 cells (over 10 passages) and this recombinant virus could replicate in various cell lines. Compared to the wild type EBIV, the rEBIV/eGFP/S caused the smaller plaques and its peak titers were lower, suggesting attenuation due to the eGFP insertion. Through the high-content screening (HCS) system, ribavirin showed an inhibitory effect on the rEBIV/eGFP/S with an EC50 of 21.91 M, while favipiravir did not inhibit, even at high concentrations. In addition, five of ninety-six natural compounds had antiviral against EBIV. The robust reverse genetics system for EBIV will facilitate investigation into replication and assembly mechanisms and assist drug and vaccine development.

category: molecular biology


title: Large-scale analysis and computer modeling reveal hidden regularities behind variability of cell division patterns in Arabidopsis thaliana embryogenesis

abstract: Noise plays a major role in cellular processes and in the development of tissues and organs. Several studies have examined the origin, the integration or the accommodation of noise in gene expression, cell growth and elaboration of organ shape. By contrast, much less is known about variability in cell division plane positioning, its origin and links with cell geometry, and its impact on tissue organization. Taking advantage of the first-stereotyped-then-variable division patterns in the embryo of the model plant Arabidopsis thaliana, we combined 3D imaging and quantitative cell shape and cell lineage analysis together with mathematical and computer modeling to perform a large scale, systematic analysis of variability in division plane orientation. Our results reveal that, paradoxically, variability in cell division patterns of Arabidopsis embryos is accompanied by a progressive reduction of cell shape heterogeneity. The paradox is solved by showing that variability operates within a reduced repertoire of possible division plane orientations that is related to cell geometry. We show that in several domains of the embryo, a recently proposed geometrical division rule recapitulates observed variable patterns, thus suggesting that variable patterns emerge from deterministic principles operating in a variable geometrical context. Our work highlights the importance of emerging patterns in the plant embryo under iterated division principles, but also reveal domains where deviations between rule predictions and experimental observations point to additional regulatory mechanisms.

category: developmental biology


title: Dynamics in vertical transmission of viruses in naturally selected and traditionally managed honey bee colonies across Europe

abstract: The  suppressed in-ovo virus infection trait (SOV) was the first trait applied in honey bee breeding programs aimed to increase resilience to virus infections, a major threat for colony survival. By screening drone eggs for viruses, the SOV trait scores the antiviral resistance of queens and its implications for vertical transmission. In this study, queens from both naturally surviving and traditionally managed colonies from across Europe were screened using a two-fold improved SOV phenotyping protocol. First, a gel-based RT-PCR was replaced by a RT-qPCR. This not only allowed quantification of the infection load but also increased the test sensitivity. Second, a genotype specific primer set was replaced by a primer set that covered all known deformed wing virus (DWV) genotypes, which resulted in higher virus loads and fewer false negative results. It was demonstrated that incidences of vertical transmission of DWV were more frequent in naturally surviving populations than in traditionally managed colonies, although the virus load in the eggs remained the same. Dynamics in vertical transmission were further emphasized when comparing virus infections with queen age. Interestingly, older queens showed significantly lower infection loads of DWV in both traditionally managed and naturally surviving colonies, as well as reduced DWV infection frequencies in traditionally managed colonies when compared with younger queens. Seasonal variation in vertical transmission was found with lower infection frequencies in spring compared to summer for DWV and black queen cell virus. Together, these patterns in vertical transmission suggest an adaptive antiviral response of queens aimed at reducing vertical transmission over time.

category: ecology


title: Population dynamics of free-roaming dogs in Europe and implications for population control

abstract: Changes in free-roaming dog population size are important indicators of the effectiveness of dog population management. Assessing the effectiveness of different management methods also requires estimating the processes that change population size, such as the rates of recruitment into and removal from a population. This is one of the first studies to quantify the size, rates of recruitment and removal, and health and welfare status of free-roaming dog populations in Europe. We determined the size, dynamics, and health status of free-roaming dog populations in Pescara, Italy, and Lviv, Ukraine, over a 15-month study period. Both study populations had ongoing dog population management through catch-neuter-release and sheltering programmes. Average monthly apparent survival probability was 0.93 (95% CI 0.81-1.00) in Pescara and 0.93 (95% CI 0.84-0.99) in Lviv. An average of 7 dogs km-2 were observed in Pescara and 40 dogs km-2 in Lviv. Per capita entry probabilities varied between 0.09 and 0.20 in Pescara, and 0.12 and 0.42 in Lviv. In Lviv, detection probability was lower on weekdays (odds ratio: 0.74, 95% CI 0.53-0.96) and higher on market days (odds ratio: 2.58, 95% CI 1.28-4.14), and apparent survival probability was lower in males (odds ratio: 0.25, 95% CI 0.03-0.59). Few juveniles were observed in the study populations, indicating that recruitment may be occurring by movement between dog subpopulations (e.g. from local owned or neighbouring free-roaming dog populations), with important consequences for population control. This study provides important data for planning effective dog population management and for informing population and infectious disease modelling. The results emphasise the importance of dog population management efforts that encourage responsible ownership through reduced abandonment and cover the entire city to reduce the risk of dog movement between sections of the city hindering management efforts.


title: Prevalence Estimates of Putatively Pathogenic Leptin Variants in the gnomAD Database

abstract: Homozygosity for pathogenic variants in the leptin gene leads to congenital leptin deficiency causing early-onset extreme obesity. This monogenic form of obesity has mainly been detected in patients from consanguineous families. Prevalence estimates for the general population using the Exome Aggregation Consortium (ExAC) database reported a low frequency of leptin mutations. One in approximately 15 million individuals will be homozygous for a deleterious leptin variant. With the present study, we aimed to extend these findings utilizing the augmented Genome Aggregation Database (gnomAD) v2.1.1 including more than 140,000 samples. In total, 68 non-synonymous and 7 loss-of-function (LoF) leptin variants were deposited in gnomAD. By predicting functional implications with the help of in silico tools, like SIFT, PolyPhen2 and MutationTaster2021, the prevalence of hetero- and homozygosity for putatively pathological variants (n = 32; pathogenic prediction by at least two tools) in the leptin gene were calculated. Across all populations, the estimated prevalence for heterozygosity for functionally relevant variants was approximately 1:2,100 and 1:17,860,000 for homozygosity. This prevalence deviated between the individual populations. Accordingly, people from South Asia were at greater risk to carry a possibly damaging leptin variant than individuals of other ancestries. Generally, this study emphasises the scarcity of deleterious leptin variants in the general population with varying prevalence for distinct study groups.

category: genetics


title: Genetics of Fertility Restoration in A1 Cytoplasmic Genetic Male Sterility (CGMS) Systems in Sorghum (Sorghum bicolor L. Moench)

abstract: The cause of Cytoplasmic Genetic Male Sterility (CGMS) is specific nuclear and mitochondrial interactions. Almost all commercial sorghum hybrids were developed using the A1 cytoplasmic genetic male sterility system. Understanding the inheritance of fertility restoration in sorghum for A1 cytoplasm, for example, can improve the selection efficiency of restorer lines for increased seed production. In a cross of male sterile line 296A with A1 cytoplasm and restorer lines comprised of a set of Recombinant Inbred Lines (RILs), the inheritance pattern of fertility restoration of sorghum was studied. The F1 hybrid was completely fertile, revealing the dominant nature of fertility restoration, which is controlled by one or two major genes with modifiers. In this study, the genetics of fertility restoration of the A1 cytoplasmic nuclear male sterility system (CGMS) in sorghum were investigated in segregating F2 and BC1 populations of A1 cytoplasm crosses. Fertility restoration was governed by a monogenic inheritance (3F:1S) mechanism represented by a single dominant gene responsible for fertility restoration in all of the crosses studied.


title: Influence of residual microbial nucleic acid in instruments for hip joint arthroplasty on false-positive results of metagenomic sequencing

abstract: The purpose of this study was to analyze the effect of residual microbial nucleic acid in instruments for hip joint arthroplasty on false-positive sequencing results. Samples were taken from 3 different acetabular reamer for hip arthroplasty in 7 different hospitals. The whole process was strictly aseptic, metagenomic next-generation sequencing (mNGS) was performed according to standard operating procedures. The sterility of instruments was confirmed by culture method. The sequencing results of specimens from different hospitals were compared to analyze the difference of background bacteria. Bioinformatics analysis and visualization were presented through R language. A total of 26 samples were processed by mNGS, including 24 instrument swab samples, 1 blank swab control, and 1 blank water control. 254,314,707 reads were sequenced in all samples. The results showed that 1.13% of Clean Reads can be matched to pathogenic microorganism genomes, of which bacterial sequences account for 87.48%, fungal sequences account for 11.18%, parasite sequences account for 1.26%, and virus sequences account for 0.06%. The results of PCA (Principal Component Analysis) demonstrated that the distribution of bacteria on the surface of instruments was significantly different between medical institutions. Through the Venn diagram, it was found that 465 species of bacteria in all region hospitals, Liaocheng Peoples Hospital had a maximum of 340 species of bacteria, followed by Guanxian County Peoples Hospital with 169 species. The clustering heat map illustrated that the distribution of bacterial groups in three different instrument samples in the same hospital was basically the same, and the bacterial genera varied significantly among hospitals. The residual microbial nucleic acid fragments are mainly bacterial DNA and represent differences in different medical institutions, The establishment of independent background bacterial libraries in different medical institutions can effectively improve the accuracy of mNGS diagnosis and help to exclude background microorganisms interference.

category: genomics


title: The DNA-based global positioning system--a theoretical framework for large-scale spatial genomics

abstract: We present GPS-seq, a theoretical framework that enables massively scalable, optics-free spatial transcriptomics. GPS-seq combines data from high-throughput sequencing with manifold learning to obtain the spatial transcriptomic landscape of a given tissue section without optical microscopy. In this framework, similar to technologies like Slide-seq and 10X Visium, tissue samples are stamped on a surface of randomly-distributed DNA-barcoded spots (or beads). The transcriptomic sequences of proximal cells are fused to DNA barcodes, enabling the recovery of a transcriptomic pixel image by high-throughput sequencing. The barcode spots serve as "anchors" which also capture spatially diffused "satellite" barcodes, and therefore allow computational reconstruction of spot positions without optical sequencing or depositing barcodes to pre-specified positions. In theory, it could generate 100 mm x 100 mm spatial transcriptomic images with 10-20 m resolution by localizing 108 DNA-barcoded pixels with a single Illumina NovaSeq run. The general framework of GPS-seq is also compatible with standard single-cell (or single-nucleus) capture methods, and any modality of single-cell genomics, such as sci-ATAC-seq, could be transformed into spatial genomics in this strategy. We envision that GPS-seq will lead to breakthrough discoveries in diverse areas of biology by enabling organ-scale imaging of multiple genomic statuses at single-cell resolution for the first time.


title: Large-bodied ornithomimosaurs inhabited Appalachia during the Late Cretaceous of North America

abstract: Reconstructing the evolution, diversity, and paleobiogeography of North Americas Late Cretaceous dinosaur assemblages requires spatiotemporally contiguous data; however, there remains a spatial and temporal disparity in dinosaur data on the continent. The rarity of vertebrate-bearing sedimentary deposits representing Turonian-Santonian ecosystems, and the relatively sparse record of dinosaurs from the eastern portion of the continent, present persistent challenges for studies of North American dinosaur evolution. Here we describe an assemblage of ornithomimosaurian materials from the Santonian Eutaw Formation of Mississippi. Morphological data coupled with osteohistological growth markers suggest the presence of two taxa of different body sizes, including one of the largest ornithomimosaurians known worldwide. The regression predicts a femoral circumference and a body mass of the Eutaw individuals similar to or greater than that of large-bodied ornithomimosaurs, Beishanlong grandis and Gallimimus bullatus. The paleohistology of MMNS VP-6332 demonstrates that the individual was at least 11 years of age (similar to B. grandis [[~]375 kg, 13-14 years old at death]). Additional pedal elements share some intriguing features with ornithomimosaurs yet suggest a larger-body size closer to Deinocheirus mirificus. The presence of a large-bodied ornithomimosaur in this region during this time is consistent with the relatively recent discoveries of early-diverging, large-bodied ornithomimosaurs from mid-Cretaceous strata of Laurasia (Arkansaurus fridayi and B. grandis). The smaller Eutaw taxon is represented by a tibia preserving seven growth cycles, with osteohistological indicators of decreasing growth, yet belongs to an individual with near reaching somatic maturity of the larger taxon, suggesting the co-existence of medium- and large-bodied ornithomimosaur taxa during the Late Cretaceous Santonian of North America. The Eutaw ornithomimosaur materials provide key information on the diversity and distribution of North American ornithomimosaurs and Appalachian dinosaurs and fit with broader evidence of multiple cohabiting species of ornithomimosaurian dinosaurs in Late Cretaceous ecosystems of Laurasia.

category: paleontology


title: Evidence for the role of music on the growth and signal response in duckweed

abstract: Sound vibration, an external mechanical force, has been proved to modulate plant growth and development like rain, wind, and vibration. However, the role of music on plants, especially on signal response, has been usually neglected in research. Herein, we investigated the growth state, gene expression, and signal response in duckweed treated with soft music. The protein content in duckweed after music treatment for 7 days was about 1.6 times that in duckweed without music treatment. Additionally, Fv/Fm in duckweed treated with music was 0.78, which was significantly higher in comparison with the control group (P < 0.01). Interestingly, music promoted the Glu and Ca signaling response. To further explore the global molecular mechanism, we performed transcriptome analysis. A total of 1296 DEGs were found for all these investigated genes in duckweed treated with music compared to the control group. Among these, up-regulation of the expression of metabolism-related genes related to glycolysis, cell wall biosynthesis, oxidative phosphorylation, and pentose phosphate pathways were found. Overall, these results provided a molecular basis to music-triggered signal response, transcriptomic, and growth changes in duckweed, which also highlighted the potential of music as an environmentally friendly stimulus to promote improved protein production in duckweed.

category: plant biology


title: Differential Hebbian learning with time-continuous signals for active noise reduction

abstract: Spike timing-dependent plasticity, related to differential Hebb-rules, has become a leading paradigm in neuronal learning, because weights can grow or shrink depending on the timing of pre- and post-synaptic signals. Here we use this paradigm to reduce unwanted (acoustic) noise. Our system relies on heterosynaptic differential Hebbian learning and we show that it can efficiently eliminate noise by up to -140 dB in multi-microphone setups under various conditions. The system quickly learns, most often within a few seconds, and it is robust with respect to different geometrical microphone configurations, too. Hence, this theoretical study demonstrates that it is possible to successfully transfer differential Hebbian learning, derived from the neurosciences, into a technical domain.


title: Cardio PyMEA: A user-friendly, open-source Python application for cardiomyocyte microelectrode array analysis

abstract: Open source analytical software for the analysis of electrophysiological cardiomyocyte data offers a variety of new functionalities to complement closed-source, proprietary tools. Here, we present the Cardio PyMEA application, a free, modifiable, and open source program for the analysis of microelectrode array (MEA) data obtained from cardiomyocyte cultures. Cardio PyMEA was written entirely in Python 3 to provide an accessible, integrated workflow that possesses a user-friendly graphical user interface (GUI) written in PyQt5 to allow for performant, cross-platform utilization. This application makes use of object-oriented programming (OOP) principles to facilitate the relatively straightforward incorporation of custom functionalities, e.g. power law analysis, that suit the needs of the user. Cardio PyMEA is available as an open source application under the terms of the GNU General Public License (GPL). The source code for Cardio PyMEA can be downloaded from Github at the following repository: https://github.com/csdunhamUC/cardio_pymea.

category: cell biology


title: The Graphical R2D2 Estimator for the Precision Matrices

abstract: Biological networks are important for the analysis of human diseases, which summarize the regulatory interactions and other relationships between different molecules. Understanding and constructing networks for molecules, such as DNA, RNA and proteins, can help elucidate the mechanisms of complex biological systems. The Gaussian Graphical Models (GGMs) are popular tools for the estimation of gene regulatory networks because of their biological interpretability. Nonetheless, reconstructing GGMs from high-dimensional datasets is still challenging and current methods cannot handle the sparsity and high-dimensionality issues arising from datasets very well. Here we developed a new GGM, called the graphical R2D2 (R2-induced Dirichlet Decomposition), based on the R2D2 priors for linear models. The proposed graphical R2D2 estimator has attractive properties in estimating the precision matrices, such as greater concentration near the origin and heavier tails than current shrinkage priors. When the true precision matrix is sparse and of high dimension, the graphical R2D2 provides the estimates with smallest information divergence from the sampling model. Besides, we also provide a full Gibbs sampler for implementing the graphical R2D2 estimator. R code is available for download from https://github.com/RavenGan/GR2D2. The graphical R2D2 estimator shows superior performance to existing techniques in simulations. We also provide breast cancer gene network analysis using the graphical R2D2 estimator and the important genes recognized from the inferred gene regulatory networks are consistent with biological ground truth.


title: Examination of genome-wide ortholog variation in clinical and environmental isolates of the fungal pathogen Aspergillus fumigatus

abstract:
Aspergillus fumigatus is both an environmental saprobe and an opportunistic human fungal pathogen. Knowledge of genomic variation across A. fumigatus isolates is essential for understanding the evolution of pathogenicity, virulence, and resistance to antifungal drugs. Here, we investigated 206 A. fumigatus isolates (133 clinical and 73 environmental isolates) aiming to identify genes with variable presence across isolates and test whether this variation was related to the clinical or environmental origin of isolates. The PanCore genome of A. fumigatus constitutes 13,085 ortholog groups, of which 7,773 (59.4%) are shared by all isolates (CORE) and 5,312 (40.6%) vary in their gene presence across isolates (ACCESSORY). Despite differences in the distribution of orthologs across all isolates, no significant differences were observed among clinical vs. environmental isolates when accounting for phylogeny. Orthologs that differ in their distribution across isolates tend to occur in low frequency and/or be restricted to specific isolates; thus, the degree of genomic conservation between orthologs of A. fumigatus is high. These results suggest that differences in the distribution of orthologs within A. fumigatus cannot be associated with the clinical or environmental origin of isolates.

ImportanceAspergillus fumigatus is a cosmopolitan species of fungi responsible for thousands of cases of invasive disease. Clinical and environmental isolates of A. fumigatus exhibit extensive phenotypic differences, including differences related to virulence and antifungal drug resistance. A comprehensive survey of the genomic diversity present in A. fumigatus and its relationship to the clinical or environmental origin of isolates can contribute to the prediction of the mechanisms of evolution and infection of the species. Our results suggest that there is no significant variation in ortholog distribution between clinical and environmental isolates when accounting for evolutionary history. The work supports the hypothesis that environmental and clinical isolates of A. fumigatus do not differ in their gene contents.


title: NAb-seq: an accurate, rapid and cost-effective method for antibody long-read sequencing in hybridoma cell lines and single B cells

abstract: Despite their common use in research, monoclonal antibodies are currently not systematically sequenced. This can lead to issues with reproducibility and the occasional loss of antibodies with loss of cell lines. Hybridoma cell lines have been the primary means of generating monoclonal antibodies from immunized animals including mice, rats, rabbits and alpacas. Excluding therapeutic antibodies, few hybridoma-derived antibody sequences are known. Sanger sequencing has been "the gold standard" for antibody gene sequencing but relies on the availability of species-specific degenerate primer sets for amplification of light and heavy antibody genes, in addition to lengthy and expensive cDNA preparation. Here we leveraged recent improvements in long-read Oxford Nanopore Technologies (ONT) sequencing to develop NAb-seq: a three-day, species-independent, and cost-effective workflow to characterize paired full-length immunoglobulin light and heavy chain genes from hybridoma cell lines. When compared to Sanger sequencing of two hybridoma cell lines, long-read ONT sequencing was highly accurate, reliable, and amenable to high throughput. We further show that the method is applicable to single cells, allowing efficient antibody discovery in rare populations such as memory B cells. In summary, NAb-seq promises to accelerate identification and validation of hybridoma antibodies as well as antibodies from single B cells used in research, diagnostics and therapeutics.


title: Rad51 is a druggable target that sustains replication fork progression upon DNA replication stress

abstract: Solving the problems that replication forks encounter when synthesizing DNA is essential to prevent genomic instability. Besides their role in DNA repair in the G2 phase, several homologous recombination proteins, specifically Rad51, have prominent roles in the S phase. Using different cellular models, Rad51 has been shown not only to be present at ongoing and arrested replication forks but also to be involved in nascent DNA protection and replication fork restart. Through pharmacological inhibition, here we study the specific role of Rad51 in the S phase. Rad51 inhibition in non-transformed cell lines did not have a major effect on replication fork progression under non-perturbed conditions, but when the same cells were subjected to replication stress, Rad51 became necessary to maintain replication fork progression. Notably, the inhibition or depletion of Rad51 did not compromise fork integrity when subjected to hydroxyurea treatment. Rad51 inhibition also did not decrease the ability to restart, but rather compromised, fork progression during reinitiation. In agreement with the presence of basal replication stress in human colorectal cancer cells, Rad51 inhibition reduced replication fork speed in these cells and increased {gamma}H2Ax foci under control conditions. These alterations could have resulted from the reduced association of DNA polymerase  to chromatin, as observed when inhibiting Rad51. It may be possible to exploit the differential dependence of non-transformed cells versus colorectal cancer cells on Rad51 activity under basal conditions to design new therapies that specifically target cancer cells.


title: Neural responses to biological motion distinguish autistic and schizotypal traits in the general population

abstract: Difficulties in social interactions are common to both autism and schizophrenia, and contribute to correlated autistic and schizotypal traits in the neurotypical population. It remains unresolved whether this represents a shared etiology or a superficial phenotypic overlap. Both conditions are associated with atypical neural activity in response to the perception of social stimuli, and also decreased neural synchronization between individuals that may prohibit establishing shared experiences. This study sought to establish if neural activity and neural synchronization associated with biological motion perception are differentially associated with autistic and schizotypal traits in the neurotypical population. Participants watched an audiovisual montage of naturalistic social interactions whilst hemodynamic brain activity was measured with fMRI. A separate sample of participants provided a continuous measure of the extent of biological motion, which was used to predict hemodynamic activity. General Linear Model analysis revealed that biological motion perception was associated with neural activity across the action-observation network. However, inter-subject phase synchronization analysis revealed that neural activity synchronized between individuals in occipital and parietal areas, but de-synchronized in temporal and frontal regions. Autistic traits were associated with a decrease in neural activity (precuneus, middle cingulate gyrus) and schizotypal traits were associated with a decrease in neural synchronization (middle and inferior frontal gyri). Biological motion perception elicits convergent and divergent patterns of neural activity and neural synchronization, and are differentially associated with shared traits related with autism and schizophrenia in the general population, suggesting that they originate from different neural mechanisms.


title: Individual behavioral trajectories shape whole-brain connectivity in mice

abstract: It is widely assumed that our actions shape our brains and that the resulting connections determine who we are. To test this idea in a reductionist setting, in which genes and environment are controlled, we investigated differences in neuroanatomy and structural covariance by ex vivo structural magnetic resonance imaging (MRI) in mice whose behavioral activity was continuously tracked for 3 months in a large, enriched environment. We confirmed that environmental enrichment increases mouse hippocampal volumes. Stratifying the enriched group according to individual longitudinal behavioral trajectories, however, revealed striking differences in mouse brain structural covariance in continuously highly active mice compared to those whose trajectories showed signs of habituating activity. Network-based statistics identified distinct sub-networks of murine structural covariance underlying these differences in behavioral activity. Together, these results reveal that differentiated behavioral trajectories of mice in an enriched environment are associated with differences in brain connectivity.


title: Successor-like representation guides the prediction of future events in human visual cortex and hippocampus

abstract: Human agents build models of their environment, which enable them to anticipate and plan upcoming events. However, little is known about the properties of such predictive models. Recently, it has been proposed that hippocampal representations take the form of a predictive map-like structure, the so-called successor representation. Here we used human fMRI to probe whether activity in the early visual cortex (V1) and hippocampus adhere to the postulated properties of the successor representation after visual sequence learning. Participants were exposed to an arbitrary spatiotemporal sequence consisting of four items (A-B-C-D). We found that after repeated exposure to the sequence, merely presenting single sequence items (e.g., - B - -) resulted in V1 activation at the successor locations of the full sequence (e.g., C-D), but not at the predecessor locations (e.g., A). This highlights that visual representations are skewed toward future states, in line with the successor representation. Similar results were also found in the hippocampus. Moreover, the hippocampus developed a tuning profile that showed sensitivity to the temporal distance in sequence-space, with fading representations for sequence events in the more distant past and future. V1, in contrast, showed a tuning profile that was only sensitive to spatial distance in stimulus-space. Together, these results provide empirical evidence for the proposition that both visual and hippocampal cortex represent a predictive map of the visual world akin to the successor representation.


title: Neural Correlates of Linguistic Collocations During Continuous Speech Perception

abstract: Language is fundamentally predictable, both on a higher schematic level as well as low-level lexical items. Regarding predictability on a lexical level, collocations are frequent co-occurrences of words that are often characterized by high strength of association. So far, psycho-and neurolin guistic studies have mostly employed highly artificial experimental paradigms in the investigation of collocations by focusing on the processing of single words or isolated sentences. In contrast, here we analyze EEG brain responses recorded during stimulation with continuous speech, i.e. audio books. We find that the N400 response to collocations is significantly different from that of non-collocations, whereas the effect varies with respect to cortical region (anterior/ posterior) and laterality (left/right). Our results are in line with studies using continuous speech, and they mostly contradict those using artificial paradigms and stimuli. To the best of our knowledge, this is the first neurolinguistic study on collocations using continuous speech stimulation.


title: At-home blood collection and stabilization in high temperature climates using homeRNA

abstract: Expanding whole blood sample collection for transcriptome analysis beyond traditional phlebotomy clinics will open new frontiers for remote immune research and telemedicine. Determining the stability of RNA in blood samples exposed to high ambient temperatures (>30{degrees}C) is necessary for deploying home-sampling in settings with elevated temperatures (e.g., studying physiological response to natural disasters that occur in warm locations or in the summer). Recently, we have developed homeRNA, a technology that allows for self-blood sampling and RNA stabilization remotely. homeRNA consists of a lancet-based blood collection device, the Tasso-SST which collects up to 0.5 mL of blood from the upper arm, and a custom-built stabilization transfer tube containing RNAlater. In this study, we investigated the robustness of our homeRNA kit in high temperature settings via two small pilot studies in Doha, Qatar (no. participants = 8), and the Western USA during the summer of 2021, which included a heatwave of unusually high temperatures in some locations (no. participants = 11). Samples collected from participants in Doha were subjected to rapid external temperature fluctuations from being moved to and from air-conditioned areas and extreme heat environments (up to 41{degrees}C external temperature during brief temperature spikes). In the USA pilot study, regions varied in outdoor temperature highs (between 25{degrees}C and 43.4{degrees}C). All samples that returned a RNA integrity number (RIN) value from the Doha, Qatar group had a RIN [≥]7.0, a typical integrity threshold for downstream transcriptomics analysis. RIN values for the Western USA samples (n=12 samples) ranged from 6.9-8.7 with 9 out of 12 samples reporting RINs [≥]7.0. Overall, our pilot data suggest that homeRNA can be used in some regions that experience elevated temperatures, opening up new geographical frontiers in disseminated transcriptome analysis for applications critical to telemedicine, global health, and expanded clinical research. Further studies, including our ongoing work in Qatar, USA, and Thailand, will continue to test the robustness of homeRNA.

category: bioengineering


title: Multifactorial remodeling of the cancer immunopeptidome by interferon gamma

abstract: IFN{gamma} alters the immunopeptidome presented on HLA class I (HLA-I), and its activity on cancer cells is known to be important for effective immunotherapy responses. We performed proteomic analyses of untreated and IFN{gamma}-treated colorectal cancer patient-derived organoids (PDOs) and combined this with transcriptomic and HLA-I immunopeptidomics data to dissect mechanisms that lead to remodeling of the immunopeptidome through IFN{gamma}. IFN{gamma}-induced changes in the abundance of source proteins, switching from the constitutive- to the immunoproteasome, and differential upregulation of different HLA alleles explained some, but not all, observed peptide abundance changes. By selecting for peptides which increased or decreased the most in abundance, but originated from proteins with limited abundance changes, we discovered that the amino acid composition of presented peptides also influences whether a peptide is up- or downregulated on HLA-I through IFN{gamma}. The presence of proline within the peptide core was most strongly associated with peptide downregulation. This was validated in an independent dataset. Proline substitution in relevant core positions did not influence the predicted HLA-I binding affinity or stability, indicating that proline effects on peptide processing may be most relevant. Understanding the multiple factors that influence the abundance of peptides presented on HLA-I in the absence or presence of IFN{gamma} is important to identify the best targets for antigen-specific cancer immunotherapies such as vaccines or T-cell receptor engineered therapeutics.

category: cancer biology


title: Keratose sponges in ancient carbonates - a problem of interpretation

abstract: Increasing current interest in sponge fossils includes numerous reports of diverse vermicular and peloidal structures interpreted as keratose sponges in Neoproterozoic to Mesozoic carbonates and in various open marine to peritidal and restricted settings. Reports of their occurrence are fundamental and far-reaching for understanding microfacies and diagenesis where they occur; and fossil biotic assemblages, as well as wider aspects of origins of animals, sponge evolution/ecology and the systemic recovery from mass extinctions. Keratose sponges: 1) have elaborate spongin skeletons but no spicules, thus lack mineral parts and therefore have poor preservation potential so that determining their presence in rocks requires interpretation; and 2) are presented in publications as interpreted fossil structures almost entirely in two-dimensional (thin section) studies, where structures claimed as sponges comprise diverse layered, network, particulate and amalgamated fabrics involving calcite sparite in a micritic groundmass. There is no verification of sponges in these cases and almost all of them can be otherwise explained; some are certainly not correctly identified. The diversity of structures seen in thin sections may be reinterpreted to include: a) meiofaunal activity; b) layered, possibly microbial (spongiostromate) accretion; c) sedimentary peloidal to clotted micrites; d) fluid escape and capture resulting in birdseye to vuggy porosities; and e) molds of siliceous sponge spicules. Without confirmation of keratose sponges in ancient carbonates, interpretations of their role in ancient carbonate systems, including facies directly after mass extinctions, are unsafe, and alternative explanations for such structures should be considered. This study calls for greater critical appraisal of evidence, to seek confirmation or not, of keratose sponge presence. (259/300 max, for Sedimentology)


title: Hypomorphism of a novel long ERα isoform causes severe reproductive dysfunctions in female mice

abstract: Estrogen receptor alpha (ER)-mediated estrogen signaling play a pivotal role in both reproductive and non-reproductive functions. Transcriptional regulation of ER gene is highly complex, with multiple transcript variants being differentially produced across the tissues. However, tissue-specific variation and physiological specificity of the ER variants are not yet fully understood. In an attempt to generate a Cre-dependently restorable ER-null mice for functional investigation of the genetic sufficiency, we unexpectedly produced ER hypomorphic mice with biased downregulation of a previously unappreciated long ER isoform that is enriched in the female reproductive organs (uterus and ovaries) and the pituitary but minimally expressed in the brain. Female homozygous mutant mice were capable of pregnancy but displayed irregular estrus cycle and rarely maintained alive newborns without significant morphological and pathological changes in reproductive system and the disruption of body weight homeostasis, indicating the vital role of this long isoform in female reproductive function. Collectively, our results define a tissue-specifically enriched long ER isoform and its preferential role in female reproductive function over body weight homeostasis.

category: physiology


title: A new monoclonal antibody enables BAR analysis of subcellular importin β1 interactomes

abstract: Importin {beta}1 (KPNB1) is a nucleocytoplasmic transport factor with critical roles in both cytoplasmic and nucleocytoplasmic transport, hence there is keen interest in the characterization of its subcellular interactomes. We found limited efficiency of BioID in detection of importin complex cargos, and therefore generated a highly specific and sensitive anti-KPNB1 monoclonal antibody to enable Biotinylation by Antibody Recognition (BAR) analysis of importin {beta}1 interactomes. The monoclonal antibody recognizes an epitope comprising residues 301-320 of human KPBN1, and strikingly is highly specific for cytoplasmic KPNB1 in diverse applications, with little or no reaction with KPNB1 in the nucleus. BAR with this novel antibody revealed numerous new interactors of importin {beta}1, expanding the KPNB1 interactome to cytoplasmic and signaling complexes that highlight potential new functions for the importins complex beyond nucleocytoplasmic transport. Data are available via ProteomeXchange with identifier PXD032728.


title: RNA-targeting CRISPR-Cas13 Provides Broad-spectrum Phage Immunity

abstract: CRISPR-Cas13 proteins are RNA-guided RNA nucleases that defend against invasive phages through general, non-specific RNA degradation upon complementary target transcript binding. Despite being RNA nucleases, Cas13 effectors are capable of inhibiting the infection of dsDNA phages but have only been investigated across a relatively small sampling of phage diversity. Here, we employ a systematic, phage-centric approach to determine the anti-phage capacity of Cas13 and find LbuCas13a to be a remarkably potent phage inhibitor. LbuCas13a confers robust, consistent antiviral activity regardless of gene essentiality, gene expression timing or target sequence location. Furthermore, after challenging LbuCas13a with eight diverse E. coli phages distributed across E. coli phage phylogenetic groups, we find no apparent phage-encoded limits to its potent antiviral activity. In contrast to other Class 2 CRISPR-Cas proteins, these results suggest that DNA phages are generally vulnerable to Cas13a targeting. Leveraging this effective anti-phage activity, LbuCas13a can be used seamlessly as a counter-selection agent for broad-spectrum phage editing. Using a two-step phage editing and enrichment approach, we show that LbuCas13a enables markerless genome edits in phages with exceptionally high efficiency and precision, including edits as small as a single codon. By taking advantage of the broad vulnerability of RNA during viral infection, Cas13a enables a generalizable strategy for editing the most abundant and diverse biological entities on Earth.


title: Conserved GYXLI motif of FlhA is involved in dynamic domain motions of FlhA required for flagellar protein export

abstract:
Flagellar structural subunits are transported via the flagellar type III secretion system (fT3SS) and assemble at the distal end of the growing flagellar structure. The C-terminal cytoplasmic domain of FlhA (FlhAC) serves as a docking platform for export substrates and flagellar chaperones and plays an important role in hierarchical protein targeting and export. FlhAC consists of domains D1, D2, D3, and D4 and adopts open and closed conformations. Gly-368 of Salmonella FlhA is located within the highly conserved GYXLI motif and is critical for the dynamic domain motions of FlhAC. However, it remains unclear how it works. Here, we report that periodic conformational changes of the GYXLI motif induce a remodeling of hydrophobic side-chain interaction networks in FlhAC and promotes the cyclic open-close domain motions of FlhAC. The temperature-sensitive flhA(G368C) mutation stabilized a completely closed conformation at 42{degrees}C through strong hydrophobic interactions between Gln-498 of domain D1 and Pro-667 of domain D4 and between Phe-459 of domain D2 and Pro-646 of domain D4, thereby inhibiting flagellar protein export by the fT3SS. Its intragenic suppressor mutations reorganized the hydrophobic interaction networks in the closed FlhAC structure, restoring the protein export activity of the fT3SS to a significant degree. Furthermore, the conformational flexibility of the GYXLI motif was critical for flagellar protein export. We propose that the conserved GYXLI motif acts as a structural switch to induce the dynamic domain motions of FlhAC required for efficient and rapid protein export by the fT3SS.

IMPORTANCEMany motile bacteria employ the flagellar type III secretion system (fT3SS) to construct flagella beyond the cytoplasmic membrane. The C-terminal cytoplasmic domain of FlhA (FlhAC), a transmembrane subunit of the fT3SS, provides binding-sites for export substrates and flagellar export chaperones to coordinate flagellar protein export with assembly. FlhAC undergoes cyclic open-close domain motions. The highly conserved Gly-368 residue of FlhA is postulated to be critical for dynamic domain motions of FlhAC. However, it remains unknown how it works. Here, we carried out mutational analysis of FlhAC combined with molecular dynamics simulation and provide evidence that the conformational flexibility of FlhAC by Gly-368 is important for remodeling hydrophobic side-chain interaction networks in FlhAC to facilitate its cyclic open-close domain motions, allowing the fT3SS to transport flagellar structural subunits for efficient and rapid flagellar assembly.


title: Candida albicans filamentation does not require the cAMP-PKA pathway in vivo

abstract:
Candida albicans is one of the most prevalent human fungal pathogens. Its ability to transition between budding yeast and filamentous morphological forms (pseudohyphae and hyphae) is tightly associated with its pathogenesis. Based on in vitro studies, the cAMP-Protein Kinase A (PKA) pathway is a key regulator of C. albicans morphogenesis. Using an intravital imaging approach, we investigated the role of the cAMP-PKA pathway during infection. Consistent with their roles in vitro, the downstream effectors of the cAMP-PKA pathway Efg1 and Nrg1 function, respectively, as an activator and a repressor of in vivo filamentation. Surprisingly, strains lacking the adenylyl cyclase, CYR1, showed only slightly reduced filamentation in vivo despite being completely unable to filament in RPMI+10% serum at 37{degrees}C. Consistent with these findings, deletion of the catalytic subunits of PKA (Tpk1 and Tpk2), either singly or in combination, generated strains that also filamented in vivo but not in vitro. In vivo transcription profiling of C. albicans isolated from both ear and kidney tissue showed that the expression of a set of 184 environmentally responsive correlated well with in vitro filamentation (R2 0.62-0.68) genes. This concordance suggests that the in vivo and in vitro transcriptional responses are similar but that the upstream regulatory mechanisms are distinct. As such, these data emphatically emphasize that C. albicans filamentation is a complex phenotype that occurs in different environments through an intricate network of distinct regulatory mechanisms.

ImportanceThe fungus Candida albicans causes a wide range of disease in humans from common diaper rash to life-threatening infections in patients with compromised immune systems. As such, the mechanisms for its ability to cause disease are of wide interest. An intensely studied virulence property of C. albicans is its ability to switch from a round yeast form to filament-like forms (hyphae and pseudohyphae). Surprisingly, we have found that a key signaling pathway that regulates this transition in vitro, the protein kinase A pathway, is not required for filamentation during infection of the host. Our work not only demonstrates that the regulation of filamentation depends upon the specific environment C. albicans inhabits but also underscores the importance of studying these mechanisms during infection.


title: Transcriptome of epibiont Saccharibacteria Nanosynbacter lyticus strain TM7x during establishment of symbiosis

abstract:
Saccharibacteria Nanosynbacter lyticus strain TM7x is a member of the broadly distributed Candidate Phylum Radiation. These bacteria have ultrasmall cell size, reduced genomes and live as epibionts on the surface of other bacteria. The mechanisms by which they establish and maintain this relationship are not yet fully understood. The transcriptomes of the epibiont TM7x and its host bacteria Schaalia odontolyticus strain XH001 were captured across the establishment of symbiosis during both the initial interaction and stable symbiosis. The results showed a dynamic interaction with large shifts in gene expression for both species between the initial encounter and stable symbiosis, notably transporter genes. During stable symbiosis, the host XH001 showed higher gene expression for peptidoglycan biosynthesis, mannosylation, cell cycle and stress related genes, but lower expression of chromosomal partitioning genes. This was consistent with the elongated cell shape seen in XH001 infected with TM7x and our discovery that infection resulted in thickened cell walls. Within TM7x, increased pili, type IV effector gene, and arginine catabolism/biosynthesis gene expression during stable symbiosis implied a key role for these functions in the interaction. Consistent with its survival and persistence in the human microbiome as an obligate epibiont with reduced de novo biosynthetic capacities, TM7x also showed higher levels for energy production and peptidoglycan biosynthesis but lower expression of stress related genes during stable symbiosis. These results imply that TM7x and its host bacteria keep a delicate balance in order to sustain an episymbiotic lifestyle.

IMPORTANCENanosynbacter lyticus type strain TM7x is the first cultivated member of the Saccharibacteria and the Candidate Phyla Radiation (CPR). It was discovered to have ultrasmall cell size with a highly reduced genome that establishes an obligate epibiotic relationship with its host bacterium. The CPR, now formally proposed as the Patescibacteria super-phylum, is a large monophyletic radiation of diverse bacteria with reduced genomes that includes Saccharibacteria. The vast majority of the CPR have yet to be cultivated in the laboratory and our insights into these unique organisms to date has been derived from only a few Saccharibacteria species. It is unknown however how these small obligate parasitic Saccharibacteria, that are missing many de novo biosynthetic pathways, are maintained at high prevalence within the human microbiome as well as in the environment. When TM7x infects its host bacterium there are distinct temporal phases, including an initial interaction, a killing phase, recovery phase, and finally stable symbiosis. Here we captured the gene expression of the host bacterium and epibiont during this dynamic interaction which represents the initial insights into the mechanisms of how these unique microbes may survive and persist.


title: Comparative Genome Analysis of Lactobacillus acidophilus Isolates from Different Ecological Niches

abstract: Lactobacillus acidophilus has been extensively applied in plentiful probiotic products and is mostly found in the gastrointestinal tract, vagina, and oral cavity of human and animal, and fermented foods. Although several studies have been performed to investigate the beneficial characteristics and genome function of L. acidophilus, comparative genomic analysis remains scarce. In this study, we collected 74 L. acidophilus genomes from our gut bacterial genome collection and the public database and conducted a comprehensive comparative genomic analysis. The analysis of average nucleotide identity (ANI), phylogenetic, gene distribution of COG and KEGG database, carbohydrates utilization, and secondary metabolites revealed the potential correlation of the genomic diversity and niches adaptation of L. acidophilus from different perspectives. In addition, the pan-genome of L. acidophilus was found to be open, with metabolism, information storage and processing genes mainly distributed in the core genome. Phage- and peptidase-associated genes were found in the genome of the specificity of animal-derived strains, which were related to adaptation of animal gut. SNP analysis showed the differences of the utilization of Vitamin B12 in cellular of L. acidophilus strains from animal gut and others. This work provides new insights for the genomic diversity analysis of Lactobacillus acidophilus and uncovers the ecological adaptation of the specific strains.


title: Structure-guided synthesis of FK506 and FK520 analogs with increased selectivity exhibit in vivo therapeutic efficacy against C. neoformans

abstract:
Calcineurin is an essential virulence factor that is conserved across human fungal pathogens including Cryptococcus neoformans, Aspergillus fumigatus, and Candida albicans. Although an excellent target for antifungal drug development, the serine-threonine phosphatase activity of calcineurin is conserved in mammals and inhibition of this activity results in immunosuppression. FK506 (tacrolimus) is a naturally produced macrocyclic compound that inhibits calcineurin by binding to the immunophilin FKBP12. Previously, our fungal calcineurin-FK506-FKBP12 structure-based approaches identified a non-conserved region of FKBP12 that can be exploited for fungal-specific targeting. These studies led to the design of an FK506 analog, APX879, modified at the C22 position that was less immunosuppressive yet maintained antifungal activity. We now report high resolution protein crystal structures of fungal FKBP12 and a human truncated calcineurin-FKBP12 bound to a natural FK506 analog, FK520 (ascomycin). Based on information from these structures and the success of APX879, we synthesized and screened a novel panel of C22-modified compounds derived from both FK506 and FK520. One compound, JH-FK-05 demonstrates broad-spectrum antifungal activity in vitro and is non-immunosuppressive in vivo. In murine models of pulmonary and disseminated C. neoformans infection, JH-FK-05 treatment significantly reduced fungal burden and extended animal survival alone and in combination with fluconazole. Furthermore, molecular dynamic simulations performed with JH-FK-05 binding to fungal and human FKBP12 identified additional residues outside of the C22 and C21 positions that could be modified to generate novel FK506 analogs with improved antifungal activity.

SignificanceDue to rising rates of antifungal drug resistance and a limited armamentarium of antifungal treatments, there is a paramount need for novel antifungal drugs to treat systemic fungal infections. Calcineurin has been established as an essential and conserved virulence factor in several fungi, making it an attractive antifungal target. However, due to the immunosuppressive action of calcineurin inhibitors, they have not been successfully utilized clinically for antifungal treatment in humans. Recent availability of crystal structures of fungal calcineurin bound inhibitor complexes have enabled the structure-guided design of FK506 analogs and led to a breakthrough in the development of a compound increased for fungal specificity. The development of a calcineurin inhibitor with reduced immunosuppressive activity and therapeutic antifungal activity would add a significant tool to the treatment options for these invasive fungal infections with exceedingly high rates of mortality.


title: Structure-guided synthesis of FK506 and FK520 analogs with increased selectivity exhibit in vivo therapeutic efficacy against Cryptococcus


title: Drug design and repurposing with a sequence-to-drug paradigm

abstract: Drug development based on target proteins has been a successful approach in recent decades. A conventional structure-based drug design pipeline is a complex, human-engineered pipeline with multiple independently optimized steps. Advances in end-to-end differentiable learning suggest the potential benefits of similarly reformulating drug design. Here, we proposed a new sequence-to-drug paradigm that discovers drug-like small-molecule modulators directly from protein sequences and validated this concept for the first time in three stages. First, we designed TransformerCPI2.0 as a core tool for the sequence-to-drug paradigm, which exhibited competitive performance with conventional structure-based drug design approaches. Second, we validated the binding knowledge that TransformerCPI2.0 has learned. Third, we applied a sequence-to-drug paradigm to discover new hits for E3 ubiquitin-protein ligases: speckle-type POZ protein (SPOP), ring finger protein 130 (RNF130) which does not have a 3D structure, and repurposed proton pump inhibitors (PPIs) for ADP-ribosylation factor 1 (ARF1). This first proof of concept shows that the sequence-to-drug paradigm is a promising direction for drug development.


title: A Conformation-specific Nanobody Targeting the NMN-activated State of SARM1

abstract: Upon axonal injury, Sterile alpha (SAM) and Toll/interleukin-1 receptor (TIR) motif containing 1 (SARM1) is activated by nicotinamide mononucleotide (NMN) to deplete NAD and consequently promote the process of axon degeneration (AxD). Currently, only the inactive form of SARM1 in its auto-inhibitory conformation has been resolved. The flexibility of the enzymatically active form of SARM1 has so far precluded its structural determination. To solve the problem, we generated a stabilizing nanobody, Nb-C6, that specifically recognized only the NMN-activated form of SARM1. The conformation specificity was verified by immunoprecipitation and surface plasmon resonance. Fluorescently labeled Nb-C6 could immunostain only the activated SARM1 in cells stimulated with CZ-48, a permeant mimetic of NMN. Expression of Nb-C6 in live cells resulted in stabilization of the active form of the endogenous and exogenous SARM1, producing and elevating cellular levels of cyclic ADP-ribose, a calcium messenger. Cryo-EM of the NMN-activated SARM1 complexed by Nb-C6 showed an octameric structure resembling a "blooming lotus" with the ARM domains bending significantly inward and swinging out together with the TIR domains to form the "petals of the lotus". Nb-C6 bound to the SAM domain of the activated SARM1 and stabilized its Armadillo repeat motif domain. Analyses using hydrogen-deuterium exchange mass spectrometry (HDX-MS), and cross-linking MS (XL-MS) indicate that the activated SARM1 is highly dynamic and flexible and the neighboring TIRs form dimers via the surface close to one BB loop. The Nanobody is thus a valuable tool for delineating the mechanism of activation of SARM1 in AxD and other cellular processes.


title: Telomerase-independent maintenance of telomere length in a vertebrate

abstract:
Telomere shortening places a key limitation on cell proliferation1. In all vertebrates explored to date, this limitation is overcome by telomerase-dependent telomere extension. Failure to maintain telomere length results in premature ageing and functional impairments in highly replicative cell populations as telomeres erode2. Alternative lengthening of telomeres (ALT), a telomerase-independent mechanism, compensates for telomere loss in a subset of human cancer cell lines 2. Here, we demonstrate that the highly regenerative newt Pleurodeles waltl lacks telomerase activity, contains telomeres distinct from all known vertebrates in both sequence and structure, and deploys ALT for physiological telomere maintenance. This constitutes the first report of telomerase-independent resolution of the end-replication problem at the whole-organism level within Chordata.

One-Sentence SummaryP. waltl telomere biology is distinct amongst vertebrates and uses ALT at the whole-organism level.


title: Chemical inhibition of prolyl hydroxylases impairs angiogenic competence of human vascular endothelium through metabolic reprogramming

abstract: Endothelial cell (EC) metabolism has emerged as a driver of angiogenesis. While hypoxia inactivates prolyl-4 hydroxylase domain containing proteins 1-3 (PHD1-3) and stabilizes hypoxia inducible factors (HIFs) stimulating angiogenesis, the effects of PHDs on EC functions remain unclear. Here, we investigated the impact of PHD inhibition by dimethyloxalylglycine (DMOG) on angiogenic competence and metabolism of human vascular ECs. PHD inhibition reduced EC proliferation, migration, and tube formation capacities. Furthermore, transcriptomic and metabolomic analyses revealed an unfavorable metabolic reprogramming for angiogenesis following treatment with DMOG. Despite the induction of glycolytic genes and high levels of lactate, multiple genes encoding sub-units of mitochondrial complex I were suppressed with concurrent decline in nicotinamide adenine dinucleotide (NAD+) levels. Importantly, defective EC migration due to DMOG could be partially restored by augmenting NAD+ levels. Combined, our data provide metabolic insights into the mechanism by which chemical PHD inhibition impairs angiogenic competence of human vascular ECs.


title: Structural basis of transcription antitermination by Q{lambda}: NusA induces refolding of Q{lambda} to form nozzle for RNA polymerase exit channel

abstract:
Lambdoid bacteriophage Q proteins are transcription antipausing and antitermination factors that enable RNA polymerase (RNAP) to read through pause and termination sites. Q proteins load onto RNAP engaged in promoter-proximal pausing at a Q binding element (QBE) and adjacent sigma-dependent pause element to yield a Q-loading complex, and translocate with RNAP as a pausing-deficient, termination-deficient Q-loaded complex. In previous work, we showed that the Q protein from bacteriophage 21 (Q21) functions by forming a "nozzle" that narrows and extends the RNAP RNA-exit channel, preventing formation of pause and termination RNA hairpins. Here, we report atomic structures of four states on the pathway of antitermination by the Q protein from bacteriophage {lambda} (Q{lambda}), a Q protein that shows no sequence similarity to Q21 and that, unlike Q21, requires the transcription elongation factor NusA for efficient antipausing and antitermination. We report structures of Q{lambda}, the Q{lambda}-QBE complex, the NusA-free "pre-engaged" Q{lambda}-loading complex, and the NusA-containing "engaged" Q{lambda}-loading complex. The results show that Q{lambda}, like Q21, forms a nozzle that narrows and extends the RNAP RNA-exit channel, preventing formation of RNA hairpins. However, the results show that Q{lambda} has no three-dimensional structural similarity to Q21, employs a different mechanism of QBE recognition than Q21, and employs a more complex process for loading onto RNAP than Q21, involving recruitment of Q{lambda} to form a "pre-engaged" loading complex, followed by NusA-facilitated refolding of Q{lambda} to form an "engaged" loading complex. The results establish Q{lambda} and Q21 are not structural homologs and are solely functional analogs.

SIGNIFICANCE STATEMENTBacteriophage Q proteins are textbook examples of regulators of gene expression that function at the level of transcription antitermination. Here, we report structures defining the mechanism of antitermination by the Q protein of bacteriophage {lambda} (Q{lambda}). The results show Q{lambda} forms a "nozzle" that narrows and extends the RNA polymerase RNA-exit channel, precluding the formation of terminator RNA hairpins. The results show Q{lambda} exhibits no structural similarity to the Q protein of bacteriophage 21 (Q21), employs a different mechanism for DNA binding than Q21, and employs a more complex process of loading onto RNA polymerase than Q21. We conclude Q{lambda} and Q21 are not structural homologs and are solely functional analogs, akin to a bird wing and a bat wing.


title: Reciprocity and interaction effectiveness in generalised mutualisms among free-living species

abstract: Mutualistic interactions among free-living species generally involve weak links and highly asymmetric dependence among partners, yet our understanding of factors beyond their emergence is still limited. Using individual-based interactions of a super-generalist fleshy-fruited plant with its frugivore assemblage we estimate the Resource Provisioning Effectiveness (RPE) and Seed Dispersal Effectiveness (SDE) to assess the balance in the exchange of resources. Plants were highly dependent on a few super-generalist frugivore species, while these interacted with most individual plants, resulting in strong asymmetries in mutual dependence. Both RPE and SDE were mainly driven by interaction frequency. Despite highly asymmetric dependences, the strong reliance on quantity largely determined high reciprocity in rewards between partners (i.e., higher energy provided, more seedlings recruited), not obscured by minor variations in the quality of animal or plant service. We anticipate reciprocity will emerge in low-intimacy mutualisms where the mutualistic outcome largely relies upon interaction frequency.


title: Unearthing the global impact of mining construction minerals on biodiversity

abstract: Construction minerals - sand, gravel, limestone - are the most extracted solid raw materials1 and account for most of the worlds anthropogenic mass, which as of 2020 outweighed all of Earths living biomass2. However, knowledge about the magnitude, geography, and profile of this widespread threat to biodiversity remains scarce and scattered3-6. Combining long-term data from the IUCN Red List and new species descriptions we provide the first systematic evaluation of species threatened by mining of construction minerals globally. We found 1,047 species in the Red List impacted by this type of mining, of which 58.5% are threatened with extinction and four species already went extinct. We also identified 234 new species descriptions in 20 biodiversity hotspots reporting impacts from mining. Temporal trends in the assessments highlight the increased saliency of this threat to biodiversity, whose full extent may well reach over 24,000 animal and plant species. While rock quarrying mostly threatens karst biodiversity and narrow-ranged species, sand and gravel extraction is a more prominent threat to freshwater and coastal systems. This study provides the first evidence base to support a global strategy to limit the biodiversity impacts of construction mineral extraction.


title: On-farm Implementation of Midseason Drainage to Decrease Greenhouse Gas Emissions and Grain Arsenic Concentration in Rice Systems

abstract: Rice (Oryza sativa L.) is an important staple crop throughout much of the world, however, it is also a significant source of agricultural methane (CH4) emissions and exposure to arsenic (As). Introduction of soil aerobic events through practices such as alternate wetting and drying or midseason drainage, in flooded rice systems can significantly decrease grain As concentration and seasonal CH4 emissions. Previous on small research plots research has shown that a single midseason drain accomplishes these goals without yield reduction, but the degree of benefit depends on soil-drying severity. A midseason drain also has the potential to fit in well within current management practices of California rice systems, however, it has not been tested across a wide range of soil types or at a scale that farmers typically manage in this region. Therefore, in this three-year study, we aimed to determine if the results from previous small plot research are similar to what can be expected on-farm. At seven on-farm trials we implemented a single midseason drain and compared the grain yields, GHG emissions, and As concentration to the traditional farmer practice (FP) practiced in an adjacent part of the field. Soil moisture parameters [perched water table, volumetric water content, gravimetric water content (GWC), and soil water potential], CH4 and nitrous oxide (N2O) emissions, grain As and cadmium concentration, and grain yield were quantified. Across site-years, midseason drainage reduced seasonal CH4 emissions by 20-77%, compared to the FP control with the magnitude of reduction related to the soil-drying severity. For every 1% reduction in soil GWC during the drainage period, CH4 emissions were reduced by approximately 3.2%, compared to 2.5% in previous on-station research using small plots. Midseason drainage increased N2O emissions (average = 0.248 kg N2O-N ha-1) compared to the control but this accounted for only 3% of the seasonal global warming potential across all drainage treatments. Drainage also decreased grain As concentration by approximately 20%, on average, but was not related to the degree of soil-drying. Importantly, midseason drainage had no significant impact on grain yields. Overall, these results confirm findings from previous on-station research, indicating that midseason drainage may be a viable on-farm management practice for GHG mitigation and for reducing grain As concentration in flooded rice fields with limited risk of yield reduction.


title: HIV-1 CD4-binding site germline antibody-Env structures inform vaccine design

abstract: BG24, a VRC01-class broadly neutralizing antibody (bNAb) against HIV-1 Env with relatively few somatic hypermutations (SHMs), represents a promising target for vaccine strategies to elicit CD4-binding site (CD4bs) bNAbs. To understand how SHMs correlate with BG24 neutralization of HIV-1, we solved 4.1 [A] and 3.4 [A] single-particle cryo-EM structures of two inferred germline (iGL) BG24 precursors complexed with engineered Env-based immunogens lacking CD4bs N-glycans. Structures revealed critical Env contacts by BG24iGL and identified antibody light chain structural features that impede Env recognition. In addition, biochemical data and cryo-EM structures of BG24iGL variants bound to Envs with CD4bs glycans present provided insights into N-glycan accommodation, including structural modes of light chain adaptations in the presence of the N276gp120 glycan. Together, these findings revealed Env regions critical for germline antibody recognition and potential sites to alter in immunogen design.

category: biochemistry


title: Structural visualization of the tubulin folding pathway directed by eukaryotic chaperonin TRiC

abstract: SUMMARYThe ATP-dependent ring-shaped chaperonin TRiC/CCT is essential for cellular proteostasis. To uncover why some eukaryotic proteins can only fold with TRiC assistance, we reconstituted the folding of {beta}-tubulin using human Prefoldin and TRiC. We find unstructured {beta}-tubulin is delivered by Prefoldin to the open TRiC chamber followed by ATP-dependent chamber closure. CryoEM resolves four near-atomic resolution structures containing progressively folded {beta}-tubulin intermediates within the closed TRiC chamber, culminating in native tubulin. This substrate folding pathway appears closely guided by site-specific interactions with conserved regions in the TRiC chamber. Initial electrostatic interactions between the TRiC interior wall and both the folded tubulin N-domain and its C-terminal E-hook tail establish the native substrate topology, thus enabling C-domain folding. Disordered CCT C-termini within the chamber promote subsequent folding of tubulin Core and Middle domains and GTP-binding. Thus, TRiCs chamber provides chemical and topological directives that shape the folding landscape of its obligate substrates.


title: Sub-cellular dynamic investigation of the multi-component drug on the gastric cancer cell BGC823 using Raman spectroscopy

abstract: The potential of Raman spectroscopy in anticancer drug study has been demonstrated, yet its ability to character systematic cellular changes caused by multi-component drugs has not been explored. Here we used micro-Raman spectroscopy combined with bright field imaging to study Compound Kushen injection (CKI) at a sub-cellular level including intracellular vesicles(IVs). In our report, CKI caused dysfunction of DNA replication and repair was displayed by Raman spectrum (RS) from the cell nucleus. Meanwhile, the dynamics of CKI induced intracellular vesicles and cell component deconstruction was delineated by RS from the cytoplasm and IVs. The lipids-related biomolecular changes were also presented by the cytoplasm RS: the lipids level in the cytoplasm first descended then uprising. In conclusion, this study validated the mechanism and displayed the dynamics of CKI in treating cancer cells. We proved the capability of subcellular micro-Raman spectroscopy for detecting systematic cellular changes and its application for multi-component drug evaluation.


title: Glioblastoma states are defined by cohabitating cellular populations with progression-, imaging- and sex-distinct patterns

abstract: Glioblastomas (GBMs) are biologically heterogeneous within and between patients. Many previous attempts to characterize this heterogeneity have classified tumors according to their omics similarities. These discrete classifications have predominantly focused on characterizing malignant cells, neglecting the immune and other cell populations that are known to be present. We leverage a manifold learning algorithm to define a low-dimensional transcriptional continuum along which heterogeneous GBM samples organize. This reveals three polarized states: invasive, immune/inflammatory, and proliferative. The location of each sample along this continuum correlates with the abundance of eighteen malignant, immune, and other cell populations. We connect these cell abundances with magnetic resonance imaging and find that the relationship between contrast enhancement and tumor composition varies with patient sex and treatment status. These findings suggest that GBM transcriptional biology is a predictably constrained continuum that contains a limited spectrum of viable cell cohabitation ecologies. Since the relationships between this ecological continuum and imaging vary with patient sex and tumor treatment status, studies that integrate imaging features with tumor biology should incorporate these variables in their design.


title: The effect of Septoria glycines and fungicide application on the soybean phyllosphere mycobiome

abstract: Septoria brown spot (caused by Septoria glycines) is the most prevalent soybean disease in Illinois. It is common to use a foliar fungicide to control Septoria brown spot and other foliar diseases. The effects of fungicide on non-target organisms in the phyllosphere are unknown. To study the effect of S. glycines and fungicide application on the soybean phyllosphere mycobiome we sequenced full-length ITS and partial LSU region using oxford nanopore technologies. Sequencing produced 3,342 operational taxonomic units (OTUs). The richness of the fungal community significantly increased with the developmental stage. The soybean lines significantly affected the mycobiome diversity at the early developmental stage but not at the reproductive stages. S. glycines did not significantly affect the alpha diversity but some significant changes were observed for the beta diversity. At the R5 stage, fungicide application significantly changed the composition of the fungal community. The fungicide treatment significantly decreased the proportion of several fungal reads, but it increased the proportion of Septoria. The core microbiome in soybean leaves was composed of Gibberella, Alternaria, Didymella, Cladosporium, Plectosphaerella, Colletotrichum, and Bipolaris. Network analysis identified significant interactions between Septoria reads and Diaporthe, Bipolaris and two other taxonomic units. In this study, we set Septoria as the target organism and demonstrated that metabarcoding could be a tool to quantify the effect of multiple treatments on the fungal phyllosphere community. Better understanding the dynamics of the phyllosphere microbiome is necessary to start untangling late-season diseases of soybean.


title: Differentiation of airway cholinergic neurons from human pluripotent stem cells for airway neurobiology studies

abstract: Airway cholinergic nerves play a key role in airway physiology and disease. In asthma and other diseases of the respiratory tract, airway cholinergic neurons undergo plasticity and contribute to airway hyperresponsiveness and mucus secretion. We currently lack mechanistic understanding of airway cholinergic neuroplasticity due to the absence of human in vitro models. Here, we developed the first human in vitro model for airway cholinergic neurons using human pluripotent stem cell (hPSC) technology. hPSCs were differentiated towards mature and functional airway cholinergic neurons via a vagal precursor. Airway cholinergic neurons were characterized by ChAT and VAChT expression, and responded to chemical stimulation with changes in Ca2+ mobilization. Co-culture of hPSC-derived airway cholinergic neurons with airway smooth muscle cells enhanced phenotypic and functional characteristics of these neurons. The differentiation protocol we developed for human airway cholinergic neurons from hPSCs allows for studies into airway neurobiology and airway neuroplasticity in disease.


title: EPHierStats: a statistical tool to model the hierarchical relationships in electrophysiological data

abstract: Electrophysiological datasets are typically analyzed under the assumption that repeated measurements of the same unit of analysis (i.e. neuron or animal) can be treated as statistically independent. Recently, this assumption has been questioned and our data confirms and quantifies this skepticism using ex vivo slice recordings of synaptic currents in D1R+ medium spiny neurons in the nucleus accumbens. We therefore present EPHierStats as a statistical framework to analyze electrophysiological datasets with large numbers of measurements (>100) per unit of analysis. This novel analysis framework enables encoding of the full hierarchical relationships between measurements in a mixed-effects general linear model while also analyzing the distribution of values in assessed variables. Our method can easily be adapted to analyze a wide range of repeated-measures electrophysiological experiments. Implementation of the EPHierStats tool will aid the adaption of modern statistical approaches that prevent pseudoreplication and its associated false discovery rate while enabling statistical assessments of the complex relationships inherent to the field of neuroscience.


title: GASEOUS TRANSMITTER MODULATION OF HYPOGLOSSAL MOTONEURON ACTIVITY

abstract: Obstructive sleep apnea (OSA) is characterized by periodic disruption of breathing due to sporadic collapse of upper airway. Mice with targeted deletion of the gene Hmox2, encoding the enzyme heme oxygenase-2 (HO-2), exhibit severe OSA. Upper airway patency is governed by genioglossal nerve activity originating from the hypoglossal motor nucleus. We tested the hypothesis that disrupted HO-2 signaling weakens hypoglossal motoneuron activity using a combination of genetic, biochemical and electrophysiological approaches. Hypoglossal neurons are positive for HO-2. Disrupting HO-2 function caused transmission failure of inspiratory rhythm from preBotzinger complex (preBotC) to the hypoglossal nucleus. Failed transmission coincided with smaller inspiratory drive currents and a reduced intrinsic excitability among hypoglossal neurons. Hypoglossal neurons also expressed cystathionine {square}-lyase (CSE), an H2S producing enzyme. H2S abundance was higher in hypoglossal neurons of HO-2 null mice than wild-type controls. Disrupting CSE function normalized transmission from preBotC to hypoglossal neurons in HO-2 null mice. An H2S donor mimicked the effects of disrupted HO-2 function in wild type mice. These findings demonstrate a hitherto uncharacterized modulation of hypoglossal motoneurons by HO-2 and CSE-derived H2S.


title: Coregistration of heading to visual cues in retrosplenial cortex

abstract: Spatial cognition depends on an accurate representation of orientation within an environment. Head direction cells in distributed brain regions receive a range of sensory inputs, but visual input is particularly important for aligning their responses to environmental landmarks. To investigate how population-level heading responses are aligned to visual input, we recorded from retrosplenial cortex (RSC) of head-fixed mice in a moving environment using two-photon calcium imaging. We show that RSC neurons are tuned to the animals relative orientation in the environment, even in the absence of head movement. Next, we found that RSC receives functionally distinct projections from visual and thalamic areas, and contains several functional classes of neurons. While some functional classes mirror RSC inputs, a newly discovered class coregisters visual and thalamic signals. Finally, decoding analyses reveal unique contributions to heading from each class. Our results suggest an RSC circuit for anchoring heading representations to environmental visual landmarks.