Mapping interindividual dynamics of innate immune response at single-cell resolution.

Common genetic variants across individuals modulate the cellular response to pathogens and are implicated in diverse immune pathologies, yet how they dynamically alter the response upon infection is not well understood. Here, we triggered antiviral responses in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-sequencing. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), a statistical approach designed to identify nonlinear dynamic genetic effects across transcriptional trajectories of cells. This approach identified 1,275 expression quantitative trait loci (local false discovery rate 10%) that manifested during the responses, many of which were colocalized with susceptibility loci identified by genome-wide association studies of infectious and autoimmune diseases, including the OAS1 splicing quantitative trait locus in a COVID-19 susceptibility locus. In summary, our analytical approach provides a unique framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

A systematic CRISPR screen defines mutational mechanisms underpinning signatures caused by replication errors and endogenous DNA damage.

Mutational signatures are imprints of pathophysiological processes arising through tumorigenesis. We generated isogenic CRISPR-Cas9 knockouts (Δ) of 43 genes in human induced pluripotent stem cells, cultured them in the absence of added DNA damage, and performed whole-genome sequencing of 173 subclones. ΔOGG1, ΔUNG, ΔEXO1, ΔRNF168, ΔMLH1, ΔMSH2, ΔMSH6, ΔPMS1, and ΔPMS2 produced marked mutational signatures indicative of being critical mitigators of endogenous DNA modifications. Detailed analyses revealed mutational mechanistic insights, including how 8-oxo-dG elimination is sequence-context-specific while uracil clearance is sequence-context-independent. Mismatch repair (MMR) deficiency signatures are engendered by oxidative damage (C>A transversions), differential misincorporation by replicative polymerases (T>C and C>T transitions), and we propose a 'reverse template slippage' model for T>A transversions. ΔMLH1, ΔMSH6, and ΔMSH2 signatures were similar to each other but distinct from ΔPMS2. Finally, we developed a classifier, MMRDetect, where application to 7,695 WGS cancers showed enhanced detection of MMR-deficient tumors, with implications for responsiveness to immunotherapies.

A Compendium of Mutational Signatures of Environmental Agents.

Whole-genome-sequencing (WGS) of human tumors has revealed distinct mutation patterns that hint at the causative origins of cancer. We examined mutational signatures in 324 WGS human-induced pluripotent stem cells exposed to 79 known or suspected environmental carcinogens. Forty-one yielded characteristic substitution mutational signatures. Some were similar to signatures found in human tumors. Additionally, six agents produced double-substitution signatures and eight produced indel signatures. Investigating mutation asymmetries across genome topography revealed fully functional mismatch and transcription-coupled repair pathways. DNA damage induced by environmental mutagens can be resolved by disparate repair and/or replicative pathways, resulting in an assortment of signature outcomes even for a single agent. This compendium of experimentally induced mutational signatures permits further exploration of roles of environmental agents in cancer etiology and underscores how human stem cell DNA is directly vulnerable to environmental agents. VIDEO ABSTRACT.

Inhibition of ß-catenin signalling in dermal fibroblasts enhances hair follicle regeneration during wound healing.

New hair follicles (HFs) do not form in adult mammalian skin unless epidermal Wnt signalling is activated genetically or within large wounds. To understand the postnatal loss of hair forming ability we monitored HF formation at small circular (2 mm) wound sites. At P2, new HFs formed in back skin, but HF formation was markedly decreased by P21. Neonatal tail also formed wound-associated HFs, albeit in smaller numbers. Postnatal loss of HF neogenesis did not correlate with wound closure rate but with a reduction in Lrig1-positive papillary fibroblasts in wounds. Comparative gene expression profiling of back and tail dermis at P1 and dorsal fibroblasts at P2 and P50 showed a correlation between loss of HF formation and decreased expression of genes associated with proliferation and Wnt/ß-catenin activity. Between P2 and P50, fibroblast density declined throughout the dermis and clones of fibroblasts became more dispersed. This correlated with a decline in fibroblasts expressing a TOPGFP reporter of Wnt activation. Surprisingly, between P2 and P50 there was no difference in fibroblast proliferation at the wound site but Wnt signalling was highly upregulated in healing dermis of P21 compared with P2 mice. Postnatal ß-catenin ablation in fibroblasts promoted HF regeneration in neonatal and adult mouse wounds, whereas ß-catenin activation reduced HF regeneration in neonatal wounds. Our data support a model whereby postnatal loss of hair forming ability in wounds reflects elevated dermal Wnt/ß-catenin activation in the wound bed, increasing the abundance of fibroblasts that are unable to induce HF formation.

