Loci for prediction of penicillin and tetracycline susceptibility in Neisseria gonorrhoeae: a genome wide association study BackgroundWhile Neisseria gonorrhoeae poses an urgent public health threat because of increasing antimicrobial resistance, much of the circulating population remains susceptible to historical treatment regimens. Point-of-care diagnostics that report susceptibility could allow for reintroduction of these regimens, but development of such diagnostics has been limited to ciprofloxacin, for which susceptibility can be predicted from a single locus.

MethodsWe assembled a dataset of 12,045 N. gonorrhoeae genomes with phenotypic resistance data for tetracycline (n = 3,611) and penicillin (n = 6,935). Using conditional genome wide association studies (GWAS), we sought to define genetic variants associated with susceptibility to penicillin and tetracycline. We evaluated the sensitivity and specificity of these variants for predicting susceptibility and non-resistance in our collection of gonococcal genomes.

FindingsIn our conditional penicillin GWAS, the presence of a genetic variant defined by a non-mosaic penA allele without an insertion at codon 345 was significantly associated with penicillin susceptibility and had the highest negative effect size of significant variants (p = 5.0 x 10-14, {beta} = -2.5). In combination with the absence of blaTEM, this variant predicted penicillin susceptibility with high specificity (99.8%) and modest sensitivity (36.7%). For tetracycline, the wild type allele at rpsJ codon 57, encoding valine, was significantly associated with tetracycline susceptibility (p = 5.6 x 10-16, {beta} = -1.6) after conditioning on the presence of tetM. The combination of rpsJ codon 57 allele and tetM absence predicted tetracycline susceptibility with high specificity (97.2%) and sensitivity (88.7%).

InterpretationAs few as two genetic loci can predict susceptibility to penicillin and tetracycline in N. gonorrhoeae with high specificity. Molecular point-of-care diagnostics targeting these loci have the potential to increase available treatments for gonorrhea.

FundingNational Institute of Allergy and Infectious Diseases, the National Science Foundation, and the Smith Family Foundation

Research in ContextO_ST_ABSEvidence before this studyC_ST_ABSWe searched PubMed with the terms "Neisseria gonorrhoeae" and "diagnostic" or "assay" plus "penicillin" or "tetracycline" for reports in any language published up to July 1, 2021. We additionally searched for "Neisseria gonorrhoeae" and "genome wide association study". We found that previously proposed molecular diagnostics for penicillin and tetracycline susceptibility either exclusively focused on plasmid-mediated resistance (i.e., targeting blaTEM or tetM) or did not include variants in genes encoding antibiotic targets (e.g., did not include penA or rpsJ). Targets for molecular surveillance have focused on resistance-associated alleles rather than susceptibility-associated alleles. We did not find any previous penicillin or tetracycline conditional genome wide association studies (GWAS) in N. gonorrhoeae.

Added value of this studyTo identify targets for molecular diagnostics that predict penicillin and tetracycline susceptibility, we conducted GWAS conditioning on the presence of plasmid-mediated resistance determinants to detect chromosomal loci with the highest association with susceptibility. We discovered a sequence (penA_01) that differentiates susceptible isolates from those with a resistance-associated insertion at codon 345 and from those with mosaic penA alleles, which is associated with penicillin susceptibility. We also found that rpsJ codon 57 was the chromosomal locus contributing the most to tetracycline susceptibility. The combination of these chromosomal loci and the absence of plasmid encoded determinants predicts penicillin and tetracycline susceptibility with high specificity in both a large global collection of N. gonorrhoeae and a validation dataset consisting of recently published genomes from CDCs Gonococcal Isolate Surveillance Program (GISP) surveillance collected in 2018.

