PolyG-DS: An ultrasensitive polyguanine tract-profiling method to detect clonal expansions and trace cell lineage.

Polyguanine tracts (PolyGs) are short guanine homopolymer repeats that are prone to accumulating mutations when cells divide. This feature makes them especially suitable for cell lineage tracing, which has been exploited to detect and characterize precancerous and cancerous somatic evolution. PolyG genotyping, however, is challenging because of the inherent biochemical difficulties in amplifying and sequencing repetitive regions. To overcome this limitation, we developed PolyG-DS, a next-generation sequencing (NGS) method that combines the error-correction capabilities of duplex sequencing (DS) with enrichment of PolyG loci using CRISPR-Cas9-targeted genomic fragmentation. PolyG-DS markedly reduces technical artifacts by comparing the sequences derived from the complementary strands of each original DNA molecule. We demonstrate that PolyG-DS genotyping is accurate, reproducible, and highly sensitive, enabling the detection of low-frequency alleles (<0.01) in spike-in samples using a panel of only 19 PolyG markers. PolyG-DS replicated prior results based on PolyG fragment length analysis by capillary electrophoresis, and exhibited higher sensitivity for identifying clonal expansions in the nondysplastic colon of patients with ulcerative colitis. We illustrate the utility of this method for resolving the phylogenetic relationship among precancerous lesions in ulcerative colitis and for tracing the metastatic dissemination of ovarian cancer. PolyG-DS enables the study of tumor evolution without prior knowledge of tumor driver mutations and provides a tool to perform cost-effective and easily scalable ultra-accurate NGS-based PolyG genotyping for multiple applications in biology, genetics, and cancer research.

poly(UG)-tailed RNAs in genome protection and epigenetic inheritance.

Mobile genetic elements threaten genome integrity in all organisms. RDE-3 (also known as MUT-2) is a ribonucleotidyltransferase that is required for transposon silencing and RNA interference in Caenorhabditis elegans1-4. When tethered to RNAs in heterologous expression systems, RDE-3 can add long stretches of alternating non-templated uridine (U) and guanosine (G) ribonucleotides to the 3' termini of these RNAs (designated poly(UG) or pUG tails)5. Here we show that, in its natural context in C. elegans, RDE-3 adds pUG tails to targets of RNA interference, as well as to transposon RNAs. RNA fragments attached to pUG tails with more than 16 perfectly alternating 3' U and G nucleotides become gene-silencing agents. pUG tails promote gene silencing by recruiting RNA-dependent RNA polymerases, which use pUG-tailed RNAs (pUG RNAs) as templates to synthesize small interfering RNAs (siRNAs). Our results show that cycles of pUG RNA-templated siRNA synthesis and siRNA-directed pUG RNA biogenesis underlie double-stranded-RNA-directed transgenerational epigenetic inheritance in the C. elegans germline. We speculate that this pUG RNA-siRNA silencing loop enables parents to inoculate progeny against the expression of unwanted or parasitic genetic elements.

Structure Validation of G-Rich RNAs in Noncoding Regions of the Human Genome.

We present the rapid biophysical characterization of six previously reported putative G-quadruplex-forming RNAs from the 5'-untranslated region (5'-UTR) of silvestrol-sensitive transcripts for investigation of their secondary structures. By NMR and CD spectroscopic analysis, we found that only a single sequence-[AGG]2 [CGG]2 C-folds into a single well-defined G-quadruplex structure. Sequences with longer poly-G strands form unspecific aggregates, whereas CGG-repeat-containing sequences exhibit a temperature-dependent equilibrium between a hairpin and a G-quadruplex structure. The applied experimental strategy is fast and provides robust readout for G-quadruplex-forming capacities of RNA oligomers.

Poly-GR dipeptide repeat polymers correlate with neurodegeneration and Clinicopathological subtypes in C9ORF72-related brain disease.

