Alu element in the RNA binding motif protein, X-linked 2 (RBMX2) gene found to be linked to bipolar disorder.

OBJECTIVE: We have used long-read single molecule, real-time (SMRT) sequencing to fully characterize a ~12Mb genomic region on chromosome Xq24-q27, significantly linked to bipolar disorder (BD) in an extended family from a genetic sub-isolate. This family segregates BD in at least four generations with 24 affected individuals. METHODS: We selected 16 family members for targeted sequencing. The selected individuals either carried the disease haplotype, were non-carriers of the disease haplotype, or served as married-in controls. We designed hybrid capture probes enriching for 5-9Kb fragments spanning the entire 12Mb region that were then sequenced to screen for candidate structural variants (SVs) that could explain the increased risk for BD in this extended family. RESULTS: Altogether, 201 variants were detected in the critically linked region. Although most of these represented common variants, three variants emerged that showed near-perfect segregation among all BD type I affected individuals. Two of the SVs were identified in or near genes belonging to the RNA Binding Motif Protein, X-Linked (RBMX) gene family-a 330bp Alu (subfamily AluYa5) deletion in intron 3 of the RBMX2 gene and an intergenic 27bp tandem repeat deletion between the RBMX and G protein-coupled receptor 101 (GPR101) genes. The third SV was a 50bp tandem repeat insertion in intron 1 of the Coagulation Factor IX (F9) gene. CONCLUSIONS: Among the three genetically linked SVs, additional evidence supported the Alu element deletion in RBMX2 as the leading candidate for contributing directly to the disease development of BD type I in this extended family.

Cross-oncopanel study reveals high sensitivity and accuracy with overall analytical performance depending on genomic regions.

BACKGROUND: Targeted sequencing using oncopanels requires comprehensive assessments of accuracy and detection sensitivity to ensure analytical validity. By employing reference materials characterized by the U.S. Food and Drug Administration-led SEquence Quality Control project phase2 (SEQC2) effort, we perform a cross-platform multi-lab evaluation of eight Pan-Cancer panels to assess best practices for oncopanel sequencing. RESULTS: All panels demonstrate high sensitivity across targeted high-confidence coding regions and variant types for the variants previously verified to have variant allele frequency (VAF) in the 5-20% range. Sensitivity is reduced by utilizing VAF thresholds due to inherent variability in VAF measurements. Enforcing a VAF threshold for reporting has a positive impact on reducing false positive calls. Importantly, the false positive rate is found to be significantly higher outside the high-confidence coding regions, resulting in lower reproducibility. Thus, region restriction and VAF thresholds lead to low relative technical variability in estimating promising biomarkers and tumor mutational burden. CONCLUSION: This comprehensive study provides actionable guidelines for oncopanel sequencing and clear evidence that supports a simplified approach to assess the analytical performance of oncopanels. It will facilitate the rapid implementation, validation, and quality control of oncopanels in clinical use.

A verified genomic reference sample for assessing performance of cancer panels detecting small variants of low allele frequency.

BACKGROUND: Oncopanel genomic testing, which identifies important somatic variants, is increasingly common in medical practice and especially in clinical trials. Currently, there is a paucity of reliable genomic reference samples having a suitably large number of pre-identified variants for properly assessing oncopanel assay analytical quality and performance. The FDA-led Sequencing and Quality Control Phase 2 (SEQC2) consortium analyze ten diverse cancer cell lines individually and their pool, termed Sample A, to develop a reference sample with suitably large numbers of coding positions with known (variant) positives and negatives for properly evaluating oncopanel analytical performance. RESULTS: In reference Sample A, we identify more than 40,000 variants down to 1% allele frequency with more than 25,000 variants having less than 20% allele frequency with 1653 variants in COSMIC-related genes. This is 5-100× more than existing commercially available samples. We also identify an unprecedented number of negative positions in coding regions, allowing statistical rigor in assessing limit-of-detection, sensitivity, and precision. Over 300 loci are randomly selected and independently verified via droplet digital PCR with 100% concordance. Agilent normal reference Sample B can be admixed with Sample A to create new samples with a similar number of known variants at much lower allele frequency than what exists in Sample A natively, including known variants having allele frequency of 0.02%, a range suitable for assessing liquid biopsy panels. CONCLUSION: These new reference samples and their admixtures provide superior capability for performing oncopanel quality control, analytical accuracy, and validation for small to large oncopanels and liquid biopsy assays.

