Endomicrobiome of in vitro and natural plants deciphering the endophytes-associated secondary metabolite biosynthesis in Picrorhiza kurrooa, a Himalayan medicinal herb

IMPORTANCE

Picrorhiza kurrooa is a major source of picrosides, potent hepatoprotective molecules. Due to the ever-increasing demands, overexploitation has caused an extensive decline in its population in the wild and placed it in the endangered plants' category. At present plant in-vitro systems are widely used for the sustainable generation of P. kurrooa plants, and also for the conservation of other commercially important, rare, endangered, and threatened plant species. Furthermore, the in-vitro-generated plants had reduced content of therapeutic secondary metabolites compared to their wild counterparts, and the reason behind, not well-explored. Here, we revealed the loss of plant-associated endophytic communities during in-vitro propagation of P. kurrooa plants which also correlated to in-planta secondary metabolite biosynthesis. Therefore, this study emphasized to consider the essential role of plant-associated endophytic communities in in-vitro practices which may be the possible reason for reduced secondary metabolites in in-vitro plants.

Technical feasibility of newborn screening for spinal muscular atrophy by next-generation DNA sequencing

Newborn screening (NBS) assays for spinal muscular atrophy (SMA) typically use a polymerase chain reaction (PCR) based assay to identify individuals with homozygous deletion in exon 7 of the SMN1 gene. Due to high DNA sequence homology between SMN1 and SMN2, it has previously been difficult to accurately bioinformatically map short reads from next-generation DNA sequencing (NGS) to SMN1, resulting in low analytical performance and preventing NGS being used for SMA screening. Advances in bioinformatics have allowed NGS to be used in diagnostic settings, but to date these assays have not reached the scale required for high volume population newborn screening and have not been performed on the dried blood spot samples that NBS programs currently use. Here we integrate an NGS assay using hybridisation-based capture with a customised bioinformatics algorithm and purpose designed high throughput reporting software into an existing NBS program to achieve a laboratory workflow for population SMA screening. We tested the NGS assay on over 2500 newborns born over 2 weeks in a NBS program in a technical feasibility study and show high sensitivity and specificity. Our results suggest NGS may be an alternate method for SMA screening by NBS programs, providing a multiplex testing platform on which potentially hundreds of inherited conditions could be simultaneously tested.

Proteomics and AQP2 regulation

The advent of modern quantitative protein mass spectrometry techniques around the turn of the 21st century has contributed to a revolution in biology referred to as 'systems biology'. These methods allow identification and quantification of thousands of proteins in a biological specimen, as well as detection and quantification of post-translational protein modifications including phosphorylation. Here, we discuss these methodologies and show how they can be applied to understand the effects of the peptide hormone vasopressin to regulate the molecular water channel aquaporin-2. The emerging picture provides a detailed framework for understanding the molecular mechanisms involved in water balance disorders.

The impact of an enhanced infection control protocol on the microbial community profile in molar root canal treatment-an in vivo NGS molecular study

INTRODUCTION

Recent findings demonstrated that one-year CBCT-based outcomes of molar root canal treatment were improved through an Enhanced infection protocol (EnP), when compared to a current best-practice standard infection control protocol (StP). The EnP comprised measures to reduce iatrogenic contamination from direct and indirect contact surfaces, including the replacement of the rubber dams, gloves, files, all instruments, and surface barriers before root canal obturation. The aim of this study was to investigate the effect of such an enhanced infection control protocol on resident microbiome present after chemomechanical instrumentation and the protocol ability in reducing iatrogenic contamination in molar teeth during root canal treatment.

METHODS

Molar teeth were block-randomized to receive either treatment under EnP or StP. To compare the differential effect of the protocol on the identity of bacteria present, one hundred and fifty, matched DNA extracts from 75 molar teeth samples (StP, n=39; EnP, n=36), were evaluated. Samples were taken before (S1) and after (S2) chemomechanical preparation and were subjected to next-generation sequencing of the V3-V4 region of the 16S rRNA gene, prior to bioinformatical identification using the HOMD oral microbiome database and downstream taxonomic processing, providing measures of richness and diversity of bacteria and significant bacterial taxa during chemomechanical instrumentation and the effect of the two treatment groups.

RESULTS

88 microbial taxa were significantly more abundant in StP S2 samples, including endodontically relevant contaminants taxa as Actinomyces, Cutibacterium, and Haemophilus. The S2 samples demonstrated fewer residual bacterial species in the EnP group, with 26.8 observed species compared to 38.3 in the StP. Reduced diversity and richness measures were noted in the EnP pre-obturation samples compared to the StP in OTU, Chao1 and ACE indices (p≤0.05). Differential microbial identities between S1 and S2 samples and protocols demonstrated that the previously observed increased effectiveness of the EnP protocol was likely to prevent recontamination or de novo contamination of the root canal space during treatment.

