Search

Telomere length (TL) regulation is an important factor in ageing, reproduction and cancer development. Genetic, hereditary and environmental factors regulating TL are currently widely investigated, however, their relative contribution to TL variability is still understudied. We have used whole genome sequencing data of 250 family trios from the Genome of the Netherlands project to perform computational measurement of TL and a series of regression and genome-wide association analyses to reveal TL inheritance patterns and associated genetic factors. Our results confirm that TL is a largely heritable trait, primarily with mother’s, and, to a lesser extent, with father’s TL having the strongest influence on the offspring. In this cohort, mother’s, but not father’s age at conception was positively linked to offspring TL. Age-related TL attrition of 40 bp/year had relatively small influence on TL variability. Finally, we have identified TL-associated variations in ribonuclease reductase catalytic subunit M1 (RRM1 gene), which is known to regulate telomere maintenance in yeast. We also highlight the importance of multivariate approach and the limitations of existing tools for the analysis of TL as a polygenic heritable quantitative trait.

Many disciplines, from human genetics and oncology to plant breeding, microbiology and virology, commonly face the challenge of analyzing rapidly increasing numbers of genomes. In case of Homo sapiens, the number of sequenced genomes will approach hundreds of thousands in the next few years. Simply scaling up established bioinformatics pipelines will not be sufficient for leveraging the full potential of such rich genomic data sets. Instead, novel, qualitatively different Computational methods and paradigms are needed.We will witness the rapid extension of Computational pan-genomics, a new sub-area of research in Computational biology. In this article, we generalize existing definitions and understand a pangenome as any collection of genomic sequences to be analyzed jointly or to be used as a reference. We examine already available approaches to construct and use pan-genomes, discuss the potential benefits of future technologies and methodologies and review open challenges from the vantage point of the above-mentioned biological disciplines. As a prominent example for a Computational paradigm shift, we particularly highlight the transition from the representation of reference genomes as strings to representations a

Background: Current normalization methods for RNA-sequencing data allow either for intersample comparison to identify differentially expressed (DE) genes or for intrasample comparison for the discovery and validation of gene signatures. Most studies on optimization of normalization methods typically use simulated data to validate methodologies. We describe a new method, GeTMM, which allows for both inter- and intrasample analyses with the same normalized data set. We used actual (i.e. not simulated) RNA-seq data from 263 colon cancers (no biological replicates) and used the same read count data to compare GeTMM with the most commonly used normalization methods (i.e. TMM (used by edgeR), RLE (used by DESeq2) and TPM) with respect to distributions, effect of RNA quality, subtype-classification, recurrence score, recall of DE genes and correlation to RT-qPCR data. Results: We observed a clear benefit for GeTMM and TPM with regard to intrasample comparison while GeTMM performed similar to TMM and RLE normalized data in intersample comparisons. Regarding DE genes, recall was found comparable among the normalization methods, while GeTMM showed the lowest number of false-positive DE genes. Remarkably, we observed limited detrimental effects in samples with low RNA quality. Conclusions: We show that GeTMM outperforms established methods with regard to intrasample comparison while performing equivalent with regard to intersample normalization using the same normalized data. These combined properties enhance the general usefulness of RNA-seq but also the comparability to the many array-based gene expression data in the public domain.

The identification of patients with high-risk stage II colon cancer who may benefit from adjuvant therapy may allow the clinical approach to be tailored for these patients based on an understanding of tumour biology. MicroRNAs have been proposed as markers of the prognosis or treatment response in colorectal cancer. Recently, a 2-microRNA signature (l et-7i and miR-10b) was proposed to identify colorectal cancer patients at risk of developing distant metastasis. We assessed the prognostic value of this signature and additional candidate microRNAs in an independent, clinically well-defined, prospectively collected cohort of primary colon cancer patients including stage I-II colon cancer without and stage III colon cancer with adjuvant treatment. The 2-microRNA signature specifically predicted hepatic recurrence in the stage I-II group, but not the overall ability to develop distant metastasis. The addition of miR-30b to the 2-microRNA signature allowed the prediction of both distant metastasis and hepatic recurrence in patients with stage I-II colon cancer who did not receive adjuvant chemotherapy. Available gene expression data allowed us to associate m iR-30b expression with axon guidance and l et-7i expression with cell adhesion, migration, and motility.

