Computational techniques for human genome resequencing using mated gapped reads

Analysis of 1386 epileptogenic brain lesions reveals association with DYRK1A and EGFR

Lesional focal epilepsy (LFE) is a common and severe seizure disorder caused by epileptogenic lesions, including malformations of cortical development (MCD) and low-grade epilepsy-associated tumors (LEAT). Understanding the genetic etiology of these lesions can inform medical and surgical treatment. We conducted a somatic variant enrichment mega-analysis in brain tissue from 1386 individuals who underwent epilepsy surgery, including 599 previously unpublished individuals with ultra-deep ( > 1600x) targeted panel sequencing. Here we confirm four known associations (BRAF, SLC35A2, MTOR, PTPN11), support eight associations without prior statistical support (FGFR1, PIK3CA, AKT3, NF1, PTEN, RHEB, KRAS, NRAS), and identify novel associations for two genes, DYRK1A and EGFR. Both novel genes show specific histopathological phenotypes, interact with LFE genes and pathways, and may represent promising candidates as biomarkers and potentially druggable targets.

Germline intergenic duplications at Xq26.1 underlie Bazex-Dupré-Christol basal cell carcinoma susceptibility syndrome

BACKGROUND

Bazex-Dupré-Christol syndrome (BDCS; MIM301845) is a rare X-linked dominant genodermatosis characterized by follicular atrophoderma, congenital hypotrichosis and multiple basal cell carcinomas (BCCs). Previous studies have linked BDCS to an 11·4-Mb interval on chromosome Xq25-q27.1. However, the genetic mechanism of BDCS remains an open question.

OBJECTIVES

To investigate the genetic aetiology and molecular mechanisms underlying BDCS.

METHODS

We ascertained multiple individuals from eight unrelated families affected with BDCS (F1-F8). Whole-exome (F1 and F2) and genome sequencing (F3) were performed to identify putative disease-causing variants within the linkage region. Array comparative genomic hybridization and quantitative polymerase chain reaction (PCR) were used to explore copy number variations, followed by long-range gap PCR and Sanger sequencing to amplify the duplication junctions and to define the head-tail junctions. Hi-C was performed on dermal fibroblasts from two affected individuals with BDCS and one control. Public datasets and tools were used to identify regulatory elements and transcription factor binding sites within the minimal duplicated region. Immunofluorescence was performed in hair follicles, BCCs and trichoepitheliomas from patients with BDCS and sporadic BCCs. The ACTRT1 variant c.547dup (p.Met183Asnfs*17), previously proposed to cause BDCS, was evaluated with t allele frequency calculator.

RESULTS

In eight families with BDCS, we identified overlapping 18-135-kb duplications (six inherited and two de novo) at Xq26.1, flanked by ARHGAP36 and IGSF1. Hi-C showed that the duplications did not affect the topologically associated domain, but may alter the interactions between flanking genes and putative enhancers located in the minimal duplicated region. We detected ARHGAP36 expression near the control hair follicular stem cell compartment, and found increased ARHGAP36 levels in hair follicles in telogen, in BCCs and in trichoepitheliomas from patients with BDCS. ARHGAP36 was also detected in sporadic BCCs from individuals without BDCS. Our modelling showed the predicted maximum tolerated minor allele frequency of ACTRT1 variants in control populations to be orders of magnitude higher than expected for a high-penetrant ultra-rare disorder, suggesting loss of function of ACTRT1 variants to be an unlikely cause for BDCS.

CONCLUSIONS

Noncoding Xq26.1 duplications cause BDCS. The BDCS duplications most likely lead to dysregulation of ARHGAP36. ARHGAP36 is a potential therapeutic target for both inherited and sporadic BCCs. What is already known about this topic? Bazex-Dupré-Christol syndrome (BDCS) is a rare X-linked basal cell carcinoma susceptibility syndrome linked to an 11·4-Mb interval on chromosome Xq25-q27.1. Loss-of-function variants in ACTRT1 and its regulatory elements were suggested to cause BDCS. What does this study add? BDCS is caused by small tandem noncoding intergenic duplications at chromosome Xq26.1. The Xq26.1 BDCS duplications likely dysregulate ARHGAP36, the flanking centromeric gene. ACTRT1 loss-of-function variants are unlikely to cause BDCS. What is the translational message? This study provides the basis for accurate genetic testing for BDCS, which will aid precise diagnosis and appropriate surveillance and clinical management. ARHGAP36 may be a novel therapeutic target for all forms of sporadic basal cell carcinomas.

Diverse tumorigenic consequences of human papillomavirus integration in primary oropharyngeal cancers

Human papillomavirus (HPV) causes 5% of all cancers and frequently integrates into host chromosomes. The HPV oncoproteins E6 and E7 are necessary but insufficient for cancer formation, indicating that additional secondary genetic events are required. Here, we investigate potential oncogenic impacts of virus integration. Analysis of 105 HPV-positive oropharyngeal cancers by whole-genome sequencing detects virus integration in 77%, revealing five statistically significant sites of recurrent integration near genes that regulate epithelial stem cell maintenance (i.e., SOX2, TP63, FGFR, MYC) and immune evasion (i.e., CD274). Genomic copy number hyperamplification is enriched 16-fold near HPV integrants, and the extent of focal host genomic instability increases with their local density. The frequency of genes expressed at extreme outlier levels is increased 86-fold within ±150 kb of integrants. Across 95% of tumors with integration, host gene transcription is disrupted via intragenic integrants, chimeric transcription, outlier expression, gene breaking, and/or de novo expression of noncoding or imprinted genes. We conclude that virus integration can contribute to carcinogenesis in a large majority of HPV-positive oropharyngeal cancers by inducing extensive disruption of host genome structure and gene expression.

Genome-wide mutational signatures revealed distinct developmental paths for human B cell lymphomas

Both somatic hypermutation (SHM) and class switch recombination (CSR) are initiated by activation-induced cytidine deaminase (AID). Dysregulation of these processes has been linked to B cell lymphomagenesis. Here we performed an in-depth analysis of diffuse large B cell lymphoma (DLBCL) and follicular lymphoma (FL) genomes. We characterized seven genomic mutational signatures, including two B cell tumor-specific signatures, one of which is novel and associated with aberrant SHM. We further identified two major mutational signatures (K1 and K2) of clustered mutations (kataegis) resulting from the activities of AID or error-prone DNA polymerase η, respectively. K1 was associated with the immunoglobulin (Ig) switch region mutations/translocations and the ABC subtype of DLBCL, whereas K2 was related to the Ig variable region mutations and the GCB subtype of DLBCL and FL. Similar patterns were also observed in chronic lymphocytic leukemia subtypes. Thus, alterations associated with aberrant CSR and SHM activities can be linked to distinct developmental paths for different subtypes of B cell lymphomas.

Whole-Genome Sequencing of Common Salivary Gland Carcinomas: Subtype-Restricted and Shared Genetic Alterations

PURPOSE

Salivary gland carcinomas (SGCs) are pathologically classified into several widely diverse subtypes, of which adenoid cystic carcinoma (ACC), mucoepidermoid carcinoma (MEC), and salivary duct carcinoma (SDC) are the most commonly encountered. A comparative genetic analysis of these subtypes provides detailed information on the genetic alterations that are associated with their tumorigenesis and may lead to the identification of biomarkers to guide tumor-specific clinical trials.

EXPERIMENTAL DESIGN

Whole-genome sequencing of 58 common SGCs (20 ACCs, 20 SDCs, and 18 MECs) was performed to catalog structural variations, copy number, rearrangements, and driver mutations. Data were bioinformatically analyzed and correlated with clinicopathologic parameters, and selected targets were validated.

RESULTS

Novel and recurrent type-specific and shared genetic alterations were identified within and among 3 subtypes. Mutually exclusive canonical fusion and nonfusion genomic alterations were identified in both ACC and MEC. In ACCs, loss of chromosome 12q was dominant in MYB or MYBL1 fusion-positive tumors and mutations of NOTCH pathway were more common in these fusion negatives. In MECs, CRTC1-MAML2 fusion-positive tumors showed frequent BAP1 mutation, and tumors lacking this fusion were enriched with LRFN1 mutation. SDCs displayed considerable genetic instability, lacked recurrent chromosomal rearrangements, and demonstrated nonoverlapping TP53 mutation and ERBB2 amplification in a subset of tumors. Limited genetic alterations, including focal amplifications of 8q21-q23, were shared by all subtypes and were associated with poor survival.

