title: A Novel Propidium Monoazide-Based PCR Assay Can Measure Viable Uropathogenic E. coli in Vitro and in Vivo

abstract:
BackgroundPolymerase chain reaction (PCR) is an important means by which to study the urine microbiome and is emerging as possible alternative to urine cultures to identify pathogens that cause urinary tract infection (UTI). However, PCR is limited by its inability to differentiate DNA originating from viable, metabolically active versus non-viable, inactive bacteria. This drawback has led to concerns that urobiome studies and PCR-based diagnosis of UTI are confounded by the presence of relic DNA from non-viable bacteria in urine. Propidium monoazide (PMA) dye can penetrate cells with compromised cell membranes and covalently bind to DNA, rendering it inaccessible to amplification by PCR. Although PMA has been shown to differentiate between non-viable and viable bacteria in various settings, its effectiveness in urine has not been previously studied. We sought to investigate the ability of PMA to differentiate between viable and non-viable bacteria in urine.

MethodsVarying amounts of viable or non-viable uropathogenic E. coli (UTI89) or buffer control were titrated with mouse urine. The samples were centrifuged to collect urine sediment or not centrifuged. Urine samples were incubated with PMA and DNA cross-linked using blue LED light. DNA was isolated and uidA gene-specific PCR was performed. For in vivo studies, mice were inoculated with UTI89, followed by ciprofloxacin treatment or no treatment. After the completion of ciprofloxacin treatment, an aliquot of urine was plated on non-selective LB agar and another aliquot was treated with PMA and subjected to uidA-specific PCR.

ResultsPMAs efficiency in excluding DNA signal from non-viable bacteria was significantly higher in bacterial samples in phosphate-buffered saline (PBS, dCT=13.69) versus bacterial samples in unspun urine (dCT=1.58). This discrepancy was diminished by spinning down urine-based bacterial samples to collect sediment and resuspending it in PBS prior to PMA treatment. In 3 of 5 replicate groups of UTI89-infected mice, no bacteria grew in culture; however, there was PCR amplification of E. coli after PMA treatment in 2 of those 3 groups.

ConclusionWe have successfully developed PMA-based PCR methods for amplifying DNA from live bacteria in urine. Our results suggest that non-PMA bound DNA from live bacteria can be present in urine, even after antibiotic treatment.

This indicates that viable but non-culturable E. coli can be present following treatment of UTI, and may explain why some patients have persistent symptoms but negative urine cultures following UTI treatment.


title: Mitochondrial active Ras2 protein promotes apoptosis and regulated cell death in a cAMP/PKA pathway-dependent manner in budding yeast.

abstract: In previous papers, using the eGFP-RBD3 probe, which binds Ras-GTP with high affinity, we showed that activated Ras proteins are localized to the plasma membrane and in the nucleus in wild-type Saccharomyces cerevisiae cells growing exponentially on glucose, while an aberrant accumulation of activated Ras in mitochondria correlates to mitochondrial dysfunction, accumulation of ROS and an increase of apoptosis. In this paper, we show that lack of TPS1, which is known to trigger apoptosis in S. cerevisiae, induces localization of active Ras proteins in mitochondria, confirming the above-mentioned correlation. Next, by characterizing the ras1{Delta} and ras2{Delta} mutants concerning localization of active Ras proteins and propensity to undergo cell death, we show that active Ras2 proteins, which accumulate in the mitochondria following addition of acetic acid, a well-known pro-apoptotic stimulus, might be the GTPases involved in regulated cell death, while active Ras1 proteins, constitutively localized in mitochondria, might be involved in a pro-survival molecular machinery. Finally, by characterizing the gpa2{Delta} and cyr1{Delta} mutants concerning the propensity to undergo cell death, we show that active mitochondrial Ras proteins promote apoptosis through the cAMP/PKA pathway.

category: microbiology


title: Guanosine-specific single-stranded ribonuclease effectors of a phytopathogenic fungus potentiate host immune responses

abstract: Plants activate immunity upon recognition of pathogen-associated molecular patterns. Although phytopathogens have evolved a set of effector proteins to counteract plant immunity, some effectors are perceived by hosts and induce immune responses. Here, we show that two secreted ribonuclease effectors, SRN1 and SRN2, encoded in a phytopathogenic fungus, Colletotrichum orbiculare, induce cell death in a signal peptide- and catalytic residue-dependent manner, when transiently expressed in Nicotiana benthamiana. The pervasive presence of SRN genes across Colletotrichum species suggested the conserved roles. Using a transient gene expression system in cucumber (Cucumis sativus), an original host of C. orbiculare, we show that SRN1 and SRN2 potentiate host pattern-triggered immunity. Consistent with this, C. orbiculare SRN1 and SRN2 deletion mutants exhibited increased virulence on the host. In vitro analysis revealed that SRN1 specifically cleaves single-stranded RNAs at guanosine, leaving a 3'-end phosphate. This activity has not been reported in plants. Importantly, the potentiation of C. sativus responses by SRN1 and SRN2 depends on the signal peptide and ribonuclease catalytic residues, suggesting that secreted SRNs cleave RNAs in apoplast and are detected by the host. We propose that the pathogen-derived apoplastic guanosine-specific single-stranded endoribonucleases lead to immunity potentiation in plants.


title: Taxonomic classification methods reveal a new subgenus in the paramyxovirus subfamily Orthoparamyxovirinae

abstract: As part of a broad One Health surveillance effort to detect novel viruses in wildlife and people, we report several paramyxoviruses sequenced primarily from bats during 2013 and 2014 in Brazil and Malaysia, including seven from which we recovered full-length genomes. Of these, six represent the first full-length paramyxovirus genomes sequenced from the Americas, including two sequences which are the first full-length bat morbillivirus genomes published to date. Our findings add to the vast number of viral sequences in public repositories that have been increasing considerably in recent years due to the rising accessibility of metagenomics. Taxonomic classification of these sequences in the absence of phenotypic data has been a significant challenge, particularly in the paramyxovirus subfamily Orthoparamyxovirinae, where the rate of discovery of novel sequences has been substantial. Using pairwise amino acid sequence classification (PASC), we describe a novel genus within this subfamily tentatively named Jeishaanvirus, which we propose should include as subgenera Jeilongvirus, Shaanvirus, and a novel South American subgenus Cadivirus. We also highlight inconsistencies in the classification of Tupaia virus and Mojiang virus using the same demarcation criteria and show that members of the proposed subgenus Shaanvirus are paraphyletic. Importantly, this study underscores the critical importance of sequence length in PASC analysis as well as the importance of biological characteristics such as genome organization in the taxonomic classification of viral sequences.


title: A pseudorabies virus serine/threonine kinase, US3, promotes retrograde transport in axons via Akt/mToRC1

abstract:
Infection of peripheral axons by alpha herpesviruses (AHVs) is a critical stage in establishing a life-long infection in the host. Upon entering the cytoplasm of axons, AHV nucleocapsids and associated inner-tegument proteins must engage the cellular retrograde transport machinery to promote the long-distance movement of virion components to the nucleus. The current model outlining this process is incomplete and further investigation is required to discover all viral and cellular determinants involved as well as the temporality of the events. Using a modified tri-chamber system, we have discovered a novel role of the pseudorabies virus (PRV) serine/threonine kinase, US3, in promoting efficient retrograde transport of nucleocapsids. We discovered that transporting nucleocapsids move at similar velocities both in the presence and absence of a functional US3 kinase; however fewer nucleocapsids are moving when US3 is absent and move for shorter periods of time before stopping, suggesting US3 is required for efficient nucleocapsid engagement with the retrograde transport machinery. This led to fewer nucleocapsids reaching the cell bodies to produce a productive infection 12hr later. Furthermore, US3 was responsible for the induction of local translation in axons as early as 1hpi through the stimulation of a PI3K/Akt-mToRC1. These data describe a novel role for US3 in the induction of local translation in axons during AHV infection, a critical step in transport of nucleocapsids to the cell body.

ImportanceNeurons are highly polarized cells with axons that can reach centimeters in length. Communication between axons at the periphery and the distant cell body is a relatively slow process involving the active transport of chemical messengers. Theres a need for axons to respond rapidly to extracellular stimuli. Translation of repressed mRNAs present within the axon occurs to enable rapid, localized responses independently of the cell body. AHVs have evolved a way to hijack local translation in the axons to promote their transport to the nucleus. We have determined the cellular mechanism and viral components involved in the induction of axonal translation. The US3 serine/threonine kinase of PRV activates Akt-mToRC1 signaling pathways early during infection to promote axonal translation. When US3 is not present, the number of moving nucleocapsids and their processivity are reduced, suggesting that US3 activity is required for efficient engagement of nucleocapsids with the retrograde transport machinery.

category: molecular biology


title: Differentiation of mouse embryonic stem cells into cells with spermatogonia-like morphology with chemical intervention-dependent increased gene expression of LIM homeobox 1 (Lhx1).

abstract: During the development of the male germline, spermatogonial stem cells (SSCs) originate from gonocytes that differentiate from primordial germ cells (PGCs). In the developing and regenerating mouse testis, expression of the gene LIM homeobox 1 (Lhx1) marks the most undifferentiated SSCs. However, an enrichment of Lhx1 expression in spermatogonia-like cells generated in vitro has not been reported so far. Previously, it was shown that a chemical intervention in male mouse embryonic stem (ES) cells in serum culture, including a timed combination of the SIRT1 inhibitor Ex-527, the DNA methyltransferase inhibitor RG-108 and the electrophilic redox cycling compound tert-butylhydroquinone (tBHQ), was associated with molecular markers of the PGC to gonocyte differentiation process. Here, we report the in vitro differentiation of male mouse ES cells, cultured under dual chemical inhibition of GSK3{beta} and MEK (termed 2i) with leukemia inhibitory factor (LIF) (termed 2iL) and serum, into cells with spermatogonia-like morphology (CSMs) and population-averaged expression of spermatogonia-specific genes. This was achieved by the removal of 2iL and a specific schedule of 2 partial medium replacements per day with alternating 8-hour and 16-hour intervals over a period of 32 days. Combination of this new cell culture protocol with the previously reported chemical intervention in ES cells changed the population-averaged expression of spermatogonia- and gonocyte-specific genes during the differentiation process and increased the population-averaged gene expression of Lhx1 in the resulting CSMs compared to CSMs without chemical intervention. Furthermore, we detected single CSMs with a strong nuclear LHX1/5 protein signal only in the chemical intervention group. Our results provide the first experimental evidence for the generation of CSMs with an enrichment of Lhx1 expression in vitro. We propose that further investigation of the CSMs generated with this in vitro system may provide new insights into male germline and stem cell development.