Shift in precipitation regime promotes interspecific hybridization of introduced Coffea species.

The frequency of plant species introductions has increased in a highly connected world, modifying species distribution patterns to include areas outside their natural ranges. These introductions provide the opportunity to gain new insight into the importance of flowering phenology as a component of adaptation to a new environment. Three Coffea species, C. arabica, C. canephora (Robusta), and C. liberica, native to intertropical Africa have been introduced to New Caledonia. On this archipelago, a secondary contact zone has been characterized where these species coexist, persist, and hybridize spontaneously. We investigated the impact of environmental changes undergone by each species following its introduction in New Caledonia on flowering phenology and overcoming reproductive barriers between sister species. We developed species distribution models and compared both environmental envelopes and climatic niches between native and introduced hybrid zones. Flowering phenology was monitored in a population in the hybrid zone along with temperature and precipitation sequences recorded at a nearby weather station. The extent and nature of hybridization events were characterized using chloroplast and nuclear microsatellite markers. The three Coffea species encountered weak environmental suitability compared to their native ranges when introduced to New Caledonia, especially C. arabica and C. canephora. The niche of the New Caledonia hybrid zone was significantly different from all three species' native niches based on identity tests (I Similarity and D Schoener's Similarity Indexes). This area appeared to exhibit intermediate conditions between the native conditions of the three species for temperature-related variables and divergent conditions for precipitation-related ones. Flowering pattern in these Coffea species was shown to have a strong genetic component that determined the time between the triggering rain and anthesis (flower opening), specific to each species. However, a precipitation regime different from those in Africa was directly involved in generating partial flowering overlap between species and thus in allowing hybridization and interspecific gene flow. Interspecific hybrids accounted for 4% of the mature individuals in the sympatric population and occurred between each pair of species with various level of introgression. Adaptation to new environmental conditions following introduction of Coffea species to New Caledonia has resulted in a secondary contact between three related species, which would not have happened in their native ranges, leading to hybridization and gene flow.

Species distribution models of tropical deep-sea snappers.

Deep-sea fisheries provide an important source of protein to Pacific Island countries and territories that are highly dependent on fish for food security. However, spatial management of these deep-sea habitats is hindered by insufficient data. We developed species distribution models using spatially limited presence data for the main harvested species in the Western Central Pacific Ocean. We used bathymetric and water temperature data to develop presence-only species distribution models for the commercially exploited deep-sea snappers Etelis Cuvier 1828, Pristipomoides Valenciennes 1830, and Aphareus Cuvier 1830. We evaluated the performance of four different algorithms (CTA, GLM, MARS, and MAXENT) within the BIOMOD framework to obtain an ensemble of predicted distributions. We projected these predictions across the Western Central Pacific Ocean to produce maps of potential deep-sea snapper distributions in 32 countries and territories. Depth was consistently the best predictor of presence for all species groups across all models. Bathymetric slope was consistently the poorest predictor. Temperature at depth was a good predictor of presence for GLM only. Model precision was highest for MAXENT and CTA. There were strong regional patterns in predicted distribution of suitable habitat, with the largest areas of suitable habitat (> 35% of the Exclusive Economic Zone) predicted in seven South Pacific countries and territories (Fiji, Matthew & Hunter, Nauru, New Caledonia, Tonga, Vanuatu and Wallis & Futuna). Predicted habitat also varied among species, with the proportion of predicted habitat highest for Aphareus and lowest for Etelis. Despite data paucity, the relationship between deep-sea snapper presence and their environments was sufficiently strong to predict their distribution across a large area of the Pacific Ocean. Our results therefore provide a strong baseline for designing monitoring programs that balance resource exploitation and conservation planning, and for predicting future distributions of deep-sea snappers.

