Search

Germline pathogenic variants in the DNA mismatch (MMR) repair genes (Lynch syndrome) predispose to colorectal (CRC) and endometrial (EC) cancer. However, mosaic variants in the MMR genes have been rarely described. We identified a likely de novo mosaic MSH6:c.1135_1139del p.Arg379* pathogenic variant in a patient diagnosed with suspected Lynch syndrome/Lynch-like syndrome. The patient developed MSH6-deficient EC and CRC at 54 and 58 years of age, respectively, without a detectable germline MMR pathogenic variant. Multigene panel sequencing of tumor and blood-derived DNA identified an MSH6 somatic mutation (MSH6:c.1135_1139del p.Arg379*) common to both the EC and CRC, raising suspicion of mosaicism. A droplet digital polymerase chain reaction (ddPCR) assay detected the MSH6 variant at 5.34% frequency in normal colonic tissue, 3.49% in saliva and 1.64% in blood DNA, demonstrating the presence of the MSH6 variant in all three germ layers. This study highlights the utility of tumor sequencing to guide sensitive ddPCR testing to detect low-level mosaicism in the MMR genes. Further investigation of the prevalence of MMR mosaicism is needed to inform routine diagnostic approaches and genetic counselling.

Background and AimsThe microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain ofEscherichia coliharboring thepksisland that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown.MethodsSBS88-positive CRCs were identified from targeted sequencing data from 5,292 CRCs from 17 studies and tested for their association with clinico-pathological features, oncogenic pathways, genomic characteristics and survival.ResultsIn total, 7.5% (398/5,292) of the CRCs were SBS88-positive, of which 98.7% (392/398) were microsatellite stable/microsatellite instability low (MSS/MSI-L), compared with 80% (3916/4894) of SBS88 negative tumors (p=1.5×10-28). Analysis of MSS/MSI-L CRCs demonstrated that SBS88 positive CRCs were associated with the distal colon (OR=1.84, 95% CI=1.40-2.42, p=1×10-5) and rectum (OR=1.90, 95% CI=1.44-2.51, p=6×10-6) tumor sites compared with the proximal colon. The top seven recurrent somatic mutations associated with SBS88-positive CRCs demonstrated mutational contexts associated with colibactin-induced DNA damage, the strongest of which was theAPC:c.835-8A>G mutation (OR=65.5, 95%CI=39.0-110.0, p=3×10-80). Large copy number alterations (CNAs) including CNA loss on 14q and gains on 13q, 16q and 20p were significantly enriched in SBS88- positive CRCs. SBS88-positive CRCs were associated with better CRC-specific survival (p=0.007; hazard ratio of 0.69, 95% CI=0.52-0.90) when stratified by age, sex, study, and by stage.ConclusionSBS88-positivity, a biomarker of colibactin-induced DNA damage, can identify a novel subtype of CRC characterized by recurrent somatic mutations, copy number alterations and better survival. These findings provide new insights for treatment and prevention strategies for this subtype of CRC.

Germline loss-of-function variants in AXIN2 are associated with oligodontia and ectodermal dysplasia. The association between colorectal cancer (CRC) and colonic polyposis is less clear despite this gene now being included in multi-gene panels for CRC. Study participants were people with genetically unexplained colonic polyposis recruited to the Genetics of Colonic Polyposis Study who had a rare germline AXIN2 gene variant identified from either clinical multi-gene panel testing (n=2) or from whole genome/exome sequencing (n=2). Variant segregation in relatives and characterisation of tumour tissue were performed where possible. Four different germline pathogenic variants in AXIN2 were identified in four families. Five of the seven carriers of the c.1049delC, p.Pro350Leufs*13 variant, two of the six carriers of the c.1994dupG, p.Asn666Glnfs*41 variant, all three carriers of c.1972delA, p.Ser658Alafs*31 variant and the single proband carrier of the c.2405G>C, p.Arg802Thr variant, which creates an alternate splice form resulting in a frameshift mutation (p.Glu763Ilefs*42), were affected by CRC and/or polyposis. Carriers had a mean age at diagnosis of CRC/polyposis of 52.5 ± 9.2 years. Colonic polyps were typically pan colonic with counts ranging from 5 to >100 (median 12.5) comprising predominantly adenomatous polyps but also serrated polyps. Two CRCs from carriers displayed evidence of a second hit via loss of heterozygosity. Oligodontia was observed in carriers from two families. Germline AXIN2 pathogenic variants from four families were associated with CRC and/or polyposis in multiple family members. These findings support the inclusion of AXIN2 in CRC and polyposis multigene panels for clinical testing.

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri. We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most “missing” orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Background Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease. Methods We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables. Results Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control). Conclusions In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