CONCLUSIONS

This study delineates type-specific and shared genetic alterations that are associated with early phenotypic commitment and the biologic progression of common SGCs. These alterations, upon validation, could serve as biomarkers in tumor-specific clinical trials.

Genomic and transcriptomic landscape of conjunctival melanoma

Conjunctival melanoma (CJM) is a rare but potentially lethal and highly-recurrent cancer of the eye. Similar to cutaneous melanoma (CM), it originates from melanocytes. Unlike CM, however, CJM is relatively poorly characterized from a genomic point of view. To fill this knowledge gap and gain insight into the genomic nature of CJM, we performed whole-exome (WES) or whole-genome sequencing (WGS) of tumor-normal tissue pairs in 14 affected individuals, as well as RNA sequencing in a subset of 11 tumor tissues. Our results show that, similarly to CM, CJM is also characterized by a very high mutation load, composed of approximately 500 somatic mutations in exonic regions. This, as well as the presence of a UV light-induced mutational signature, are clear signs of the role of sunlight in CJM tumorigenesis. In addition, the genomic classification of CM proposed by TCGA seems to be well-applicable to CJM, with the presence of four typical subclasses defined on the basis of the most frequently mutated genes: BRAF, NF1, RAS, and triple wild-type. In line with these results, transcriptomic analyses revealed similarities with CM as well, namely the presence of a transcriptomic subtype enriched for immune genes and a subtype enriched for genes associated with keratins and epithelial functions. Finally, in seven tumors we detected somatic mutations in ACSS3, a possible new candidate oncogene. Transfected conjunctival melanoma cells overexpressing mutant ACSS3 showed higher proliferative activity, supporting the direct involvement of this gene in the tumorigenesis of CJM. Altogether, our results provide the first unbiased and complete genomic and transcriptomic classification of CJM.

Conjunctival melanoma is an extremely rare form of cancer of the eye that arises from melanocytes–the cells producing the protective pigment melanin–in the outmost layer of the eye: the conjunctiva. This tissue, similarly to the skin, can also be exposed to UV light radiation from the sun. We investigated the genetic background of this rare form of cancer in samples from fourteen patients, by global DNA and RNA sequencing. Our results showed that conjunctival melanoma is genetically very similar to cutaneous melanoma. More precisely, in tumor DNA we detected signs of damage caused by UV light, as well as mutations in the genes BRAF, NF1 and NRAS/HRAS, previously described to be involved in cutaneous melanoma. Analysis of tumor gene expression also revealed similarities between these two types of cancer, some of which could be used as prognostic factors or as indicators of a patients’ response to therapy. In addition, we identified frequent somatic mutations in ACSS3, a gene not yet associated with either conjunctival or cutaneous melanoma, which represents a potential key player in oncogenesis of conjunctival melanoma.

A novel nicastrin mutation in a three-generation Dutch family with hidradenitis suppurativa: a search for functional significance

Background

Mutations in the γ‐secretase enzyme subunits have been described in multiple kindreds with familial hidradenitis suppurativa (HS).

Objective

In this study, we report a novel nicastrin (NCSTN) mutation causing HS in a Dutch family. We sought to explore the immunobiological function of NCSTN mutations using data of the Immunological Genome Project.

Methods

Blood samples of three affected and two unaffected family members were collected. Whole‐genome sequencing was performed using genomic DNA isolated from peripheral blood leucocytes. Sanger sequencing was done to confirm the causative NCSTN variant and the familial segregation. The microarray data set of the Immunological Genome Project was used for thorough dissection of the expression and function of wildtype NCSTN in the immune system.

Results

In a family consisting of 23 members, we found an autosomal dominant inheritance pattern of HS and detected a novel splice site mutation (c.1912_1915delCAGT) in the NCSTN gene resulting in a frameshift and subsequent premature stop. All affected individuals had HS lesions on non‐flexural and atypical locations. Wildtype NCSTN appears to be upregulated in myeloid cells like monocytes and macrophages, and in mesenchymal cells such as fibroblastic reticular cells and fibroblasts. In addition, within the 25 highest co‐expressed genes with NCSTN we identified CAPNS1,ARNT and PPARD.

Conclusion

This study reports the identification a novel NCSTN gene splice site mutation which causes familial HS. The associated immunobiological functions of NCSTN and its co‐expressed genes ARNT and PPARD link genetics to the most common environmental and metabolic HS risk factors which are smoking and obesity.

Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay

The diversity of disease presentations warrants one single assay for detection and delineation of various genomic disorders. Herein, we describe a gel-free and biotin-capture-free mate-pair method through coupling Controlled Polymerizations by Adapter-Ligation (CP-AL). We first demonstrated the feasibility and ease-of-use in monitoring DNA nick translation and primer extension by limiting the nucleotide input. By coupling these two controlled polymerizations by a reported non-conventional adapter-ligation reaction 3′ branch ligation, we evidenced that CP-AL significantly increased DNA circularization efficiency (by 4-fold) and was applicable for different sequencing methods but at a faction of current cost. Its advantages were further demonstrated by fully elimination of small-insert-contaminated (by 39.3-fold) with a ∼50% increment of physical coverage, and producing uniform genome/exome coverage and the lowest chimeric rate. It achieved single-nucleotide variants detection with sensitivity and specificity up to 97.3 and 99.7%, respectively, compared with data from small-insert libraries. In addition, this method can provide a comprehensive delineation of structural rearrangements, evidenced by a potential diagnosis in a patient with oligo-atheno-terato-spermia. Moreover, it enables accurate mutation identification by integration of genomic variants from different aberration types. Overall, it provides a potential single-integrated solution for detecting various genomic variants, facilitating a genetic diagnosis in human diseases.

Genome Sequencing Explores Complexity of Chromosomal Abnormalities in Recurrent Miscarriage

Recurrent miscarriage (RM) affects millions of couples globally, and half of them have no demonstrated etiology. Genome sequencing (GS) is an enhanced and novel cytogenetic tool to define the contribution of chromosomal abnormalities in human diseases. In this study we evaluated its utility in RM-affected couples. We performed low-pass GS retrospectively for 1,090 RM-affected couples, all of whom had routine chromosome analysis. A customized sequencing and interpretation pipeline was developed to identify chromosomal rearrangements and deletions/duplications with confirmation by fluorescence in situ hybridization, chromosomal microarray analysis, and PCR studies. Low-pass GS yielded results in 1,077 of 1,090 couples (98.8%) and detected 127 chromosomal abnormalities in 11.7% (126/1,077) of couples; both members of one couple were identified with inversions. Of the 126 couples, 39.7% (50/126) had received former diagnostic results by karyotyping characteristic of normal human male or female karyotypes. Low-pass GS revealed additional chromosomal abnormalities in 50 (4.0%) couples, including eight with balanced translocations and 42 inversions. Follow-up studies of these couples showed a higher miscarriage/fetal-anomaly rate of 5/10 (50%) compared to 21/93 (22.6%) in couples with normal GS, resulting in a relative risk of 2.2 (95% confidence interval, 1.1 to 4.6). In these couples, this protocol significantly increased the diagnostic yield of chromosomal abnormalities per couple (11.7%) in comparison to chromosome analysis (8.0%, chi-square test p = 0.000751). In summary, low-pass GS identified underlying chromosomal aberrations in 1 in 9 RM-affected couples, enabling identification of a subgroup of couples with increased risk of subsequent miscarriage who would benefit from a personalized intervention.

Leveraging known genomic variants to improve detection of variants, especially close-by Indels

Motivation

The detection of genomic variants has great significance in genomics, bioinformatics, biomedical research and its applications. However, despite a lot of effort, Indels and structural variants are still under-characterized compared to SNPs. Current approaches based on next-generation sequencing data usually require large numbers of reads (high coverage) to be able to detect such types of variants accurately. However Indels, especially those close to each other, are still hard to detect accurately.

Results