Wildfire risk for main vegetation units in a biodiversity hotspot: modeling approach in New Caledonia, South Pacific.

Wildfire has been recognized as one of the most ubiquitous disturbance agents to impact on natural environments. In this study, our main objective was to propose a modeling approach to investigate the potential impact of wildfire on biodiversity. The method is illustrated with an application example in New Caledonia where conservation and sustainable biodiversity management represent an important challenge. Firstly, a biodiversity loss index, including the diversity and the vulnerability indexes, was calculated for every vegetation unit in New Caledonia and mapped according to its distribution over the New Caledonian mainland. Then, based on spatially explicit fire behavior simulations (using the FLAMMAP software) and fire ignition probabilities, two original fire risk assessment approaches were proposed: a one-off event model and a multi-event burn probability model. The spatial distribution of fire risk across New Caledonia was similar for both indices with very small localized spots having high risk. The patterns relating to highest risk are all located around the remaining sclerophyll forest fragments and are representing 0.012% of the mainland surface. A small part of maquis and areas adjacent to dense humid forest on ultramafic substrates should also be monitored. Vegetation interfaces between secondary and primary units displayed high risk and should represent priority zones for fire effects mitigation. Low fire ignition probability in anthropogenic-free areas decreases drastically the risk. A one-off event associated risk allowed localizing of the most likely ignition areas with potential for extensive damage. Emergency actions could aim limiting specific fire spread known to have high impact or consist of on targeting high risk areas to limit one-off fire ignitions. Spatially explicit information on burning probability is necessary for setting strategic fire and fuel management planning. Both risk indices provide clues to preserve New Caledonia hot spot of biodiversity facing wildfires.

The subtle role of climate change on population genetic structure in Canada lynx.

Anthropogenically driven climatic change is expected to reshape global patterns of species distribution and abundance. Given recent links between genetic variation and environmental patterns, climate change may similarly impact genetic population structure, but we lack information on the spatial and mechanistic underpinnings of genetic-climate associations. Here, we show that current genetic variability of Canada lynx (Lynx canadensis) is strongly correlated with a winter climate gradient (i.e. increasing snow depth and winter precipitation from west-to-east) across the Pacific-North American (PNO) to North Atlantic Oscillation (NAO) climatic systems. This relationship was stronger than isolation by distance and not explained by landscape variables or changes in abundance. Thus, these patterns suggest that individuals restricted dispersal across the climate boundary, likely in the absence of changes in habitat quality. We propose habitat imprinting on snow conditions as one possible explanation for this unusual phenomenon. Coupling historical climate data with future projections, we also found increasingly diverging snow conditions between the two climate systems. Based on genetic simulations using projected climate data (2041-2070), we predicted that this divergence could lead to a threefold increase in genetic differentiation, potentially leading to isolated east-west populations of lynx in North America. Our results imply that subtle genetic structure can be governed by current climate and that substantive genetic differentiation and related ecological divergence may arise from changing climate patterns.

Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia.

