Search

Lynch syndrome (LS) is the most common inherited cancer syndrome. It is inherited via a monoallelic germline variant in one of the DNA mismatch repair (MMR) genes. LS carriers have a broad 30-80% risk of developing various malignancies, and more precise, individual risk estimations would be of high clinical value, allowing tailored cancer prevention and surveillance. Due to MMR deficiency, LS cancers are characterized by the accumulation of frameshift mutations leading to highly immunogenic frameshift peptides (FSPs). Thus, immune surveillance is proposed to inhibit the outgrowth of MMR-deficient cell clones. Recent studies have shown that immunoediting during the evolution of MMR-deficient cancers leads to a counter-selection of highly immunogenic antigens. The immunogenicity of FSPs is dependent on the antigen presentation. One crucial factor determining antigen presentation is the HLA genotype. Hence, a LS carrier's HLA genotype plays an important role in the presentation of FSP antigens to the immune system, and may influence the likelihood of progression from pre-cancerous lesions to cancer. To address the challenge of clarifying this possibility including diverse populations with different HLA types, we have established the INDICATE initiative (Individual cancer risk by HLA type, http://indicate-lynch.org/), an international network aiming at a systematic evaluation of the HLA genotype as a possible cancer risk modifier in LS. Here we summarize the current knowledge on the role of HLA type in cancer risk and outline future research directions to delineate possible association in the scenario of LS with genetically defined risk population and highly immunogenic tumors. peerReviewed