Frontotemporal lobar degeneration (FTLD) is heterogeneous in clinical presentation, neuropathological characteristics and genetics. An expanded GGGGCC hexanucleotide repeat in C9ORF72 is the most common genetic cause of both FTLD and motor neuron disease (MND). Dipeptide repeat polymers (DPR) are generated through repeat-associated non-ATG translation, and they aggregate in neuronal inclusions with a distribution distinct from that of TDP-43 pathology. Recent studies from animal and cell culture models suggest that DPR might be toxic, but that toxicity may differ for specific DPR. Arginine containing DPR (poly-GR and poly-PR) have the greatest toxicity and are less frequent than other DPR (poly-GP, poly-GA). A unique feature of arginine-containing DPR is their potential for post-translational modification by methyl-transferases, which produces methylarginine DPR. In this report, we explored the relationship of DPR and methylarginine to markers of neurodegeneration using quantitative digital microscopic methods in 40 patients with C9ORF72 mutations and one of three different clinicopathologic phenotypes, FTLD, FTLD-MND or MND. We find that density and distribution of poly-GR inclusions are different from poly-GA and poly-GP inclusions. We also demonstrate colocalization of poly-GR with asymmetrical dimethylarginine (aDMA) immunoreactivity in regions with neurodegeneration. Differences in aDMA were also noted by clinical phenotype. FTLD-MND had the highest burden of poly-GR pathology compared to FTLD and MND, while FTLD-MND had higher burden of aDMA than FTLD. The results suggest that poly-GR pathology is associated with toxicity and neurodegeneration. It remains to be determined if dimethylarginine modification of poly-GR could contribute to its toxicity.

No effect of AR polyG polymorphism on spinal and bulbar muscular atrophy phenotype.

BACKGROUND AND PURPOSE: Disease severity varies considerably among patients with Spinal and Bulbar Muscular Atrophy (SBMA). Our aim was to investigate the role of androgen receptor (AR) polymorphic repeats in SBMA phenotype. METHODS: We analyzed the length of AR polyQ and polyG tracts in 159 SBMA patients. RESULTS: No relationship between polyG size or polyG/polyQ haplotypes and clinical phenotype was found. An independent negative correlation between polyQ-length and onset of weakness was confirmed (P < 0.001). CONCLUSIONS: The negative results of our study prompt to continue the search for potential disease modifiers in SBMA outside the AR gene.

Binding of Metallated Porphyrin-Imidazophenazine Conjugate to Tetramolecular Quadruplex Formed by Poly(G): a Spectroscopic Investigation.

The binding of telomerase inhibitor ZnTMPyP(3+)-ImPzn, Zn(II) derivative of tricationic porphyrin-imidazophenazine conjugate, to tetramolecular quadruplex structure formed by poly(G) was studied in aqueous solutions at neutral pH and near physiological ionic strength using absorption and polarized fluorescent spectroscopy techniques. Three binding modes were determined from the dependences of the fluorescence intensity and polarization degree for the porphyrin and phenazine moieties of the conjugate on molar polymer-to-dye ratio (P/D). The first one is outside electrostatic binding of positively charged porphyrin fragments to anionic phosphate groups of the polymer which prevails only at very low P/D values and manifests itself by substantial fluorescence quenching. It is suggested that the formation of externally bound porphyrin dimers occurs. The other two binding modes observed at high P/D are embedding of the ZnTMPyP(3+) moiety into the groove of poly(G) quadruplex accompanied by more than 3-fold enhancement of the conjugate emission, and simultaneous intercalation of the phenazine fragment between the guanine bases accompanied by the increase of its fluorescence polarization degree up to 0.25. Thus Zn(II) conjugate seems to be promising ligand for the stabilization of G-quadruplex structures since porphyrin binding to poly(G) is strengthened by additional intercalation of phenazine moiety.

Guanine polynucleotides are self-antigens for human natural autoantibodies and are significantly reduced in the human genome.