We introduce a novel approach that leverages known variant information, e.g. provided by dbSNP, dbVar, ExAC or the 1000 Genomes Project, to improve sensitivity of detecting variants, especially close-by Indels. In our approach, the standard reference genome and the known variants are combined to build a meta-reference, which is expected to be probabilistically closer to the subject genomes than the standard reference. An alignment algorithm, which can take into account known variant information, is developed to accurately align reads to the meta-reference. This strategy resulted in accurate alignment and variant calling even with low coverage data. We showed that compared to popular methods such as GATK and SAMtools, our method significantly improves the sensitivity of detecting variants, especially Indels that are close to each other. In particular, our method was able to call these close-by Indels at a 15-20% higher sensitivity than other methods at low coverage, and still get 1-5% higher sensitivity at high coverage, at competitive precision. These results were validated using simulated data with variant profiles extracted from the 1000 Genomes Project data, and real data from the Illumina Platinum Genomes Project and ExAC database. Our finding suggests that by incorporating known variant information in an appropriate manner, sensitive variant calling is possible at a low cost.

Availability and implementation

Implementation can be found in our public code repository https://github.com/namsyvo/IVC.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genome sequencing uncovers phenocopies in primary progressive multiple sclerosis

OBJECTIVE

Primary progressive multiple sclerosis (PPMS) causes accumulation of neurological disability from disease onset without clinical attacks typical of relapsing multiple sclerosis (RMS). However, whether genetic variation influences the disease course remains unclear. We aimed to determine whether mutations causative of neurological disorders that share features with multiple sclerosis (MS) contribute to risk for developing PPMS.

METHODS

We examined whole-genome sequencing (WGS) data from 38 PPMS and 81 healthy subjects of European ancestry. We selected pathogenic variants exclusively found in PPMS patients that cause monogenic neurological disorders and performed two rounds of replication genotyping in 746 PPMS, 3,049 RMS, and 1,000 healthy subjects. To refine our findings, we examined the burden of rare, potentially pathogenic mutations in 41 genes that cause hereditary spastic paraplegias (HSPs) in PPMS (n = 314), secondary progressive multiple sclerosis (SPMS; n = 587), RMS (n = 2,248), and healthy subjects (n = 987) genotyped using the MS replication chip.

RESULTS

WGS and replication studies identified three pathogenic variants in PPMS patients that cause neurological disorders sharing features with MS: KIF5A p.Ala361Val in spastic paraplegia 10; MLC1 p.Pro92Ser in megalencephalic leukodystrophy with subcortical cysts, and REEP1 c.606 + 43G>T in Spastic Paraplegia 31. Moreover, we detected a significant enrichment of HSP-related mutations in PPMS patients compared to controls (risk ratio [RR] = 1.95; 95% confidence interval [CI], 1.27-2.98; p = 0.002), as well as in SPMS patients compared to controls (RR = 1.57; 95% CI, 1.18-2.10; p = 0.002). Importantly, this enrichment was not detected in RMS.

INTERPRETATION

This study provides evidence to support the hypothesis that rare Mendelian genetic variants contribute to the risk for developing progressive forms of MS. Ann Neurol 2018;83:51-63.

The landscape of somatic mutation in sporadic Chinese colorectal cancer

Colorectal cancer is the fifth prevalent cancer in China. Nevertheless, a large-scale characterization of Chinese colorectal cancer mutation spectrum has not been carried out. In this study, we have performed whole exome-sequencing analysis of 98 patients’ tumor samples with matched pairs of normal colon tissues using Illumina and Complete Genomics high-throughput sequencing platforms. Canonical CRC somatic gene mutations with high prevalence (>10%) have been verified, including TP53, APC, KRAS, SMAD4, FBXW7 and PIK3CA. PEG3 is identified as a novel frequently mutated gene (10.6%). APC and Wnt signaling exhibit significantly lower mutation frequencies than those in TCGA data. Analysis with clinical characteristics indicates that APC gene and Wnt signaling display lower mutation rate in lymph node positive cancer than negative ones, which are not observed in TCGA data. APC gene and Wnt signaling are considered as the key molecule and pathway for colorectal cancer initiation, and these findings greatly undermine their importance in tumor progression for Chinese patients. Taken together, the application of next-generation sequencing has led to the determination of novel somatic mutations and alternative disease mechanisms in colorectal cancer progression, which may be useful for understanding disease mechanism and personalizing treatment for Chinese patients.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

Advanced Whole-Genome Sequencing and Analysis of Fetal Genomes from Amniotic Fluid

BACKGROUND

Amniocentesis is a common procedure, the primary purpose of which is to collect cells from the fetus to allow testing for abnormal chromosomes, altered chromosomal copy number, or a small number of genes that have small single- to multibase defects. Here we demonstrate the feasibility of generating an accurate whole-genome sequence of a fetus from either the cellular or cell-free DNA (cfDNA) of an amniotic sample.

METHODS

cfDNA and DNA isolated from the cell pellet of 31 amniocenteses were sequenced to approximately 50× genome coverage by use of the Complete Genomics nanoarray platform. In a subset of the samples, long fragment read libraries were generated from DNA isolated from cells and sequenced to approximately 100× genome coverage.

RESULTS

Concordance of variant calls between the 2 DNA sources and with parental libraries was >96%. Two fetal genomes were found to harbor potentially detrimental variants in chromodomain helicase DNA binding protein 8 (CHD8) and LDL receptor-related protein 1 (LRP1), variations of which have been associated with autism spectrum disorder and keratosis pilaris atrophicans, respectively. We also discovered drug sensitivities and carrier information of fetuses for a variety of diseases.

CONCLUSIONS

We were able to elucidate the complete genome sequence of 31 fetuses from amniotic fluid and demonstrate that the cfDNA or DNA from the cell pellet can be analyzed with little difference in quality. We believe that current technologies could analyze this material in a highly accurate and complete manner and that analyses like these should be considered for addition to current amniocentesis procedures.

Whole-genome sequencing reveals principles of brain retrotransposition in neurodevelopmental disorders

Neural progenitor cells undergo somatic retrotransposition events, mainly involving L1 elements, which can be potentially deleterious. Here, we analyze the whole genomes of 20 brain samples and 80 non-brain samples, and characterized the retrotransposition landscape of patients affected by a variety of neurodevelopmental disorders including Rett syndrome, tuberous sclerosis, ataxia-telangiectasia and autism. We report that the number of retrotranspositions in brain tissues is higher than that observed in non-brain samples and even higher in pathologic vs normal brains. The majority of somatic brain retrotransposons integrate into pre-existing repetitive elements, preferentially A/T rich L1 sequences, resulting in nested insertions. Our findings document the fingerprints of encoded endonuclease independent mechanisms in the majority of L1 brain insertion events. The insertions are "non-classical" in that they are truncated at both ends, integrate in the same orientation as the host element, and their target sequences are enriched with a CCATT motif in contrast to the classical endonuclease motif of most other retrotranspositions. We show that L1Hs elements integrate preferentially into genes associated with neural functions and diseases. We propose that pre-existing retrotransposons act as "lightning rods" for novel insertions, which may give fine modulation of gene expression while safeguarding from deleterious events. Overwhelmingly uncontrolled retrotransposition may breach this safeguard mechanism and increase the risk of harmful mutagenesis in neurodevelopmental disorders.

VarCards: an integrated genetic and clinical database for coding variants in the human genome

