Integrated spatial and multimodal single-cell transcriptomics reveal patient-dependent cell heterogeneity in splenic marginal zone lymphoma.

Biological hallmarks of splenic marginal zone lymphoma (SMZL) remain poorly described. Herein, we performed in-depth SMZL characterization through multimodal single-cell analyses of paired blood/spleen samples. The 3'-single-cell RNA-sequencing, Cellular Indexing of Transcriptomes and Epitopes by sequencing, and 5'-V(D)J single-cell RNA-sequencing datasets were integrated to characterize SMZL transcriptome profiles, including B-cell receptor and T-cell receptor repertoires. Hyperexpanded B-cell clones in the spleen were at a memory-like stage, whereas recirculating tumor B-cells in blood encompassed multiple differentiation stages, indicating an unexpected desynchronization of the B-cell maturation program in SMZL cells. Spatial transcriptomics showed the enrichment of T-effector and T-follicular helper (TFH) signatures in the nodular subtype of SMZL. This latter also exhibited gene-based cell-cell interactions suggestive of dynamic crosstalk between TFH and cancer cells in transcriptomics, further substantiated by using imaging mass cytometry. Our findings provide a comprehensive high-resolution description of SMZL biological hallmarks and characterize, for the first time in situ, inter- and intra-patient heterogeneity at both transcriptomic and protein levels. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Cell specification and functional interactions in the pig blastocyst inferred from single-cell transcriptomics and uterine fluids proteomics.

The embryonic development of the pig comprises a long in utero pre- and peri-implantation development, which dramatically differs from mice and humans. During this peri-implantation period, a complex series of paracrine signals establishes an intimate dialogue between the embryo and the uterus. To better understand the biology of the pig blastocyst during this period, we generated a large dataset of single-cell RNAseq from early and hatched blastocysts, spheroid and ovoid conceptus and proteomic datasets from corresponding uterine fluids. Our results confirm the molecular specificity and functionality of the three main cell populations. We also discovered two previously unknown subpopulations of the trophectoderm, one characterised by the expression of LRP2, which could represent progenitor cells, and the other, expressing pro-apoptotic markers, which could correspond to the Rauber's layer. Our work provides new insights into the biology of these populations, their reciprocal functional interactions, and the molecular dialogue with the maternal uterine environment.

Single-cell RNA sequencing identifies senescence as therapeutic target in rhabdomyolysis-induced acute kidney injury.

BACKGROUND: The role of macrophages in the development of rhabdomyolysis-induced acute kidney injury (RM-AKI) has been established, but an in-depth understanding of the changes in the immune landscape could help to improve targeted strategies. Whereas senescence is usually associated with chronic kidney processes, we also wished to explore whether senescence could also occur in AKI and whether senolytics could act on immune cells. METHODS: Single-cell RNA sequencing was used in the murine glycerol-induced RM-AKI model to dissect the transcriptomic characteristics of CD45+ live cells sorted from kidneys 2 days after injury. Public datasets from murine AKI models were reanalysed to explore cellular senescence signature in tubular epithelial cells (TECs). A combination of senolytics (dasatinib and quercetin, DQ) was administered to mice exposed or not to RM-AKI. RESULTS: Unsupervised clustering of nearly 17 000 single-cell transcriptomes identified seven known immune cell clusters. Sub-clustering of the mononuclear phagocyte cells revealed nine distinct cell sub-populations differently modified with RM. One macrophage cluster was particularly interesting since it behaved as a critical node in a trajectory connecting one major histocompatibility complex class IIhigh (MHCIIhigh) cluster only present in Control to two MHCIIlow clusters only present in RM-AKI. This critical cluster expressed a senescence gene signature, that was very different from that of the TECs. Senolytic DQ treatment blocked the switch from a F4/80highCD11blow to F4/80lowCD11bhigh phenotype, which correlated with prolonged nephroprotection in RM-AKI. CONCLUSIONS: Single-cell RNA sequencing unmasked novel transitional macrophage subpopulation associated with RM-AKI characterized by the activation of cellular senescence processes. This work provides a proof-of-concept that senolytics nephroprotective effects may rely, at least in part, on subtle immune modulation.