Identifying tumor DNA mismatch repair deficiency (dMMR) is important for precision medicine. We assessed tumor features, individually and in combination, in whole-exome sequenced (WES) colorectal cancers (CRCs) and in panel sequenced CRCs, endometrial cancers (ECs) and sebaceous skin tumors (SSTs) for their accuracy in detecting dMMR. CRCs (n=300) with WES, where MMR status was determined by immunohistochemistry, were assessed for microsatellite instability (MSMuTect, MANTIS, MSIseq, MSISensor), COSMIC tumor mutational signatures (TMS) and somatic mutation counts. A 10-fold cross-validation approach (100 repeats) evaluated the dMMR prediction accuracy for 1) individual features, 2) Lasso statistical model and 3) an additive feature combination approach. Panel sequenced tumors (29 CRCs, 22 ECs, 20 SSTs) were assessed for the top performing dMMR predicting features/models using these three approaches. For WES CRCs, 10 features provided >80% dMMR prediction accuracy, with MSMuTect, MSIseq, and MANTIS achieving ≥99% accuracy. The Lasso model achieved 98.3%. The additive feature approach with ≥3/6 of MSMuTect, MANTIS, MSIseq, MSISensor, INDEL count or TMS ID2+ID7 achieved 99.7% accuracy. For the panel sequenced tumors, the additive feature combination approach of ≥3/6 achieved accuracies of 100%, 95.5% and 100%, for CRCs, ECs, and SSTs, respectively. The microsatellite instability calling tools performed well in WES CRCs, however, an approach combining tumor features may improve dMMR prediction in both WES and panel sequenced data across tissue types.

Patients in whom mismatch repair (MMR)-deficient cancer develops in the absence of pathogenic variants of germline MMR genes or somatic hypermethylation of the MLH1 gene promoter are classified as having suspected Lynch syndrome (SLS). Germline whole-genome sequencing (WGS) and targeted and genome-wide tumor sequencing were applied to identify the underlying cause of tumor MMR deficiency in SLS. Germline WGS was performed on samples from 14 cancer-affected patients with SLS, including two sets of first-degree relatives. MMR genes were assessed for germline pathogenic variants, including complex structural rearrangements and noncoding variants. Tumor tissue was assessed for somatic MMR gene mutations using targeted, whole-exome sequencing or WGS. Germline WGS identified pathogenic MMR variants in 3 of the 14 cases (21.4%), including a 9.5-megabase inversion disrupting MSH2 in a mother and daughter. Excluding these 3 MMR carriers, tumor sequencing identified at least two somatic MMR gene mutations in 8 of 11 tumors tested (72.7%). In a second mother–daughter pair, a somatic cause of tumor MMR deficiency was supported by the presence of double somatic MSH2 mutations in their respective tumors. More than 70% of SLS cases had double somatic MMR mutations in the absence of germline pathogenic variants in the MMR or other DNA repair–related genes on WGS, and, therefore, were confidently assigned a noninherited cause of tumor MMR deficiency.

