Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Identification of population-informative markers from high-density genotyping data through combined feature selection and machine learning algorithms: Application to European autochthonous and cosmopolitan pig breeds.

Large genotyping datasets, obtained from high-density single nucleotide polymorphism (SNP) arrays, developed for different livestock species, can be used to describe and differentiate breeds or populations. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this study, we applied the Boruta algorithm, a wrapper of the machine learning random forest algorithm, on a database of 23 European pig breeds (20 autochthonous and three cosmopolitan breeds) genotyped with a 70k SNP chip, to pre-select informative SNPs. To identify different sets of SNPs, these pre-selected markers were then ranked with random forest based on their mean decrease accuracy and mean decrease gene indexes. We evaluated the efficiency of these subsets for breed classification and the usefulness of this approach to detect candidate genes affecting breed-specific phenotypes and relevant production traits that might differ among breeds. The lowest overall classification error (2.3%) was reached with a subpanel including only 398 SNPs (ranked based on their mean decrease accuracy), with no classification error in seven breeds using up to 49 SNPs. Several SNPs of these selected subpanels were in genomic regions in which previous studies had identified signatures of selection or genes associated with morphological or production traits that distinguish the analysed breeds. Therefore, even if these approaches have not been originally designed to identify signatures of selection, the obtained results showed that they could potentially be useful for this purpose.

Genome-wide association studies for diarrhoea outcomes identified genomic regions affecting resistance to a severe enteropathy in suckling rabbits.

Selection and breeding strategies to improve resistance to enteropathies are essential to reaching the sustainability of the rabbit production systems. However, disease heterogeneity (having only as major visible symptom diarrhoea) and low disease heritability are two barriers for the implementation of these strategies. Diarrhoea condition can affect rabbits at different life stages, starting from the suckling period, with large negative economic impacts. In this study, from a commercial population of suckling rabbits (derived from 133 litters) that experienced an outbreak of enteropathy, we first selected a few animals that died with severe symptoms of diarrhoea and characterized their microbiota, using 16S rRNA gene sequencing data. Clostridium genus was consistently present in all affected specimens. In addition, with the aim to identify genetic markers in the rabbit genome that could be used as selection tools, we performed genome-wide association studies for symptoms of diarrhoea in the same commercial rabbit population. These studies were also complemented with FST analyses between the same groups of rabbits. A total of 332 suckling rabbits (151 with severe symptoms of diarrhoea, 42 with mild symptoms and 129 without any symptoms till the weaning period), derived from 45 different litters (a subset of the 133 litters) were genotyped with the Affymetrix Axiom OrcunSNP Array. In both genomic approaches, rabbits within litters were paired to constitute two groups (susceptible and resistant, including the mildly affected in one or the other group) and run case and control genome-wide association analyses. Genomic heritability estimated in the designed experimental structure integrated in a commercial breeding scheme was 0.19-0.21 (s.e. 0.09-0.10). A total of eight genomic regions on rabbit chromosome 2 (OCU2), OCU3, OCU7, OCU12, OCU13, OCU16 and in an unassembled scaffold had significant single nucleotide polymorphisms (SNPs) and/or markers that trespassed the FST percentile distribution. Among these regions, three main peaks of SNPs were identified on OCU12, OCU13 and OCU16. The QTL region on OCU13 encompasses several genes that encode members of a family of immunoglobulin Fc receptors (FCER1G, FCRLA, FCRLB and FCGR2A) involved in the immune innate system, which might be important candidate genes for this pathogenic condition. The results obtained in this study demonstrated that resistance to an enteropathy occurring in suckling rabbits is in part genetically determined and can be dissected at the genomic level, providing DNA markers that could be used in breeding programmes to increase resistance to enteropathies in meat rabbits.

Persistence of autozygosity in crossbreds between autochthonous and cosmopolitan breeds of swine: a simulation study.

Crossbreeding might be a valid strategy to valorize local pig breeds. Crossbreeding should reduce homozygosity and, as a consequence, yield hybrid vigor for fitness and production traits. This study aimed to quantify the persistence of autozygosity in terminal crossbred pigs compared with purebreds and, in turn, identify genomic regions where autozygosity's persistence would not be found. The study was based on genotyping data from 20 European local pig breeds and three cosmopolitan pig breeds used to simulate crossbred offspring. This study consisted of two steps. First, one hundred matings were simulated for each pairwise combination of the 23 considered breeds (for a total of 276 combinations), ignoring the sex of the parent individuals in order to generate purebred and crossbred matings leveraging all the germplasm available. Second, a few preselected terminal-maternal breed pairs were used to mimic a realistic terminal crossbreeding system: (i) Mora Romagnola (boars) or Cinta Senese (boars) crossed with Large White (sows) or Landrace (sows); (ii) Duroc (boars) crossed with Mora Romagnola (sows) or Cinta Senese (sows). Runs of homozygosity was used to estimate genome-wide autozygosity (FROH). Observed FROH was higher in purebreds than in crossbreds, although some crossbred combinations showed higher FROH than other purebred combinations. Among the purebreds, the highest FROH values were observed in Mora Romagnola and Turopolje (0.50 and 0.46, respectively). FROH ranged from 0.04 to 0.16 in the crossbreds Alentejana × Large White and Alentejana × Iberian, respectively. Persistence of autozygosity was found in several genomic segments harboring regions where quantitative trait loci (QTLs) were found in the literature. The regions were enriched in QTLs involved in fatty acid metabolism and associated with performance traits. This simulation shows that autozygosity persists in most breed combinations of terminal crosses. Results suggest that a strategy for crossbreeding is implemented when leveraging autochthonous and cosmopolitan breeds to obtain most of the hybrid vigor.

Investigation of ABO Gene Variants across More Than 60 Pig Breeds and Populations and Other Suidae Species Using Whole-Genome Sequencing Datasets.

Polymorphisms in the human ABO gene determine the major blood classification system based on the three well-known forms: A; B; and O. In pigs that carry only two main alleles in this gene (A and O), we still need to obtain a more comprehensive distribution of variants, which could also impact its function. In this study, we mined more than 500 whole-genome sequencing datasets to obtain information on the ABO gene in different Suidae species, pig breeds, and populations and provide (i) a comprehensive distribution of the A and O alleles, (ii) evolutionary relationships of ABO gene sequences across Suidae species, and (iii) an exploratory evaluation of the effect of the different ABO gene variants on production traits and blood-related parameters in Italian Large White pigs. We confirmed that allele O is likely under balancing selection, present in all Sus species investigated, without being fixed in any of them. We reported a novel structural variant in perfect linkage disequilibrium with allele O that made it possible to estimate the evolutionary time window of occurrence of this functional allele. We also identified two single nucleotide polymorphisms that were suggestively associated with plasma magnesium levels in pigs. Other studies can also be constructed over our results to further evaluate the effect of this gene on economically relevant traits and basic biological functions.