Identification of two novel mutations in POU4F3 gene associated with autosomal dominant hearing loss in Chinese families

Autosomal dominant non‐syndromic hearing loss is genetically heterogeneous with 47 genes identified to date, including POU4F3. In this study, by using a next‐generation sequencing panel targeting 127 deafness genes, we identified a pathogenic frameshift mutation c.704_705del and a missense mutation c.593G>A in two three‐generation Chinese families with late‐onset progressive ADNSHL, respectively. The novel mutations of POU4F3 co‐segregated with the deafness phenotype in these two families. c.704_705del caused a frameshift p.T235fs and c.593G>A caused an amino acid substitution of p.R198H. Both mutations led to an abnormal and incomplete protein structure. POU4F3 with either of the two mutations was transiently transfected into HEI‐OC1 and HEK 293 cell lines and immunofluorescence assay was performed to investigate the subcellular localization of mutated protein. The results indicated that both c.704_705del (p.T235fs) and c.593G>A (p.R198H) could impair the nuclear localization function of POU4F3. The p.R198H POU4F3 protein was detected as a weak band of the correct molecular weight, indicating that the stability of p.R198H POU4F3 differed from that of the wild‐type protein. While, the p.T235fs POU4F3 protein was expressed with a smaller molecular weight, implying this mutation result in a frameshift and premature termination of the POU4F3 protein. In summary, we report two novel mutations of POU4F3 associated with progressive ADNSHL and explored their effects on POU4F3 nuclear localization. These findings expanded the mutation spectrum of POU4F3 and provided new knowledge for the pathogenesis of POU4F3 in hearing loss.     

Assessing structural variation in a personal genome—towards a human reference diploid genome

Background

Characterizing large genomic variants is essential to expanding the research and clinical applications of genome sequencing. While multiple data types and methods are available to detect these structural variants (SVs), they remain less characterized than smaller variants because of SV diversity, complexity, and size. These challenges are exacerbated by the experimental and computational demands of SV analysis. Here, we characterize the SV content of a personal genome with Parliament, a publicly available consensus SV-calling infrastructure that merges multiple data types and SV detection methods.

Results

We demonstrate Parliament’s efficacy via integrated analyses of data from whole-genome array comparative genomic hybridization, short-read next-generation sequencing, long-read (Pacific BioSciences RSII), long-insert (Illumina Nextera), and whole-genome architecture (BioNano Irys) data from the personal genome of a single subject (HS1011). From this genome, Parliament identified 31,007 genomic loci between 100 bp and 1 Mbp that are inconsistent with the hg19 reference assembly. Of these loci, 9,777 are supported as putative SVs by hybrid local assembly, long-read PacBio data, or multi-source heuristics. These SVs span 59 Mbp of the reference genome (1.8%) and include 3,801 events identified only with long-read data. The HS1011 data and complete Parliament infrastructure, including a BAM-to-SV workflow, are available on the cloud-based service DNAnexus.

Conclusions

HS1011 SV analysis reveals the limits and advantages of multiple sequencing technologies, specifically the impact of long-read SV discovery. With the full Parliament infrastructure, the HS1011 data constitute a public resource for novel SV discovery, software calibration, and personal genome structural variation analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1479-3) contains supplementary material, which is available to authorized users.     

Targeted massive parallel sequencing: the effective detection of novel causative mutations associated with hearing loss in small families

ABSTRACT: BACKGROUND: Hereditary hearing loss is one of the most common heterogeneous disorders, and genetic variants that can cause hearing loss have been identified in over fifty genes. Most of these hearing loss genes have been detected using classical genetic methods, typically starting with linkage analysis in large families with hereditary hearing loss. However, these classical strategies are not well suited for mutation analysis in smaller families who have insufficient genetic information. METHODS: Eighty known hearing loss genes were selected and simultaneously sequenced by targeted next-generation sequencing (NGS) in 8 Korean families with autosomal dominant non-syndromic sensorineural hearing loss. RESULTS: Five mutations in known hearing loss genes, including 1 nonsense and 4 missense mutations, were identified in 5 different genes (ACTG1, MYO1F, DIAPH1, POU4F3 and EYA4), and the genotypes for these mutations were consistent with the autosomal dominant inheritance pattern of hearing loss in each family. No mutational hot-spots were revealed in these Korean families. CONCLUSION: Targeted NGS allowed for the detection of pathogenic mutations in affected individuals who were not candidates for classical genetic studies. This report is the first documenting the effective use of an NGS technique to detect pathogenic mutations that underlie hearing loss in an East Asian population. Using this NGS technique to establish a database of common mutations in Korean patients with hearing loss and further data accumulation will contribute to the early diagnosis and fundamental therapies for hereditary hearing loss.     

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6^Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.     

