Whole-genome Sequence Variants Research Articles

PSAP-Genomic-Regions: A Method Leveraging Population Data to Prioritize Coding and Non-Coding Variants in Whole Genome Sequencing for Rare Disease Diagnosis.

The introduction of Next-Generation Sequencing technologies in the clinics has improved rare disease diagnosis. Nonetheless, for very heterogeneous or very rare diseases, more than half of cases still lack molecular diagnosis. Novel strategies are needed to prioritize variants within a single individual. The Population Sampling Probability (PSAP) method was developed to meet this aim but only for coding variants in exome data. Here, we propose an extension of the PSAP method to the non-coding genome called PSAP-genomic-regions. In this extension, instead of considering genes as testing units (PSAP-genes strategy), we use genomic regions defined over the whole genome that pinpoint potential functional constraints. We conceived an evaluation protocol for our method using artificially generated disease exomes and genomes, by inserting coding and non-coding pathogenic ClinVar variants in large data sets of exomes and genomes from the general population. PSAP-genomic-regions significantly improves the ranking of these variants compared to using a pathogenicity score alone. Using PSAP-genomic-regions, more than 50% of non-coding ClinVar variants were among the top 10 variants of the genome. On real sequencing data from six patients with Cerebral Small Vessel Disease and nine patients with male infertility, all causal variants were ranked in the top 100 variants with PSAP-genomic-regions. By revisiting the testing units used in the PSAP method to include non-coding variants, we have developed PSAP-genomic-regions, an efficient whole-genome prioritization tool which offers promising results for the diagnosis of unresolved rare diseases.

ACMGA: a reference-free multiple-genome alignment pipeline for plant species

BackgroundThe short-read whole-genome sequencing (WGS) approach has been widely applied to investigate the genomic variation in the natural populations of many plant species. With the rapid advancements in long-read sequencing and genome assembly technologies, high-quality genome sequences are available for a group of varieties for many plant species. These genome sequences are expected to help researchers comprehensively investigate any type of genomic variants that are missed by the WGS technology. However, multiple genome alignment (MGA) tools designed by the human genome research community might be unsuitable for plant genomes.ResultsTo fill this gap, we developed the AnchorWave-Cactus Multiple Genome Alignment (ACMGA) pipeline, which improved the alignment of repeat elements and could identify long (> 50 bp) deletions or insertions (INDELs). We conducted MGA using ACMGA and Cactus for 8 Arabidopsis (Arabidopsis thaliana) and 26 Maize (Zea mays) de novo assembled genome sequences and compared them with the previously published short-read variant calling results. MGA identified more single nucleotide variants (SNVs) and long INDELs than did previously published WGS variant callings. Additionally, ACMGA detected significantly more SNVs and long INDELs in repetitive regions and the whole genome than did Cactus. Compared with the results of Cactus, the results of ACMGA were more similar to the previously published variants called using short-read. These two MGA pipelines identified numerous multi-allelic variants that were missed by the WGS variant calling pipeline.ConclusionsAligning denovo assembled genome sequences could identify more SNVs and INDELs than mapping short-read. ACMGA combines the advantages of AnchorWave and Cactus and offers a practical solution for plant MGA by integrating global alignment, a 2-piece-affine-gap cost strategy, and the progressive MGA algorithm.

Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries.

We develop a method, SBayesRC, that integrates genome-wide association study (GWAS) summary statistics with functional genomic annotations to improve polygenic prediction of complex traits. Our method is scalable to whole-genome variant analysis and refines signals from functional annotations by allowing them to affect both causal variant probability and causal effect distribution. We analyze 50 complex traits and diseases using ∼7 million common single-nucleotide polymorphisms (SNPs) and 96 annotations. SBayesRC improves prediction accuracy by 14% in European ancestry and up to 34% in cross-ancestry prediction compared to the baseline method SBayesR, which does not use annotations, and outperforms other methods, including LDpred2, LDpred-funct, MegaPRS, PolyPred-S and PRS-CSx. Investigation of factors affecting prediction accuracy identifies a significant interaction between SNP density and annotation information, suggesting whole-genome sequence variants with annotations may further improve prediction. Functional partitioning analysis highlights a major contribution of evolutionary constrained regions to prediction accuracy and the largest per-SNP contribution from nonsynonymous SNPs.

