Whole-exome sequencing reveals a C-terminal germline variant in CEBPA-associated acute myeloid leukemia: 45-year follow-up of a large family.

Familial acute myeloid leukemia is rare and linked to germline mutations in RUNX1, GATA2 or CCAAT/enhancer binding protein-α (CEBPA). We re-evaluated a large family with acute myeloid leukemia originally seen at NIH in 1969. We utilized whole-exome sequencing to study this family, and conducted in silico bioinformatics analysis, protein structural modeling and laboratory experiments to assess the impact of the identified CEBPA Q311P mutation. Unlike most previously identified germline mutations in CEBPA, which were N-terminal frameshift mutations, we identified a novel Q311P variant that was located in the C-terminal bZip domain of C/EBPα. Protein structural modeling suggested that the Q311P mutation alters the ability of the CEBPA dimer to bind DNA. Electrophoretic mobility shift assays showed that the Q311P mutant had attenuated binding to DNA, as predicted by the protein modeling. Consistent with these findings, we found that the Q311P mutation has reduced transactivation, consistent with a loss-of-function mutation. From 45 years of follow-up, we observed incomplete penetrance (46%) of CEBPA Q311P. This study of a large multi-generational pedigree reveals that a germline mutation in the C-terminal bZip domain can alter the ability of C/EBP-α to bind DNA and reduces transactivation, leading to acute myeloid leukemia.

New insights into the performance of human whole-exome capture platforms.

Whole exome sequencing (WES) is increasingly used in research and diagnostics. WES users expect coverage of the entire coding region of known genes as well as sufficient read depth for the covered regions. It is, however, unknown which recent WES platform is most suitable to meet these expectations. We present insights into the performance of the most recent standard exome enrichment platforms from Agilent, NimbleGen and Illumina applied to six different DNA samples by two sequencing vendors per platform. Our results suggest that both Agilent and NimbleGen overall perform better than Illumina and that the high enrichment performance of Agilent is stable among samples and between vendors, whereas NimbleGen is only able to achieve vendor- and sample-specific best exome coverage. Moreover, the recent Agilent platform overall captures more coding exons with sufficient read depth than NimbleGen and Illumina. Due to considerable gaps in effective exome coverage, however, the three platforms cannot capture all known coding exons alone or in combination, requiring improvement. Our data emphasize the importance of evaluation of updated platform versions and suggest that enrichment-free whole genome sequencing can overcome the limitations of WES in sufficiently covering coding exons, especially GC-rich regions, and in characterizing structural variants.

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

Whole-Genome Sequencing of a Canine Family Trio Reveals a FAM83G Variant Associated with Hereditary Footpad Hyperkeratosis

Over 250 Mendelian traits and disorders, caused by rare alleles have been mapped in the canine genome. Although each disease is rare in the dog as a species, they are collectively common and have major impact on canine health. With SNP-based genotyping arrays, genome-wide association studies (GWAS) have proven to be a powerful method to map the genomic region of interest when 10-20 cases and 10-20 controls are available. However, to identify the genetic variant in associated regions, fine-mapping and targeted re-sequencing is required. Here we present a new approach using whole-genome sequencing (WGS) of a family trio without prior GWAS. As a proof-of-concept, we chose an autosomal recessive disease known as hereditary footpad hyperkeratosis (HFH) in Kromfohrl änder dogs. To our knowledge, this is the first time this family trio WGS-approach, has successfully been used to identify a genetic variant that perfectly segregates with a canine disorder. The sequencing of three Kromfohrl änder dogs from a family trio (an affected offspring and both its healthy parents) resulted in an average genome coverage of 9.2X per individual. After applying stringent filtering criteria for candidate causative coding variants, 527 single nucleotide variants (SNVs) and 15 indels were found to be homozygous in the affected offspring and heterozygous in the parents. Using the computer software packages ANNOVAR and SIFT to functionally annotate coding sequence differences and to predict their functional effect, resulted in seven candidate variants located in six different genes. Of these, only FAM83G:c155G>C (p.R52P) was found to be concordant in eight additional cases and 16 healthy Kromfohrl änder dogs.

ENU-induced phenovariance in mice: inferences from 587 mutations