Effects of single and combined genotypes of MC4R and POU1F1 genes on two production traits in Langshan chicken

The objective of this study was to analyze the effects of single and combined genotypes of MC4R and POU1F1 genes in Chinese well-known indigenous chicken (Langshan chicken) population. Genetic variants within MC4R gene and POU1F1 gene were screened through PCR-SSCP and DNA sequencing methods. A C/T mutation at nt 944 in MC4R gene (NC_006089.2:g. 944C>T) and a G/A mutation at nt 3109 in POU1F1 gene (NC_006088.2:g. 3109 G>A) were identified. Associations between the mutations of the two genes with two production traits were analyzed. The results showed that, at MC4R locus, individuals with BB and AB genotypes had highly significantly higher body weight at 16 weeks (p < 0.01) than did those with the AA genotype. And, individuals within AA and AB genotypes had significantly higher egg numbers at 300 days (p < 0.05). At POU1F1 locus, individuals with CD genotype had higher body weight at 16 weeks and egg numbers at 300 days (p < 0.05). Furthermore, combined genotypes from these two loci were found to be associated with egg numbers at 300 days (p < 0.05). The individuals within combined genotype AB/CD had higher egg production. Therefore, variations identified within the MC4R and POU1F1genes are suitable for future use in identifying chickens with the genetic potential of higher body weight and reproductive traits, at least in the population of Langshan chickens.     

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type—which included inversions, duplications, deletions, translocations, and mobile element insertions—was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.     

Targeted Exome Sequencing Integrated with Clinicopathological Information Reveals Novel and Rare Mutations in Atypical,Suspected and Unknown Cases of Alport Syndrome or Proteinuria

We applied customized targeted next-generation exome sequencing (NGS) to determine if mutations in genes associated with renal malformations, Alport syndrome (AS) or nephrotic syndrome are a potential cause of renal abnormalities in patients with equivocal or atypical presentation. We first sequenced 4,041 exons representing 292 kidney disease genes in a Caucasian woman with a history of congenital vesicoureteral reflux (VUR), recurrent urinary tract infections and hydronephrosis who presented with nephrotic range proteinuria at the age of 45. Her biopsy was remarkable for focal segmental glomerulosclerosis (FSGS), a potential complication of longstanding VUR. She had no family history of renal disease. Her proteinuria improved initially, however, several years later she presented with worsening proteinuria and microhematuria. NGS analysis revealed two deleterious COL4A3 mutations, one novel and the other previously reported in AS, and a novel deleterious SALL2 mutation, a gene linked to renal malformations. Pedigree analysis confirmed that COL4A3 mutations were nonallelic and compound heterozygous. The genomic results in conjunction with subsequent abnormal electron microscopy, Collagen IV minor chain immunohistochemistry and progressive sensorineural hearing loss confirmed AS. We then modified our NGS approach to enable more efficient discovery of variants associated with AS or a subset of FSGS by multiplexing targeted exome sequencing of 19 genes associated with AS or FSGS in 14 patients. Using this approach, we found novel or known COL4A3 or COL4A5 mutations in a subset of patients with clinically diagnosed or suspected AS, APOL1 variants associated with FSGS in African Americans and novel mutations in genes associated with nephrotic syndrome. These studies demonstrate the successful application of targeted capture-based exome sequencing to simultaneously evaluate genetic variations in many genes in patients with complex renal phenotypes and provide insights into etiology of conditions with equivocal clinical and pathologic presentations.     