Nuclear HMGB1 protects from nonalcoholic fatty liver disease through negative regulation of liver X receptor.

Dysregulations of lipid metabolism in the liver may trigger steatosis progression, leading to potentially severe clinical consequences such as nonalcoholic fatty liver diseases (NAFLDs). Molecular mechanisms underlying liver lipogenesis are very complex and fine-tuned by chromatin dynamics and multiple key transcription factors. Here, we demonstrate that the nuclear factor HMGB1 acts as a strong repressor of liver lipogenesis. Mice with liver-specific Hmgb1 deficiency display exacerbated liver steatosis, while Hmgb1-overexpressing mice exhibited a protection from fatty liver progression when subjected to nutritional stress. Global transcriptome and functional analysis revealed that the deletion of Hmgb1 gene enhances LXRα and PPARγ activity. HMGB1 repression is not mediated through nucleosome landscape reorganization but rather via a preferential DNA occupation in a region carrying genes regulated by LXRα and PPARγ. Together, these findings suggest that hepatocellular HMGB1 protects from liver steatosis development. HMGB1 may constitute a new attractive option to therapeutically target the LXRα-PPARγ axis during NAFLD.

Tumor-associated high endothelial venules mediate lymphocyte entry into tumors and predict response to PD-1 plus CTLA-4 combination immunotherapy.

Recruitment of lymphocytes into tumors is critical for anti-tumor immunity and efficacious immunotherapy. We show in murine models that tumor-associated high endothelial venules (TA-HEVs) are major sites of lymphocyte entry into tumors at baseline and upon treatment with anti-PD-1/anti-CTLA-4 immune checkpoint blockade (ICB). TA-HEV endothelial cells (TA-HECs) derive from post-capillary venules, co-express MECA-79+ HEV sialomucins and E/P-selectins, and are associated with homing and infiltration into tumors of various T cell subsets. Intravital microscopy further shows that TA-HEVs are the main sites of lymphocyte arrest and extravasation into ICB-treated tumors. Increasing TA-HEC frequency and maturation increases the proportion of tumor-infiltrating stem-like CD8+ T cells, and ameliorates ICB efficacy. Analysis of tumor biopsies from 93 patients with metastatic melanoma reveals that TA-HEVs are predictive of better response and survival upon treatment with anti-PD-1/anti-CTLA-4 combination. These studies provide critical insights into the mechanisms governing lymphocyte trafficking in cancer immunity and immunotherapy.

Identification of Loci Enabling Stable and High-Level Heterologous Gene Expression.

Efficient and reliable genome engineering technologies have yet to be developed for diatoms. The delivery of DNA in diatoms results in the random integration of multiple copies, quite often leading to heterogeneous gene activity, as well as host instability. Transgenic diatoms are generally selected on the basis of transgene expression or high enzyme activity, without consideration of the copy number or the integration locus. Here, we propose an integrated pipeline for the diatom, Phaeodactylum tricornutum, that accurately quantifies transgene activity using a ß-glucuronidase assay and the number of transgene copies integrated into the genome through Droplet Digital PCR (ddPCR). An exhaustive and systematic analysis performed on 93 strains indicated that 42% of them exhibited high ß-glucuronidase activity. Though most were attributed to high transgene copy numbers, we succeeded in isolating single-copy clones, as well as sequencing the integration loci. In addition to demonstrating the impact of the genomic integration site on gene activity, this study identifies integration sites for stable transgene expression in Phaeodactylum tricornutum.

Single-cell RNA sequencing unveils the shared and the distinct cytotoxic hallmarks of human TCRVδ1 and TCRVδ2 γδ T lymphocytes.

