NTRK fusions and Trk proteins: what are they and how to test for them.

The NTRK genes include a family of three genes, NTRK1, NTRK2, and NTRK3, which are associated with fusions with a variety of partner genes, leading to upregulation of three proteins, TrkA, TrkB, and TrkC. NTRK fusions occur in a variety of solid tumors: at high incidence in secretory carcinoma of the breast and salivary glands, congenital mesoblastic nephroma, and infantile fibrosarcoma; at intermediate incidence in thyroid carcinoma, particularly postradiation carcinomas and a subset of aggressive papillary carcinomas, Spitzoid melanocytic neoplasms, pediatric midline gliomas (particularly pontine glioma), and KIT/PDGFRA/RAS negative gastrointestinal stromal sarcomas; and at a low incidence in many other solid tumors. With new FDA-approved treatments available and effective in treating patients whose tumors harbor NTRK fusions, testing for these fusions has become important. A variety of technologies can be used for testing, including FISH, PCR, DNA, and RNA-based next-generation sequencing, and immunohistochemistry. RNA-based next-generation sequencing represents the gold standard for the identification of NTRK fusions, but FISH using break-apart probes and DNA-based next-generation sequencing also represent adequate approaches. Immunohistochemistry to detect increased levels of Trk protein may be very useful as a screening technology to reduce costs, although it alone does not represent a definitive diagnostic methodology.

A role for BATF3 in TH9 differentiation and T-cell-driven mucosal pathologies.

T helper 9 (TH9) cells are important for the development of inflammatory and allergic diseases. The TH9 transcriptional network converges signals from cytokines and antigen presentation but is incompletely understood. Here, we identified TL1A, a member of the TNF superfamily, as a strong inducer of mouse and human TH9 differentiation. Mechanistically, TL1A induced the expression of the transcription factors BATF and BATF3 and facilitated their binding to the Il9 promoter leading to enhanced secretion of IL-9. BATF- and BATF3-deficiencies impaired IL-9 secretion under TH9 and TH9-TL1A-polarizing conditions. In vivo, using a T-cell transfer model, we demonstrated that TL1A promoted IL-9-dependent, TH9 cell-induced intestinal and lung inflammation. Neutralizing IL-9 antibodies attenuated TL1A-driven mucosal inflammation. Batf3-/- TH9-TL1A cells induced reduced inflammation and cytokine expression in vivo compared to WT cells. Our results demonstrate that TL1A promotes TH9 cell differentiation and function and define a role for BATF3 in T-cell-driven mucosal inflammation.

Genome-wide analysis suggests a differential microRNA signature associated with normal and diabetic human corneal limbus.

Small non-coding RNAs, in particular microRNAs (miRNAs), regulate fine-tuning of gene expression and can impact a wide range of biological processes. However, their roles in normal and diseased limbal epithelial stem cells (LESC) remain unknown. Using deep sequencing analysis, we investigated miRNA expression profiles in central and limbal regions of normal and diabetic human corneas. We identified differentially expressed miRNAs in limbus vs. central cornea in normal and diabetic (DM) corneas including both type 1 (T1DM/IDDM) and type 2 (T2DM/NIDDM) diabetes. Some miRNAs such as miR-10b that was upregulated in limbus vs. central cornea and in diabetic vs. normal limbus also showed significant increase in T1DM vs. T2DM limbus. Overexpression of miR-10b increased Ki-67 staining in human organ-cultured corneas and proliferation rate in cultured corneal epithelial cells. MiR-10b transfected human organ-cultured corneas showed downregulation of PAX6 and DKK1 and upregulation of keratin 17 protein expression levels. In summary, we report for the first time differential miRNA signatures of T1DM and T2DM corneal limbus harboring LESC and show that miR-10b could be involved in the LESC maintenance and/or their early differentiation. Furthermore, miR-10b upregulation may be an important mechanism of corneal diabetic alterations especially in the T1DM patients.