In the course of investigating anti-DNA autoantibodies, we examined IgM and IgG antibodies to poly-G and other oligonucleotides in the sera of healthy persons and those diagnosed with systemic lupus erythematosus (SLE), scleroderma (SSc), or pemphigus vulgaris (PV); we used an antigen microarray and informatic analysis. We now report that all of the 135 humans studied, irrespective of health or autoimmune disease, manifested relatively high amounts of IgG antibodies binding to the 20-mer G oligonucleotide (G20); no participants entirely lacked this reactivity. IgG antibodies to homo-nucleotides A20, C20 or T20 were present only in the sera of SLE patients who were positive for antibodies to dsDNA. The prevalence of anti-G20 antibodies led us to survey human, mouse and Drosophila melanogaster (fruit fly) genomes for runs of T20 and G20 or more: runs of T20 appear > 170,000 times compared with only 93 runs of G20 or more in the human genome; of these runs, 40 were close to brain-associated genes. Mouse and fruit fly genomes showed significantly lower T20/G20 ratios than did human genomes. Moreover, sera from both healthy and SLE mice contained relatively little or no anti-G20 antibodies; so natural anti-G20 antibodies appear to be characteristic of humans. These unexpected observations invite investigation of the immune functions of anti-G20 antibodies in human health and disease and of runs of G20 in the human genome.

The frequency of poly(G) tracts in the human genome and their use as a sensor of DNA damage.

Tandem repeats of short DNA sequences are commonly found in human DNA. These simple sequence repeats or microsatellites are highly polymorphic in the human genome. Since the anti-tumour agent cisplatin preferentially forms DNA adducts at runs of consecutive guanine nucleotides (poly(G)), the position and frequency of occurrence of poly(G) sequences in the updated human genome was investigated. There are more runs of consecutive guanines than would be expected by random chance. This especially true for poly(G) sequences longer than approximately n=9. A plot of poly(G) length against log(observed/expected) frequency produced a straight line for n>9. A similar observation was also found for poly(A) DNA sequence repeats. This data implied that the increase in observed/expected frequency is directly related to length of DNA repeat. It was proposed that long runs of consecutive guanine nucleotides could be a sensitive sensor of cellular DNA damage since a number of DNA damaging agents cause lesions at poly(G) sequences.

Identification of functional, short-lived isoform of linker for activation of T cells (LAT).

Linker for activation of T cells (LAT) is a transmembrane adaptor protein playing a key role in the development, activation and maintenance of peripheral homeostasis of T cells. In this study we identified a functional isoform of LAT. It originates from an intron 6 retention event generating an in-frame splice variant of LAT mRNA denoted as LATi6. Comparison of LATi6 expression in peripheral blood leukocytes of human and several other mammalian species revealed that it varied from being virtually absent in the mouse to being predominant in the cow. Analysis of LAT isoform frequency expressed from minigene splicing reporters carrying loss- or gain-of-function point mutations within intronic polyguanine sequences showed that these elements are critical for controlling the intron 6 removal. The protein product of LATi6 isoform (LATi6) ectopically expressed in LAT-deficient JCam 2.5 cell line localized correctly to subcellular compartments and supported T-cell receptor signaling but differed from the canonical LAT protein by displaying a shorter half-life and mediating an increased interleukin-2 secretion upon prolonged CD3/CD28 crosslinking. Altogether, our data suggest that the appearance of LATi6 isoform is an evolutionary innovation that may contribute to a more efficient proofreading control of effector T-cell response.

Mathematical and live meningococcal models for simple sequence repeat dynamics - coherent predictions and observations.

Evolvability by means of simple sequence repeat (SSR) instability is a feature under the constant influence of opposing selective pressures to expand and compress the repeat tract and is mechanistically influenced by factors that affect genetic instability. In addition to direct selection for protein expression and structural integrity, other factors that influence tract length evolution were studied. The genetic instability of SSRs that switch the expression of antibiotic resistance ON and OFF was modelled mathematically and monitored in a panel of live meningococcal strains. The mathematical model showed that the SSR length of a theoretical locus in an evolving population may be shaped by direct selection of expression status (ON or OFF), tract length dependent (α) and tract length independent factors (ß). According to the model an increase in α drives the evolution towards shorter tracts. An increase in ß drives the evolution towards a normal distribution of tract lengths given that an upper and a lower limit are set. Insertion and deletion biases were shown to skew allelic distributions in both directions. The meningococcal SSR model was tested in vivo by monitoring the frequency of spectinomycin resistance OFF→ON switching in a designed locus. The instability of a comprehensive panel of the homopolymeric SSRs, constituted of a range of 5-13 guanine nucleotides, was monitored in wildtype and mismatch repair deficient backgrounds. Both the repeat length itself and mismatch repair deficiency were shown to influence the genetic instability of the homopolymeric tracts. A possible insertion bias was observed in tracts ≤G10. Finally, an inverse correlation between the number of tract-encoded amino acids and growth in the presence of ON-selection illustrated a limitation to SSR expansion in an essential gene associated with the designed model locus and the protein function mediating antibiotic resistance.