γδ T lymphocytes represent ∼1% of human peripheral blood mononuclear cells and even more cells in most tissues of vertebrates. Although they have important anticancer functions, most current single-cell RNA sequencing (scRNA-seq) studies do not identify γδ T lymphocytes because their transcriptomes at the single-cell level are unknown. Here we show that high-resolution clustering of large scRNA-seq datasets and a combination of gene signatures allow the specific detection of human γδ T lymphocytes and identification of their T cell receptor (TCR)Vδ1 and TCRVδ2 subsets in large datasets from complex cell mixtures. In t-distributed stochastic neighbor embedding plots from blood and tumor samples, the few γδ T lymphocytes appear collectively embedded between cytotoxic CD8 T and NK cells. Their TCRVδ1 and TCRVδ2 subsets form close yet distinct subclusters, respectively neighboring NK and CD8 T cells because of expression of shared and distinct cytotoxic maturation genes. Similar pseudotime maturation trajectories of TCRVδ1 and TCRVδ2 γδ T lymphocytes were discovered, unveiling in both subsets an unattended pool of terminally differentiated effector memory cells with preserved proliferative capacity, a finding confirmed by in vitro proliferation assays. Overall, the single-cell transcriptomes of thousands of individual γδ T lymphocytes from different CMV+ and CMV- donors reflect cytotoxic maturation stages driven by the immunological history of donors. This landmark study establishes the rationale for identification, subtyping, and deep characterization of human γδ T lymphocytes in further scRNA-seq studies of complex tissues in physiological and disease conditions.

Single-Cell Analysis Reveals Heterogeneity of High Endothelial Venules and Different Regulation of Genes Controlling Lymphocyte Entry to Lymph Nodes.

High-endothelial venules (HEVs) are specialized blood vessels allowing recirculation of naive lymphocytes through lymphoid organs. Here, using full-length, single-cell RNA sequencing, RNA fluorescence in situ hybridization (FISH), flow cytometry, and immunohistofluorescence, we reveal the heterogeneity of HEVs in adult mouse peripheral lymph nodes (PLNs) under conditions of homeostasis, antigenic stimulation, and after inhibition of lymphotoxin-ß receptor (LTßR) signaling. We demonstrate that HEV endothelial cells are in an activated state during homeostasis, and we identify the genes characteristic of the differentiated HEV phenotype. We show that LTßR signaling regulates many HEV genes and pathways in resting PLNs and that immune stimulation induces a global and temporary inflammatory phenotype in HEVs without compromising their ability to recruit naive lymphocytes. Most importantly, we uncover differences in the regulation of genes controlling lymphocyte trafficking, Glycam1, Fut7, Gcnt1, Chst4, B3gnt3, and Ccl21a, that have implications for HEV function and regulation in health and disease.

Niacin induces miR-502-3p expression which impairs insulin sensitivity in human adipocytes.

MicroRNAs have been involved in insulin resistance (IR). As the mechanism whereby niacin, an anti-dyslipidemic agent, leads to IR remains elusive, we sought to identify differentially expressed microRNAs in adipose tissue (AT) of individuals receiving niacin and to explore the link between microRNAs, niacin and IR in human adipocytes.In a double-blind controlled study, 22 obese men received extended-release niacin or placebo over 8 weeks. Bioclinical data and subcutaneous AT biopsies were obtained before and after treatment. AT microRNA expression profiles were determined using RTqPCR for 758 human-specific microRNAs. hMADS adipocytes were treated with niacin, or acipimox (a niacin-like drug without effect on IR), or transfected with miR-502-3p. Glucose uptake and Western blotting were performed.In obese men, insulin sensitivity decreased after niacin treatment. In AT, the expression of 6 microRNAs including miR-502-3p was up-regulated. Treatment of hMADS adipocytes with niacin specifically increased miR-502-3p expression. Acipimox had no effect. Overexpression of miR-502-3p in adipocytes led to reduced insulin-induced glucose uptake and lower insulin-stimulated AKT phosphorylation.Long term niacin treatment altered microRNA expression levels in human AT. Increased miR-502-3p expression may play a role in the mediation of IR due to niacin in adipocytes.The study is registered in Clinical Trials NCT01083329 and EudraCT 2009-012124-85.