Genomic targets for high-resolution inference of kinship, ancestry and disease susceptibility in orang-utans (genus: Pongo).

BACKGROUND: Orang-utans comprise three critically endangered species endemic to the islands of Borneo and Sumatra. Though whole-genome sequencing has recently accelerated our understanding of their evolutionary history, the costs of implementing routine genome screening and diagnostics remain prohibitive. Capitalizing on a tri-fold locus discovery approach, combining data from published whole-genome sequences, novel whole-exome sequencing, and microarray-derived genotype data, we aimed to develop a highly informative gene-focused panel of targets that can be used to address a broad range of research questions. RESULTS: We identified and present genomic co-ordinates for 175,186 SNPs and 2315 Y-chromosomal targets, plus 185 genes either known or presumed to be pathogenic in cardiovascular (N = 109) or respiratory (N = 43) diseases in humans - the primary and secondary causes of captive orang-utan mortality - or a majority of other human diseases (N = 33). As proof of concept, we designed and synthesized 'SeqCap' hybrid capture probes for these targets, demonstrating cost-effective target enrichment and reduced-representation sequencing. CONCLUSIONS: Our targets are of broad utility in studies of orang-utan ancestry, admixture and disease susceptibility and aetiology, and thus are of value in addressing questions key to the survival of these species. To facilitate comparative analyses, these targets could now be standardized for future orang-utan population genomic studies. The targets are broadly compatible with commercial target enrichment platforms and can be utilized as published here to synthesize applicable probes.

Representative Sequencing: Unbiased Sampling of Solid Tumor Tissue.

Although thousands of solid tumors have been sequenced to date, a fundamental under-sampling bias is inherent in current methodologies. This is caused by a tissue sample input of fixed dimensions (e.g., 6 mm biopsy), which becomes grossly under-powered as tumor volume scales. Here, we demonstrate representative sequencing (Rep-Seq) as a new method to achieve unbiased tumor tissue sampling. Rep-Seq uses fixed residual tumor material, which is homogenized and subjected to next-generation sequencing. Analysis of intratumor tumor mutation burden (TMB) variability shows a high level of misclassification using current single-biopsy methods, with 20% of lung and 52% of bladder tumors having at least one biopsy with high TMB but low clonal TMB overall. Misclassification rates by contrast are reduced to 2% (lung) and 4% (bladder) when a more representative sampling methodology is used. Rep-Seq offers an improved sampling protocol for tumor profiling, with significant potential for improved clinical utility and more accurate deconvolution of clonal structure.

Structural variation and its potential impact on genome instability: Novel discoveries in the EGFR landscape by long-read sequencing.

Structural variation (SV) is typically defined as variation within the human genome that exceeds 50 base pairs (bp). SV may be copy number neutral or it may involve duplications, deletions, and complex rearrangements. Recent studies have shown SV to be associated with many human diseases. However, studies of SV have been challenging due to technological constraints. With the advent of third generation (long-read) sequencing technology, exploration of longer stretches of DNA not easily examined previously has been made possible. In the present study, we utilized third generation (long-read) sequencing techniques to examine SV in the EGFR landscape of four haplotypes derived from two human samples. We analyzed the EGFR gene and its landscape (+/- 500,000 base pairs) using this approach and were able to identify a region of non-coding DNA with over 90% similarity to the most common activating EGFR mutation in non-small cell lung cancer. Based on previously published Alu-element genome instability algorithms, we propose a molecular mechanism to explain how this non-coding region of DNA may be interacting with and impacting the stability of the EGFR gene and potentially generating this cancer-driver gene. By these techniques, we were also able to identify previously hidden structural variation in the four haplotypes and in the human reference genome (hg38). We applied previously published algorithms to compare the relative stabilities of these five different EGFR gene landscape haplotypes to estimate their relative potentials to generate the EGFR exon 19, 15 bp canonical deletion. To our knowledge, the present study is the first to use the differences in genomic architecture between targeted cancer-linked phased haplotypes to estimate their relative potentials to form a common cancer-linked driver mutation.

Tracking Cancer Evolution Reveals Constrained Routes to Metastases: TRACERx Renal.