Background We present a compendium of N-ethyl-N-nitrosourea (ENU)-induced mouse mutations, identified in our laboratory over a period of 10 years either on the basis of phenotype or whole genome and/or whole exome sequencing, and archived in the Mutagenetix database. Our purpose is threefold: 1) to formally describe many point mutations, including those that were not previously disclosed in peer-reviewed publications; 2) to assess the characteristics of these mutations; and 3) to estimate the likelihood that a missense mutation induced by ENU will create a detectable phenotype. Findings In the context of an ENU mutagenesis program for C57BL/6J mice, a total of 185 phenotypes were tracked to mutations in 129 genes. In addition, 402 incidental mutations were identified and predicted to affect 390 genes. As previously reported, ENU shows strand asymmetry in its induction of mutations, particularly favoring T to A rather than A to T in the sense strand of coding regions and splice junctions. Some amino acid substitutions are far more likely to be damaging than others, and some are far more likely to be observed. Indeed, from among a total of 494 non-synonymous coding mutations, ENU was observed to create only 114 of the 182 possible amino acid substitutions that single base changes can achieve. Based on differences in overt null allele frequencies observed in phenotypic vs. non-phenotypic mutation sets, we infer that ENU-induced missense mutations create detectable phenotype only about 1 in 4.7 times. While the remaining mutations may not be functionally neutral, they are, on average, beneath the limits of detection of the phenotypic assays we applied. Conclusions Collectively, these mutations add to our understanding of the chemical specificity of ENU, the types of amino acid substitutions it creates, and its efficiency in causing phenovariance. Our data support the validity of computational algorithms for the prediction of damage caused by amino acid substitutions, and may lead to refined predictions as to whether specific amino acid changes are responsible for observed phenotypes. These data form the basis for closer in silico estimations of the number of genes mutated to a state of phenovariance by ENU within a population of G3 mice.

Disease evolution and outcomes in familial AML with germline CEBPA mutations

In-depth molecular investigation of familial leukemia has been limited by the rarity of recognized cases. This study examines the genetic events initiating leukemia and details the clinical progression of disease across multiple families harboring germ-line CEBPA mutations. Clinical data were collected from 10 CEBPA-mutated families, representing 24 members with acute myeloid leukemia (AML). Whole-exome (WES) and deep sequencing were performed to genetically profile tumors and define patterns of clonal evolution. Germline CEBPA mutations clustered within the N-terminal and were highly penetrant, with AML presenting at a median age of 24.5 years (range, 1.75-46 years). In all diagnostic tumors tested (n = 18), double CEBPA mutations (CEBPAdm) were detected, with acquired (somatic) mutations preferentially targeting the C-terminal. Somatic CEBPA mutations were unstable throughout the disease course, with different mutations identified at recurrence. Deep sequencing of diagnostic and relapse paired samples confirmed that relapse-associated CEBPA mutations were absent at diagnosis, suggesting recurrence was triggered by novel, independent clones. Integrated WES and deep sequencing subsequently revealed an entirely new complement of mutations at relapse, verifying the presentation of a de novo leukemic episode. The cumulative incidence of relapse in familial AML was 56% at 10 years (n = 11), and 3 patients experienced ≥3 disease episodes over a period of 17 to 20 years. Durable responses to secondary therapies were observed, with prolonged median survival after relapse (8 years) and long-term overall survival (10-year overall survival, 67%). Our data reveal that familial CEBPA-mutated AML exhibits a unique model of disease progression, associated with favorable long-term outcomes.

A Staphylococcus xylosus isolate with a new mecC allotype

Recently, a novel variant of mecA known as mecC (mecA(LGA251)) was identified in Staphylococcus aureus isolates from both humans and animals. In this study, we identified a Staphylococcus xylosus isolate that harbors a new allotype of the mecC gene, mecC1. Whole-genome sequencing revealed that mecC1 forms part of a class E mec complex (mecI-mecR1-mecC1-blaZ) located at the orfX locus as part of a likely staphylococcal cassette chromosome mec element (SCCmec) remnant, which also contains a number of other genes present on the type XI SCCmec.

Novel pseudo-staphylococcal cassette chromosome mec element (ψSCCmec57395) in methicillin-resistant Staphylococcus pseudintermedius CC45