category: cell biology


title: Archaeal origins of gamete fusion

abstract:
Sexual reproduction consists of genome reduction by meiosis and subsequent gamete fusion. Presence of meiotic genes in prokaryotes suggests that DNA repair mechanisms evolved toward meiotic recombination; however, fusogenic proteins resembling those found in eukaryotes were not identified in prokaryotes. Here, we identify archaeal proteins that are homologs of fusexins, a superfamily of fusogens that mediate eukaryotic gamete and somatic cell fusion, as well as virus entry. The crystal structure of a trimeric archaeal Fusexin1 reveals novel features such as a six-helix bundle and an additional globular domain. Ectopically expressed Fusexin1 can fuse mammalian cells, and this process involves the additional domain and a conserved fusion loop. Archaeal fusexin genes exist within integrated mobile elements, potentially linking ancient archaeal gene exchanges and eukaryotic sex.

One-Sentence SummaryCell membrane fusion proteins of viruses and eukaryotes are also present in archaea.

category: evolutionary biology


title: Discovery of archaeal Fusexins homologous to eukaryotic HAP2/GCS1 gamete fusion proteins


title: The Nup2 meiotic-autonomous region relieves autoinhibition of Nup60 to promote progression of meiosis and sporulation in Saccharomyces cerevisiae

abstract: During meiosis, chromosomes undergo dramatic changes in structural organization, nuclear positioning, and motion. Although the nuclear pore complex has been shown to affect genome organization and function in vegetative cells, its role in meiotic chromosome dynamics has remained largely unexplored. Recent work in the budding yeast Saccharomyces cerevisiae demonstrated that the mobile nucleoporin Nup2 is required for normal progression through meiosis I prophase and sporulation in strains where telomere-led chromosome movement has been compromised. The meiotic autonomous region (MAR), a short fragment of Nup2 responsible for its role in meiosis, was shown to localize to the nuclear envelope via Nup60 and to bind to meiotic chromosomes. To understand the relative contribution these two activities have on MAR function, we first carried out a screen for MAR mutants defective in sporulation and found that all the mutations disrupt interaction with both Nup60 and meiotic chromosomes. Moreover, nup60 mutants phenocopy nup2 mutants, exhibiting similar nuclear division kinetics, sporulation efficiencies, and genetic interactions with mutations that affect the telomere bouquet. Although full-length Nup60 requires Nup2 for function, removal of Nup60s C-terminus allows Nup60 to bind meiotic chromosomes and promote sporulation without Nup2. In contrast, binding of the MAR to meiotic chromosomes is completely dependent on Nup60. Our findings uncover an inhibitory function for the Nup60 C-terminus and suggest that Nup60 mediates recruitment of meiotic chromosomes to the nuclear envelope, while Nup2 plays a secondary role counteracting Nup60s autoinhibition.

category: genetics


title: Single-cell transcriptomics reveals expression profiles of Trypanosoma brucei sexual stages

abstract: Early diverging lineages such as trypanosomes can provide clues to the evolution of sexual reproduction in eukaryotes. In Trypanosoma brucei, the pathogen that causes Human African Trypanosomiasis, sexual reproduction occurs in the salivary glands of the insect host, but analysis of the molecular signatures that define these sexual forms is complicated because they mingle with more numerous, mitotically-dividing developmental stages. We used single-cell RNA-sequencing (scRNAseq) to profile 388 individual trypanosomes from midgut, proventriculus, and salivary glands of infected tsetse flies allowing us to identify tissue-specific cell types. Further investigation of salivary gland parasite transcriptomes revealed fine-scale changes in gene expression over a developmental progression from putative sexual forms through metacyclics expressing variant surface glycoprotein genes. The cluster of cells potentially containing sexual forms was characterized by high level transcription of the gamete fusion protein HAP2, together with an array of surface proteins and several genes of unknown function. We linked these expression patterns to distinct morphological forms using immunofluorescence assays and reporter gene expression to demonstrate that the kinetoplastid-conserved gene Tb927.10.12080 is exclusively expressed at high levels by meiotic intermediates and gametes. We speculate that this protein, currently of unknown function, plays a role in gamete formation and/or fusion.

category: genomics


title: Characterisation of a highly diverged mitochondrial ATP synthase peripheral stalk subunit b in Trypanosoma brucei

abstract: The mitochondrial F1Fo ATP synthase of Trypanosoma brucei has been studied in detail. Whereas its F1 moiety is relatively highly conserved in structure and composition, the same is not the case for the Fo part and the peripheral stalk. A core subunit of the latter, the normally conserved subunit b, could not be identified in trypanosomes suggesting that it might be absent. Here we have identified a 17 kDa mitochondrial protein of the inner membrane that is essential for normal growth, efficient oxidative phosphorylation and membrane potential maintenance. Pulldown experiments and native PAGE analysis indicate that the protein is associated with the F1Fo ATP synthase. Its ablation reduces the levels of Fo subunits, but not those of F1, and disturbs the cell cycle. HHpred analysis showed that the protein has structural similarities to subunit b of other species, indicating that the Fo part of the trypanosomal ATP synthase does contain a highly diverged subunit b. Thus, the Fo part of the trypanosomal ATPase synthase may be more widely conserved than initially thought.

category: biochemistry


title: Truncation of poly-N-acetylglucosamine (PNAG) polymerization with N-acetylglucosamine analogues

abstract: Bacteria require polysaccharides for structure, survival, and virulence. Despite the central role these structures play in microbiology few tools are available to manipulate their production. In E. coli the glycosyltransferase complex PgaCD produces poly-N-acetylglucosamine (PNAG), an extracellular matrix polysaccharide required for biofilm formation. We report that C6-substituted (H, F, N3, SH, NH2) UDP-GlcNAc substrate analogues are inhibitors of PgaCD. In vitro the inhibitors cause PNAG chain termination; consistent with the mechanism of PNAG polymerization from the non-reducing terminus. In vivo, expression of the GlcNAc-1-kinase NahK in E. coli provided a non-native GlcNAc salvage pathway that produced the UDP-GlcNAc analogue inhibitors in situ. The 6-fluoro and 6-deoxy derivatives were potent inhibitors of biofilm formation in the transformed strain, providing a tool to manipulate this key exopolysaccharide. Characterization of the UDP-GlcNAc pool and quantification of PNAG generation support PNAG termination as the primary in vivo mechanism of biofilm inhibition by 6-fluoro UDP-GlcNAc.


title: Ribosomal protein S1 plays a critical role in horizontal gene transfer by mediating the expression of foreign mRNAs

abstract: The emergence of multi-antibiotic resistant bacteria is one of the largest threats to global heath. This rise is due to the genomic plasticity of bacteria, allowing rapid acquisition of antibiotic resistance through the uptake of foreign DNA (i.e. horizontal gene transfer, HGT). This genomic plasticity is not limited to DNA from bacteria, highly divergent (trans-kingdom) mRNA have been reported to drive translation in E. coli. Trans-kingdom activity has been attributed to mRNA tertiary structure suggesting the bacterial translation machinery bottle-necks HGT, restricting the expression of foreign DNA. However, here we show that tertiary structure is not responsible for ribosome recruitment and that the translation efficiency is dependent on ribosomal protein S1 and an A-rich Shine-Dalgarno-like element. The S1-facilitated ability of ribosomes to identify and exploit A-rich sequences in foreign RNA highlights the important role that S1 plays in horizontal gene transfer, the robustness of canonical prokaryotic translation, and bacterial evolution.


title: Replication Fork Uncoupling Causes Nascent Strand Degradation and Fork Reversal

abstract: Genotoxins cause nascent strand degradation (NSD) and fork reversal during DNA replication. NSD and fork reversal are crucial for genome stability and exploited by chemotherapeutic approaches. However, it is unclear how NSD and fork reversal are triggered. Additionally, the fate of the replicative helicase during these processes is unknown. We developed a biochemical approach to study synchronous, localized NSD and fork reversal using Xenopus egg extracts. We show that replication fork uncoupling stimulates NSD of both nascent strands and progressive conversion of uncoupled forks to reversed forks. The replicative helicase remains bound during NSD and fork reversal, indicating that both processes take place behind the helicase. Unexpectedly, NSD occurs before and after fork reversal, indicating that multiple degradation steps take place. Overall, our data show that uncoupling causes NSD and fork reversal and identify key steps involved in these processes.


title: Flow Cytometric Evaluation of Yeast-Bacterial Cell-Cell Interactions

abstract: Synthetic cell-cell interaction systems can be useful for understanding multicellular communities or for screening binding molecules. We adapt a previously characterized set of synthetic cognate nanobody-antigen pairs to a yeast-bacteria coincubation format and use flow cytometry to evaluate cell-cell interactions mediated by binding between surface-displayed molecules. We further use fluorescence-activated cell sorting (FACS) to enrich for a specific yeast-displayed nanobody within a mixed yeast-display population. Finally, we demonstrate that this system supports characterization of a therapeutically relevant nanobody-antigen interaction: a previously discovered nanobody that binds to the intimin protein expressed on the surface of enterohemorrhagic E. coli. Overall, our findings indicate that the yeast-bacteria format supports efficient evaluation of ligand-target interactions. With further development, this format may facilitate systematic characterization and high throughput discovery of bacterial surface-binding molecules.

category: synthetic biology


title: Obesity influences composition of salivary and fecal microbiota and impacts the interactions between bacterial taxa