Clear-cell renal cell carcinoma (ccRCC) exhibits a broad range of metastatic phenotypes that have not been systematically studied to date. Here, we analyzed 575 primary and 335 metastatic biopsies across 100 patients with metastatic ccRCC, including two cases sampledat post-mortem. Metastatic competence was afforded by chromosome complexity, and we identify 9p loss as a highly selected event driving metastasis and ccRCC-related mortality (p = 0.0014). Distinct patterns of metastatic dissemination were observed, including rapid progression to multiple tissue sites seeded by primary tumors of monoclonal structure. By contrast, we observed attenuated progression in cases characterized by high primary tumor heterogeneity, with metastatic competence acquired gradually and initial progression to solitary metastasis. Finally, we observed early divergence of primitive ancestral clones and protracted latency of up to two decades as a feature of pancreatic metastases.

Targeted Bisulfite Sequencing Using the SeqCap Epi Enrichment System.

Cytosine methylation has been shown to have a role in a host of biological processes. In mammalian biology these include stem cell differentiation, embryonic development, genomic imprinting, inflammation, and silencing of transposable elements. Given the central importance of these processes, it is not surprising to find aberrant cytosine methylation patterns associated with many disorders in humans, including cancer, cardiovascular disease, and neurological disease. While whole genome shotgun bisulfite sequencing (WGBS) has recently become feasible, generating high sequence coverage data for the entire genome is expensive, both in terms of money and analysis time, when generally only a small subset of the genome is of interest to most researchers. This report details a procedure for the targeted enrichment of bisulfite treated DNA via SeqCap Epi, allowing high resolution focus of next generation sequencing onto a subset of the genome for high resolution cytosine methylation analysis. Regions ranging in size from only a few kb up to over 200 Mb may be targeted, including the use of the SeqCap Epi CpGiant design which is designed to target 5.5 million CpGs in the human genome. Finally, multiple samples may be multiplexed and sequenced together to provide an inexpensive method of generating methylation data for a large number of samples in a high throughput fashion.

Targeted LncRNA Sequencing with the SeqCap RNA Enrichment System.

Sequencing-based whole-transcriptome analysis (i.e., RNA-Seq) can be a powerful tool when used to measure gene expression, detect novel transcripts, characterize transcript isoforms, and identify sequence polymorphisms. However, this method can be inefficient when the goal is to study only one component of the transcriptome, such as long noncoding RNAs (lncRNAs), which constitute only a small fraction of transcripts in a total RNA sample. Here, we describe a target enrichment method where a total RNA sample is converted to a sequencing-ready cDNA library and hybridized to a complex pool of lncRNA-specific biotinylated long oligonucleotide capture probes prior to sequencing. The resulting sequence data are highly enriched for the targets of interest, dramatically increasing the efficiency of next-generation sequencing approaches for the analysis of lncRNAs.

Dynamic Reorganization of Extremely Long-Range Promoter-Promoter Interactions between Two States of Pluripotency.

Serum-to-2i interconversion of mouse embryonic stem cells (mESCs) is a valuable in vitro model for early embryonic development. To assess whether 3D chromatin organization changes during this transition, we established Capture Hi-C with target-sequence enrichment of DNase I hypersensitive sites. We detected extremely long-range intra- and inter-chromosomal interactions between a small subset of H3K27me3 marked bivalent promoters involving the Hox clusters in serum-grown cells. Notably, these promoter-mediated interactions are not present in 2i ground-state pluripotent mESCs but appear upon their further development into primed-like serum mESCs. Reverting serum mESCs to ground-state 2i mESCs removes these promoter-promoter interactions in a spatiotemporal manner. H3K27me3, which is largely absent at bivalent promoters in ground-state 2i mESCs, is necessary, but not sufficient, to establish these interactions, as confirmed by Capture Hi-C on Eed(-/-) serum mESCs. Our results implicate H3K27me3 and PRC2 as critical players in chromatin alteration during priming of ESCs for differentiation.

LINE-1 Mediated Insertion into Poc1a (Protein of Centriole 1 A) Causes Growth Insufficiency and Male Infertility in Mice.

Skeletal dysplasias are a common, genetically heterogeneous cause of short stature that can result from disruptions in many cellular processes. We report the identification of the lesion responsible for skeletal dysplasia and male infertility in the spontaneous, recessive mouse mutant chagun. We determined that Poc1a, encoding protein of the centriole 1a, is disrupted by the insertion of a processed Cenpw cDNA, which is flanked by target site duplications, suggestive of a LINE-1 retrotransposon-mediated event. Mutant fibroblasts have impaired cilia formation and multipolar spindles. Male infertility is caused by defective spermatogenesis early in meiosis and progressive germ cell loss. Spermatogonial stem cell transplantation studies revealed that Poc1a is essential for normal function of both Sertoli cells and germ cells. The proliferative zone of the growth plate is small and disorganized because chondrocytes fail to re-align after cell division and undergo increased apoptosis. Poc1a and several other genes associated with centrosome function can affect the skeleton and lead to skeletal dysplasias and primordial dwarfisms. This mouse mutant reveals how centrosome dysfunction contributes to defects in skeletal growth and male infertility.