CONCLUSIONS

The implemented enhanced infection control protocol resulted in a specific reduction of microbial taxa often associated with recontamination or iatrogenic contamination, suggesting the basis for improved infection control measures during root canal treatment.

Systems Biology of the Vasopressin V2 Receptor: New Tools for Discovery of Molecular Actions of a GPCR

Systems biology can be defined as the study of a biological process in which all of the relevant components are investigated together in parallel to discover the mechanism. Although the approach is not new, it has come to the forefront as a result of genome sequencing projects completed in the first few years of the current century. It has elements of large-scale data acquisition (chiefly next-generation sequencing-based methods and protein mass spectrometry) and large-scale data analysis (big data integration and Bayesian modeling). Here we discuss these methodologies and show how they can be applied to understand the downstream effects of GPCR signaling, specifically looking at how the neurohypophyseal peptide hormone vasopressin, working through the V2 receptor and PKA activation, regulates the water channel aquaporin-2. The emerging picture provides a detailedframework for understanding the molecular mechanisms involved in water balance disorders, pointing the way to improved treatment of both polyuric disorders and water-retention disorders causing dilutional hyponatremia.

Phage display broadly identifies inhibitor-reactive regions in von Willebrand factor

BACKGROUND

Correction of von Willebrand factor (VWF) deficiency with replacement products containing VWF can lead to the development of anti-VWF alloantibodies (i.e., VWF inhibitors) in patients with severe von Willebrand disease (VWD).

OBJECTIVE

Locate inhibitor-reactive regions within VWF using phage display.

METHODS

We screened a phage library displaying random, overlapping fragments covering the full-length VWF protein sequence for binding to a commercial anti-VWF antibody or to immunoglobulins from three type 3 VWD patients who developed VWF inhibitors in response to treatment with plasma-derived VWF. Immunoreactive phage clones were identified and quantified by next-generation DNA sequencing (NGS).

RESULTS

Next-generation DNA sequencing markedly increased the number of phages analyzed for locating immunoreactive regions within VWF following a single round of selection and identified regions not recognized in previous reports using standard phage display methods. Extending this approach to characterize VWF inhibitors from three type 3 VWD patients (including two siblings homozygous for the same VWF gene deletion) revealed patterns of immunoreactivity distinct from the commercial antibody and between unrelated patients, though with notable areas of overlap. Alloantibody reactivity against the VWF propeptide is consistent with incomplete removal of the propeptide from plasma-derived VWF replacement products.

CONCLUSION

These results demonstrate the utility of phage display and NGS to characterize diverse anti-VWF antibody reactivities.

Characterization of Oral Microbiota in Cats: Novel Insights on the Potential Role of Fungi in Feline Chronic Gingivostomatitis

Previous studies have suggested the involvement of viral and bacterial components in the initiation and progression of feline chronic gingivostomatitis (FCGS), but the role of fungi remains entirely unknown. This pilot study aimed to investigate the bacteriome and mycobiome in feline oral health and disease. Physical exams, including oral health assessment, of privately owned, clinically healthy (CH) cats (n = 14) and cats affected by FCGS (n = 14) were performed. Using a sterile swab, oral tissue surfaces of CH and FCGS cats were sampled and submitted for 16S rRNA and ITS-2 next-generation DNA sequencing. A high number of fungal species (n = 186) was detected, with Malassezia restricta, Malassezia arunalokei, Cladosporium penidielloides/salinae, and Aspergillaceae sp. being significantly enriched in FCGS samples, and Saccharomyces cerevisiae in CH samples. The bacteriome was significantly distinct between groups, and significant inter-kingdom interactions were documented. Bergeyella zoohelcum was identified as a potential biomarker of a healthy feline oral microbiome. These data suggest that fungi might play a role in the etiology and pathogenesis of FCGS, and that oral health should not simply be regarded as the absence of microbial infections. Instead, it may be viewed as the biological interactions between bacterial and fungal populations that coexist to preserve a complex equilibrium in the microenvironment of the mouth. Additional investigations are needed to improve our understanding of the feline oral ecosystem and the potential interactions between viruses, bacteria, and fungi in FCGS.

Identifying Optimal Loci for the Molecular Diagnosis of Microsatellite Instability

BACKGROUND

Microsatellite instability (MSI) predicts oncological response to checkpoint blockade immunotherapies. Although microsatellite mutation is pathognomonic for the condition, loci have unequal diagnostic value for predicting MSI within and across cancer types.

METHODS

To better inform molecular diagnosis of MSI, we examined 9438 tumor-normal exome pairs and 901 whole genome sequence pairs from 32 different cancer types and cataloged genome-wide microsatellite instability events. Using a statistical framework, we identified microsatellite mutations that were predictive of MSI within and across cancer types. The diagnostic accuracy of different subsets of maximally informative markers was estimated computationally using a dedicated validation set.