A comparative study of structural variant calling in WGS from Alzheimer's disease families.

Detecting structural variants (SVs) in whole-genome sequencing poses significant challenges. We present a protocol for variant calling, merging, genotyping, sensitivity analysis, and laboratory validation for generating a high-quality SV call set in whole-genome sequencing from the Alzheimer's Disease Sequencing Project comprising 578 individuals from 111 families. Employing two complementary pipelines, Scalpel and Parliament, for SV/indel calling, we assessed sensitivity through sample replicates (N = 9) with in silico variant spike-ins. We developed a novel metric, D-score, to evaluate caller specificity for deletions. The accuracy of deletions was evaluated by Sanger sequencing. We generated a high-quality call set of 152,301 deletions of diverse sizes. Sanger sequencing validated 114 of 146 detected deletions (78.1%). Scalpel excelled in accuracy for deletions ≤100 bp, whereas Parliament was optimal for deletions >900 bp. Overall, 83.0% and 72.5% of calls by Scalpel and Parliament were validated, respectively, including all 11 deletions called by both Parliament and Scalpel between 101 and 900 bp. Our flexible protocol successfully generated a high-quality deletion call set and a truth set of Sanger sequencing-validated deletions with precise breakpoints spanning 1-17,000 bp.

Ancestral alleles defined for 70 million cattle variants using a population-based likelihood ratio test

BackgroundThe study of ancestral alleles provides insights into the evolutionary history, selection, and genetic structures of a population. In cattle, ancestral alleles are widely used in genetic analyses, including the detection of signatures of selection, determination of breed ancestry, and identification of admixture. Having a comprehensive list of ancestral alleles is expected to improve the accuracy of these genetic analyses. However, the list of ancestral alleles in cattle, especially at the whole genome sequence level, is far from complete. In fact, the current largest list of ancestral alleles (~ 42 million) represents less than 28% of the total number of detected variants in cattle. To address this issue and develop a genomic resource for evolutionary studies, we determined ancestral alleles in cattle by comparing prior derived whole-genome sequence variants to an out-species group using a population-based likelihood ratio test.ResultsOur study determined and makes available the largest list of ancestral alleles in cattle to date (70.1 million) and includes 2.3 million on the X chromosome. There was high concordance (97.6%) of the determined ancestral alleles with those from previous studies when only high-probability ancestral alleles were considered (29.8 million positions) and another 23.5 million high-confidence ancestral alleles were novel, expanding the available reference list to improve the accuracies of genetic analyses involving ancestral alleles. The high concordance of the results with previous studies implies that our approach using genomic sequence variants and a likelihood ratio test to determine ancestral alleles is appropriate.ConclusionsConsidering the high concordance of ancestral alleles across studies, the ancestral alleles determined in this study including those not previously listed, particularly those with high-probability estimates, may be used for further genetic analyses with reasonable accuracy. Our approach that used predetermined variants in species and the likelihood ratio test to determine ancestral alleles is applicable to other species for which sequence level genotypes are available.

RT-PCR genotyping assays to identify SARS-CoV-2 variants in England in 2021: a design and retrospective evaluation study

