Motile cilia are essential components of the mucociliary escalator and are central to respiratory-tract host defenses. Abnormalities in these evolutionarily conserved organelles cause primary ciliary dyskinesia (PCD). Despite recent strides characterizing the ciliome and sensory ciliopathies through exploration of the phenotype-genotype associations in model organisms, the genetic bases of most cases of PCD remain elusive. We identified nine related subjects with PCD from geographically dispersed Amish communities and performed exome sequencing of two affected individuals and their unaffected parents. A single autosomal-recessive nonsynonymous missense mutation was identified in HEATR2, an uncharacterized gene that belongs to a family not previously associated with ciliary assembly or function. Airway epithelial cells isolated from PCD-affected individuals had markedly reduced HEATR2 levels, absent dynein arms, and loss of ciliary beating. MicroRNA-mediated silencing of the orthologous gene in Chlamydomonas reinhardtii resulted in absent outer dynein arms, reduced flagellar beat frequency, and decreased cell velocity. These findings were recapitulated by small hairpin RNA-mediated knockdown of HEATR2 in airway epithelial cells from unaffected donors. Moreover, immunohistochemistry studies in human airway epithelial cells showed that HEATR2 was localized to the cytoplasm and not in cilia, which suggests a role in either dynein arm transport or assembly. The identification of HEATR2 contributes to the growing number of genes associated with PCD identified in both individuals and model organisms and shows that exome sequencing in family studies facilitates the discovery of novel disease-causing gene mutations.

Acute lymphoblastic leukemia (ALL) is the most common cancer diagnosed in children under the age of 15. In addition to genetic aberrations, epigenetic modifications such as DNA methylation are altered in cancer and impact gene expression. To identify epigenetic alterations in ALL, genome-wide methylation profiles were generated using the methylated CpG island recovery assay followed by next-generation sequencing. More than 25,000 differentially methylated regions (DMR) were observed in ALL patients with »90% present within intronic or intergenic regions. To determine the regulatory potential of the DMR, whole-transcriptome analysis was performed and integrated with methylation data. Aberrant promoter methylation was associated with the altered expression of genes involved in transcriptional regulation, apoptosis, and proliferation. Novel enhancer-like sequences were identified within intronic and intergenic DMR. Aberrant methylation in these regions was associated with the altered expression of neighboring genes involved in cell cycle processes, lymphocyte activation and apoptosis. These genes include potential epi-driver genes, such as SYNE1, PTPRS, PAWR, HDAC9, RGCC, MCOLN2, LYN, TRAF3, FLT1, and MELK, which may provide a selective advantage to leukemic cells. In addition, the differential expression of epigenetic modifier genes, pseudogenes, and non-coding RNAs was also observed accentuating the role of erroneous epigenetic gene regulation in ALL.

BackgroundRare genetic variation in the human population is a major source of pathophysiological variability and has been implicated in a host of complex phenotypes and diseases. Finding disease-related genes harboring disparate functional rare variants requires sequencing of many individuals across many genomic regions and comparing against unaffected cohorts. However, despite persistent declines in sequencing costs, population-based rare variant detection across large genomic target regions remains cost prohibitive for most investigators. In addition, DNA samples are often precious and hybridization methods typically require large amounts of input DNA. Pooled sample DNA sequencing is a cost and time-efficient strategy for surveying populations of individuals for rare variants. We set out to 1) create a scalable, multiplexing method for custom capture with or without individual DNA indexing that was amenable to low amounts of input DNA and 2) expand the functionality of the SPLINTER algorithm for calling substitutions, insertions and deletions across either candidate genes or the entire exome by integrating the variant calling algorithm with the dynamic programming aligner, Novoalign.ResultsWe report methodology for pooled hybridization capture with pre-enrichment, indexed multiplexing of up to 48 individuals or non-indexed pooled sequencing of up to 92 individuals with as little as 70 ng of DNA per person. Modified solid phase reversible immobilization bead purification strategies enable no sample transfers from sonication in 96-well plates through adapter ligation, resulting in 50% less library preparation reagent consumption. Custom Y-shaped adapters containing novel 7 base pair index sequences with a Hamming distance of ≥2 were directly ligated onto fragmented source DNA eliminating the need for PCR to incorporate indexes, and was followed by a custom blocking strategy using a single oligonucleotide regardless of index sequence. These results were obtained aligning raw reads against the entire genome using Novoalign followed by variant calling of non-indexed pools using SPLINTER or SAMtools for indexed samples. With these pipelines, we find sensitivity and specificity of 99.4% and 99.7% for pooled exome sequencing. Sensitivity, and to a lesser degree specificity, proved to be a function of coverage. For rare variants (≤2% minor allele frequency), we achieved sensitivity and specificity of ≥94.9% and ≥99.99% for custom capture of 2.5 Mb in multiplexed libraries of 22–48 individuals with only ≥5-fold coverage/chromosome, but these parameters improved to ≥98.7 and 100% with 20-fold coverage/chromosome.ConclusionsThis highly scalable methodology enables accurate rare variant detection, with or without individual DNA sample indexing, while reducing the amount of required source DNA and total costs through less hybridization reagent consumption, multi-sample sonication in a standard PCR plate, multiplexed pre-enrichment pooling with a single hybridization and lesser sequencing coverage required to obtain high sens...

