Re-evaluating the genotypes of patients with adenomatous polyposis of unknown etiology: a nationwide study.

In the Danish Polyposis Register, patients with over 100 cumulative colorectal adenomas of unknown genetic etiology, named in this study colorectal polyposis (CP), is registered and treated as familial adenomatous polyposis (FAP). In this study, we performed genetic analyses, including whole genome sequencing (WGS), of all Danish patients registered with CP and estimated the detection rate of pathogenic variants (PV). We identified 231 families in the Polyposis Register, 31 of which had CP. A polyposis-associated gene panel was performed and, if negative, patients were offered WGS and screening for mosaicism in blood and/or adenomas. Next-generation sequencing (NGS) was carried out for 27 of the families (four declined). PVs were detected in 11 families, and WGS revealed three additional structural variants in APC. Mosaicism of a PV in APC was detected in two families. As the variant detection rate of eligible families was 60%, 93% of families in the register now have a known genetic etiology.

Clinical implications of genetic testing in familial intermediate and late-onset colorectal cancer.

The genetic background of familial, late-onset colorectal cancer (CRC) (i.e., onset > age 50 years) has not been studied as thoroughly as other subgroups of familial CRC, and the proportion of families with a germline genetic predisposition to CRC remains to be defined. To define the contribution of known or suggested CRC predisposition genes to familial late-onset CRC, we analyzed 32 well-established or candidate CRC predisposition genes in 75 families with late-onset CRC. We identified pathogenic or likely pathogenic variants in five patients in MSH6 (n = 1), MUTYH (monoallelic; n = 2) and NTHL1 (monoallelic; n = 2). In addition, we identified a number of variants of unknown significance in particular in the lower penetrant Lynch syndrome-associated mismatch repair (MMR) gene MSH6 (n = 6). In conclusion, screening using a comprehensive cancer gene panel in families with accumulation of late-onset CRC appears not to have a significant clinical value due to the low level of high-risk pathogenic variants detected. Our data suggest that only patients with abnormal MMR immunohistochemistry (IHC) or microsatellite instability (MSI) analyses, suggestive of Lynch syndrome, or a family history indicating another cancer predisposition syndrome should be prioritized for such genetic evaluations. Variants in MSH6 and MUTYH have previously been proposed to be involved in digenic or oligogenic hereditary predisposition to CRC. Accumulation of variants in MSH6 and monoallelic, pathogenic variants in MUTYH in our study indicates that digenic or oligogenic inheritance might be involved in late-onset CRC and warrants further studies of complex types of inheritance.

A Comparison of Tools for Copy-Number Variation Detection in Germline Whole Exome and Whole Genome Sequencing Data.

Copy-number variations (CNVs) have important clinical implications for several diseases and cancers. Relevant CNVs are hard to detect because common structural variations define large parts of the human genome. CNV calling from short-read sequencing would allow single protocol full genomic profiling. We reviewed 50 popular CNV calling tools and included 11 tools for benchmarking in a reference cohort encompassing 39 whole genome sequencing (WGS) samples paired current clinical standard-SNP-array based CNV calling. Additionally, for nine samples we also performed whole exome sequencing (WES), to address the effect of sequencing protocol on CNV calling. Furthermore, we included Gold Standard reference sample NA12878, and tested 12 samples with CNVs confirmed by multiplex ligation-dependent probe amplification (MLPA). Tool performance varied greatly in the number of called CNVs and bias for CNV lengths. Some tools had near-perfect recall of CNVs from arrays for some samples, but poor precision. Several tools had better performance for NA12878, which could be a result of overfitting. We suggest combining the best tools also based on different methodologies: GATK gCNV, Lumpy, DELLY, and cn.MOPS. Reducing the total number of called variants could potentially be assisted by the use of background panels for filtering of frequently called variants.

A rare missense variant in APC interrupts splicing and causes AFAP in two Danish families.

BACKGROUND: We report the first case of a missense variant in the APC gene that interrupts splicing by creating a new cryptic acceptor site. The variant, c.289G>A, p.(Gly97Arg), is located in exon 3, and qualitative and semi-quantitative RNA splicing analysis reveal that the variant results in skipping of the last 70 nucleotides of the exon, which leads to the introduction of a frameshift and a premature stop codon. CASE PRESENTATION: The variant was detected in two, apparently unrelated, Danish families with an accumulation of colorectal cancers, colonic adenomas and other cancers. The families both have an attenuated familial adenomatous polyposis phenotype, which is consistent with the association of pathogenic variants in the 5' end of the gene.One variant-carrier also had Caroli Disease and a Caroli Disease associated hepatic mucinous cystadenocarcinoma. This is the first description of a person with both Caroli Disease and a pathogenic APC variant, and although the APC variant is not known to be connected to the development of the hepatic malformations in Caroli Disease, it remains unclear whether the variant could have contributed to the carcinogenesis of the liver tumour. CONCLUSIONS: Based on functional and co-segregation data we classify the APC c.289G>A, p.(Gly97Arg) variant as pathogenic (class 5). Our findings emphasize the importance of a functional evaluation of missense variants although located far from the exon-intron boundaries.