abstract:
Obesity is an increasing global health concern and is associated with a broad range of morbidities. The gut microbiota are increasingly recognized as important contributors to obesity and cardiometabolic health. This study aimed to characterize oral and gut microbial communities, and evaluate host:microbiota interactions between clinical obesity classifications. We performed 16S rDNA sequencing on fecal and salivary samples, global metabolomics profiling on plasma and stool samples, and dietary profiling in 135 healthy individuals. We grouped individuals by obesity status, based on body mass index (BMI), including lean (BMI 18-24.9), overweight (BMI 25-29.9), or obese (BMI [≥]30). We analyzed differences in microbiome composition, community inter-relationships, and predicted microbial function by obesity status. We found that salivary bacterial communities of lean and obese individuals were compositionally and phylogenetically distinct. An increase in obesity status was positively associated with strong correlations between bacterial taxa, particularly with bacterial groups implicated in metabolic disorders including Fretibacterium, and Tannerella. Consumption of sweeteners, especially xylitol, significantly influenced compositional and phylogenetic diversities of salivary and fecal bacterial communities. In addition, obesity groups exhibited differences in predicted bacterial metabolic activity, which was correlated with hosts metabolite concentrations. Overall, obesity was associated with distinct changes in bacterial community dynamics, particularly in saliva. Consideration of microbiome community structure, and inclusion of salivary samples may improve our ability to understand pathways linking microbiota to obesity and cardiometabolic disease.

IMPORTANCEObesity is a worldwide epidemic that is associated with a wide range of health issues. Microbiota were shown to influence metabolism and obesity development. Our study aimed to evaluate the interactions between obesity, salivary and fecal microbiota, and metabolite concentrations in healthy individuals. The oral bacterial community was more impacted by the obesity status of the host than fecal microbiota. Consistently for oral and fecal microbiota, the number of strong interactions between bacteria increased with the increase in the obesity status. Several predicted microbial metabolic pathways that were shown to be associated with metabolic health were uniquely enriched between obesity groups. In addition, these metabolic pathways were correlated with plasma and stool metabolites. Our results suggest that oral microbiota might better reflect the obesity status of the host than fecal microbiota, and that correlations between microbial taxa are altered during obesity.


title: Co-cultivation is a powerful approach to produce a robust functionally designed synthetic consortium as a live biotherapeutic product (LBP)

abstract: The successes of fecal microbiota transplants (FMT) have provided the necessary proof-of-concept for microbiome therapeutics. Because of the many risks and uncertainties associated with feces-based therapies, defined microbial consortia that modify the microbiome in a targeted manner have emerged as a promising safer alternative to FMT. The development of such live biotherapeutic products has important challenges, including the selection of appropriate strains and the production of the consortia at scale. Here, we report on an ecology and biotechnology-based approach to microbial consortium design that overcomes these issues. We designed a nine-strain consortium that emulates the central metabolic pathways of carbohydrate fermentation in the healthy human gut microbiota. We show that continuous co-culturing the bacteria produce a stable consortium whose activity is distinct from an equivalent mix of individually cultured strains. Further, we showed that our function-based consortium is as effective as FMT in counteracting dysbiosis in a dextran sodium sulfate mouse model of acute colitis. We propose that combining a bottom-up functional design with continuous co-cultivation is a powerful strategy to produce robust, functionally designed synthetic consortia for therapeutic use.


title: 3D Structure Determination of Protein Complexes using Matrix-Landing Mass Spectrometry

abstract: Native mass spectrometry (MS) is an emerging technology that can provide complementary data to electron microscopy (EM) for protein structure characterization. Beyond the ability to provide mass measurements of gas-phase biomolecular ions, MS instruments offer the ability to purify, select, and precisely control the spatial location of these ions. Here we present a modified Orbitrap MS system capable of depositing a native MS ion beam onto EM grids. We further describe use of a chemical landing matrix that both preserves and protects the structural integrity of the deposited particles. With this system we obtained the first 3D reconstructed structure of gas-phase, deposited biomolecular ions - the 800 KDa protein complex GroEL. These data provide direct evidence that non-covalent protein complexes can indeed retain their condensed-phase structures following ionization and vaporization. Finally, we describe how further developments of this technology could pave the way to an integrated MS-EM technology with promise to provide improved cryo-EM sample preparation over conventional plunge-freezing techniques.


title: Early life stress-induced vulnerability to postpartum mental disturbance: prolonged dysregulation of the HPA axis and behavior

abstract: Early life stress (ELS) increases the risk for postpartum depression (PPD). Patients with depression who experienced ELS tend to be treatment refractory. Nevertheless, it remains elusive how ELS underlies the pathophysiology of PPD at the mechanistic level. Here we observed that adolescent social isolation in mice resulted in an aberrantly sustained elevation of glucocorticoids via hypothalamic-pituitary-adrenal (HPA) axis dysregulation and long-lasting behavioral deficits during the postpartum period. The postpartum behavioral changes elicited by ELS were not ameliorated by the medicines currently used for PPD in behavioral assays that are frequently used in drug discovery for human depression. However, post-delivery treatment with a glucocorticoid receptor antagonist effectively ameliorated the deficits. We also demonstrated a significant impact of ELS on the HPA axis dysregulation and PPD in humans. In summary, we show the validity and utility of a novel model and a novel option of treatment for a subset of treatment-refractory PPD.


title: Local and global chemical shaping of bacterial communities by redox potential

abstract: Thermodynamics predicts a positive correlation between environmental redox potential and oxidation state of molecules; if found for microbial communities it would imply a new kind of deterministic eco-evolutionary process. This study examines evidence for local- and global-scale correlations between oxidation-reduction potential (ORP or Eh) in environmental samples and carbon oxidation state (ZC) of estimated bacterial and archaeal community proteomes. Seventy-nine public datasets for seven environment types (river & seawater, lake & pond, alkaline spring, hot spring, groundwater, sediment, and soil) were analyzed. Taxonomic abundances inferred from high-throughput 16S rRNA gene sequences were combined with NCBI Reference Sequence (RefSeq) proteomes to estimate the amino acid compositions and chemical formulas (CcHhNnOoSs) of community proteomes, which yield ZC. Alkaline hot springs have the lowest ZC for both bacterial and archaeal domains of any environment. Positive global correlations between redox potential and ZC are found for bacterial communities in lake & pond, groundwater, and soil environments, but not archaeal communities, suggesting a broad ecological signal of chemical shaping in Bacteria.

category: ecology


title: Epigenetic clock and methylation studies in gray short-tailed opossums

abstract: The opossum (Monodelphis domestica), with its sequenced genome, ease of laboratory care and experimental manipulation, and unique biology, is the most used laboratory marsupial. Using the mammalian methylation array, we generated DNA methylation data from n=100 opossum tissues including blood, liver, and tail. We contrast age-related changes in the opossum methylome to those of C57BL/6J mice. We present several epigenetic clocks for opossums that are distinguished by their compatibility with tissue type (pan-tissue and blood clock) and species (opossum and human). Two dual-species human-opossum pan-tissue clocks accurately measure chronological age and relative age, respectively. These human-opossum epigenetic clocks are expected to provide a significant boost to the attractiveness of opossum as a biological model.


title: scINSIGHT for interpreting single-cell gene expression from biologically heterogeneous data

abstract: The increasing number of scRNA-seq data emphasizes the need for integrative analysis to interpret similarities and differences between single-cell samples. Even though different batch effect removal methods have been developed, none of the existing methods is suitable for het-erogeneous single-cell samples coming from multiple biological conditions. To address this challenge, we propose a method named scINSIGHT to learn coordinated gene expression patterns that are common among or specific to different biological conditions, offering a unique chance to identify cellular identities and key biological processes across single-cell samples. We have evaluated scINSIGHT in comparison with state-of-the-art methods using simulated and real data, which consistently demonstrate its improved performance. In addition, our results show the applicability of scINSIGHT in diverse biomedical and clinical problems.

category: bioinformatics


title: The R346K Mutation in the Mu Variant of SARS-CoV-2 Alter the Interactions with Monoclonal Antibodies from Class 2: A Free Energy of Perturbation Study

abstract: The Mu variant of SARS-CoV-2 has been recently classified as a variant of interest (VOI) by the world health organization (WHO) but limited data are available at the moment. In particular, a special attention was given to the R346K mutation located in the receptor binding domain (RBD). In the current study we performed Free energy of perturbation (FEP) calculations to elucidate it possible impact on a set of neutralizing monoclonal antibodies (mAbs) which have been shown to be strong inhibitors of the most other known COVID-19 variants. Our results show that R346K affects the class 2 antibodies but its effect is not so significant (0.66 kcal/mol); i.e. reduces the binding with RBD about 3 times. An identical value was calculated also in the presence of both class 1 and class 2 antibodies (BD-812/836). Further, a similar reduction in the binding (0.4 kcal/mol) was obtained for BD-821/771 pair of mAbs. For comparison, the addition of K417N mutation, present in the newly registered Mu variant in July 2021 in UK, affected the class 1 mAbs by 1.29 kcal/mol reducing stronger the binding by about 10 times. Thus, the resistance effect of R346K mutation in the Mu variant is possible but not so significant and is due to the additional decrease of antibody neutralization based on the reduced binding of class 2 antibodies.


title: Quantitative parameters of bacterial RNA polymerase open-complex formation, stabilization and disruption on a consensus promoter

abstract: Transcription initiation is the first step in gene expression, and is therefore strongly regulated in all domains of life. The RNA polymerase (RNAP) first associates with the initiation factor{sigma} to form a holoenzyme, which binds, bends and opens the promoter in a succession of reversible states. These states are critical for transcription regulation, but remain poorly understood. Here, we addressed the mechanism of open complex formation by monitoring its assembly/disassembly kinetics on individual consensus lacUV5 promoters using high-throughput single-molecule magnetic tweezers. We probed the key protein-DNA interactions governing the open-complex formation and dissociation pathway by modulating the dynamics at different concentrations of monovalent salts and varying temperatures. Consistent with ensemble studies, we observed that RPO is a stable, slowly reversible state that is preceded by a kinetically significant open intermediate (RPI), from which the holoenzyme dissociates. A strong anion concentration and type dependence indicates that the RPO stabilization may involve sequence-independent interactions between the DNA and the holoenzyme, driven by a non-Coulombic effect consistent with the non-template DNA strand interacting with{sigma} and the RNAP {beta} subunit. The temperature dependence provides the energy scale of open-complex formation and further supports the existence of additional intermediates.