A growing number of genomic tools and databases were developed to facilitate the interpretation of genomic variants, particularly in coding regions. However, these tools are separately available in different online websites or databases, making it challenging for general clinicians, geneticists and biologists to obtain the first-hand information regarding some particular variants and genes of interest. Starting with coding regions and splice sties, we artificially generated all possible single nucleotide variants (n = 110 154 363) and cataloged all reported insertion and deletions (n = 1 223 370). We then annotated these variants with respect to functional consequences from more than 60 genomic data sources to develop a database, named VarCards (http://varcards.biols.ac.cn/), by which users can conveniently search, browse and annotate the variant- and gene-level implications of given variants, including the following information: (i) functional effects; (ii) functional consequences through different in silico algorithms; (iii) allele frequencies in different populations; (iv) disease- and phenotype-related knowledge; (v) general meaningful gene-level information; and (vi) drug-gene interactions. As a case study, we successfully employed VarCards in interpretation of de novo mutations in autism spectrum disorders. In conclusion, VarCards provides an intuitive interface of necessary information for researchers to prioritize candidate variations and genes.

Assessment of the incorporation of CNV surveillance into gene panel next-generation sequencing testing for inherited retinal diseases

BACKGROUND

Diagnostic use of gene panel next-generation sequencing (NGS) techniques is commonplace for individuals with inherited retinal dystrophies (IRDs), a highly genetically heterogeneous group of disorders. However, these techniques have often failed to capture the complete spectrum of genomic variation causing IRD, including CNVs. This study assessed the applicability of introducing CNV surveillance into first-tier diagnostic gene panel NGS services for IRD.

METHODS

Three read-depth algorithms were applied to gene panel NGS data sets for 550 referred individuals, and informatics strategies used for quality assurance and CNV filtering. CNV events were confirmed and reported to referring clinicians through an accredited diagnostic laboratory.

RESULTS

We confirmed the presence of 33 deletions and 11 duplications, determining these findings to contribute to the confirmed or provisional molecular diagnosis of IRD for 25 individuals. We show that at least 7% of individuals referred for diagnostic testing for IRD have a CNV within genes relevant to their clinical diagnosis, and determined a positive predictive value of 79% for the employed CNV filtering techniques.

CONCLUSION

Incorporation of CNV analysis increases diagnostic yield of gene panel NGS diagnostic tests for IRD, increases clarity in diagnostic reporting and expands the spectrum of known disease-causing mutations.

Analysis and Annotation of Whole-Genome or Whole-Exome Sequencing Derived Variants for Clinical Diagnosis

Over the last 10 years, next-generation sequencing (NGS) has transformed genomic research through substantial advances in technology and reduction in the cost of sequencing, and also in the systems required for analysis of these large volumes of data. This technology is now being used as a standard molecular diagnostic test in some clinical settings. The advances in sequencing have come so rapidly that the major bottleneck in identification of causal variants is no longer the sequencing or analysis (given access to appropriate tools), but rather clinical interpretation. Interpretation of genetic findings in a complex and ever changing clinical setting is scarcely a new challenge, but the task is increasingly complex in clinical genome-wide sequencing given the dramatic increase in dataset size and complexity. This increase requires application of appropriate interpretation tools, as well as development and application of appropriate methodologies and standard procedures. This unit provides an overview of these items. Specific challenges related to implementation of genome-wide sequencing in a clinical setting are discussed. © 2017 by John Wiley & Sons, Inc.

Tools for annotation and comparison of structural variation

The impact of structural variants (SVs) on a variety of organisms and diseases like cancer has become increasingly evident. Methods for SV detection when studying genomic differences across cells, individuals or populations are being actively developed. Currently, just a few methods are available to compare different SVs callsets, and no specialized methods are available to annotate SVs that account for the unique characteristics of these variant types. Here, we introduce SURVIVOR_ant, a tool that compares types and breakpoints for candidate SVs from different callsets and enables fast comparison of SVs to genomic features such as genes and repetitive regions, as well as to previously established SV datasets such as from the 1000 Genomes Project. As proof of concept we compared 16 SV callsets generated by different SV calling methods on a single genome, the Genome in a Bottle sample HG002 (Ashkenazi son), and annotated the SVs with gene annotations, 1000 Genomes Project SV calls, and four different types of repetitive regions. Computation time to annotate 134,528 SVs with 33,954 of annotations was 22 seconds on a laptop.

Mutations of conserved non-coding elements of PITX2 in patients with ocular dysgenesis and developmental glaucoma

Mutations in FOXC1 and PITX2 constitute the most common causes of ocular anterior segment dysgenesis (ASD), and confer a high risk for secondary glaucoma. The genetic causes underlying ASD in approximately half of patients remain unknown, despite many of them being screened by whole exome sequencing. Here, we performed whole genome sequencing on DNA from two affected individuals from a family with dominantly inherited ASD and glaucoma to identify a 748-kb deletion in a gene desert that contains conserved putative PITX2 regulatory elements. We used CRISPR/Cas9 to delete the orthologous region in zebrafish in order to test the pathogenicity of this structural variant. Deletion in zebrafish reduced pitx2 expression during development and resulted in shallow anterior chambers. We screened additional patients for copy number variation of the putative regulatory elements and found an overlapping deletion in a second family and in a potentially-ancestrally-related index patient with ASD and glaucoma. These data suggest that mutations affecting conserved non-coding elements of PITX2 may constitute an important class of mutations in patients with ASD for whom the molecular cause of their disease have not yet been identified. Improved functional annotation of the human genome and transition to sequencing of patient genomes instead of exomes will be required before the magnitude of this class of mutations is fully understood.

Evidence for genetic association between chromosome 1q loci and predisposition to colorectal neoplasia

BACKGROUND

A substantial fraction of familial colorectal cancer (CRC) and polyposis heritability remains unexplained. This study aimed to identify predisposing loci in patients with these disorders.

METHODS

Homozygosity mapping was performed using 222 563 SNPs in 302 index patients with various colorectal neoplasms and 3367 controls. Linkage analysis, exome and whole-genome sequencing were performed in a family affected by microsatellite stable CRCs. Candidate variants were genotyped in 10 554 cases and 21 480 controls. Gene expression was assessed at the mRNA and protein level.

RESULTS

Homozygosity mapping revealed a disease-associated region at 1q32.3 which was part of the linkage region 1q32.2-42.2 identified in the CRC family. This includes a region previously associated with risk of CRC. Sequencing identified the p.Asp1432Glu variant in the MIA3 gene (known as TANGO1 or TANGO) and 472 additional rare, shared variants within the linkage region. In both cases and controls the population frequency was 0.02% for this MIA3 variant. The MIA3 mutant allele showed predominant mRNA expression in normal, cancer and precancerous tissues. Furthermore, immunohistochemistry revealed increased expression of MIA3 in adenomatous tissues.

CONCLUSIONS

Taken together, our two independent strategies associate genetic variations in chromosome 1q loci and predisposition to familial CRC and polyps, which warrants further investigation.

A novel MERTK mutation causing retinitis pigmentosa

PURPOSE

Retinitis pigmentosa (RP) is a genetically heterogeneous inherited retinal dystrophy. To date, over 80 genes have been implicated in RP. However, the disease demonstrates significant locus and allelic heterogeneity not entirely captured by current testing platforms. The purpose of the present study was to characterize the underlying mutation in a patient with RP without a molecular diagnosis after initial genetic testing.

METHODS

Whole-exome sequencing of the affected proband was performed. Candidate gene mutations were selected based on adherence to expected genetic inheritance pattern and predicted pathogenicity. Sanger sequencing of MERTK was completed on the patient's unaffected mother, affected brother, and unaffected sister to determine genetic phase.

RESULTS

Eight sequence variants were identified in the proband in known RP-associated genes. Sequence analysis revealed that the proband was a compound heterozygote with two independent mutations in MERTK, a novel nonsense mutation (c.2179C > T) and a previously reported missense variant (c.2530C > T). The proband's affected brother also had both mutations. Predicted phase was confirmed in unaffected family members.

CONCLUSION

Our study identifies a novel nonsense mutation in MERTK in a family with RP and no prior molecular diagnosis. The present study also demonstrates the clinical value of exome sequencing in determining the genetic basis of Mendelian diseases when standard genetic testing is unsuccessful.

A population-specific reference panel empowers genetic studies of Anabaptist populations

Genotype imputation is a powerful strategy for achieving the large sample sizes required for identification of variants underlying complex phenotypes, but imputation of rare variants remains problematic. Genetically isolated populations offer one solution, however population-specific reference panels are needed to assure optimal imputation accuracy and allele frequency estimation. Here we report the Anabaptist Genome Reference Panel (AGRP), the first whole-genome catalogue of variants and phased haplotypes in people of Amish and Mennonite ancestry. Based on high-depth whole-genome sequence (WGS) from 265 individuals, the AGRP contains >12 M high-confidence single nucleotide variants and short indels, of which ~12.5% are novel. These Anabaptist-specific variants were more deleterious than variants with comparable frequencies observed in the 1000 Genomes panel. About 43,000 variants showed enriched allele frequencies in AGRP, consistent with drift. When combined with the 1000 Genomes Project reference panel, the AGRP substantially improved imputation, especially for rarer variants. The AGRP is freely available to researchers through an imputation server.

Genome sequencing as a platform for pharmacogenetic genotyping: a pediatric cohort study

Whole-genome sequencing and whole-exome sequencing have proven valuable for diagnosing inherited diseases, particularly in children. However, usage of sequencing data as a pharmacogenetic screening tool to ensure medication safety and effectiveness remains to be explored. Sixty-seven variants in 19 genes with known effects on drug response were compared between genome sequencing and targeted genotyping data for coverage and concordance in 98 pediatric patients. We used targeted genotyping data as a benchmark to assess accuracy of variant calling, and to identify copy number variations of the CYP2D6 gene. We then predicted clinical impact of these variants on drug therapy. We find genotype concordance across those panels to be > 97%. Concordance of CYP2D6 predicted phenotype between estimates of whole-genome sequencing and targeted genotyping panel were 90%; a result from a lower coverage depth or variant calling difficulties in our whole-genome sequencing data when copy number variation and/or the CYP2D6*4 haplotype were present. Importantly, 95 children had at least one clinically actionable pharmacogenetic variant. Diagnostic genomic sequencing data can be used for pre-emptive pharmacogenetic screening. However, concordance between genome-wide sequencing and target genotyping needs to be characterized for each of the pharmacologically important genes.

Genome sequencing, in addition to helping diagnose disease, can inform medication decisions and improve drug safety for children. Ronald Cohn, Shinya Ito and colleagues at the Hospital for Sick Children in Toronto, Canada, studied a cohort of 98 pediatric patients who had undergone whole genome sequencing to help diagnose their unexplained congenital malformations or neurodevelopmental disorders. The researchers looked for 67 DNA variants found in 19 genes with known effects on drug responses. They used targeted genotyping to assess the accuracy of the sequence data. Sequencing proved to be more than 99% accurate for all but one of the pharmacologically important genes, showing the power of diagnostic genomic sequencing to identify DNA variants in children that affect medication safety and effectiveness. However, the accuracy of the method may need to be validated for each relevant gene.

Principles and Recommendations for Standardizing the Use of the Next-Generation Sequencing Variant File in Clinical Settings

A national workgroup convened by the Centers for Disease Control and Prevention identified principles and made recommendations for standardizing the description of sequence data contained within the variant file generated during the course of clinical next-generation sequence analysis for diagnosing human heritable conditions. The specifications for variant files were initially developed to be flexible with regard to content representation to support a variety of research applications. This flexibility permits variation with regard to how sequence findings are described and this depends, in part, on the conventions used. For clinical laboratory testing, this poses a problem because these differences can compromise the capability to compare sequence findings among laboratories to confirm results and to query databases to identify clinically relevant variants. To provide for a more consistent representation of sequence findings described within variant files, the workgroup made several recommendations that considered alignment to a common reference sequence, variant caller settings, use of genomic coordinates, and gene and variant naming conventions. These recommendations were considered with regard to the existing variant file specifications presently used in the clinical setting. Adoption of these recommendations is anticipated to reduce the potential for ambiguity in describing sequence findings and facilitate the sharing of genomic data among clinical laboratories and other entities.

Validation of copy number variation analysis for next-generation sequencing diagnostics

Although a common cause of disease, copy number variants (CNVs) have not routinely been identified from next-generation sequencing (NGS) data in a clinical context. This study aimed to examine the sensitivity and specificity of a widely used software package, ExomeDepth, to identify CNVs from targeted NGS data sets. We benchmarked the accuracy of CNV detection using ExomeDepth v1.1.6 applied to targeted NGS data sets by comparison to CNV events detected through whole-genome sequencing for 25 individuals and determined the sensitivity and specificity of ExomeDepth applied to these targeted NGS data sets to be 100% and 99.8%, respectively. To define quality assurance metrics for CNV surveillance through ExomeDepth, we undertook simulation of single-exon (n=1000) and multiple-exon heterozygous deletion events (n=1749), determining a sensitivity of 97% (n=2749). We identified that the extent of sequencing coverage, the inter- and intra-sample variability in the depth of sequencing coverage and the composition of analysis regions are all important determinants of successful CNV surveillance through ExomeDepth. We then applied these quality assurance metrics during CNV surveillance for 140 individuals across 12 distinct clinical areas, encompassing over 500 potential rare disease diagnoses. All 140 individuals lacked molecular diagnoses after routine clinical NGS testing, and by application of ExomeDepth, we identified 17 CNVs contributing to the cause of a Mendelian disorder. Our findings support the integration of CNV detection using ExomeDepth v1.1.6 with routine targeted NGS diagnostic services for Mendelian disorders. Implementation of this strategy increases diagnostic yields and enhances clinical care.

Differences between the genomes of lymphoblastoid cell lines and blood-derived samples

Lymphoblastoid cell lines (LCLs) represent a convenient research tool for expanding the amount of biologic material available from an individual. LCLs are commonly used as reference materials, most notably from the Genome in a Bottle Consortium. However, the question remains how faithfully LCL-derived genome assemblies represent the germline genome of the donor individual as compared to the genome assemblies derived from peripheral blood mononuclear cells. We present an in-depth comparison of a large collection of LCL- and peripheral blood mononuclear cell-derived genomes in terms of distributions of coverage and copy number alterations. We found significant differences in the depth of coverage and copy number calls, which may be driven by differential replication timing. Importantly, these copy number changes preferentially affect regions closer to genes and with higher GC content. This suggests that genomic studies based on LCLs may display locus-specific biases, and that conclusions based on analysis of depth of coverage and copy number variation may require further scrutiny.

Whole genome sequencing of a dizygotic twin suggests a role for the serotonin receptor HTR7 in autism spectrum disorder

Whole genome sequencing of a severely affected dizygotic twin with an autism spectrum disorder and intellectual disability revealed a compound heterozygous mutation in the HTR7 gene as the only variation not detected in control databases. Each parent carries one allele of the mutation, which is not present in an unaffected stepsister. The HTR7 gene encodes the 5-HT₇ serotonin receptor that is involved in brain development, synaptic transmission, and plasticity. The paternally inherited p.W60C variant is situated at an evolutionary conserved nucleotide and predicted damaging by Polyphen2. A mutation akin to the maternally inherited pV286I mutation has been reported to significantly affect the binding characteristics of the receptor. Therefore, the observed sequence alterations provide a first suggestive link between a genetic abnormality in the HTR7 gene and a neurodevelopmental disorder. © 2016 Wiley Periodicals, Inc.

Comprehensive Genetic Landscape of Uveal Melanoma by Whole-Genome Sequencing

Uveal melanoma (UM) is a rare intraocular tumor that, similar to cutaneous melanoma, originates from melanocytes. To gain insights into its genetics, we performed whole-genome sequencing at very deep coverage of tumor-control pairs in 33 samples (24 primary and 9 metastases). Genome-wide, the number of coding mutations was rather low (only 17 variants per tumor on average; range 7-28), thus radically different from cutaneous melanoma, where hundreds of exonic DNA insults are usually detected. Furthermore, no UV light-induced mutational signature was identified. Recurrent coding mutations were found in the known UM drivers GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. Other genes, i.e., TP53BP1, CSMD1, TTC28, DLK2, and KTN1, were also found to harbor somatic mutations in more than one individual, possibly indicating a previously undescribed association with UM pathogenesis. De novo assembly of unmatched reads from non-coding DNA revealed peculiar copy-number variations defining specific UM subtypes, which in turn could be associated with metastatic transformation. Mutational-driven comparison with other tumor types showed that UM is very similar to pediatric tumors, characterized by very few somatic insults and, possibly, important epigenetic changes. Through the analysis of whole-genome sequencing data, our findings shed new light on the molecular genetics of uveal melanoma, delineating it as an atypical tumor of the adult for which somatic events other than mutations in exonic DNA shape its genetic landscape and define its metastatic potential.

The whole genome sequences and experimentally phased haplotypes of over 100 personal genomes

Background

Since the completion of the Human Genome Project in 2003, it is estimated that more than 200,000 individual whole human genomes have been sequenced. A stunning accomplishment in such a short period of time. However, most of these were sequenced without experimental haplotype data and are therefore missing an important aspect of genome biology. In addition, much of the genomic data is not available to the public and lacks phenotypic information.

Findings

As part of the Personal Genome Project, blood samples from 184 participants were collected and processed using Complete Genomics’ Long Fragment Read technology. Here, we present the experimental whole genome haplotyping and sequencing of these samples to an average read coverage depth of 100X. This is approximately three-fold higher than the read coverage applied to most whole human genome assemblies and ensures the highest quality results. Currently, 114 genomes from this dataset are freely available in the GigaDB repository and are associated with rich phenotypic data; the remaining 70 should be added in the near future as they are approved through the PGP data release process. For reproducibility analyses, 20 genomes were sequenced at least twice using independent LFR barcoded libraries. Seven genomes were also sequenced using Complete Genomics’ standard non-barcoded library process. In addition, we report 2.6 million high-quality, rare variants not previously identified in the Single Nucleotide Polymorphisms database or the 1000 Genomes Project Phase 3 data.

Conclusions

These genomes represent a unique source of haplotype and phenotype data for the scientific community and should help to expand our understanding of human genome evolution and function.

Electronic supplementary material

The online version of this article (doi:10.1186/s13742-016-0148-z) contains supplementary material, which is available to authorized users.

Premalignant SOX2 overexpression in the fallopian tubes of ovarian cancer patients: Discovery and validation studies

Current screening methods for ovarian cancer can only detect advanced disease. Earlier detection has proved difficult because the molecular precursors involved in the natural history of the disease are unknown. To identify early driver mutations in ovarian cancer cells, we used dense whole genome sequencing of micrometastases and microscopic residual disease collected at three time points over three years from a single patient during treatment for high-grade serous ovarian cancer (HGSOC). The functional and clinical significance of the identified mutations was examined using a combination of population-based whole genome sequencing, targeted deep sequencing, multi-center analysis of protein expression, loss of function experiments in an in-vivo reporter assay and mammalian models, and gain of function experiments in primary cultured fallopian tube epithelial (FTE) cells. We identified frequent mutations involving a 40kb distal repressor region for the key stem cell differentiation gene SOX2. In the apparently normal FTE, the region was also mutated. This was associated with a profound increase in SOX2 expression (p<2(-16)), which was not found in patients without cancer (n=108). Importantly, we show that SOX2 overexpression in FTE is nearly ubiquitous in patients with HGSOCs (n=100), and common in BRCA1-BRCA2 mutation carriers (n=71) who underwent prophylactic salpingo-oophorectomy. We propose that the finding of SOX2 overexpression in FTE could be exploited to develop biomarkers for detecting disease at a premalignant stage, which would reduce mortality from this devastating disease.

Parent-of-origin-specific signatures of de novo mutations

De novo mutations (DNMs) originating in gametogenesis are an important source of genetic variation. We use a data set of 7,216 autosomal DNMs with resolved parent of origin from whole-genome sequencing of 816 parent-offspring trios to investigate differences between maternally and paternally derived DNMs and study the underlying mutational mechanisms. Our results show that the number of DNMs in offspring increases not only with paternal age, but also with maternal age, and that some genome regions show enrichment for maternally derived DNMs. We identify parent-of-origin-specific mutation signatures that become more pronounced with increased parental age, pointing to different mutational mechanisms in spermatogenesis and oogenesis. Moreover, we find DNMs that are spatially clustered to have a unique mutational signature with no significant differences between parental alleles, suggesting a different mutational mechanism. Our findings provide insights into the molecular mechanisms that underlie mutagenesis and are relevant to disease and evolution in humans.

Extensive sequencing of seven human genomes to characterize benchmark reference materials

The Genome in a Bottle Consortium, hosted by the National Institute of Standards and Technology (NIST) is creating reference materials and data for human genome sequencing, as well as methods for genome comparison and benchmarking. Here, we describe a large, diverse set of sequencing data for seven human genomes; five are current or candidate NIST Reference Materials. The pilot genome, NA12878, has been released as NIST RM 8398. We also describe data from two Personal Genome Project trios, one of Ashkenazim Jewish ancestry and one of Chinese ancestry. The data come from 12 technologies: BioNano Genomics, Complete Genomics paired-end and LFR, Ion Proton exome, Oxford Nanopore, Pacific Biosciences, SOLiD, 10X Genomics GemCode WGS, and Illumina exome and WGS paired-end, mate-pair, and synthetic long reads. Cell lines, DNA, and data from these individuals are publicly available. Therefore, we expect these data to be useful for revealing novel information about the human genome and improving sequencing technologies, SNP, indel, and structural variant calling, and de novo assembly.

Genomic Characterization of Esophageal Squamous Cell Carcinoma Reveals Critical Genes Underlying Tumorigenesis and Poor Prognosis

The genetic mechanisms underlying the poor prognosis of esophageal squamous cell carcinoma (ESCC) are not well understood. Here, we report somatic mutations found in ESCC from sequencing 10 whole-genome and 57 whole-exome matched tumor-normal sample pairs. Among the identified genes, we characterized mutations in VANGL1 and showed that they accelerated cell growth in vitro. We also found that five other genes, including three coding genes (SHANK2, MYBL2, FADD) and two non-coding genes (miR-4707-5p, PCAT1), were involved in somatic copy-number alterations (SCNAs) or structural variants (SVs). A survival analysis based on the expression profiles of 321 individuals with ESCC indicated that these genes were significantly associated with poorer survival. Subsequently, we performed functional studies, which showed that miR-4707-5p and MYBL2 promoted proliferation and metastasis. Together, our results shed light on somatic mutations and genomic events that contribute to ESCC tumorigenesis and prognosis and might suggest therapeutic targets.

Mutations in C8ORF37 cause Bardet Biedl syndrome (BBS21)

Bardet Biedl syndrome (BBS) is a multisystem genetically heterogeneous ciliopathy that most commonly leads to obesity, photoreceptor degeneration, digit anomalies, genito-urinary abnormalities, as well as cognitive impairment with autism, among other features. Sequencing of a DNA sample from a 17-year-old female affected with BBS did not identify any mutation in the known BBS genes. Whole-genome sequencing identified a novel loss-of-function disease-causing homozygous mutation (K102*) in C8ORF37, a gene coding for a cilia protein. The proband was overweight (body mass index 29.1) with a slowly progressive rod-cone dystrophy, a mild learning difficulty, high myopia, three limb post-axial polydactyly, horseshoe kidney, abnormally positioned uterus and elevated liver enzymes. Mutations in C8ORF37 were previously associated with severe autosomal recessive retinal dystrophies (retinitis pigmentosa RP64 and cone-rod dystrophy CORD16) but not BBS. To elucidate the functional role of C8ORF37 in a vertebrate system, we performed gene knockdown in Danio rerio and assessed the cardinal features of BBS and visual function. Knockdown of c8orf37 resulted in impaired visual behavior and BBS-related phenotypes, specifically, defects in the formation of Kupffer's vesicle and delays in retrograde transport. Specificity of these phenotypes to BBS knockdown was shown with rescue experiments. Over-expression of human missense mutations in zebrafish also resulted in impaired visual behavior and BBS-related phenotypes. This is the first functional validation and association of C8ORF37 mutations with the BBS phenotype, which identifies BBS21. The zebrafish studies hereby show that C8ORF37 variants underlie clinically diagnosed BBS-related phenotypes as well as isolated retinal degeneration.

Recurrent viral infections associated with a homozygous CORO1A mutation that disrupts oligomerization and cytoskeletal association

BACKGROUND

Coronin-1A (CORO1A) is a regulator of actin dynamics important for T-cell homeostasis. CORO1A deficiency causes T(-)B(+) natural killer-positive severe combined immunodeficiency or T-cell lymphopenia with severe viral infections. However, because all known human mutations in CORO1A abrogate protein expression, the role of the protein's functional domains in host immunity is unknown.

OBJECTIVE

We sought to identify the cause of the primary immunodeficiency in 2 young adult siblings with a history of disseminated varicella, cutaneous warts, and CD4(+) T-cell lymphopenia.

METHODS

We performed immunologic, genetic, and biochemical studies in the patients, family members, and healthy control subjects.

RESULTS

Both patients had CD4(+) T-cell lymphopenia and decreased lymphocyte proliferation to mitogens. IgG, IgM, IgA, and specific antibody responses were normal. Whole-genome sequencing identified a homozygous frameshift mutation in CORO1A disrupting the last 2 C-terminal domains by replacing 61 amino acids with a novel 91-amino-acid sequence. The CORO1A(S401fs) mutant was expressed in the patients' lymphocytes at a level comparable with that of wild-type CORO1A in normal lymphocytes but did not oligomerize and had impaired cytoskeletal association. CORO1A(S401fs) was associated with increased filamentous actin accumulation in T cells, severely defective thymic output, and impaired T-cell survival but normal calcium flux and cytotoxicity, demonstrating the importance of CORO1A oligomerization and subcellular localization in T-cell homeostasis.

CONCLUSIONS

We describe a truncating mutation in CORO1A that permits protein expression and survival into young adulthood. Our studies demonstrate the importance of intact CORO1A C-terminal domains in thymic egress and T-cell survival, as well as in defense against viral pathogens.

Mutational Landscape of Aggressive Prostate Tumors in African American Men

Prostate cancer is the most frequently diagnosed and second most fatal nonskin cancer among men in the United States. African American men are two times more likely to develop and die of prostate cancer compared with men of other ancestries. Previous whole genome or exome tumor-sequencing studies of prostate cancer have primarily focused on men of European ancestry. In this study, we sequenced and characterized somatic mutations in aggressive (Gleason ≥7, stage ≥T2b) prostate tumors from 24 African American patients. We describe the locations and prevalence of small somatic mutations (up to 50 bases in length), copy number aberrations, and structural rearrangements in the tumor genomes compared with patient-matched normal genomes. We observed several mutation patterns consistent with previous studies, such as large copy number aberrations in chromosome 8 and complex rearrangement chains. However, TMPRSS2-ERG gene fusions and PTEN losses occurred in only 21% and 8% of the African American patients, respectively, far less common than in patients of European ancestry. We also identified mutations that appeared specific to or more common in African American patients, including a novel CDC27-OAT gene fusion occurring in 17% of patients. The genomic aberrations reported in this study warrant further investigation of their biologic significant role in the incidence and clinical outcomes of prostate cancer in African Americans. Cancer Res; 76(7); 1860-8. ©2016 AACR.

Whole genome sequence analysis of BT-474 using complete Genomics' standard and long fragment read technologies

BACKGROUND

The cell line BT-474 is a popular cell line for studying the biology of cancer and developing novel drugs. However, there is no complete, published genome sequence for this highly utilized scientific resource. In this study we sought to provide a comprehensive and useful data set for the scientific community by generating a whole genome sequence for BT-474.

FINDINGS

Five μg of genomic DNA, isolated from an early passage of the BT-474 cell line, was used to generate a whole genome sequence (114X coverage) using Complete Genomics' standard sequencing process. To provide additional variant phasing and structural variation data we also processed and analyzed two separate libraries of 5 and 6 individual cells to depths of 99X and 87X, respectively, using Complete Genomics' Long Fragment Read (LFR) technology.

CONCLUSIONS

BT-474 is a highly aneuploid cell line with an extremely complex genome sequence. This ~300X total coverage genome sequence provides a more complete understanding of this highly utilized cell line at the genomic level.

Whole Genome Sequencing Increases Molecular Diagnostic Yield Compared with Current Diagnostic Testing for Inherited Retinal Disease

PURPOSE

To compare the efficacy of whole genome sequencing (WGS) with targeted next-generation sequencing (NGS) in the diagnosis of inherited retinal disease (IRD).

DESIGN

Case series.

PARTICIPANTS

A total of 562 patients diagnosed with IRD.

METHODS

We performed a direct comparative analysis of current molecular diagnostics with WGS. We retrospectively reviewed the findings from a diagnostic NGS DNA test for 562 patients with IRD. A subset of 46 of 562 patients (encompassing potential clinical outcomes of diagnostic analysis) also underwent WGS, and we compared mutation detection rates and molecular diagnostic yields. In addition, we compared the sensitivity and specificity of the 2 techniques to identify known single nucleotide variants (SNVs) using 6 control samples with publically available genotype data.

MAIN OUTCOME MEASURES

Diagnostic yield of genomic testing.

RESULTS

Across known disease-causing genes, targeted NGS and WGS achieved similar levels of sensitivity and specificity for SNV detection. However, WGS also identified 14 clinically relevant genetic variants through WGS that had not been identified by NGS diagnostic testing for the 46 individuals with IRD. These variants included large deletions and variants in noncoding regions of the genome. Identification of these variants confirmed a molecular diagnosis of IRD for 11 of the 33 individuals referred for WGS who had not obtained a molecular diagnosis through targeted NGS testing. Weighted estimates, accounting for population structure, suggest that WGS methods could result in an overall 29% (95% confidence interval, 15-45) uplift in diagnostic yield.

CONCLUSIONS

We show that WGS methods can detect disease-causing genetic variants missed by current NGS diagnostic methodologies for IRD and thereby demonstrate the clinical utility and additional value of WGS.

Whole Genome Sequencing Expands Diagnostic Utility and Improves Clinical Management in Pediatric Medicine

The standard of care for first-tier clinical investigation of the etiology of congenital malformations and neurodevelopmental disorders is chromosome microarray analysis (CMA) for copy number variations (CNVs), often followed by gene(s)-specific sequencing searching for smaller insertion-deletions (indels) and single nucleotide variant (SNV) mutations. Whole genome sequencing (WGS) has the potential to capture all classes of genetic variation in one experiment; however, the diagnostic yield for mutation detection of WGS compared to CMA, and other tests, needs to be established. In a prospective study we utilized WGS and comprehensive medical annotation to assess 100 patients referred to a paediatric genetics service and compared the diagnostic yield versus standard genetic testing. WGS identified genetic variants meeting clinical diagnostic criteria in 34% of cases, representing a 4-fold increase in diagnostic rate over CMA (8%) (p-value = 1.42e-05) alone and >2-fold increase in CMA plus targeted gene sequencing (13%) (p-value = 0.0009). WGS identified all rare clinically significant CNVs that were detected by CMA. In 26 patients, WGS revealed indel and missense mutations presenting in a dominant (63%) or a recessive (37%) manner. We found four subjects with mutations in at least two genes associated with distinct genetic disorders, including two cases harboring a pathogenic CNV and SNV. When considering medically actionable secondary findings in addition to primary WGS findings, 38% of patients would benefit from genetic counseling. Clinical implementation of WGS as a primary test will provide a higher diagnostic yield than conventional genetic testing and potentially reduce the time required to reach a genetic diagnosis.

Intrahaplotypic Variants Differentiate Complex Linkage Disequilibrium within Human MHC Haplotypes

Distinct regions of long-range genetic fixation in the human MHC region, known as conserved extended haplotypes (CEHs), possess unique genomic characteristics and are strongly associated with numerous diseases. While CEHs appear to be homogeneous by SNP analysis, the nature of fine variations within their genomic structure is unknown. Using multiple, MHC-homozygous cell lines, we demonstrate extensive sequence conservation in two common Asian MHC haplotypes: A33-B58-DR3 and A2-B46-DR9. However, characterization of phase-resolved MHC haplotypes revealed unique intra-CEH patterns of variation and uncovered 127 single nucleotide variants (SNVs) which are missing from public databases. We further show that the strong linkage disequilibrium structure within the human MHC that typically confounds precise identification of genetic features can be resolved using intra-CEH variants, as evidenced by rs3129063 and rs448489, which affect expression of ZFP57, a gene important in methylation and epigenetic regulation. This study demonstrates an improved strategy that can be used towards genetic dissection of diseases.

Compound heterozygous mutations in the noncoding RNU4ATAC cause Roifman Syndrome by disrupting minor intron splicing

Roifman Syndrome is a rare congenital disorder characterized by growth retardation, cognitive delay, spondyloepiphyseal dysplasia and antibody deficiency. Here we utilize whole-genome sequencing of Roifman Syndrome patients to reveal compound heterozygous rare variants that disrupt highly conserved positions of the RNU4ATAC small nuclear RNA gene, a minor spliceosome component that is essential for minor intron splicing. Targeted sequencing confirms allele segregation in six cases from four unrelated families. RNU4ATAC rare variants have been recently reported to cause microcephalic osteodysplastic primordial dwarfism, type I (MOPD1), whose phenotype is distinct from Roifman Syndrome. Strikingly, all six of the Roifman Syndrome cases have one variant that overlaps MOPD1-implicated structural elements, while the other variant overlaps a highly conserved structural element not previously implicated in disease. RNA-seq analysis confirms extensive and specific defects of minor intron splicing. Available allele frequency data suggest that recessive genetic disorders caused by RNU4ATAC rare variants may be more prevalent than previously reported.

Roifman Syndrome is a rare disorder whose disease manifestations include growth retardation, spondyloepiphyseal dysplasia and immunodeficiency. Here, the authors use whole-genome sequencing to discover that rare compound heterozygous variants disrupting the small nuclear RNA gene RNU4ATAC cause Roifman Syndrome.

Spaced Seed Data Structures for De Novo Assembly

De novo assembly of the genome of a species is essential in the absence of a reference genome sequence. Many scalable assembly algorithms use the de Bruijn graph (DBG) paradigm to reconstruct genomes, where a table of subsequences of a certain length is derived from the reads, and their overlaps are analyzed to assemble sequences. Despite longer subsequences unlocking longer genomic features for assembly, associated increase in compute resources limits the practicability of DBG over other assembly archetypes already designed for longer reads. Here, we revisit the DBG paradigm to adapt it to the changing sequencing technology landscape and introduce three data structure designs for spaced seeds in the form of paired subsequences. These data structures address memory and run time constraints imposed by longer reads. We observe that when a fixed distance separates seed pairs, it provides increased sequence specificity with increased gap length. Further, we note that Bloom filters would be suitable to implicitly store spaced seeds and be tolerant to sequencing errors. Building on this concept, we describe a data structure for tracking the frequencies of observed spaced seeds. These data structure designs will have applications in genome, transcriptome and metagenome assemblies, and read error correction.

Whole-Genome Sequencing Suggests Schizophrenia Risk Mechanisms in Humans with 22q11.2 Deletion Syndrome

Chromosome 22q11.2 microdeletions impart a high but incomplete risk for schizophrenia. Possible mechanisms include genome-wide effects of DGCR8 haploinsufficiency. In a proof-of-principle study to assess the power of this model, we used high-quality, whole-genome sequencing of nine individuals with 22q11.2 deletions and extreme phenotypes (schizophrenia, or no psychotic disorder at age >50 years). The schizophrenia group had a greater burden of rare, damaging variants impacting protein-coding neurofunctional genes, including genes involved in neuron projection (nominal P = 0.02, joint burden of three variant types). Variants in the intact 22q11.2 region were not major contributors. Restricting to genes affected by a DGCR8 mechanism tended to amplify between-group differences. Damaging variants in highly conserved long intergenic noncoding RNA genes also were enriched in the schizophrenia group (nominal P = 0.04). The findings support the 22q11.2 deletion model as a threshold-lowering first hit for schizophrenia risk. If applied to a larger and thus better-powered cohort, this appears to be a promising approach to identify genome-wide rare variants in coding and noncoding sequence that perturb gene networks relevant to idiopathic schizophrenia. Similarly designed studies exploiting genetic models may prove useful to help delineate the genetic architecture of other complex phenotypes.

Gain-of-Function Mutations in ZIC1 Are Associated with Coronal Craniosynostosis and Learning Disability

Human ZIC1 (zinc finger protein of cerebellum 1), one of five homologs of the Drosophila pair-rule gene odd-paired, encodes a transcription factor previously implicated in vertebrate brain development. Heterozygous deletions of ZIC1 and its nearby paralog ZIC4 on chromosome 3q25.1 are associated with Dandy-Walker malformation of the cerebellum, and loss of the orthologous Zic1 gene in the mouse causes cerebellar hypoplasia and vertebral defects. We describe individuals from five families with heterozygous mutations located in the final (third) exon of ZIC1 (encoding four nonsense and one missense change) who have a distinct phenotype in which severe craniosynostosis, specifically involving the coronal sutures, and variable learning disability are the most characteristic features. The location of the nonsense mutations predicts escape of mutant ZIC1 transcripts from nonsense-mediated decay, which was confirmed in a cell line from an affected individual. Both nonsense and missense mutations are associated with altered and/or enhanced expression of a target gene, engrailed-2, in a Xenopus embryo assay. Analysis of mouse embryos revealed a localized domain of Zic1 expression at embryonic days 11.5-12.5 in a region overlapping the supraorbital regulatory center, which patterns the coronal suture. We conclude that the human mutations uncover a previously unsuspected role for Zic1 in early cranial suture development, potentially by regulating engrailed 1, which was previously shown to be critical for positioning of the murine coronal suture. The diagnosis of a ZIC1 mutation has significant implications for prognosis and we recommend genetic testing when common causes of coronal synostosis have been excluded.

DNAseq Workflow in a Diagnostic Context and an Example of a User Friendly Implementation

Over recent years next generation sequencing (NGS) technologies evolved from costly tools used by very few, to a much more accessible and economically viable technology. Through this recently gained popularity, its use-cases expanded from research environments into clinical settings. But the technical know-how and infrastructure required to analyze the data remain an obstacle for a wider adoption of this technology, especially in smaller laboratories. We present GensearchNGS, a commercial DNAseq software suite distributed by Phenosystems SA. The focus of GensearchNGS is the optimal usage of already existing infrastructure, while keeping its use simple. This is achieved through the integration of existing tools in a comprehensive software environment, as well as custom algorithms developed with the restrictions of limited infrastructures in mind. This includes the possibility to connect multiple computers to speed up computing intensive parts of the analysis such as sequence alignments. We present a typical DNAseq workflow for NGS data analysis and the approach GensearchNGS takes to implement it. The presented workflow goes from raw data quality control to the final variant report. This includes features such as gene panels and the integration of online databases, like Ensembl for annotations or Cafe Variome for variant sharing.

The Theory and Practice of Genome Sequence Assembly

The current genomic revolution was made possible by joint advances in genome sequencing technologies and computational approaches for analyzing sequence data. The close interaction between biologists and computational scientists is perhaps most apparent in the development of approaches for sequencing entire genomes, a feat that would not be possible without sophisticated computational tools called genome assemblers (short for genome sequence assemblers). Here, we survey the key developments in algorithms for assembling genome sequences since the development of the first DNA sequencing methods more than 35 years ago.

Copy number variants encompassing Mendelian disease genes in a large multigenerational family segregating bipolar disorder

BACKGROUND

Bipolar affective disorder (BP) is a common, highly heritable psychiatric disorder characterized by periods of depression and mania. Using dense SNP genotype data, we characterized CNVs in 388 members of an Old Order Amish Pedigree with bipolar disorder. We identified CNV regions arising from common ancestral mutations by utilizing the pedigree information. By combining this analysis with whole genome sequence data in the same individuals, we also explored the role of compound heterozygosity.

RESULTS

Here we describe 541 inherited CNV regions, of which 268 are rare in a control population of European origin but present in a large number of Amish individuals. In addition, we highlight a set of CNVs found at higher frequencies in BP individuals, and within genes known to play a role in human development and disease. As in prior reports, we find no evidence for an increased burden of CNVs in BP individuals, but we report a trend towards a higher burden of CNVs in known Mendelian disease loci in bipolar individuals (BPI and BPII, p = 0.06).

CONCLUSIONS

We conclude that CNVs may be contributing factors in the phenotypic presentation of mood disorders and co-morbid medical conditions in this family. These results reinforce the hypothesis of a complex genetic architecture underlying BP disorder, and suggest that the role of CNVs should continue to be investigated in BP data sets.