Inherited genetic variation profoundly influences cancer risk and outcome. While the impact of germline single nucleotide polymorphisms has been well-studied in several cancer types, the effects of germline structural variants (gSVs) on cancer biology and clinical outcomes is largely unknown. From our cohort of 300 men with localized, intermediate risk prostate cancer, we identified 6,003 gSVs present in at least 3% of patients; 48 were associated with recurrent somatic alterations or clinical outcome. Of these, approximately 50% were associated with expression of nearby genes or intersected with exons or regulatory regions. Using external cohorts, we validated three gSVs that were strongly associated with poor clinical outcomes, including an inversion at chr14q24.1 present in ~20% of patients. Notably, a strong synergistic effect on outcome was observed in patients with somatic TP53 alterations or high genomic instability, defining a new aggressive prostate cancer subtype with chr14INV as a novel, recurrent biomarker. Citation Format: Nicholas K. Wang, Alexandre Rouette, Kathleen E. Houlahan, Takafumi N. Yamaguchi, Julie Livingstone, Chol-Hee Jung, Peter Georgeson, Michael Fraser, Yu-Jia Shiah, Cindy Q. Yao, Vincent Huang, Natalie S. Fox, Natalie Kurganovs, Katayoon Kasaian, Veronica Y. Sabelnykova, Jay Jayalath, Kenneth Weke, Helen Zhu, Theodorus van der Kwast, Tony Papenfuss, Housheng H. He, Niall M. Corcoran, Robert G. Bristow, Alexandre R. Zlotta, Christopher Hovens, Paul C. Boutros. Germline structural variants shape prostate cancer clinical and molecular evolution. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4305.

Carriers of germline biallelic pathogenic variants in the MUTYH gene have a high risk of colorectal cancer. We test 5649 colorectal cancers to evaluate the discriminatory potential of a tumor mutational signature specific to MUTYH for identifying biallelic carriers and classifying variants of uncertain clinical significance (VUS). Using a tumor and matched germline targeted multi-gene panel approach, our classifier identifies all biallelic MUTYH carriers and all known non-carriers in an independent test set of 3019 colorectal cancers (accuracy = 100% (95% confidence interval 99.87–100%)). All monoallelic MUTYH carriers are classified with the non-MUTYH carriers. The classifier provides evidence for a pathogenic classification for two VUS and a benign classification for five VUS. Somatic hotspot mutations KRAS p.G12C and PIK3CA p.Q546K are associated with colorectal cancers from biallelic MUTYH carriers compared with non-carriers (p = 2 × 10−23 and p = 6 × 10−11, respectively). Here, we demonstrate the potential application of mutational signatures to tumor sequencing workflows to improve the identification of biallelic MUTYH carriers.

We investigated aberrant DNA methylation (DNAm) changes and the contribution of ageing-associated methylomic drift and age acceleration to early-onset colorectal cancer (EOCRC) carcinogenesis. Genome-wide DNAm profiling using the Infinium HM450K on 97 EOCRC tumour and 54 normal colonic mucosa samples was compared with: (1) intermediate-onset CRC (IOCRC, diagnosed between 50–70 years, 343 tumour and 35 normal), and (2) late-onset CRC (LOCRC, &gt, 70 years, 318 tumour and 40 normal). CpGs associated with age-related methylation drift were identified using a public dataset of 231 normal mucosa samples from people without CRC. DNAm-age was estimated using epiTOC2. Common to all three age-of-onset groups, 88,385 (20% of all CpGs) CpGs were differentially methylated between tumour and normal mucosa. We identified 234 differentially methylated genes that were unique to the EOCRC group, 13 of these DMRs/genes were replicated in EOCRC compared with LOCRCs from TCGA. In normal mucosa from people without CRC, we identified 28,154 CpGs that undergo ageing-related DNAm drift, and of those, 65% were aberrantly methylated in EOCRC tumours. Based on the mitotic-based DNAm clock epiTOC2, we identified age acceleration in normal mucosa of people with EOCRC compared with normal mucosa from the IOCRC, LOCRC groups (p = 3.7 × 10−16) and young people without CRC (p = 5.8 × 10−6). EOCRC acquires unique DNAm alterations at 234 loci. CpGs associated with ageing-associated drift were widely affected in EOCRC without needing the decades-long accrual of DNAm drift as commonly seen in intermediate- and late-onset CRCs. Accelerated ageing in normal mucosa from people with EOCRC potentially underlies the earlier age of diagnosis in CRC carcinogenesis.