Implications of all the available evidenceThe chromosomal loci penA_01 and rpsJ codon 57 in combination with plasmid loci blaTEM and tetM are candidates for the development of point-of-care molecular diagnostics for penicillin and tetracycline susceptibility. The loci may be combined with the currently available ciprofloxacin susceptibility diagnostics to predict susceptibility to multiple antibiotics. Additionally, our study suggests that conditional GWAS focused on variants associated with susceptibility may be a promising approach to identify minimal sets of loci for molecular diagnostics and surveillance.
REGEN-COV protects against viral escape in preclinical and human studies Monoclonal antibodies against SARS-CoV-2 are a clinically validated therapeutic option against COVID-19. As rapidly emerging virus mutants are becoming the next major concern in the fight against the global pandemic, it is imperative that these therapeutic treatments provide coverage against circulating variants and do not contribute to development of treatment emergent resistance. To this end, we investigated the sequence diversity of the spike protein and monitored emergence of minor virus variants in SARS-COV-2 isolates found in COVID-19 patients or identified from preclinical in vitro and in vivo studies. This study demonstrates that a combination of non-competing antibodies, REGEN-COV, not only provides full coverage against current variants of concern/interest but also protects against emergence of new such variants and their potential seeding into the population in a clinical setting.
Strong immobilization of C. elegans en masse on cultivation plates by thermoelectric cooling for high-throughput submicron-resolution microscopy Despite its profound impact on biology, most high-resolution in vivo microscopy approaches remain low throughput because current immobilization techniques require significant manual effort. We greatly accelerate imaging of the nematode Caenorhabditis elegans by implementing a simple cooling approach to easily immobilize entire populations directly on their cultivation plates. We optimize and characterize cooling immobilization. Counterintuitively, relatively warmer temperatures immobilize animals significantly more effectively than colder temperatures utilized in prior studies. This enhanced immobilization enables clear submicron-resolution fluorescence imaging, which is challenging to achieve with most current immobilization techniques. We demonstrate 64x magnification 3D imaging and timelapse recording of neurons in adults and embryos without motion blur. Compared to standard azide immobilization, cooling immobilization reduces the animal preparation and recovery time by >98%, significantly increasing experimental speed. By obviating individual animal manipulation, our approach could also empower automated imaging of large C. elegans populations within standard experimental setups and workflows. F1F0 ATP Hydrolysis is a Determinant of Metabolic Rate, a Correlate of Lifespan, and a Weakness of Cancer Dramatic rectal temperature decrease, in mice administered with a drug that specifically inhibits F1F0 ATP hydrolysis, suggests that F1F0 ATP hydrolysis is a major determinant of metabolic rate in vivo. Across twelve investigated species, less F1F0 ATP hydrolysis correlates with greater maximal lifespan. Specific drug inhibition of F1F0 ATP hydrolysis exerts potent anti-cancer activity in vitro. A novel two-generation approach to the population dynamics of gregarious parasitoids Parasitoids, as important natural enemies, help maintain balance in natural ecosystems. Their population dynamics is generally predicted from the number of individuals. Here, using gregarious parasitoids as models, we show that this traditional approach omits one important parameter: the mothers manipulation of offspring fertility due to the clutch size-body size- fertility correlation. As a result of this correlation, when females deliberately adjust the number of offspring laid in a host, they determine not only the number but also the body sizes and reproductive potentials of those offspring. For the first time in the model species Anaphes flavipes, we determined the parasitoids offspring fertility from clutch size. Using this, we experimentally clarified the advantage of specific clutch size combinations and we show that identical fertility in the F1 generation can lead to distinctly different fertility values in the F2 generation. Even with the same number of hosts, lower fertility in the F1 generation can cause higher fertility in the F2 generation. Based on these results, we propose a novel two-generation approach which includes the clutch size-body size-fertility correlation. Our novel approach provides a new perspective for determining the individual fitness levels of gregarious parasitoids with new options for the modelling of parasitoid population dynamics. Genetic instability from a single S-phase after whole genome duplication Diploid and stable karyotypes are associated with health and fitness in animals. In contrast, whole genome duplications (WGDs) - doubling full chromosome content - are linked to genetic instability (GIN) and frequently found in human cancers 1-3. It has been established that WGDs fuel chromosome instability through abnormal mitosis 4-8, however, the immediate consequences of tetraploidy in the first interphase are not known. This is an essential question because single WGD events such as cytokinesis failure can promote tumorigenesis 9. Here, we found that newly born tetraploid human cells undergo high rates of DNA damage during DNA replication in the first S-phase. Using DNA combing and single cell sequencing, we show that DNA replication dynamics is perturbed, generating under- and over-replicated regions. Mechanistically, we found that these defects result from the lack of protein mass scaling up at the G1/S transition, which impairs the fidelity of DNA replication. This work shows that within a single interphase, unscheduled tetraploid cells can acquire highly abnormal karyotypes. These findings provide an explanation for the GIN landscape that favors tumorigenesis after tetraploidization. CRCM5484: A BET- BDII Selective Compound With Differential Anti-Leukemic Drug Modulation Differentially screening the Fr-PPIChem chemical library on the BET BRD4-BDII versus -BDI bromodomains led to the discovery of a BDII selective tetrahydropyridothienopyrimidinone (THPTP)-based compound. Structure-activity relationship (SAR) and hit-to-lead approaches allowed us to develop CRCM5484, a potent inhibitor of BET proteins with a preferential and 475-fold selectivity for the second bromodomain of the BRD3 protein (BRD3-BDII) over its first bromodomain (BRD3-BDI). Its very low activity was demonstrated in various cell-based assays, corresponding with recent data describing other selective BDII compounds. However, screening on a drug sensitivity and resistance-profiling platform revealed its ability to modulate the anti-leukemic activity in combination with various FDA-approved and/or in-development drugs in a cell- and context-dependent differential manner. Altogether, the results confirm the originality of the THPTP molecular mode of action in the BD cavity and its potential as starting scaffold for the development of potent and selective bromodomain inhibitors. Probing Patterning in Microbial Consortia with picCASO: a Cellular Automaton for Spatial Organisation Microbial consortia exhibit spatial patterning across diverse environments. Since probing the self-organization of natural microbial communities is limited by their inherent complexity, synthetic models have emerged as attractive alternatives. In this study, we develop novel frameworks of bacterial communication and explore the emergent spatiotemporal organization of microbes. Specifically, we build quorum sensing-mediated models of microbial growth that are utilized to characterize the dynamics of communities from arbitrary initial configurations and establish the effectiveness of our communication strategies in coupling the growth rates of microbes. Our simulations indicate that the behavior of quorum sensing-coupled consortia can be most effectively modulated by the rates of secretion of AHLs. Such a mechanism of control enables the construction of desired relative populations of constituent species in spatially organized populations. Our models accurately recapitulate previous experiments that have investigated pattern formation in synthetic multi-cellular systems. Additionally, our software tool enables the easy implementation and analysis of our frameworks for a variety of initial configurations and simplifies the development of sophisticated gene circuits facilitating distributed computing. Overall, we demonstrate the potential of spatial organization as a tunable parameter in synthetic biology by introducing a communication paradigm based on the location and strength of coupling of microbial strains.