Whole-genome sequencing (WGS) is the gold standard diagnostic tool to identify and genetically characterise emerging pathogen mutations (variants), but cost, capacity, and timeliness limit its use when large populations need rapidly assessing. We assessed the potential of genotyping assays to provide accurate and timely variant information at scale by retrospectively examining surveillance for SARS-CoV-2 variants in England between March and September, 2021, when genotyping assays were used widely for variant detection. We chose a panel of four RT-PCR genotyping assays to detect circulating variants of SARS-COV-2 in England and developed a decision algorithm to assign a probable SARS-CoV-2 variant to samples using the assay results. We extracted surveillance data from the UK Health Security Agency databases for 115 934 SARS-CoV-2-positive samples (March 1-Sept 6, 2021) when variant information was available from both genotyping and WGS. By comparing the genotyping and WGS variant result, we calculated accuracy metrics (ie, sensitivity, specificity, and positive predictive value [PPV]) and the time difference between the sample collection date and the availability of variant information. We assessed the number of samples with a variant assigned from genotyping or WGS, or both, over time. Genotyping and an initial decision algorithm (April 10-May 11, 2021 data) were accurate for key variant assignment: sensitivities and PPVs were 0·99 (95% CI 0·99-0·99) for the alpha, 1·00 (1·00-1·00) for the beta, and 0·91(0·80-1·00) for the gamma variants; specificities were 0·97 (0·96-0·98), 1·00 (1·00-1·00), and 1·00 (1·00-1·00), respectively. A subsequent decision algorithm over a longer time period (May 27-Sept 6, 2021 data) remained accurate for key variant assignment: sensitivities were 0·91 (95% CI 0·74-1·00) for the beta, 0·98 (0·98-0·99) for the delta, and 0·93(0·81-1·00) for the gammavariants; specificities were 1·00 (1·00-1·00), 0·96 (0·96-0·97), and 1·00 (1·00-1·00), respectively; and PPVs were 0·83 (0·62-1·00), 1·00 (1·00-1·00), and 0·78 (0·59-0·97), respectively. Genotyping produced variant information a medianof 3days (IQR 2-4) after the sample collection date, which was faster than with WGS (9days [8-11]). The flexibility ofgenotyping enabled a nine-times increase in the quantity of samples tested for variants by this method (from5000 to45 000). RT-PCR genotyping assays are suitable for high-throughput variant surveillance and could complement WGS, enabling larger scale testing for known variants and timelier results, with important implications for effective public health responses and disease control globally, especially in settings with low WGS capacity. However, the choice of panels of RT-PCR assays is highly dependent on database information on circulating variants generated by WGS, which could limit the use of genotyping assays when new variants are emerging and spreading rapidly. UK Health Security Agency and National Institute for Health Research Health Protection Research Unit in Emergency Preparedness and Response.

Sequence Analysis of Six Candidate Genes in Miniature Schnauzers with Primary Hypertriglyceridemia.

Miniature Schnauzers are predisposed to primary hypertriglyceridemia (HTG). In this study, we performed whole genome sequencing (WGS) of eight Miniature Schnauzers with primary HTG and screened for risk variants in six HTG candidate genes: LPL, APOC2, APOA5, GPIHBP1, LMF1, and APOE. Variants were filtered to identify those present in ≥2 Miniature Schnauzers with primary HTG and uncommon (<10% allele frequency) in a WGS variant database including 613 dogs from 61 other breeds. Three variants passed filtering: an APOE TATA box deletion, an LMF1 intronic SNP, and a GPIHBP1 missense variant. The APOE and GPIHBP1 variants were genotyped in a cohort of 108 Miniature Schnauzers, including 68 with primary HTG and 40 controls. A multivariable regression model, including age and sex, did not identify an effect of APOE (estimate = 0.18, std. error = 0.14; p = 0.20) or GPIHBP1 genotypes (estimate = -0.26, std. error = 0.42; p = 0.54) on triglyceride concentration. In conclusion, we did not identify a monogenic cause for primary HTG in Miniature Schnauzers in the six genes evaluated. However, if HTG in Miniature Schnauzers is a complex disease resulting from the cumulative effects of multiple variants and environment, the identified variants cannot be ruled out as contributing factors.

GENETIC PLEIOTROPY OF KIDNEY FUNCTION AND SOLUBLE RECEPTOR FOR AGE: THE LONG LIFE FAMILY STUDY