category: biophysics


title: Chasing the fitness optimum: temporal variation in the genetic and environmental expression of life-history traits for a perennial plant

abstract: The ability of plants to track shifting fitness optima is crucial within the context of global change, where increasing environmental extremes may have dramatic consequences to life history, fitness, and ultimately species persistence. However, to track changing conditions relies upon the complex relationship between genetic and environmental variance, where selection may favor plasticity, the evolution of genetic differences, or both depending on the spatial and temporal scale of environmental heterogeneity. Over three years, we compared the genetic and environmental components of phenological and life-history variation in a common environment for the spring perennial Geum triflorum. Populations were sourced from alvar habitats that exhibit extreme, but predictable annual flood-desiccation cycles and prairie habitats that exhibit similar, but less predictable variation in water availability. Narrow-sense heritabilities were generally higher for early life history (emergence probability) relative to later life history traits (total seed mass), indicating that traits associated with establishment within an environment are under stronger genetic control relative to later life-history fitness expressions, where plasticity may play a larger role. This pattern was particularly notable in seeds sourced from environmentally extreme, but predictable alvar habitats relative to less predictable prairie seed sources. Fitness landscapes based on seed source origin, largely characterized by varying water availability and flower production, described selection as the degree of maladaptation to the prairie common garden environment relative to seed source environment. Plants from alvar populations were consistently closer to the fitness optimum across all years. Annually, the breadth of the fitness optimum expanded primarily along a moisture gradient, with inclusion of more populations onto the expanding optimum. These results highlight the importance of temporally and spatially varying selection for the evolution of life history, indicating plasticity within perennial systems may over time become the primary mechanism to track fitness for later life history events.


title: Phylogenomic analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving the diversifications, with insights into the complex biogeographic history in the Northern Hemisphere

abstract: Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, hundreds of single-copy nuclear (SCN) genes (785) and complete plastomes (75) were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, it was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event best explains the cytonuclear discordance that occurred in the Eocene in western North America. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. We detected one deep hybridization event (Malus doumeri) involving the ancestor of pome-bearing species and Docynia delavayi, and one recent hybridization event (Malus coronaria) between M. sieversii and a combined clade of M. ioensis and M. angustifolia. Furthermore, our historical biogeographic analysis combining living and fossil species supported a widespread East Asian-western North American origin of Malus s.l., followed by a series of extinction events in the Eocene in northern East Aisa and western North America. This study provides a valuable evolutionary framework for the breeding and crop improvement of apples and their close relatives.


title: Diversity of reptile sex chromosome evolution revealed by cytogenetic and linked-read sequencing

abstract: Reptile sex determination is attracting much attention because the great diversity of sex-determination and dosage compensation mechanisms permits us to approach fundamental questions about sex chromosome turnover and evolution. However, reptile sex chromosome variation remains largely uncharacterized and no reptile master sex determination genes have yet been identified. Here we describe a powerful and cost-effective "chromosomics" approach, combining probes generated from the microdissected sex chromosomes with transcriptome sequencing to explore this diversity in non-model Australian reptiles with heteromorphic or cryptic sex chromosomes. We tested the pipeline on a turtle, a gecko, and a worm-lizard, and we also identified sequences located on sex chromosomes in a monitor lizard using linked-read sequencing. Genes identified on sex chromosomes were compared to the chicken genome to identify homologous regions among the four species. We identified candidate sex determining genes within these regions, including conserved vertebrate sex-determining genes pdgfa, pdgfra amh and wt1, and demonstrated their testis or ovary-specific expression. All four species showed gene-by-gene rather than chromosome-wide dosage compensation. Our results imply that reptile sex chromosomes originated by independent acquisition of sex-determining genes on different autosomes, as well as translocations between different ancestral macro- and micro-chromosomes. We discuss the evolutionary drivers of the slow differentiation, but rapid turnover, of reptile sex chromosomes.


title: Screening of CHO-K1 endogenous promoters for expressing recombinant proteins in mammalian cell cultures

abstract: For the production of recombinant protein therapeutics in mammalian cells, a high rate of gene expression is desired and hence strong viral-derived promoters are commonly used. However, they usually induce cellular stress and can be susceptible to epigenetic silencing. Endogenous promoters, which coordinates their activity with cellular and bioprocess dynamics while at the same time they maintain high expression levels, may help to avoid such drawbacks. In this work, endogenous promoters were identified based on high expression levels in RNA-seq data of CHO-K1 cells cultured in high density. The promoters of Actb, Ctsz, Hmox1, Hspa5, Vim and Rps18 genes were selected for generating new expression vectors for the production of recombinant proteins in mammalian cells. The in silico-derived promoter regions were experimentally verified and the majority showed transcriptional activity comparable or higher than CMV. Also, stable expression following a reduction of culture temperature was investigated. The characterized endogenous promoters (excluding Rps18) constitute a promising alternative to CMV promoter due to their high strength, long-term expression stability and integration into the regulatory network of the host cell. These promoters may also comprise an initial panel for designing cell engineering strategies and synthetic promoters, as well as for industrial cell line development.


title: Top-down predictions of visual features dynamically reverse their bottom-up processing in the occipito-ventral pathway to facilitate stimulus disambiguation and behavior

abstract: The prevalent conception of vision-for-categorization suggests an interplay of two dynamic flows of information within the occipito-ventral pathway. The bottom-up flow progressively reduces the high-dimensional input into a lower-dimensional representation that is compared with memory to produce categorization behavior. The top-down flow predicts category information (i.e. features) from memory that propagates down the same hierarchy to facilitate input processing and behavior. However, the neural mechanisms that support such dynamic feature propagation up and down the visual hierarchy and how they facilitate behavior remain unclear. Here, we studied them using a prediction experiment that cued participants (N = 11) to the spatial location (left vs. right) and spatial frequency (SF, Low, LSF, vs. High, HSF) contents of an upcoming Gabor patch. Using concurrent MEG recordings of each participants neural activity, we compared the top-down flow of representation of the predicted Gabor contents (i.e. left vs. right; LSF vs. HSF) to their bottom-up flow. We show (1) that top-down prediction improves speed of categorization in all participants, (2) the top-down flow of prediction reverses the bottom-up representation of the Gabor stimuli, going from deep right fusiform gyrus sources down to occipital cortex sources contra-lateral to the expected Gabor location and (3) that predicted Gabors are better represented when the stimulus is eventually shown, leading to faster categorizations. Our results therefore trace the dynamic top-down flow of a predicted visual content that chronologically and hierarchically reversed bottom-up processing, further facilitates visual representations in early visual cortex and subsequent categorization behavior.

category: neuroscience


title: Honey bees can store and retrieve independent memory traces after complex experiences that combine appetitive and aversive associations

abstract: Real-world experiences do often mix appetitive and aversive events. Understanding the ability of animals to extract, store and use this information is an important issue in neurobiology. We used honey bees as model to study learning and memory after a differential conditioning that combines appetitive and aversive training trials. First of all, we describe an aversive conditioning paradigm that constitutes a clear opposite of the well known appetitive olfactory conditioning of the proboscis extension response. A neutral odour is presented paired with the bitter substance quinine. Aversive memory is evidenced later as an odour-specific impairment in appetitive conditioning. Then we tested the effect of mixing appetitive and aversive conditioning trials distributed along the same training session. Differential conditioning protocols like this were used before to study the ability to discriminate odours, however they were not focused on whether appetitive and aversive memories are formed. We found that after a differential conditioning, honey bees establish independent appetitive and aversive memories that do not interfere with each other during acquisition or storage. Finally, we moved the question forward to retrieval and memory expression to evaluate what happens when appetitive and the aversive learned odours are mixed during test. Interestingly, opposite memories compete in a way that they do not cancel each other out. Honey bees showed the ability to switch from expressing appetitive to aversive memory depending on their satiation level.

category: animal behavior and cognition


title: PoSTcode: Probabilistic image-based spatial transcriptomics decoder

abstract: Realising the full potential of novel image-based spatial transcriptomic (IST) technologies requires robust and accurate algorithms for decoding the hundreds of thousand fluorescent signals each derived from single molecules of mRNA. In this paper, we introduce PoSTcode, a probabilistic method for transcript decoding from cyclic multi-channel images, whose effectiveness is demonstrated on multiple large-scale datasets generated using different versions of the in situ sequencing protocols. PoSTcode is based on a re-parametrised matrix-variate Gaussian mixture model designed to account for correlated noise across fluorescence channels and imaging cycles. PoSTcode is shown to recover up to 50% more confidently decoded molecules while simultaneously decreasing transcript mislabeling when compared to existing decoding techniques. In addition, we demonstrate its increased stability to various types of noise and tuning parameters, which makes this new approach reliable and easy to use in practice. Lastly, we show that PoSTcode produces fewer doublet signals compared to a pixel-based decoding algorithm.


title: RNA-seq, bioinformatic identification of potential MicroRNA-Like Small RNAs in the edible mushroom Agaricus bisporus and experimental approach for their validation