Hypermutable DNA chronicles the evolution of human colon cancer.

Intratumor genetic heterogeneity reflects the evolutionary history of a cancer and is thought to influence treatment outcomes. Here we report that a simple PCR-based assay interrogating somatic variation in hypermutable polyguanine (poly-G) repeats can provide a rapid and reliable assessment of mitotic history and clonal architecture in human cancer. We use poly-G repeat genotyping to study the evolution of colon carcinoma. In a cohort of 22 patients, we detect poly-G variants in 91% of tumors. Patient age is positively correlated with somatic mutation frequency, suggesting that some poly-G variants accumulate before the onset of carcinogenesis during normal division in colonic stem cells. Poorly differentiated tumors have fewer mutations than well-differentiated tumors, possibly indicating a shorter mitotic history of the founder cell in these cancers. We generate poly-G mutation profiles of spatially separated samples from primary carcinomas and matched metastases to build well-supported phylogenetic trees that illuminate individual patients' path of metastatic progression. Our results show varying degrees of intratumor heterogeneity among patients. Finally, we show that poly-G mutations can be found in other cancers than colon carcinoma. Our approach can generate reliable maps of intratumor heterogeneity in large numbers of patients with minimal time and cost expenditure.

Frequency and distribution of simple and compound microsatellites in forty-eight Human papillomavirus (HPV) genomes.

Simple sequence repeats (SSRs) are tandem-repeated sequences ubiquitously present but differentially distributed across genomes. Present study is a systematic analysis for incidence, composition and complexity of different microsatellites in 48 representative Human papillomavirus (HPV) genomes. The analysis revealed a total of 1868 SSRs and 120 cSSRs. However, four genomes (HPV-60, HPV-92, HPV-112 and HPV-136) lacked any cSSR content; while HPV-31 accounted for a maximum of 10 cSSRs. An overall increase in cSSR% with higher dMAX was observed. The SSRs and cSSRs were prevalent in coding regions. Poly(A/T) repeats were significantly more abundant than poly(G/C) repeats possibly due to high (A/T) content of the HPV genomes. Further, higher prevalence of di-nucleotide repeats over tri-nucleotide repeats may be attributed to instability of former because of higher slippage rate. An in-depth study of the satellite sequences would provide an insight into the imperfections and evolution of microsatellites.

Clonal expansions and short telomeres are associated with neoplasia in early-onset, but not late-onset, ulcerative colitis.

BACKGROUND: Patients with ulcerative colitis (UC) are at risk of developing colorectal cancer. We have previously reported that cancer progression is associated with the presence of clonal expansions and shorter telomeres in nondysplastic mucosa. We sought to validate these findings in an independent case-control study. METHODS: This study included 33 patients with UC: 14 progressors (patients with high-grade dysplasia or cancer) and 19 nonprogressors. For each patient, a mean of 5 nondysplastic biopsies from proximal, mid, and distal colon were assessed for clonal expansions, as determined by clonal length altering mutations in polyguanine tracts, and telomere length, as measured by quantitative PCR. Both parameters were compared with individual clinicopathological characteristics. RESULTS: Clonal expansions and shorter telomeres were more frequent in nondysplastic biopsies from UC progressors than nonprogressors, but only for patients with early-onset of UC (diagnosis at younger than 50 years of age). Late-onset progressor patients had very few or no clonal expansions and longer telomeres. A few nonprogressors exhibited clonal expansions, which were associated with older age and shorter telomeres. In progressors, clonal expansions were associated with proximity to dysplasia. The mean percentage of clonally expanded mutations distinguished early-onset progressors from nonprogressors with 100% sensitivity and 80% specificity. CONCLUSIONS: Early-onset progressors develop cancer in a field of clonally expanded epithelium with shorter telomeres. The detection of these clones in a few random nondysplastic colon biopsies is a promising cancer biomarker in early-onset UC. Curiously, patients with late-onset UC seem to develop cancer without the involvement of such fields.