Significant Improvement in Detecting BRAF, KRAS, and EGFR Mutations Using Next-Generation Sequencing as Compared with FDA-Cleared Kits.

INTRODUCTION: We compared mutations detected in EGFR, KRAS, and BRAF genes using next-generation sequencing (NGS) and confirmed by Sanger sequencing with mutations that could be detected by FDA-cleared testing kits. METHODS: Paraffin-embedded tissue from 822 patients was tested for mutations in EGFR, KRAS, and BRAF by NGS. Sanger sequencing of hot spots was used with locked nucleic acid to increase sensitivity for specific hot-spot mutations. This included 442 (54%) lung cancers, 168 (20%) colorectal cancers, 29 (4%) brain tumors, 33 (4%) melanomas, 14 (2%) thyroid cancers, and 16% others (pancreas, head and neck, and cancer of unknown origin). Results were compared with the approved list of detectable mutations in FDA kits for EGFR, KRAS, and BRAF. RESULTS: Of the 101 patients with EGFR abnormalities as detected by NGS, only 58 (57%) were detectable by cobas v2 and only 35 (35%) by therascreen. Therefore, 42 and 65%, respectively, more mutations were detected by NGS, including two patients with EGFR amplification. Of the 117 patients with BRAF mutation detected by NGS, 62 (53%) mutations were within codon 600, detectable by commercial kits, but 55 (47%) of the mutations were outside codon V600, detected by NGS only. Of the 321 patients with mutations in KRAS detected by NGS, 284 (88.5%) had mutations detectable by therascreen and 300 (93.5%) had mutations detectable by cobas. Therefore, 11.5 and 6.5% additional KRAS mutations were detected by NGS, respectively. CONCLUSION: NGS provides significantly more comprehensive testing for mutations as compared with FDA-cleared kits currently available commercially.

Single-cell RNA sequencing identifies diverse roles of epithelial cells in idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a lethal interstitial lung disease characterized by airway remodeling, inflammation, alveolar destruction, and fibrosis. We utilized single-cell RNA sequencing (scRNA-seq) to identify epithelial cell types and associated biological processes involved in the pathogenesis of IPF. Transcriptomic analysis of normal human lung epithelial cells defined gene expression patterns associated with highly differentiated alveolar type 2 (AT2) cells, indicated by enrichment of RNAs critical for surfactant homeostasis. In contrast, scRNA-seq of IPF cells identified 3 distinct subsets of epithelial cell types with characteristics of conducting airway basal and goblet cells and an additional atypical transitional cell that contributes to pathological processes in IPF. Individual IPF cells frequently coexpressed alveolar type 1 (AT1), AT2, and conducting airway selective markers, demonstrating "indeterminate" states of differentiation not seen in normal lung development. Pathway analysis predicted aberrant activation of canonical signaling via TGF-ß, HIPPO/YAP, P53, WNT, and AKT/PI3K. Immunofluorescence confocal microscopy identified the disruption of alveolar structure and loss of the normal proximal-peripheral differentiation of pulmonary epithelial cells. scRNA-seq analyses identified loss of normal epithelial cell identities and unique contributions of epithelial cells to the pathogenesis of IPF. The present study provides a rich data source to further explore lung health and disease.

Immunological Consequences of Intestinal Fungal Dysbiosis.

Compared to bacteria, the role of fungi within the intestinal microbiota is poorly understood. In this study we investigated whether the presence of a "healthy" fungal community in the gut is important for modulating immune function. Prolonged oral treatment of mice with antifungal drugs resulted in increased disease severity in acute and chronic models of colitis, and also exacerbated the development of allergic airway disease. Microbiota profiling revealed restructuring of fungal and bacterial communities. Specifically, representation of Candida spp. was reduced, while Aspergillus, Wallemia, and Epicoccum spp. were increased. Oral supplementation with a mixture of three fungi found to expand during antifungal treatment (Aspergillus amstelodami, Epicoccum nigrum, and Wallemia sebi) was sufficient to recapitulate the exacerbating effects of antifungal drugs on allergic airway disease. Taken together, these results indicate that disruption of commensal fungal populations can influence local and peripheral immune responses and enhance relevant disease states.