Genetic characterization of methicillin-resistant Staphylococcus pseudintermedius (MRSP) from Thailand and Israel revealed the presence of a predominant atypical clonal lineage which was not typeable by SmaI-PFGE and SCCmec typing. All the atypical isolates (n = 34) belonged to CC45 (30 ST45 and 2 ST179 isolates, 1 ST57 isolate, and 1 ST85 isolate). The isolates originated from healthy and diseased dogs and cats, as well as from the environment of one clinic. Cfr9I-pulsed-field gel electrophoresis (Cfr9I-PFGE) and dru typing permitted the further distinction of CC45 isolates from the two different countries. Microarray analysis identified genes that confer resistance to β-lactams (mecA; blaZ), aminoglycosides [aac(6')-Ie-aph(2')-Ia; aph(3')-III; ant(6)-Ia], macrolides and lincosamides [erm(B)], tetracyclines [tet(M)], trimethoprim [dfr(G)], streptothricin (sat4), and chloramphenicol (catpC221). Fluoroquinolone resistance was attributed to specific amino acid substitutions, i.e., Ser84Leu in GyrA and Ser80Ile and Asp84Asn in GrlA. A novel pseudo-staphylococcal cassette chromosome (ΨSCCmec57395) element was identified in MRSP strain 57395 (sequence type ST45) by whole-genome sequencing. The 12,282-bp ΨSCCmec57395 element contained a class C1 mec gene complex but no ccr genes. In addition to the methicillin resistance gene mecA, ΨSCCmec57395 also carried determinants of resistance to heavy metals, such as arsenic, cadmium, and copper. Bsu36I restriction analysis of the ΨSCCmec57395 element amplified by long-range PCR revealed the presence of ΨSCCmec57395 in the 33 additional isolates of MRSP CC45. The ΨSCCmec57395 element represents a new class of SCCmec and has been identified in MRSP of CC45, which is a predominant clonal lineage in Israel and Thailand.

Identification of a SIRT1 mutation in a family with type 1 diabetes.

Type 1 diabetes is caused by autoimmune-mediated β cell destruction leading to insulin deficiency. The histone deacetylase SIRT1 plays an essential role in modulating several age-related diseases. Here we describe a family carrying a mutation in the SIRT1 gene, in which all five affected members developed an autoimmune disorder: four developed type 1 diabetes, and one developed ulcerative colitis. Initially, a 26-year-old man was diagnosed with the typical features of type 1 diabetes, including lean body mass, autoantibodies, T cell reactivity to β cell antigens, and a rapid dependence on insulin. Direct and exome sequencing identified the presence of a T-to-C exchange in exon 1 of SIRT1, corresponding to a leucine-to-proline mutation at residue 107. Expression of SIRT1-L107P in insulin-producing cells resulted in overproduction of nitric oxide, cytokines, and chemokines. These observations identify a role for SIRT1 in human autoimmunity and unveil a monogenic form of type 1 diabetes.

NGS Nominated CELA1, HSPG2, and KCNK5 as Candidate Genes for Predisposition to Balkan Endemic Nephropathy

Balkan endemic nephropathy (BEN) is a familial chronic tubulointerstitial disease with insidious onset and slow progression leading to terminal renal failure. The results of molecular biological investigations propose that BEN is a multifactorial disease with genetic predisposition to environmental risk agents. Exome sequencing of 22 000 genes with Illumina Nextera Exome Enrichment Kit was performed on 22 DNA samples (11 Bulgarian patients and 11 Serbian patients). Software analysis was performed via NextGene, Provean, and PolyPhen. The frequency of all annotated genetic variants with deleterious/damaging effect was compared with those of European populations. Then we focused on nonannotated variants (with no data available about them and not found in healthy Bulgarian controls). There is no statistically significant difference between annotated variants in BEN patients and European populations. From nonannotated variants with more than 40% frequency in both patients' groups, we nominated 3 genes with possible deleterious/damaging variants-CELA1, HSPG2, and KCNK5. Mutant genes (CELA1, HSPG2, and KCNK5) in BEN patients encode proteins involved in basement membrane/extracellular matrix and vascular tone, tightly connected to process of angiogenesis. We suggest that an abnormal process of angiogenesis plays a key role in the molecular pathogenesis of BEN.

A point mutation in cpsE renders Streptococcus pneumoniae nonencapsulated and enhances its growth, adherence and competence

BackgroundThe polysaccharide capsule is a major virulence factor of the important human pathogen Streptococcus pneumoniae. However, S. pneumoniae strains lacking capsule do occur.ResultsHere, we report a nasopharyngeal isolate of Streptococcus pneumoniae composed of a mixture of two phenotypes; one encapsulated (serotype 18C) and the other nonencapsulated, determined by serotyping, electron microscopy and fluorescence isothiocyanate dextran exclusion assay.By whole genome sequencing, we demonstrated that the phenotypes differ by a single nucleotide base pair in capsular gene cpsE (C to G change at gene position 1135) predicted to result in amino acid change from arginine to glycine at position 379, located in the cytoplasmic, enzymatically active, region of this transmembrane protein. This SNP is responsible for loss of capsule production as the phenotype is transferred with the capsule operon. The nonencapsulated variant is superior in growth in vitro and is also 117-fold more adherent to and more invasive into Detroit 562 human epithelial cells than the encapsulated variant.Expression of six competence pathway genes and one competence-associated gene was 11 to 34-fold higher in the nonencapsulated variant than the encapsulated and transformation frequency was 3.7-fold greater.ConclusionsWe identified a new single point mutation in capsule gene cpsE of a clinical S. pneumoniae serotype 18C isolate sufficient to cause loss of capsule expression resulting in the co-existence of the encapsulated and nonencapsulated phenotype. The mutation caused phenotypic changes in growth, adherence to epithelial cells and transformability. Mutation in capsule gene cpsE may be a way for S. pneumoniae to lose its capsule and increase its colonization potential.

Global phylogenomic analysis of nonencapsulated Streptococcus pneumoniae reveals a deep-branching classic lineage that is distinct from multiple sporadic lineages.

The surrounding capsule of Streptococcus pneumoniae has been identified as a major virulence factor and is targeted by pneumococcal conjugate vaccines (PCV). However, nonencapsulated Streptococcus pneumoniae (Non-Ec-Sp) have also been isolated globally, mainly in carriage studies. It is unknown if Non-Ec-Sp evolve sporadically, if they have high antibiotic non-susceptiblity rates and a unique, specific gene content. Here, whole genome sequencing of 131 Non-Ec-Sp isolates sourced from 17 different locations around the world was performed. Results revealed a deep-branching classic lineage that is distinct from multiple sporadic lineages. The sporadic lineages clustered with a previously sequenced, global collection of encapsulated S. pneumoniae (Ec-Sp) isolates while the classic lineage is comprised mainly of the frequently identified multi-locus sequences types ST344 (n=39) and ST448 (n=40). All ST344 and nine ST448 isolates had high non-susceptiblity rates to β-lactams and other antimicrobials. Analysis of the accessory genome reveals that the classic Non-Ec-Sp contained an increased number of mobile elements, than Ec-Sp and sporadic Non-Ec-Sp. Performing adherence assays to human epithelial cells for selected classic and sporadic Non-Ec-Sp revealed that the presence of a integrative conjugative element (ICE) results in increased adherence to human epithelial cells (P=0.005). In contrast, sporadic Non-Ec-Sp lacking the ICE had greater growth in vitro possibly resulting in improved fitness. In conclusion, Non-Ec-Sp isolates from the classic lineage have evolved separately. They have spread globally, are well adapted to nasopharyngeal carriage and are able to coexist with Ec-Sp. Due to continued use of pneumococcal conjugate vaccines, Non-Ec-Sp may become more prevalent.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Looking the cow in the eye: deletion in the NID1 gene is associated with recessive inherited cataract in Romagnola cattle

Cataract is a known condition leading to opacification of the eye lens causing partial or total blindness. Mutations are known to cause autosomal dominant or recessive inherited forms of cataracts in humans, mice, rats, guinea pigs and dogs. The use of large-sized animal models instead of those using mice for the study of this condition has been discussed due to the small size of rodent lenses. Four juvenile-onset cases of bilateral incomplete immature nuclear cataract were recently observed in Romagnola cattle. Pedigree analysis suggested a monogenic autosomal recessive inheritance. In addition to the cataract, one of the cases displayed abnormal head movements. Genome-wide association and homozygosity mapping and subsequent whole genome sequencing of a single case identified two perfectly associated sequence variants in a critical interval of 7.2 Mb on cattle chromosome 28: a missense point mutation located in an uncharacterized locus and an 855 bp deletion across the exon 19/intron 19 border of the bovine nidogen 1 (NID1) gene (c.3579_3604+829del). RT-PCR showed that NID1 is expressed in bovine lenses while the transcript of the second locus was absent. The NID1 deletion leads to the skipping of exon 19 during transcription and is therefore predicted to cause a frameshift and premature stop codon (p.1164fs27X). The truncated protein lacks a C-terminal domain essential for binding with matrix assembly complexes. Nidogen 1 deficient mice show neurological abnormalities and highly irregular crystal lens alterations. This study adds NID1 to the list of candidate genes for inherited cataract in humans and is the first report of a naturally occurring mutation leading to non-syndromic catarct in cattle provides a potential large animal model for human cataract.