The experiments described in this research article were designed to test the effect of rare variants into genomic prediction in dairy cattle. Common polymorphisms are able to explain only a small proportion of the underlying genetic variation of complex phenotypes. Variants representing functional mutations with large effects on complex phenotypes are expected to be rare due to natural (humans) or artificial (livestock) selection pressure. Therefore, it is important to check whether the use of rare variants could increase the accuracy of ranking of animals by providing the tool for more precise differentiation among the bulls with high additive genetic merit. The goal of our study was to verify whether including rare variants in a genomic selection model allows for a more accurate description of the additive genetic background of traits under selection in dairy cattle. We used the linear mixed model for comparison SNP estimates for Holstein-Friesian cattle of the two data sets – a set containing only single nucleotide polymorphisms defined by minor allele frequency ≥ 0.01, which is routinely used in the Polish genomic evaluation system (46,216 SNPs), and a set containing SNPs selected based only on the call rate (54,378 SNPs). Based on the SNP estimates we also calculated DGV and GEBV and compared them between both data sets. In all the analyses we used production, fertility, conformation and udder health traits. We also assessed the time required for the two most computationally demanding components of genomic selection: preparing genotype data, and estimation of SNP effects between those two data sets. The results of our study indicated that the analysis including rare variants resulted in changes in the individual ranking of the top 100 male and female candidates, but had no effect on the outcome of the quality of EBV prediction as expressed by the Interbull validation test.

(1) Background: The goal of our study was to identify SNPs, metabolic pathways (KEGG), and gene ontology (GO) terms significantly associated with calving and workability; (2) Methods: Based on the EuroGenomics reference data set, we analyzed direct (DCE) and maternal (MCE) calving ease, direct (DSB), and maternal (MSB) stillbirth, milking speed (MSP), and temperament (TEM). We estimated SNP effects using a multi-SNP mixed-model. Further, SNP positions were mapped to genes, and GO terms/KEGG pathways of the corresponding genes were assigned. The estimation of GO term/KEGG pathway effects was based on a mixed-model using SNP effects as dependent variables; (3) Results: The number of significant SNPs comprised 59 for DCE, 25 for DSB and MSP, 17 for MCE and MSB, and 7 for TEM. Significant KEGG pathways were found for MSB (2), TEM (2), and MSP (1), while 11 GO terms were significant for MSP, 10 for DCE, 8 for DSB and TEM, 5 for MCE, and 3 for MSB.; (4) Conclusions: From the perspective of a better understanding of the genomic background of the phenotypes, traits with low heritabilities suggest that the focus should be moved from single genes to metabolic pathways or gene ontologies significant for the phenotype.

The genome sequencing technologies reveal the molecular mechanisms of differentiated thyroid cancer (DTC). Unlike somatic mutation analysis from thyroidectomy samples, germline mutations showing genetic susceptibility to DTC are less understood. The study aimed to assess the prevalence of germline mutations predisposing to DTC in a cohort of Polish individuals based on their whole genome sequencing (WGS) data. We analyzed sequencing data from 1076 unrelated individuals, released openly for academic and clinical research as The Thousand Polish Genomes database (https://1000polishgenomes.com). The list of genes chosen for further analysis was based on the review of previous studies. The cohort contained 104 variants located within the coding and noncoding DNA sequences of 90 genes selected by ClinVar classification as pathogenic and potentially pathogenic. The frequency of variants in the Polish cohort (our study) was compared to the frequency estimated for the non-Finnish European population, obtained from the gnomAD database (gnomad.broadinstitute.org). Statistically significant variants included 23 genes. Even though the Polish population is genetically similar to other European populations, there are significant differences in variant frequencies contributing to the disease development and progression, such as RET, CHECK2, BRCA1, SLC26A4 or TERT. Further studies are needed to identify genomic variants associated directly with DTC.