Comprehensive assessment of expression of insulin signaling pathway components in subcutaneous adipose tissue of women with and without polycystic ovary syndrome.

OBJECTIVE: Insulin resistance is a common feature of polycystic ovary syndrome (PCOS). The insulin signaling pathway consists of two major pathways, the metabolic and the mitogenic cascades. The many components of these pathways have not been comprehensively analyzed for differential expression in insulin-responsive tissues in PCOS. The goal of this study was to determine whether the core elements of the insulin signal transduction cascade were differentially expressed in subcutaneous adipose tissue (SAT) between PCOS and controls. MATERIALS/METHODS: Quantitative real-time PCR for 36 insulin signaling pathway genes was performed subcutaneous adipose tissue from 22 white PCOS and 13 healthy controls. RESULTS: Genes in the insulin signaling pathway were not differentially expressed in subcutaneous adipose tissue between PCOS and controls (P>0.05 for all). Components mainly of the mitogenic pathway were correlated with both androgens and metabolic phenotypes. Expression levels of five genes (MKNK1, HRAS, NRAS, KRAS, and GSK3A) were positively correlated with total testosterone level (ρ>0, P<0.05). Inverse correlation was found between expression of six genes (HRAS, MAP2K2, NRAS, MAPK3, GRB2, and SHC1) and metabolic traits (body mass index, fasting glucose, fasting insulin, and HOMA-IR) (ρ<0, P<0.05). CONCLUSIONS: Differential expression of core insulin signaling pathway components in subcutaneous adipose tissue is not a major contributor to the pathogenesis of PCOS. Correlation between clinical phenotypes and expression of several genes in the mitogenic limb of the insulin signaling pathway suggests mitogenic signaling by insulin may regulate steroidogenesis and glucose homeostasis.

The development of innate lymphoid cells requires TOX-dependent generation of a common innate lymphoid cell progenitor.

Diverse innate lymphoid cell (ILC) subtypes have been defined on the basis of effector function and transcription factor expression. ILCs derive from common lymphoid progenitors, although the transcriptional pathways that lead to ILC-lineage specification remain poorly characterized. Here we found that the transcriptional regulator TOX was required for the in vivo differentiation of common lymphoid progenitors into ILC lineage-restricted cells. In vitro modeling demonstrated that TOX deficiency resulted in early defects in the survival or proliferation of progenitor cells, as well as ILC differentiation at a later stage. In addition, comparative transcriptome analysis of bone marrow progenitors revealed that TOX-deficient cells failed to upregulate many genes of the ILC program, including genes that are targets of Notch, which indicated that TOX is a key determinant of early specification to the ILC lineage.

Mycobiome: Approaches to analysis of intestinal fungi.

Massively parallel sequencing (MPSS) of bacterial 16S rDNA has been widely used to characterize the microbial makeup of the human and mouse gastrointestinal tract. However, techniques for fungal microbiota (mycobiota) profiling remain relatively under-developed. Compared to 16S profiling, the size and sequence context of the fungal Internal Transcribed Spacer 1 (ITS1), the most common target for mycobiota profiling, are highly variable. Using representative gastrointestinal tract fungi to build a known "mock" library, we examine how this sequence variability affects data quality derived from Illumina Miseq and Ion Torrent PGM sequencing pipelines. Also, while analysis of bacterial 16S profiles is facilitated by the presence of high-quality well-accepted databases of bacterial 16S sequences, such an accepted database has not yet emerged to facilitate fungal ITS sequence characterization, and we observe that redundant and inconsistent ITS1 sequence representation in publically available fungal reference databases affect quantitation and annotation of species in the gut. To address this problem, we have constructed a manually curated reference database optimized for annotation of gastrointestinal fungi. This targeted host-associated fungi (THF) database contains 1817 ITS1 sequences representing sequence diversity in genera previously identified in human and mouse gut. We observe that this database consistently outperforms three common ITS database alternatives on comprehensiveness, taxonomy assignment accuracy and computational efficiency in analyzing sequencing data from the mouse gastrointestinal tract.