abstract: Although genomes from many edible mushrooms are sequenced, studies on fungal miRNAs are scarce. Most of the bioinformatic tools are designed for plants or animals but fungal miRNAs processing and expression share similarities and differences with both kingdoms. Moreover, since mushroom species such as Agaricus bisporus (white button mushroom) are frequently consumed as food, controversial discussions are still evaluating whether their miRNAs might or might not be assimilated, perhaps within extracellular vesicles (i.e exosomes). Therefore, the A. bisporus RNA-seq was studied in order to identify potential de novo miRNA-like small RNAs (milRNAs) that might allow their later detection in the diet. Results pointed to 1 already known and 37 de novo milRNAss. Three milRNAss were selected for RT-qPCR experiments. Precursors and mature milRNAs were found in the edible parts (caps and stipes) validating the predictions carried out in silico. When their potential gene targets were investigated, results pointed that mostly were involved in primary and secondary metabolic regulation. However, when human transcriptome is used as target, the results suggest that they might interfere with important biological processes related with cancer, infectious process and neurodegenerative diseases.


title: Investigating the Role of Chromatin Remodeler FOXA1 in Ferroptotic Cell Death

abstract:
Ferroptosis is a lipid peroxidation-dependent mechanism of regulated cell death known to suppress tumor proliferation and progression. Although several genetic and protein hallmarks have been identified in ferroptotic cell death, it remains challenging to fully characterize ferroptosis signaling pathways and to find suitable biomarkers. Moreover, changes taking place in the epigenome of ferroptotic cells remain poorly studied. In this context, we aimed to investigate the role of chromatin remodeler forkhead box protein A1 (FOXA1) in RSL3-treated multiple myeloma cells because, similar to ferroptosis, this transcription factor has been associated with changes in the lipid metabolism, DNA damage, and epithelial-to-mesenchymal transition (EMT). RNA sequencing and Western blot analysis revealed that FOXA1 expression is consistently upregulated upon ferroptosis induction in different in vitro and in vivo disease models. In silico motif analysis and transcription factor enrichment analysis further suggested that ferroptosis-mediated FOXA1 expression is orchestrated by specificity protein 1 (Sp1), a transcription factor known to be influenced by lipid peroxidation. Remarkably, FOXA1 upregulation in ferroptotic myeloma cells did not alter hormone signaling or EMT, two key downstream signaling pathways of FOXA1. CUT&RUN genome-wide transcriptional binding site profiling showed that GPX4-inhibition by RSL3 triggered loss of binding of FOXA1 to pericentromeric regions in multiple myeloma cells, suggesting that this transcription factor is possibly involved in genomic instability, DNA damage, or cellular senescence under ferroptotic conditions.

O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=152 SRC="FIGDIR/small/461056v1_ufig1.gif" ALT="Figure 1">
View larger version (35K):
org.highwire.dtl.DTLVardef@1fe44acorg.highwire.dtl.DTLVardef@4ea970org.highwire.dtl.DTLVardef@560fd0org.highwire.dtl.DTLVardef@283874_HPS_FORMAT_FIGEXP  M_FIG C_FIG

category: cancer biology


title: Impact of baseline culture conditions of mouse-derived cancer organoids when determining therapeutic response and tumor heterogeneity

abstract: Representative models are needed to screen new therapies for patients with cancer. Cancer organoids are a leap forward as a culture model that faithfully represents the disease. Mouse-derived cancer organoids (MDCOs) are becoming increasingly popular, however there has yet to be a standardized method to assess therapeutic response and identify subpopulation heterogeneity. There are multiple factors unique to organoid culture that could affect how therapeutic response and MDCO heterogeneity are assessed. Here we describe an analysis of nearly 3,500 individual MDCOs where individual organoid morphologic tracking was performed. Change in MDCO diameter was assessed in the presence of control media or targeted therapies. Individual organoid tracking was identified to be more sensitive to treatment response than well-level assessment. The impact of different generations of mice of the same genotype, different regions of the colon, and organoid specific characteristics including baseline size, passage number, plating density, and location within the matrix were examined. Only the starting size of the MDCO altered the subsequent growth. Here we establish organoid culture parameters for individual organoid morphologic tracking to determine therapeutic response and growth/response heterogeneity for translational studies using murine colorectal cancer organoids.


title: Creatine transporter deficient rat show motor dysfunction linked with cerebellar alterations and muscle creatine deficiency without muscle atrophy

abstract: Creatine (Cr) is a nitrogenous organic acid and plays roles as fast phosphate energy buffer to replenish ATP, osmolyte, antioxidant, neuromodulator, and as a compound with anabolic and ergogenic properties in muscle. Cr is taken from the diet or endogenously synthetized by the enzymes AGAT and GAMT, and specifically taken up by the transporter SLC6A8. Loss-of-function mutations in the genes encoding for the enzymes or the transporter cause Cerebral Creatine Deficiency Syndromes (CCDS). CCDS are characterized by brain Cr deficiency, intellectual disability with severe speech delay, behavioral troubles, epilepsy and motor dysfunction. Among CCDS, the X-linked Cr transporter deficiency (CTD) is the most prevalent with no efficient treatment so far. Different animal models of CTD show reduced brain Cr levels, cognitive deficiencies and together they cover other traits similar to those of patients. However, motor function was poorly explored in CTD models and some controversies in the phenotype exist in comparison with CTD patients. Our recently described Slc6a8Y389C knock-in (KI) rat model of CTD showed mild impaired motor function, morphological alterations in cerebellum, reduced muscular mass, Cr deficiency and increased guanidinoacetate content in muscle, although no consistent signs of muscle atrophy. Our results indicate that such motor dysfunction co-occurred with both nervous and muscle dysfunction, suggesting that muscle strength and performance as well as neuronal connectivity might be affected by this Cr deficiency in muscle and brain.

category: pathology


title: Creatine transporter deficient rat model show motor dysfunction, cerebellar alterations and muscle creatine deficiency without muscle atrophy


title: Prion-like α-synuclein pathology in the brains of infants: Krabbe disease as a novel seed-competent α-synucleinopathy

abstract: Krabbe disease (KD) is an infantile neurodegenerative disorder resulting from pathogenic variants in the GALC gene which causes accumulation of the toxic sphingolipid psychosine. GALC variants are associated with increased risk of Lewy body diseases (LBD), an umbrella term for age-associated neurodegenerative diseases in which the protein -synuclein aggregates into Lewy bodies. To explore whether -synuclein in KD has pathological similarities to that in LBD, we compared post-mortem KD tissue to that of infant control cases and identified alterations to -synuclein localisation and expression of modifications associated with LBD. To determine whether -synuclein in KD displayed pathogenic properties associated with LBD we evaluated its seeding capacity using the real-time quaking-induced conversion assay. Strikingly, seeded aggregation of -synuclein resulted in the formation of fibrillar aggregates similar to those observed in LBD, confirming the prion-like capacity of KD-derived -synuclein. These observations constitute the first report of prion-like -synuclein in the brain tissue of infants and challenge the putative view that -synuclein pathology is merely an age-associated phenomenon, instead suggesting it can result from alterations to biological processes such as sphingolipid homeostasis. Our findings have important implications for understanding the mechanisms underlying Lewy body formation in LBD.


title: Microscopic Quantification of Oxygen Consumption across Cortical Layers

abstract: The cerebral cortex is organized in cortical layers that differ in their cellular density, composition, and wiring. Cortical laminar architecture is also readily revealed by staining for cytochrome oxidase - the last enzyme in the respiratory electron transport chain located in the inner mitochondrial membrane. It has been hypothesized that a high-density band of cytochrome oxidase in cortical layer IV reflects higher oxygen consumption under baseline (unstimulated) conditions. Here, we tested the above hypothesis using direct measurements of the partial pressure of O2 (pO2) in cortical tissue by means of 2-photon phosphorescence lifetime microscopy (2PLM). We revisited our previously developed method for extraction of the cerebral metabolic rate of O2 (CMRO2) based on 2-photon pO2 measurements around diving arterioles and applied this method to estimate baseline CMRO2 in awake mice across cortical layers. To our surprise, our results revealed a decrease in baseline CMRO2 from layer I to layer IV. This decrease of CMRO2 with cortical depth was paralleled by an increase in tissue oxygenation. Higher baseline oxygenation and cytochrome density in layer IV may serve as an O2 reserve during surges of neuronal activity or certain metabolically active brain states rather than baseline energy needs. Our study provides the first quantification of microscopically resolved CMRO2 across cortical layers as a step towards informed interpretation and modeling of cortical-layer-specific Blood Oxygen Level Dependent (BOLD) functional Magnetic Resonance Imaging (fMRI) signals.


title: Deep learning models of cognitive processes constrained by human brain connectomes

abstract:
Brain decoding aims to infer cognitive states from recordings of brain activity. Current literature has mainly focused on isolated brain regions engaged in specific experimental conditions, but ignored the integrative nature of cognitive processes recruiting distributed brain networks. To tackle this issue, we propose a connectome-based graph neural network to integrate distributed patterns of brain activity in a multiscale manner, ranging from localized brain regions, to a specific brain circuit/network and towards the full brain. We evaluate the decoding model using a large task-fMRI database from the human connectome project. By implementing connectomic constraints and multiscale interactions in deep graph convolutions, the model achieves high accuracy of decoding 21 cognitive states (93%, chancel level: 4.8%) and shows high robustness against adversarial attacks on the graph architecture. Our study bridges human connectomes with deep learning techniques and provides new avenues to study the underlying neural substrates of human cognition at scale.

HighlightsO_LIA full-brain integrative model is critical for cognitive decoding
C_LIO_LIHuman connectomes and long-range connections accelerate the propagation of activity
C_LIO_LIChebNet decoding is robust to random attacks on brain connectomes and regions
C_LI


title: PyMVPD: A toolbox for multivariate pattern dependence

abstract: AO_SCPLOWBSTRACTC_SCPLOWCognitive tasks engage multiple brain regions. Studying how these regions interact is key to understand the neural bases of cognition. Standard approaches to model the interactions between brain regions rely on univariate statistical dependence. However, newly developed methods can capture multivariate dependence. Multivariate Pattern Dependence (MVPD) is a powerful and flexible approach that trains and tests multivariate models of the interactions between brain regions using independent data. In this article, we introduce PyMVPD: an open source toolbox for Multivariate Pattern Dependence. The toolbox includes pre-implemented linear regression models and artificial neural network models of the interactions between regions. It is designed to be easily customizable. We demonstrate example applications of PyMVPD using well-studied seed regions such as the fusiform face area (FFA) and the parahippocampal place area (PPA). Next, we compare the performance of different model architectures. Overall, artificial neural networks outperform linear regression. Importantly, the best performing architecture is region-dependent: MVPD subdivides cortex in distinct, contiguous regions whose interaction with FFA and PPA is best captured by different models.


