Search

As biological sequence databases continue growing, so do the insight that they promise to shed on the shape of the genetic diversity of life. However, to fulfil this promise the software must remain usable, be able to accommodate a large amount of data and allow use of modern high performance computing infrastructure. In this study we present a reimplementation as well as an extension of a technique using indicator vectors to compute and visualize similarities between sets of nucleotide sequences. We have a flexible and easy to use python program relying on standard and open-source libraries. Our tool allows analysis of very large complement of sequences using code parallelization, as well as by providing routines to split a computational task in smaller and manageable subtasks whose results are then merged. This implementation also facilitates adding new sequences into an indicator vector-based representation without re-computing the whole set. The efficient synthesis of data into knowledge is no trivial matter given the size and rapid growth of biological sequence databases. Based on previous results regarding the properties of indicator vectors, the open-source approach proposed here efficiently and flexibly supports comparative analysis of genetic diversity at a large scale. Our software is freely available at: https://github.com/WandrilleD/pyKleeBarcode., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Duchemin, Thaler. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages., (© 2023 Society for Experimental Biology and John Wiley & Sons Ltd.)

Growing evidence links COVID-19 with acute and long-term neurological dysfunction. However, the pathophysiological mechanisms resulting in central nervous system involvement remain unclear, posing both diagnostic and therapeutic challenges. Here we show outcomes of a cross-sectional clinical study (NCT04472013) including clinical and imaging data and corresponding multidimensional characterization of immune mediators in the cerebrospinal fluid (CSF) and plasma of patients belonging to different Neuro-COVID severity classes. The most prominent signs of severe Neuro-COVID are blood-brain barrier (BBB) impairment, elevated microglia activation markers and a polyclonal B cell response targeting self-antigens and non-self-antigens. COVID-19 patients show decreased regional brain volumes associating with specific CSF parameters, however, COVID-19 patients characterized by plasma cytokine storm are presenting with a non-inflammatory CSF profile. Post-acute COVID-19 syndrome strongly associates with a distinctive set of CSF and plasma mediators. Collectively, we identify several potentially actionable targets to prevent or intervene with the neurological consequences of SARS-CoV-2 infection., (© 2022. The Author(s).)