BackgroundThe role of DNA methylation (DNAm) in the carcinogenesis of colorectal cancer (CRC) diagnosed MethodsGenome-wide DNAm profiling using the Infinium HM450K on 97 EOCRC tumour and 54 normal colonic mucosa samples was compared with: 1) intermediate-onset CRC (IOCRC; diagnosed between 50-70 years; 343 tumour and 35 normal); and 2) late-onset CRC (LOCRC; >70 years; 318 tumour and 40 normal). CpGs associated with age-related methylation drift were identified using a public dataset of 231 normal mucosa samples from people without CRC. DNAm-age was estimated using epiTOC2.ResultsCommon to all three age-of-onset groups, 88,385 (20% of all CpGs) CpGs were differentially methylated between tumour and normal mucosa. We identified 234 differentially methylated genes that were unique to the EOCRC group. In normal mucosa from people without CRC, we identified 28,154 CpGs that undergo ageing-related DNAm drift and of those, 65% were aberrantly methylated in EOCRC tumours. Based on the mitotic-based DNAm clock epiTOC2, we identified age acceleration in normal mucosa of people with EOCRC compared with normal mucosa from the IOCRC, LOCRC groups (p=3.7×10−16) and young people without CRC (p=5.8×10−6).ConclusionEOCRC acquires unique DNAm alterations at 234 loci. CpGs associated with ageing-associated drift were widely affected in EOCRC without needing the decades-long accrual of DNAm drift as commonly seen in intermediate- and late-onset CRCs. We found accelerated ageing in normal mucosa from people with EOCRC, as evidenced by a faster stem-cell division rate, potentially contributing to EOCRC carcinogenesis.

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri. We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most “missing” orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Contains fulltext : 228713.pdf (Publisher’s version ) (Open Access)

Recent publications have shown patients with defects in the DNA mismatch repair (MMR) pathway driven by either MSH2 or MSH6 loss experience a significant increase in the incidence of prostate cancer. Moreover, this increased incidence of prostate cancer is accompanied by rapid disease progression and poor clinical outcomes. We show that androgen-receptor activation, a key driver of prostate carcinogenesis, can disrupt the MSH2 gene in prostate cancer. We screened tumours from two cohorts (recurrent/non-recurrent) of prostate cancer patients to confirm the loss of MSH2 protein expression and identified decreased MSH2 expression in recurrent cases. Stratifying the independent TCGA prostate cancer cohort for MSH2/6 expression revealed that patients with lower levels of MSH2/6 had significant worse outcomes, in contrast, endometrial and colorectal cancer patients with lower MSH2/6 levels. MMRd endometrial and colorectal tumours showed the expected increase in mutational burden, microsatellite instability and enhanced immune cell mobilisation but this was not evident in prostate tumours. We have shown that loss or reduced levels of MSH2/MSH6 protein in prostate cancer is associated with poor outcome. However, our data indicate that this is not associated with a statistically significant increase in mutational burden, microsatellite instability or immune cell mobilisation in a cohort of primary prostate cancers.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21–0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21–1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features., Large scale sequencing study is of paramount importance to unravel the heterogeneity of colorectal cancer. Here, the authors sequenced 205 cancer genes in more than 2000 tumours and identified additional mutated driver genes, determined that mutational burden and specific mutations in TP53 are associated with survival odds.

