Variant calling and genotyping accuracy of ddRAD-seq: Comparison with 20X WGS in layers.

Whole Genome Sequencing (WGS) remains a costly or unsuitable method for routine genotyping of laying hens. Until now, breeding companies have been using or developing SNP chips. Nevertheless, alternatives methods based on sequencing have been developed. Among these, reduced representation sequencing approaches can offer sequencing quality and cost-effectiveness by reducing the genomic regions covered by sequencing. The aim of this study was to evaluate the ability of double digested Restriction site Associated DNA sequencing (ddRAD-seq) to identify and genotype SNPs in laying hens, by comparison with a presumed reliable WGS approach. Firstly, the sensitivity and precision of variant calling and the genotyping reliability of ddRADseq were determined. Next, the SNP Call Rate (CRSNP) and mean depth of sequencing per SNP (DPSNP) were compared between both methods. Finally, the effect of multiple combinations of thresholds for these parameters on genotyping reliability and amount of remaining SNPs in ddRAD-seq was studied. In raw form, the ddRAD-seq identified 349,497 SNPs evenly distributed on the genome with a CRSNP of 0.55, a DPSNP of 11X and a mean genotyping reliability rate per SNP of 80%. Considering genomic regions covered by expected enzymatic fragments (EFs), the sensitivity of the ddRAD-seq was estimated at 32.4% and its precision at 96.4%. The low CRSNP and DPSNP values were explained by the detection of SNPs outside the EFs theoretically generated by the ddRAD-seq protocol. Indeed, SNPs outside the EFs had significantly lower CRSNP (0.25) and DPSNP (1X) values than SNPs within the EFs (0.7 and 17X, resp.). The study demonstrated the relationship between CRSNP, DPSNP, genotyping reliability and the number of SNPs retained, to provide a decision-support tool for defining filtration thresholds. Severe quality control over ddRAD-seq data allowed to retain a minimum of 40% of the SNPs with a CcR of 98%. Then, ddRAD-seq was defined as a suitable method for variant calling and genotyping in layers.

Enriched atlas of lncRNA and protein-coding genes for the GRCg7b chicken assembly and its functional annotation across 47 tissues.

Gene atlases for livestock are steadily improving thanks to new genome assemblies and new expression data improving the gene annotation. However, gene content varies across databases due to differences in RNA sequencing data and bioinformatics pipelines, especially for long non-coding RNAs (lncRNAs) which have higher tissue and developmental specificity and are harder to consistently identify compared to protein coding genes (PCGs). As done previously in 2020 for chicken assemblies galgal5 and GRCg6a, we provide a new gene atlas, lncRNA-enriched, for the latest GRCg7b chicken assembly, integrating "NCBI RefSeq", "EMBL-EBI Ensembl/GENCODE" reference annotations and other resources such as FAANG and NONCODE. As a result, the number of PCGs increases from 18,022 (RefSeq) and 17,007 (Ensembl) to 24,102, and that of lncRNAs from 5789 (RefSeq) and 11,944 (Ensembl) to 44,428. Using 1400 public RNA-seq transcriptome representing 47 tissues, we provided expression evidence for 35,257 (79%) lncRNAs and 22,468 (93%) PCGs, supporting the relevance of this atlas. Further characterization including tissue-specificity, sex-differential expression and gene configurations are provided. We also identified conserved miRNA-hosting genes with human counterparts, suggesting common function. The annotated atlas is available at gega.sigenae.org.

Restriction site-associated DNA sequencing technologies as an alternative to low-density SNP chips for genomic selection: a simulation study in layer chickens.