Amborella trichopoda Baill. (Amborellaceae, Amborellales), the sole living member of the sister group to all other extant angiosperms, is endemic to New Caledonia. We addressed the intraspecific phylogeography of Amborella by investigating whether its present population genetic structure could be related to its current and past habitats. We found moderate range-wide genetic diversity based on nuclear microsatellite data and detected four well-differentiated, geographically distinct genetic groups using Bayesian clustering analyses. We modelled the ecological niche of Amborella based on the current climatic and environmental conditions. The predictive ability of the model was very good throughout the Central East mainland zone, but Amborella was predicted in the northern part of the island where this plant has not been reported. Furthermore, no significant barrier was detected based on habitat suitability that could explain the genetic differentiation across the area. Conversely, we found that the main genetic clusters could be related to the distribution of the suitable habitat at the last glacial maximum (LGM, c. 21,000 years BP), when Amborella experienced a dramatic 96.5% reduction in suitable area. At least two lineages survived in distinct putative refugia located in the Massif des Lèvres and in the vicinity of Mount Aoupinié. Our findings finally confirmed the importance of LGM rainforest refugia in shaping the current intra- and interspecific diversity in New Caledonian plants and revealed the possibility of an as yet unreported refugium. The combination of niche modelling and population genetics thereby offered novel insight into the biogeographical history of an emblematic taxon.

The interfollicular epidermis of adult mouse tail comprises two distinct cell lineages that are differentially regulated by Wnt, Edaradd, and Lrig1.

Current models of how mouse tail interfollicular epidermis (ife) is maintained overlook the coexistence of two distinct terminal differentiation programs: parakeratotic (scale) and orthokeratotic (interscale). lineage tracing and clonal analysis revealed that scale and interscale are maintained by unipotent cells in the underlying basal layer, with scale progenitors dividing more rapidly than interscale progenitors. Although scales are pigmented and precisely aligned with hair follicles, melanocytes and follicles were not necessary for scale differentiation. Epidermal Wnt signaling was required for scale enlargement during development and for postnatal maintenance of scale-interscale boundaries. Loss of Edaradd inhibited ventral scale formation, whereas loss of Lrig1 led to scale enlargement and fusion. In wild-type skin, Lrig1 was not expressed in IFE but was selectively upregulated in dermal fibroblasts underlying the interscale. We conclude that the different IFE differentiation compartments are maintained by distinct stem cell populations and are regulated by epidermal and dermal signals.

Microsatellite markers for Amborella (Amborellaceae), a monotypic genus endemic to New Caledonia.

PREMISE OF THE STUDY: Informative markers are required for assessing the diversity of Amborella trichopoda, the only species of its order, endemic to New Caledonia and considered to be the sister species to all flowering plants. Therefore, expressed sequence tag (EST)-based microsatellite markers were developed. • METHODS AND RESULTS: Fifty-five microsatellite loci were characterized in 14896 putative unigenes, which were generated by assembling A. trichopoda ESTs from the public sequence database. Seventeen markers revealed polymorphism in 80 adult shrubs from three populations. The number of alleles per locus ranged from two to 12, with a total of 132 alleles scored. The mean expected heterozygosity per population ranged from 0.336 to 0.567. • CONCLUSIONS: These markers offer an appropriate amount of variation to investigate genetic diversity structure, gene flow, and other conservation issues.

Current genetic differentiation of Coffea canephora Pierre ex A. Froehn in the Guineo-Congolian African zone: cumulative impact of ancient climatic changes and recent human activities.

BACKGROUND: Among Coffea species, C. canephora has the widest natural distribution area in tropical African forests. It represents a good model for analyzing the geographical distribution of diversity in relation to locations proposed as part of the "refuge theory". In this study, we used both microsatellite (simple sequence repeat, SSR) and restriction fragment length polymorphism (RFLP) markers to investigate the genetic variation pattern of C. canephora in the Guineo-Congolean distribution zone. RESULTS: Both markers were first compared in terms of their informativeness and efficiency in a study of genetic diversity and relationships among wild C. canephora genotypes. As expected, SSR markers were found to have a higher genetic distance detection capacity than RFLP. Nevertheless, similarity matrices showed significant correlations when Mantel's test was carried out (r = 0.66, p < 0.0001). Finally, both markers were equally effective for group discrimination and phylogenetic studies, but SSR markers tended to outperform RFLP markers in discriminating the source of an individual among diversity groups and in putative hybrid detection. Five well defined genetic groups, one in the Upper Guinean forests, the four others in the Lower Guinean forests, were identified, corresponding to geographical patterning in the individuals. CONCLUSION: Our data suggested that the Dahomey Gap, a biogeographical barrier, played a role in wild C. canephora differentiation. Climatic variations during the Pleistocene and/or Holocene probably caused the subgroup differentiation in the Congolese zone through the presence of a mosaic of putative refugia. Recent hybridization between C. canephora diversity groups, both for spontaneous individuals and cultivars, was further characterised according to their geographic dissemination or breeding history as a consequence of human activities.