Abstract Chronic kidney disease (CKD) patients have increased inflammation and oxidative stress that may escalate the production of advanced glycation end-products (AGEs), a biomarker implicated in aging. Circulate soluble receptor for AGE (sRAGE) is associated with CKD and inversely associated with the estimated glomerular filtration rate (eGFR). In the Long Life Family Study, eGFR measured by creatinine (eGFRcr) and cystatin-C (eGFRcys) correlated with sRAGE (r~ -0.27, p&amp;lt; 1x10-30). We investigated genetic pleiotropy conducting GWAS for eGFRcr, eGFRcys, and sRAGE and correlated meta-analysis (CMA) on GWAS p-values from whole-genome sequence variants in 4,182 subjects (age range: 24–110). Of the 59 loci identified (p&amp;lt; 5x10-8), 42 were novel eGFR discoveries. We performed TWAS and CMA-TWAS for eGFRcr, eGFRcys, and sRAGE on blood RNA-seq in 1,209 subjects and identified 17 genes with pleiotropic expressions for eGFRs and sRAGE (Bonferroni-p&amp;lt; 2.73x10-6). Several CMA-GWAS gene variants (e.g., SH3GLB1/SELENOF, CENPP) were eQTLs in glomeruli and tubule kidney (Human Kidney eQTL Atlas, HK-Atlas) and other tissues (GTEx-v8). The bioinformatics analyses identified (i)genes harboring eQTLs in glomeruli and tubule kidney (HK-Atlas), (ii)enhancer-promoter gene variants associated with kidney function-related phenotypes at genome-wide significance (GeneHancer/GeneCards/GWAS-catalog), and (iii)enhancer regions predicting gene expressions in kidney related-function (EPRI, enhancer-promoter-RNA maps). Some functional genes were found among bioinformatics (e.g., MTX1, TOPORS, SNTB1, LSP1, SNTB1, ATG2A, TP53INP2) and others were unique (e.g., TEC, CCN1, and ERGIC1, GeneHancer/GeneCards/GWAS-catalog). Genomic, transcriptomic, and bioinformatics results demonstrated evidence that some novel variants and genes have functional regulatory features in the kidney and are likely involved in biological and clinical aging.

AGIDB: a versatile database for genotype imputation and variant decoding across species.