Author SummaryInteracting microorganisms that coexist in a given environment tend to have well-defined spatial arrangements. While the emergence of such organization is seen across different microbiomes in nature, it is hitherto not well understood. Decoding the inherent spatial patterning of microbes is constrained by the complexity of their natural habitats. Here, we take advantage of synthetic models of microbial communities to study the dynamics of emergent spatial organization. Our framework of bacterial communication utilizes modular synthetic devices to couple growth rates. In addition to uncovering potential principles of spatial organization, this work enables the construction of complex genetic circuits distributed across communicating strains. Additionally, we have developed a software tool, picCASO, that streamlines the investigation of microbial systems communicating through such frameworks. Orientation dependent proton transverse relaxation in the human brain white matter: The magic angle effect on a cylindrical helix PurposeTo overcome some limitations of prior orientation-dependent proton transverse relaxation formalisms in white matter (WM) with a novel framework based on the generalized magic angle effect function.

MethodsA cylindrical helix model was developed embracing anisotropic rotational and translational diffusion of restricted molecules in human brain WM, with the former characterized by an axially symmetric system. Transverse relaxation rates R2 and [Formula] were divided into isotropic [Formula] and anisotropic parts, [Formula], with  denoting an open angle and{varepsilon} 0 an orientation ({Phi}) offset from DTI-derived primary diffusivity direction. The proposed framework (Fit A) was compared with prior models without{varepsilon} 0 on previously published water and methylene proton transverse relaxation rates from developing, healthy, and pathological WM at 3T. Goodness of fit was represented by root-mean-square error (RMSE). F-test and linear correlation were used with statistical significance set to P [&le;] 0.05.

