The primacy of NF1 loss as the driver of tumorigenesis in neurofibromatosis type 1-associated plexiform neurofibromas.

Neurofibromatosis type 1 (NF1) is a common tumor-predisposition disorder due to germline mutations in the tumor suppressor gene NF1. A virtually pathognomonic finding of NF1 is the plexiform neurofibroma (PN), a benign, likely congenital tumor that arises from bi-allelic inactivation of NF1. PN can undergo transformation to a malignant peripheral nerve sheath tumor, an aggressive soft-tissue sarcoma. To better understand the non-NF1 genetic contributions to PN pathogenesis, we performed whole-exome sequencing, RNASeq profiling and genome-wide copy-number determination for 23 low-passage Schwann cell cultures established from surgical PN material with matching germline DNA. All resected tumors were derived from routine debulking surgeries. None of the tumors were considered at risk for malignant transformation at the time; for example, there was no pain or rapid growth. Deep (~500X) NF1 exon sequencing was also conducted on tumor DNA. Non-NF1 somatic mutation verification was performed using the Ampliseq/IonTorrent platform. We identified 100% of the germline NF1 mutations and found somatic NF1 inactivation in 74% of the PN. One individual with three PNs had different NF1 somatic mutations in each tumor. The median number of somatic mutations per sample, including NF1, was one (range 0-8). NF1 was the only gene that was recurrently somatically inactivated in multiple tumors. Gene Set Enrichment Analysis of transcriptome-wide tumor RNA sequencing identified five significant (FDR<0.01) and seven trending (0.01⩽FDR<0.02) gene sets related to DNA replication, telomere maintenance and elongation, cell cycle progression, signal transduction and cell proliferation. We found no recurrent non-NF1 locus copy-number variation in PN. This is the first multi-sample whole-exome and whole-transcriptome sequencing study of NF1-associated PN. Taken together with concurrent copy-number data, our comprehensive genetic analysis reveals the primacy of NF1 loss as the driver of PN tumorigenesis.

Exome Sequencing and High-Density Microarray Testing in Monozygotic Twin Pairs Discordant for Features of VACTERL Association.

Exome sequencing offers an efficient and affordable method to interrogate genetic factors involved in human disease. Performing exome sequencing of monozygotic twins discordant for VACTERL (Vertebral anomalies, Anal atresia, Cardiac malformations, Tracheo-Esophageal fistula, Renal anomalies, and Limb abnormalities) association-type congenital malformations was hypothesized to potentially reveal discordant variants that could demonstrate disease cause(s). After demonstrating monozygosity, we applied high-density microarrays and exome sequencing to 2 twin pairs in which 1 twin had features of VACTERL association while the other was phenotypically normal (demonstrated through comprehensive clinical and radiological evaluation). No obvious discordant genotypic results were found that would explain phenotypic discordance. We conclude that VACTERL association is a complex disease, and while performing microarray analysis and exome sequencing on phenotypically discordant monozygotic twins may hypothetically reveal genetic causes of disorders, challenges remain in applying these methods in this circumstance.

Applying Genomic Analysis to Newborn Screening.

Large-scale genomic analysis such as whole-exome and whole-genome sequencing is becoming increasingly prevalent in the research arena. Clinically, many potential uses of this technology have been proposed. One such application is the extension or augmentation of newborn screening. In order to explore this application, we examined data from 3 children with normal newborn screens who underwent whole-exome sequencing as part of research participation. We analyzed sequence information for 151 selected genes associated with conditions ascertained by newborn screening. We compared findings with publicly available databases and results from over 500 individuals who underwent whole-exome sequencing at the same facility. Novel variants were confirmed through bidirectional dideoxynucleotide sequencing. High-density microarrays (Illumina Omni1-Quad) were also performed to detect potential copy number variations affecting these genes. We detected an average of 87 genetic variants per individual. After excluding artifacts, 96% of the variants were found to be reported in public databases and have no evidence of pathogenicity. No variants were identified that would predict disease in the tested individuals, which is in accordance with their normal newborn screens. However, we identified 6 previously reported variants and 2 novel variants that, according to published literature, could result in affected offspring if the reproductive partner were also a mutation carrier; other specific molecular findings highlight additional means by which genomic testing could augment newborn screening.

