Publications

Lowland tropical bryophytes have been perceived as excellent dispersers. In such groups, the inverse isolation hypothesis proposes that spatial genetic structure is erased beyond the limits of short-distance dispersal. Here, we determine the influence of environmental variation and geographic barriers on the spatial genetic structure of a widely dispersed and phylogenetically independent sample of Amazonian bryophytes. Single nucleotide polymorphism data were produced from a restriction site-associated DNA sequencing protocol for 10 species and analyzed through F-statistics and Mantel tests. Neither isolation-by-environment nor the impact of geographic barriers were recovered from the analyses. However, significant isolation-by-distance patterns were observed for 8 out of the 10 investigated species beyond the scale of short-distance dispersal (> 1 km), offering evidence contrary to the inverse isolation hypothesis. Despite a cadre of life-history traits and distributional patterns suggesting that tropical bryophytes are highly vagile, our analyses reveal spatial genetic structures comparable to those documented for angiosperms, whose diaspores are orders of magnitude larger. Dispersal limitation for tropical bryophytes flies in the face of traditional assumptions regarding their dispersal potential, and suggests that the plight of this component of cryptic biodiversity is more dire than previously considered in light of accelerated forest fragmentation in the Amazon.

Drosophila circadian behavior relies on the network of heterogeneous groups of clock neurons. Short- and long-range signaling within the pacemaker circuit coordinates molecular and neural rhythms of clock neurons to generate coherent behavioral output. The neurochemistry of circadian behavior is complex and remains incompletely understood. Here we demonstrate that the gaseous messenger nitric oxide (NO) is a signaling molecule linking circadian pacemaker to rhythmic locomotor activity. We show that mutants lacking nitric oxide synthase (NOS) have behavioral arrhythmia in constant darkness, although molecular clocks in the main pacemaker neurons are unaffected. Behavioral phenotypes of mutants are due in part to the malformation of neurites of the main pacemaker neurons, s-LNvs. Using cell-type selective and stage-specific gain- and loss-of-function of NOS, we also demonstrate that NO secreted from diverse cellular clusters affect behavioral rhythms. Furthermore, we identify the perineurial glia, one of the two glial subtypes that form the blood-brain barrier, as the major source of NO that regulates circadian locomotor output. These results reveal for the first time the critical role of NO signaling in the Drosophila circadian system and highlight the importance of neuro-glial interaction in the neural circuit output.