Key innovations enable access to new adaptive zones and are often linked to increased species diversification. As such, innovations have attracted much attention, yet their concrete consequences on the subsequent evolutionary trajectory and diversification of the bearing lineages remain unclear. Water striders and relatives (Hemiptera: Heteroptera: Gerromorpha) represent a monophyletic lineage of insects that transitioned to live on the water-air interface and that diversified to occupy ponds, puddles, streams, mangroves and even oceans. This lineage offers an excellent model to study the patterns and processes underlying species diversification following the conquest of new adaptive zones. However, such studies require a reliable and comprehensive phylogeny of the infraorder. Based on whole transcriptomic datasets of 97 species and fossil records, we reconstructed a new phylogeny of the Gerromorpha that resolved inconsistencies and uncovered strong support for previously unknown relationships between some important taxa. We then used this phylogeny to reconstruct the ancestral state of a set of adaptations associated with water surface invasion (fluid locomotion, dispersal and transition to saline waters) and sexual dimorphism. Our results uncovered important patterns and dynamics of phenotypic evolution, revealing how the initial event of water surface invasion enabled multiple subsequent transitions to new adaptive zones on the water surfaces. This phylogeny and the associated transcriptomic datasets constitute highly valuable resources, making Gerromorpha an attractive model lineage to study phenotypic evolution., (© The Author(s) 2022. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

A patient-tailored, ex vivo drug response platform for glioblastoma (GBM) would facilitate therapy planning, provide insights into treatment-induced mechanisms in the immune tumor microenvironment (iTME), and enable the discovery of biomarkers of response. We cultured regionally annotated GBM explants in perfusion bioreactors to assess iTME responses to immunotherapy. Explants were treated with anti-CD47, anti-PD-1, or their combination, and analyzed by multiplexed microscopy [CO-Detection by indEXing (CODEX)], enabling the spatially resolved identification of >850,000 single cells, accompanied by explant secretome interrogation. Center and periphery explants differed in their cell type and soluble factor composition, and responses to immunotherapy. A subset of explants displayed increased interferon-γ levels, which correlated with shifts in immune cell composition within specified tissue compartments. Our study demonstrates that ex vivo immunotherapy of GBM explants enables an active antitumoral immune response within the tumor center and provides a framework for multidimensional personalized assessment of tumor response to immunotherapy.

Insecticide resistance in Anopheles mosquitoes is seriously threatening the success of insecticide-based malaria vector control. Surveillance of insecticide resistance in mosquito populations and identifying the underlying mechanisms enables optimisation of vector control strategies. Here, we investigated the molecular mechanisms of insecticide resistance in three Anopheles coluzzii field populations from southern Côte d'Ivoire, including Agboville, Dabou and Tiassalé. All three populations were resistant to bendiocarb, deltamethrin and DDT, but not or only very weakly resistant to malathion. The absence of malathion resistance is an unexpected result because we found the acetylcholinesterase mutation Ace1-G280S at high frequencies, which would typically confer cross-resistance to carbamates and organophosphates, including malathion. Notably, Tiassalé was the most susceptible population to malathion while being the most resistant one to the pyrethroid deltamethrin. The resistance ratio to deltamethrin between Tiassalé and the laboratory reference colony was 1,800 fold. By sequencing the transcriptome of individual mosquitoes, we found numerous cytochrome P450-dependent monooxygenases - including CYP6M2, CYP6P2, CYP6P3, CYP6P4 and CYP6P5 - overexpressed in all three field populations. This could be an indication for negative cross-resistance caused by overexpression of pyrethroid-detoxifying cytochrome P450s that may activate pro-insecticides, thereby increasing malathion susceptibility. In addition to the P450s, we found several overexpressed carboxylesterases, glutathione S-transferases and other candidates putatively involved in insecticide resistance., Competing Interests: The authors have declared that no competing interests exist.

The cell cycle is a fundamental process that has been extensively studied in bacteria. However, many of its components and their interactions with machineries involved in other cellular processes are poorly understood. Furthermore, most knowledge relies on the study of a few models, but the real diversity of the cell division apparatus and its evolution are largely unknown. Here, we present a massive in-silico analysis of cell division and associated processes in around 1,000 genomes of the Firmicutes, a major bacterial phylum encompassing models (i.e. Bacillus subtilis, Streptococcus pneumoniae, and Staphylococcus aureus), as well as many important pathogens. We analyzed over 160 proteins by using an original approach combining phylogenetic reconciliation, phylogenetic profiles, and gene cluster survey. Our results reveal the presence of substantial differences among clades and pinpoints a number of evolutionary hotspots. In particular, the emergence of Bacilli coincides with an expansion of the gene repertoires involved in cell wall synthesis and remodeling. We also highlight major genomic rearrangements at the emergence of Streptococcaceae. We establish a functional network in Firmicutes that allows identifying new functional links inside one same process such as between FtsW (peptidoglycan polymerase) and a previously undescribed Penicilin-Binding Protein or between different processes, such as replication and cell wall synthesis. Finally, we identify new candidates involved in sporulation and cell wall synthesis. Our results provide a previously undescribed view on the diversity of the bacterial cell cycle, testable hypotheses for further experimental studies, and a methodological framework for the analysis of any other biological system., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background: Barley is one of the founder crops of Neolithic agriculture and is among the most-grown cereals today. The only trait that universally differentiates the cultivated and wild subspecies is 'non-brittleness' of the rachis (the stem of the inflorescence), which facilitates harvesting of the crop. Other phenotypic differences appear to result from facultative or regional selective pressures. The population structure resulting from these regional events has been interpreted as evidence for multiple domestications or a mosaic ancestry involving genetic interaction between multiple wild or proto-domesticated lineages. However, each of the three mutations that confer non-brittleness originated in the western Fertile Crescent, arguing against multiregional origins for the crop., Results: We examined exome data for 310 wild, cultivated and hybrid/feral barley accessions and showed that cultivated barley is structured into six genetically-defined groups that display admixture, resulting at least in part from two or more significant passages of gene flow with distinct wild populations. The six groups are descended from a single founding population that emerged in the western Fertile Crescent. Only a few loci were universally targeted by selection, the identity of these suggesting that changes in seedling emergence and pathogen resistance could represent crucial domestication switches. Subsequent selection operated on a regional basis and strongly contributed to differentiation of the genetic groups., Conclusions: Identification of genetically-defined groups provides clarity to our understanding of the population history of cultivated barley. Inference of population splits and mixtures together with analysis of selection sweeps indicate descent from a single founding population, which emerged in the western Fertile Crescent. This founding population underwent relatively little genetic selection, those changes that did occur affecting traits involved in seedling emergence and pathogen resistance, indicating that these phenotypes should be considered as 'domestication traits'. During its expansion out of the western Fertile Crescent, the crop underwent regional episodes of gene flow and selection, giving rise to a modern genetic signature that has been interpreted as evidence for multiple domestications, but which we show can be rationalized with a single origin.

Motivation: A reconciliation is an annotation of the nodes of a gene tree with evolutionary events-for example, speciation, gene duplication, transfer, loss, etc.-along with a mapping onto a species tree. Many algorithms and software produce or use reconciliations but often using different reconciliation formats, regarding the type of events considered or whether the species tree is dated or not. This complicates the comparison and communication between different programs., Results: Here, we gather a consortium of software developers in gene tree species tree reconciliation to propose and endorse a format that aims to promote an integrative-albeit flexible-specification of phylogenetic reconciliations. This format, named recPhyloXML, is accompanied by several tools such as a reconciled tree visualizer and conversion utilities., Availability and Implementation: http://phylariane.univ-lyon1.fr/recphyloxml/.

Background: Genomes rearrangements carry valuable information for phylogenetic inference or the elucidation of molecular mechanisms of adaptation. However, the detection of genome rearrangements is often hampered by current deficiencies in data and methods: Genomes obtained from short sequence reads have generally very fragmented assemblies, and comparing multiple gene orders generally leads to computationally intractable algorithmic questions., Results: We present a computational method, ADSEQ, which, by combining ancestral gene order reconstruction, comparative scaffolding and de novo scaffolding methods, overcomes these two caveats. ADSEQ provides simultaneously improved assemblies and ancestral genomes, with statistical supports on all local features. Compared to previous comparative methods, it runs in polynomial time, it samples solutions in a probabilistic space, and it can handle a significantly larger gene complement from the considered extant genomes, with complex histories including gene duplications and losses. We use ADSEQ to provide improved assemblies and a genome history made of duplications, losses, gene translocations, rearrangements, of 18 complete Anopheles genomes, including several important malaria vectors. We also provide additional support for a differentiated mode of evolution of the sex chromosome and of the autosomes in these mosquito genomes., Conclusions: We demonstrate the method's ability to improve extant assemblies accurately through a procedure simulating realistic assembly fragmentation. We study a debated issue regarding the phylogeny of the Gambiae complex group of Anopheles genomes in the light of the evolution of chromosomal rearrangements, suggesting that the phylogenetic signal they carry can differ from the phylogenetic signal carried by gene sequences, more prone to introgression.

Background: Forming a new species through the merger of two or more divergent parent species is increasingly seen as a key phenomenon in the evolution of many biological systems. However, little is known about how expression of parental gene copies (homeologs) responds following genome merger. High throughput RNA sequencing now makes this analysis technically feasible, but tools to determine homeolog expression are still in their infancy., Results: Here we present HyLiTE - a single-step analysis to obtain tables of homeolog expression in a hybrid or allopolyploid and its parent species directly from raw mRNA sequence files. By implementing on-the-fly detection of diagnostic parental polymorphisms, HyLiTE can perform SNP calling and read classification simultaneously, thus allowing HyLiTE to be run as parallelized code. HyLiTE accommodates any number of parent species, multiple data sources (including genomic DNA reads to improve SNP detection), and implements a statistical framework optimized for genes with low to moderate expression., Conclusions: HyLiTE is a flexible and easy-to-use program designed for bench biologists to explore patterns of gene expression following genome merger. HyLiTE offers practical advantages over manual methods and existing programs, has been designed to accommodate a wide range of genome merger systems, can identify SNPs that arose following genome merger, and offers accurate performance on non-model organisms.

Background: We propose the computational reconstruction of a whole bacterial ancestral genome at the nucleotide scale, and its validation by a sequence of ancient DNA. This rare possibility is offered by an ancient sequence of the late middle ages plague agent. It has been hypothesized to be ancestral to extant Yersinia pestis strains based on the pattern of nucleotide substitutions. But the dynamics of indels, duplications, insertion sequences and rearrangements has impacted all genomes much more than the substitution process, which makes the ancestral reconstruction task challenging., Results: We use a set of gene families from 13 Yersinia species, construct reconciled phylogenies for all of them, and determine gene orders in ancestral species. Gene trees integrate information from the sequence, the species tree and gene order. We reconstruct ancestral sequences for ancestral genic and intergenic regions, providing nearly a complete genome sequence for the ancestor, containing a chromosome and three plasmids., Conclusion: The comparison of the ancestral and ancient sequences provides a unique opportunity to assess the quality of ancestral genome reconstruction methods. But the quality of the sequencing and assembly of the ancient sequence can also be questioned by this comparison.