Caloric Restriction and Diet-Induced Weight Loss Do Not Induce Browning of Human Subcutaneous White Adipose Tissue in Women and Men with Obesity.

Caloric restriction (CR) is standard lifestyle therapy in obesity management. CR-induced weight loss improves the metabolic profile of individuals with obesity. In mice, occurrence of beige fat cells in white fat depots favors a metabolically healthy phenotype, and CR promotes browning of white adipose tissue (WAT). Here, human subcutaneous abdominal WAT samples were analyzed in 289 individuals with obesity following a two-phase dietary intervention consisting of an 8 week very low calorie diet and a 6-month weight-maintenance phase. Before the intervention, we show sex differences and seasonal variation, with higher expression of brown and beige markers in women with obesity and during winter, respectively. The very low calorie diet resulted in decreased browning of subcutaneous abdominal WAT. During the whole dietary intervention, evolution of body fat and insulin resistance was independent of changes in brown and beige fat markers. These data suggest that diet-induced effects on body fat and insulin resistance are independent of subcutaneous abdominal WAT browning in people with obesity.

Muscle transcriptome analysis reveals molecular pathways and biomarkers involved in extreme ultimate pH and meat defect occurrence in chicken.

The processing ability and sensory quality of chicken breast meat are highly related to its ultimate pH (pHu), which is mainly determined by the amount of glycogen in the muscle at death. To unravel the molecular mechanisms underlying glycogen and meat pHu variations and to identify predictive biomarkers of these traits, a transcriptome profiling analysis was performed using an Agilent custom chicken 8 × 60 K microarray. The breast muscle gene expression patterns were studied in two chicken lines experimentally selected for high (pHu+) and low (pHu-) pHu values of the breast meat. Across the 1,436 differentially expressed (DE) genes found between the two lines, many were involved in biological processes related to muscle development and remodelling and carbohydrate and energy metabolism. The functional analysis showed an intensive use of carbohydrate metabolism to produce energy in the pHu- line, while alternative catabolic pathways were solicited in the muscle of the pHu+ broilers, compromising their muscle development and integrity. After a validation step on a population of 278 broilers using microfluidic RT-qPCR, 20 genes were identified by partial least squares regression as good predictors of the pHu, opening new perspectives of screening broilers likely to present meat quality defects.

Circulating miR-155, miR-145 and let-7c as diagnostic biomarkers of the coronary artery disease.

Coronary artery disease (CAD) is the most prevalent cause of mortality and morbidity worldwide and the number of individuals at risk is increasing. To better manage cardiovascular diseases, improved tools for risk prediction including the identification of novel accurate biomarkers are needed. MicroRNA (miRNA) are essential post-transcriptional modulators of gene expression leading to mRNA suppression or translational repression. Specific expression profiles of circulating miRNA have emerged as potential noninvasive diagnostic biomarkers of diseases. The aim of this study was to identify the potential diagnostic value of circulating miRNA with CAD. Circulating miR-145, miR-155, miR-92a and let-7c were selected and validated by quantitative PCR in 69 patients with CAD and 30 control subjects from the cross-sectional study GENES. The expression of miR-145, miR-155 and let-7c showed significantly reduced expression in patients with CAD compared to controls. Multivariate logistic regression analysis revealed that low levels of circulating let-7c, miR-145 and miR-155 were associated with CAD. Receiver operating curves analysis showed that let-7c, miR-145 or miR-155 were powerful markers for detecting CAD. Furthermore, we demonstrated that the combination of the three circulating miRNA managed to deliver a specific signature for diagnosing CAD.

Genetic diversity, linkage disequilibrium and power of a large grapevine (Vitis vinifera L) diversity panel newly designed for association studies.