Splice site strength-dependent activity and genetic buffering by poly-G runs.

Pre-mRNA splicing is regulated through the combinatorial activity of RNA motifs, including splice sites and splicing regulatory elements. Here we show that the activity of the G-run (polyguanine sequence) class of splicing enhancer elements is approximately 4-fold higher when adjacent to intermediate strength 5' splice sites (ss) than when adjacent to weak 5' ss, and approximately 1.3-fold higher relative to strong 5' ss. We observed this dependence on 5' ss strength in both splicing reporters and in global microarray and mRNA-Seq analyses of splicing changes following RNA interference against heterogeneous nuclear ribonucleoprotein (hnRNP) H, which cross-linked to G-runs adjacent to many regulated exons. An exon's responsiveness to changes in hnRNP H levels therefore depends in a complex way on G-run abundance and 5' ss strength. This pattern of activity enables G-runs and hnRNP H to buffer the effects of 5' ss mutations, augmenting both the frequency of 5' ss polymorphism and the evolution of new splicing patterns. Certain other splicing factors may function similarly.

Poly-G/poly-C tracts in the genomes of Caenorhabditis.

BACKGROUND: In the genome of Caenorhabditis elegans, homopolymeric poly-G/poly-C tracts (G/C tracts) exist at high frequency and are maintained by the activity of the DOG-1 protein. The frequency and distribution of G/C tracts in the genomes of C. elegans and the related nematode, C. briggsae were analyzed to investigate possible biological roles for G/C tracts. RESULTS: In C. elegans, G/C tracts are distributed along every chromosome in a non-random pattern. Most G/C tracts are within introns or are close to genes. Analysis of SAGE data showed that G/C tracts correlate with the levels of regional gene expression in C. elegans. G/C tracts are over-represented and dispersed across all chromosomes in another Caenorhabditis species, C. briggsae. However, the positions and distribution of G/C tracts in C. briggsae differ from those in C. elegans. Furthermore, the C. briggsae dog-1 ortholog CBG19723 can rescue the mutator phenotype of C. elegans dog-1 mutants. CONCLUSION: The abundance and genomic distribution of G/C tracts in C. elegans, the effect of G/C tracts on regional transcription levels, and the lack of positional conservation of G/C tracts in C. briggsae suggest a role for G/C tracts in chromatin structure but not in the transcriptional regulation of specific genes.

Identification and functional analysis of common sequence variants in the DFNA15 gene, Brn-3c.

A rare mutation in Brn-3c (Brn3.1, POU4F3) underlies adult onset hearing loss (DFNA15) and targeted deletion of this gene in mice leads to complete deafness due to loss of sensory hair cells from the cochlea. Therefore the aim of our study was to identify and characterise common functional variation in the Brn-3c gene, which could potentially be a genetic risk for more common forms of adult onset hearing loss. We identified seven sequence variants at the Brn-3c locus. One of these, a novel, common variant at position -3432 was extremely complex consisting of a variable guanine repeat that also exhibited single nucleotide substitutions within the poly-guanine repeat: -3432 poly-G polymorphism. In-vitro studies show that this polymorphism modifies binding affinity for the SP1 transcription factor. Furthermore, reporter constructs of the Brn-3c 5'-flanking region containing different -3432 poly-G alleles show altered transcriptional activity when endogenous SP1 levels are reduced using a dominant negative approach. Results also indicate that this effect is influenced by the length of a novel polymorphic (GT)(n) repeat at position -566 in the Brn-3c 5'-flanking region. In summary, our data show there are common sequence variants in the Brn-3c 5'-flanking region that affect transcriptional regulation in vitro; these variants are candidates for large-scale population based association analysis as they could potentially affect the genetic risk for more common types of adult onset hearing loss.