BACKGROUND: To reduce the cost of genomic selection, a low-density (LD) single nucleotide polymorphism (SNP) chip can be used in combination with imputation for genotyping selection candidates instead of using a high-density (HD) SNP chip. Next-generation sequencing (NGS) techniques have been increasingly used in livestock species but remain expensive for routine use for genomic selection. An alternative and cost-efficient solution is to use restriction site-associated DNA sequencing (RADseq) techniques to sequence only a fraction of the genome using restriction enzymes. From this perspective, use of RADseq techniques followed by an imputation step on HD chip as alternatives to LD chips for genomic selection was studied in a pure layer line. RESULTS: Genome reduction and sequencing fragments were identified on reference genome using four restriction enzymes (EcoRI, TaqI, AvaII and PstI) and a double-digest RADseq (ddRADseq) method (TaqI-PstI). The SNPs contained in these fragments were detected from the 20X sequence data of the individuals in our population. Imputation accuracy on HD chip with these genotypes was assessed as the mean correlation between true and imputed genotypes. Several production traits were evaluated using single-step GBLUP methodology. The impact of imputation errors on the ranking of the selection candidates was assessed by comparing a genomic evaluation based on ancestry using true HD or imputed HD genotyping. The relative accuracy of genomic estimated breeding values (GEBVs) was investigated by considering the GEBVs estimated on offspring as a reference. With AvaII or PstI and ddRADseq with TaqI and PstI, more than 10 K SNPs were detected in common with the HD SNP chip, resulting in an imputation accuracy greater than 0.97. The impact of imputation errors on genomic evaluation of the breeders was reduced, with a Spearman correlation greater than 0.99. Finally, the relative accuracy of GEBVs was equivalent. CONCLUSIONS: RADseq approaches can be interesting alternatives to low-density SNP chips for genomic selection. With more than 10 K SNPs in common with the SNPs of the HD SNP chip, good imputation and genomic evaluation results can be obtained. However, with real data, heterogeneity between individuals with missing data must be considered.

Watch Out for a Second SNP: Focus on Multi-Nucleotide Variants in Coding Regions and Rescued Stop-Gained.

Most single-nucleotide polymorphisms (SNPs) are located in non-coding regions, but the fraction usually studied is harbored in protein-coding regions because potential impacts on proteins are relatively easy to predict by popular tools such as the Variant Effect Predictor. These tools annotate variants independently without considering the potential effect of grouped or haplotypic variations, often called "multi-nucleotide variants" (MNVs). Here, we used a large RNA-seq dataset to survey MNVs, comprising 382 chicken samples originating from 11 populations analyzed in the companion paper in which 9.5M SNPs- including 3.3M SNPs with reliable genotypes-were detected. We focused our study on in-codon MNVs and evaluate their potential mis-annotation. Using GATK HaplotypeCaller read-based phasing results, we identified 2,965 MNVs observed in at least five individuals located in 1,792 genes. We found 41.1% of them showing a novel impact when compared to the effect of their constituent SNPs analyzed separately. The biggest impact variation flux concerns the originally annotated stop-gained consequences, for which around 95% were rescued; this flux is followed by the missense consequences for which 37% were reannotated with a different amino acid. We then present in more depth the rescued stop-gained MNVs and give an illustration in the SLC27A4 gene. As previously shown in human datasets, our results in chicken demonstrate the value of haplotype-aware variant annotation, and the interest to consider MNVs in the coding region, particularly when searching for severe functional consequence such as stop-gained variants.

RNA-Seq Data for Reliable SNP Detection and Genotype Calling: Interest for Coding Variant Characterization and Cis-Regulation Analysis by Allele-Specific Expression in Livestock Species.

In addition to their common usages to study gene expression, RNA-seq data accumulated over the last 10 years are a yet-unexploited resource of SNPs in numerous individuals from different populations. SNP detection by RNA-seq is particularly interesting for livestock species since whole genome sequencing is expensive and exome sequencing tools are unavailable. These SNPs detected in expressed regions can be used to characterize variants affecting protein functions, and to study cis-regulated genes by analyzing allele-specific expression (ASE) in the tissue of interest. However, gene expression can be highly variable, and filters for SNP detection using the popular GATK toolkit are not yet standardized, making SNP detection and genotype calling by RNA-seq a challenging endeavor. We compared SNP calling results using GATK suggested filters, on two chicken populations for which both RNA-seq and DNA-seq data were available for the same samples of the same tissue. We showed, in expressed regions, a RNA-seq precision of 91% (SNPs detected by RNA-seq and shared by DNA-seq) and we characterized the remaining 9% of SNPs. We then studied the genotype (GT) obtained by RNA-seq and the impact of two factors (GT call-rate and read number per GT) on the concordance of GT with DNA-seq; we proposed thresholds for them leading to a 95% concordance. Applying these thresholds to 767 multi-tissue RNA-seq of 382 birds of 11 chicken populations, we found 9.5 M SNPs in total, of which ∼550,000 SNPs per tissue and population with a reliable GT (call rate ≥ 50%) and among them, ∼340,000 with a MAF ≥ 10%. We showed that such RNA-seq data from one tissue can be used to (i) detect SNPs with a strong predicted impact on proteins, despite their scarcity in each population (16,307 SIFT deleterious missenses and 590 stop-gained), (ii) study, on a large scale, cis-regulations of gene expression, with ∼81% of protein-coding and 68% of long non-coding genes (TPM ≥ 1) that can be analyzed for ASE, and with ∼29% of them that were cis-regulated, and (iii) analyze population genetic using such SNPs located in expressed regions. This work shows that RNA-seq data can be used with good confidence to detect SNPs and associated GT within various populations and used them for different analyses as GTEx studies.