The high cost of large-scale, high-coverage whole-genome sequencing has limited its application in genomics and genetics research. The common approach has been to impute whole-genome sequence variants obtained from a few individuals for a larger population of interest individually genotyped using SNP chip. An alternative involves low-coverage whole-genome sequencing (lcWGS) of all individuals in the larger population, followed by imputation to sequence resolution. To overcome limitations of processing lcWGS data and meeting specific genotype imputation requirements, we developed AGIDB (https://agidb.pro), a website comprising tools and database with an unprecedented sample size and comprehensive variant decoding for animals. AGIDB integrates whole-genome sequencing and chip data from 17360 and 174945 individuals, respectively, across 89 species to identify over one billion variants, totaling a massive 688.57 TB of processed data. AGIDB focuses on integrating multiple genotype imputation scenarios. It also provides user-friendly searching and data analysis modules that enable comprehensive annotation of genetic variants for specific populations. To meet a wide range of research requirements, AGIDB offers downloadable reference panels for each species in addition to its extensive dataset, variant decoding and utility tools. We hope that AGIDB will become a key foundational resource in genetics and breeding, providing robust support to researchers.

Genomic prediction based on selective linkage disequilibrium pruning of low-coverage whole-genome sequence variants in a pure Duroc population

BackgroundAlthough the accumulation of whole-genome sequencing (WGS) data has accelerated the identification of mutations underlying complex traits, its impact on the accuracy of genomic predictions is limited. Reliable genotyping data and pre-selected beneficial loci can be used to improve prediction accuracy. Previously, we reported a low-coverage sequencing genotyping method that yielded 11.3 million highly accurate single-nucleotide polymorphisms (SNPs) in pigs. Here, we introduce a method termed selective linkage disequilibrium pruning (SLDP), which refines the set of SNPs that show a large gain during prediction of complex traits using whole-genome SNP data.ResultsWe used the SLDP method to identify and select markers among millions of SNPs based on genome-wide association study (GWAS) prior information. We evaluated the performance of SLDP with respect to three real traits and six simulated traits with varying genetic architectures using two representative models (genomic best linear unbiased prediction and BayesR) on samples from 3579 Duroc boars. SLDP was determined by testing 180 combinations of two core parameters (GWAS P-value thresholds and linkage disequilibrium r2). The parameters for each trait were optimized in the training population by five fold cross-validation and then tested in the validation population. Similar to previous GWAS prior-based methods, the performance of SLDP was mainly affected by the genetic architecture of the traits analyzed. Specifically, SLDP performed better for traits controlled by major quantitative trait loci (QTL) or a small number of quantitative trait nucleotides (QTN). Compared with two commercial SNP chips, genotyping-by-sequencing data, and an unselected whole-genome SNP panel, the SLDP strategy led to significant improvements in prediction accuracy, which ranged from 0.84 to 3.22% for real traits controlled by major or moderate QTL and from 1.23 to 11.47% for simulated traits controlled by a small number of QTN.ConclusionsThe SLDP marker selection method can be incorporated into mainstream prediction models to yield accuracy improvements for traits with a relatively simple genetic architecture, however, it has no significant advantage for traits not controlled by major QTL. The main factors that affect its performance are the genetic architecture of traits and the reliability of GWAS prior information. Our findings can facilitate the application of WGS-based genomic selection.

Multi-line ssGBLUP evaluation using preselected markers from whole-genome sequence data in pigs.

Genomic evaluations in pigs could benefit from using multi-line data along with whole-genome sequencing (WGS) if the data are large enough to represent the variability across populations. The objective of this study was to investigate strategies to combine large-scale data from different terminal pig lines in a multi-line genomic evaluation (MLE) through single-step GBLUP (ssGBLUP) models while including variants preselected from whole-genome sequence (WGS) data. We investigated single-line and multi-line evaluations for five traits recorded in three terminal lines. The number of sequenced animals in each line ranged from 731 to 1,865, with 60k to 104k imputed to WGS. Unknown parent groups (UPG) and metafounders (MF) were explored to account for genetic differences among the lines and improve the compatibility between pedigree and genomic relationships in the MLE. Sequence variants were preselected based on multi-line genome-wide association studies (GWAS) or linkage disequilibrium (LD) pruning. These preselected variant sets were used for ssGBLUP predictions without and with weights from BayesR, and the performances were compared to that of a commercial porcine single-nucleotide polymorphisms (SNP) chip. Using UPG and MF in MLE showed small to no gain in prediction accuracy (up to 0.02), depending on the lines and traits, compared to the single-line genomic evaluation (SLE). Likewise, adding selected variants from the GWAS to the commercial SNP chip resulted in a maximum increase of 0.02 in the prediction accuracy, only for average daily feed intake in the most numerous lines. In addition, no benefits were observed when using preselected sequence variants in multi-line genomic predictions. Weights from BayesR did not help improve the performance of ssGBLUP. This study revealed limited benefits of using preselected whole-genome sequence variants for multi-line genomic predictions, even when tens of thousands of animals had imputed sequence data. Correctly accounting for line differences with UPG or MF in MLE is essential to obtain predictions similar to SLE; however, the only observed benefit of an MLE is to have comparable predictions across lines. Further investigation into the amount of data and novel methods to preselect whole-genome causative variants in combined populations would be of significant interest.

GWAS on Imputed Whole-Genome Sequence Variants Reveal Genes Associated with Resistance to Piscirickettsia salmonis in Rainbow Trout (Oncorhynchus mykiss).

Genome-wide association studies (GWAS) allow the identification of associations between genetic variants and important phenotypes in domestic animals, including disease-resistance traits. Whole Genome Sequencing (WGS) data can help increase the resolution and statistical power of association mapping. Here, we conduced GWAS to asses he facultative intracellular bacterium Piscirickettsia salmonis, which affects farmed rainbow trout, Oncorhynchus mykiss, in Chile using imputed genotypes at the sequence level and searched for candidate genes located in genomic regions associated with the trait. A total of 2130 rainbow trout were intraperitoneally challenged with P. salmonis under controlled conditions and genotyped using a 57K single nucleotide polymorphism (SNP) panel. Genotype imputation was performed in all the genotyped animals using WGS data from 102 individuals. A total of 488,979 imputed WGS variants were available in the 2130 individuals after quality control. GWAS revealed genome-wide significant quantitative trait loci (QTL) in Omy02, Omy03, Omy25, Omy26 and Omy27 for time to death and in Omy26 for binary survival. Twenty-four (24) candidate genes associated with P. salmonis resistance were identified, which were mainly related to phagocytosis, innate immune response, inflammation, oxidative response, lipid metabolism and apoptotic process. Our results provide further knowledge on the genetic variants and genes associated with resistance to intracellular bacterial infection in rainbow trout.

Exploring the optimal strategy of imputation from SNP array to whole-genome sequencing data in farm animals.

Genotype imputation from BeadChip to whole-genome sequencing (WGS) data is a cost-effective method of obtaining genotypes of WGS variants. Beagle, one of the most popular imputation software programs, has been widely used for genotype inference in humans and non-human species. A few studies have systematically and comprehensively compared the performance of beagle versions and parameter settings of farm animals. Here, we investigated the imputation performance of three representative versions of Beagle (Beagle 4.1, Beagle 5.0, and Beagle 5.4), and the effective population size (Ne) parameter setting for three species (cattle, pig, and chicken). Six scenarios were investigated to explore the impact of certain key factors on imputation performance. The results showed that the default Ne (1,000,000) is not suitable for livestock and poultry in small reference or low-density arrays of target panels, with 2.47%–10.45% drops in accuracy. Beagle 5 significantly reduced the computation time (4.66-fold–13.24-fold) without an accuracy loss. In addition, using a large combined-reference panel or high-density chip provides greater imputation accuracy, especially for low minor allele frequency (MAF) variants. Finally, a highly significant correlation in the measures of imputation accuracy can be obtained with an MAF equal to or greater than 0.05.

Improving Genomic Selection for Heat Tolerance in Dairy Cattle: Current Opportunities and Future Directions.

Heat tolerance is the ability of an animal to maintain production and reproduction levels under hot and humid conditions and is now a trait of economic relevance in dairy systems worldwide because of an escalating warming climate. The Australian dairy population is one of the excellent study models for enhancing our understanding of the biology of heat tolerance because they are predominantly kept outdoors on pastures where they experience direct effects of weather elements (e.g., solar radiation). In this article, we focus on evidence from recent studies in Australia that leveraged large a dataset [∼40,000 animals with phenotypes and 15 million whole-genome sequence variants] to elucidate the genetic basis of thermal stress as a critical part of the strategy to breed cattle adapted to warmer environments. Genotype-by-environment interaction (i.e., G × E) due to temperature and humidity variation is increasing, meaning animals are becoming less adapted (i.e., more sensitive) to changing environments. There are opportunities to reverse this trend and accelerate adaptation to warming climate by 1) selecting robust or heat-resilient animals and 2) including resilience indicators in breeding goals. Candidate causal variants related to the nervous system and metabolic functions are relevant for heat tolerance and, therefore, key for improving this trait. This could include adding these variants in the custom SNP panels used for routine genomic evaluations or as the basis to design specific agonist or antagonist compounds for lowering core body temperature under heat stress conditions. Indeed, it was encouraging to see that adding prioritized functionally relevant variants into the 50k SNP panel (i.e., the industry panel used for genomic evaluation in Australia) increased the prediction accuracy of heat tolerance by up to 10% units. This gain in accuracy is critical because genetic improvement has a linear relationship with prediction accuracy. Overall, while this article used data mainly from Australia, this could benefit other countries that aim to develop breeding values for heat tolerance, considering that the warming climate is becoming a topical issue worldwide.