The A645D mutation in the hinge region of the human androgen receptor (AR) gene modulates AR activity, depending on the context of the polymorphic glutamine and glycine repeats.

BACKGROUND: Sufficient androgen receptor (AR) activity is crucial for normal male sexual differentiation. Here we report on two unrelated 46, XY patients suffering from undervirilization and genital malformations. Both patients had a short polyglycine (polyG) repeat of 10 residues and a relatively long polyglutamine (polyQ) repeat of 28 and 30 residues within the transactivation domain of the AR. In addition, they also harbor a rare A645D substitution. OBJECTIVE: We made a set of AR expression plasmid constructs with varying polyQ and polyG tract sizes in context with or without the A645D substitution and analyzed their in vitro transactivation capacity in transfected CHO cells. RESULTS: We found that a short polyG repeat downmodulated AR activity to approximately 60-65% of the wild-type receptor. This effect was aggravated by A645D in context of a long polyQ repeat to less than 50% activity. In contrast, in the context of a short polyQ and a short polyG repeat, the A645D mutation rescues AR activity to almost wild-type levels, demonstrating a contradictory effect of this mutation, depending on the size of the polymorphic repeats. CONCLUSIONS: A combination of a short polyG repeat with a long polyQ repeat and an A645D substitution might contribute to the development of virilization disorders and explain the observed phenotypes of our patients as a form of androgen insensitivity. The whole recreation of AR sequence variations including individual polymorphic repeat sizes could unravel possible interference of mutations and variations on AR activity by in vitro transfection.

A mutational analysis of Caenorhabditis elegans in space.

The International Caenorhabditis elegans Experiment First Flight (ICE-First) was a project using C. elegans as a model organism to study the biological effects of short duration spaceflight (11 days in the International Space Station). As a member of the ICE-First research team, our group focused on the mutational effects of spaceflight. Several approaches were taken to measure mutational changes that occurred during the spaceflight including measurement of the integrity of poly-G/poly-C tracts, determination of the mutation frequency in the unc-22 gene, analysis of lethal mutations captured by the genetic balancer eT1(III;V), and identification of alterations in telomere length. By comparing the efficiency, sensitivity, and convenience of these methods, we deduced that the eT1 balancer system is well-suited for capturing, maintaining and recovering mutational events that occur over several generations during spaceflight. In the course of this experiment, we have extended the usefulness of the eT1 balancer system by identifying the physical breakpoints of the eT1 translocation and have developed a PCR assay to follow the eT1 chromosomes. C. elegans animals were grown in a defined liquid media during the spaceflight. This is the first analysis of genetic changes in C. elegans grown in the defined media. Although no significant difference in mutation rate was detected between spaceflight and control samples, which is not surprising given the short duration of the spaceflight, we demonstrate here the utility of worms as an integrating biological dosimeter for spaceflight.

Direct molecular haplotyping of the IVS-8 poly(TG) and polyT repeat tracts in the cystic fibrosis gene by melting curve analysis of hybridization probes.

BACKGROUND: Molecular haplotyping is a developing technology with great potential for use in clinical diagnostics. We describe a haplotyping method that uses PCR combined with hybridization probes. METHODS: We designed a LightCycler assay that uses fluorescence resonance energy transfer hybridization probes to haplotype the poly(TG) and polyT (TG-T) tract in the IVS-8 region of the CFTR gene. The reporter probe was designed as a perfect match to the TG12-5T allele. RESULTS: Analysis of 132 samples revealed 9 unique derivative melting temperatures (Tms); the lowest was 42.4 degrees C and the highest was 63.6 degrees C. The lowest Tms were in the TGn-9T group, the intermediate Tms in the TGn-7T group, and the highest Tms in the TGn-5T group. Haplotype frequencies were highest (39%) for TG11-7T and lowest (0.4%) for TG13-5T. CONCLUSIONS: Different combinations of polymorphisms under the reporter hybridization probe had unique and characteristic Tms. This property enables genotyping as well as determination of the phase of multiple variants under the probe, a principle we demonstrated by haplotyping the TG-T repeat tract in the IVS-8 region of the CFTR gene.