ObjectiveGermline pathogenic variants (PVs) in the DNA mismatch repair (MMR) genes and in the base excision repair gene MUTYH underlie hereditary colorectal cancer (CRC) and polyposis syndromes. We evaluated the robustness and discriminatory potential of tumour mutational signatures in CRCs for identifying germline PV carriers.DesignWhole exome sequencing of formalin-fixed paraffin embedded (FFPE) CRC tissue was performed on 33 MMR germline PV carriers, 12 biallelic MUTYH germline PV carriers, 25 sporadic MLH1 methylated MMR-deficient CRCs (MMRd controls) and 160 sporadic MMR-proficient CRCs (MMRp controls) and included 498 TCGA CRC tumours. COSMIC V3 single base substitution (SBS) and indel (ID) mutational signatures were assessed for their ability to differentiate CRCs that developed in carriers from non-carriers.ResultsThe combination of mutational signatures SBS18 and SBS36 contributing >30% of a CRC’s signature profile was able to discriminate biallelic MUTYH carriers from all other non-carrier control CRCs with 100% accuracy (area under the curve (AUC) 1.0). SBS18 and SBS36 were associated with specific MUTYH variants p.Gly396Asp (p=0.025) and p.Tyr179Cys (p=5×10−5), respectively. The combination of ID2 and ID7 could discriminate the 33 MMR PV carrier CRCs from the MMRp control CRCs (AUC 0.99), however, SBS and ID signatures, alone or in combination, could not provide complete discrimination (AUC 0.79) between CRCs from MMR PV carriers and sporadic MMRd controls.ConclusionAssessment of SBS and ID signatures can discriminate CRCs from biallelic MUTYH carriers and MMR PV carriers from non-carriers with high accuracy, demonstrating utility as a potential diagnostic and variant classification tool.SIGNIFICANCE OF THE STUDYWhat is already known about this subject?Identifying carriers of pathogenic variants (PVs) in moderate/high-risk colorectal cancer (CRC) and polyposis susceptibility genes has clinical relevance for diagnosis, targeted screening and prevention strategies, prognosis, and treatment options. However, challenges still remain in the identification of carriers and the classification of rare variants in these genes.Previous studies have identified tumour mutational signatures that result from defective DNA repair including DNA mismatch repair (MMR) deficiency and base excision repair defects, DNA repair mechanisms that underlie the common hereditary CRC and polyposis syndromes but their diagnostic utility in CRC is unknown.What are the new findings?Single base substitution (SBS)-related mutational signatures derived from whole exome sequencing of formalin-fixed paraffin embedded (FFPE)-derived CRC tissue DNA can effectively discriminate CRCs that developed in biallelic MUTYH PV carriers from CRC-affected non-carriers.CRCs that develop in MMR PV carriers (Lynch syndrome) can be effectively differentiated from sporadic MMR-proficient CRC by a combination of indel (ID) signatures, but the SBS and ID tumour mutational signatures are less effective at discriminating Lynch syndrome-related CRC from sporadic MMR-deficient CRC resulting from MLH1 gene promoter hypermethylation.The SBS and ID mutational signatures associated with biallelic MUTYH PV carriers and MMR PV carriers are robust to changes in experimental settings.We demonstrate the optimal experimental settings for calculating mutational signatures and define thresholds that optimise sensitivity and specificity for classifying CRC aetiology as hereditary or non-hereditary.How might it impact on clinical practice in the foreseeable future?Deriving SBS- and ID-related mutational signatures from CRCs can identify carriers of PVs in hereditary CRC and polyposis susceptibility genes.The application of mutational signatures has the potential to improve the diagnosis of hereditary CRC and aid in variant classification, leading to improved clinical management and CRC prevention.

Background Bioinformatics software tools are often created ad hoc, frequently by people without extensive training in software development. In particular, for beginners, the barrier to entry in bioinformatics software development is high, especially if they want to adopt good programming practices. Even experienced developers do not always follow best practices. This results in the proliferation of poorer-quality bioinformatics software, leading to limited scalability and inefficient use of resources; lack of reproducibility, usability, adaptability, and interoperability; and erroneous or inaccurate results. Findings We have developed Bionitio, a tool that automates the process of starting new bioinformatics software projects following recommended best practices. With a single command, the user can create a new well-structured project in 1 of 12 programming languages. The resulting software is functional, carrying out a prototypical bioinformatics task, and thus serves as both a working example and a template for building new tools. Key features include command-line argument parsing, error handling, progress logging, defined exit status values, a test suite, a version number, standardized building and packaging, user documentation, code documentation, a standard open source software license, software revision control, and containerization. Conclusions Bionitio serves as a learning aid for beginner-to-intermediate bioinformatics programmers and provides an excellent starting point for new projects. This helps developers adopt good programming practices from the beginning of a project and encourages high-quality tools to be developed more rapidly. This also benefits users because tools are more easily installed and consistent in their usage. Bionitio is released as open source software under the MIT License and is available at https://github.com/bionitio-team/bionitio.