BACKGROUND & AIMS: Constitutional mismatch repair deficiency (CMMRD) is a rare recessive childhood cancer predisposition syndrome caused by germline mismatch repair variants. Constitutional microsatellite instability (cMSI) is a CMMRD diagnostic hallmark and may associate with cancer risk. We quantified cMSI in a large CMMRD patient cohort to explore genotype-phenotype correlations using novel MSI markers selected for instability in blood.METHODS: Three CMMRD, 1 Lynch syndrome, and 2 control blood samples were genome sequenced to >120× depth. A pilot cohort of 8 CMMRD and 38 control blood samples and a blinded cohort of 56 CMMRD, 8 suspected CMMRD, 40 Lynch syndrome, and 43 control blood samples were amplicon sequenced to 5000× depth. Sample cMSI score was calculated using a published method comparing microsatellite reference allele frequencies with 80 controls.RESULTS: Thirty-two mononucleotide repeats were selected from blood genome and pilot amplicon sequencing data. cMSI scoring using these MSI markers achieved 100% sensitivity (95% CI, 93.6%-100.0%) and specificity (95% CI 97.9%-100.0%), was reproducible, and was superior to an established tumor MSI marker panel. Lower cMSI scores were found in patients with CMMRD with MSH6 deficiency and patients with at least 1 mismatch repair missense variant, and patients with biallelic truncating/copy number variants had higher scores. cMSI score did not correlate with age at first tumor.CONCLUSIONS: We present an inexpensive and scalable cMSI assay that enhances CMMRD detection relative to existing methods. cMSI score is associated with mismatch repair genotype but not phenotype, suggesting it is not a useful predictor of cancer risk.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed by computational detection tools. During the last decade, a plethora of such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools were used and detected over 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were empirically validated using three orthogonal methods. Generally, tool-specific precision values are high and similar (median of 98.8%, 96.3%, and 95.5% for qPCR, RNase R, and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant tool differentiators. Furthermore, we demonstrate the complementarity of tools through the increase in detection sensitivity by considering the union of highly-precise tool combinations while keeping the number of false discoveries low. Finally, based on the benchmarking results, recommendations are put forward for circRNA detection and validation.

Identification of mismatch repair (MMR)-deficient colorectal cancers (CRCs) is recommended for Lynch syndrome (LS) screening, and supports targeting of immune checkpoint inhibitors. Microsatellite instability (MSI) analysis is commonly used to test for MMR deficiency. Testing biopsies prior to tumour resection can inform surgical and therapeutic decisions, but can be limited by DNA quantity. MSI analysis of voided urine could also provide much needed surveillance for genitourinary tract cancers in LS. Here, we reconfigure an existing molecular inversion probe-based MSI and BRAF c.1799T > A assay to a multiplex PCR (mPCR) format, and demonstrate that it can sample >140 unique molecules per marker from

Purpose Biallelic germline mismatch repair (MMR) gene pathogenic variants (PVs) cause constitutional MMR deficiency (CMMRD), a highly penetrant childhood cancer syndrome phenotypically overlapping with neurofibromatosis type 1 (NF1). CMMRD testing in suspected NF1 children without NF1/SPRED1 PVs enables inclusion of CMMRD positives into monitoring programs prior to tumor onset. However, testing is associated with potential harms and the prevalence of CMMRD among these children is unknown. Methods Using a simple and scalable microsatellite instability (MSI) assay of non-neoplastic leukocyte DNA to detect CMMRD, we retrospectively screened >700 children suspected of sporadic NF1 but lacking NF1/SPRED1 PVs. Results For three of seven MSI-positive patients germline MMR gene PVs confirmed the diagnosis of CMMRD. Founder variants NM_000535.5(PMS2):c.736_741delinsTGTGTGTGAAG, prevalent in Europe and North America, and NM_000179.2(MSH6):c.10C>G, affecting 1:400 French Canadians, represented two of five PVs. The prevalence of CMMRD was 3/735 (0.41%, 95% confidence interval [CI]: 0.08–1.19%). Conclusion Our empirical data provide reliable numbers for genetic counseling and confirm previous prevalence estimations, on which Care for CMMRD consortium guidelines are based. These advocate CMMRD testing of preselected patients rather than offering reflex testing to all suspected sporadic NF1 children lacking NF1/SPRED1 PVs. The possibility of founder effects should be considered alongside these testing guidelines.

Interactions between the products of the nuclear and mitochondrial genomes are critical for the function of most eukaryotic cells. Recently the introduction of mitochondrial replacement therapy has raised the question of incompatibilities between mitochondrial and nuclear variants, and their potential influence on the genetic makeup of human populations. Such interactions could also contribute to the variability of the penetrance of pathogenic DNA variants. This led us to investigate the frequencies of combinations of nuclear and mitochondrial SNP alleles (mitonuclear combinations) in healthy individuals (n=5375) and in a cohort of patients with Parkinson’s disease (PD, n=2210). In the unaffected population, we were not able to find associations between nuclear and mitochondrial variants with a false discovery rate below 0.05 after accounting for multiple testing (i.e., the number of combinations examined). However, in the PD cohort, five combinations surpassed this threshold. Next, combining both cohorts, we investigated whether these associations were modulated by disease status. All five combinations were significant (p−3 for all tests).These combinations also showed significant evidence for an effect of the interaction between the mitochondrial and nuclear variants on disease risk. Their nuclear components mapped to TBCA, NIBAN3 and GLT25D1 and an uncharacterised intergenic region. In summary, starting from a single cohort design we identified combinations of nuclear and mitochondrial variants affecting PD disease risk.

Microsatellite instability (MSI) testing of colorectal cancers (CRCs) is used to screen for Lynch syndrome (LS), a hereditary cancer‐predisposition, and can be used to predict response to immunotherapy. Here, we present a single‐molecule molecular inversion probe and sequencing‐based MSI assay and demonstrate its clinical validity according to existing guidelines. We amplified 24 microsatellites in multiplex and trained a classifier using 98 CRCs, which accommodates marker specific sensitivities to MSI. Sample classification achieved 100% concordance with the MSI Analysis System v1.2 (Promega) in three independent cohorts, totaling 220 CRCs. Backward–forward stepwise selection was used to identify a 6‐marker subset of equal accuracy to the 24‐marker panel. Assessment of assay detection limits showed that the 24‐marker panel is marginally more robust to sample variables than the 6‐marker subset, detecting as little as 3% high levels of MSI DNA in sample mixtures, and requiring a minimum of 10 template molecules to be sequenced per marker for >95% accuracy. BRAF c.1799 mutation analysis was also included to streamline LS testing, with all c.1799T>A variants being correctly identified. The assay, therefore, provides a cheap, robust, automatable, and scalable MSI test with internal quality controls, suitable for clinical cancer diagnostics.

IntroductionDe novo mutations (DNMs) are known to play a prominent role in sporadic disorders with reduced fitness1. We hypothesize that DNMs play an important role in male infertility and explain a significant fraction of the genetic causes of this understudied disorder. To test this hypothesis, we performed trio-based exome-sequencing in a unique cohort of 185 infertile males and their unaffected parents. Following a systematic analysis, 29 of 145 rare protein altering DNMs were classified as possibly causative of the male infertility phenotype. We observed a significant enrichment of Loss-of-Function (LoF) DNMs in LoF-intolerant genes (p-value=1.00×10-5) as well as predicted pathogenic missense DNMs in missense-intolerant genes (p-value=5.01×10-4). One DNM gene identified, RBM5, is an essential regulator of male germ cell pre-mRNA splicing2. In a follow-up study, 5 rare pathogenic missense mutations affecting this gene were observed in a cohort of 2,279 infertile patients, with no such mutations found in a cohort of 5,784 fertile men (p-value=0.009). Our results provide the first evidence for the role of DNMs in severe male infertility and point to many new candidate genes affecting fertility.

Several studies have linked circulating cell-free mitochondrial DNA (ccf-mtDNA) to human disease. In particular, reduced ccf-mtDNA levels in the cerebrospinal fluid (CSF) of both Alzheimer’s and Parkinson’s disease (PD) patients have raised the hypothesis that ccf-mtDNA could be used as a biomarker for neurodegenerative disease onset and progression. However, how a reduction of CSF ccf-mtDNA levels relates to neurodegeneration remains unclear. Many factors are likely to influence ccf-mtDNA levels, such as concomitant therapeutic treatment and comorbidities. In this study we aimed to investigate these factors, quantifying CSF ccf-mtDNA from the Parkinson’s Progression Markers Initiative in 372 PD patients and 159 matched controls at two time points. We found that ccf-mtDNA levels appear significantly reduced in PD cases when compared to matched controls and are associated with cognitive impairment. However, our data indicate that this reduction in ccf-mtDNA is also associated with the commencement, type and duration of treatment. Additionally, we found that ccf-mtDNA levels are associated with comorbidities such as depression and insomnia, however this was only significant if measured in the absence of treatment. We conclude that in PD, similar to reports in HIV and sepsis, comorbidities and treatment can both influence ccf-mtDNA homeostasis, raising the possibility that ccf-mtDNA may be useful as a biomarker for treatment response or the development of secondary phenotypes. Given that, clinically, PD manifests often decades after neurodegeneration begins, predicting who will develop disease is important. Also, identifying patients who will respond to existing treatments or develop secondary phenotypes will have increased clinical importance as PD incidence rises.

Background Circular RNAs (circRNAs) are predominantly derived from protein coding genes, and some can act as microRNA sponges or transcriptional regulators. Changes in circRNA levels have been identified during human development which may be functionally important, but lineage-specific analyses are currently lacking. To address this, we performed RNAseq analysis of human embryonic stem (ES) cells differentiated for 90 days towards 3D laminated retina. Results A transcriptome-wide increase in circRNA expression, size, and exon count was observed, with circRNA levels reaching a plateau by day 45. Parallel statistical analyses, controlling for sample and locus specific effects, identified 239 circRNAs with expression changes distinct from the transcriptome-wide pattern, but these all also increased in abundance over time. Surprisingly, circRNAs derived from long non-coding RNAs (lncRNAs) were found to account for a significantly larger proportion of transcripts from their loci of origin than circRNAs from coding genes. The most abundant, circRMST:E12-E6, showed a > 100X increase during differentiation accompanied by an isoform switch, and accounts for > 99% of RMST transcripts in many adult tissues. The second most abundant, circFIRRE:E10-E5, accounts for > 98% of FIRRE transcripts in differentiating human ES cells, and is one of 39 FIRRE circRNAs, many of which include multiple unannotated exons. Conclusions Our results suggest that during human ES cell differentiation, changes in circRNA levels are primarily globally controlled. They also suggest that RMST and FIRRE, genes with established roles in neurogenesis and topological organisation of chromosomal domains respectively, are processed as circular lncRNAs with only minor linear species.

Prescription of direct oral anticoagulants (DOAC) compared to warfarin for treating atrial fibrillation patients have increased substantially since their introduction in the England's National Health Service. Assessment of the risk of strokes and bleeds in relation to the large-scale uptake in DOACs compared to warfarin at the clinical commissioning group (CCG) level needs to be carried out. Publicly available- aggregated, CCG level, multi-source health and prescription records data were interrogated to investigate the association between prescription rate of DOACs and stroke/ bleed events during the period of 2013 to 2016. Variability of prescription rates and patient numbers across 208 CCGs were used to infer the effect of DOACs on stroke and bleed risk. Relative risk (RR) and 95% credible intervals (CI) were estimated using Markov chain Monte Carlo approach in JAGS. During the study period, the proportion of DOAC prescriptions increased at an average rate of 122% per annum. DOAC prescription was association with a 50% reduction in ischaemic (RR = 0.48, 95% CI = 0.39, 0.57) and haemorrhagic stroke (RR = 0.50, 95% CI = 0.26-0.77). In contrast, DOAC prescription reached significant association with reduction in gastrointestinal bleeds (RR = 0.86, 95% CI = 0.73-0.98) but not clinically relevant bleeds (RR = 0.95, 95% CI = 0.85-1.05). Sex stratified data showed significant association between DOAC prescription and reduction in haemorrhagic stroke risk (RR = 0.40, 95% CI = 0.28-0.52) and gastrointestinal bleeds (RR = 0.76, 95% CI = 0.63-0.93) in males only. Age stratified data suggested significant association with reduction in risk of both ischaemic and haemorrhagic strokes in patients aged 70 years and above, and reduction in risk of clinically relevant and gastrointestinal bleeds in patients aged 70-79 years only. Publicly available health and prescription data for the English population indicates reduction in stroke and bleed risk in specific age and sex sub-groups with the uptake of DOACs compared to warfarin between 2013 and 2016.

Introduction There is strong evidence from familial recurrence studies for a genetic predisposition to sporadic, non-syndromic Tetralogy of Fallot (TOF). TOF is the most common, cyanotic congenital heart disease (CHD) phenotype yet the cause for the majority of cases remains elusive. Rare genetic variants have been identified as important contributors to the risk of CHD, but relatively small numbers of TOF cases have been studied to date. Methods and Results 829 TOF patients underwent whole exome sequencing (WES), the largest cohort of non-syndromic TOF patients reported to date. The prevalence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database (gnomAD) and a scaled combined annotation-dependent depletion (CADD) score of ≥20. Clustering analysis of variants revealed that two genes, NOTCH1 and FLT4, surpassed thresholds for genome-wide significance (assigned as P Conclusion In summary, the NOTCH1 locus is the most frequent site of genetic variants predisposing to non-syndromic TOF, followed by FLT4. Together, variants in these genes are found in almost 7% of TOF patients. Conflict of Interest None

Rationale: Familial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, nonsyndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease phenotype. Rare genetic variants have been identified as important contributors to the risk of congenital heart disease, but relatively small numbers of TOF cases have been studied to date. Objective: We used whole exome sequencing to assess the prevalence of unique, deleterious variants in the largest cohort of nonsyndromic TOF patients reported to date. Methods and Results: Eight hundred twenty-nine TOF patients underwent whole exome sequencing. The presence of unique, deleterious variants was determined; defined by their absence in the Genome Aggregation Database and a scaled combined annotation-dependent depletion score of ≥20. The clustering of variants in 2 genes, NOTCH1 and FLT4 , surpassed thresholds for genome-wide significance (assigned as P −8 ) after correction for multiple comparisons. NOTCH1 was most frequently found to harbor unique, deleterious variants. Thirty-one changes were observed in 37 probands (4.5%; 95% CI, 3.2%–6.1%) and included 7 loss-of-function variants 22 missense variants and 2 in-frame indels. Sanger sequencing of the unaffected parents of 7 cases identified 5 de novo variants. Three NOTCH1 variants (p.G200R, p.C607Y, and p.N1875S) were subjected to functional evaluation, and 2 showed a reduction in Jagged1-induced NOTCH signaling. FLT4 variants were found in 2.4% (95% CI, 1.6%–3.8%) of TOF patients, with 21 patients harboring 22 unique, deleterious variants. The variants identified were distinct to those that cause the congenital lymphoedema syndrome Milroy disease. In addition to NOTCH1 , FLT4 and the well-established TOF gene, TBX1 , we identified potential association with variants in several other candidates, including RYR1 , ZFPM1 , CAMTA2 , DLX6 , and PCM1 . Conclusions: The NOTCH1 locus is the most frequent site of genetic variants predisposing to nonsyndromic TOF, followed by FLT4 . Together, variants in these genes are found in almost 7% of TOF patients.

Despite the known importance of zinc for human immunity, molecular insights into its roles have remained limited. Here we report a novel autosomal recessive disease characterized by absent B cells, agammaglobulinemia and early onset infections in five unrelated families. The immunodeficiency results from hypomorphic mutations of SLC39A7, which encodes the endoplasmic reticulum-to-cytoplasm zinc transporter ZIP7. Using CRISPR-Cas9 mutagenesis we have precisely modeled ZIP7 deficiency in mice. Homozygosity for a null allele caused embryonic death, but hypomorphic alleles reproduced the block in B cell development seen in patients. B cells from mutant mice exhibited a diminished concentration of cytoplasmic free zinc, increased phosphatase activity and decreased phosphorylation of signaling molecules downstream of the pre-B cell and B cell receptors. Our findings highlight a specific role for cytosolic Zn2+ in modulating B cell receptor signal strength and positive selection.

In many repeat diseases, such as Huntington's disease (HD), ongoing repeat expansions in affected tissues contribute to disease onset, progression and severity. Inducing contractions of expanded repeats by exogenous agents is not yet possible. Traditional approaches would target proteins driving repeat mutations. Here we report a compound, naphthyridine-azaquinolone (NA), that specifically binds slipped-CAG DNA intermediates of expansion mutations, a previously unsuspected target. NA efficiently induces repeat contractions in HD patient cells as well as en masse contractions in medium spiny neurons of HD mouse striatum. Contractions are specific for the expanded allele, independently of DNA replication, require transcription across the coding CTG strand and arise by blocking repair of CAG slip-outs. NA-induced contractions depend on active expansions driven by MutSβ. NA injections in HD mouse striatum reduce mutant HTT protein aggregates, a biomarker of HD pathogenesis and severity. Repeat-structure-specific DNA ligands are a novel avenue to contract expanded repeats.

Context: The pathogenesis of autoimmune Addison's disease (AAD) is thought to be due to interplay of genetic, immune, and environmental factors. A month-of-birth effect, with increased risk for those born in autumn/winter months, has been described in autoimmune conditions such as type 1 diabetes and autoimmune thyroid disease. Objective: Month-of-birth effect was investigated in 2 independent cohorts of AAD subjects. Design, Setting, and Patients: The monthly distribution of birth in AAD patients was compared with that of the general population using the cosinor test. A total of 415 AAD subjects from the United Kingdom cohort were compared with 8 180 180 United Kingdom births, and 231 AAD subjects from the Polish cohort were compared with 2 421 384 Polish births. Main Outcome Measures: Association between month of birth and the susceptibility to AAD. Results: In the entire cohort of AAD subjects, month-of-birth distribution analysis showed significant periodicity with peak of births in December and trough in May (P = .028). Analysis of the odds ratio distribution based on month of birth in 2 cohorts of patients with AAD versus the general population revealed a December peak and May trough, and January peak and July trough, in the United Kingdom and Polish cohorts, respectively. Conclusion: For the first time, we demonstrate that month of birth exerts an effect on the risk of developing AAD, with excess risk in individuals born in winter months and a protective effect when born in the summer. Exposure to seasonal viral infections in the perinatal period, coupled with vitamin D deficiency, could lead to dysregulation of innate immunity affecting the risk of developing AAD., We show for the first time that month of birth exerts an effect on the risk of developing AAD, with excess risk in individuals born in winter months and a protective effect when born in the summer.

PURPOSE\ud \ud To investigate the clinical and functional aspects of MST1 (STK4) deficiency in a profoundly CD4-lymphopenic kindred with a novel homozygous nonsense mutation in STK4. Although recent studies have described the cellular effects of murine Mst1 deficiency, the phenotype of MST1-deficient human lymphocytes has yet to be fully explored. Patient lymphocytes were therefore investigated in the context of current knowledge of murine Mst1 deficiency.\ud \ud METHODS\ud \ud Genetic etiology was identified by whole exome sequencing of genomic DNA from two siblings, combined with linkage analysis in the wider family. MST1 protein expression was assessed by immunoblotting. The ability of patient lymphocytes to adhere to ICAM-1 under flow conditions was measured, and transwell assays were used to assess chemotaxis. Chemokine receptor expression was examined by flow cytometry and receptor signalling by immunoblotting.\ud \ud RESULTS\ud \ud A homozygous nonsense mutation in STK4 (c.442C > T, p.Arg148Stop) was found in the patients, leading to a lack of MST1 protein expression. Patient leukocytes exhibited deficient chemotaxis after stimulation with CXCL11, despite preserved expression of CXCR3. Patient lymphocytes were also unable to bind effectively to immobilised ICAM-1 under flow conditions, in keeping with a failure to develop high affinity binding.\ud \ud CONCLUSION\ud \ud The observed abnormalities of adhesion and migration imply a profound trafficking defect among human MST1-deficient lymphocytes. By analogy with murine Mst1 deficiency and other defects of leucocyte trafficking, this is likely to contribute to immunodeficiency by impairing key aspects of T-cell development and function such as positive selection in the thymus, thymic egress and immune synapse formation in the periphery.

BackgroundClinical guidelines recommend microsatellite instability (MSI) and BRAF V600E testing of all colorectal cancers (CRCs) to screen for Lynch syndrome (LS), a hereditary predisposition to cancer. MSI is also associated with response to immunotherapy. However, uptake of MSI testing is poor and current assays are not suitable for high throughput diagnostics.We aimed to develop a cheap and scalable sequencing assay for MSI classification, which is robust to variables in clinical samples and simultaneously tests for BRAF V600E to streamline the LS screening pipeline.Methods24 short (7-12bp) microsatellites and the BRAF V600E locus were amplified in multiplex using single molecule molecular inversion probes (smMIPs) and sequenced using the Illumina MiSeq platform. Reads were aligned to reference genome hg19. An MSI classifier was trained from 98 CRCs and validated in 99 independent CRCs collected in pathology laboratories in Edinburgh, Spain and Newcastle.ResultsThe smMIP-based MSI assay has 100% accuracy for MSI status relative to MSI Analysis System (Promega). MSI classification is reproducible (100% concordance) and is robust to sample variables, detecting less than 5% MSI-high content in template DNA and giving reliable classification from sequencing only 75 DNA molecules per marker. BRAF V600E was detected with mutant allele frequencies down to 1.7%.ConclusionsOur short microsatellite, smMIP-based, MSI assay provides a cheap and fully automatable assessment of MSI status and BRAF mutation. It is readily scalable to high throughput cancer diagnostics, and is suitable both as a companion diagnostic for immunotherapy and for streamlined LS screening.

Somatic mutations in mononucleotide repeats are commonly used to assess the mismatch repair status of tumours. Current tests focus on repeats with a length above 15bp, which tend to be somatically more unstable than shorter ones. These longer repeats also have a substantially higher PCR error rate and tests that use capillary electrophoresis for fragment size analysis often require expert interpretation. In this communication, we present a panel of 17 short repeats (length 7–12bp) for sequence-based microsatellite instability (MSI) testing. Using a simple scoring procedure that incorporates the allelic distribution of the mutant repeats, and analysis of two cohort of tumours totalling 209 samples, we show that this panel is able to discriminate between MMR proficient and deficient tumours, even when constitutional DNA is not available. In the training cohort, the method achieved 100% concordance with fragment analysis, while in the testing cohort, 4 discordant samples were observed (corresponding to 97% concordance). Of these, 2 showed discrepancies between fragment analysis and immunohistochemistry and one was reclassified after re-testing using fragment analysis. These results indicate that our approach offers the option of a reliable, scalable routine test for MSI.

AimsFamilial recurrence studies provide strong evidence for a genetic component to the predisposition to sporadic, non-syndromic Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease (CHD) phenotype. Rare genetic variants have been identified as important contributors to the risk of CHD, but relatively small numbers of TOF cases have been studied to date. Here, we use whole exome sequencing to assess the prevalence of rare, potentially deleterious variants in candidate genes previously associated with both syndromic and non-syndromic TOF, in the largest cohort of non-syndromic TOF patients reported to date.Methods & Results829 non-syndromic TOF patients underwent whole exome sequencing. A systematic review of the literature was conducted which revealed 77 genes in which mutations had been reported in patients with TOF. The presence of rare, deleterious variants in the 77 candidate genes was determined, defined by a minor allele frequency of ≤ 0.001 and scaled combined annotation-dependent depletion (CADD) score of ≥ 20. We found a clustering of heterozygous rare, deleterious variants inNOTCH1(P=1.89E-15),DOCK6(P=2.93E-07),MYOM2(P= 7.35E-05),TTC37(P=0.016),MESP1(P=0.024) andTBX1(P=0.039), after correcting for multiple testing.NOTCH1was most frequently found to harbour deleterious variants. Changes were observed in 49 patients (6%; 95% confidence interval [CI]: 4.5% - 7.8%) and included six truncating/frameshift variants and forty missense variants. Sanger sequencing of the unaffected parents of thirteen cases identified fivede novovariants. Variants were not confined to a single functional domain of the NOTCH1 protein but significant clustering of variants was evident in the EGF-like repeats (P=0.018). ThreeNOTCH1missense variants (p.G200R, p.C607Y andde novop.N1875S) were subjected to functional evaluation and showed a reduction in Jagged1 ligand-induced NOTCH signalling. p.C607Y, which exhibited the most significant reduction in signalling, also perturbed S1 cleavage of the NOTCH1 receptor in the Golgi.ConclusionTheNOTCH1locus is a frequent site of genetic variants predisposing to non-syndromic TOF with 6% of patients exhibiting rare, deleterious variants. Our data supports the polygenic origin of TOF and suggests larger studies may identify additional loci.

Despite the known importance of zinc for human immunity, molecular insights into its roles have remained limited. Here we report a novel autosomal recessive disease characterized by absent B cells, agammaglobulinemia and early onset infections in five unrelated families. The immunodeficiency results from hypomorphic mutations of SLC39A7, which encodes the endoplasmic reticulum-to-cytoplasm zinc transporter ZIP7. Using CRISPR-Cas9 mutagenesis we have precisely modeled ZIP7 deficiency in mice. Homozygosity for a null allele caused embryonic death, but hypomorphic alleles reproduced the block in B cell development seen in patients. B cells from mutant mice exhibited a diminished concentration of cytoplasmic free zinc, increased phosphatase activity and decreased phosphorylation of signaling molecules downstream of the pre-B cell and B cell receptors. Our findings highlight a specific role for cytosolic Zn

Somatic mutations in mononucleotide repeats are commonly used to assess the mismatch repair status of tumours. Current tests focus on repeats with a length above 15bp, which tend to be somatically more unstable than shorter ones. These longer repeats also have a substantially higher PCR error rate, and tests that use capillary electrophoresis for fragment size analysis often require expert interpretation. In this communication, we present a panel of 17 short repeats (length 7-12bp) for sequence-based microsatellite instability (MSI) testing. Using a simple scoring procedure that incorporates the allelic distribution of the mutant repeats, and analysis of two cohort of tumours totalling 209 samples, we show that this panel is able to discriminate between MMR proficient and deficient tumours, even when constitutional DNA is not available. In the training cohort, the method achieved 100% concordance with fragment analysis, while in the testing cohort, 4 discordant samples were observed (corresponding to 97% concordance). Of these, 2 showed discrepancies between fragment analysis and immunohistochemistry and one was reclassified after re-testing using fragment analysis. These results indicate that our approach offers the option of a reliable, scalable routine test for MSI.

Background Circular RNAs (circRNAs) are predominantly derived from protein coding genes, and some can act as microRNA sponges or transcriptional regulators. Changes in circRNA levels have been identified during human development which may be functionally important, but lineage-specific analyses are currently lacking. To address this, we performed RNAseq analysis of human embryonic stem (ES) cells differentiated for 90 days towards 3D laminated retina. Results A transcriptome-wide increase in circRNA expression, size, and exon count was observed, with circRNA levels reaching a plateau by day 45. Parallel statistical analyses, controlling for sample and locus specific effects, identified 239 circRNAs with expression changes distinct from the transcriptome-wide pattern, but these all also increased in abundance over time. Surprisingly, circRNAs derived from long non-coding RNAs (lncRNAs) were found to account for a significantly larger proportion of transcripts from their loci of origin than circRNAs from coding genes. The most abundant, circRMST:E12-E6, showed a > 100X increase during differentiation accompanied by an isoform switch, and accounts for > 99% of RMST transcripts in many adult tissues. The second most abundant, circFIRRE:E10-E5, accounts for > 98% of FIRRE transcripts in differentiating human ES cells, and is one of 39 FIRRE circRNAs, many of which include multiple unannotated exons. Conclusions Our results suggest that during human ES cell differentiation, changes in circRNA levels are primarily globally controlled. They also suggest that RMST and FIRRE, genes with established roles in neurogenesis and topological organisation of chromosomal domains respectively, are processed as circular lncRNAs with only minor linear species. Electronic supplementary material The online version of this article (10.1186/s12864-018-4660-7) contains supplementary material, which is available to authorized users.

Genome-wide association studies (GWAS) have identified genetic variants in a number of chromosomal regions that are associated with atrial fibrillation (AF). The mechanisms underlying these associations are unknown, but are likely to involve effects of the risk haplotypes on expression of neighbouring genes. To investigate the association between genetic variants at AF-associated loci and expression of nearby candidate genes in human atrial tissue and peripheral blood. Right atrial appendage (RAA) samples were collected from 122 patients undergoing cardiac surgery, of these, 12 patients also had left atrial appendage samples taken. 22 patients had a history of AF. Peripheral blood samples were collected from 405 patients undergoing diagnostic cardiac catheterisation. In order to tag genetic variation at each of nine loci, a total of 367 single nucleotide polymorphisms (SNPs) were genotyped using the Sequenom platform. Total expression of 16 candidate genes in the nine AF-associated regions was measured by quantitative PCR. The relative expression of each allele of the candidate genes was measured on the Sequenom platform using one or more transcribed SNPs to distinguish between alleles in heterozygotes. We tested association between the SNPs of interest and gene expression using total gene expression (integrating cis and trans acting sources of variation), and allelic expression ratios (specific for cis acting influences), in atrial tissue and peripheral blood. We adjusted for multiple comparisons using a Bonferroni approach. In subsidiary analyses, we compared the expression of candidate genes between patients with and without a history of AF. Total expression of 15 transcripts of 14 genes and allelic expression ratio of 14 transcripts of 14 genes in genomic regions associated with AF were measured in right atrial appendage tissue. 8 of these transcripts were also expressed in peripheral blood. Risk alleles at AF-associated SNPs were associated in cis with an increased expression of PITX2a (2.01-fold, p=6.5×10(-4)); and with decreased expression of MYOZ1 (0.39 fold; p=5.5×10(-15)), CAV1 (0.89 fold; p=5.9×10(-8)), C9orf3 (0.91 fold; 1.5×10(-5)), and FANCC (0.94-fold; p=8.9×10(-8)) in right atrial appendage. Of these five genes, only CAV1 was expressed in peripheral blood; association between the same AF risk alleles and lower expression of CAV1 was confirmed (0.91 fold decrease; p=4.2×10(-5)). A history of AF was also associated with a decrease in expression of CAV1 in both right and left atria (0.84 and 0.85 fold, respectively; p=0.03), congruent with the magnitude of the effect of the risk SNP on expression, and independent of genotype. The analyses in peripheral blood showed association between AF risk SNPs and decreased expression of KCNN3 (0.85-fold; p=2.1×10(-4)); and increased expression of SYNE2 (1.12-fold; p=7.5×10(-24)); however, these associations were not detectable in atrial tissue. We identified novel cis-acting associations in atrial tissue between AF risk SNPs and increased expression of PITX2a/b; and decreased expression of CAV1 (an association also seen in peripheral blood), C9orf3 and FANCC. We also confirmed a previously described association between AF risk variants and MYOZ1 expression. Analyses of peripheral blood illustrated tissue-specificity of cardiac eQTLs and highlight the need for larger-scale genome-wide eQTL studies in cardiac tissue. Our results suggest novel aetiological roles for genes in four AF-associated genomic regions.

End-stage renal disease (ESRD) presenting in a familial autosomal dominant pattern points to an underlying monogenic cause. Nail-patella syndrome (NPS) is an autosomal dominant disorder that may lead to ESRD caused by mutations in the transcription factor LMX1B. Renal-limited forms of this disease, termed nail-patella-like renal disease (NPLRD), and LMX1B nephropathy have recently been described. We report a large family, from the North East of England, with seven affected members with varying phenotypes of renal disease, ranging from ESRD at 28 years of age to microscopic haematuria and proteinuria and relatively preserved renal function. In this family, there were no extra-renal manifestations to suggest NPS. Genome-wide linkage studies and inheritance by descent (IBD) suggested disease loci on Chromosome 1 and 9. Whole exome sequencing (WES) analysis identified a novel sequence variant (p.R249Q) in the LMX1B gene in each of the three samples submitted, which was confirmed using Sanger sequencing. The variant segregated with the disease in all affected individuals. In silico modelling revealed that R249 is putatively located in close proximity to the DNA phosphoskeleton, supporting a role for this residue in the interaction between the LMX1B homeodomain and its target DNA. WES and analysis of potential target genes, including CD2AP, NPHS2, COL4A3, COL4A4 and COL4A5, did not reveal any co-inherited pathogenic variants. In conclusion, we confirm a novel LMX1B mutation in a large family with an autosomal dominant pattern of nephropathy. This report confirms that LMX1B mutations may cause a glomerulopathy without extra-renal manifestations. A molecular genetic diagnosis of LMX1B nephropathy thus provides a definitive diagnosis, prevents the need for renal biopsies and allows at risk family members to be screened.

PURPOSE: Mitochondrial disorders are a common cause of inherited metabolic disease and can be due to mutations affecting mitochondrial DNA or nuclear DNA. The current diagnostic approach involves the targeted resequencing of mitochondrial DNA and candidate nuclear genes, usually proceeds step by step, and is time consuming and costly. Recent evidence suggests that variations in mitochondrial DNA sequence can be obtained from whole-exome sequence data, raising the possibility of a comprehensive single diagnostic test to detect pathogenic point mutations. METHODS: We compared the mitochondrial DNA sequence derived from off-target exome reads with conventional mitochondrial DNA Sanger sequencing in 46 subjects. RESULTS: Mitochondrial DNA sequences can be reliably obtained using three different whole-exome sequence capture kits. Coverage correlates with the relative amount of mitochondrial DNA in the original genomic DNA sample, heteroplasmy levels can be determined using variant and total read depths, and-providing there is a minimum read depth of 20-fold-rare sequencing errors occur at a rate similar to that observed with conventional Sanger sequencing. CONCLUSION: This offers the prospect of using whole-exome sequence in a diagnostic setting to screen not only all protein coding nuclear genes but also all mitochondrial DNA genes for pathogenic mutations. Off-target mitochondrial DNA reads can also be used to assess quality control and maternal ancestry, inform on ethnic origin, and allow genetic disease association studies not previously anticipated with existing whole-exome data sets.

Recurrent rearrangements of chromosome 1q21.1 that occur via non-allelic homologous recombination have been associated with variable phenotypes exhibiting incomplete penetrance, including congenital heart disease (CHD). However, the gene or genes within the ~1 Mb critical region responsible for each of the associated phenotypes remains unknown. We examined the 1q21.1 locus in 948 patients with tetralogy of Fallot (TOF), 1488 patients with other forms of CHD and 6760 ethnically matched controls using single nucleotide polymorphism genotyping arrays (Illumina 660W and Affymetrix 6.0) and multiplex ligation-dependent probe amplification. We found that duplication of 1q21.1 was more common in cases of TOF than in controls [odds ratio (OR) 30.9, 95% confidence interval (CI) 8.9-107.6); P = 2.2 × 10(-7)], but deletion was not. In contrast, deletion of 1q21.1 was more common in cases of non-TOF CHD than in controls [OR 5.5 (95% CI 1.4-22.0); P = 0.04] while duplication was not. We also detected rare (n = 3) 100-200 kb duplications within the critical region of 1q21.1 in cases of TOF. These small duplications encompassed a single gene in common, GJA5, and were enriched in cases of TOF in comparison to controls [OR = 10.7 (95% CI 1.8-64.3), P = 0.01]. These findings show that duplication and deletion at chromosome 1q21.1 exhibit a degree of phenotypic specificity in CHD, and implicate GJA5 as the gene responsible for the CHD phenotypes observed with copy number imbalances at this locus.

Vesicoureteric reflux (VUR) is the commonest urological anomaly in children. Despite treatment improvements, associated renal lesions - congenital dysplasia, acquired scarring or both - are a common cause of childhood hypertension and renal failure. Primary VUR is familial, with transmission rate and sibling risk both approaching 50%, and appears highly genetically heterogeneous. It is often associated with other developmental anomalies of the urinary tract, emphasising its etiology as a disorder of urogenital tract development. We conducted a genome-wide linkage and association study in three European populations to search for loci predisposing to VUR. Family-based association analysis of 1098 parent-affected-child trios and case/control association analysis of 1147 cases and 3789 controls did not reveal any compelling associations, but parametric linkage analysis of 460 families (1062 affected individuals) under a dominant model identified a single region, on 10q26, that showed strong linkage (HLOD = 4.90; ZLRLOD = 4.39) to VUR. The ~9Mb region contains 69 genes, including some good biological candidates. Resequencing this region in selected individuals did not clearly implicate any gene but FOXI2, FANK1 and GLRX3 remain candidates for further investigation. This, the largest genetic study of VUR to date, highlights the 10q26 region as a major genetic contributor to VUR in European populations.

Purpose We aimed to achieve a retrospective molecular diagnosis by applying state-of-the-art genomic sequencing methods to past patients with T-B+NK+ severe combined immunodeficiency (SCID). We included identification of copy number variations (CNVs) by whole exome sequencing (WES) using the CNV calling method ExomeDepth to detect gene alterations for which routine Sanger sequencing analysis is not suitable, such as large heterozygous deletions. Methods Of a total of 12 undiagnosed patients with T-B+NK+ SCID, we analyzed eight probands by WES, using GATK to detect single nucleotide variants (SNVs) and small insertions and deletions (INDELs) and ExomeDepth to detect CNVs. Results We found heterozygous single- or multi-exon deletions in IL7R, a known disease gene for autosomal recessive T-B+NK+ SCID, in four families (seven patients). In three families (five patients), these deletions coexisted with a heterozygous splice site or nonsense mutation elsewhere in the same gene, consistent with compound heterozygosity. In our cohort, about a quarter of T-B+NK+ SCID patients (26%) had such compound heterozygous IL7R deletions. Conclusions We show that heterozygous IL7R exon deletions are common in T-B+NK+ SCID and are detectable by WES. They should be considered if Sanger sequencing fails to detect homozygous or compound heterozygous IL7R SNVs or INDELs. Electronic supplementary material The online version of this article (doi:10.1007/s10875-016-0343-9) contains supplementary material, which is available to authorized users.

The demonstration of impaired C regulation in the thrombotic microangiopathy (TMA) atypical hemolytic uremic syndrome (aHUS) resulted in the successful introduction of the C inhibitor eculizumab into clinical practice. C abnormalities account for approximately 50% of aHUS cases; however, mutations in the non-C gene diacylglycerol kinase-ε have been described recently in individuals not responsive to eculizumab. We report here a family in which the proposita presented with aHUS but did not respond to eculizumab. Her mother had previously presented with a post-renal transplant TMA. Both the proposita and her mother also had Charcot-Marie-Tooth disease. Using whole-exome sequencing, we identified a mutation in the inverted formin 2 gene (INF2) in the mutational hotspot for FSGS. Subsequent analysis of the Newcastle aHUS cohort identified another family with a functionally-significant mutation in INF2 In this family, renal transplantation was associated with post-transplant TMA. All individuals with INF2 mutations presenting with a TMA also had aHUS risk haplotypes, potentially accounting for the genetic pleiotropy. Identifying individuals with TMAs who may not respond to eculizumab will avoid prolonged exposure of such individuals to the infectious complications of terminal pathway C blockade.

Objective To quantify any risk of atrial fibrillation (AF) associated with ivabradine treatment by meta-analysis of clinical trial data. Methods Medline, Embase, Web of Knowledge and the Cochrane central register of controlled trials were searched for double-blinded randomised controlled trials of ivabradine with a minimum follow-up period of 4 weeks. For studies where AF data were unpublished, safety data were obtained from the European Medicines Agency (EMeA) website and personal communications. Studies were appraised for risk of bias using components recommended by the Cochrane Collaboration. Meta-analyses were performed of relative risk of AF and absolute risk difference of AF per year of treatment. The main outcome measure was incident AF during the follow-up period. Results AF data were available from 11 studies: one from the published report, six from the EMeA and four from personal communications. Ivabradine treatment was associated with a relative risk of AF of 1.15 (95% CI 1.07 to 1.24, p=0.0027) among 21 571 patients in the meta-analysis. From this we estimated that the number needed to harm for ivabradine would be 208 (95% CI 122 to 667) per year of treatment. Conclusions AF is a substantially more common side effect of ivabradine treatment than one patient in 10 000, the risk presently reported in the product literature. The incidence of AF has not routinely been reported in clinical trials of ivabradine.

Context: The pathogenesis of autoimmune Addison's disease (AAD) is thought to be due to interplay of genetic, immune, and environmental factors. A month-of-birth effect, with increased risk for those born in autumn/winter months, has been described in autoimmune conditions such as type 1 diabetes and autoimmune thyroid disease. Objective: Month-of-birth effect was investigated in 2 independent cohorts of AAD subjects. Design, Setting, and Patients: The monthly distribution of birth in AAD patients was compared with that of the general population using the cosinor test. A total of 415 AAD subjects from the United Kingdom cohort were compared with 8 180 180 United Kingdom births, and 231 AAD subjects from the Polish cohort were compared with 2 421 384 Polish births. Main Outcome Measures: Association between month of birth and the susceptibility to AAD. Results: In the entire cohort of AAD subjects, month-of-birth distribution analysis showed significant periodicity with peak of births in December and troug...

The molecular diagnosis of adult-onset autosomal recessive cerebellar ataxias remains challenging because of genetic heterogeneity. However, recently developed molecular genetic techniques will potentially revolutionize the diagnostic approach. Here we set out to define the genetic basis of the ataxia in two brothers with no molecular diagnosis. Clinical evaluation was followed by whole-exome second-generation sequencing and comparative genomic hybridization to determine the diagnosis. Whole-exome sequencing identified a hemizygous novel spastic ataxia of Charlevoix-Saguenay (SACS) stop-codon mutation in both brothers (c.13048G→T, p.E4350*) that was present in the mother, but not the father. Comparative genomic hybridization revealed a 0.7-Mb deletion on chromosome 13q12.12 in both brothers, which included SACS and was heterozygous in the asymptomatic father. The milder phenotype, caused by a whole-gene deletion and a stop-codon mutation in SACS, indicates a loss-of-function mechanism in late-onset autosomal recessive spastic ataxia of Charlevoix-Saguenay (ARSACS), and illustrates the importance of chromosomal rearrangements in the investigation of adult-onset ataxia.

Objective Resident macrophages play an important role in atheromatous plaque rupture. The macrophage gene expression signature associated with plaque rupture is incompletely defined due to the complex cellular heterogeneity in the plaque. We aimed to characterise differential gene expression in resident plaque macrophages from ruptured and stable human atheromatous lesions. Methods and results We performed genome-wide expression analyses of isolated macrophage-rich regions of stable and ruptured human atherosclerotic plaques. Plaques present in carotid endarterectomy specimens were designated as stable or ruptured using clinical, radiological and histopathological criteria. Macrophage-rich regions were excised from 5 ruptured and 6 stable plaques by laser micro-dissection. Transcriptional profiling was performed using Affymetrix microarrays. The profiles were characteristic of activated macrophages. At a false discovery rate of 10%, 914 genes were differentially expressed between stable and ruptured plaques. The findings were confirmed in fourteen further stable and ruptured samples for a subset of eleven genes with the highest expression differences (p, Highlights ► We examined gene expression in macrophages from stable and ruptured plaques. ► The PPAR/adipocytokine signaling pathway was upregulated in ruptured plaques. ► FABP4 and Leptin were highly expressed in ruptured atheromatous plaque macrophages. ► Down-regulation of PPAR/adipocytokine signaling may have therapeutic potential.

Previous studies have shown that copy-number variants (CNVs) contribute to the risk of complex developmental phenotypes. However, the contribution of global CNV burden to the risk of sporadic congenital heart disease (CHD) remains incompletely defined. We generated genome-wide CNV data by using Illumina 660W-Quad SNP arrays in 2,256 individuals with CHD, 283 trio CHD-affected families, and 1,538 controls. We found association of rare genic deletions with CHD risk (odds ratio [OR] = 1.8, p = 0.0008). Rare deletions in study participants with CHD had higher gene content (p = 0.001) with higher haploinsufficiency scores (p = 0.03) than they did in controls, and they were enriched with Wnt-signaling genes (p = 1 × 10(-5)). Recurrent 15q11.2 deletions were associated with CHD risk (OR = 8.2, p = 0.02). Rare de novo CNVs were observed in ∼5% of CHD trios; 10 out of 11 occurred on the paternally transmitted chromosome (p = 0.01). Some of the rare de novo CNVs spanned genes known to be involved in heart development (e.g., HAND2 and GJA5). Rare genic deletions contribute ∼4% of the population-attributable risk of sporadic CHD. Second to previously described CNVs at 1q21.1, deletions at 15q11.2 and those implicating Wnt signaling are the most significant contributors to the risk of sporadic CHD. Rare de novo CNVs identified in CHD trios exhibit paternal origin bias. ispartof: American Journal of Human Genetics vol:91 issue:3 pages:489-501 ispartof: location:United States status: published