RESULTS

Twenty-five cancer types exhibited hypermutated states consistent with MSI. Recurrently mutated microsatellites associated with MSI were identifiable in 15 cancer types, but were largely specific to individual cancer types. Cancer-specific microsatellite panels of 1 to 7 loci were needed to attain ≥95% diagnostic sensitivity and specificity for 11 cancer types, and in 8 of the cancer types, 100% sensitivity and specificity were achieved. Breast cancer required 800 loci to achieve comparable performance. We were unable to identify recurrent microsatellite mutations supporting reliable MSI diagnosis in ovarian tumors. Features associated with informative microsatellites were cataloged.

CONCLUSIONS

Most microsatellites informative for MSI are specific to particular cancer types, requiring the use of tissue-specific loci for optimal diagnosis. Limited numbers of markers are needed to provide accurate MSI diagnosis in most tumor types, but it is challenging to diagnose breast and ovarian cancers using predefined microsatellite locus panels.

Congenital defects of glycosylation: Novel presentations with mainly neurological involvement and variable dysmorphic features

The pathophysiology of congenital defects of glycosylation (CDG) is complex and the diagnosis has been a challenge because of the overlapping clinical signs and symptoms as well as a large number of disorders. Isoelectric focusing of transferrin has been used as a screening method but has limitations. Individual enzyme or molecular genetic tests have been difficult to perform. In this study, we aimed to describe CDG patients who were referred to from different departments either without a preliminary diagnosis or suspected to have a genetic disorder other than CDG. The patients were diagnosed mainly with a 450 gene next-generation DNA sequencing panel for inborn errors of metabolism, which also included 25 genes for CDG. A total of 862 patients were investigated with the panel, whereby homozygous (10) or compound heterozygous (4) mutations were found in a total of 14 (1.6%) patients. A total of 13 different mutations were discovered, 10 of them being novel. Interestingly, none of the patients was suspected to have a CDG before referral. This report expands the clinical/laboratory findings in patients with CDG and stresses on the fact that CDG should be in the differential list for pediatric patients presented with nonspecific dysmorphic features and neurological delays/regression. Also, next-generation DNA sequencing with panel approach was noticed to have a significant diagnostic potential in patients presented with nonspecific neurologic and dysmorphic findings.

EGFRvIII expression and isocitrate dehydrogenase mutations in patients with glioma

Molecular pathology and personalized medicine are still being evolved in Saudi Arabia, and genetic testing for the detection of mutations as cancer markers have not been established in the diagnostics laboratories in Saudi Arabia. The aim of the present study was to determine the prevalence of isocitrate dehydrogenase (IDH1 and IDH2) mutations and epidermal growth factor receptor variant (EGFRv)III transcript expression in Saudi Arabian patients with glioma. Out of 117 brain tumors tested by reverse transcription-quantitative PCR for EGFRvIII, 41 cases tested positive. In the glioblastoma (GBM) category, 28/55 tumors were positive, in astrocytoma tumors 5/22, and in oligodendrogliomas 4/13 cases were positive respectively. EGFRvIII transcript was sequenced by capillary electrophoresis to demonstrate the presence of EGFRvIII-specific junction where exons 2–7 were deleted. In the present study 106 tumors were sequenced for IDH1 exon-4 mutations using the capillary sequencing method. The most common substitution missense mutation c.395G>A was found in 16 tumors. In the case of adamantinomatous craniopharyngioma, a novel missense mutation in c.472C>T was detected in IDH2 gene. Using next-generation sequencing (NGS), 74 tumors were sequenced for the IDH1 gene, and a total of 8 missense variants were identified in 36 tumors in a population of Saudi Arabia. The missense mutation (c.395G>A) was detected in 29/36 of tumors. A novel intronic mutation in c.414+9T>A was found in 13 cases in the IDH1 gene. In addition, one case exhibited a novel synonymous mutation in c.369A>G. Eleven tumors were found to have compound mutations in the IDH1 gene. In IDH2 gene, out of a total of 16 variants found in 6 out of 45 tumors, nine were missense, five were synonymous and one was intronic. This is the first report from Saudi Arabian laboratories analyzing glioma tumors for EGFRvIII expression, and the first study from Saudi Arabia to analyze IDH mutations in gliomas using the capillary and NGS methods.

Cultivating DNA Sequencing Technology After the Human Genome Project

