Genetic predisposition to differentiated thyroid cancer in the Polish population.

INTRODUCTION: Genome sequencing technologies reveal molecular mechanisms of differentiated thyroid cancer (DTC). Unlike somatic mutation analysis from thyroidectomy samples, germline mutations showing genetic susceptibility to DTC are less understood. OBJECTIVES: 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 data. PATIENTS AND METHODS: We analyzed sequencing data from 1076 unrelated individuals totaling over 1018 billion read pairs and yielding an average 35.26 × read depth per genome, 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. RESULTS: 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 was compared with the frequency estimated for the non­Finnish European population obtained from the gnomAD database (gnomad.broadinstitute.org). Significant differences in variant frequency were found for the APC, ARSB, ATM, BRCA1, CHEK2, DICER1, GPD1L, INSR, KCNJ10, MYH9, PALB2, PLCB1, PLEKHG5, PTEN, RET, SEC23B, SERPINA1, SLC26A4, SMAD3, STK11, TERT, TOE1, and WRN genes. CONCLUSIONS: Even though the Polish population is genetically similar to the other European populations, there are significant differences in variant frequencies contributing to the disease development and progression, such as those in the RET, CHEK2, BRCA1, SLC26A4, or TERT genes. Further studies are needed to identify genomic variants associated directly with DTC.

A Study of 41 Canine Orthologues of Human Genes Involved in Monogenic Obesity Reveals Marker in the ADCY3 for Body Weight in Labrador Retrievers.

Obesity and overweight are common conditions in dogs, but individual susceptibility varies with numerous risk factors, including diet, age, sterilization, and gender. In addition to environmental and biological factors, genetic and epigenetic risk factors can influence predisposition to canine obesity, however, they remain unknown. Labrador Retrievers are one of the breeds that are prone to obesity. The purpose of this study was to analyse 41 canine orthologues of human genes linked to monogenic obesity in humans to identify genes associated with body weight in Labrador Retriever dogs. We analysed 11,520 variants from 50 dogs using a linear mixed model with sex, age, and sterilization as covariates and population structure as a random effect. Estimates obtained from the model were subjected to a maxT permutation procedure to adjust p-values for FWER < 0.05. Only the ADCY3 gene showed statistically significant association: TA>T deletion located at 17:19,222,459 in 1/20 intron (per allele effect of 5.56 kg, SE 0.018, p-value = 5.83 × 10-5, TA/TA: 11 dogs; TA/T: 32 dogs; T/T: 7 dogs). Mutations in the ADCY3 gene have already been associated with obesity in mice and humans, making it a promising marker for canine obesity research. Our results provide further evidence that the genetic makeup of obesity in Labrador Retriever dogs contains genes with large effect sizes.

The cancer-risk variant frequency among Polish population reported by the first national whole-genome sequencing study.

Introduction: Population-based cancer screening has raised many controversies in recent years, not only regarding the costs but also regarding the ethical nature and issues related to variant interpretation. Nowadays, genetic cancer screening standards are different in every country and usually encompass only individuals with a personal or family history of relevant cancer. Methods: Here we performed a broad genetic screening for cancer-related rare germline variants on population data from the Thousand Polish Genomes database based on 1076 Polish unrelated individuals that underwent whole genome sequencing (WGS). Results: We identified 19 551 rare variants in 806 genes related to oncological diseases, among them 89% have been located in non-coding regions. The combined BRCA1/BRCA2 pathogenic/likely pathogenic according to ClinVar allele frequency in the unselected population of 1076 Poles was 0.42%, corresponding to nine carriers. Discussion: Altogether, on the population level, we found especially problematic the assessment of the pathogenicity of variants and the relation of ACMG guidelines to the population frequency. Some of the variants may be overinterpreted as disease-causing due to their rarity or lack of annotation in the databases. On the other hand, some relevant variants may have been overseen given that there is little pooled population whole genome data on oncology. Before population WGS screening will become a standard, further studies are needed to assess the frequency of the variants suspected to be pathogenic on the population level and with reporting of likely benign variants.

Better safe than sorry-Whole-genome sequencing indicates that missense variants are significant in susceptibility to COVID-19.