Characterization of functional methylomes by next-generation capture sequencing identifies novel disease-associated variants.

Most genome-wide methylation studies (EWAS) of multifactorial disease traits use targeted arrays or enrichment methodologies preferentially covering CpG-dense regions, to characterize sufficiently large samples. To overcome this limitation, we present here a new customizable, cost-effective approach, methylC-capture sequencing (MCC-Seq), for sequencing functional methylomes, while simultaneously providing genetic variation information. To illustrate MCC-Seq, we use whole-genome bisulfite sequencing on adipose tissue (AT) samples and public databases to design AT-specific panels. We establish its efficiency for high-density interrogation of methylome variability by systematic comparisons with other approaches and demonstrate its applicability by identifying novel methylation variation within enhancers strongly correlated to plasma triglyceride and HDL-cholesterol, including at CD36. Our more comprehensive AT panel assesses tissue methylation and genotypes in parallel at ∼4 and ∼3 M sites, respectively. Our study demonstrates that MCC-Seq provides comparable accuracy to alternative approaches but enables more efficient cataloguing of functional and disease-relevant epigenetic and genetic variants for large-scale EWAS.

Post-conversion targeted capture of modified cytosines in mammalian and plant genomes.

We present a capture-based approach for bisulfite-converted DNA that allows interrogation of pre-defined genomic locations, allowing quantitative and qualitative assessments of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) at CG dinucleotides and in non-CG contexts (CHG, CHH) in mammalian and plant genomes. We show the technique works robustly and reproducibly using as little as 500 ng of starting DNA, with results correlating well with whole genome bisulfite sequencing data, and demonstrate that human DNA can be tested in samples contaminated with microbial DNA. This targeting approach will allow cell type-specific designs to maximize the value of 5mC and 5hmC sequencing.

Genetic perturbation of the maize methylome.

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.

Characterization of the peptidylglycine α-amidating monooxygenase (PAM) from the venom ducts of neogastropods, Conus bullatus and Conus geographus.

Cone snails, genus Conus, are predatory marine snails that use venom to capture their prey. This venom contains a diverse array of peptide toxins, known as conotoxins, which undergo a diverse set of posttranslational modifications. Amidating enzymes modify peptides and proteins containing a C-terminal glycine residue, resulting in loss of the glycine residue and amidation of the preceding residue. A significant fraction of peptides present in the venom of cone snails contain C-terminal amidated residues, which are important for optimizing biological activity. This study describes the characterization of the amidating enzyme, peptidylglycine α-amidating monooxygenase (PAM), present in the venom duct of cone snails, Conus bullatus and Conus geographus. PAM is known to carry out two functions, peptidyl α-hydroxylating monooxygenase (PHM) and peptidylamido-glycolate lyase (PAL). In some animals, such as Drosophila melanogaster, these two functions are present in separate polypeptides, working as individual enzymes. In other animals, such as mammals and in Aplysia californica, PAM activity resides in a single, bifunctional polypeptide. Using specific oligonucleotide primers and reverse transcription-polymerase chain reaction we have identified and cloned from the venom duct cDNA library, a cDNA with 49% homology to PAM from A. californica. We have determined that both the PHM and PAL activities are encoded in one mRNA polynucleotide in both C. bullatus and C. geographus. We have directly demonstrated enzymatic activity catalyzing the conversion of dansyl-YVG-COOH to dansyl-YV-NH2 in cloned cDNA expressed in Drosophila S2 cells.

Spry1 and Spry2 are necessary for lens vesicle separation and corneal differentiation.

