Single-locus enrichment without amplification for sequencing and direct detection of epigenetic modifications.

A gene-level targeted enrichment method for direct detection of epigenetic modifications is described. The approach is demonstrated on the CGG-repeat region of the FMR1 gene, for which large repeat expansions, hitherto refractory to sequencing, are known to cause fragile X syndrome. In addition to achieving a single-locus enrichment of nearly 700,000-fold, the elimination of all amplification steps removes PCR-induced bias in the repeat count and preserves the native epigenetic modifications of the DNA. In conjunction with the single-molecule real-time sequencing approach, this enrichment method enables direct readout of the methylation status and the CGG repeat number of the FMR1 allele(s) for a clonally derived cell line. The current method avoids potential biases introduced through chemical modification and/or amplification methods for indirect detection of CpG methylation events.

Differential increases of specific FMR1 mRNA isoforms in premutation carriers.

BACKGROUND: Over 40% of male and ∼16% of female carriers of a premutation FMR1 allele (55-200 CGG repeats) will develop fragile X-associated tremor/ataxia syndrome, an adult onset neurodegenerative disorder, while about 20% of female carriers will develop fragile X-associated primary ovarian insufficiency. Marked elevation in FMR1 mRNA transcript levels has been observed with premutation alleles, and RNA toxicity due to increased mRNA levels is the leading molecular mechanism proposed for these disorders. However, although the FMR1 gene undergoes alternative splicing, it is unknown whether all or only some of the isoforms are overexpressed in premutation carriers and which isoforms may contribute to the premutation pathology. METHODS: To address this question, we have applied a long-read sequencing approach using single-molecule real-time (SMRT) sequencing and qRT-PCR. RESULTS: Our SMRT sequencing analysis performed on peripheral blood mononuclear cells, fibroblasts and brain tissue samples derived from premutation carriers and controls revealed the existence of 16 isoforms of 24 predicted variants. Although the relative abundance of all mRNA isoforms was significantly increased in the premutation group, as expected based on the bulk increase in mRNA levels, there was a disproportionate (fourfold to sixfold) increase, relative to the overall increase in mRNA, in the abundance of isoforms spliced at both exons 12 and 14, specifically Iso10 and Iso10b, containing the complete exon 15 and differing only in splicing in exon 17. CONCLUSIONS: These findings suggest that RNA toxicity may arise from a relative increase of all FMR1 mRNA isoforms. Interestingly, the Iso10 and Iso10b mRNA isoforms, lacking the C-terminal functional sites for fragile X mental retardation protein function, are the most increased in premutation carriers relative to normal, suggesting a functional relevance in the pathology of FMR1-associated disorders.

Sequencing the unsequenceable: expanded CGG-repeat alleles of the fragile X gene.

The human fragile X mental retardation 1 (FMR1) gene contains a (CGG)(n) trinucleotide repeat in its 5' untranslated region (5'UTR). Expansions of this repeat result in a number of clinical disorders with distinct molecular pathologies, including fragile X syndrome (FXS; full mutation range, greater than 200 CGG repeats) and fragile X-associated tremor/ataxia syndrome (FXTAS; premutation range, 55-200 repeats). Study of these diseases has been limited by an inability to sequence expanded CGG repeats, particularly in the full mutation range, with existing DNA sequencing technologies. Single-molecule, real-time (SMRT) sequencing provides an approach to sequencing that is fundamentally different from other "next-generation" sequencing platforms, and is well suited for long, repetitive DNA sequences. We report the first sequence data for expanded CGG-repeat FMR1 alleles in the full mutation range that reveal the confounding effects of CGG-repeat tracts on both cloning and PCR. A unique feature of SMRT sequencing is its ability to yield real-time information on the rates of nucleoside addition by the tethered DNA polymerase; for the CGG-repeat alleles, we find a strand-specific effect of CGG-repeat DNA on the interpulse distance. This kinetic signature reveals a novel aspect of the repeat element; namely, that the particular G bias within the CGG/CCG-repeat element influences polymerase activity in a manner that extends beyond simple nearest-neighbor effects. These observations provide a baseline for future kinetic studies of repeat elements, as well as for studies of epigenetic and other chemical modifications thereof.

A flexible and efficient template format for circular consensus sequencing and SNP detection.