Undoubtedly, genetic factors play an important role in susceptibility and resistance to COVID-19. In this study, we conducted the GWAS analysis. Out of 15,489,173 SNPs, we identified 18,191 significant SNPs for severe and 11,799 SNPs for resistant phenotype, showing that a great number of loci were significant in different COVID-19 representations. The majority of variants were synonymous (60.56% for severe, 58.46% for resistant phenotype) or located in introns (55.77% for severe, 59.83% for resistant phenotype). We identified the most significant SNPs for a severe outcome (in AJAP1 intron) and for COVID resistance (in FIG4 intron). We found no missense variants with a potential causal function on resistance to COVID-19; however, two missense variants were determined as significant a severe phenotype (in PM20D1 and LRP4 exons). None of the aforementioned SNPs and missense variants found in this study have been previously associated with COVID-19.

Population WGS-based spinal muscular atrophy carrier screening in a cohort of 1076 healthy Polish individuals.

Spinal muscular atrophy is a severe neuromuscular disorder with an autosomal recessive inheritance pattern. The disease-causing gene is SMN1, and its paralogue, SMN2, is a disease course modifier. Both genes SMN1 and SMN2 show over 99.9% sequence identity and a high rate of crossing over in the genomic region. Due to this reason, SMN1/SMN2 is usually excluded from the whole-genome sequencing (WGS) analysis and investigated with traditional methods, such as MLPA and qPCR. Recently, novel bioinformatic algorithms dedicated to analyzing this particular genomic region have been developed. Here, we analyze the SMN1/SMN2 genomic region with a dedicated program, SMNCopyNumberCaller. We report a similar prevalence of SMN1 gene deletion carrier status (1 per 41 people) to published data from the Polish population (1 per 35 people). Additionally, SMNCopyNumberCaller can identify SMN2 CNVs and SMN2Δ7-8 present in 153 healthy Polish individuals. Two other programs for the CNV analysis in standard genomic regions were not able to provide reliable results. Using WGS-based tools for SMN1/2 genomic region analysis is not only an efficient method in terms of time but will also enable more complex analysis such screening for markers related with a silent carrier status and identification of further genetic modifiers. Although still an experimental method, soon WGS-based SMN1/SMN2 carrier identification may become a standard method for patients screened with WGS for other purposes.

Gene Variants Related to Cardiovascular and Pulmonary Diseases May Correlate with Severe Outcome of COVID-19.

Background: Severe outcomes of COVID-19 account for up to 15% of all cases. The study aims to check if any gene variants related to cardiovascular (CVD) and pulmonary diseases (PD) are correlated with a severe outcome of COVID-19 in a Polish cohort of COVID-19 patients. Methods: In this study, a subset of 747 samples from unrelated individuals collected across Poland in 2020 and 2021 was used and whole-genome sequencing was performed. Results: The GWAS analysis of SNPs and short indels located in genes related to CVD identified one variant significant in COVID-19 severe outcome in the HADHA gene, while for the PD gene panel, we found two significant variants in the DRC1 gene. In this study, both potentially protective and risk variants were identified, of which variants in the HADHA gene deserve the most attention. Conclusions: This is the first study reporting the association between the HADHA and DRC1 genetic variants and COVID-19 severe outcome based on the cohort WGS analysis. Although all the identified variants are localised in introns, they may be correlated and therefore inherited along with other risk variants, potentially causative to severe outcome of COVID-19 but not discovered yet.

Beyond GWAS-Could Genetic Differentiation within the Allograft Rejection Pathway Shape Natural Immunity to COVID-19?

COVID-19 infections pose a serious global health concern so it is crucial to identify the biomarkers for the susceptibility to and resistance against this disease that could help in a rapid risk assessment and reliable decisions being made on patients' treatment and their potential hospitalisation. Several studies investigated the factors associated with severe COVID-19 outcomes that can be either environmental, population based, or genetic. It was demonstrated that the genetics of the host plays an important role in the various immune responses and, therefore, there are different clinical presentations of COVID-19 infection. In this study, we aimed to use variant descriptive statistics from GWAS (Genome-Wide Association Study) and variant genomic annotations to identify metabolic pathways that are associated with a severe COVID-19 infection as well as pathways related to resistance to COVID-19. For this purpose, we applied a custom-designed mixed linear model implemented into custom-written software. Our analysis of more than 12.5 million SNPs did not indicate any pathway that was significant for a severe COVID-19 infection. However, the Allograft rejection pathway (hsa05330) was significant (p = 0.01087) for resistance to the infection. The majority of the 27 SNP marking genes constituting the Allograft rejection pathway were located on chromosome 6 (19 SNPs) and the remainder were mapped to chromosomes 2, 3, 10, 12, 20, and X. This pathway comprises several immune system components crucial for the self versus non-self recognition, but also the components of antiviral immunity. Our study demonstrated that not only single variants are important for resistance to COVID-19, but also the cumulative impact of several SNPs within the same pathway matters.