Developmental control of segment numbers in vertebrates.

Segmentation or metamery in vertebrates is best illustrated by the repetition of the vertebrae and ribs, their associated skeletal muscles and blood vessels, and the spinal nerves and ganglia. The segment number varies tremendously among the different vertebrate species, ranging from as few as six vertebrae in some frogs to as many as several hundred in some snakes and fish. In vertebrates, metameric segments or somites form sequentially during body axis formation. This results in the embryonic axis becoming entirely segmented into metameric units from the level of the otic vesicle almost to the very tip of the tail. The total segment number mostly depends on two parameters: (1) the control of the posterior growth of the body axis during somitogenesis-more same-size segments can be formed in a longer axis and (2) segment size--more smaller--size segments can be formed in a same-size body axis. During evolution, independent variations of these parameters could explain the huge diversity in segment numbers observed among vertebrate species. These variations in segment numbers are accompanied by diversity in the regionalization of the vertebral column. For example, amniotes can exhibit up to five different types of vertebrae: cervical, thoracic, lumbar, sacral and caudal, the number of which varies according to the species. This regionalization of the vertebral column is controlled by the Hox family of transcription factors. We propose that during development, dissociation of the Hox- and segmentation-clock-dependent vertebral patterning systems explains the enormous diversity of vertebral formulae observed in vertebrates.

Control of segment number in vertebrate embryos.

The vertebrate body axis is subdivided into repeated segments, best exemplified by the vertebrae that derive from embryonic somites. The number of somites is precisely defined for any given species but varies widely from one species to another. To determine the mechanism controlling somite number, we have compared somitogenesis in zebrafish, chicken, mouse and corn snake embryos. Here we present evidence that in all of these species a similar 'clock-and-wavefront' mechanism operates to control somitogenesis; in all of them, somitogenesis is brought to an end through a process in which the presomitic mesoderm, having first increased in size, gradually shrinks until it is exhausted, terminating somite formation. In snake embryos, however, the segmentation clock rate is much faster relative to developmental rate than in other amniotes, leading to a greatly increased number of smaller-sized somites.

FGF8, Wnt8 and Myf5 are target genes of Tbx6 during anteroposterior specification in Xenopus embryo.

The T-box transcription factor Tbx6 is required for somite formation and loss-of-function or reduced activity of Tbx6 result in absence of posterior paraxial mesoderm or disorganized somites, but how it is involved in a regulatory hierarchy during Xenopus early development is not clear. We show here that Tbx6 is expressed in the lateral and ventral mesoderm of early gastrula, and it is necessary and sufficient to directly and indirectly regulate the expression of a subset of early mesodermal and endodermal genes. Ectopic expression of Tbx6 inhibits early neuroectodermal gene expression by strongly inducing the expression of posterior mesodermal genes, and expands the mesoderm territory at the expense of neuroectoderm. Conversely, overexpression of a dominant negative Tbx6 mutant in the ventral mesoderm inhibits the expression of several mesodermal genes and results in neural induction in a dose-dependent manner. Using a hormone-inducible form of Tbx6, we have identified FGF8, Xwnt8 and XMyf5 as immediate early responsive genes of Tbx6, and the induction of these genes by Tbx6 is independent of Xbra and VegT. These target genes act downstream and mediate the function of Tbx6 in anteroposterior specification. Our results therefore identify a regulatory cascade governed by Tbx6 in the specification of posterior mesoderm during Xenopus early development.