SNP Linkage Analysis and Whole Exome Sequencing Identify a Novel POU4F3 Mutation in Autosomal Dominant Late-Onset Nonsyndromic Hearing Loss (DFNA15)

Autosomal dominant non-syndromic hearing loss (AD-NSHL) is one of the most common genetic diseases in human and is well-known for the considerable genetic heterogeneity. In this study, we utilized whole exome sequencing (WES) and linkage analysis for direct genetic diagnosis in AD-NSHL. The Korean family had typical AD-NSHL running over 6 generations. Linkage analysis was performed by using genome-wide single nucleotide polymorphism (SNP) chip and pinpointed a genomic region on 5q31 with a significant linkage signal. Sequential filtering of variants obtained from WES, application of the linkage region, bioinformatic analyses, and Sanger sequencing validation identified a novel missense mutation Arg326Lys (c.977G>A) in the POU homeodomain of the POU4F3 gene as the candidate disease-causing mutation in the family. POU4F3 is a known disease gene causing AD-HSLH (DFNA15) described in 5 unrelated families until now each with a unique mutation. Arg326Lys was the first missense mutation affecting the 3^rd alpha helix of the POU homeodomain harboring a bipartite nuclear localization signal sequence. The phenotype findings in our family further supported previously noted intrafamilial and interfamilial variability of DFNA15. This study demonstrated that WES in combination with linkage analysis utilizing bi-allelic SNP markers successfully identified the disease locus and causative mutation in AD-NSHL.     

Building a Robust Tumor Profiling Program: Synergy between Next-Generation Sequencing and Targeted Single-Gene Testing

Next-generation sequencing (NGS) is a powerful platform for identifying cancer mutations. Routine clinical adoption of NGS requires optimized quality control metrics to ensure accurate results. To assess the robustness of our clinical NGS pipeline, we analyzed the results of 304 solid tumor and hematologic malignancy specimens tested simultaneously by NGS and one or more targeted single-gene tests (EGFR, KRAS, BRAF, NPM1, FLT3, and JAK2). For samples that passed our validated tumor percentage and DNA quality and quantity thresholds, there was perfect concordance between NGS and targeted single-gene tests with the exception of two FLT3 internal tandem duplications that fell below the stringent pre-established reporting threshold but were readily detected by manual inspection. In addition, NGS identified clinically significant mutations not covered by single-gene tests. These findings confirm NGS as a reliable platform for routine clinical use when appropriate quality control metrics, such as tumor percentage and DNA quality cutoffs, are in place. Based on our findings, we suggest a simple workflow that should facilitate adoption of clinical oncologic NGS services at other institutions.     

Structural Variants Selected during Yak Domestication Inferred from Long-Read Whole-Genome Sequencing

Structural variants (SVs) represent an important genetic resource for both natural and artificial selection. Here we present a chromosome-scale reference genome for domestic yak (Bos grunniens) that has longer contigs and scaffolds (N50 44.72 and 114.39 Mb, respectively) than reported for any other ruminant genome. We further obtained long-read resequencing data for 6 wild and 23 domestic yaks and constructed a genetic SV map of 372,220 SVs that covers the geographic range of the yaks. The majority of the SVs contains repetitive sequences and several are in or near genes. By comparing SVs in domestic and wild yaks, we identified genes that are predominantly related to the nervous system, behavior, immunity, and reproduction and may have been targeted by artificial selection during yak domestication. These findings provide new insights in the domestication of animals living at high altitude and highlight the importance of SVs in animal domestication.     

The Accuracy,Feasibility and Challenges of Sequencing Short Tandem Repeats Using Next-Generation Sequencing Platforms