Background Muir‐Torre syndrome is defined by the development of sebaceous skin lesions in individuals who carry a germline mismatch repair (MMR) gene mutation. Loss of expression of MMR proteins is frequently observed in sebaceous skin lesions, but MMR‐deficiency alone is not diagnostic for carrying a germline MMR gene mutation. Methods Whole exome sequencing was performed on three MMR‐deficient sebaceous lesions from individuals with MSH2 gene mutations (Lynch syndrome) and three MMR‐proficient sebaceous lesions from individuals without Lynch syndrome with the aim of characterizing the tumor mutational signatures, somatic mutation burden, and microsatellite instability status. Thirty predefined somatic mutational signatures were calculated for each lesion. Results Signature 1 was ubiquitous across the six lesions tested. Signatures 6 and 15, associated with defective DNA MMR, were significantly more prevalent in the MMR‐deficient lesions from the MSH2 carriers compared with the MMR‐proficient non‐Lynch sebaceous lesions (mean ± SD=41.0 ± 8.2% vs. 2.3 ± 4.0%, p = 0.0018). Tumor mutation burden was, on average, significantly higher in the MMR‐deficient lesions compared with the MMR‐proficient lesions (23.3 ± 11.4 vs. 1.8 ± 0.8 mutations/Mb, p = 0.03). All four sebaceous lesions observed in sun exposed areas of the body demonstrated signature 7 related to ultraviolet light exposure. Conclusion Tumor mutational signatures 6 and 15 and somatic mutation burden were effective in differentiating Lynch‐related from non‐Lynch sebaceous lesions.

DNA methylation influences predisposition, development and prognosis for many diseases, including cancer. However, it is not uncommon to encounter samples with incorrect sex labelling or atypical sex chromosome arrangement. Sex is one of the strongest influencers of the genomic distribution of DNA methylation and, therefore, correct assignment of sex and filtering of abnormal samples are essential for the quality control of study data. Differences in sex chromosome copy numbers between sexes and X-chromosome inactivation in females result in distinctive sex-specific patterns in the distribution of DNA methylation levels. In this study, we present a software tool, sEst, which incorporates clustering analysis to infer sex and to detect sex-chromosome abnormalities from DNA methylation microarray data. Testing with two publicly available datasets demonstrated that sEst not only correctly inferred the sex of the test samples, but also identified mislabelled samples and samples with potential sex-chromosome abnormalities, such as Klinefelter syndrome and Turner syndrome, the latter being a feature not offered by existing methods. Considering that sex and the sex-chromosome abnormalities can have large effects on many phenotypes, including diseases, our method can make a significant contribution to DNA methylation studies that are based on microarray platforms.

Background Large-scale computational prediction of protein structures represents a cost-effective alternative to empirical structure determination with particular promise for non-model organisms and neglected pathogens. Conventional sequence-based tools are insufficient to annotate the genomes of such divergent biological systems. Conversely, protein structure tolerates substantial variation in primary amino acid sequence and is thus a robust indicator of biochemical function. Structural proteomics is poised to become a standard part of pathogen genomics research; however, informatic methods are now required to assign confidence in large volumes of predicted structures. Aims Our aim was to predict the proteome of a neglected human pathogen, Giardia duodenalis, and stratify predicted structures into high- and lower-confidence categories using a variety of metrics in isolation and combination. Methods We used the I-TASSER suite to predict structural models for ∼5,000 proteins encoded in G. duodenalis and identify their closest empirically-determined structural homologues in the Protein Data Bank. Models were assigned to high- or lower-confidence categories depending on the presence of matching protein family (Pfam) domains in query and reference peptides. Metrics output from the suite and derived metrics were assessed for their ability to predict the high-confidence category individually, and in combination through development of a random forest classifier. Results We identified 1,095 high-confidence models including 212 hypothetical proteins. Amino acid identity between query and reference peptides was the greatest individual predictor of high-confidence status; however, the random forest classifier outperformed any metric in isolation (area under the receiver operating characteristic curve = 0.976) and identified a subset of 305 high-confidence-like models, corresponding to false-positive predictions. High-confidence models exhibited greater transcriptional abundance, and the classifier generalized across species, indicating the broad utility of this approach for automatically stratifying predicted structures. Additional structure-based clustering was used to cross-check confidence predictions in an expanded family of Nek kinases. Several high-confidence-like proteins yielded substantial new insight into mechanisms of redox balance in G. duodenalis-a system central to the efficacy of limited anti-giardial drugs. Conclusion Structural proteomics combined with machine learning can aid genome annotation for genetically divergent organisms, including human pathogens, and stratify predicted structures to promote efficient allocation of limited resources for experimental investigation.