An integrative atlas of chicken long non-coding genes and their annotations across 25 tissues.

Long non-coding RNAs (LNC) regulate numerous biological processes. In contrast to human, the identification of LNC in farm species, like chicken, is still lacunar. We propose a catalogue of 52,075 chicken genes enriched in LNC ( http://www.fragencode.org/ ), built from the Ensembl reference extended using novel LNC modelled here from 364 RNA-seq and LNC from four public databases. The Ensembl reference grew from 4,643 to 30,084 LNC, of which 59% and 41% with expression ≥ 0.5 and ≥ 1 TPM respectively. Characterization of these LNC relatively to the closest protein coding genes (PCG) revealed that 79% of LNC are in intergenic regions, as in other species. Expression analysis across 25 tissues revealed an enrichment of co-expressed LNC:PCG pairs, suggesting co-regulation and/or co-function. As expected LNC were more tissue-specific than PCG (25% vs. 10%). Similarly to human, 16% of chicken LNC hosted one or more miRNA. We highlighted a new chicken LNC, hosting miR155, conserved in human, highly expressed in immune tissues like miR155, and correlated with immunity-related PCG in both species. Among LNC:PCG pairs tissue-specific in the same tissue, we revealed an enrichment of divergent pairs with the PCG coding transcription factors, as for example LHX5, HXD3 and TBX4, in both human and chicken.

Long noncoding RNAs in lipid metabolism: literature review and conservation analysis across species.

BACKGROUND: Lipids are important for the cell and organism life since they are major components of membranes, energy reserves and are also signal molecules. The main organs for the energy synthesis and storage are the liver and adipose tissue, both in humans and in more distant species such as chicken. Long noncoding RNAs (lncRNAs) are known to be involved in many biological processes including lipid metabolism. RESULTS: In this context, this paper provides the most exhaustive list of lncRNAs involved in lipid metabolism with 60 genes identified after an in-depth analysis of the bibliography, while all "review" type articles list a total of 27 genes. These 60 lncRNAs are mainly described in human or mice and only a few of them have a precise described mode-of-action. Because these genes are still named in a non-standard way making such a study tedious, we propose a standard name for this list according to the rules dictated by the HUGO consortium. Moreover, we identified about 10% of lncRNAs which are conserved between mammals and chicken and 2% between mammals and fishes. Finally, we demonstrated that two lncRNA were wrongly considered as lncRNAs in the literature since they are 3' extensions of the closest coding gene. CONCLUSIONS: Such a lncRNAs catalogue can participate to the understanding of the lipid metabolism regulators; it can be useful to better understand the genetic regulation of some human diseases (obesity, hepatic steatosis) or traits of economic interest in livestock species (meat quality, carcass composition). We have no doubt that this first set will be rapidly enriched in coming years.

Canine neuropathies: powerful spontaneous models for human hereditary sensory neuropathies.

In humans, hereditary sensory neuropathies (HSN), also known as hereditary sensory and autonomic neuropathies (HSAN), constitute a clinically and genetically heterogeneous group of disorders characterized by progressive sensory loss, often accompanied by chronic skin ulcerations and nail dystrophic changes. To date, although around 20 genes have already been discovered, they do not explain the genetic causes of all patients. In dogs, similar neuropathies are also diagnosed, several breeds being predisposed to specific forms of the disease. Indeed, the breed specificity of most canine genetic diseases is due to the small numbers of founders and high levels of inbreeding. Recent knowledge and tools developed to study the canine genome efficiently allows deciphering the genetic bases of such diseases. To date, a dozen breeds are recognized to develop specific HSN. For the Border collie and hunting dog breeds, the genes involved have recently been discovered. Other affected breeds thus constitute potential genetic models, with new genes to be found in dogs that can be considered as candidate genes for human HSAN/HSN. Here, we review the different forms of human and canine HSAN/HSN and we present a novel form in Fox terrier cases, highlighting the advantages of the dog model for such rare human diseases.

Natural models for retinitis pigmentosa: progressive retinal atrophy in dog breeds.

Retinitis pigmentosa (RP) is a heterogeneous group of inherited retinal disorders eventually leading to blindness with different ages of onset, progression and severity. Human RP, first characterized by the progressive degeneration of rod photoreceptor cells, shows high genetic heterogeneity with more than 90 genes identified. However, about one-third of patients have no known genetic causes. Interestingly, dogs are also severely affected by similar diseases, called progressive retinal atrophy (PRA). Indeed, RP and PRA have comparable clinical signs, physiopathology and outcomes, similar diagnosis methods and most often, orthologous genes are involved. The many different dog PRAs often segregate in specific breeds. Indeed, undesired alleles have been selected and amplified through drastic selection and excessive use of inbreeding. Out of the 400 breeds, nearly 100 have an inherited form of PRA, which are natural animal models that can be used to investigate the genetics, disease progression and therapies in dogs for the benefit of both dogs and humans. Recent knowledge on the canine genome and access to new genotyping and sequencing technologies now efficiently allows the identification of mutations involved in canine genetic diseases. To date, PRA genes identified in dog breeds correspond to the same genes in humans and represent relevant RP models, and new genes found in dogs represent good candidate for still unknown human RP. We present here a review of the main advantages of the dog models for human RP with the genes already identified and an X-linked PRA in the Border collie as a model for orphan X-linked RPs in human.

Characterisation and functional predictions of canine long non-coding RNAs.

Long non-coding RNAs (lncRNAs) are a family of heterogeneous RNAs that play major roles in multiple biological processes. We recently identified an extended repertoire of more than 10,000 lncRNAs of the domestic dog however, predicting their biological functionality remains challenging. In this study, we have characterised the expression profiles of 10,444 canine lncRNAs in 26 distinct tissue types, representing various anatomical systems. We showed that lncRNA expressions are mainly clustered by tissue type and we highlighted that 44% of canine lncRNAs are expressed in a tissue-specific manner. We further demonstrated that tissue-specificity correlates with specific families of canine transposable elements. In addition, we identified more than 900 conserved dog-human lncRNAs for which we show their overall reproducible expression patterns between dog and human through comparative transcriptomics. Finally, co-expression analyses of lncRNA and neighbouring protein-coding genes identified more than 3,400 canine lncRNAs, suggesting that functional roles of these lncRNAs act as regulatory elements. Altogether, this genomic and transcriptomic integrative study of lncRNAs constitutes a major resource to investigate genotype to phenotype relationships and biomedical research in the dog species.

Discovery of Human-Similar Gene Fusions in Canine Cancers.

Canine cancers represent a tremendous natural resource due to their incidence and striking similarities to human cancers, sharing similar clinical and pathologic features as well as oncogenic events, including identical somatic mutations. Considering the importance of gene fusions as driver alterations, we explored their relevance in canine cancers. We focused on three distinct human-comparable canine cancers representing different tissues and embryonic origins. Through RNA-Seq, we discovered similar gene fusions as those found in their human counterparts: IGK-CCND3 in B-cell lymphoma, MPB-BRAF in glioma, and COL3A1-PDGFB in dermatofibrosarcoma protuberans-like. We showed not only similar partner genes but also identical breakpoints leading to oncogene overexpression. This study demonstrates similar gene fusion partners and mechanisms in human-dog corresponding tumors and allows for selection of targeted therapies in preclinical and clinical trials with pet dogs prior to human trials, within the framework of personalized medicine. Cancer Res; 77(21); 5721-7. ©2017 AACR.

FEELnc: a tool for long non-coding RNA annotation and its application to the dog transcriptome.

Whole transcriptome sequencing (RNA-seq) has become a standard for cataloguing and monitoring RNA populations. One of the main bottlenecks, however, is to correctly identify the different classes of RNAs among the plethora of reconstructed transcripts, particularly those that will be translated (mRNAs) from the class of long non-coding RNAs (lncRNAs). Here, we present FEELnc (FlExible Extraction of LncRNAs), an alignment-free program that accurately annotates lncRNAs based on a Random Forest model trained with general features such as multi k-mer frequencies and relaxed open reading frames. Benchmarking versus five state-of-the-art tools shows that FEELnc achieves similar or better classification performance on GENCODE and NONCODE data sets. The program also provides specific modules that enable the user to fine-tune classification accuracy, to formalize the annotation of lncRNA classes and to identify lncRNAs even in the absence of a training set of non-coding RNAs. We used FEELnc on a real data set comprising 20 canine RNA-seq samples produced by the European LUPA consortium to substantially expand the canine genome annotation to include 10 374 novel lncRNAs and 58 640 mRNA transcripts. FEELnc moves beyond conventional coding potential classifiers by providing a standardized and complete solution for annotating lncRNAs and is freely available at https://github.com/tderrien/FEELnc.

Amy2B copy number variation reveals starch diet adaptations in ancient European dogs.

Extant dog and wolf DNA indicates that dog domestication was accompanied by the selection of a series of duplications on the Amy2B gene coding for pancreatic amylase. In this study, we used a palaeogenetic approach to investigate the timing and expansion of the Amy2B gene in the ancient dog populations of Western and Eastern Europe and Southwest Asia. Quantitative polymerase chain reaction was used to estimate the copy numbers of this gene for 13 ancient dog samples, dated to between 15 000 and 4000 years before present (cal. BP). This evidenced an increase of Amy2B copies in ancient dogs from as early as the 7th millennium cal. BP in Southeastern Europe. We found that the gene expansion was not fixed across all dogs within this early farming context, with ancient dogs bearing between 2 and 20 diploid copies of the gene. The results also suggested that selection for the increased Amy2B copy number started 7000 years cal. BP, at the latest. This expansion reflects a local adaptation that allowed dogs to thrive on a starch rich diet, especially within early farming societies, and suggests a biocultural coevolution of dog genes and human culture.

A Point Mutation in a lincRNA Upstream of GDNF Is Associated to a Canine Insensitivity to Pain: A Spontaneous Model for Human Sensory Neuropathies.

Human Hereditary Sensory Autonomic Neuropathies (HSANs) are characterized by insensitivity to pain, sometimes combined with self-mutilation. Strikingly, several sporting dog breeds are particularly affected by such neuropathies. Clinical signs appear in young puppies and consist of acral analgesia, with or without sudden intense licking, biting and severe self-mutilation of the feet, whereas proprioception, motor abilities and spinal reflexes remain intact. Through a Genome Wide Association Study (GWAS) with 24 affected and 30 unaffected sporting dogs using the Canine HD 170K SNP array (Illumina), we identified a 1.8 Mb homozygous locus on canine chromosome 4 (adj. p-val = 2.5x10-6). Targeted high-throughput sequencing of this locus in 4 affected and 4 unaffected dogs identified 478 variants. Only one variant perfectly segregated with the expected recessive inheritance in 300 sporting dogs of known clinical status, while it was never present in 900 unaffected dogs from 130 other breeds. This variant, located 90 kb upstream of the GDNF gene, a highly relevant neurotrophic factor candidate gene, lies in a long intergenic non-coding RNAs (lincRNA), GDNF-AS. Using human comparative genomic analysis, we observed that the canine variant maps onto an enhancer element. Quantitative RT-PCR of dorsal root ganglia RNAs of affected dogs showed a significant decrease of both GDNF mRNA and GDNF-AS expression levels (respectively 60% and 80%), as compared to unaffected dogs. We thus performed gel shift assays (EMSA) that reveal that the canine variant significantly alters the binding of regulatory elements. Altogether, these results allowed the identification in dogs of GDNF as a relevant candidate for human HSAN and insensitivity to pain, but also shed light on the regulation of GDNF transcription. Finally, such results allow proposing these sporting dog breeds as natural models for clinical trials with a double benefit for human and veterinary medicine.

Evidence of coat color variation sheds new light on ancient canids.

We have used a paleogenetics approach to investigate the genetic landscape of coat color variation in ancient Eurasian dog and wolf populations. We amplified DNA fragments of two genes controlling coat color, Mc1r (Melanocortin 1 Receptor) and CBD103 (canine-ß-defensin), in respectively 15 and 19 ancient canids (dogs and wolf morphotypes) from 14 different archeological sites, throughout Asia and Europe spanning from ca. 12 000 B.P. (end of Upper Palaeolithic) to ca. 4000 B.P. (Bronze Age). We provide evidence of a new variant (R301C) of the Melanocortin 1 receptor (Mc1r) and highlight the presence of the beta-defensin melanistic mutation (CDB103-K locus) on ancient DNA from dog-and wolf-morphotype specimens. We show that the dominant K(B) allele (CBD103), which causes melanism, and R301C (Mc1r), the variant that may cause light hair color, are present as early as the beginning of the Holocene, over 10,000 years ago. These results underline the genetic diversity of prehistoric dogs. This diversity may have partly stemmed not only from the wolf gene pool captured by domestication but also from mutations very likely linked to the relaxation of natural selection pressure occurring in-line with this process.

PNPLA1 mutations cause autosomal recessive congenital ichthyosis in golden retriever dogs and humans.

Ichthyoses comprise a heterogeneous group of genodermatoses characterized by abnormal desquamation over the whole body, for which the genetic causes of several human forms remain unknown. We used a spontaneous dog model in the golden retriever breed, which is affected by a lamellar ichthyosis resembling human autosomal recessive congenital ichthyoses (ARCI), to carry out a genome-wide association study. We identified a homozygous insertion-deletion (indel) mutation in PNPLA1 that leads to a premature stop codon in all affected golden retriever dogs. We subsequently found one missense and one nonsense mutation in the catalytic domain of human PNPLA1 in six individuals with ARCI from two families. Further experiments highlighted the importance of PNPLA1 in the formation of the epidermal lipid barrier. This study identifies a new gene involved in human ichthyoses and provides insights into the localization and function of this yet uncharacterized member of the PNPLA protein family.

Comparison of buccal and blood-derived canine DNA, either native or whole genome amplified, for array-based genome-wide association studies.

BACKGROUND: The availability of array-based genotyping platforms for single nucleotide polymorphisms (SNPs) for the canine genome has expanded the opportunities to undertake genome-wide association (GWA) studies to identify the genetic basis for Mendelian and complex traits. Whole blood as the source of high quality DNA is undisputed but often proves impractical for collection of the large numbers of samples necessary to discover the loci underlying complex traits. Further, many countries prohibit the collection of blood from dogs unless medically necessary thereby restricting access to critical control samples from healthy dogs. Alternate sources of DNA, typically from buccal cytobrush extractions, while convenient, have been suggested to have low yield and perform poorly in GWA. Yet buccal cytobrushes provide a cost-effective means of collecting DNA, are readily accepted by dog owners, and represent a large resource base in many canine genetics laboratories. To increase the DNA quantities, whole genome amplification (WGA) can be performed. Thus, the present study assessed the utility of buccal-derived DNA as well as whole genome amplification in comparison to blood samples for use on the most recent iteration of the canine HD SNP array (Illumina). FINDINGS: In both buccal and blood samples, whether whole genome amplified or not, 97% of the samples had SNP call rates in excess of 80% indicating that the vast majority of the SNPs would be suitable to perform association studies regardless of the DNA source. Similarly, there were no significant differences in marker intensity measurements between buccal and blood samples for copy number variations (CNV) analysis. CONCLUSIONS: All DNA samples assayed, buccal or blood, native or whole genome amplified, are appropriate for use in array-based genome-wide association studies. The concordance between subsets of dogs for which both buccal and blood samples, or those samples whole genome amplified, was shown to average >99%. Thus, the two DNA sources were comparable in the generation of SNP genotypes and intensity values to estimate structural variation indicating the utility for the use of buccal cytobrush samples and the reliability of whole genome amplification for genome-wide association and CNV studies.

Identification of a gene expression profile associated with operational tolerance among a selected group of stable kidney transplant patients.

Despite their utility, immunosuppressive treatments have numerous side effects, including infectious complications, malignancies and metabolic disorders, all of which contribute to long-term graft loss. In addition to the development of new pharmaceutical products with reduced toxicity and more comfortable modes of administration, tailoring immunosuppression according to the immune status of each patient would represent a significant breakthrough. Gene expression profiling has been shown to be a clinically relevant monitoring tool. In this paper, we have assessed the overall long-term kidney transplant outcome and attempted to identify operationally tolerant-like patients among recipients with stable clinical status at least 5 years post-transplantation. We thus measured a combination of noninvasive blood biomarkers of operational tolerance in a cohort of 144 stable patients and showed that only 3.5% exhibited a gene expression profile of operational tolerance, suggesting that such a profile can be detected under immunosuppressive therapy but that its frequency is low in kidney transplant recipients when compared with liver transplant recipients. We suggest that a rational approach to patient selection, based on a combination of clinical and biological characteristics, may help to provide a safer method for identification of patients potentially suitable for immunosuppressive drug weaning procedures.