When the Human Genome Project was completed in 2003, automated Sanger DNA sequencing with fluorescent dye labels was the dominant technology. Several nascent alternative methods based on older ideas that had not been fully developed were the focus of technical researchers and companies. Funding agencies recognized the dynamic nature of technology development and that, beyond the Human Genome Project, there were growing opportunities to deploy DNA sequencing in biological research. Consequently, the National Human Genome Research Institute of the National Institutes of Health created a program-widely known as the Advanced Sequencing Technology Program-that stimulated all stages of development of new DNA sequencing methods, from innovation to advanced manufacturing and production testing, with the goal of reducing the cost of sequencing a human genome first to $100,000 and then to $1,000. The events of this period provide a powerful example of how judicious funding of academic and commercial partners can rapidly advance core technology developments that lead to profound advances across the scientific landscape.

Oxford nanopore sequencing enables rapid discovery of single-domain antibodies from phage display libraries

Antibody (Ab) repertoire sequencing using high-throughput massively parallel technologies has contributed substantially to the understanding of Ab responses following infection, vaccination and autoimmunity. Because individual B-cell receptors are recombined and diversified somatically, genomic comparisons are limited, and distinguishing rare variants from sequencing errors is a major challenge. Oxford Nanopore Technologies' MinION is a highly portable and cost-effective third-generation sequencing instrument, but has not been used for Ab repertoire sequencing due to its high error rate (approximately 1/10 bases). Here, we applied nanopore sequencing to single-domain Ab (sdAb) repertoires and phage-displayed sdAb libraries. We show that despite low overall data fidelity, sdAb sequences could be reconstructed above a frequency threshold (∼100 copies); however, distinguishing clonal sdAb variants was not always possible. The data quality was sufficient to enable rapid identification of antigen-specific sdAb sequences enriched during panning of phage display libraries, obviating the need for screening single clones.

polyRAD: Genotype Calling with Uncertainty from Sequencing Data in Polyploids and Diploids

Low or uneven read depth is a common limitation of genotyping-by-sequencing (GBS) and restriction site-associated DNA sequencing (RAD-seq), resulting in high missing data rates, heterozygotes miscalled as homozygotes, and uncertainty of allele copy number in heterozygous polyploids. Bayesian genotype calling can mitigate these issues, but previously has only been implemented in software that requires a reference genome or uses priors that may be inappropriate for the population. Here we present several novel Bayesian algorithms that estimate genotype posterior probabilities, all of which are implemented in a new R package, polyRAD. Appropriate priors can be specified for mapping populations, populations in Hardy-Weinberg equilibrium, or structured populations, and in each case can be informed by genotypes at linked markers. The polyRAD software imports read depth from several existing pipelines, and outputs continuous or discrete numerical genotypes suitable for analyses such as genome-wide association and genomic prediction.

Deep sequencing of the T cell receptor β repertoire reveals signature patterns and clonal drift in atherosclerotic plaques and patients

The T cell receptor (TCR) β repertoire directly reflects the status of T cell function. Meanwhile, the immune/inflammatory responses regulated by T cells are the critical determinants of atherosclerosis development. However, due to technical limitations, the composition and molecular characteristics of the TCR repertoire in atherosclerotic patients have not been fully elucidated. In the present study, we use powerful immune repertoire sequencing technology to study this issue. Results show that the utilization of V and/or J genes and the diversity of TCRβ repertoire in atherosclerotic plaques are significantly reduced compared to those in the peripheral blood of normal subjects and atherosclerotic patients. The frequencies of the common T cell clones with certain lengths of the complement determining region 3 regions are notably different among all groups. The high-frequency common clones are also increased in the atherosclerotic plaques compared to that in the other two groups. The expansion of several T cell clonotypes (V29-1J2-1, V20-1J1-6, V6-3J2-7 and V11-2J2-2) is validated in atherosclerotic patients. In short, this study reveals that the diversity of TCR β repertoire significantly decreases in atherosclerotic plaques, probably because of the reduced utilization of VJ genes and marked expansion of some T cell subclones. It provides the basis for understanding the roles of T lymphocytes in the pathogenesis of atherosclerosis.

Studying the human oral microbiome: challenges and the evolution of solutions

Since the pioneering work of van Leeuwenhoek in 1684, subsequently built upon by other renowned microbiologists Robert Koch, Willoughby Miller and GV Black, oral microbiology has developed innovative techniques to study the oral microflora (now termed the 'oral microbiome'). The advent of molecular techniques such as DNA-DNA hybridization, polymerase chain reaction and DNA sequencing has created an array of opportunities to construct a comprehensive picture of the diversity and composition of the oral microbiome. Approximately 700 oral bacterial species have been identified, of which 50% have yet to be cultivated, and some of these are known only by their signature DNA sequences. The synergism of ever-evolving culture-based and state-of-the-art culture-independent molecular techniques has facilitated in-depth understanding of the dynamics, acquisition and transfer of oral bacteria, along with their role in oral and general health and disease. Further research is needed to not only analyse but also to make sense of the ever-increasing volumes of data which these molecular techniques (especially high-throughput DNA sequencing) are generating, as well as why particular bacteria are present and what they are 'actually doing' there. This review presents a comprehensive literature search of oral microbiology-related methods currently used to study the oral microbiome.