ResultsFit A significantly (P<0.01) outperformed prior models as demonstrated by reduced RMSEs, e.g., 0.349 vs. 0.724 in myelin water. Fitted{varepsilon} 0 was in good agreement with calculated{varepsilon} 0 from directional diffusivities. Compared with those from healthy adult, the fitted [Formula], and  from neonates were substantially reduced but{varepsilon} 0 increased, consistent with incomplete myelination. Significant positive [Formula] and negative ( and [Formula]) correlations were found with aging (demyelination) in elderly.

ConclusionThe developed framework can better characterize orientation dependences from a wide range of proton transverse relaxation measurements in human brain WM, shedding new light on myelin microstructural alterations at the molecular level. An Extended Motif in the SARS-CoV-2 Spike Modulates Binding and Release of Host Coatomer in Retrograde Trafficking {beta}-Coronaviruses such as SARS-CoV-2 hijack coatomer protein-I (COPI) for spike protein retrograde trafficking to the progeny assembly site in endoplasmic reticulum-Golgi intermediate compartment (ERGIC). However, limited residue-level details are available into how the spike interacts with COPI. Here we identify an extended COPI binding motif in the spike that encompasses the canonical K-x-H dibasic sequence. This motif demonstrates selectivity for COPI subunit. Guided by an in silico analysis of dibasic motifs in the human proteome, we employ mutagenesis and binding assays to show that the spike motif terminal residues are critical modulators of complex dissociation, which is essential for spike release in ERGIC. COPI residues critical for spike motif binding are elucidated by mutagenesis and crystallography and found to be conserved in the zoonotic reservoirs, bats, pangolins, camels, and in humans. Collectively, our investigation on the spike motif identifies key COPI binding determinants with implications for retrograde trafficking. Clonal evolution in serially passaged Cryptococcus neoformans x deneoformans hybrids reveals a heterogenous landscape of genomic change Cryptococcus neoformans x deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high ([~]10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in subtelomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments. Predicting the zoonotic capacity of mammals to transmit SARS-CoV-2 Back and forth transmission of SARS-CoV-2 between humans and animals may lead to wild reservoirs of virus that can endanger efforts toward long-term control of COVID-19 in people, and protecting vulnerable animal populations that are particularly susceptible to lethal disease. Predicting high risk host species is key to targeting field surveillance and lab experiments that validate host zoonotic potential. A major bottleneck to predicting animal hosts is the small number of species with available molecular information about the structure of ACE2, a key cellular receptor required for viral cell entry. We overcome this bottleneck by combining species ecological and biological traits with 3D modeling of virus and host cell protein interactions using machine learning methods. This approach enables predictions about the zoonotic capacity of SARS-CoV-2 for over 5,000 mammals -- an order of magnitude more species than previously possible. The high accuracy predictions achieved by this approach are strongly corroborated by in vivo empirical studies. We identify numerous common mammal species whose predicted zoonotic capacity and close proximity to humans may further enhance the risk of spillover and spillback transmission of SARS-CoV-2. Our results reveal high priority areas of geographic overlap between global COVID-19 hotspots and potential new mammal hosts of SARS-CoV-2. With molecular sequence data available for only a small fraction of potential host species, predictive modeling integrating data across multiple biological scales offers a conceptual advance that may expand our predictive capacity for zoonotic viruses with similarly unknown and potentially broad host ranges. β-catenin and canonical Wnts in Hydra pattern formation: Insights from mathematical modelling into two distinct systems Formation of the body axes and the subsequent formation of the apical termini are fundamental steps during animal development. Nuclear {beta}-catenin and canonical Wnt3 have been identified as major players in both these processes in Hydra. Based on molecular knowledge of canonical Wnt signalling directly linking nuclear {beta}-catenin and Wnt3 activity, de novo axis formation and head formation have frequently been assumed to be part of the same pattern formation system. In this study, we investigated this relationship by combining new model simulations with available experimental results. We demonstrated that nuclear {beta}-catenin and Wnt3 likely contribute to two distinct de novo pattern formation systems in Hydra, organizing development and differentiation on two different spatial scales. In particular, the results suggested that nuclear {beta}-catenin acts at the scale of the whole body to control axis formation, whereas canonical Wnt3 signalling is involved in a downstream pathway responsible for small-scale patterning of the head. These results suggest that the underlying mechanisms may also be more complex than previously assumed in other animals where axis formation has been ascribed to canonical Wnt signalling. Model sharing in the human medial temporal lobe Effective planning involves knowing where different actions will take us. However natural environments are rich and complex, leading to an exponential increase in memory demand as a plan grows in depth. One potential solution to this problem is to share the neural state transition functions used for planning between similar contexts. Here, we asked human participants to perform a sequential decision making task designed so that knowledge could be shared between some contexts but not others. Computational modelling showed that participants shared a model of state transitions between contexts where appropriate. fMRI data identified the medial temporal lobe as a locus for learning of state transitions, and within the same region, correlated BOLD patterns were observed in contexts where state transition information were shared. Finally, we show that the transition model is updated more strongly following the receipt of positive compared to negative outcomes, a finding that challenges conventional theories of planning which assume knowledge about our environment is updated independently of outcomes received. Together, these findings propose a computational and neural account of how information relevant for planning can be shared between contexts. Ribosome actuator via the SpoT/ppGpp pathway to mitigate gene overexpression burden Gene overexpression causes a burden on ribosomes, which decreases the expression of other genes and growth rate. Here, we engineer a genetic ribosome actuator that increases gene expression rate and growth rate by over 150% and 80%, respectively. The actuator expresses the hydrolysis domain of the SpoT enzyme (SpoTH) in a bacterial strain with high basal levels of ppGpp, thus increasing ribosome availability. We use this actuator in a feedforward control arrangement, wherein SpoTH is co-expressed with a gene of interest (GOI). Without SpoTH co-expression, activation of the GOI decreases the expression of a constitutive gene and growth rate by over 80% and 40%, respectively. By contrast, with SpoTH co-expression, the ribosome binding site (RBS) of SpoTH can be tuned to achieve inappreciable change of either the expression of the constitutive gene or growth rate when the GOI is activated to the same level. We finally employ the actuator to keep growth rate constant while expressing dCas9 to a level that would otherwise decrease growth rate by more than 40% due to toxicity. This actuator will be relevant in applications where high cellular protein concentrations are desired without growth defects or global reduction of gene expression. Diverse stem-chondrichthyan oral structures and evidence for an independently acquired acanthodid dentition The teeth of sharks famously form a series of transversely-organised files with conveyor-belt replacement that are borne directly on the jaw cartilages, in contrast to the dermal plate-borne dentition of bony fishes that undergoes site-specific replacement. A major obstacle in understanding how this system evolved is the poorly understood relationships of the earliest chondrichthyans and the profusion of morphologically and terminologically diverse bones, cartilages, splints and whorls that they possess. Here we use tomographic methods to investigate mandibular structures in several early branching  acanthodian-grade stem-chondrichthyans. We show that the dentigerous jaw bones of disparate genera of ischnacanthids are united by a common construction, being growing bones with non-shedding dentition. Mandibular splints, which support the ventro-lateral edge of the Meckels cartilage in some taxa, are formed from dermal bone and may be an acanthodid synapomorphy. We demonstrate that the teeth of Acanthodopsis are borne directly on the mandibular cartilage and that this taxon is deeply nested within an edentulous radiation, representing an unexpected independent origin of teeth. Many or even all of the range of unusual oral structures may be apomorphic, but they should nonetheless be considered when building hypotheses of tooth and jaw evolution, both in chondrichthyans and more broadly. A revised view on the evolution of glutamine synthetase isoenzymes in plants Glutamine synthetase (GS) is a key enzyme responsible for the incorporation of inorganic nitrogen in the form of ammonium into the amino acid glutamine. The genes encoding GS are among the oldest existing genes in living organisms. In plants, two groups of functional GS enzymes are found: eubacterial GSIIb (GLN2) and eukaryotic GSIIe (GLN1/GS). Phylogenetic analyses have shown that the GLN2 group originated from bacteria following horizontal gene transfer. Only GLN1/GS genes are found in vascular plants, which suggests that they are involved in the final adaptation of plants to terrestrial life. The present phylogenetic study reclassifies the different GS of seed plants into three clusters: GS1a, GS1b and GS2. The presence of genes encoding GS2 has been expanded to Cycadopsida gymnosperms, which suggests the origin of this gene in a common ancestor of Cycadopsida, Ginkgoopsida and angiosperms. GS1a genes have been identified in all gymnosperms, basal angiosperms and some Magnoliidae species. Previous studies in conifers and the gene expression profiles obtained in ginkgo and magnolia in the present work could explain the absence of GS1a in more recent angiosperm species (e.g., monocots and eudicots) due to the redundant roles of GS1a and GS2 in photosynthetic cells. Altogether, the results provide a better understanding of the evolution of plant GS isoenzymes and their physiological roles, which is valuable for improving crop nitrogen use efficiency and productivity. Simulation of foraging behavior using a decision-making agent with Bayesian and inverse Bayesian inference: Temporal correlations and power laws in displacement patterns It has been stated that in human migratory behavior, the step length series may have temporal correlation and that there is some relationship between this time dependency and the fact that the frequency distribution of step length follows a power-law distribution. Furthermore, the frequency of occurrence of the step length in some large marine organisms has been found to switch between power-law and exponential distributions, depending on the difficulty of prey acquisition. However, to date it has not been clarified how the aforementioned three phenomena arise: the positive correlation created in the step length series, the relation between the positive correlation of the step length series and the form of an individuals step length distribution, and the switching between power-law and exponential distributions depending on the abundance of prey. This study simulated foraging behavior using the Bayesian decision-making agent simultaneously performing both knowledge learning and knowledge-based inference to analyze how the aforementioned three phenomena arise. In the agent with learning and inference, past experiences were stored as hypotheses (knowledge) and they were used in current foraging behavior; at the same time, the hypothesis continued to be updated based on new experiences. The simulation results show that the agent with both learning and inference has a mechanism that simultaneously causes all the phenomena. Nuclear NAD+-biosynthetic enzyme NMNAT1 facilitates survival of developing retinal neurons Despite mounting evidence that the mammalian retina is exceptionally reliant on proper NAD+ homeostasis for health and function, the specific roles of subcellular NAD+ pools in retinal development, maintenance, and disease remain obscure. Here, we show that deletion of the nuclear-localized NAD+ synthase nicotinamide mononucleotide adenylyltransferase-1 (NMNAT1) in the developing murine retina causes early and severe degeneration of photoreceptors and select inner retinal neurons via multiple distinct cell death pathways. This severe phenotype is associated with disruptions to retinal central carbon metabolism, purine nucleotide synthesis, and amino acid pathways. Furthermore, large-scale transcriptomics reveals dysregulation of a collection of photoreceptor and synapse-specific genes in NMNAT1 knockout retinas prior to detectable morphological or metabolic alterations. Collectively, our study reveals previously unrecognized complexity in NMNAT1-associated retinal degeneration and suggests a yet-undescribed role for NMNAT1 in gene regulation during photoreceptor terminal differentiation. Protein sequence sampling and prediction from structural data The increasing integration between protein engineering and machine learning has led to many interesting results. A problem still to solve is to evaluate the likelihood that a sequence will fold into a target structure. This problem can be also viewed as sequence prediction from a known structure.

In the current work, we propose improvements in the recent architecture of Geometric Vector Perceptrons [1] in order to optimize the sampling of sequences from a known backbone structure. The proposed model differs from the original in that there is: (i) no updating in the vectorial embedding, only in the scalar one, (ii) only one layer of decoding. The first aspect improves the accuracy of the model and reduces the use of memory, the second allows for training of the model with several tasks without incurring data leakage.

We treat the trained classifier as an Energy-Based Model and sample sequences by sampling amino acids in a non-autoreggresive manner in the empty positions of the sequence using energy-guided criteria and followed by random mutation optimization. We improve the median identity of samples from 40.2% to 44.7%.

An additional question worth investigating is whether sampled and original sequences fold into similar structures independent of their identity. We chose proteins in our test set whose sampled sequences show low identity (under 30%) but for which our model predicted favorable energies. We used AlphaFold [2, 3] and observed that the predicted structures for sampled sequences highly resemble the predicted structures for original sequences, with an average TM-score of 0.84. Neuroinvasion and neurotropism by SARS-CoV-2 variants in the K18-hACE2 mouse Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) not only affects the respiratory tract but also causes neurological symptoms such as loss of smell and taste, headache, fatigue or severe cerebrovascular complications. Using transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) we investigated the spatiotemporal distribution and pathomorphological features in the CNS following intranasal infection with SARS-CoV-2 variants, also after prior influenza A virus infection. Apart from Omicron, we found all variants to frequently spread to and within the CNS. Infection was restricted to neurons and appeared to spread from the olfactory bulb mainly in basally orientated regions in the brain and into the spinal cord, independent of ACE2 expression and without evidence of neuronal cell death, axonal damage or demyelination. However, microglial activation, microgliosis and a mild macrophage and T cell dominated inflammatory response was consistently observed, accompanied by apoptotic death of endothelial, microglial and immune cells, without their apparent infection. Microgliosis and immune cell apoptosis indicate a potential role of microglia for pathogenesis and viral effect in COVID-19 and possible impairment of neurological functions, especially in long COVID. These data may also be informative for the selection of therapeutic candidates, and broadly support investigation of agents with adequate penetration into relevant regions of the CNS. Enrichment analyses identify shared associations for 25 quantitative traits in over 600,000 individuals from seven diverse ancestries Since 2005, genome-wide association (GWA) datasets have been largely biased toward sampling European ancestry individuals, and recent studies have shown that GWA results estimated from self-identified European individuals are not transferable to non-European individuals due to various confounding challenges. Here, we demonstrate that enrichment analyses which aggregate SNP-level association statistics at multiple genomic scales--from genes to genomic regions and pathways--have been underutilized in the GWA era and can generate biologically interpretable hypotheses regarding the genetic basis of complex trait architecture. We illustrate examples of the robust associations generated by enrichment analyses while studying 25 continuous traits assayed in 566,786 individuals from seven diverse self-identified human ancestries in the UK Biobank and the Biobank Japan, as well as 44,348 admixed individuals from the PAGE consortium including cohorts of African-American, Hispanic and Latin American, Native Hawaiian, and American Indian/Alaska Native individuals. We identify 1,000 gene-level associations that are genome-wide significant in at least two ancestry cohorts across these 25 traits, as well as highly conserved pathway associations with triglyceride levels in European, East Asian, and Native Hawaiian cohorts. Single-domain antibodies efficiently neutralize SARS-CoV-2 variants of concern, including Omicron variant Virus-neutralizing antibodies are one of the few treatment options for COVID-19. The evolution of SARS-CoV-2 virus has led to the emergence of virus variants with reduced sensitivity to some antibody-based therapies. The development of potent antibodies with a broad spectrum of neutralizing activity is urgently needed. Here we isolated a panel of single-domain antibodies that specifically bind to the receptor-binding domain of SARS-CoV-2 S glycoprotein. Three of the selected antibodies exhibiting most robust neutralization potency were used to generate dimeric molecules. We observed that these modifications resulted in up to a 200-fold increase in neutralizing activity. The most potent heterodimeric molecule efficiently neutralized each of SARS-CoV-2 variant of concern, including Alpha, Beta, Gamma, Delta and Omicron variants. This heterodimeric molecule could be a promising drug candidate for a treatment for COVID-19 caused by virus variants of concern. Repeated subfunctionalization of a modular antimicrobial peptide gene for neural function Antimicrobial peptides (AMPs) are host-encoded antibiotics that combat invading pathogens. These genes commonly encode multiple products as post-translationally cleaved polypeptides. Recent studies have highlighted roles for AMPs in neurological contexts suggesting functions for these defence molecules beyond infection. During our immune study characterizing the antimicrobial peptide gene Baramicin, we recovered multiple Baramicin paralogs in Drosophila melanogaster and other species, united by their N-terminal IM24 domain. Not all paralogs were immune-induced. Here, through careful dissection of the Baramicin familys evolutionary history, we find that these non-immune paralogs result from repeated events of duplication and subsequent truncation of the coding sequence from an immune-inducible ancestor. These truncations leave only the IM24 domain as the prominent gene product. Surprisingly, using mutation and targeted gene silencing we demonstrate that two such genes are adapted for function in neural contexts in D. melanogaster. We also show enrichment in the head for independent Baramicin genes in other species. The Baramicin evolutionary history reveals that the IM24 Baramicin domain is not strictly useful in an immune context. We thus provide a case study for how an AMP-encoding gene might play dual roles in both immune and non-immune processes via its multiple peptide products. We reflect on these findings to highlight a blind spot in the way researchers approach AMP research in in vivo contexts.