BACKGROUND: As for many crops, new high-quality grapevine varieties requiring less pesticide and adapted to climate change are needed. In perennial species, breeding is a long process which can be speeded up by gaining knowledge about quantitative trait loci linked to agronomic traits variation. However, due to the long juvenile period of these species, establishing numerous highly recombinant populations for high resolution mapping is both costly and time-consuming. Genome wide association studies in germplasm panels is an alternative method of choice, since it allows identifying the main quantitative trait loci with high resolution by exploiting past recombination events between cultivars. Such studies require adequate panel design to represent most of the available genetic and phenotypic diversity. Assessing linkage disequilibrium extent and panel power is also needed to determine the marker density required for association studies. RESULTS: Starting from the largest grapevine collection worldwide maintained in Vassal (France), we designed a diversity panel of 279 cultivars with limited relatedness, reflecting the low structuration in three genetic pools resulting from different uses (table vs wine) and geographical origin (East vs West), and including the major founders of modern cultivars. With 20 simple sequence repeat markers and five quantitative traits, we showed that our panel adequately captured most of the genetic and phenotypic diversity existing within the entire Vassal collection. To assess linkage disequilibrium extent and panel power, we genotyped single nucleotide polymorphisms: 372 over four genomic regions and 129 distributed over the whole genome. Linkage disequilibrium, measured by correlation corrected for kinship, reached 0.2 for a physical distance between 9 and 458 Kb depending on genetic pool and genomic region, with varying size of linkage disequilibrium blocks. This panel achieved reasonable power to detect associations between traits with high broad-sense heritability (> 0.7) and causal loci with intermediate allelic frequency and strong effect (explaining > 10 % of total variance). CONCLUSIONS: Our association panel constitutes a new, highly valuable resource for genetic association studies in grapevine, and deserves dissemination to diverse field and greenhouse trials to gain more insight into the genetic control of many agronomic traits and their interaction with the environment.

Salivary MicroRNA in Pancreatic Cancer Patients.

BACKGROUND: Pancreatic cancer is the fourth leading cause of cancer death in Western countries, with the lowest 1-year survival rate among commonly diagnosed cancers. Reliable biomarkers for pancreatic cancer diagnosis are lacking and are urgently needed to allow for curative surgery. As microRNA (miRNA) recently emerged as candidate biomarkers for this disease, we explored in the present pilot study the differences in salivary microRNA profiles between patients with pancreatic tumors that are not eligible for surgery, precancerous lesions, inflammatory disease or cancer-free patients as a potential early diagnostic tool. METHODS: Whole saliva samples from patients with pancreatic cancer (n = 7), pancreatitis (n = 4), IPMN (n = 2), or healthy controls (n = 4) were obtained during endoscopic examination. After total RNA isolation, expression of 94 candidate miRNAs was screened by q(RT)PCR using Biomark Fluidgm. Human-derived pancreatic cancer cells were xenografted in athymic mice as an experimental model of pancreatic cancer. RESULTS: We identified hsa-miR-21, hsa-miR-23a, hsa-miR-23b and miR-29c as being significantly upregulated in saliva of pancreatic cancer patients compared to control, showing sensitivities of 71.4%, 85.7%, 85,7% and 57%, respectively and excellent specificity (100%). Interestingly, hsa-miR-23a and hsa-miR23b are overexpressed in the saliva of patients with pancreatic cancer precursor lesions. We found that hsa-miR-210 and let-7c are overexpressed in the saliva of patients with pancreatitis as compared to the control group, with sensitivity of 100% and 75%, and specificity of 100% and 80%, respectively. Last hsa-miR-216 was upregulated in cancer patients as compared to patients diagnosed with pancreatitis, with sensitivity of 50% and specificity of 100%. In experimental models of PDAC, salivary microRNA detection precedes systemic detection of cancer cells markers. CONCLUSIONS: Our novel findings indicate that salivary miRNA are discriminatory in pancreatic cancer patients that are not eligible for surgery. In addition, we demonstrate in experimental models that salivary miRNA detection precedes systemic detection of cancer cells markers. This study stems for the use of salivary miRNA as biomarker for the early diagnosis of patients with unresectable pancreatic cancer.