Clinical, immunologic, and genetic spectrum of 696 patients with combined immunodeficiency

BACKGROUND

Combined immunodeficiencies (CIDs) are diseases of defective adaptive immunity with diverse clinical phenotypes. Although CIDs are more prevalent in the Middle East than Western countries, the resources for genetic diagnosis are limited.

OBJECTIVES

This study aims to characterize the categories of patients with CIDs in Iran clinically and genetically.

METHODS

Clinical and laboratory data were obtained from 696 patients with CIDs. Patients were subdivided into those with syndromic (344 patients) and nonsyndromic (352 patients) CIDs. Targeted DNA sequencing was performed on 243 (34.9%) patients.

RESULTS

The overall diagnostic yield of the 243 sequenced patients was 77.8% (189 patients). The clinical diagnosis of hyper-IgE syndrome (P < .001), onset of disease at greater than 5 years (P = .02), and absence of multiple affected family members (P = .04) were significantly more frequent in the patients without a genetic diagnosis. An autosomal recessive disease was found in 62.9% of patients, reflecting the high rate of consanguinity in this cohort. Mutations impairing VDJ recombination and DNA repair were the most common underlying causes of CIDs. However, in patients with syndromic CIDs, autosomal recessive mutations in ataxia-telangiectasia mutated (ATM), autosomal dominant mutations in signal transducer and activator of transcription 3 (STAT3), and microdeletions in 22q11.21 were the most commonly affected genomic loci. Patients with syndromic CIDs had a significantly lower 5-year survival rate rather than those with nonsyndromic CIDs.

CONCLUSIONS

This study provides proof of principle for the application of targeted next-generation sequencing panels in countries with limited diagnostic resources. The effect of genetic diagnosis on clinical care requires continued improvements in therapeutic resources for these patients.

Hyaluronan-binding peptide for targeting peritoneal carcinomatosis

Peritoneal carcinomatosis results from dissemination of solid tumors in the peritoneal cavity, and is a common site of metastasis in patients with carcinomas of gastrointestinal or gynecological origin. Peritoneal carcinomatosis treatment is challenging as poorly vascularized, disseminated peritoneal micro-tumors are shielded from systemic anticancer drugs and drive tumor regrowth. Here, we describe the identification and validation of a tumor homing peptide CKRDLSRRC (IP3), which upon intraperitoneal administration delivers payloads to peritoneal metastases. IP3 peptide was identified by in vivo phage display on a mouse model of peritoneal carcinomatosis of gastric origin (MKN-45P), using high-throughput sequencing of the peptide-encoding region of phage genome as a readout. The IP3 peptide contains a hyaluronan-binding motif, and fluorescein-labeled IP3 peptide bound to immobilized hyaluronan in vitro. After intraperitoneal administration in mice bearing peritoneal metastases of gastric and colon origin, IP3 peptide homed robustly to macrophage-rich regions in peritoneal tumors, including poorly vascularized micro-tumors. Finally, we show that IP3 functionalization conferred silver nanoparticles the ability to home to peritoneal tumors of gastric and colonic origin, suggesting that it could facilitate targeted delivery of nanoscale payloads to peritoneal tumors. Collectively, our study suggests that the IP3 peptide has potential applications for targeting drugs, nanoparticles, and imaging agents to peritoneal tumors.

Identification of novel proteins binding the AU-rich element of α-prothymosin mRNA through the selection of open reading frames (RIDome)

We describe here a platform for high-throughput protein expression and interaction analysis aimed at identifying the RNA-interacting domainome. This approach combines the selection of a phage library displaying "filtered" open reading frames with next-generation DNA sequencing. The method was validated using an RNA bait corresponding to the AU-rich element of α-prothymosin, an RNA motif that promotes mRNA stability and translation through its interaction with the RNA-binding protein ELAVL1. With this strategy, we not only confirmed known RNA-binding proteins that specifically interact with the target RNA (such as ELAVL1/HuR and RBM38) but also identified proteins not previously known to be ARE-binding (R3HDM2 and RALY). We propose this technology as a novel approach for studying the RNA-binding proteome.

Peptide Synthesis on a Next-Generation DNA Sequencing Platform