PURPOSE: The studies reported here were performed to analyze the roles of Sproutys (Sprys), downstream targets and negative feedback regulators of the fibroblast growth factor (FGF) signaling pathway, in lens and corneal differentiation. METHODS: Spry1 and -2 were conditionally deleted in the lens and corneal epithelial precursors using the Le-Cre transgene and floxed alleles of Spry1 and -2. Alterations in lens and corneal development were assessed by hematoxylin and eosin staining, in situ hybridization, and immunohistochemistry. RESULTS: Spry1 and -2 were upregulated in the lens fibers at the onset of fiber differentiation. FGF signaling was both necessary and sufficient for induction of Spry1 and -2 in the lens fiber cells. Spry1 and -2 single- or double-null lenses failed to separate from the overlying ectoderm and showed persistent keratolenticular stalks. Apoptosis of stalk cells, normally seen during lens vesicle detachment from the ectoderm, was inhibited in Spry mutant lenses, with concomitant ERK activation. Prox1 and p57(KIP2), normally upregulated at the onset of fiber differentiation were prematurely induced in the Spry mutant lens epithelial cells. However, terminal differentiation markers such as ß- or γ-crystallin were not induced. Corneal epithelial precursors in Spry1 and -2 double mutants showed increased proliferation with elevated expression of Erm and DUSP6 and decreased expression of the corneal differentiation marker K12. CONCLUSIONS: Collectively, the results indicate that Spry1 and -2 (1) through negative modulation of ERKs allow lens vesicle separation, (2) are targets of FGF signaling in the lens during initiation of fiber differentiation and (3) function redundantly in the corneal epithelial cells to suppress proliferation.

Targeted enrichment beyond the consensus coding DNA sequence exome reveals exons with higher variant densities.

BACKGROUND: Enrichment of loci by DNA hybridization-capture, followed by high-throughput sequencing, is an important tool in modern genetics. Currently, the most common targets for enrichment are the protein coding exons represented by the consensus coding DNA sequence (CCDS). The CCDS, however, excludes many actual or computationally predicted coding exons present in other databases, such as RefSeq and Vega, and non-coding functional elements such as untranslated and regulatory regions. The number of variants per base pair (variant density) and our ability to interrogate regions outside of the CCDS regions is consequently less well understood. RESULTS: We examine capture sequence data from outside of the CCDS regions and find that extremes of GC content that are present in different subregions of the genome can reduce the local capture sequence coverage to less than 50% relative to the CCDS. This effect is due to biases inherent in both the Illumina and SOLiD sequencing platforms that are exacerbated by the capture process. Interestingly, for two subregion types, microRNA and predicted exons, the capture process yields higher than expected coverage when compared to whole genome sequencing. Lastly, we examine the variation present in non-CCDS regions and find that predicted exons, as well as exonic regions specific to RefSeq and Vega, show much higher variant densities than the CCDS. CONCLUSIONS: We show that regions outside of the CCDS perform less efficiently in capture sequence experiments. Further, we show that the variant density in computationally predicted exons is more than 2.5-times higher than that observed in the CCDS.

Whole exome capture in solution with 3 Gbp of data.

We have developed a solution-based method for targeted DNA capture-sequencing that is directed to the complete human exome. Using this approach allows the discovery of greater than 95% of all expected heterozygous singe base variants, requires as little as 3 Gbp of raw sequence data and constitutes an effective tool for identifying rare coding alleles in large scale genomic studies.

The endocytic recycling regulator EHD1 is essential for spermatogenesis and male fertility in mice.

BACKGROUND: The C-terminal Eps15 homology domain-containing protein 1 (EHD1) is ubiquitously expressed and regulates the endocytic trafficking and recycling of membrane components and several transmembrane receptors. To elucidate the function of EHD1 in mammalian development, we generated Ehd1-/- mice using a Cre/loxP system. RESULTS: Both male and female Ehd1-/- mice survived at sub-Mendelian ratios. A proportion of Ehd1-/- mice were viable and showed smaller size at birth, which continued into adulthood. Ehd1-/- adult males were infertile and displayed decreased testis size, whereas Ehd1-/- females were fertile. In situ hybridization and immunohistochemistry of developing wildtype mouse testes revealed EHD1 expression in most cells of the seminiferous epithelia. Histopathology revealed abnormal spermatogenesis in the seminiferous tubules and the absence of mature spermatozoa in the epididymides of Ehd1-/- males. Seminiferous tubules showed disruption of the normal spermatogenic cycle with abnormal acrosomal development on round spermatids, clumping of acrosomes, misaligned spermatids and the absence of normal elongated spermatids in Ehd1-/- males. Light and electron microscopy analyses indicated that elongated spermatids were abnormally phagocytosed by Sertoli cells in Ehd1-/- mice. CONCLUSIONS: Contrary to a previous report, these results demonstrate an important role for EHD1 in pre- and post-natal development with a specific role in spermatogenesis.