Significance statementAntimicrobial peptides are immune proteins recently implicated in neurological roles. To date little attention has been paid to the contributions of different gene products in this function. Here we show that an antimicrobial peptide gene encodes multiple products with either immune-specific or neurological roles. Effect of a low concentration in chloride, sodium and potassium on oocyte maturation, oocyte hydration, ovulation and egg quality in rainbow trout. Water salinity is an important environmental factor known to have detrimental effects on salmonid reproduction, mostly when migrating female broodfish are held in sea water. In contrast, data obtained in freshwater are scarce and the impact of low water salinity, and more specifically of low water concentrations in sodium, chloride and potassium, during reproduction in freshwater is currently unknown. For this reason, and because ion and water fluxes are critical for the final steps of the female gamete formation, including oocyte hydration and ovulation, the aim of the present study was to investigate the impact of low salinity water on final oocyte maturation, ovulation and, ultimately, on egg quality, using rainbow trout as a physiological model and relevant aquaculture species.

Fish from the same commercial strain were raised either in a site characterized by low concentrations of Na+, K+, and Cl- ions in the water or in a closely located control site exhibiting higher concentration in these elements. Egg quality and duration of final oocyte maturation were investigated using innovative phenotyping tools such as automatic assessment of egg viability using the VisEgg system and non-invasive echograph-based monitoring of final oocyte maturation duration, respectively. Oocyte hydration during final oocyte maturation and after ovulation was also investigated. Finally, molecular phenotyping was performed using real-time PCR-based monitoring of several key players of final oocyte maturation and ovulation associated with ion and water transport, inflammation, proteolytic activity, and coagulation. Oocyte hydration and gene expression data were analyzed in the light of the duration of final oocyte maturation.