Methods for displaying large numbers of peptides on solid surfaces are essential for high-throughput characterization of peptide function and binding properties. Here we describe a method for converting the >10(7) flow cell-bound clusters of identical DNA strands generated by the Illumina DNA sequencing technology into clusters of complementary RNA, and subsequently peptide clusters. We modified the flow-cell-bound primers with ribonucleotides thus enabling them to be used by poliovirus polymerase 3D(pol) . The primers hybridize to the clustered DNA thus leading to RNA clusters. The RNAs fold into functional protein- or small molecule-binding aptamers. We used the mRNA-display approach to synthesize flow-cell-tethered peptides from these RNA clusters. The peptides showed selective binding to cognate antibodies. The methods described here provide an approach for using DNA clusters to template peptide synthesis on an Illumina flow cell, thus providing new opportunities for massively parallel peptide-based assays.

Identification of cross-reactive single-domain antibodies against serum albumin using next-generation DNA sequencing

Antibodies that cross-react with multiple isoforms or homologue of a given protein are often desirable, especially in therapeutic applications. Here, we report the identification of several unique, clonally unrelated, single-domain antibodies (sdAbs) that bind to multiple serum albumin orthologues (human, rhesus, rat and mouse) with low-to-medium nanomolar affinity from a llama immunized only with human serum albumin. Using single-round panning of a phage-displayed sdAb library against serum albumins and next-generation DNA sequencing, we were able to predict patterns of sdAb reactivity to the albumins of different species with ∼90% accuracy. We expect this strategy to be generally applicable for identifying cross-reactive sdAbs, particularly when these exist at low frequency and/or are poorly enriched by panning. Moreover, the sdAbs identified here are of potential interest for serum half-life extension of biologics.

A rigorous approach for selection of optimal variant sets for carrier screening with demonstration of clinical utility

Carrier screening for certain diseases is recommended by major medical and Ashkenazi Jewish (AJ) societies. Most carrier screening panels test only for common, ethnic-specific variants. However, with formerly isolated ethnic groups becoming increasingly intermixed, this approach is becoming inadequate. Our objective was to develop a rigorous process to curate all variants, for relevant genes, into a database and then apply stringent clinical validity classification criteria to each in order to retain only those with clear evidence for pathogenicity. The resulting variant set, in conjunction with next-generation DNA sequencing (NGS), then affords the capability for an ethnically diverse, comprehensive, highly specific carrier-screening assay. The clinical utility of our approach was demonstrated by screening a pan-ethnic population of 22,864 individuals for Bloom syndrome carrier status using a BLM variant panel comprised of 50 pathogenic variants. In addition to carriers of the common AJ founder variant, we identified 57 carriers of other pathogenic BLM variants. All variants reported had previously been curated and their clinical validity documented, or were of a type that met our stringent, preassigned validity criteria. Thus, it was possible to confidently report an increased number of Bloom’s syndrome carriers compared to traditional, ethnicity-based screening, while not reducing the specificity of the screening due to reporting variants of unknown clinical significance.

Fast and cost-effective genetic mapping in apple using next-generation sequencing

Next-generation DNA sequencing (NGS) produces vast amounts of DNA sequence data, but it is not specifically designed to generate data suitable for genetic mapping. Recently developed DNA library preparation methods for NGS have helped solve this problem, however, by combining the use of reduced representation libraries with DNA sample barcoding to generate genome-wide genotype data from a common set of genetic markers across a large number of samples. Here we use such a method, called genotyping-by-sequencing (GBS), to produce a data set for genetic mapping in an F1 population of apples (Malus × domestica) segregating for skin color. We show that GBS produces a relatively large, but extremely sparse, genotype matrix: over 270,000 SNPs were discovered but most SNPs have too much missing data across samples to be useful for genetic mapping. After filtering for genotype quality and missing data, only 6% of the 85 million DNA sequence reads contributed to useful genotype calls. Despite this limitation, using existing software and a set of simple heuristics, we generated a final genotype matrix containing 3967 SNPs from 89 DNA samples from a single lane of Illumina HiSeq and used it to create a saturated genetic linkage map and to identify a known QTL underlying apple skin color. We therefore demonstrate that GBS is a cost-effective method for generating genome-wide SNP data suitable for genetic mapping in a highly diverse and heterozygous agricultural species. We anticipate future improvements to the GBS analysis pipeline presented here that will enhance the utility of next-generation DNA sequence data for the purposes of genetic mapping across diverse species.

SSR_pipeline: a bioinformatic infrastructure for identifying microsatellites from paired-end Illumina high-throughput DNA sequencing data