Amplicon-Based NGS Panels for Actionable Cancer Target Identification in Follicular Cell-Derived Thyroid Neoplasia.

Follicular cell-derived thyroid cancers are heterogenous and morphological classification is a complex and highly specialized task. Hence, identification of somatic alterations could provide insights to tumor biology and serve as an add-on diagnostic tool. Furthermore, results from these add-on tools could point in the direction of a more personalized treatment strategy. In the present study we set out to identify and validate the somatic mutation profile in a sample-set of follicular cell-derived thyroid neoplasia. One-hundred-and-one archived formalin fixed paraffin embedded (FFPE) tissue samples from patients diagnosed with follicular cell-derived thyroid neoplasia were included, and upon DNA-extraction and qualitative measurements 99 samples were eligible for amplicon-based next-generation-sequencing. Libraries were generated using the TruSeq Amplicon Cancer Panel, followed by sequencing using a MiSeq. Upon data processing and variant filtering all variants were manually assessed to exclude false positive mutations in the final curated list. Moreover, hot-spot mutations were validated using an independent platform from Agilent. Each diagnostic group were correlated to mutation burden and individual mutations were classified according to recent guidelines for somatic mutation classification. Close to 100% of the archived FFPE samples were eligible for DNA-library preparation and amplicon sequencing based on DNA quality criterion. The distribution of mutations in the specific diagnostic groups resulted in a higher mutation frequency among the most dedifferentiated than in the groups with a more differentiated cell profile. Based on the distribution mutations across the samples and using hierarchical clustering, we generated four tentative mutational signatures; highly mutated tumors; tumors with mainly NRAS and TP53 mutations; BRAF mutated tumors and tumors with none or single sporadic mutations. Future studies including more samples and follow-up data may amend these signatures, however our results imply that morphological classification of follicular cell derived thyroid neoplasia could be supplemented with a somatic mutational signature. Taken together, broad screening of the somatic alterations in FFPE tissue of thyroid neoplasia is comprehensible and essential for future identification of possible treatment targets and personalized medicine.

WISH-R- a fast and efficient tool for construction of epistatic networks for complex traits and diseases.

BACKGROUND: Genetic epistasis is an often-overlooked area in the study of the genomics of complex traits. Genome-wide association studies are a useful tool for revealing potential causal genetic variants, but in this context, epistasis is generally ignored. Data complexity and interpretation issues make it difficult to process and interpret epistasis. As the number of interaction grows exponentially with the number of variants, computational limitation is a bottleneck. Gene Network based strategies have been successful in integrating biological data and identifying relevant hub genes and pathways related to complex traits. In this study, epistatic interactions and network-based analysis are combined in the Weighted Interaction SNP hub (WISH) method and implemented in an efficient and easy to use R package. RESULTS: The WISH R package (WISH-R) was developed to calculate epistatic interactions on a genome-wide level based on genomic data. It is easy to use and install, and works on regular genomic data. The package filters data based on linkage disequilibrium and calculates epistatic interaction coefficients between SNP pairs based on a parallelized efficient linear model and generalized linear model implementations. Normalized epistatic coefficients are analyzed in a network framework, alleviating multiple testing issues and integrating biological signal to identify modules and pathways related to complex traits. Functions for visualizing results and testing runtimes are also provided. CONCLUSION: The WISH-R package is an efficient implementation for analyzing genome-wide epistasis for complex diseases and traits. It includes methods and strategies for analyzing epistasis from initial data filtering until final data interpretation. WISH offers a new way to analyze genomic data by combining epistasis and network based analysis in one method and provides options for visualizations. This alleviates many of the existing hurdles in the analysis of genomic interactions.

Carboxylesterase 1 genes: systematic review and evaluation of existing genotyping procedures.

The carboxylesterase 1 gene (CES1) encodes a hydrolase that metabolizes commonly used drugs. The CES1-related pseudogene, carboxylesterase 1 pseudogene 1 (CES1P1), has been implicated in gene exchange with CES1 and in the formation of hybrid genes including the carboxylesterase 1A2 gene (CES1A2). Hence, the CES1 region is complex. Using in silico PCR and alignment, we assessed the specificity of PCR-assisted procedures for genotyping CES1, CES1A2 and CES1P1 in studies identified in PubMed. We identified 33 such studies and excluded those that were not the first to use a procedure or lacked sequence information. After this 17 studies remained. Ten of these used haplotype-specific amplification, restriction enzyme treatment or amplicon sequencing, and included five that were predicted to lack specificity. All procedures for genotyping of single nucleotide polymorphisms in eight studies lacked specificity. One of these studies also used amplicon sequencing, thus being present in the group above. Some primers and their intended targets were mismatched. We provide experimental evidence that one of the procedures lacked specificity. Additionally, a complex pattern of segmental duplications in the CES1 region was revealed. In conclusion, many procedures for CES1, CES1A2 and CES1P1 genotyping appear to lack specificity. Knowledge about the segmental duplications may improve the typing of these genes.