First-in-man phase 1 clinical trial of gene therapy for advanced pancreatic cancer: safety, biodistribution, and preliminary clinical findings.

This phase 1 trial was aimed to determine the safety, pharmacokinetics, and preliminary clinical activity of CYL-02, a nonviral gene therapy product that sensitizes pancreatic cancer cells to chemotherapy. CYL-02 was administrated using endoscopic ultrasound in 22 patients with pancreatic cancer that concomitantly received chemotherapy (gemcitabine). The maximum-tolerated dose (MTD) exceeded the maximal feasible dose of CYL-02 and was not identified. Treatment-related toxicities were mild, without serious adverse events. Pharmacokinetic analysis revealed a dose-dependent increase in CYL-02 DNA exposure in blood and tumors, while therapeutic RNAs were detected in tumors. No objective response was observed, but nine patients showed stable disease up to 6 months following treatment and two of these patients experienced long-term survival. Panels of plasmatic microRNAs and proteins were identified as predictive of gene therapy efficacy. We demonstrate that CYL-02 nonviral gene therapy has a favorable safety profile and is well tolerated in patients. We characterize CYL-02 biodistribution and demonstrate therapeutic gene expression in tumors. Treated patients experienced stability of disease and predictive biomarkers of response to treatment were identified. These promising results warrant further evaluation in phase 2 clinical trial.

Specific macrophage subtypes influence the progression of rhabdomyolysis-induced kidney injury.

Rhabdomyolysis can be life threatening if complicated by AKI. Macrophage infiltration has been observed in rat kidneys after glycerol-induced rhabdomyolysis, but the role of macrophages in rhabdomyolysis-induced AKI remains unknown. Here, in a patient diagnosed with rhabdomyolysis, we detected substantial macrophage infiltration in the kidney. In a mouse model of rhabdomyolysis-induced AKI, diverse renal macrophage phenotypes were observed depending on the stage of the disease. Two days after rhabdomyolysis, F4/80(low)CD11b(high)Ly6b(high)CD206(low) kidney macrophages were dominant, whereas by day 8, F4/80(high)CD11b(+)Ly6b(low)CD206(high) cells became the most abundant. Single-cell gene expression analyses of FACS-sorted macrophages revealed that these subpopulations were heterogeneous and that individual cells simultaneously expressed both M1 and M2 markers. Liposomal clodronate-mediated macrophage depletion significantly reduced the early infiltration of F4/80(low)CD11b(high)Ly6b(high)CD206(low) macrophages. Furthermore, transcriptionally regulated targets potentially involved in disease progression, including fibronectin, collagen III, and chemoattractants that were identified via single-cell analysis, were verified as macrophage-dependent in situ. In vitro, myoglobin treatment induced proximal tubular cells to secrete chemoattractants and macrophages to express proinflammatory markers. At day 30, liposomal clodronate-mediated macrophage depletion reduced fibrosis and improved both kidney repair and mouse survival. Seven months after rhabdomyolysis, histologic lesions were still present but were substantially reduced with prior depletion of macrophages. These results suggest an important role for macrophages in rhabdomyolysis-induced AKI progression and advocate the utility of long-term follow-up for patients with this disease.

Gene-based single nucleotide polymorphism discovery in bovine muscle using next-generation transcriptomic sequencing.