Performance analysis of a novel artificial intelligence-based computational method on published ex vivo drug sensitivity data to support targeted treatment decisions in precision oncology.

e13618 Background: Comprehensive molecular profiling is readily available for clinical practice. An extensive amount of published evidence provides information about the potential functional relevance of many driver genes and genetic alterations. But, due to the large number of driver genes and alterations, and the long-tail frequency distribution of these alterations and their possible combinations, very few single driver alterations have reached high-enough-level evidence alone to support clinical decisions. We hypothesized that by aggregating evidence following logical principles of reasoning by a computational system, this vast molecular information can be used to improve personalized treatment decisions. In the present study, we used previously published ex vivo drug sensitivity data to analyze the performance of an artificial intelligence-based computational method, the digital drug assignment (DDA), which has shown utility by improving treatment decisions in case of complex molecular profiles acquired in the SHIVA01 trial. Methods: We selected 111 cases with whole-genome sequencing (WGS) and ex vivo drug sensitivity data of a previously published acute myeloid leukemia study (Tyner et al, 2018). WGS variants were filtered for a preselected hematology-related panel of 446 genes and uploaded to a DDA-based software system to calculate the aggregated evidence level (AEL) values of associated molecularly targeted agents. DDA-predicted sensitivity (or resistance) was defined as AEL &gt; 0 (or AEL &lt; 0) in the presence of at least one actionable driver. Area under the curve (AUC) values were used for determining the ex vivo sensitivity or resistance of leukemia cells to 40 approved drugs and 53 developmental compounds. Results: The AUC values were significantly different in the drug-sensitive and -resistant groups forecasted by the DDA (167.1 and 205.5, respectively, p &lt; 0.0001) and differed significantly from the average AUC value (194.6). Overall, sensitivity was correctly predicted in 66% of compound-sample pairs (n = 671). 88 approved drugs had AEL value over 1000; of these 73% were effective according to the ex vivo results. While forecasted resistance was confirmed in 63% of the cases. With only the most sensitive/resistant 20% of cases considered from the ex vivo data, sensitivity was accurately predicted in 75% of approved compound-sample pairs (n = 173). 37 approved drugs had AEL value over 1000; of these 81% were confirmed as sensitive. Conclusions: The DDA-based computational reasoning has a promising performance in forecasting sensitivity and resistance to a broad spectrum of targeted agents based on molecular information. Therefore, it has the potential to automate, standardize and improve complex molecular profile-based targeted treatment decisions.