A novel template design for single-molecule sequencing is introduced, a structure we refer to as a SMRTbell template. This structure consists of a double-stranded portion, containing the insert of interest, and a single-stranded hairpin loop on either end, which provides a site for primer binding. Structurally, this format resembles a linear double-stranded molecule, and yet it is topologically circular. When placed into a single-molecule sequencing reaction, the SMRTbell template format enables a consensus sequence to be obtained from multiple passes on a single molecule. Furthermore, this consensus sequence is obtained from both the sense and antisense strands of the insert region. In this article, we present a universal method for constructing these templates, as well as an application of their use. We demonstrate the generation of high-quality consensus accuracy from single molecules, as well as the use of SMRTbell templates in the identification of rare sequence variants.

Accurate single molecule FRET efficiency determination for surface immobilized DNA using maximum likelihood calculated lifetimes.

Single molecule fluorescent lifetime trajectories of surface immobilized double-stranded DNA coupled with a tetramethylrhodmaine and Cy5 FRET pair were directly measured using time-tagged single-photon counting and scanning confocal microscopy. A modified maximum likelihood estimator (MLE) was developed to compensate for localized background fluorescence and instrument response. With this algorithm, we were able to robustly extract fluorescent lifetimes from their respective decays with as few as 20 photons. Fluorescent lifetimes extracted using an MLE were found to be highly dependent on background fluorescence. We show that appropriate factors are required to extract true lifetime trajectories from single fluorophores.

Long time scale blinking kinetics of cyanine fluorophores conjugated to DNA and its effect on Förster resonance energy transfer.

The blinking kinetics of individual Cy5 fluorophores conjugated to DNA are directly measured using single-molecule spectroscopy. Under deoxygenated aqueous conditions, Cy5 fluorescence exhibits spontaneous and reversible on/off fluctuations with a period lasting seconds. This blinking is observed when directly exciting Cy5 with 640 nm light and by Forster resonance energy transfer (FRET). We find that Cy5 blinking is influenced by the proximity of the donor, the structure of the donor, the presence of 514 nm excitation, and FRET. In the context of single-molecule FRET, blinking of the acceptor produces anticorrelated donor-acceptor intensity fluctuations, which can be difficult to discern from variations in the interdye distance. Slow blinking is, in particular, problematic because it overlaps with biologically relevant time scales. By employing an alternating 514640 nm laser excitation scheme, we show that the dark states can be readily resolved and discriminated from FRET distance fluctuations.

Using fluorescence resonance energy transfer to measure distances along individual DNA molecules: corrections due to nonideal transfer.

Single molecule fluorescence resonance energy transfer has been extensively used to measure distance changes and kinetics in various biomolecular systems. However, due to complications involving multiple de-excitation pathways of the dyes, the absolute inter-dye distance information has seldom been recovered. To circumvent this we directly probe the relative variations in the quantum yield of individual fluorophores. B-DNA was used as a scaffold to position the donor (Cy3 or TMR) at precise distances from the acceptor (Cy5) within the Forster radius. We found that the variation in the Cy3 quantum yield is approximately 5 times larger than that of TMR. By taking into account the molecule-to-molecule variability in the acceptor/donor quantum yield ratio, the apparent fluorescence resonance energy transfer efficiencies were scaled to yield the theoretical values. We obtained very good agreement with a physical model that predicts distances along B-DNA.

Ligand binding to somatostatin receptors induces receptor-specific oligomer formation in live cells.

Heptahelical receptors (HHRs) are generally thought to function as monomeric entities. Several HHRs such as somatostatin receptors (SSTRs), however, form homo- and heterooligomers when activated by ligand binding. By using dual fluorescent ligands simultaneously applied to live cells monotransfected with SSTR5 (R5) or SSTR1 (R1), or cotransfected with R5 and R1, we have analyzed the ligand receptor stoichiometry and aggregation states for the three receptor systems by fluorescence resonance energy transfer and fluorescence correlation spectroscopy. Both homo- and heterooligomeric receptors are occupied by two ligand molecules. We find that monomeric, homooligomeric, and heterooligomeric receptor species occur in the same cell cotransfected with two SSTRs, and that oligomerization of SSTRs is regulated by ligand binding by a selective process that is restricted to some (R5) but not other (R1) SSTR subtypes. We propose that induction by ligand of different oligomeric states of SSTRs represents a unique mechanism for generating signaling specificity not only within the SSTR family but more generally in the HHR family.