Exome sequencing identified rare recurrent copy number variants and hereditary breast cancer susceptibility.

Copy number variants (CNVs) are a major source of genetic variation and can disrupt genes or affect gene dosage. They are known to be causal or underlie predisposition to various diseases. However, the role of CNVs in inherited breast cancer susceptibility has not been thoroughly investigated. To address this, we performed whole-exome sequencing based analysis of rare CNVs in 98 high-risk Northern Finnish breast cancer cases. After filtering, selected candidate alleles were validated and characterized with a combination of orthogonal methods, including PCR-based approaches, optical genome mapping and long-read sequencing. This revealed three recurrent alterations: a 31 kb deletion co-occurring with a retrotransposon insertion (delins) in RAD52, a 13.4 kb deletion in HSD17B14 and a 64 kb partial duplication of RAD51C. Notably, all these genes encode proteins involved in pathways previously identified as essential for breast cancer development. Variants were genotyped in geographically matched cases and controls (altogether 278 hereditary and 1983 unselected breast cancer cases, and 1229 controls). The RAD52 delins and HSD17B14 deletion both showed significant enrichment among cases with indications of hereditary disease susceptibility. RAD52 delins was identified in 7/278 cases (2.5%, P = 0.034, OR = 2.86, 95% CI = 1.10-7.45) and HSD17B14 deletion in 8/278 cases (2.9%, P = 0.014, OR = 3.28, 95% CI = 1.31-8.23), the frequency of both variants in the controls being 11/1229 (0.9%). This suggests a role for RAD52 and HSD17B14 in hereditary breast cancer susceptibility. The RAD51C duplication was very rare, identified only in 2/278 of hereditary cases and 2/1229 controls (P = 0.157, OR = 4.45, 95% CI = 0.62-31.70). The identification of recurrent CNVs in these genes, and especially the relatively high frequency of RAD52 and HSD17B14 alterations in the Finnish population, highlights the importance of studying CNVs alongside single nucleotide variants when searching for genetic factors underlying hereditary disease predisposition.

Strong Purifying Selection in Haploid Tissue-Specific Genes of Scots Pine Supports the Masking Theory.

The masking theory states that genes expressed in a haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in a diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intraspecific competition. Evidence from other plant groups is scarce, and results show no support for the masking theory. Here, we have used a gymnosperm Scots pine megagametophyte, a maternally derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue and that this signal of strong selection is not an artifact driven by high expression levels.

Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

Within- and between-population comparisons suggest independently acting selection maintaining parallel clines in Scots pine (Pinus sylvestris).

Parallel clines in traits related to adaptation in a species can be due to independent selection on a pair of traits, or due to selection in one trait resulting in a parallel cline in a correlated trait. To distinguish between the mechanisms giving rise to parallel adaptive population divergence of multiple traits along an environmental gradient we need to study variation, correlations, and selective forces within individual populations along the gradient. In many tree species, budset timing (BST) forms a latitudinal cline, and parallel clinal variation is also found in other seedling traits, such as first-year height (FYH) and fall frost injury (FFI). In this study, we set up a common garden experiment with open pollinated progeny from natural populations of Scots pine (Pinus sylvestris), with one large sample from single population (500 families) and smaller samples from across a latitudinal gradient. BST, FYH and induced FFI were first measured in a greenhouse. The seedlings were then planted in the field, where survival and height were measured at the age of 9 years as fitness proxies. We compared between- and within-population variation and genetic correlations of these three seedling traits, and estimated selection gradients at the family level in our main population, taking into account the potential effects of seed weight. Between-population genetic correlations between seedling traits were high (0.76-0.95). Within-population genetic correlations in the main population were lower (0.14-0.35), as in other populations (0.10-0.39). Within population, extensive adaptive variation persists in the seedling traits, in line with rather weak selection gradients, yet maintaining the clines. Although our sampling does not cover the whole cline equally, the results suggest that the individual clines in these traits are maintained by largely independently acting selection, which results in fewer constraints in adaptation under changing climate.

Optical Genome Mapping as an Alternative to FISH-Based Cytogenetic Assessment in Chronic Lymphocytic Leukemia.

The fluorescence in situ hybridization (FISH) technique plays an important role in the risk stratification and clinical management of patients with chronic lymphocytic leukemia (CLL). For genome-wide analysis, FISH needs to be complemented with other cytogenetic methods, including karyotyping and/or chromosomal microarrays. However, this is often not feasible in a diagnostic setup. Optical genome mapping (OGM) is a novel technique for high-resolution genome-wide detection of structural variants (SVs), and previous studies have indicated that OGM could serve as a generic cytogenetic tool for hematological malignancies. Herein, we report the results from our study evaluating the concordance of OGM and standard-of-care FISH in 18 CLL samples. The results were fully concordant between these two techniques in the blinded comparison. Using in silico dilution series, the lowest limit of detection with OGM was determined to range between 3 and 9% variant allele fractions. Genome-wide analysis by OGM revealed additional (>1 Mb) aberrations in 78% of the samples, including both unbalanced and balanced SVs. Importantly, OGM also enabled the detection of clinically relevant complex karyotypes, undetectable by FISH, in three samples. Overall, this study demonstrates the potential of OGM as a first-tier cytogenetic test for CLL and as a powerful tool for genome-wide SV analysis.

Evaluating the role of CHEK2 p.(Asp438Tyr) allele in inherited breast cancer predisposition.

CHEK2 is a well-established breast cancer susceptibility gene. The most frequent pathogenic CHEK2 variant is 1100delC, a loss-of-function mutation conferring 2-fold risk for breast cancer. This gene also harbors other rare variants encountered in the clinical gene panels for hereditary cancer. One of these is CHEK2 c.1312 G > T, p.(Asp438Tyr) in the kinase domain of the protein, but due to its rarity its clinical significance for breast cancer predisposition has remained unclear. Here, we tested the prevalence of CHEK2 p.(Asp438Tyr) allele showing enrichment in the Northern Finnish population, in a total of 2284 breast cancer patients from this geographical region. Genotyping was performed for DNA samples extracted from peripheral blood using high-resolution melt analysis. Fourteen CHEK2 p.(Asp438Tyr) carriers were identified (14/2284, 0.6%, P = 0.67): two in the cohort of breast cancer cases with the indication of inherited disease susceptibility (2/281, 0.7%, P = 1.00) and twelve in the breast cancer cohort unselected for the family history of disease and age at disease onset (12/2003, 0.6%, P = 0.66). This frequency did not differ from the frequency in the general population (10/1299, 0.8%). No CHEK2 p.(Asp438Tyr) homozygotes were identified. Our results indicate that CHEK2 p.(Asp438Tyr) carriers do not have an increased risk for breast cancer and the classification of the CHEK2 p.(Asp438Tyr) variant can be changed from the variant of uncertain significance (VUS) to likely benign for breast cancer.