Contribution of rare whole-genome sequencing variants to plasma protein levels and the missing heritability

Despite the success of genome-wide association studies, much of the genetic contribution to complex traits remains unexplained. Here, we analyse high coverage whole-genome sequencing data, to evaluate the contribution of rare genetic variants to 414 plasma proteins. The frequency distribution of genetic variants is skewed towards the rare spectrum, and damaging variants are more often rare. We estimate that less than 4.3% of the narrow-sense heritability is expected to be explained by rare variants in our cohort. Using a gene-based approach, we identify Cis-associations for 237 of the proteins, which is slightly more compared to a GWAS (N = 213), and we identify 34 associated loci in Trans. Several associations are driven by rare variants, which have larger effects, on average. We therefore conclude that rare variants could be of importance for precision medicine applications, but have a more limited contribution to the missing heritability of complex diseases.

Functionally prioritised whole-genome sequence variants improve the accuracy of genomic prediction for heat tolerance

BackgroundHeat tolerance is a trait of economic importance in the context of warm climates and the effects of global warming on livestock production, reproduction, health, and well-being. This study investigated the improvement in prediction accuracy for heat tolerance when selected sets of sequence variants from a large genome-wide association study (GWAS) were combined with a standard 50k single nucleotide polymorphism (SNP) panel used by the dairy industry.MethodsOver 40,000 dairy cattle with genotype and phenotype data were analysed. The phenotypes used to measure an individual’s heat tolerance were defined as the rate of decline in milk production traits with rising temperature and humidity. We used Holstein and Jersey cows to select sequence variants linked to heat tolerance. The prioritised sequence variants were the most significant SNPs passing a GWAS p-value threshold selected based on sliding 100-kb windows along each chromosome. We used a bull reference set to develop the genomic prediction equations, which were then validated in an independent set of Holstein, Jersey, and crossbred cows. Prediction analyses were performed using the BayesR, BayesRC, and GBLUP methods.ResultsThe accuracy of genomic prediction for heat tolerance improved by up to 0.07, 0.05, and 0.10 units in Holstein, Jersey, and crossbred cows, respectively, when sets of selected sequence markers from Holstein cows were added to the 50k SNP panel. However, in some scenarios, the prediction accuracy decreased unexpectedly with the largest drop of − 0.10 units for the heat tolerance fat yield trait observed in Jersey cows when 50k plus pre-selected SNPs from Holstein cows were used. Using pre-selected SNPs discovered on a combined set of Holstein and Jersey cows generally improved the accuracy, especially in the Jersey validation. In addition, combining Holstein and Jersey bulls in the reference set generally improved prediction accuracy in most scenarios compared to using only Holstein bulls as the reference set.ConclusionsInformative sequence markers can be prioritised to improve the genomic prediction of heat tolerance in different breeds. In addition to providing biological insight, these variants could also have a direct application for developing customized SNP arrays or can be used via imputation in current industry SNP panels.

Large-scale association study on daily weight gain in pigs reveals overlap of genetic factors for growth in humans