Molecular assessment of differences in the duodenal microbiome in subjects with irritable bowel syndrome.

OBJECTIVE: Breath testing and duodenal culture studies suggest that a significant proportion of irritable bowel syndrome (IBS) patients have small intestinal bacterial overgrowth. In this study, we extended these data through 16S rDNA amplicon sequencing and quantitative PCR (qPCR) analyses of duodenal aspirates from a large cohort of IBS, non-IBS and control subjects. MATERIALS AND METHODS: Consecutive subjects presenting for esophagogastroduodenoscopy only and healthy controls were recruited. Exclusion criteria included recent antibiotic or probiotic use. Following extensive medical work-up, patients were evaluated for symptoms of IBS. DNAs were isolated from duodenal aspirates obtained during endoscopy. Microbial populations in a subset of IBS subjects and controls were compared by 16S profiling. Duodenal microbes were then quantitated in the entire cohort by qPCR and the results compared with quantitative live culture data. RESULTS: A total of 258 subjects were recruited (21 healthy, 163 non-healthy non-IBS, and 74 IBS). 16S profiling in five IBS and five control subjects revealed significantly lower microbial diversity in the duodenum in IBS, with significant alterations in 12 genera (false discovery rate < 0.15), including overrepresentation of Escherichia/Shigella (p = 0.005) and Aeromonas (p = 0.051) and underrepresentation of Acinetobacter (p = 0.024), Citrobacter (p = 0.031) and Microvirgula (p = 0.036). qPCR in all 258 subjects confirmed greater levels of Escherichia coli in IBS and also revealed increases in Klebsiella spp, which correlated strongly with quantitative culture data. CONCLUSIONS: 16S rDNA sequencing confirms microbial overgrowth in the small bowel in IBS, with a concomitant reduction in diversity. qPCR supports alterations in specific microbial populations in IBS.

TOX3 is expressed in mammary ER(+) epithelial cells and regulates ER target genes in luminal breast cancer.

BACKGROUND: A breast cancer susceptibility locus has been mapped to the gene encoding TOX3. Little is known regarding the expression pattern or biological role of TOX3 in breast cancer or in the mammary gland. Here we analyzed TOX3 expression in murine and human mammary glands and in molecular subtypes of breast cancer, and assessed its ability to alter the biology of breast cancer cells. METHODS: We used a cell sorting strategy, followed by quantitative real-time PCR, to study TOX3 gene expression in the mouse mammary gland. To study the expression of this nuclear protein in human mammary glands and breast tumors, we generated a rabbit monoclonal antibody specific for human TOX3. In vitro studies were performed on MCF7, BT474 and MDA-MB-231 cell lines to study the effects of TOX3 modulation on gene expression in the context of breast cancer cells. RESULTS: We found TOX3 expression in estrogen receptor-positive mammary epithelial cells, including progenitor cells. A subset of breast tumors also highly expresses TOX3, with poor outcome associated with high expression of TOX3 in luminal B breast cancers. We also demonstrate the ability of TOX3 to alter gene expression in MCF7 luminal breast cancer cells, including cancer relevant genes TFF1 and CXCR4. Knockdown of TOX3 in a luminal B breast cancer cell line that highly expresses TOX3 is associated with slower growth. Surprisingly, TOX3 is also shown to regulate TFF1 in an estrogen-independent and tamoxifen-insensitive manner. CONCLUSIONS: These results demonstrate that high expression of this protein likely plays a crucial role in breast cancer progression. This is in sharp contrast to previous studies that indicated breast cancer susceptibility is associated with lower expression of TOX3. Together, these results suggest two different roles for TOX3, one in the initiation of breast cancer, potentially related to expression of TOX3 in mammary epithelial cell progenitors, and another role for this nuclear protein in the progression of cancer. In addition, these results can begin to shed light on the reported association of TOX3 expression and breast cancer metastasis to the bone, and point to TOX3 as a novel regulator of estrogen receptor-mediated gene expression.