Contains fulltext : 190288.pdf (publisher's version ) (Open Access)

Cancer affects millions of people worldwide. With the advent of novel DNA sequencing technologies,whole genome sequencing (WGS) is becoming an integral part of cancer diagnostics that can potentially enable tailored treatments of individual patients. To alleviate these problems,a user-friendly, portable and extendable SV calling workflow, sv-callers developed, that includes four state-of-the-art tools to detect SVs in cancer genomes using on-premises HPC systems. The workflow's parallel execution environment enables to scale from a single computer to high-performance compute clusters with minimal effort. The workflow supports SV analysis in either germline or somatic mode, and requires a list of (paired) WGS samples including a reference genome as input. Users may change the workflow parameters and/or software versions using the YAML configuration files. We performed SV analyses on single and paired (tumor/normal) WGS samples, and report on the results obtained using different academic HPC systems.

Contains fulltext : 177761.pdf (publisher's version ) (Open Access), OBJECTIVE: Overall and splice specific expression of Spleen Tyrosine Kinase (SYK) has been posed as a marker predicting both poor and favorable outcome in various epithelial malignancies. However, its role in colorectal cancer is largely unknown. The aim of this study was to explore the prognostic role of SYK in three cohorts of colon cancer patients. METHODS: Total messenger RNA (mRNA) expression of SYK, SYK(T), and mRNA expression of its two splice variants SYK short (S) and SYK long (L) were measured using quantitative reverse transcriptase (RT-qPCR) in 240 primary colon cancer patients (n = 160 patients with chemonaive lymph node negative [LNN] and n = 80 patients with adjuvant treated lymph node positive [LNP] colon cancer) and related to microsatellite instability (MSI), known colorectal cancer mutations, and disease-free (DFS), hepatic metastasis-free (HFS) and overall survival (OS). Two independent cohorts of patients with respectively 48 and 118 chemonaive LNN colon cancer were used for validation. RESULTS: Expression of SYK and its splice variants was significantly lower in tumors with MSI, and in KRAS wild type, BRAF mutant and PTEN mutant tumors. In a multivariate Cox regression analysis, as a continuous variable, increasing SYK(S) mRNA expression was associated with worse HFS (Hazard Ratio[HR] = 1.83; 95% Confidence Interval[CI] = 1.08-3.12; p = 0.026) in the LNN group, indicating a prognostic role for SYK(S) mRNA in patients with chemonaive LNN colon cancer. However, only a non-significant trend between SYK(S) and HFS in one of the two validation cohorts was observed (HR = 4.68; 95%CI = 0.75-29.15; p = 0.098). CONCLUSION: In our cohort, we discovered SYK(S) as a significant prognostic marker for HFS for patients with untreated LNN colon cancer. This association could however not be confirmed in two independent smaller cohorts, suggesting that further extensive validation is needed to confirm the prognostic value of SYK(S) expression in chemonaive LNN c

Genomic rearrangements that give rise to oncogenic gene fusions can offer actionable targets for cancer therapy. Here we present a systematic analysis of oncogenic gene fusions among a clinically well-characterized, prospectively collected set of 278 primary colon cancers spanning diverse tumor stages and clinical outcomes. Gene fusions and somatic genetic variations were identified in fresh frozen clinical specimens by Illumina RNA-sequencing, the STAR fusion gene detection pipeline, and GATK RNA-seq variant calling. We considered gene fusions to be pathogenically relevant when recurrent, producing divergent gene expression (outlier analysis), or as functionally important (e.g., kinase fusions). Overall, 2.5% of all specimens were defined as harboring a relevant gene fusion (kinase fusions 1.8%). Novel configurations of BRAF, NTRK3, and RET gene fusions resulting from chromosomal translocations were identified. An R-spondin fusion was found in only one tumor (0.35%), much less than an earlier reported frequency of 10% in colorectal cancers. We also found a novel fusion involving USP9X-ERAS formed by chromothripsis and leading to high expression of ERAS, a constitutively active RAS protein normally expressed only in embryonic stem cells. This USP9X–ERAS fusion appeared highly oncogenic on the basis of its ability to activate AKT signaling. Oncogenic fusions were identified only in lymph node–negative tumors that lacked BRAF or KRAS mutations. In summary, we identified several novel oncogenic gene fusions in colorectal cancer that may drive malignant development and offer new targets for personalized therapy.

Despite the clinical significance of balanced chromosomal abnormalities (BCAs), their characterization has largely been restricted to cytogenetic resolution. We explored the landscape of BCAs at nucleotide resolution in 273 subjects with a spectrum of congenital anomalies. Whole-genome sequencing revised 93% of karyotypes and demonstrated complexity that was cryptic to karyotyping in 21% of BCAs, highlighting the limitations of conventional cytogenetic approaches. At least 33.9% of BCAs resulted in gene disruption that likely contributed to the developmental phenotype, 5.2% were associated with pathogenic genomic imbalances, and 7.3% disrupted topologically associated domains (TADs) encompassing known syndromic loci. Remarkably, BCA br

OBJECTIVE: Overall and splice specific expression of Spleen Tyrosine Kinase (SYK) has been posed as a marker predicting both poor and favorable outcome in various epithelial malignancies. However, its role in colorectal cancer is largely unknown. The aim of this study was to explore the prognostic role of SYK in three cohorts of colon cancer patients.METHODS: Total messenger RNA (mRNA) expression of SYK, SYK(T), and mRNA expression of its two splice variants SYK short (S) and SYK long (L) were measured using quantitative reverse transcriptase (RT-qPCR) in 240 primary colon cancer patients (n = 160 patients with chemonaive lymph node negative [LNN] and n = 80 patients with adjuvant treated lymph node positive [LNP] colon cancer) and related to microsatellite instability (MSI), known colorectal cancer mutations, and disease-free (DFS), hepatic metastasis-free (HFS) and overall survival (OS). Two independent cohorts of patients with respectively 48 and 118 chemonaive LNN colon cancer were used for validation.RESULTS: Expression of SYK and its splice variants was significantly lower in tumors with MSI, and in KRAS wild type, BRAF mutant and PTEN mutant tumors. In a multivariate Cox regression analysis, as a continuous variable, increasing SYK(S) mRNA expression was associated with worse HFS (Hazard Ratio[HR] = 1.83; 95% Confidence Interval[

Germline mutation detection from human DNA sequence data is challenging due to the rarity of such events relative to the intrinsic error rates of sequencing technologies and the uneven coverage across the genome. We developed PhaseByTransmission (PBT) to identify de novo single nucleotide variants and short insertions and deletions (indels) from sequence data collected in parent-offspring trios. We compute the joint probability of the data given the genotype likelihoods in the individual family members, the known familial relationships and a prior probability for the mutation rate. Candidate de novo mutations (DNMs) are reported along with their posterior probability, providing a systematic way to prioritize them for validation. Our tool is integrated in the Genome Analysis Toolkit and can be used together with the ReadBackedPhasing module to infer the parental origin of DNMs based on phase-informative reads. Using simulated data, we show that PBT outperforms existing tools, especially in low coverage data and on the X chromosome. We further show that PBT displays high validation rates on empirical parent-offspring sequencing data for whole-exome data from 104 trios and X-chromosome data from 249 parent-offspring families. Finally, we demonstrate an association between father's age at conception and the number of DNMs in female offspring's X chromosome, consistent with previous literature reports.

Structural variation (SV) represents a major source of differences between individual human genomes and has been linked to disease phenotypes. However, the majority of studies provide neither a global view of the full spectrum of these variants nor integrate them into reference panels of genetic variation. Here, we analyse whole genome sequencing data of 769 individuals from 250 Dutch families, and provide a haplotype-resolved map of 1.9 million genome variants across 9 different variant classes, including novel forms of complex indels, and retrotransposition-mediated insertions of mobile elements and processed RNAs. A large proportion are previously under reported variants sized between 21 and 100 bp. We detect 4 megabases of novel sequence, encoding 11 new transcripts. Finally, we show 191 known, trait-associated SNPs to be in strong linkage disequilibrium with SVs and demonstrate that our panel facilitates accurate imputation of SVs in unrelated individuals.

Mutations create variation in the population, fuel evolution and cause genetic diseases. Current knowledge about de novo mutations is incomplete and mostly indirect. Here we analyze 11,020 de novo mutations from the whole genomes of 250 families. We show that de novo mutations in the offspring of older fathers are not only more numerous but also occur more frequently in early-replicating, genic regions. Functional regions exhibit higher mutation rates due to CpG dinucleotides and show signatures of transcription-coupled repair, whereas mutation clusters with a unique signature point to a new mutational mechanism. Mutation and recombination rates independently associate with nucleotide diversity, and regional variation in human-chimpanzee divergence is only partly explained by heterogeneity in mutation rate. Finally, we provide a genome-wide mutation rate map for medical and population genetics applications. Our results provide new insights and refine long-standing hypotheses about human mutagenesis.

1 juni 2015, Contains fulltext : 154750.pdf (publisher's version ) (Open Access), Small insertions and deletions (indels) and large structural variations (SVs) are major contributors to human genetic diversity and disease. However, mutation rates and characteristics of de novo indels and SVs in the general population have remained largely unexplored. We report 332 validated de novo structural changes identified in whole genomes of 250 families, including complex indels, retrotransposon insertions, and interchromosomal events. These data indicate a mutation rate of 2.94 indels (1-20 bp) and 0.16 SVs (>20 bp) per generation. De novo structural changes affect on average 4.1 kbp of genomic sequence and 29 coding bases per generation, which is 91 and 52 times more nucleotides than de novo substitutions, respectively. This contrasts with the equal genomic footprint of inherited SVs and substitutions. An excess of structural changes originated on paternal haplotypes. Additionally, we observed a nonuniform distribution of de novo SVs across offspring. These results reveal the importance of different mutational mechanisms to changes in human genome structure across generations.

Small insertions and deletions (indels) and large structural variations (SVs) are major contributors to human genetic diversity and disease. However, mutation rates and characteristics of de novo indels and SVs in the general population have remained largely unexplored. We report 332 validated de novo structural changes identified in whole genomes of 250 families, including complex indels, retrotransposon insertions, and interchromosomal events. These data indicate a mutation rate of 2.94 indels (120 bp) and 0.16 SVs (>20 bp) per generation. De novo structural changes affect on average 4.1 kbp of genomic sequence and 29 coding bases per generation, which is 91 and 52 times more nucleotides than de novo substitutions, respectively. This contrasts with the equal genomic footprint of inherited SVs and substitutions. An excess of structural changes originated on paternal haplotypes. Additionally, we observed a nonuniform distribution of de novo SVs across offspring. These results reveal the importance of different mutational mechanisms to changes in human genome structure across generations.

Item does not contain fulltext, OBJECTIVES: Autoinflammatory disorders are disorders of the innate immune system. Standard genetic testing provided no correct diagnosis in a female patient from a non-consanguineous family of British descent with a colchicine-responsive autosomal dominant periodic fever syndrome. We aimed to unravel the genetic cause of the symptoms. METHODS: Whole exome sequencing was used to screen for novel sequence variants, which were validated by direct Sanger sequencing. Ex vivo stimulation with peripheral blood mononuclear cells was performed to study the functional consequences of the mutation. mRNA and cytokine levels were measured by quantitative PCR and ELISA, respectively. RESULTS: Whole exome sequencing revealed a novel missense sequence variant, not seen in around 6800 controls, mapping to exon 8 of the MEFV gene (c.1730C>A; p.T577N), co-segregating perfectly with disease in this family. Other mutations at the same amino acid (c.1730C>G; p.T577S and c.1729A>T; p.T577S) were found in a family of Turkish descent, with autosomal dominant inheritance of familial Mediterranean fever (FMF)-like phenotype, and a Dutch patient, respectively. Moreover, a mutation (c.1729A>G; p.T577A) was detected in two Dutch siblings, who had episodes of inflammation of varying severity not resembling FMF. Peripheral blood mononuclear cells from one patient of the index family showed increased basal interleukin 1beta mRNA levels and cytokine responses after lipopolysaccharide stimulation. Responses normalised with colchicine treatment. CONCLUSIONS: Heterozygous mutations at amino acid position 577 of pyrin can induce an autosomal dominant autoinflammatory syndrome. This suggests that T577, located in front of the C-terminal B30.2/SPRY domain, is crucial for pyrin function.

Contains fulltext : 137213.pdf (publisher's version ) (Closed access), Whole-genome sequencing enables complete characterization of genetic variation, but geographic clustering of rare alleles demands many diverse populations be studied. Here we describe the Genome of the Netherlands (GoNL) Project, in which we sequenced the whole genomes of 250 Dutch parent-offspring families and constructed a haplotype map of 20.4 million single-nucleotide variants and 1.2 million insertions and deletions. The intermediate coverage ( approximately 13x) and trio design enabled extensive characterization of structural variation, including midsize events (30-500 bp) previously poorly catalogued and de novo mutations. We demonstrate that the quality of the haplotypes boosts imputation accuracy in independent samples, especially for lower frequency alleles. Population genetic analyses demonstrate fine-scale structure across the country and support multiple ancient migrations, consistent with historical changes in sea level and flooding. The GoNL Project illustrates how single-population whole-genome sequencing can provide detailed characterization of genetic variation and may guide the design of future population studies.

Contains fulltext : 139109.pdf (publisher's version ) (Open Access), Genomic rearrangements are a common cause of human congenital abnormalities. However, their origin and consequences are poorly understood. We performed molecular analysis of two patients with congenital disease who carried de novo genomic rearrangements. We found that the rearrangements in both patients hit genes that are recurrently rearranged in cancer (ETV1, FOXP1, and microRNA cluster C19MC) and drive formation of fusion genes similar to those described in cancer. Subsequent analysis of a large set of 552 de novo germline genomic rearrangements underlying congenital disorders revealed enrichment for genes rearranged in cancer and overlap with somatic cancer breakpoints. Breakpoints of common (inherited) germline structural variations also overlap with cancer breakpoints but are depleted for cancer genes. We propose that the same genomic positions are prone to genomic rearrangements in germline and soma but that timing and context of breakage determines whether developmental defects or cancer are promoted.

BACKGROUND: Paired-tag sequencing approaches are commonly used for the analysis of genome structure. However, mammalian genomes have a complex organization with a variety of repetitive elements that complicate comprehensive genome-wide analyses. RESULTS: Here, we systematically assessed the utility of paired-end and mate-pair (MP) next-generation sequencing libraries with insert sizes ranging from 170 bp to 25 kb, for genome coverage and for improving scaffolding of a mammalian genome (Rattus norvegicus). Despite a lower library complexity, large insert MP libraries (20 or 25 kb) provided very high physical genome coverage and were found to efficiently span repeat elements in the genome. Medium-sized (5, 8 or 15 kb) MP libraries were much more efficient for genome structure analysis than the more commonly used shorter insert paired-end and 3 kb MP libraries. Furthermore, the combination of medium- and large insert libraries resulted in a 3-fold increase in N50 in scaffolding processes. Finally, we show that our data can be used to evaluate and improve contig order and orientation in the current rat reference genome assembly. CONCLUSIONS: We conclude that applying combinations of mate-pair libraries with insert sizes that match the distributions of repetitive elements improves contig scaffolding and can contribute to the finishing of draft genomes., BACKGROUND: Paired-tag sequencing approaches are commonly used for the analysis of genome structure. However, mammalian genomes have a complex organization with a variety of repetitive elements that complicate comprehensive genome-wide analyses. RESULTS: Here, we systematically assessed the utility of paired-end and mate-pair (MP) next-generation sequencing libraries with insert sizes ranging from 170 bp to 25 kb, for genome coverage and for improving scaffolding of a mammalian genome (Rattus norvegicus). Despite a lower library complexity, large insert MP libraries (20 or 25 kb) provided very high physical genome coverage and were found to efficiently span repeat elements in the genome. Medium-sized (5, 8 or 15 kb) MP libraries were much more efficient for genome structure analysis than the more commonly used shorter insert paired-end and 3 kb MP libraries. Furthermore, the combination of medium- and large insert libraries resulted in a 3-fold increase in N50 in scaffolding processes. Finally, we show that our data can be used to evaluate and improve contig order and orientation in the current rat reference genome assembly. CONCLUSIONS: We conclude that applying combinations of mate-pair libraries with insert sizes that match the distributions of repetitive elements improves contig scaffolding and can contribute to the finishing of draft genomes.

Genomic rearrangements may give rise to congenital disease and contribute to cancer development. Recent evidence has shown that very complex genomic rearrangements in cancer cells can result from a single catastrophic event of massive DNA breakage and repair, termed chromothripsis. This results in heavily rearranged chromosomes comprising frequent sequence losses. A very similar process of chromosome shattering is found for complex chromosome rearrangements in the germline of patients with congenital disorders. Here, we review the literature on chromothripsis in cancer and congenital disease. We describe differences and similarities for chromothripsis rearrangements in somatic tissue and the germ line and we discuss the cellular origin and molecular mechanisms of chromothripsis., Genomic rearrangements may give rise to congenital disease and contribute to cancer development. Recent evidence has shown that very complex genomic rearrangements in cancer cells can result from a single catastrophic event of massive DNA breakage and repair, termed chromothripsis. This results in heavily rearranged chromosomes comprising frequent sequence losses. A very similar process of chromosome shattering is found for complex chromosome rearrangements in the germline of patients with congenital disorders. Here, we review the literature on chromothripsis in cancer and congenital disease. We describe differences and similarities for chromothripsis rearrangements in somatic tissue and the germ line and we discuss the cellular origin and molecular mechanisms of chromothripsis.

1 januari 2012, Item does not contain fulltext, OBJECTIVE: To evaluate the significance of NR5A1 mutations in a large, well-phenotyped cohort of women with primary ovarian insufficiency (POI). Mutations in the NR5A1 gene (SF-1) were previously described in disorders of sexual development and adrenal insufficiency. Recently, a high frequency of NR5A1 gene mutations was reported in a small group of women with POI. DESIGN: Cross-sectional cohort study. SETTING: University hospital. PATIENT(S): Well-phenotyped women (n = 386) with secondary amenorrhea and diagnosed with POI, including women with familial POI (n = 77). INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): The entire coding region and splice sites of the NR5A1 gene were PCR-amplified and sequenced. The pathogenicity of identified mutations was predicted in silico by assessing Align-GVGD class and Grantham score. RESULT(S): Sequencing was successful in 356 patients with POI. In total, 9 mutations were identified in 10 patients. Five of these mutations concerned novel nonconservative mutations occurring in 5 patients. Prediction of effect on protein function showed low to intermediate pathogenicity for all nonconservative mutations. The overall NR5A1 gene mutation rate was 1.4%. CONCLUSION(S): The current study demonstrates that mutations in the NR5A1 gene are rare in women with POI. Primary ovarian insufficiency remains unexplained in the great majority of patients; therefore, continued efforts are needed to elucidate its underlying genetic factors.

Structural variations are genomic rearrangements that contribute significantly to evolution, natural variation between humans, and are often involved in genetic disorders. Cellular stresses and errors in repair mechanisms can lead to a large variety of structural variation events throughout the genome. Traditional microscopy- and array-based methods are used for the detection of larger events or copy number variations. Next generation sequencing has in theory enabled the detection of all types of structural variants in the human genome at unprecedented accuracy. In practice, a significant challenge lies in the development of computational methods that are able to identify these structural variants based on the generated data. In the last several years, many tools have been developed based on four different categories of information that can be obtained from sequencing experiments: read pairs, read depths, split reads and assembled sequences. In this thesis, I first introduce the topic of structural variation by discussing its impact in various areas, what mechanisms can lead to its formation, and the types of structural variation that can occur. Subsequently, I describe the array-based and sequencing-based methods that can be used to detect structural variation. Finally, I give an overview of the tools that are currently available to detect signatures of structural variants in NGS data and their properties, and conclude by discussing the current capabilities of these tools, possible future directions and expectations for the future

Item does not contain fulltext, Cantu syndrome is characterized by congenital hypertrichosis, distinctive facial features, osteochondrodysplasia and cardiac defects. By using family-based exome sequencing, we identified a de novo mutation in ABCC9. Subsequently, we discovered novel dominant missense mutations in ABCC9 in 14 of the 16 individuals with Cantu syndrome examined. The ABCC9 protein is part of an ATP-dependent potassium (K(ATP)) channel that couples the metabolic state of a cell with its electrical activity. All mutations altered amino acids in or close to the transmembrane domains of ABCC9. Using electrophysiological measurements, we show that mutations in ABCC9 reduce the ATP-mediated potassium channel inhibition, resulting in channel opening. Moreover, similarities between the phenotype of individuals with Cantu syndrome and side effects from the K(ATP) channel agonist minoxidil indicate that the mutations in ABCC9 result in channel opening. Given the availability of ABCC9 antagonists, our findings may have direct implications for the treatment of individuals with Cantu syndrome.

Array-based genome-wide segmental aneuploidy screening detects both de novo and inherited copy number variations (CNVs). In sporadic patients de novo CNVs are interpreted as potentially pathogenic. However, a deletion, transmitted from a healthy parent, may be pathogenic if it overlaps with a mutated second allele inherited from the other healthy parent. To detect such events, we performed multiplex enrichment and next-generation sequencing of the entire coding sequence of all genes within unique hemizygous deletion regions in 20 patients (1.53 Mb capture footprint). Out of the detected 703 non-synonymous single-nucleotide variants (SNVs), 8 represented variants being unmasked by a hemizygous deletion. Although evaluation of inheritance patterns, Grantham matrix scores, evolutionary conservation and bioinformatic predictions did not consistently indicate pathogenicity of these variants, no definitive conclusions can be drawn without functional validation. However, in one patient with severe mental retardation, lack of speech, microcephaly, cheilognathopalatoschisis and bilateral hearing loss, we discovered a second smaller deletion, inherited from the other healthy parent, resulting in loss of both alleles of the highly conserved heat shock factor binding protein 1 (HSBP1) gene. Conceivably, inherited deletions may unmask rare pathogenic variants that may exert a phenotypic impact through a recessive mode of gene action.European Journal of Human Genetics advance online publication, 18 January 2012; doi:10.1038/ejhg.2011.263., Array-based genome-wide segmental aneuploidy screening detects both de novo and inherited copy number variations (CNVs). In sporadic patients de novo CNVs are interpreted as potentially pathogenic. However, a deletion, transmitted from a healthy parent, may be pathogenic if it overlaps with a mutated second allele inherited from the other healthy parent. To detect such events, we performed multiplex enrichment and next-generation sequencing of the entire coding sequence of all genes within unique hemizygous deletion regions in 20 patients (1.53 Mb capture footprint). Out of the detected 703 non-synonymous single-nucleotide variants (SNVs), 8 represented variants being unmasked by a hemizygous deletion. Although evaluation of inheritance patterns, Grantham matrix scores, evolutionary conservation and bioinformatic predictions did not consistently indicate pathogenicity of these variants, no definitive conclusions can be drawn without functional validation. However, in one patient with severe mental retardation, lack of speech, microcephaly, cheilognathopalatoschisis and bilateral hearing loss, we discovered a second smaller deletion, inherited from the other healthy parent, resulting in loss of both alleles of the highly conserved heat shock factor binding protein 1 (HSBP1) gene. Conceivably, inherited deletions may unmask rare pathogenic variants that may exert a phenotypic impact through a recessive mode of gene action.European Journal of Human Genetics advance online publication, 18 January 2012; doi:10.1038/ejhg.2011.263.

Chromothripsis represents a novel phenomenon in the structural variation landscape of cancer genomes. Here, we analyze the genomes of ten patients with congenital disease who were preselected to carry complex chromosomal rearrangements with more than two breakpoints. The rearrangements displayed unanticipated complexity resembling chromothripsis. We find that eight of them contain hallmarks of multiple clustered double-stranded DNA breaks (DSBs) on one or more chromosomes. In addition, nucleotide resolution analysis of 98 breakpoint junctions indicates that break repair involves nonhomologous or microhomology-mediated end joining. We observed that these eight rearrangements are balanced or contain sporadic deletions ranging in size between a few hundred base pairs and several megabases. The two remaining complex rearrangements did not display signs of DSBs and contain duplications, indicative of rearrangement processes involving template switching. Our work provides detailed insight into the characteristics of chromothripsis and supports a role for clustered DSBs driving some constitutional chromothripsis rearrangements., Chromothripsis represents a novel phenomenon in the structural variation landscape of cancer genomes. Here, we analyze the genomes of ten patients with congenital disease who were preselected to carry complex chromosomal rearrangements with more than two breakpoints. The rearrangements displayed unanticipated complexity resembling chromothripsis. We find that eight of them contain hallmarks of multiple clustered double-stranded DNA breaks (DSBs) on one or more chromosomes. In addition, nucleotide resolution analysis of 98 breakpoint junctions indicates that break repair involves nonhomologous or microhomology-mediated end joining. We observed that these eight rearrangements are balanced or contain sporadic deletions ranging in size between a few hundred base pairs and several megabases. The two remaining complex rearrangements did not display signs of DSBs and contain duplications, indicative of rearrangement processes involving template switching. Our work provides detailed insight into the characteristics of chromothripsis and supports a role for clustered DSBs driving some constitutional chromothripsis rearrangements.

Research on DNA has evolved from the discovery of the double-helix structure in 1953 to structural variations today. Structural variations are all genomic rearrangements bigger than one base pair. This definition includes deletions, insertions, translocations, inversions, and duplications. Genomic rearrangements can have an influence on phenotype, and are thus associated with diseases. A Structural variation in a somatic cell might change susceptibility to cancer while a de novo rearrangement in a germ cell might result in congenital defects. Sequencing the break point can aid in relating the variant to a phenotypic effect and may help identifying a mutational mechanism. Three major mechanisms have currently been suggested. NAHR and NHEJ are double strand DNA break repair mechanisms. FoSTeS (or MMBIR) is a replication based mechanism. Chromothripsis, retrotransposition, alternative FoSTeS and alternative end-joining (MMEJ) are also suggested mechanisms, resulting in structural variations. Finding and defining both pathogenic and non pathogenic structural variations is important, since we will then be able to establish the cause for some diseases. In the project described in this article, the occurrence of four recurrent non pathogenic deletions in the population was determined. This experiment shows that non-pathogenic rearrangements are quite common in the population. The deletions in chromosomes 1, 5, 22, and the X chromosome are present in 35% to 93% of the population. Furthermore, a second experiment was performed in which structural variations of two children with congenital defects were sequenced by capillary sequencing. The goal of this experiment was to identify a possible cause for their abnormalities and to establish which mutational mechanism could have led to the structural variation. No de novo mutations were found in one of the patients. Two mutations that he inherited from his mother were caused by MMEJ and retrotransposition. In the other patient, tw

The costs and efforts for sample preparation of hundreds of individuals, their genomic enrichment for regions of interest, and sufficient deep sequencing bring a significant burden to next-generation sequencing-based experiments. We investigated whether pooling of samples at the level of genomic DNA would be a more versatile strategy for lowering the costs and efforts for common disease-associated rare variant detection in candidate genes or associated loci in a substantial patient cohort. We performed a pilot experiment using five pools of 20 abdominal aortic aneurysm (AAA) patients that were enriched on separate microarrays for the reported 9p21.3 associated locus and 42 additional AAA candidate genes, and sequenced on the SOLiD platform. Here, we discuss challenges and limitations connected to this approach and show that the high number of novel variants detected per pool and allele frequency deviations to the usually highly false positive cut-off region for variant detection in non-pooled samples can be limiting factors for successful variant prioritization and confirmation. We conclude that barcode indexing of individual samples before pooling followed by a multiplexed enrichment strategy should be preferred for detection of rare genetic variants in larger sample sets rather than a genomic DNA pooling strategy., The costs and efforts for sample preparation of hundreds of individuals, their genomic enrichment for regions of interest, and sufficient deep sequencing bring a significant burden to next-generation sequencing-based experiments. We investigated whether pooling of samples at the level of genomic DNA would be a more versatile strategy for lowering the costs and efforts for common disease-associated rare variant detection in candidate genes or associated loci in a substantial patient cohort. We performed a pilot experiment using five pools of 20 abdominal aortic aneurysm (AAA) patients that were enriched on separate microarrays for the reported 9p21.3 associated locus and 42 additional AAA candidate genes, and sequenced on the SOLiD platform. Here, we discuss challenges and limitations connected to this approach and show that the high number of novel variants detected per pool and allele frequency deviations to the usually highly false positive cut-off region for variant detection in non-pooled samples can be limiting factors for successful variant prioritization and confirmation. We conclude that barcode indexing of individual samples before pooling followed by a multiplexed enrichment strategy should be preferred for detection of rare genetic variants in larger sample sets rather than a genomic DNA pooling strategy.

A variety of mutational mechanisms shape the dynamic architecture of human genomes and occasionally result in congenital defects and disease. Here, we used genome-wide long mate-pair sequencing to systematically screen for inherited and de novo structural variation in a trio including a child with severe congenital abnormalities. We identified 4321 inherited structural variants and 17 de novo rearrangements. We characterized the de novo structural changes to the base-pair level revealing a complex series of balanced inter- and intra-chromosomal rearrangements consisting of 12 breakpoints involving chromosomes 1, 4 and 10. Detailed inspection of breakpoint regions indicated that a series of simultaneous double-stranded DNA breaks caused local shattering of chromosomes. Fusion of the resulting chromosomal fragments involved non-homologous end joining, since junction points displayed limited or no homology and small insertions and deletions. The pattern of random joining of chromosomal fragments that we observe here strongly resembles the somatic rearrangement patterns--termed chromothripsis--that have recently been described in deranged cancer cells. We conclude that a similar mechanism may also drive the formation of de novo structural variation in the germline. [KEYWORDS: Base Sequence, Child, Chromosome Aberrations, Chromosome Breakage, Chromosomes, Human, Pair 1/genetics, Chromosomes, Human, Pair 10/genetics, Chromosomes, Human, Pair 4/genetics, Computational Biology, Female, Gene Order, Gene Rearrangement/ genetics, Germ Cells, Humans, Male, Models, Genetic, Molecular Sequence Data, Sequence Analysis, DNA], A variety of mutational mechanisms shape the dynamic architecture of human genomes and occasionally result in congenital defects and disease. Here, we used genome-wide long mate-pair sequencing to systematically screen for inherited and de novo structural variation in a trio including a child with severe congenital abnormalities. We identified 4321 inherited structural variants and 17 de novo rearrangements. We characterized the de novo structural changes to the base-pair level revealing a complex series of balanced inter- and intra-chromosomal rearrangements consisting of 12 breakpoints involving chromosomes 1, 4 and 10. Detailed inspection of breakpoint regions indicated that a series of simultaneous double-stranded DNA breaks caused local shattering of chromosomes. Fusion of the resulting chromosomal fragments involved non-homologous end joining, since junction points displayed limited or no homology and small insertions and deletions. The pattern of random joining of chromosomal fragments that we observe here strongly resembles the somatic rearrangement patterns--termed chromothripsis--that have recently been described in deranged cancer cells. We conclude that a similar mechanism may also drive the formation of de novo structural variation in the germline. [KEYWORDS: Base Sequence, Child, Chromosome Aberrations, Chromosome Breakage, Chromosomes, Human, Pair 1/genetics, Chromosomes, Human, Pair 10/genetics, Chromosomes, Human, Pair 4/genetics, Computational Biology, Female, Gene Order, Gene Rearrangement/ genetics, Germ Cells, Humans, Male, Models, Genetic, Molecular Sequence Data, Sequence Analysis, DNA]

The unprecedented increase in the throughput of DNA sequencing driven by next-generation technologies now allows efficient analysis of the complete protein-coding regions of genomes (exomes) for multiple samples in a single sequencing run. However, sample preparation and targeted enrichment of multiple samples has become a rate-limiting and costly step in high-throughput genetic analysis. Here we present an efficient protocol for parallel library preparation and targeted enrichment of pooled multiplexed bar-coded samples. The procedure is compatible with microarray-based and solution-based capture approaches. The high flexibility of this method allows multiplexing of 3-5 samples for whole-exome experiments, 20 samples for targeted footprints of 5 Mb and 96 samples for targeted footprints of 0.4 Mb. From library preparation to post-enrichment amplification, including hybridization time, the protocol takes 5-6 d for array-based enrichment and 3-4 d for solution-based enrichment. Our method provides a cost-effective approach for a broad range of applications, including targeted resequencing of large sample collections (e.g., follow-up genome-wide association studies), and whole-exome or custom mini-genome sequencing projects. This protocol gives details for a single-tube procedure, but scaling to a manual or automated 96-well plate format is possible and discussed., The unprecedented increase in the throughput of DNA sequencing driven by next-generation technologies now allows efficient analysis of the complete protein-coding regions of genomes (exomes) for multiple samples in a single sequencing run. However, sample preparation and targeted enrichment of multiple samples has become a rate-limiting and costly step in high-throughput genetic analysis. Here we present an efficient protocol for parallel library preparation and targeted enrichment of pooled multiplexed bar-coded samples. The procedure is compatible with microarray-based and solution-based capture approaches. The high flexibility of this method allows multiplexing of 3-5 samples for whole-exome experiments, 20 samples for targeted footprints of 5 Mb and 96 samples for targeted footprints of 0.4 Mb. From library preparation to post-enrichment amplification, including hybridization time, the protocol takes 5-6 d for array-based enrichment and 3-4 d for solution-based enrichment. Our method provides a cost-effective approach for a broad range of applications, including targeted resequencing of large sample collections (e.g., follow-up genome-wide association studies), and whole-exome or custom mini-genome sequencing projects. This protocol gives details for a single-tube procedure, but scaling to a manual or automated 96-well plate format is possible and discussed.

BACKGROUND: Structural rearrangements form a major class of somatic variation in cancer genomes. Local chromosome shattering, termed chromothripsis, is a mechanism proposed to be the cause of clustered chromosomal rearrangements and was recently described to occur in a small percentage of tumors. The significance of these clusters for tumor development or metastatic spread is largely unclear. RESULTS: We used genome-wide long mate-pair sequencing and SNP array profiling to reveal that chromothripsis is a widespread phenomenon in primary colorectal cancer and metastases. We find large and small chromothripsis events in nearly every colorectal tumor sample and show that several breakpoints of chromothripsis clusters and isolated rearrangements affect cancer genes, including NOTCH2, EXO1 and MLL3. We complemented the structural variation studies by sequencing the coding regions of a cancer exome in all colorectal tumor samples and found somatic mutations in 24 genes, including APC, KRAS, SMAD4 and PIK3CA. A pairwise comparison of somatic variations in primary and metastatic samples indicated that many chromothripsis clusters, isolated rearrangements and point mutations are exclusively present in either the primary tumor or the metastasis and may affect cancer genes in a lesion-specific manner. CONCLUSIONS: We conclude that chromothripsis is a prevalent mechanism driving structural rearrangements in colorectal cancer and show that a complex interplay between point mutations, simple copy number changes and chromothripsis events drive colorectal tumor development and metastasis., BACKGROUND: Structural rearrangements form a major class of somatic variation in cancer genomes. Local chromosome shattering, termed chromothripsis, is a mechanism proposed to be the cause of clustered chromosomal rearrangements and was recently described to occur in a small percentage of tumors. The significance of these clusters for tumor development or metastatic spread is largely unclear. RESULTS: We used genome-wide long mate-pair sequencing and SNP array profiling to reveal that chromothripsis is a widespread phenomenon in primary colorectal cancer and metastases. We find large and small chromothripsis events in nearly every colorectal tumor sample and show that several breakpoints of chromothripsis clusters and isolated rearrangements affect cancer genes, including NOTCH2, EXO1 and MLL3. We complemented the structural variation studies by sequencing the coding regions of a cancer exome in all colorectal tumor samples and found somatic mutations in 24 genes, including APC, KRAS, SMAD4 and PIK3CA. A pairwise comparison of somatic variations in primary and metastatic samples indicated that many chromothripsis clusters, isolated rearrangements and point mutations are exclusively present in either the primary tumor or the metastasis and may affect cancer genes in a lesion-specific manner. CONCLUSIONS: We conclude that chromothripsis is a prevalent mechanism driving structural rearrangements in colorectal cancer and show that a complex interplay between point mutations, simple copy number changes and chromothripsis events drive colorectal tumor development and metastasis.

MicroRNAs (miRNAs) are noncoding RNA molecules of 21 to 24 nt that regulate the expression of target genes in a post-transcriptional manner. Although evidence indicates that miRNAs play essential roles in embryogenesis, cell differentiation and pathogenesis of human diseases, extensive miRNA profiling in cells or tissues has been hampered by the lack of sensitive cloning methods. Here we describe a highly efficient profiling method, termed miRNA amplification profiling (mRAP), as well as its application both to mouse embryos at various developmental stages and to adult mouse organs. A total of 77,436 Small-RNA species was sequenced, with 11,776 of these sequences found to match previously described miRNAs. With the use of a newly developed computational prediction algorithm, we further identified 229 independent candidates for previously unknown miRNAs. The expression of some of these candidate miRNAs was confirmed by northern blot analysis and whole-mount in situ hybridization. Our data thus indicate that the total number of miRNAs in vertebrates is larger than previously appreciated and that the expression of these molecules is tightly controlled in a tissue- and developmental stage-specific manner., MicroRNAs (miRNAs) are noncoding RNA molecules of 21 to 24 nt that regulate the expression of target genes in a post-transcriptional manner. Although evidence indicates that miRNAs play essential roles in embryogenesis, cell differentiation and pathogenesis of human diseases, extensive miRNA profiling in cells or tissues has been hampered by the lack of sensitive cloning methods. Here we describe a highly efficient profiling method, termed miRNA amplification profiling (mRAP), as well as its application both to mouse embryos at various developmental stages and to adult mouse organs. A total of 77,436 Small-RNA species was sequenced, with 11,776 of these sequences found to match previously described miRNAs. With the use of a newly developed computational prediction algorithm, we further identified 229 independent candidates for previously unknown miRNAs. The expression of some of these candidate miRNAs was confirmed by northern blot analysis and whole-mount in situ hybridization. Our data thus indicate that the total number of miRNAs in vertebrates is larger than previously appreciated and that the expression of these molecules is tightly controlled in a tissue- and developmental stage-specific manner.

MicroRNAs (miRNAs) control gene expression by translational inhibition and destabilization of mRNAs. While hundreds of miRNAs have been found, only a few have been studied in detail. miRNAs have been implicated in tissue morphogenesis, cellular processes like apoptosis, and major signaling pathways. Emerging evidence suggests a direct link between miRNAs and disease, and miRNA expression signatures are associated with various types of cancer. In addition, the gain and loss of miRNA target sites appears to be causal to some genetic disorders. Here, we discuss the current literature on the role of miRNAs in animal development and disease., MicroRNAs (miRNAs) control gene expression by translational inhibition and destabilization of mRNAs. While hundreds of miRNAs have been found, only a few have been studied in detail. miRNAs have been implicated in tissue morphogenesis, cellular processes like apoptosis, and major signaling pathways. Emerging evidence suggests a direct link between miRNAs and disease, and miRNA expression signatures are associated with various types of cancer. In addition, the gain and loss of miRNA target sites appears to be causal to some genetic disorders. Here, we discuss the current literature on the role of miRNAs in animal development and disease.

MicroRNAs (miRNAs) play an important role in development and regulate the expression of many animal genes by post-transcriptional gene silencing. Here we describe the cloning and expression of new miRNAs from zebrafish. By high-throughput sequencing of small-RNA cDNA libraries from 5-day-old zebrafish larvae and adult zebrafish brain we found 139 known miRNAs and 66 new miRNAs. For 65 known miRNAs and for 11 new miRNAs we also cloned the miRNA star sequence. We analyzed the temporal and spatial expression patterns for 35 new miRNAs and for 32 known miRNAs in the zebrafish by whole mount in situ hybridization and northern blotting. Overall, 23 of the 35 new miRNAs and 30 of the 32 known miRNAs could be detected. We found that most miRNAs were expressed during later stages of development. Some were expressed ubiquitously, but many of the miRNAs were expressed in a tissue-specific manner. Most newly discovered miRNAs have low expression levels and are less conserved in other vertebrate species. Our cloning and expression analysis indicates that most abundant and conserved miRNAs in zebrafish are now known., MicroRNAs (miRNAs) play an important role in development and regulate the expression of many animal genes by post-transcriptional gene silencing. Here we describe the cloning and expression of new miRNAs from zebrafish. By high-throughput sequencing of small-RNA cDNA libraries from 5-day-old zebrafish larvae and adult zebrafish brain we found 139 known miRNAs and 66 new miRNAs. For 65 known miRNAs and for 11 new miRNAs we also cloned the miRNA star sequence. We analyzed the temporal and spatial expression patterns for 35 new miRNAs and for 32 known miRNAs in the zebrafish by whole mount in situ hybridization and northern blotting. Overall, 23 of the 35 new miRNAs and 30 of the 32 known miRNAs could be detected. We found that most miRNAs were expressed during later stages of development. Some were expressed ubiquitously, but many of the miRNAs were expressed in a tissue-specific manner. Most newly discovered miRNAs have low expression levels and are less conserved in other vertebrate species. Our cloning and expression analysis indicates that most abundant and conserved miRNAs in zebrafish are now known.

MicroRNAs (miRNAs) attenuate gene expression by means of translational inhibition and mRNA degradation. They are abundant, highly conserved, and predicted to regulate a large number of transcripts. Several hundred miRNA classes are known, and many are associated with cell proliferation and differentiation. Many exhibit tissue-specific expression, which aids in evaluating their functions, and it has been assumed that their high level of sequence conservation implies a high level of expression conservation. A limited amount of data supports this, although discrepancies do exist. By comparing the expression of approximately 100 miRNAs in medaka and chicken with existing data for zebrafish and mouse, we conclude that the timing and location of miRNA expression is not strictly conserved. In some instances, differences in expression are associated with changes in miRNA copy number, genomic context, or both between species. Variation in miRNA expression is more pronounced the greater the differences in physiology, and it is enticing to speculate that changes in miRNA expression may play a role in shaping the physiological differences produced during animal development., MicroRNAs (miRNAs) attenuate gene expression by means of translational inhibition and mRNA degradation. They are abundant, highly conserved, and predicted to regulate a large number of transcripts. Several hundred miRNA classes are known, and many are associated with cell proliferation and differentiation. Many exhibit tissue-specific expression, which aids in evaluating their functions, and it has been assumed that their high level of sequence conservation implies a high level of expression conservation. A limited amount of data supports this, although discrepancies do exist. By comparing the expression of approximately 100 miRNAs in medaka and chicken with existing data for zebrafish and mouse, we conclude that the timing and location of miRNA expression is not strictly conserved. In some instances, differences in expression are associated with changes in miRNA copy number, genomic context, or both between species. Variation in miRNA expression is more pronounced the greater the differences in physiology, and it is enticing to speculate that changes in miRNA expression may play a role in shaping the physiological differences produced during animal development.