BackgroundImputation from genotyping array to whole-genome sequence variants using resequencing of representative reference populations enhances our ability to map genetic factors affecting complex phenotypes in livestock species. The accumulation of knowledge about gene function in human and laboratory animals can provide substantial advantage for genomic research in livestock species.ResultsIn this study, 201,388 pigs from three commercial Danish breeds genotyped with low to medium (8.5k to 70k) SNP arrays were imputed to whole genome sequence variants using a two-step approach. Both imputation steps achieved high accuracies, and in total this yielded 26,447,434 markers on 18 autosomes. The average estimated imputation accuracy of markers with minor allele frequency ≥ 0.05 was 0.94. To overcome the memory consumption of running genome-wide association study (GWAS) for each breed, we performed within-breed subpopulation GWAS then within-breed meta-analysis for average daily weight gain (ADG), followed by a multi-breed meta-analysis of GWAS summary statistics. We identified 15 quantitative trait loci (QTL). Our post-GWAS analysis strategy to prioritize of candidate genes including information like gene ontology, mammalian phenotype database, differential expression gene analysis of high and low feed efficiency pig and human GWAS catalog for height, obesity, and body mass index, we proposed MRAP2, LEPROT, PMAIP1, ENSSSCG00000036234, BMP2, ELFN1, LIG4 and FAM155A as the candidate genes with biological support for ADG in pigs.ConclusionOur post-GWAS analysis strategy helped to identify candidate genes not just by distance to the lead SNP but also by multiple sources of biological evidence. Besides, the identified QTL overlap with genes which are known for their association with human growth-related traits. The GWAS with this large data set showed the power to map the genetic factors associated with ADG in pigs and have added to our understanding of the genetics of growth across mammalian species.

Genome-wide association study with imputed whole-genome sequence variants including large deletions for female fertility in 3 Nordic dairy cattle breeds

Fertility is an economically important trait in livestock. Poor fertility in dairy cattle can be due to loss-of-function variants affecting any essential gene that causes early embryonic mortality in homozygotes. To identify fertility-associated quantitative trait loci, we performed single-marker association analyses for 8 fertility traits in Holstein, Jersey, and Nordic Red Dairy cattle using imputed whole-genome sequence variants including SNPs, indels, and large deletion. We then performed stepwise selection of independent markers from GWAS loci using conditional and joint association analyses. From single-marker analyses for fertility traits, we reported genome-wide significant associations of 30,384 SNPs, 178 indels, and 3 deletions in Holstein; 23,481 SNPs, 189 indels, and 13 deletions in Nordic Red; and 17 SNPs in Jersey cattle. Conditional and joint association analyses identified 37 and 23 independent associations in Holstein and Nordic Red Dairy cattle, respectively. Fertility-associated GWAS loci were enriched for developmental and cellular processes (Gene Ontology enrichment, false discovery rate < 0.05). For these quantitative trait loci regions (top marker and 500 kb of surrounding regions), we proposed several candidate genes with functional annotations corresponding to embryonic lethality and various fertility-related phenotypes in mouse and cattle. The inclusion of these top markers in future releases of the custom SNP chip used for genomic evaluations will enable their validation in independent populations and improve the accuracy of genomic predictions.

A Next Generation Sequencing-Based Protocol for Screening of Variants of Concern in Autism Spectrum Disorder.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with strong genetic influences. There is an increasing demand for ASD genetic testing beyond the traditionally recommended microarray and syndromic autism testing; however, the current whole genome sequencing (WGS) and whole exome sequencing (WES) methods are lacking an academic standard for WGS variant annotation, reporting, and interpretation, tailored towards patients with ASD and offer very limited interpretation for clinical significance. Using WGS data from six family trios, we demonstrate the clinical feasibility and technical implementation of an evidence-based, fully transparent bioinformatics pipeline and report framework for an ASD-focused WGS genetic report. We confirmed a portion of the key variants with Sanger sequencing and provided interpretation with consideration of patients’ clinical symptoms and detailed literature review. Furthermore, we showed that identification of the genetic contributions of ASD core symptoms and comorbidities may promote a better understanding of the ASD pathophysiology, lead to early detection of associated comorbidities, and facilitate pharmacologic intervention based on pathological pathways inferred from the genetic information. We will make the bioinformatics pipeline and interpretation framework publicly available, in an easily accessible format, after validation with a larger cohort. We hope that the present proposed protocol can serve as a starting point to invite discourse and debate to further improve approaches in WGS-based genetic consultation for patients with ASD.