Differentiation of human limbal-derived induced pluripotent stem cells into limbal-like epithelium.

Limbal epithelial stem cell (LESC) deficiency (LSCD) leads to corneal abnormalities resulting in compromised vision and blindness. LSCD can be potentially treated by transplantation of appropriate cells, which should be easily expandable and bankable. Induced pluripotent stem cells (iPSCs) are a promising source of transplantable LESCs. The purpose of this study was to generate human iPSCs and direct them to limbal differentiation by maintaining them on natural substrata mimicking the native LESC niche, including feederless denuded human amniotic membrane (HAM) and de-epithelialized corneas. These iPSCs were generated with nonintegrating vectors from human primary limbal epithelial cells. This choice of parent cells was supposed to enhance limbal cell differentiation from iPSCs by partial retention of parental epigenetic signatures in iPSCs. When the gene methylation patterns were compared in iPSCs to parental LESCs using Illumina global methylation arrays, limbal-derived iPSCs had fewer unique methylation changes than fibroblast-derived iPSCs, suggesting retention of epigenetic memory during reprogramming. Limbal iPSCs cultured for 2 weeks on HAM developed markedly higher expression of putative LESC markers ABCG2, ΔNp63α, keratins 14, 15, and 17, N-cadherin, and TrkA than did fibroblast iPSCs. On HAM culture, the methylation profiles of select limbal iPSC genes (including NTRK1, coding for TrkA protein) became closer to the parental cells, but fibroblast iPSCs remained closer to parental fibroblasts. On denuded air-lifted corneas, limbal iPSCs even upregulated differentiated corneal keratins 3 and 12. These data emphasize the importance of the natural niche and limbal tissue of origin in generating iPSCs as a LESC source with translational potential for LSCD treatment.

Differentially expressed wound healing-related microRNAs in the human diabetic cornea.

MicroRNAs are powerful gene expression regulators, but their corneal repertoire and potential changes in corneal diseases remain unknown. Our purpose was to identify miRNAs altered in the human diabetic cornea by microarray analysis, and to examine their effects on wound healing in cultured telomerase-immortalized human corneal epithelial cells (HCEC) in vitro. Total RNA was extracted from age-matched human autopsy normal (n=6) and diabetic (n=6) central corneas, Flash Tag end-labeled, and hybridized to Affymetrix® GeneChip® miRNA Arrays. Select miRNAs associated with diabetic cornea were validated by quantitative RT-PCR (Q-PCR) and by in situ hybridization (ISH) in independent samples. HCEC were transfected with human pre-miR™miRNA precursors (h-miR) or their inhibitors (antagomirs) using Lipofectamine 2000. Confluent transfected cultures were scratch-wounded with P200 pipette tip. Wound closure was monitored by digital photography. Expression of signaling proteins was detected by immunostaining and Western blot. Using microarrays, 29 miRNAs were identified as differentially expressed in diabetic samples. Two miRNA candidates showing the highest fold increased in expression in the diabetic cornea were confirmed by Q-PCR and further characterized. HCEC transfection with h-miR-146a or h-miR-424 significantly retarded wound closure, but their respective antagomirs significantly enhanced wound healing vs. controls. Cells treated with h-miR-146a or h-miR-424 had decreased p-p38 and p-EGFR staining, but these increased over control levels close to the wound edge upon antagomir treatment. In conclusion, several miRNAs with increased expression in human diabetic central corneas were found. Two such miRNAs inhibited cultured corneal epithelial cell wound healing. Dysregulation of miRNA expression in human diabetic cornea may be an important mediator of abnormal wound healing.