A high-resolution 15,000(Rad) radiation hybrid panel for the domestic cat.

The current genetic and recombination maps of the cat have fewer than 3,000 markers and a resolution limit greater than 1 Mb. To complement the first-generation domestic cat maps, support higher resolution mapping studies, and aid genome assembly in specific areas as well as in the whole genome, a 15,000(Rad) radiation hybrid (RH) panel for the domestic cat was generated. Fibroblasts from the female Abyssinian cat that was used to generate the cat genomic sequence were fused to a Chinese hamster cell line (A23), producing 150 hybrid lines. The clones were initially characterized using 39 short tandem repeats (STRs) and 1,536 SNP markers. The utility of whole-genome amplification in preserving and extending RH panel DNA was also tested using 10 STR markers; no significant difference in retention was observed. The resolution of the 15,000(Rad) RH panel was established by constructing framework maps across 10 different 1-Mb regions on different feline chromosomes. In these regions, 2-point analysis was used to estimate RH distances, which compared favorably with the estimation of physical distances. The study demonstrates that the 15,000(Rad) RH panel constitutes a powerful tool for constructing high-resolution maps, having an average resolution of 40.1 kb per marker across the ten 1-Mb regions. In addition, the RH panel will complement existing genomic resources for the domestic cat, aid in the accurate re-assemblies of the forthcoming cat genomic sequence, and support cross-species genomic comparisons.

SSAHA: a fast search method for large DNA databases.

We describe an algorithm, SSAHA (Sequence Search and Alignment by Hashing Algorithm), for performing fast searches on databases containing multiple gigabases of DNA. Sequences in the database are preprocessed by breaking them into consecutive k-tuples of k contiguous bases and then using a hash table to store the position of each occurrence of each k-tuple. Searching for a query sequence in the database is done by obtaining from the hash table the "hits" for each k-tuple in the query sequence and then performing a sort on the results. We discuss the effect of the tuple length k on the search speed, memory usage, and sensitivity of the algorithm and present the results of computational experiments which show that SSAHA can be three to four orders of magnitude faster than BLAST or FASTA, while requiring less memory than suffix tree methods. The SSAHA algorithm is used for high-throughput single nucleotide polymorphism (SNP) detection and very large scale sequence assembly. Also, it provides Web-based sequence search facilities for Ensembl projects.

Cancer and genomics.

Identification of the genes that cause oncogenesis is a central aim of cancer research. We searched the proteins predicted from the draft human genome sequence for paralogues of known tumour suppressor genes, but no novel genes were identified. We then assessed whether it was possible to search directly for oncogenic sequence changes in cancer cells by comparing cancer genome sequences against the draft genome. Apparently chimaeric transcripts (from oncogenic fusion genes generated by chromosomal translocations, the ends of which mapped to different genomic locations) were detected to the same degree in both normal and neoplastic tissues, indicating a significant level of false positives. Our experiment underscores the limited amount and variable quality of DNA sequence from cancer cells that is currently available.

A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms.

We describe a map of 1.42 million single nucleotide polymorphisms (SNPs) distributed throughout the human genome, providing an average density on available sequence of one SNP every 1.9 kilobases. These SNPs were primarily discovered by two projects: The SNP Consortium and the analysis of clone overlaps by the International Human Genome Sequencing Consortium. The map integrates all publicly available SNPs with described genes and other genomic features. We estimate that 60,000 SNPs fall within exon (coding and untranslated regions), and 85% of exons are within 5 kb of the nearest SNP. Nucleotide diversity varies greatly across the genome, in a manner broadly consistent with a standard population genetic model of human history. This high-density SNP map provides a public resource for defining haplotype variation across the genome, and should help to identify biomedically important genes for diagnosis and therapy.

An SNP map of human chromosome 22.