Carboxylesterase 1A2 encoding gene with increased transcription and potential rapid drug metabolism in Asian populations.

The carboxylesterase 1 gene (CES1) encodes a hydrolase implicated in the metabolism of commonly used drugs. CES1A2, a hybrid of CES1 and a CES1-like pseudogene, has a promoter that is weak in most individuals. However, some individuals harbor a promoter haplotype of this gene with two overlapping Sp1 sites that confer significantly increased transcription potentially leading to rapid drug metabolism. This CES1A2 haplotype has previously been reported to be common among Asians. Using polymerase chain reaction followed by sequencing, the present study examined variation in the promoter and 5' untranslated region of CES1A2 in 120 Han Chinese and 120 Japanese people enrolled in the 1000 Genomes Project. We identified 11 single nucleotide variations, two of which were novel, in 145 of the individuals who were found to carry CES1A2. Alignment analysis indicated that the CES1A2 haplotype with the overlapping Sp1 sites has been generated by incorporation of a segment of CES1. All minor allele frequencies were equal to or below 0.022 and the frequencies of the minor haplotypes were up to 40-fold lower than previously reported, including that of the haplotype with the overlapping Sp1 sites. This information is novel and suggests that the pharmacogenetic relevance of CES1A2 is limited in Asians.

Reappraisal of the genetic diversity and pharmacogenetic assessment of CES1.

The CES1 gene encodes a hydrolase that metabolizes important drugs. Variants generated by exchange of segments with CES1P1 complicate genotyping of CES1. Using a highly specific procedure we examined DNA samples from 200 Caucasians and identified 46 single nucleotide variants (SNVs) in CES1 and 21 SNVs in CES1A2, a hybrid composed of CES1 and CES1P1. Several of these SNVs were novel. The frequencies of SNVs with a potential functional impact were below 0.02 suggesting limited pharmacogenetic potential for CES1 genotyping. In silico PCR revealed that the majority of the primer pairs for amplification of CES1 or CES1A2 in three previous studies lacked specificity, which partially explains a limited overlap with our findings.

Identification of 2 loci associated with development of myxomatous mitral valve disease in Cavalier King Charles Spaniels.

Myxomatous mitral valve disease (MMVD) is the most common heart disease in dogs. It is characterized by chronic progressive degenerative lesions of the mitral valve. The valve leaflets become thickened and prolapse into the left atrium resulting in mitral regurgitation (MR). MMVD is most prevalent in small to medium sized dog breeds, Cavalier King Charles Spaniels (CKCS) in particular. The onset of MMVD is highly age dependent, and at the age of 10 years, nearly all CKCS are affected. The incidence of a similar disease in humans-mitral valve prolapse-is 1-5%. By defining CKCSs with an early onset of MMVD as cases and old dogs with no or mild signs of MMVD as controls, we conducted a genome-wide association study (GWAS) to identify loci associated with development of MMVD. We have identified a 1.58 Mb region on CFA13 (P(genome) = 4.0 × 10(-5)) and a 1.68 Mb region on CFA14 (P(genome) = 7.9 × 10(-4)) associated with development of MMVD. This confirms the power of using the dog as a model to uncover potential candidate regions involved in the molecular mechanisms behind complex traits.

Expression studies of the obesity candidate gene FTO in pig.

Obesity is an increasing problem worldwide and research on candidate genes in good animal models is highly needed. The pig is an excellent model as its metabolism, organ size, and eating habits resemble that of humans. The present study is focused on the characterization of the fat mass and obesity associated gene (FTO) in pig. This gene has recently been associated with increased body mass index in several human populations. To establish information on the expression profile of FTO in the pig we performed quantitative PCR in a panel of adult pig tissues and in tissues sampled at different developmental stages. Expression of the FTO transcript was detected in all tissues tested with significantly higher levels in brain tissues (cortex, cerebellum, and hippocampus; P < 0.001). These levels varied through the development and between the specific parts of the brain studied (i.e., frontal cortex and cerebellum). Additionally, in order to see the involvement of the FTO gene in obesity, the changes in expression level were investigated in a nutritional study in brain of Gottingen minipigs under a high cholesterol diet. Significantly higher (P < 0.01) levels of FTO transcript were found in cortex when compared with cerebellum of the high-cholesterol fed pigs. Furthermore, SNPs were investigated in the coding sequence of the FTO in the Gottingen minipig and in the Danish commercial pig. Eleven synonymous SNPs and a two bp insertion were found between the two pig lines.