title: A multi-dataset evaluation of frame censoring for task-based fMRI

abstract: Subject motion during fMRI can affect our ability to accurately measure signals of interest. In recent years, frame censoring--that is, statistically excluding motion-contaminated data within the general linear model using nuisance regressors--has appeared in several task-based fMRI studies as a mitigation strategy. However, there have been few systematic investigations quantifying its efficacy. In the present study, we compared the performance of frame censoring to several other common motion correction approaches for task-based fMRI using open data and reproducible workflows. We analyzed eight publicly-available datasets representing eleven distinct tasks in child, adolescent, and adult participants. Performance was quantified using maximum t-values in group analyses, and ROI-based mean activation and split-half reliability in single subjects. We compared frame censoring across several thresholds to the use of 6 and 24 canonical motion regressors, wavelet despiking, robust weighted least squares, and untrained ICA-based denoising, for a total of 240 separate analyses. Thresholds used to identify censored frames were based on both motion estimates (FD) and image intensity changes (DVARS). Relative to standard motion regressors, we found consistent improvements for modest amounts of frame censoring (e.g., 1-2% data loss), although these gains were frequently comparable to what could be achieved using other techniques. Importantly, no single approach consistently outperformed the others across all datasets and tasks. These findings suggest that the choice of a motion mitigation strategy depends on both the dataset and the outcome metric of interest.


title: A multi-dataset evaluation of frame censoring for motion correction in task-based fMRI


title: Aperiodic EEG activity masks the dynamics of neural oscillations during loss of consciousness from propofol

abstract: EEGs are known to provide biomarkers for consciousness. Although EEG correlates of loss of consciousness (LOC) are often ascribed to changes in neural synchrony, mounting evidence suggests that some changes result from asynchronous neural activity. By combining EEG recordings of humans undergoing propofol administration with biophysical modelling, we present here a principled decomposition of EEG changes during LOC into synchronous and asynchronous sources. Our results reveal that IPSP decay rate and mean spike rate shape aperiodic EEG features, and that propofols effects on these parameters largely explain the changes in EEG spectra following propofol infusion. We further show that traditional spectral EEG analysis likely conflates these effects with changes in rhythmic activity, thereby masking the true dynamics of neural synchrony. We conclude that the well-documented propofol-induced alpha rhythm in fact appears before LOC, and that the moment of LOC is uniquely correlated with the sudden appearance of a delta rhythm.


title: A Proposed Unified Mitotic Chromosome Architecture