Interactions between commensal fungi and the C-type lectin receptor Dectin-1 influence colitis.

The intestinal microflora, typically equated with bacteria, influences diseases such as obesity and inflammatory bowel disease. Here, we show that the mammalian gut contains a rich fungal community that interacts with the immune system through the innate immune receptor Dectin-1. Mice lacking Dectin-1 exhibited increased susceptibility to chemically induced colitis, which was the result of altered responses to indigenous fungi. In humans, we identified a polymorphism in the gene for Dectin-1 (CLEC7A) that is strongly linked to a severe form of ulcerative colitis. Together, our findings reveal a eukaryotic fungal community in the gut (the "mycobiome") that coexists with bacteria and substantially expands the repertoire of organisms interacting with the intestinal immune system to influence health and disease.

Dominant and recessive forms of fibrochondrogenesis resulting from mutations at a second locus, COL11A2.

Fibrochondrogenesis is a severe, recessively inherited skeletal dysplasia shown to result from mutations in the gene encoding the proα1(XI) chain of type XI collagen, COL11A1. The first of two cases reported here was the affected offspring of first cousins and sequence analysis excluded mutations in COL11A1. Consequently, whole-genome SNP genotyping was performed to identify blocks of homozygosity, identical-by-descent, wherein the disease locus would reside. COL11A1 was not within a region of homozygosity, further excluding it as the disease locus, but the gene encoding the proα2(XI) chain of type XI collagen, COL11A2, was located within a large region of homozygosity. Sequence analysis identified homozygosity for a splice donor mutation in intron 18. Exon trapping demonstrated that the mutation resulted in skipping of exon 18 and predicted deletion of 18 amino acids from the triple helical domain of the protein. In the second case, heterozygosity for a de novo 9 bp deletion in exon 40 of COL11A2 was identified, indicating that there are autosomal dominant forms of fibrochondrogenesis. These findings thus demonstrate that fibrochondrogenesis can result from either recessively or dominantly inherited mutations in COL11A2.

The importance of conventional radiography in the mutational analysis of skeletal dysplasias (the TRPV4 mutational family).

The spondylo and spondylometaphyseal dysplasias (SMDs) are characterized by vertebral changes and metaphyseal abnormalities of the tubular bones, which produce a phenotypic spectrum of disorders from the mild autosomal-dominant brachyolmia to SMD Kozlowski to autosomal-dominant metatropic dysplasia. Investigations have recently drawn on the similar radiographic features of those conditions to define a new family of skeletal dysplasias caused by mutations in the transient receptor potential cation channel vanilloid 4 (TRPV4). This review demonstrates the significance of radiography in the discovery of a new bone dysplasia family due to mutations in a single gene.

BMPER mutation in diaphanospondylodysostosis identified by ancestral autozygosity mapping and targeted high-throughput sequencing.

Diaphanospondylodysostosis (DSD) is a rare, recessively inherited, perinatal lethal skeletal disorder. The low frequency and perinatal lethality of DSD makes assembling a large set of families for traditional linkage-based genetic approaches challenging. By searching for evidence of unknown ancestral consanguinity, we identified two autozygous intervals, comprising 34 Mbps, unique to a single case of DSD. Empirically testing for ancestral consanguinity was effective in localizing the causative variant, thereby reducing the genomic space within which the mutation resides. High-throughput sequence analysis of exons captured from these intervals demonstrated that the affected individual was homozygous for a null mutation in BMPER, which encodes the bone morphogenetic protein-binding endothelial cell precursor-derived regulator. Mutations in BMPER were subsequently found in three additional DSD cases, confirming that defects in BMPER produce DSD. Phenotypic similarities between DSD and Bmper null mice indicate that BMPER-mediated signaling plays an essential role in vertebral segmentation early in human development.

Quantitative gene expression profiles in real time from expressed sequence tag databases.