BACKGROUND: Genetic information based on molecular markers has increasingly being used in cattle breeding improvement programmes, as a mean to improve conventionally phenotypic selection. Advances in molecular genetics have led to the identification of several genetic markers associated with genes affecting economic traits. Until recently, the identification of the causative genetic variants involved in the phenotypes of interest has remained a difficult task. The advent of novel sequencing technologies now offers a new opportunity for the identification of such variants. Despite sequencing costs plummeting, sequencing whole-genomes or large targeted regions is still too expensive for most laboratories. A transcriptomic-based sequencing approach offers a cheaper alternative to identify a large number of polymorphisms and possibly to discover causative variants. In the present study, we performed a gene-based single nucleotide polymorphism (SNP) discovery analysis in bovine Longissimus thoraci, using RNA-Seq. To our knowledge, this represents the first study done in bovine muscle. RESULTS: Messenger RNAs from Longissimus thoraci from three Limousin bull calves were subjected to high-throughput sequencing. Approximately 36-46 million paired-end reads were obtained per library. A total of 19,752 transcripts were identified and 34,376 different SNPs were detected. Fifty-five percent of the SNPs were found in coding regions and ~22% resulted in an amino acid change. Applying a very stringent SNP quality threshold, we detected 8,407 different high-confidence SNPs, 18% of which are non synonymous coding SNPs. To analyse the accuracy of RNA-Seq technology for SNP detection, 48 SNPs were selected for validation by genotyping. No discrepancies were observed when using the highest SNP probability threshold. To test the usefulness of the identified SNPs, the 48 selected SNPs were assessed by genotyping 93 bovine samples, representing mostly the nine major breeds used in France. Principal component analysis indicates a clear separation between the nine populations. CONCLUSIONS: The RNA-Seq data and the collection of newly discovered coding SNPs improve the genomic resources available for cattle, especially for beef breeds. The large amount of variation present in genes expressed in Limousin Longissimus thoracis, especially the large number of non synonymous coding SNPs, may prove useful to study the mechanisms underlying the genetic variability of meat quality traits.

Single-nucleotide polymorphism discovery and diversity in the model legume Medicago truncatula.

Extensive genomic resources are available in the model legume Medicago truncatula. Here, we present the discovery and design of the first array of single-nucleotide polymorphism (SNP) markers in M. truncatula through large-scale Sanger resequencing of genomic fragments spanning the genome, in a diverse panel of 16 M. truncatula accessions. Both anonymous fragments and fragments targeting candidate genes for flowering phenology and symbiosis were surveyed for nucleotide variation in almost 230 kb of unique genomic regions. A set of 384 SNP markers was designed for an Illumina's GoldenGate assay, genotyped on a collection of 192 inbred lines (CC192) representing the geographical range of the species and used to survey the diversity of two natural populations. Finally, 86% of the tested SNPs were of high quality and exhibited polymorphism in the CC192 collection. Even at the population level, we detected polymorphism for more than 50% of the selected SNPs. Analysis of the allele frequency spectrum in the CC192 showed a reduced ascertainment bias, mostly limited to very rare alleles (frequency <0.01). The substantial polymorphism detected at the species and population levels, the high marker quality and the potential to survey large samples of individuals make this set of SNP markers a valuable tool to improve our understanding of the effect of demographic and selective factors that shape the natural genetic diversity within the selfing species Medicago truncatula.

DNA fingerprinting and new tools for fine-scale discrimination of Arabidopsis thaliana accessions.

One of the main strengths of Arabidopsis thaliana as a model species is the impressive number of public resources available to the scientific community. Exploring species genetic diversity--and therefore adaptation--relies on collections of individuals from natural populations taken from diverse environments. Nevertheless, due to a few mislabeling events or genotype mixtures, some variants available in stock centers have been misidentified, causing inconsistencies and limiting the potential of genetic analyses. To improve the identification of natural accessions, we genotyped 1311 seed stocks from our Versailles Arabidopsis Stock Center and from other collections to determine their molecular profiles at 341 single nucleotide polymorphism markers. These profiles were used to compare genotypes at both the intra- and inter-accession levels. We confirmed previously described inconsistencies and revealed new ones, and suggest likely identities for accessions whose lineage had been lost. We also developed two new tools: a minimal fingerprint computation to quickly verify the identity of an accession, and an optimized marker set to assist in the identification of unknown or mixed accessions. These tools are available on a dedicated web interface called ANATool (https://www.versailles.inra.fr/ijpb/crb/anatool) that provides a simple and efficient means to verify or determine the identity of A. thaliana accessions in any laboratory, without the need for any specific or expensive technology.