SSR_pipeline is a flexible set of programs designed to efficiently identify simple sequence repeats (e.g., microsatellites) from paired-end high-throughput Illumina DNA sequencing data. The program suite contains 3 analysis modules along with a fourth control module that can automate analyses of large volumes of data. The modules are used to 1) identify the subset of paired-end sequences that pass Illumina quality standards, 2) align paired-end reads into a single composite DNA sequence, and 3) identify sequences that possess microsatellites (both simple and compound) conforming to user-specified parameters. The microsatellite search algorithm is extremely efficient, and we have used it to identify repeats with motifs from 2 to 25 bp in length. Each of the 3 analysis modules can also be used independently to provide greater flexibility or to work with FASTQ or FASTA files generated from other sequencing platforms (Roche 454, Ion Torrent, etc.). We demonstrate use of the program with data from the brine fly Ephydra packardi (Diptera: Ephydridae) and provide empirical timing benchmarks to illustrate program performance on a common desktop computer environment. We further show that the Illumina platform is capable of identifying large numbers of microsatellites, even when using unenriched sample libraries and a very small percentage of the sequencing capacity from a single DNA sequencing run. All modules from SSR_pipeline are implemented in the Python programming language and can therefore be used from nearly any computer operating system (Linux, Macintosh, and Windows).

Next-generation DNA sequencing of HEXA: a step in the right direction for carrier screening

Tay-Sachs disease (TSD) is the prototype for ethnic-based carrier screening, with a carrier rate of ∼1/27 in Ashkenazi Jews and French Canadians. HexA enzyme analysis is the current gold standard for TSD carrier screening (detection rate ∼98%), but has technical limitations. We compared DNA analysis by next-generation DNA sequencing (NGS) plus an assay for the 7.6 kb deletion to enzyme analysis for TSD carrier screening using 74 samples collected from participants at a TSD family conference. Fifty-one of 74 participants had positive enzyme results (46 carriers, five late-onset Tay-Sachs [LOTS]), 16 had negative, and seven had inconclusive results. NGS + 7.6 kb del screening of HEXA found a pathogenic mutation, pseudoallele, or variant of unknown significance (VUS) in 100% of the enzyme-positive or obligate carrier/enzyme-inconclusive samples. NGS detected the B1 allele in two enzyme-negative obligate carriers. Our data indicate that NGS can be used as a TSD clinical carrier screening tool. We demonstrate that NGS can be superior in detecting TSD carriers compared to traditional enzyme and genotyping methodologies, which are limited by false-positive and false-negative results and ethnically focused, limited mutation panels, respectively, but is not ready for sole use due to lack of information regarding some VUS.

Non-invasive prenatal testing using massively parallel sequencing of maternal plasma DNA: from molecular karyotyping to fetal whole-genome sequencing

The discovery of cell-free fetal DNA in maternal plasma in 1997 has stimulated a rapid development of non-invasive prenatal testing. The recent advent of massively parallel sequencing has allowed the analysis of circulating cell-free fetal DNA to be performed with unprecedented sensitivity and precision. Fetal trisomies 21, 18 and 13 are now robustly detectable in maternal plasma and such analyses have been available clinically since 2011. Fetal genome-wide molecular karyotyping and whole-genome sequencing have now been demonstrated in a number of proof-of-concept studies. Genome-wide and targeted sequencing of maternal plasma has been shown to allow the non-invasive prenatal testing of β-thalassaemia and can potentially be generalized to other monogenic diseases. It is thus expected that plasma DNA-based non-invasive prenatal testing will play an increasingly important role in future obstetric care. It is thus timely and important that the ethical, social and legal issues of non-invasive prenatal testing be discussed actively by all parties involved in prenatal care.

High intratumor genetic heterogeneity is related to worse outcome in patients with head and neck squamous cell carcinoma

BACKGROUND

Although the presence of genetic heterogeneity within the tumors of individual patients is established, it is unclear whether greater heterogeneity predicts a worse outcome. A quantitative measure of genetic heterogeneity based on next-generation sequencing (NGS) data, mutant-allele tumor heterogeneity (MATH), was previously developed and applied to a data set on head and neck squamous cell carcinoma (HNSCC). Whether this measure correlates with clinical outcome was not previously assessed.

METHODS

The authors examined the association between MATH and clinical, pathologic, and overall survival data for 74 patients with HNSCC for whom exome sequencing was completed.

RESULTS

High MATH (a MATH value above the median) was found to be significantly associated with shorter overall survival (hazards ratio, 2.5; 95% confidence interval, 1.3-4.8). MATH was similarly found to be associated with adverse outcomes in clinically high-risk patients with an advanced stage of disease, and in those with tumors classified as high risk on the basis of validated biomarkers including those that were negative for human papillomavirus or having disruptive tumor protein p53 mutations. In patients who received chemotherapy, the hazards ratio for high MATH was 4.1 (95% confidence interval, 1.6-10.2).

CONCLUSIONS

This novel measure of tumor genetic heterogeneity is significantly associated with tumor progression and adverse treatment outcomes, thereby supporting the hypothesis that higher genetic heterogeneity portends a worse clinical outcome in patients with HNSCC. The prognostic value of some known biomarkers may be the result of their association with high genetic heterogeneity. MATH provides a useful measure of that heterogeneity to be prospectively validated as NGS data from homogeneously treated patient cohorts become available.