An accumulation of expressed sequence tag (EST) data in the public domain and the availability of bioinformatic programs have made EST gene expression profiling a common practice. However, the utility and validity of using EST databases (e.g., dbEST) has been criticized, particularly for quantitative assessment of gene expression. Problems with EST sequencing errors, library construction, EST annotation, and multiple paralogs make generation of specific and sensitive qualitative arid quantitative expression profiles a concern. In addition, most EST-derived expression data exists in previously assembled databases. The Virtual Northern Blot (VNB) (http: //tlab.bu.edu/vnb.html) allows generation, evaluation, and optimization of expression profiles in real time, which is especially important for alternatively spliced, novel, or poorly characterized genes. Representative gene families with variable nucleotide sequence identity, tissue specificity, and levels of expression (bcl-xl, aldoA, and cyp2d9) are used to assess the quality of VNB's output. The profiles generated by VNB are more sensitive and specific than those constructed with ESTs listed in preindexed databases at UCSC and NCBI. Moreover, quantitative expression profiles produced by VNB are comparable to quantization obtained from Northern blots and qPCR. The VNB pipeline generates real-time gene expression profiles for single-gene queries that are both qualitatively and quantitatively reliable.

Ciliary abnormalities due to defects in the retrograde transport protein DYNC2H1 in short-rib polydactyly syndrome.

The short-rib polydactyly (SRP) syndromes are a heterogeneous group of perinatal lethal skeletal disorders with polydactyly and multisystem organ abnormalities. Homozygosity by descent mapping in a consanguineous SRP family identified a genomic region that contained DYNC2H1, a cytoplasmic dynein involved in retrograde transport in the cilium. Affected individuals in the family were homozygous for an exon 12 missense mutation that predicted the amino acid substitution R587C. Compound heterozygosity for one missense and one null mutation was identified in two additional nonconsanguineous SRP families. Cultured chondrocytes from affected individuals showed morphologically abnormal, shortened cilia. In addition, the chondrocytes showed abnormal cytoskeletal microtubule architecture, implicating an altered microtubule network as part of the disease process. These findings establish SRP as a cilia disorder and demonstrate that DYNC2H1 is essential for skeletogenesis and growth.

Improving the efficiency of genomic loci capture using oligonucleotide arrays for high throughput resequencing.

BACKGROUND: The emergence of next-generation sequencing technology presents tremendous opportunities to accelerate the discovery of rare variants or mutations that underlie human genetic disorders. Although the complete sequencing of the affected individuals' genomes would be the most powerful approach to finding such variants, the cost of such efforts make it impractical for routine use in disease gene research. In cases where candidate genes or loci can be defined by linkage, association, or phenotypic studies, the practical sequencing target can be made much smaller than the whole genome, and it becomes critical to have capture methods that can be used to purify the desired portion of the genome for shotgun short-read sequencing without biasing allelic representation or coverage. One major approach is array-based capture which relies on the ability to create a custom in-situ synthesized oligonucleotide microarray for use as a collection of hybridization capture probes. This approach is being used by our group and others routinely and we are continuing to improve its performance. RESULTS: Here, we provide a complete protocol optimized for large aggregate sequence intervals and demonstrate its utility with the capture of all predicted amino acid coding sequence from 3,038 human genes using 241,700 60-mer oligonucleotides. Further, we demonstrate two techniques by which the efficiency of the capture can be increased: by introducing a step to block cross hybridization mediated by common adapter sequences used in sequencing library construction, and by repeating the hybridization capture step. These improvements can boost the targeting efficiency to the point where over 85% of the mapped sequence reads fall within 100 bases of the targeted regions. CONCLUSIONS: The complete protocol introduced in this paper enables researchers to perform practical capture experiments, and includes two novel methods for increasing the targeting efficiency. Coupled with the new massively parallel sequencing technologies, this provides a powerful approach to identifying disease-causing genetic variants that can be localized within the genome by traditional methods.