To date we have little knowledge of how accurate next-generation sequencing (NGS) technologies are in sequencing repetitive sequences beyond known limitations to accurately sequence homopolymers. Only a handful of previous reports have evaluated the potential of NGS for sequencing short tandem repeats (microsatellites) and no empirical study has compared and evaluated the performance of more than one NGS platform with the same dataset. Here we examined yeast microsatellite variants from both long-read (454-sequencing) and short-read (Illumina) NGS platforms and compared these to data derived through Sanger sequencing. In addition, we investigated any locus-specific biases and differences that might have resulted from variability in microsatellite repeat number, repeat motif or type of mutation. Out of 112 insertion/deletion variants identified among 45 microsatellite amplicons in our study, we found 87.5% agreement between the 454-platform and Sanger sequencing in frequency of variant detection after Benjamini-Hochberg correction for multiple tests. For a subset of 21 microsatellite amplicons derived from Illumina sequencing, the results of short-read platform were highly consistent with the other two platforms, with 100% agreement with 454-sequencing and 93.6% agreement with the Sanger method after Benjamini-Hochberg correction. We found that the microsatellite attributes copy number, repeat motif and type of mutation did not have a significant effect on differences seen between the sequencing platforms. We show that both long-read and short-read NGS platforms can be used to sequence short tandem repeats accurately, which makes it feasible to consider the use of these platforms in high-throughput genotyping. It appears the major requirement for achieving both high accuracy and rare variant detection in microsatellite genotyping is sufficient read depth coverage. This might be a challenge because each platform generates a consistent pattern of non-uniform sequence coverage, which, as our study suggests, may affect some types of tandem repeats more than others.     

Measuring KRAS Mutations in Circulating Tumor DNA by Droplet Digital PCR and Next-Generation Sequencing

Measuring total cell-free DNA (cfDNA) or cancer-specific mutations herein has presented as new tools in aiding the treatment of cancer patients. Studies show that total cfDNA bears prognostic value in metastatic colorectal cancer (mCRC) and that measuring cancer-specific mutations could supplement biopsies. However, limited information is available on the performance of different methods. Blood samples from 28 patients with mCRC and known KRAS mutation status were included. cfDNA was extracted and quantified with droplet digital polymerase chain reaction (ddPCR) measuring Beta-2 Microglobulin. KRAS mutation detection was performed using ddPCR (Bio-Rad) and next-generation sequencing (NGS, Ion Torrent PGM). Comparing KRAS mutation status in plasma and tissue revealed concordance rates of 79% and 89% for NGS and ddPCR. Strong correlation between the methods was observed. Most KRAS mutations were also detectable in 10-fold diluted samples using the ddPCR. We find that for detection of KRAS mutations in ctDNA ddPCR was superior to NGS both in analysis success rate and concordance to tissue. We further present results indicating that lower amount of plasma may be used for detection of KRAS mutations in mCRC.     

High Prevalence and Clinical Relevance of Genes Affected by Chromosomal Breaks in Colorectal Cancer

Background

Cancer is caused by somatic DNA alterations such as gene point mutations, DNA copy number aberrations (CNA) and structural variants (SVs). Genome-wide analyses of SVs in large sample series with well-documented clinical information are still scarce. Consequently, the impact of SVs on carcinogenesis and patient outcome remains poorly understood. This study aimed to perform a systematic analysis of genes that are affected by CNA-associated chromosomal breaks in colorectal cancer (CRC) and to determine the clinical relevance of recurrent breakpoint genes.

Methods

Primary CRC samples of patients with metastatic disease from CAIRO and CAIRO2 clinical trials were previously characterized by array-comparative genomic hybridization. These data were now used to determine the prevalence of CNA-associated chromosomal breaks within genes across 352 CRC samples. In addition, mutation status of the commonly affected APC, TP53, KRAS, PIK3CA, FBXW7, SMAD4, BRAF and NRAS genes was determined for 204 CRC samples by targeted massive parallel sequencing. Clinical relevance was assessed upon stratification of patients based on gene mutations and gene breakpoints that were observed in >3% of CRC cases.

Results

In total, 748 genes were identified that were recurrently affected by chromosomal breaks (FDR <0.1). MACROD2 was affected in 41% of CRC samples and another 169 genes showed breakpoints in >3% of cases, indicating that prevalence of gene breakpoints is comparable to the prevalence of well-known gene point mutations. Patient stratification based on gene breakpoints and point mutations revealed one CRC subtype with very poor prognosis.

Conclusions

We conclude that CNA-associated chromosomal breaks within genes represent a highly prevalent and clinically relevant subset of SVs in CRC.