CAPN3 c.1746-20C>G variant is hypomorphic for LGMD R1 calpain 3-related.

The investigated intronic CAPN3 variant NM_000070.3:c.1746-20C>G occurs in the Central and Eastern Europe with a frequency of >1% and there are conflicting interpretations on its pathogenicity. We collected data on 14 patients carrying the CAPN3 c.1746-20C>G variant in trans position with another CAPN3 pathogenic/likely pathogenic variant. The patients compound heterozygous for the CAPN3 c.1746-20C>G variant presented a phenotype consistent with calpainopathy of mild/medium severity. This variant is most frequent in the North/West regions of Russia and may originate from that area. Molecular studies revealed that different splicing isoforms are produced in the muscle. We hypothesize that c.1746-20C>G is a hypomorphic variant with a reduction of RNA and protein expression and only individuals having a higher ratio of abnormal isoforms are affected. Reclassification of the CAPN3 variant c.1746-20C>G from variant with a conflicting interpretation of pathogenicity to hypomorphic variant explains many unidentified cases of limb girdle muscular dystrophy R1 calpain 3-related in Eastern and Central Europe.

The Thousand Polish Genomes-A Database of Polish Variant Allele Frequencies.

Although Slavic populations account for over 4.5% of world inhabitants, no centralised, open-source reference database of genetic variation of any Slavic population exists to date. Such data are crucial for clinical genetics, biomedical research, as well as archeological and historical studies. The Polish population, which is homogenous and sedentary in its nature but influenced by many migrations of the past, is unique and could serve as a genetic reference for the Slavic nations. In this study, we analysed whole genomes of 1222 Poles to identify and genotype a wide spectrum of genomic variation, such as small and structural variants, runs of homozygosity, mitochondrial haplogroups, and de novo variants. Common variant analyses showed that the Polish cohort is highly homogenous and shares ancestry with other European populations. In rare variant analyses, we identified 32 autosomal-recessive genes with significantly different frequencies of pathogenic alleles in the Polish population as compared to the non-Finish Europeans, including C2, TGM5, NUP93, C19orf12, and PROP1. The allele frequencies for small and structural variants, calculated for 1076 unrelated individuals, are released publicly as The Thousand Polish Genomes database, and will contribute to the worldwide genomic resources available to researchers and clinicians.

Genetic Variability of Methane Production and Concentration Measured in the Breath of Polish Holstein-Friesian Cattle.

The genetic architecture of methane (CH4) production remains largely unknown. We aimed to estimate its heritability and to perform genome-wide association studies (GWAS) for the identification of candidate genes associated with two phenotypes: CH4 in parts per million/day (CH4 ppm/d) and CH4 in grams/day (CH4 g/d). We studied 483 Polish Holstein-Friesian cows kept on two commercial farms in Poland. Measurements of CH4 and carbon dioxide (CO2) concentrations exhaled by cows during milking were obtained using gas analyzers installed in the automated milking system on the farms. Genomic analyses were performed using a single-step BLUP approach. The percentage of genetic variance explained by SNPs was calculated for each SNP separately and then for the windows of neighbouring SNPs. The heritability of CH4 ppm/d ranged from 0 to 0.14, with an average of 0.085. The heritability of CH4 g/d ranged from 0.13 to 0.26, with an average of 0.22. The GWAS detected potential candidate SNPs on BTA 14 which explained ~0.9% of genetic variance for CH4 ppm/d and ~1% of genetic variance for CH4 g/d. All identified SNPs were located in the TRPS1 gene. We showed that methane traits are partially controlled by genes; however, the detected SNPs explained only a small part of genetic variation-implying that both CH4 ppm/d and CH4 g/d are highly polygenic traits.