Proliferating animal cells are able to orient their mitotic spindles along their interphase cell axis, setting up the axis of cell division, despite rounding up as they enter mitosis. This has previously been attributed to molecular memory and, more specifically, to the maintenance of adhesions and retraction fibers in mitosis [1, 2, 3, 4, 5, 6], which are thought to act as local cues that pattern cortical Gαi, LGN, and nuclear mitotic apparatus protein (NuMA) [3, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. This cortical machinery then recruits and activates Dynein motors, which pull on astral microtubules to position the mitotic spindle. Here, we reveal a dynamic two-way crosstalk between the spindle and cortical motor complexes that depends on a Ran-guanosine triphosphate (GTP) signal [12], which is sufficient to drive continuous monopolar spindle motion independently of adhesive cues in flattened human cells in culture. Building on previous work [1, 12, 19, 20, 21, 22, 23], we implemented a physical model of the system that recapitulates the observed spindle-cortex interactions. Strikingly, when this model was used to study spindle dynamics in cells entering mitosis, the chromatin-based signal was found to preferentially clear force generators from the short cell axis, so that cortical motors pulling on astral microtubules align bipolar spindles with the interphase long cell axis, without requiring a fixed cue or a physical memory of interphase shape. Thus, our analysis shows that the ability of chromatin to pattern the cortex during the process of mitotic rounding is sufficient to translate interphase shape into a cortical pattern that can be read by the spindle, which then guides the axis of cell division.

The Clarion-Clipperton Zone (CCZ) occupies a vast swathe of the Pacific with extensive polymetallic nodule deposits. Eastern and central parts host diverse assemblages of xenophyophores (megafaunal agglutinated foraminifera). Here we describe xenophyophores obtained using a Remotely Operated Vehicle from the western CCZ. Eleven distinct forms include two known species, Stannophyllum zonarium Haeckel, 1888 and Aschemonella monile Gooday and Holzmann in Gooday et al., 2017b. Another four are described as new species based on morphological and genetic data. In Abyssalia foliformis gen. nov., sp. nov. and Abyssalia sphaerica sp. nov. the flattened or spherical test comprises a homogeneous framework of sponge spicules. Psammina tenuis sp. nov. has a delicate, thin, plate-like test. Moanammina semicircularis gen. nov., sp. nov. has a stalked, fan-shaped test and is genetically identical to 'Galatheammina sp. 6' of Gooday and co-workers from the eastern CCZ. Sequence data revealed a spherical 'mudball', which disintegrated and cannot be formally described, to be a novel xenophyophore. Finally, four morphospecies are represented by dead tests: Psammina spp., Reticulammina sp., and an unknown genus with a unique test structure. This collection enhances our knowledge of Pacific xenophyophore diversity and provides the first genetic confirmation of wide geographic ranges for abyssal species.

Restriction site Associated DNA Sequencing (RAD-Seq) is a technique characterized by the sequencing of specific loci along the genome that is widely employed in the field of evolutionary biology since it allows to exploit variants (mainly Single Nucleotide Polymorphism-SNPs) information from entire populations at a reduced cost. Common RAD dedicated tools, such as or , are based on all-vs-all read alignments, which require consequent time and computing resources. We present an original method, DiscoSnp-RAD, that avoids this pitfall since variants are detected by exploiting specific parts of the assembly graph built from the reads, hence preventing all-vs-all read alignments. We tested the implementation on simulated datasets of increasing size, up to 1,000 samples, and on real RAD-Seq data from 259 specimens of flies, morphologically assigned to seven species. All individuals were successfully assigned to their species using both STRUCTURE and Maximum Likelihood phylogenetic reconstruction. Moreover, identified variants succeeded to reveal a within-species genetic structure linked to the geographic distribution. Furthermore, our results show that DiscoSnp-RAD is significantly faster than state-of-the-art tools. The overall results show that DiscoSnp-RAD is suitable to identify variants from RAD-Seq data, it does not require time-consuming parameterization steps and it stands out from other tools due to its completely different principle, making it substantially faster, in particular on large datasets.

Coelacanths are iconic fishes represented today by a single marine genus. The group was a little bit more diversified in the Mesozoic, with representatives in marine and continental environments in the Late Cretaceous. Here we describe isolated skull bones of the last know freshwater coelacanths found in several fossil sites from the Early Campanian to the Early Maastrichtian of Southern France (in the Departments of Aude, Bouches-du-Rhône, Hérault, and Var). The sample does not allow distinguishing different species, and all material is referred to Axelrodichthys megadromos Cavin, Valentin, Garcia originally described from the locality of Ventabren in Southern France. A reconstruction of the skull is proposed. Previously unrecognized features are described, including parts of the postparietal portion of the skull, of the suspensorium and of the mandible. The new data confirm the assignation of the species to the mawsoniids, and more specifically to Axelrodichthys. A cladistic analysis scoring new character states provides a similar topology than a previous analysis, i.e. A. megadromos is placed in a polytomy with Axelrodichthys araripensis and Lualabaea lerichei, two species from the Early Cretaceous of Brazil and from the Late Jurassic of the Democratic Republic of the Congo, respectively. A. megadromos appears to have been restricted to freshwater environments, to the contrary of oldest Western Gondwanan representatives of the family that were able to live in brackish and marine waters. A. megadromos is the last representative of the mawsoniids and its occurrence in Europe is probably the result of a dispersal event from Western Gondwana that happened somewhen in the Cretaceous. Based on the available data, the mawsoniids went extinct in the mid-Maastrichthian, i.e. before the end-Cretaceous mass extinction. But it is possible that the fossil record of this family, which has been only recently recognized in Late Cretaceous European deposits, will geographically and stratigraphically widen with further discoveries.

Restriction site Associated DNA Sequencing (RAD-Seq) is a technique characterized by the sequencing of specific loci along the genome that is widely employed in the field of evolutionary biology since it allows to exploit variants (mainly Single Nucleotide Polymorphism-SNPs) information from entire populations at a reduced cost. Common RAD dedicated tools, such as STACKS or IPyRAD, are based on all-vs-all read alignments, which require consequent time and computing resources. We present an original method, DiscoSnp-RAD, that avoids this pitfall since variants are detected by exploiting specific parts of the assembly graph built from the reads, hence preventing all-vs-all read alignments. We tested the implementation on simulated datasets of increasing size, up to 1,000 samples, and on real RAD-Seq data from 259 specimens of Chiastocheta flies, morphologically assigned to seven species. All individuals were successfully assigned to their species using both STRUCTURE and Maximum Likelihood phylogenetic reconstruction. Moreover, identified variants succeeded to reveal a within-species genetic structure linked to the geographic distribution. Furthermore, our results show that DiscoSnp-RAD is significantly faster than state-of-the-art tools. The overall results show that DiscoSnp-RAD is suitable to identify variants from RAD-Seq data, it does not require time-consuming parameterization steps and it stands out from other tools due to its completely different principle, making it substantially faster, in particular on large datasets.

We report detailed peptide binding affinities between 438 HLA Class I and Class II proteins and complete proteomes of seven pandemic human viruses, including coronaviruses, influenza viruses and HIV-1. We contrast these affinities with HLA allele frequencies across hundreds of human populations worldwide. Statistical modelling shows that peptide binding affinities classified into four distinct categories depend on the HLA locus but that the type of virus is only a weak predictor, except in the case of HIV-1. Amongst the strong HLA binders (IC  ≤ 50), we uncovered 16 alleles (the top ones being A*02:02, B*15:03 and DRB1*01:02) binding more than 1% of peptides derived from all viruses, 9 (top ones including HLA-A*68:01, B*15:25, C*03:02 and DRB1*07:01) binding all viruses except HIV-1, and 15 (top ones A*02:01 and C*14:02) only binding coronaviruses. The frequencies of strongest and weakest HLA peptide binders differ significantly among populations from different geographic regions, with Indigenous peoples of America showing both higher frequencies of strongest and lower frequencies of weakest binders. As many HLA proteins are strong binders of peptides from distinct viral families, we discuss this result in relation to possible signatures of natural selection on HLA promiscuous alleles due to undetermined past pathogenic infections. Although highly relevant for evolutionary genetics and the development of vaccine therapies, these results should not lead to forget that individual resistance and vulnerability to diseases go beyond the sole HLA allelic affinity and depend on multiple, complex and often unknown biological, environmental and other variables. This article is protected by copyright. All rights reserved.

Cavalli-Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east-west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially-explicit computer simulations of genetic data under the following scenarios: (i) variable levels of admixture between Paleolithic and Neolithic populations, (ii) migration through long-distance dispersal (LDD), (iii) Paleolithic range contraction induced by the last glacial maximum (LGM) and, (iv) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analyzed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east-west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east-west geographic axis of this continent.

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (A) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (B) De novo bioindicator analyses; (C) Structural community metrics including inferred ecological networks; and (D) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programs that leverage recent analytical advancements, while pointing out current limitations and future research needs.