abstract:
A molecular architecture is proposed for an example mitotic chromosome, human Chromosome 10. This architecture is built on a previously described interphase chromosome structure based on Cryo-EM cellular tomography (1), thus unifying chromosome structure throughout the complete mitotic cycle. The basic organizational principle, for mitotic chromosomes, is specific coiling of the 11-nm nucleosome fiber into large scale approximately 200 nm structures (a Slinky (2, motif cited in 3) in interphase, and then further modification and subsequent additional coiling for the final structure. The final mitotic chromosome architecture accounts for the dimensional values as well as the well known cytological configurations. In addition, proof is experimentally provided, by digital PCR technology, that G1 T-cell nuclei are diploid, thus one DNA molecule per chromosome. Many nucleosome linker DNA sequences, the promotors and enhancers, are suggestive of optimal exposure on the surfaces of the large-scale coils.

Significance StatementThe significance of this proposed mitotic chromosome architecture is that a specific, sequenced chromosome, human Chromosome 10, can be built into a specific architecture that accounts for the dimensional values and cytological descriptions, a first time result. Since this molecular architecture is an extension of the interphase chromosome structure, a coiling of the 11-nm nucleosome fiber with further coiling, a unifying molecular structure motif is present throughout the entire mitotic cycle, interphase through mitosis.


title: A multifaceted ecological assessment reveals the invasion of the freshwater red macroalga Montagnia macrospora (Batrachospermales, Rhodophyta) in Taiwan

abstract: Invasive freshwater macroalgae are rarely described. Recently, a freshwater red alga, Montagnia macrospora, was introduced from South America to East Asia via the global aquarium trade. The earliest occurrence record of this alga in Taiwan is dated 2005. To determine whether M. macrospora has become an invasive species in Taiwan and to understand its traits that facilitated its invasion, we took a total-evidence approach that combines examination of ecological background and population genetic analysis. Our island-wide survey showed that M. macrospora alga was widespread in the field across Taiwan, where the climate greatly differs from that of South America. Our population genetic analysis revealed that the cox2-3 sequences of all the specimens of M. macrospora from Taiwan were identical, consistent with the hypothesis that the alga expanded through asexual reproduction. Moreover, during our long-term ecological assessments and field surveys, we observed that M. macrospora is an ecological generalist that can self-sustain for a decade and bloom. Taken together, our data suggest that M. macrospora has successfully invaded the freshwater ecosystems in Taiwan due to its ability to disperse asexually and to grow under broad environmental conditions. We hope that our study brings attention to invasive freshwater algae, which have been overlooked in conservation planning and management.


title: New cycle, same old mistakes? Overlapping vs. discrete generations in long-term recurrent selection

abstract:
BackgroundRecurrent selection is a foundational breeding method for quantitative trait improvement. It typically features rapid breeding cycles that can lead to high rates of genetic gain. In recurrent phenotypic selection, generations do not overlap, which means that breeding candidates are evaluated and considered for selection for only one cycle. With recurrent genomic selection, candidates can be evaluated based on genomic estimated breeding values indefinitely, therefore facilitating overlapping generations. Candidates with true high breeding values that were discarded in one cycle due to underestimation of breeding value could be identified and selected in subsequent cycles. The consequences of allowing generations to overlap in recurrent selection are unknown. We assessed whether maintaining overlapping and discrete generations led to differences in genetic gain for phenotypic, genomic truncation, and genomic optimum contribution recurrent selection by simulation of traits with various heritabilities and genetic architectures across fifty breeding cycles. We also assessed differences of overlapping and discrete generations in a conventional breeding scheme with multiple stages and cohorts.

ResultsWith phenotypic selection, overlapping generations led to decreased genetic gain compared to discrete generations due to increased selection error bias. Selected individuals, which were in the upper tail of the distribution of phenotypic values, tended to also have high absolute error relative to their true breeding value compared to the overall population. Without repeated phenotyping, these individuals erroneously believed to have high value were repeatedly selected across cycles, leading to decreased genetic gain. With genomic truncation selection, overlapping and discrete generations performed similarly as updating breeding values precluded repeatedly selecting individuals with inaccurately high estimates of breeding values in subsequent cycles. Overlapping generations did not outperform discrete generations in the presence of a positive genetic trend with genomic truncation selection, as past generations had lower mean genetic values than the current generation of selection candidates. With genomic optimum contribution selection, overlapping and discrete generations performed similarly, but overlapping generations slightly outperformed discrete generations in the long term if the targeted inbreeding rate was extremely low.

ConclusionsMaintaining discrete generations in recurrent phenotypic selection leads to increased genetic gain, especially at low heritabilities, by preventing selection error bias. With genomic truncation selection and genomic optimum contribution selection, genetic gain does not differ between discrete and overlapping generations assuming non-genetic effects are not present. Overlapping generations may increase genetic gain in the long term with very low targeted rates of inbreeding in genomic optimum contribution selection.


title: From compartments to gene loops: Functions of the 3D genome in the human brain

abstract: The 3D genome plays a key role in the regulation of gene expression. However, little is known about the spatiotemporal organization of chromatin during human brain development. We investigated the 3D genome in human fetal cortical plate and in adult prefrontal cortical neurons and glia. We found that neurons have weaker compartments than glia that emerge during fetal development. Furthermore, neurons form loop domains whereas glia form compartment domains. We show through CRISPRi on CNTNAP2 that transcription is coupled to loop domain insulation. Gene regulation during neural development involves increased use of enhancer-promoter and repressor-promoter loops. Finally, transcription is associated with gene loops. Altogether, we provide novel insights into the relationship between gene expression and different scales of chromatin organization in the human brain.


title: Specificity of California mouse pup vocalizations in response to olfactory cues

abstract: In rodents, young pups communicate with their parents through harmonic calls and ultrasonic vocalizations (USVs). These forms of communication can improve chances of survival, since pups rely on their parents for thermoregulation, nutrition and protection. The extent to which pups modulate calls in response to their surroundings remains unclear. In this study we examined whether olfactory stimuli influence characteristics of pup calls, and how these calls may be affected by pup sex and litter size, in the California mouse (Peromyscus californicus). Pups were isolated and audio recorded during an initial, 3-minute control period, after which they were exposed for 5 minutes to bedding containing one of 4 olfactory cues: scent from their home cage, scent from the home cage of an unfamiliar family, coyote urine, or no scent. Latency to call, call rate, call duration and call characteristics (e.g. frequency and amplitude) were compared between the control period and olfactory-exposure period as well as among olfactory conditions. Pups from 2-pup litters called more quietly (lower amplitude) when exposed to odor from a predator while pups from 3-pup litters called louder (higher amplitude). Additionally, pups tended to reduce their call rates in response to odors from their home cage, consistent with contact quieting. However, pups tended to increase their rate of calling when exposed to predator urine, in contrast to the expectations of predator-induced vocal suppression. Lastly, male pups produced higher-frequency calls and more USVs than females. These results indicate that a number of pup call characteristics in this species can be influenced by acute olfactory stimuli as well as factors such as litter size and sex. The value of these pup call variations for offspring-parent communication is unclear: whether they elicit different parental responses is unknown and would be an interesting/valuable/informative avenue for future studies.


title: A weighted constraint satisfaction approach to human goal-directed decision making

abstract:
When we plan for long-range goals, proximal information cannot be exploited in a blindly myopic way, as relevant future information must also be considered. But when a subgoal must be resolved first, irrelevant future information should not interfere with the processing of more proximal, subgoal-relevant information. We explore the idea that decision making in both situations relies on the flexible modulation of the degree to which different pieces of information under consideration are weighted, rather than explicitly decomposing a problem into smaller parts and solving each part independently. We asked participants to find the shortest goal-reaching paths in mazes and modeled their initial path choices as a noisy, weighted information integration process. In a base task where choosing the optimal initial path required weighting starting-point and goal-proximal factors equally, participants did take both constraints into account, with participants who made more accurate choices tending to exhibit more balanced weighting. The base task was then embedded as an initial subtask in a larger maze, where the same two factors constrained the optimal path to a subgoal, and the final goal position was irrelevant to the initial path choice. In this more complex task, participants choices reflected predominant consideration of the subgoal-relevant constraints, but also some influence of the initially-irrelevant final goal. More accurate participants placed much less weight on the optimality-irrelevant goal and again tended to weight the two initially-relevant constraints more equally. These findings suggest that humans may rely on a graded, task-sensitive weighting of multiple constraints to generate approximately optimal decision outcomes in both hierarchical and non-hierarchical goal-directed tasks.

Author SummaryDifferent problems require the consideration of different information sources, including often useful long-range, future information that may impact our immediate decisions. However, when future information is irrelevant to a key subgoal, it can be desirable to focus on achieving the subgoal first. We suggest that humans rely on appropriately weighting relevant information over irrelevant information to generate decision outcomes in both types of situations. We conducted behavioral experiments and fitted models of decision processes to understand to what extent people considered various task factors in choosing the initial path in different mazes, both when a simple maze occurred alone or was embedded as an initial part in a larger maze. Our results show that people approximate the optimal decision outcomes in both tasks by modulating the weighting of different factors during planning, and that people who made more accurate initial path choices modulated these weightings more successfully than those who made less accurate choices.


title: gcFront: a tool for determining a Pareto front of growth-coupled cell factory designs

abstract:
MotivationA widely applicable strategy to create cell factories is to knock out (KO) genes or reactions to redirect cell metabolism so that chemical synthesis is made obligatory when the cell grows at its maximum rate. Synthesis is thus growth-coupled, and the stronger the coupling the more deleterious any impediments in synthesis are to cell growth, making high producer phenotypes evolutionarily robust. Additionally, we desire that these strains grow and synthesise at high rates. Genome-scale metabolic models can be used to explore and identify KOs that growth-couple synthesis, but these are rare in an immense design space, making the search difficult and slow.

ResultsTo address this multi-objective optimization problem, we developed a software tool named gcFront - using a genetic algorithm it explores KOs that maximise cell growth, product synthesis, and coupling strength. Moreover, our measure of coupling strength facilitates the search so that gcFront not only finds a growth coupled design in minutes but also outputs many alternative Pareto optimal designs from a single run - granting users flexibility in selecting designs to take to the lab.

Availability and ImplementationgcFront, with documentation and a workable tutorial, is freely available at GitHub: https://github.com/lLegon/gcFront, the repository of which is archived at Zenodo, DOI: 10.5281/zenodo.6338595 (Legon et al., 2022).

Supplementary InformationSupplementary notes and data files are available at Bioinformatics online.


title: netCRS: Network-based comorbidity risk score for prediction of myocardial infarction using biobank-scaled PheWAS data

abstract: The polygenic risk score (PRS) can help to identify individuals genetic susceptibility for various diseases by combining patient genetic profiles and identified single-nucleotide polymorphisms (SNPs) from genome-wide association studies. Although multiple diseases will usually afflict patients at once or in succession, conventional PRSs fail to consider genetic relationships across multiple diseases. Even multi-trait PRSs, which take into account genetic effects for more than one disease at a time, fail to consider a sufficient number of phenotypes to accurately reflect the state of disease comorbidity in a patient, or are biased in terms of the traits that are selected. Thus, we developed novel network-based comorbidity risk scores to quantify associations among multiple phenotypes from phenome-wide association studies (PheWAS). We first constructed a disease-SNP heterogeneous multi-layered network (DS-Net), which consists of a disease network (disease-layer) and SNP network (SNP-layer). The disease-layer describes the population-level interactome from PheWAS data. The SNP-layer was constructed according to linkage disequilibrium. Both layers were attached to transform the information from a population-level interactome to individual-level inferences. Then, graph-based semi-supervised learning was applied to predict possible comorbidity scores on disease-layer for each subject. The SNP-layer serves as receiving individual genotyping data in the scoring process, and the disease-layer serves as the propagated output for an individuals multiple disease comorbidity scores. The possible comorbidity scores were combined by logistic regression, and it is denoted as netCRS. The DS-Net was constructed from UK Biobank PheWAS data, and the individual genetic profiles were collected from the Penn Medicine Biobank. As a proof-of-concept study, myocardial infarction (MI) was selected to compare netCRS with the PRS with pruning and thresholding (PRS-PT). The combined model (netCRS + PRS-PT + covariates) achieved an AUC improvement of 6.26% compared to the (PRS-PT + covariates) model. In terms of risk stratification, the combined model was able to capture the risk of MI up to approximately eight-fold higher than that of the low-risk group. The netCRS and PRS-PT complement each other in predicting high-risk groups of patients with MI. We expect that using these risk prediction models will allow for the development of prevention strategies and reduction of MI morbidity and mortality.


title: Eoulsan 2: an efficient workflow manager for reproducible bulk, long-read and single-cell transcriptomics analyses

abstract:
AO_SCPLOWBSTRACTC_SCPLOWO_ST_ABSMotivationC_ST_ABSCore sequencing facilities produce huge amounts of sequencing data that need to be analysed with automated workflows to ensure reproducibility and traceability. Eoulsan is a versatile open-source workflow engine meeting the needs of core facilities, by automating the analysis of a large number of samples. Its core design separates the description of the workflow from the actual commands to be run. This originality simplifies its usage as the user does not need to handle code, while ensuring reproducibility. Eoulsan was initially developed for bulk RNA-seq data, but the transcriptomics applications have recently widened with the advent of long-read sequencing and single-cell technologies, calling for the development of new workflows.

ResultWe present Eoulsan 2, a major update that (i) enhances the workflow manager itself, (ii) facilitates the development of new modules, and (iii) expands its applications to long reads RNA-seq (Oxford Nanopore Technologies) and scRNA-seq (Smart-seq2 and 10x Genomics). The workflow manager has been rewritten, with support for execution on a larger choice of computational infrastructure (workstations, Hadoop clusters, and various job schedulers for cluster usage). Eoulsan now facilitates the development of new modules, by reusing wrappers developed for the Galaxy platform, with support for container images (Docker or Singularity) packaging tools to execute. Finally, Eoulsan natively integrates novel modules for bulk RNA-seq, as well as others specifically designed for processing long read RNA-seq and scRNA-seq. Eoulsan 2 is distributed with ready-to-use workflows and companion tutorials.

Availability and implementationEoulsan is implemented in Java, supported on Linux systems and distributed under the LGPL and CeCILL-C licenses at: http://outils.genomique.biologie.ens.fr/eoulsan/. The source code and sample workflows are available on GitHub: https://github.com/GenomicParisCentre/eoulsan. A GitHub repository for modules using the Galaxy tool XML syntax is further provided at: https://github.com/GenomicParisCentre/eoulsan-tools

Contacteoulsan@bio.ens.psl.eu


title: FAIRly big: A framework for computationally reproducible processing of large-scale data

abstract: Large-scale datasets present unique opportunities to perform scientific investigations with un-precedented breadth. However, they also pose considerable challenges for the findability, accessibility, interoperability, and reusability (FAIR) of research outcomes due to infrastructure limitations, data usage constraints, or software license restrictions. Here we introduce a DataLad-based, domain-agnostic framework suitable for reproducible data processing in compliance with open science mandates. The framework attempts to minimize platform idiosyncrasies and performance-related complexities. It affords the capture of machine-actionable computational provenance records that can be used to retrace and verify the origins of research outcomes, as well as be re-executed independent of the original computing infrastructure. We demonstrate the frameworks performance using two showcases: one highlighting data sharing and transparency (using the studyforrest.org dataset) and another highlighting scalability (using the largest public brain imaging dataset available: the UK Biobank dataset).


title: DeepImageTranslator V2: analysis of multimodal medical images using semantic segmentation maps generated through deep learning

abstract:
IntroductionAnalysis of multimodal medical images often requires the selection of one or many anatomical regions of interest (ROIs) for extraction of useful statistics. This task can prove laborious when a manual approach is used. We have previously developed a user-friendly software tool for image-to-image translation using deep learning. Therefore, we present herein an update to the DeepImageTranslator software with the addition of a tool for multimodal medical image segmentation analysis (hereby referred to as the MMMISA).

MethodsThe MMMISA was implemented using the Tkinter library; backend computations were implemented using the Pydicom, Numpy, and OpenCV libraries. We tested our software using 4188 slices from whole-body axial 2-deoxy-2-[18F]-fluoroglucose-position emission tomography/computed tomography scans ([18F]-FDG-PET/CT) of 10 patients from the American College of Radiology Imaging Network-Head and Neck Squamous Cell Carcinoma (ACRIN-HNSCC) database. Using the deep learning software DeepImageTranslator, a model was trained with 36 randomly selected CT slices and manually labelled semantic segmentation maps. Utilizing the trained model, all the CT scans of the 10 HNSCC patients were segmented with high accuracy. Segmentation maps generated using the deep convolutional network were then used to measure organ specific [18F]-FDG uptake. We also compared measurements performed using the MMMISA and those made with manually selected ROIs.

ResultsThe MMMISA is a tool that allows user to select ROIs based on deep learning-generated segmentation maps and to compute accurate statistics for these ROIs based on coregistered multimodal images. We found that organ-specific [18F]-FDG uptake measured using multiple manually selected ROIs is concordant with whole-tissue measurements made with segmentation maps using the MMMISA tool.


title: The role of cell geometry and cell-cell communication in gradient sensing

abstract:
Cells can measure shallow gradients of external signals to initiate and accomplish a migration or a morphogenetic process. Recently, starting from mathematical models like the local-excitation global-inhibition (LEGI) model and with the support of empirical evidence, it has been proposed that cellular communication improves the measurement of an external gradient. However, the mathematical models that have been used have over-simplified geometries (e.g., they are uni-dimensional) or assumptions about cellular communication, which limit the possibility to analyze the gradient sensing ability of more complex cellular systems. Here, we generalize the existing models to study the effects on gradient sensing of cell number, geometry and of long-versus short-range cellular communication in 2D systems representing epithelial tissues. We find that increasing the cell number can be detrimental for gradient sensing when the communication is weak and limited to nearest neighbour cells, while it is beneficial when there is long-range communication. We also find that, with long-range communication, the gradient sensing ability improves for tissues with more disordered geometries; on the other hand, an ordered structure with mostly hexagonal cells is advantageous with nearest neighbour communication. Our results considerably extend the current models of gradient sensing by epithelial tissues, making a step further toward predicting the mechanism of communication and its putative mediator in many biological processes.

Author summaryGroups of cells collectively migrate in many biological processes, ranging from development to cancer metastasis. The migration is often driven by the gradient of a signaling molecule that can be shallow and noisy, raising the question of how cells can measure it reliably. Cellular communication has recently been suggested to play a key role in gradient sensing, and mathematical models with simplified cellular geometries have been developed to help interpret and design experiments. In this work, we generalize the existing mathematical models to investigate how short- and long-range cellular communication can increase gradient sensing in two-dimensional models of epithelial tissues. We analyze various cellular geometries and tissue size, and we identify the optimal setting that corresponds to different types of communication. Our findings will help pinpoint the communication mechanisms at work in a given tissue and the properties of the molecules that mediate the communication.


title: A multi-omics approach identifies pancreatic cancer cell extracellular vesicles as mediators of the unfolded protein response in normal pancreatic epithelial cells.

abstract: Although cancer-derived extracellular vesicles (cEVs) are thought to play a pivotal role in promoting cancer progression events, their precise effect on neighboring normal cells is unknown. In this study, we investigated the impact of pancreatic cancer ductal adenocarcinoma (PDAC) derived EVs on recipient non-tumorigenic pancreatic normal epithelial cells upon internalization. We show that PDAC cEVs increase the proliferation and invasive capability of treated normal cells. We further demonstrate that cEVs induce endoplasmic reticulum (ER) stress and the unfolded protein response (UPR) in treated normal pancreatic epithelial cells within 24 hours. Subsequently, these cells release several inflammatory cytokines. Leveraging a layered multi-omics approach, we analyzed EV cargo from a panel of 6 PDAC and 2 normal pancreas cell lines, using multiple EV isolation methods. We found that cEVs were enriched for an array of biomolecules which can induce or regulate ER stress and the UPR, including palmitic acid, sphingomyelins, metabolic regulators of tRNA charging and proteins which regulate trafficking and degradation. We further show that palmitic acid, at doses relevant to those found in cEVs, is sufficient to induce ER stress in normal pancreas cells. These results suggest that cEV cargo packaging may be designed to disseminate proliferative and invasive characteristics upon internalization by distant recipient normal cells, hitherto unreported. This study is among the first to highlight a major role for PDAC cEVs to induce stress in treated normal pancreas cells that may modulate a systemic response leading to altered phenotypes. For the first time, our study implicates cEV transported palmitic acid as a potential driver in this process. These findings highlight the importance of EVs in mediating disease etiology and open potential areas of investigation toward understanding the role of cEV lipids in promoting cell transformation in the surrounding microenvironment.


title: Modular, robust and extendible multicellular circuit design in yeast

abstract: Division of labor between cells is ubiquitous in biology but the use of multi-cellular consortia for engineering applications is only beginning to be explored. A significant advantage of multi-cellular circuits is their potential to be modular with respect to composition but this claim has not yet been extensively tested using experiments and quantitative modeling. Here, we construct a library of 24 yeast strains capable of sending, receiving or responding to three molecular signals, characterize them experimentally and build quantitative models of their input-output relationships. We then compose these strains into two- and three-strain cascades as well a four-strain bistable switch and show that experimentally measured consortia dynamics can be predicted from the models of the constituent parts. To further explore the achievable range of behaviors, we perform a fully automated computational search over all two-, three- and four-strain consortia to identify combinations that realize target behaviors including logic gates, band-pass filters and time pulses. Strain combinations that are predicted to map onto a target behavior are further computationally optimized and then experimentally tested. Experiments closely track computational predictions. The high reliability of these model descriptions further strengthens the feasibility and highlights the potential for distributed computing in synthetic biology.


title: Benchmarking imputation methods for network inference using a novel method of synthetic scRNA-seq data generation

abstract: Single cell RNA-sequencing (scRNA-seq) has very rapidly become the new workhorse of modern biology providing an unprecedented global view on cellular diversity and heterogeneity. In particular, the structure of gene-gene expression correlation contains information on the underlying gene regulatory networks. However, interpretation of scRNA-seq data is challenging due to specific experimental error and biases that are unique to this kind of data including drop-out (or technical zeros). To deal with this problem several methods for imputation of zeros for scRNA-seq have been developed. However, it is not clear how these processing steps affect inference of genetic networks from single cell data. Here, we introduce Biomodelling.jl, a tool for generation of synthetic scRNA-seq data using multiscale modelling of stochastic gene regulatory networks in growing and dividing cells. Our tool produces realistic transcription data with a known ground truth network topology that can be used to benchmark different approaches for gene regulatory network inference. Using this tool we investigate the impact of different imputation methods on the performance of several network inference algorithms. Biomodelling.jl provides a versatile and useful tool for future development and benchmarking of network inference approaches using scRNA-seq data.

category: systems biology


title: DDX17 is required for efficient DSB repair at DNA:RNA hybrid deficient loci

abstract: Proteins with RNA-binding activity are increasingly being implicated in DNA damage responses (DDR). Additionally, DNA:RNA-hybrids are rapidly generated around DNA double-strand breaks (DSBs), and are essential for effective repair. Here, using a meta-analysis of proteomic data, we identify novel DNA repair proteins and characterise a novel role for DDX17 in DNA repair. We found DDX17 to be required for both cell survival and DNA repair in response to numerous agents that induce DSBs. Analysis of DSB repair factor recruitment to damage sites suggested a role for DDX17 early in the DSB ubiquitin cascade. Genome-wide mapping of R-loops revealed that while DDX17 promotes the formation of DNA:RNA-hybrids around DSB sites, this role is specific to loci that are naturally deficient for DNA:RNA-hybrids. We propose that DDX17 facilitates DSB repair at loci that are inefficient at forming DNA:RNA-hybrids by catalysing the formation of DSB-induced hybrids, thereby allowing propagation of the damage response.


title: Making Genomic Surveillance Deliver: A Lineage Classification and Nomenclature System to Inform Rabies Virus Elimination

abstract:
The availability of pathogen sequence data and use of genomic surveillance is rapidly increasing. Genomic tools and classification systems need updating to reflect this. Here, rabies virus is used as an example to showcase the potential value of updated genomic tools to enhance surveillance to better understand epidemiological dynamics and improve disease control. Previous studies have described the evolutionary history of rabies virus; however, the resulting taxonomy lacks the definition necessary to identify incursions, lineage turnover and transmission routes at high resolution. Here we propose a lineage classification system based on the dynamic nomenclature used for SARS-CoV-2, defining a lineage by phylogenetic methods, for tracking virus spread and comparing sequences across geographic areas. We demonstrate this system through application to the globally distributed Cosmopolitan clade of rabies virus, defining 73 total lineages within the clade, beyond the 22 previously reported. We further show how integration of this tool with a new rabies virus sequence data resource (RABV-GLUE) enables rapid application, for example, highlighting lineage dynamics relevant to control and elimination programmes, such as identifying importations and their sources, and areas of persistence and transmission, including transboundary incursions. This system and the tools developed should be useful for coordinating and targeting control programmes and monitoring progress as we work towards eliminating dog-mediated rabies, as well as having potential for broad application to the surveillance of other viruses.

Author SummaryThe importance of the ability to track the diversity and spread of viruses in a universal way that can be clearly communicated has been highlighted during the SARS-CoV-2 pandemic. This, accompanied with the increase in the availability and use of pathogen sequence data, means the development of new genomic tools and classification systems can strengthen outbreak response and disease control. Here, we present an easy-to-use objective and transferable classification tool for tracking viruses at high resolution. We use rabies virus, a neglected zoonotic disease that causes around 59,000 human deaths each year, as an example use case of this tool. Applying our tool to a global clade of rabies virus, we find an over 200% increase in the definition at which we can classify the virus, allowing us to identify areas of persistence and transmission that were not previously apparent, and patterns of virus spread. Insights from the application of this tool should prove valuable in targeting vaccination campaigns and improving surveillance as countries work towards the elimination of dog-mediated rabies.


title: The Hair Cell Analysis Toolbox: A machine learning-based whole cochlea analysis pipeline

abstract: Auditory hair cells, the whole length of the cochlea, are routinely visualized using light microscopy techniques. It is common, therefore, for one to collect more data than is practical to analyze manually. There are currently no widely accepted tools for unsupervised, unbiased, and comprehensive analysis of cells in an entire cochlea. This represents a stark gap between image-based data and other tests of cochlear function. To close this gap, we present a machine learning-based hair cell analysis toolbox, for the analysis of whole cochleae, imaged with confocal microscopy. The software presented here allows the automation of common image analysis tasks such as counting hair cells, determining their best frequency, as well as quantifying single cell immunofluorescence intensities along the entire cochlear coil. We hope these automated tools will remove a considerable barrier in cochlear image analysis, allowing for more informative and less selective data analysis practices.