The human genome sequence will provide a reference for measuring DNA sequence variation in human populations. Sequence variants are responsible for the genetic component of individuality, including complex characteristics such as disease susceptibility and drug response. Most sequence variants are single nucleotide polymorphisms (SNPs), where two alternate bases occur at one position. Comparison of any two genomes reveals around 1 SNP per kilobase. A sufficiently dense map of SNPs would allow the detection of sequence variants responsible for particular characteristics on the basis that they are associated with a specific SNP allele. Here we have evaluated large-scale sequencing approaches to obtaining SNPs, and have constructed a map of 2,730 SNPs on human chromosome 22. Most of the SNPs are within 25 kilobases of a transcribed exon, and are valuable for association studies. We have scaled up the process, detecting over 65,000 SNPs in the genome as part of The SNP Consortium programme, which is on target to build a map of 1 SNP every 5 kilobases that is integrated with the human genome sequence and that is freely available in the public domain.

Quantitative DNA fiber mapping.

The assembly of sequence ready, high-resolution physical maps and construction of minimally overlapping contigs for the human as well as model genomes requires accurate determination of the extent of overlap between adjacent clones as well as their relative orientation. This is presently done by procedures such as clone fingerprinting, Southern blot analysis or clone end sequencing. We present a complementary analytical technique to map directly cloned DNA sequences on to individual stretched DNA molecules. This approach uses the hydrodynamic force of a receding meniscus to prepare straight high molecular weight DNA molecules that provide a linear template of approximately 2.3 kb/microns on to which the cloned probes can be mapped by in situ hybridization. This technique has numerous advantages such as a very high density of mapping templates, reproducible stretching of the mapping template providing a linear genomic scale, determination of clone orientation and direct visualization of DNA repeats. The utility and accuracy of quantitative DNA fiber mapping are illustrated through three examples: (i) mapping of lambda DNA restriction fragments along linearized approximately 49 kb long lambda phage DNA molecules with approximately 1 kb precision; (ii) localization of the overlap between a cosmid and a colinear P1 clone; and (iii) mapping of P1 clones along an approximately 490 kb yeast artificial chromosome (YAC) with approximately 5 kb precision and estimation of the approximately 25 kb gap between them.

Dicentric chromosome frequency analysis using slit-scan flow cytometry.

Slit-scan flow cytometry (SSFCM) was used to quantify the frequency of dicentric chromosomes in human lymphoblastoid cells following gamma irradiation. In this study, cultured human cells were irradiated with 0, 0.25, 0.5, 1.0, and 2.0 Gy of 0.66 MeV gamma-rays, cultured for an additional 11 h, and treated for 5 h with colcemid. Chromosomes were then isolated, stained with propidium iodide, and analyzed using SSFCM for total fluorescence and slit-scan profile. The frequency of chromosomes having DNA contents greater than once and less than twice the DNA content of the number 1 chromosome and producing trimodal profiles was determined at each dose. This frequency was used as an estimate of the relative dicentric chromosome frequency at that dose. The estimated dicentric chromosome frequency per cell, f(D), increased with dose, D, in a linear-quadratic manner according to the relation f(D) = 4.52 x 10(-5) + 5.72 x 10(-5) D + 1.19 x 10(-4) D2.

Validation of the saline-dilution method for measuring cardiac output by simultaneous measurement with a perivascular electromagnetic flowprobe.

The validity of measurements of right heart output obtained by the saline-indicator method with a catheter-mounted, tetrapolar resistivity sensor placed in the pulmonary artery was assessed. In ten anesthetized dogs, cardiac output was measured simultaneously by the saline-indicator method, using 5 ml of 3 per cent saline as the indicator, and by a perivascular electromagnetic flowprobe placed around the trunk of the pulmonary artery. Cardiac output was altered over a wide range by increasing the depth of halothane anesthesia and by bolus injections of isoproterenol. Linear regression analysis of cardiac output determined from 455 saline-indicator curves and from simultaneous recordings of the electromagnetic flow signal yielded a slope of 0.973, a Y intercept of -0.0047 1/minute, and a correlation coefficient of 0.965. Hence, this method provides accurate and precise measurement of cardiac output. The saline-indicator method with the catheter-mounted resistivity sensor is simple to use and eliminates or minimizes the disadvantages of other indicator-dilution techniques.