Random amplification and pyrosequencing for identification of novel viral genome sequences

ssRNA viruses have high levels of genomic divergence, which can lead to difficulty in genomic characterization of new viruses using traditional PCR amplification and sequencing methods. In this study, random reverse transcription, anchored random PCR amplification, and high-throughput pyrosequencing were used to identify orthobunyavirus sequences from total RNA extracted from viral cultures of acute febrile illness specimens. Draft genome sequence for the orthobunyavirus L segment was assembled and sequentially extended using de novo assembly contigs from pyrosequencing reads and orthobunyavirus sequences in GenBank as guidance. Accuracy and continuous coverage were achieved by mapping all reads to the L segment draft sequence. Subsequently, RT-PCR and Sanger sequencing were used to complete the genome sequence. The complete L segment was found to be 6936 bases in length, encoding a 2248-aa putative RNA polymerase. The identified L segment was distinct from previously published South American orthobunyaviruses, sharing 63% and 54% identity at the nucleotide and amino acid level, respectively, with the complete Oropouche virus L segment and 73% and 81% identity at the nucleotide and amino acid level, respectively, with a partial Caraparu virus L segment. The result demonstrated the effectiveness of a sequence-independent amplification and next-generation sequencing approach for obtaining complete viral genomes from total nucleic acid extracts and its use in pathogen discovery.

Methods in DNA methylation profiling

Metastable and somatically heritable patterns of DNA methylation provide an important level of genomic regulation. In this article, we review methods for analyzing these genome-wide epigenetic patterns and offer a perspective on the ever-expanding literature, which we hope will be useful for investigators who are new to this area. The historical aspects that we cover will be helpful in interpreting this literature and we hope that our discussion of the newest analytical methods will stimulate future progress. We emphasize that no single approach can provide a complete view of the overall methylome, and that combinations of several modalities applied to the same sample set will give the clearest picture. Given the unexpected epigenomic patterns and new biological principles, as well as new disease markers, that have been uncovered in recent studies, it is likely that important discoveries will continue to be made using genome-wide DNA methylation profiling.

Using Drosophila melanogaster as a Model for Genotoxic Chemical Mutational Studies with a New Program, SnpSift

This paper describes a new program SnpSift for filtering differential DNA sequence variants between two or more experimental genomes after genotoxic chemical exposure. Here, we illustrate how SnpSift can be used to identify candidate phenotype-relevant variants including single nucleotide polymorphisms, multiple nucleotide polymorphisms, insertions, and deletions (InDels) in mutant strains isolated from genome-wide chemical mutagenesis of Drosophila melanogaster. First, the genomes of two independently isolated mutant fly strains that are allelic for a novel recessive male-sterile locus generated by genotoxic chemical exposure were sequenced using the Illumina next-generation DNA sequencer to obtain 20- to 29-fold coverage of the euchromatic sequences. The sequencing reads were processed and variants were called using standard bioinformatic tools. Next, SnpEff was used to annotate all sequence variants and their potential mutational effects on associated genes. Then, SnpSift was used to filter and select differential variants that potentially disrupt a common gene in the two allelic mutant strains. The potential causative DNA lesions were partially validated by capillary sequencing of polymerase chain reaction-amplified DNA in the genetic interval as defined by meiotic mapping and deletions that remove defined regions of the chromosome. Of the five candidate genes located in the genetic interval, the Pka-like gene CG12069 was found to carry a separate pre-mature stop codon mutation in each of the two allelic mutants whereas the other four candidate genes within the interval have wild-type sequences. The Pka-like gene is therefore a strong candidate gene for the male-sterile locus. These results demonstrate that combining SnpEff and SnpSift can expedite the identification of candidate phenotype-causative mutations in chemically mutagenized Drosophila strains. This technique can also be used to characterize the variety of mutations generated by genotoxic chemicals.

Pyrosequencing as a tool for better understanding of human microbiomes

Next-generation sequencing technologies have revolutionized the analysis of microbial communities in diverse environments, including the human body. This article reviews several aspects of one of these technologies, the pyrosequencing technique, including its principles, applications, and significant contribution to the study of the human microbiome, with especial emphasis on the oral microbiome. The results brought about by pyrosequencing studies have significantly contributed to refining and augmenting the knowledge of the community membership and structure in and on the human body in healthy and diseased conditions. Because most oral infectious diseases are currently regarded as biofilm-related polymicrobial infections, high-throughput sequencing technologies have the potential to disclose specific patterns related to health or disease. Further advances in technology hold the perspective to have important implications in terms of accurate diagnosis and more effective preventive and therapeutic measures for common oral diseases.