We detected flavins in the growth medium of the methanotrophic bacterium Methylocystis species strain M. Flavin secretion correlates with growth stage and increases under iron starvation conditions. Two other methanotrophs, Methylosinus trichosporium OB3b and Methylococcus capsulatus (Bath), secrete flavins, suggesting that flavin secretion may be common to many methanotrophic bacteria.Methanotrophic bacteria utilize methane as their sole carbon source and play a key role in the global carbon cycle. The first step in methane metabolism, the oxidation of methane to methanol, is catalyzed by methane monooxygenase (MMO) metalloenzymes by either a copper-dependent membranebound particulate MMO (pMMO) (9) or an iron-dependent, soluble MMO (sMMO) (18). Therefore, metals are necessary for the growth of methanotrophic bacteria (12,13). Despite the importance of metals in methanotroph physiology, uptake of metal ions by methanotrophic bacteria is not well understood. To study metal acquisition, we purified and characterized components of the spent growth medium from Methylocystis species strain M. Unexpectedly, we detected extracellular flavins (riboflavin and flavin mononucleotide [FMN]).To perform these experiments, the type II methanotroph Methylocystis sp. strain M was grown using a procedure similar to that described for Methylosinus trichosporium OB3b (7, 8, 10) (see the supplemental material). The spent growth medium from Methylocystis sp. strain M was purified by fast protein liquid chromatography (FPLC) using a Discovery DSC-18 reverse-phase column. The eluted fractions exhibit UV-visible spectral features that suggest the presence of flavins (Fig. 1A). When the fractions were purified further using a Supelco C 18 column, the purified flavins exhibited absorption features that were identical to those of commercial riboflavin (see Fig. S1 in the supplemental material). Upon excitation at 441 nm, the fluorescence emission spectrum of the purified fractions exhibited a broad peak between 450 and 700 nm, with a maximum at 523 nm (Fig. 1B). This feature is characteristic of the isoalloxazine ring in flavins (11,21).To characterize the amount of extracellular flavins secreted throughout growth, the cells were grown under three metal availability conditions: default growth ( Fig. S2 of the supplemental material. The growth rates, from an average of biological triplicates, determined for default, iron-starved, and copper-starved growth conditions between 12 and 28 h of growth are 0.043 Ϯ 0.003 h Ϫ1 , 0.014 Ϯ 0.003 h Ϫ1 , and 0.052 Ϯ 0.002 h Ϫ1 , respectively. These values imply that the availability of iron affects the growth of Methylocystis sp. strain M.We then characterized the amount of extracellular flavins in the spent medium isolated from three different stages of growth. The samples used for quantitation were normalized based on optical density measured at 600 nm. The dry weight and cell length, normalized to the optical density, remained unaffected by growth under different metal concentrations (see Fig. S4). ...

DNA methylation is responsible for regulating gene expression and cellular differentiation and for maintaining genomic stability during normal human development. Furthermore, it plays a significant role in the regulation of hematopoiesis. In order to elucidate the influence of DNA methylation during B-cell development, genome-wide DNA methylation status of pro-B, pre-BI, pre-BII, and naïve-B-cells isolated from human umbilical cord blood was determined using the methylated CpG island recovery assay followed by next generation sequencing. On average, 182–200 million sequences were generated for each precursor B-cell subset in 10 biological replicates. An overall decrease in methylation was observed during the transition from pro-B to pre-BI, whereas no differential methylation was observed in the pre-BI to pre-BII transition or in the pre-BII to naïve B-cell transition. Most of the methylated regions were located within intergenic and intronic regions not present in a CpG island context. Putative novel enhancers were identified in these regions that were differentially methylated between pro-B and pre-BI cells. The genome-wide methylation profiles are publically available and may be used to gain a better understanding of the involvement of atypical DNA methylation in the pathogenesis of malignancies associated with precursor B-cells.