Here we show that low water salinity (i.e., low water concentration in chloride, sodium and potassium) negatively influences final oocyte maturation, ovulation and, ultimately, egg quality. Low water salinity triggered delayed ovulation and lower oocyte viability. When investigating the impact of low water salinity on final oocyte maturation duration, individuals presenting the most severe phenotypes exhibited impaired oocyte hydration and abnormally reduced gene expression levels of several key players of the ovulatory process. While the under expression of water (i.e., aquaporins) and ion (i.e., solute carriers) transporters is consistent with impaired oocyte hydration, our observations also indicate that the overall ovulatory gene expression program is disrupted. Our results raise the question of the mechanisms underlying the negative influence of low salinity water on the dynamics of the preovulatory phase, on the control of the oocyte homeostasis, including hydration, and on the overall success of the maturation-ovulation process.

HighlightsO_LILow water salinity impairs final oocyte maturation and egg quality in rainbow trout
C_LIO_LILow water salinity induces delayed ovulation and impaired oocyte hydration
C_LIO_LILow water salinity induces a dysregulation of several key ovulatory genes
C_LIO_LIMonitoring of final oocyte maturation can be performed using ultrasound staging
C_LI TeamTree analysis: a new approach to evaluate scientific production Advances in science and technology depend on the work of research teams and the publication of results through peer-reviewed articles representing a growing socio-economic resource. Current methods to mine the scientific literature regarding a field of interest focus on content, but the workforce credited by authorship remains largely unexplored, and appropriate measures of scientific production are debated. Here, a new bibliometric approach named TeamTree analysis is introduced that visualizes the development and composition of the workforce driving a field. A new citation-independent measure that scales with the H index estimates impact based on publication record, genealogical ties and collaborative connections. This author-centered approach complements existing tools to mine the scientific literature and to evaluate research across disciplines.