Constitutional biological processes involve the generation of DNA double-strand breaks (DSBs). The production of such breaks and their subsequent resolution are also highly relevant to neurodegenerative diseases and cancer, in which extensive DNA fragmentation has been described Stephens et al. (2011), Blondet et al. (2001). Tchurikov et al. Tchurikov et al. (2011, 2013) have reported previously that frequent sites of DSBs occur in chromosomal domains involved in the co-ordinated expression of genes. This group report that hot spots of DSBs in human HEK293T cells often coincide with H3K4me3 marks, associated with active transcription Kravatsky et al. (2015) and that frequent sites of DNA double-strand breakage are likely to be relevant to cancer genomics Tchurikov et al. (2013, 2016) . Recently, they applied a RAFT (rapid amplification of forum termini) protocol that selects for blunt-ended DSB sites and mapped these to the human genome within defined co-ordinate ‘windows’. In this paper, we re-analyse public RAFT data to derive sites of DSBs at the single-nucleotide level across the built genome for human HEK293T cells (https://figshare.com/s/35220b2b79eaaaf64ed8). This refined mapping, combined with accessory ENCODE data tracks and ribosomal DNA-related sequence annotations, will likely be of value for the design of clinically relevant targeted assays such as those for cancer susceptibility, diagnosis, treatment-matching and prognostication.

Importance Optimal use of whole-exome sequencing (WES) in the pediatric setting requires an understanding of who should be considered for testing and when it should be performed to maximize clinical utility and cost-effectiveness. Objectives To investigate the impact of WES in sequencing-naive children suspected of having a monogenic disorder and evaluate its cost-effectiveness if WES had been available at different time points in their diagnostic trajectory. Design, Setting, and Participants This prospective study was part of the Melbourne Genomics Health Alliance demonstration project. At the ambulatory outpatient clinics of the Victorian Clinical Genetics Services at the Royal Children’s Hospital, Melbourne, Australia, children older than 2 years suspected of having a monogenic disorder were prospectively recruited from May 1 through November 30, 2015, by clinical geneticists after referral from general and subspecialist pediatricians. All children had nondiagnostic microarrays and no prior single-gene or panel sequencing. Exposures All children underwent singleton WES with targeted phenotype-driven analysis. Main Outcomes and Measures The study examined the clinical utility of a molecular diagnosis and the cost-effectiveness of alternative diagnostic trajectories, depending on timing of WES. Results Of 61 children originally assessed, 44 (21 [48%] male and 23 [52%] female) aged 2 to 18 years (mean age at initial presentation, 28 months; range, 0-121 months) were recruited, and a diagnosis was achieved in 23 (52%) by singleton WES. The diagnoses were unexpected in 8 of 23 (35%), and clinical management was altered in 6 of 23 (26%). The mean duration of the diagnostic odyssey was 6 years, with each child having a mean of 19 tests and 4 clinical genetics and 4 nongenetics specialist consultations, and 26 (59%) underwent a procedure while under general anesthetic for diagnostic purposes. Economic analyses of the diagnostic trajectory identified that WES performed at initial tertiary presentation resulted in an incremental cost savings of A$9020 (US$6838) per additional diagnosis (95% CI, A$4304-A$15 404 [US$3263-US$11 678]) compared with the standard diagnostic pathway. Even if WES were performed at the first genetics appointment, there would be an incremental cost savings of A$5461 (US$4140) (95% CI, A$1433-A$10 557 [US$1086- US$8004]) per additional diagnosis compared with the standard diagnostic pathway. Conclusions and Relevance Singleton WES in children with suspected monogenic conditions has high diagnostic yield, and cost-effectiveness is maximized by early application in the diagnostic pathway. Pediatricians should consider early referral of children with undiagnosed syndromes to clinical geneticists.

Item does not contain fulltext The devastating clinical presentation of X-linked lissencephaly with abnormal genitalia (XLAG) is invariably caused by loss-of-function mutations in the Aristaless-related homeobox (ARX) gene. Mutations in this X-chromosome gene contribute to intellectual disability (ID) with co-morbidities including seizures and movement disorders such as dystonia in affected males. The detection of affected females with mutations in ARX is increasing. We present a family with multiple affected individuals, including two females. Two male siblings presenting with XLAG were deceased prior to full-term gestation or within the first few weeks of life. Of the two female siblings, one presented with behavioral disturbances, mild ID, a seizure disorder, and complete agenesis of the corpus callosum (ACC), similar to the mother's phenotype. A novel insertion mutation in Exon 2 of ARX was identified, c.982delCinsTTT predicted to cause a frameshift at p.(Q328Ffs* 37). Our finding is consistent with loss-of-function mutations in ARX causing XLAG in hemizygous males and extends the findings of ID and seizures in heterozygous females. We review the reported phenotypes of females with mutations in ARX and highlight the importance of screening ARX in male and female patients with ID, seizures, and in particular with complete ACC.

To prospectively evaluate the diagnostic and clinical utility of singleton whole-exome sequencing (WES) as a first-tier test in infants with suspected monogenic disease. Singleton WES was performed as a first-tier sequencing test in infants recruited from a single pediatric tertiary center. This occurred in parallel with standard investigations, including single- or multigene panel sequencing when clinically indicated. The diagnosis rate, clinical utility, and impact on management of singleton WES were evaluated. Of 80 enrolled infants, 46 received a molecular genetic diagnosis through singleton WES (57.5%) compared with 11 (13.75%) who underwent standard investigations in the same patient group. Clinical management changed following exome diagnosis in 15 of 46 diagnosed participants (32.6%). Twelve relatives received a genetic diagnosis following cascade testing, and 28 couples were identified as being at high risk of recurrence in future pregnancies. This prospective study provides strong evidence for increased diagnostic and clinical utility of singleton WES as a first-tier sequencing test for infants with a suspected monogenic disorder. Singleton WES outperformed standard care in terms of diagnosis rate and the benefits of a diagnosis, namely, impact on management of the child and clarification of reproductive risks for the extended family in a timely manner. Genet Med 18 11, 1090–1096.

Background Previously, we described ROVER, a DNA variant caller which identifies genetic variants from PCR-targeted massively parallel sequencing (MPS) datasets generated by the Hi-Plex protocol. ROVER permits stringent filtering of sequencing chemistry-induced errors by requiring reported variants to appear in both reads of overlapping pairs above certain thresholds of occurrence. ROVER was developed in tandem with Hi-Plex and has been used successfully to screen for genetic mutations in the breast cancer predisposition gene PALB2. ROVER is applied to MPS data in BAM format and, therefore, relies on sequence reads being mapped to a reference genome. In this paper, we describe an improvement to ROVER, called UNDR ROVER (Unmapped primer-Directed ROVER), which accepts MPS data in FASTQ format, avoiding the need for a computationally expensive mapping stage. It does so by taking advantage of the location-specific nature of PCR-targeted MPS data. Results The UNDR ROVER algorithm achieves the same stringent variant calling as its predecessor with a significant runtime performance improvement. In one indicative sequencing experiment, UNDR ROVER (in its fastest mode) required 8-fold less sequential computation time than the ROVER pipeline and 13-fold less sequential computation time than a variant calling pipeline based on the popular GATK tool. UNDR ROVER is implemented in Python and runs on all popular POSIX-like operating systems (Linux, OS X). It requires as input a tab-delimited format file containing primer sequence information, a FASTA format file containing the reference genome sequence, and paired FASTQ files containing sequence reads. Primer sequences at the 5′ end of reads associate read-pairs with their targeted amplicon and, thus, their expected corresponding coordinates in the reference genome. The primer-intervening sequence of each read is compared against the reference sequence from the same location and variants are identified using the same algorithm as ROVER. Specifically, for a variant to be ‘called’ it must appear at the same location in both of the overlapping reads above user-defined thresholds of minimum number of reads and proportion of reads. Conclusions UNDR ROVER provides the same rapid and accurate genetic variant calling as its predecessor with greatly reduced computational costs.