Massively parallel haplotyping on microscopic beads for the high-throughput phase analysis of single molecules

Exploring FGFR3 Mutations in the Male Germline: Implications for Clonal Germline Expansions and Paternal Age-Related Dysplasias

Delayed fatherhood results in a higher risk of inheriting a new germline mutation that might result in a congenital disorder in the offspring. In particular, some FGFR3 mutations increase in frequency with age, but there are still a large number of uncharacterized FGFR3 mutations that could be expanding in the male germline with potentially early- or late-onset effects in the offspring. Here, we used digital polymerase chain reaction to assess the frequency and spatial distribution of 10 different FGFR3 missense substitutions in the sexually mature male germline. Our functional assessment of the receptor signaling of the variants with biophysical methods showed that 9 of these variants resulted in a higher activation of the receptor´s downstream signaling, resulting in 2 different expansion behaviors. Variants that form larger subclonal expansions in a dissected postmortem testis also showed a positive correlation of the substitution frequency with the sperm donor's age, and a high and ligand-independent FGFR3 activation. In contrast, variants that measured high FGFR3 signaling and elevated substitution frequencies independent of the donor's age did not result in measurable subclonal expansions in the testis. This suggests that promiscuous signal activation might also result in an accumulation of mutations before the sexual maturation of the male gonad with clones staying relatively constant in size throughout time. Collectively, these results provide novel insights into our understanding of the mutagenesis of driver mutations and their resulting mosaicism in the male germline with important consequences for the transmission and recurrence of associated disorders.

Exploring the Micro-Mosaic Landscape of FGFR3 Mutations in the Ageing Male Germline and Their Potential Implications in Meiotic Differentiation

Advanced paternal age increases the risk of transmitting de novo germline mutations, particularly missense mutations activating the receptor tyrosine kinase (RTK) signalling pathway, as exemplified by the FGFR3 mutation, which is linked to achondroplasia (ACH). This risk is attributed to the expansion of spermatogonial stem cells carrying the mutation, forming sub-clonal clusters in the ageing testis, thereby increasing the frequency of mutant sperm and the number of affected offspring from older fathers. While prior studies proposed a correlation between sub-clonal cluster expansion in the testis and elevated mutant sperm production in older donors, limited data exist on the universality of this phenomenon. Our study addresses this gap by examining the testis-expansion patterns, as well as the increases in mutations in sperm for two FGFR3 variants-c.1138G>A (p.G380R) and c.1948A>G (p.K650E)-which are associated with ACH or thanatophoric dysplasia (TDII), respectively. Unlike the ACH mutation, which showed sub-clonal expansion events in an aged testis and a significant increase in mutant sperm with the donor's age, as also reported in other studies, the TDII mutation showed focal mutation pockets in the testis but exhibited reduced transmission into sperm and no significant age-related increase. The mechanism behind this divergence remains unclear, suggesting potential pleiotropic effects of aberrant RTK signalling in the male germline, possibly hindering differentiation requiring meiosis. This study provides further insights into the transmission risks of micro-mosaics associated with advanced paternal age in the male germline.

Discovery of an unusually high number of de novo mutations in sperm of older men using duplex sequencing

De novo mutations (DNMs) are important players in heritable diseases and evolution. Of particular interest are highly recurrent DNMs associated with congenital disorders that have been described as selfish mutations expanding in the male germline, thus becoming more frequent with age. Here, we have adapted duplex sequencing (DS), an ultradeep sequencing method that renders sequence information on both DNA strands; thus, one mutation can be reliably called in millions of sequenced bases. With DS, we examined ∼4.5 kb of the FGFR3 coding region in sperm DNA from older and younger donors. We identified sites with variant allele frequencies (VAFs) of 10-4 to 10-5, with an overall mutation frequency of the region of ∼6 × 10-7 Some of the substitutions are recurrent and are found at a higher VAF in older donors than in younger ones or are found exclusively in older donors. Also, older donors harbor more mutations associated with congenital disorders. Other mutations are present in both age groups, suggesting that these might result from a different mechanism (e.g., postzygotic mosaicism). We also observe that independent of age, the frequency and deleteriousness of the mutational spectra are more similar to COSMIC than to gnomAD variants. Our approach is an important strategy to identify mutations that could be associated with a gain of function of the receptor tyrosine kinase activity, with unexplored consequences in a society with delayed fatherhood.

Increased yields of duplex sequencing data by a series of quality control tools

Duplex sequencing is currently the most reliable method to identify ultra-low frequency DNA variants by grouping sequence reads derived from the same DNA molecule into families with information on the forward and reverse strand. However, only a small proportion of reads are assembled into duplex consensus sequences (DCS), and reads with potentially valuable information are discarded at different steps of the bioinformatics pipeline, especially reads without a family. We developed a bioinformatics toolset that analyses the tag and family composition with the purpose to understand data loss and implement modifications to maximize the data output for the variant calling. Specifically, our tools show that tags contain polymerase chain reaction and sequencing errors that contribute to data loss and lower DCS yields. Our tools also identified chimeras, which likely reflect barcode collisions. Finally, we also developed a tool that re-examines variant calls from raw reads and provides different summary data that categorizes the confidence level of a variant call by a tier-based system. With this tool, we can include reads without a family and check the reliability of the call, that increases substantially the sequencing depth for variant calling, a particular important advantage for low-input samples or low-coverage regions.

Optimized methodology for product recovery following emulsion PCR: applications for amplification of aptamer libraries and other complex templates

Bias and background issues make efficient amplification of complex template mixes such as aptamer and genomic DNA libraries via conventional PCR methods difficult; emulsion PCR is being increasingly used in such scenarios to circumvent these problems. However, before products generated via emulsion PCR can be used in downstream workflows, they need to be recovered from the water-in-oil emulsion. Often, emulsions are broken following amplification using volatile organic solvents, and product is subsequently isolated via precipitation. Unfortunately, the use of such solvents requires the implementation of special environmental controls, and the yield and purity of DNA isolated by precipitation can be highly variable. Here, we describe the optimization of a simple protocol which can be used to recover products following emulsion PCR using a 2-butanol extraction and subsequent DNA isolation via a commercially available clean-up kit. This protocol avoids the use of volatile solvents and precipitation steps, and we demonstrate that it can be used to reliably recover DNA from water-in-oil emulsions with efficiencies as high as 90%. Furthermore, we illustrate the practical applicability of this protocol by demonstrating how it can be implemented to recover a complex random aptamer library following amplification via emulsion PCR.

Challenging the proposed causes of the PCR plateau phase

Despite the wide-spread use of the polymerase chain reaction (PCR) in various life-science applications, the causes of arrested amplicon generation in late cycles have not been confidently identified. This so-called plateau phase has been attributed to depletion or thermal break-down of primers or nucleotides, thermal inactivation of the DNA polymerase, and product accumulation resulting in competition between primer annealing and product re-hybridization as well as blocking of DNA polymerase by double-stranded amplicons. In the current study, we experimentally investigate the proposed limiting factors of PCR product formation. By applying robust and validated qPCR assays, we elucidate the impact of adding non-target and target amplicons to the reactions, mimicking the high amount of products in late PCR cycles. Further, the impact of increased primer concentrations and thermal stability of reagents are explored. Our results show that high amounts of non-target amplicons inhibit amplification by binding to the DNA polymerase, but that this effect is counteracted by addition of more DNA polymerase or prolonged annealing/extension times. Adding high amounts of target amplicons that also act as templates in the reaction is far less inhibitory to amplification, although a decrease in amplification rate is seen. When primer concentrations are increased, both amplification rates and end-product yields are elevated. Taken together, our results suggest that the main cause of PCR plateau formation is primer depletion and not product accumulation or degradation of reagents. We stress that a PCR plateau caused by primer depletion is assay-dependent, i.e. dependent on the primer design and primer characteristics such as the probability of primer-dimer formation. Our findings contribute to an improved understanding of the major parameters controlling the PCR dynamics at later cycles and the limitations of continued product formation, which in the end can facilitate PCR optimization.

Discovery of Rare Haplotypes by Typing Millions of Single-Molecules with Bead Emulsion Haplotyping (BEH)

Characterizing polymorphisms on single molecules renders the phase of different alleles, and thus, haplotype information. Here, we describe a high-throughput method to genotype hundreds-of thousands single molecules in parallel using bead-emulsion haplotyping (BEH). Haplotyping via BEH is an emulsion-PCR-based method that was adapted to amplify multiple DNA fragments on paramagnetic, microscopic beads within a compartment formed by an aqueous-oil emulsion. This generates beads covered by thousands of clonal copies from several polymorphic regions of an initial DNA molecule that are then genotyped with fluorescently labeled probes. With BEH, up to three different polymorphisms (or more if several polymorphisms are within an amplicon) can be typed within a fragment of several kilobases in a singleexperiment, rendering haplotype information of a very large number of initial single molecules. The high throughput and digital nature of the method makes it ideal to quantify rare haplotypes or to assess the haplotype diversity in complex samples.

Artifactual mutations resulting from DNA lesions limit detection levels in ultrasensitive sequencing applications

The need in cancer research or evolutionary biology to detect rare mutations or variants present at very low frequencies (<10⁻⁵) poses an increasing demand on lowering the detection limits of available methods. Here we demonstrated that amplifiable DNA lesions introduce important error sources in ultrasensitive technologies such as single molecule PCR (smPCR) applications (e.g. droplet-digital PCR), or next-generation sequencing (NGS) based methods. Using templates with known amplifiable lesions (8-oxoguanine, deaminated 5-methylcytosine, uracil, and DNA heteroduplexes), we assessed with smPCR and duplex sequencing that templates with these lesions were amplified very efficiently by proofreading polymerases (except uracil), leading to G->T, and to a lesser extent, to unreported G->C substitutions at 8-oxoguanine lesions, and C->T transitions in amplified uracil containing templates. Long heat incubations common in many DNA extraction protocols significantly increased the number of G->T substitutions. Moreover, in ∼50-80% smPCR reactions we observed the random amplification preference of only one of both DNA strands explaining the known ‘PCR jackpot effect’, with the result that a lesion became indistinguishable from a true mutation or variant. Finally, we showed that artifactual mutations derived from uracil and 8-oxoguanine could be significantly reduced by DNA repair enzymes.

Pitfalls of haplotype phasing from amplicon-based long-read sequencing

The long-read sequencers from Pacific Bioscience (PacBio) and Oxford Nanopore Technologies (ONT) offer the opportunity to phase mutations multiple kilobases apart directly from sequencing reads. In this study, we used long-range PCR with ONT and PacBio sequencing to phase two variants 9 kb apart in the RET gene. We also re-analysed data from a recent paper which had apparently successfully used ONT to phase clinically important haplotypes at the CYP2D6 and HLA loci. From these analyses, we demonstrate PCR-chimera formation during PCR amplification and reference alignment bias are pitfalls that need to be considered when attempting to phase variants using amplicon-based long-read sequencing technologies. These methodological pitfalls need to be avoided if the opportunities provided by long-read sequencers are to be fully exploited.

Distinguishing low frequency mutations from RT-PCR and sequence errors in viral deep sequencing data

BACKGROUND

RNA viruses have high mutation rates and exist within their hosts as large, complex and heterogeneous populations, comprising a spectrum of related but non-identical genome sequences. Next generation sequencing is revolutionising the study of viral populations by enabling the ultra deep sequencing of their genomes, and the subsequent identification of the full spectrum of variants within the population. Identification of low frequency variants is important for our understanding of mutational dynamics, disease progression, immune pressure, and for the detection of drug resistant or pathogenic mutations. However, the current challenge is to accurately model the errors in the sequence data and distinguish real viral variants, particularly those that exist at low frequency, from errors introduced during sequencing and sample processing, which can both be substantial.

RESULTS

We have created a novel set of laboratory control samples that are derived from a plasmid containing a full-length viral genome with extremely limited diversity in the starting population. One sample was sequenced without PCR amplification whilst the other samples were subjected to increasing amounts of RT and PCR amplification prior to ultra-deep sequencing. This enabled the level of error introduced by the RT and PCR processes to be assessed and minimum frequency thresholds to be set for true viral variant identification. We developed a genome-scale computational model of the sample processing and NGS calling process to gain a detailed understanding of the errors at each step, which predicted that RT and PCR errors are more likely to occur at some genomic sites than others. The model can also be used to investigate whether the number of observed mutations at a given site of interest is greater than would be expected from processing errors alone in any NGS data set. After providing basic sample processing information and the site's coverage and quality scores, the model utilises the fitted RT-PCR error distributions to simulate the number of mutations that would be observed from processing errors alone.

CONCLUSIONS

These data sets and models provide an effective means of separating true viral mutations from those erroneously introduced during sample processing and sequencing.

Aquaporin 5 expression in mouse mammary gland cells is not driven by promoter methylation

Several studies have revealed that aquaporins play a role in tumor progression and invasion. In breast carcinomas, high levels of aquaporin 5 (AQP5), a membrane protein involved in water transport, have been linked to increased cell proliferation and migration, thus facilitating tumor progression. Despite the potential role of AQP5 in mammary oncogenesis, the mechanisms controlling mammary AQP5 expression are poorly understood. In other tissues, AQP5 expression has been correlated with its promoter methylation, yet, very little is known about AQP5 promoter methylation in the mammary gland. In this work, we used the mouse mammary gland cell line EpH4, in which we controlled AQP5 expression via the steroid hormone dexamethasone (Dex) to further investigate mechanisms regulating AQP5 expression. In this system, we observed a rapid drop of AQP5 mRNA levels with a delay of several hours in AQP5 protein, suggesting transcriptional control of AQP5 levels. Yet, AQP5 expression was independent of its promoter methylation, or to the presence of negative glucocorticoid receptor elements (nGREs) in its imminent promoter region, but was rather influenced by the cell proliferative state or cell density. We conclude that AQP5 promoter methylation is not a universal mechanism for AQP5 regulation and varies on cell and tissue type.

Crossovers are associated with mutation and biased gene conversion at recombination hotspots

Meiosis is a potentially important source of germline mutations, as sites of meiotic recombination experience recurrent double-strand breaks (DSBs). However, evidence for a local mutagenic effect of recombination from population sequence data has been equivocal, likely because mutation is only one of several forces shaping sequence variation. By sequencing large numbers of single crossover molecules obtained from human sperm for two recombination hotspots, we find direct evidence that recombination is mutagenic: Crossovers carry more de novo mutations than nonrecombinant DNA molecules analyzed for the same donors and hotspots. The observed mutations were primarily CG to TA transitions, with a higher frequency of transitions at CpG than non-CpGs sites. This enrichment of mutations at CpG sites at hotspots could predominate in methylated regions involving frequent single-stranded DNA processing as part of DSB repair. In addition, our data set provides evidence that GC alleles are preferentially transmitted during crossing over, opposing mutation, and shows that GC-biased gene conversion (gBGC) predominates over mutation in the sequence evolution of hotspots. These findings are consistent with the idea that gBGC could be an adaptation to counteract the mutational load of recombination.

Establishment of a reborn MMV-microarray technology: realization of microbiome analysis and other hitherto inaccessible technologies

BACKGROUND

With the accelerating development of bioscience, the problem of research cost has become important. We previously devised and developed a novel concept microarray with manageable volumes (MMV) using a soft gel. It demonstrated the great potential of the MMV technology with the examples of 1024-parallel-cell culture and PCR experiments. However, its full potential failed to be expressed, owing to the nature of the material used for the MMV chip.

RESULTS

In the present study, by developing plastic-based MMVs and associated technologies, we introduced novel technologies such as C2D2P (in which the cells in each well are converted from DNA to protein in 1024-parallel), NGS-non-dependent microbiome analysis, and other powerful applications.

CONCLUSIONS

The reborn MMV-microarray technology has proven to be highly efficient and cost-effective (with approximately 100-fold cost reduction) and enables us to realize hitherto unattainable technologies.

Killer-cell Immunoglobulin-like Receptor gene linkage and copy number variation analysis by droplet digital PCR

The Killer-cell Immunoglobulin-like Receptor (KIR) gene complex has considerable biomedical importance. Patterns of polymorphism in the KIR region include variability in the gene content of haplotypes and diverse structural arrangements. Droplet digital PCR (ddPCR) was used to identify different haplotype motifs and to enumerate KIR copy number variants (CNVs). ddPCR detected a variety of KIR haplotype configurations in DNA from well-characterized cell lines. Mendelian segregation of ddPCR-estimated KIR2DL5 CNVs was observed in Gambian families and CNV typing of other KIRs was shown to be accurate when compared to an established quantitative PCR method.

New evidence for positive selection helps explain the paternal age effect observed in achondroplasia

There are certain de novo germline mutations associated with genetic disorders whose mutation rates per generation are orders of magnitude higher than the genome average. Moreover, these mutations occur exclusively in the male germ line and older men have a higher probability of having an affected child than younger ones, known as the paternal age effect (PAE). The classic example of a genetic disorder exhibiting a PAE is achondroplasia, caused predominantly by a single-nucleotide substitution (c.1138G>A) in FGFR3. To elucidate what mechanisms might be driving the high frequency of this mutation in the male germline, we examined the spatial distribution of the c.1138G>A substitution in a testis from an 80-year-old unaffected man. Using a technology based on bead-emulsion amplification, we were able to measure mutation frequencies in 192 individual pieces of the dissected testis with a false-positive rate lower than 2.7 × 10⁻⁶. We observed that most mutations are clustered in a few pieces with 95% of all mutations occurring in 27% of the total testis. Using computational simulations, we rejected the model proposing an elevated mutation rate per cell division at this nucleotide site. Instead, we determined that the observed mutation distribution fits a germline selection model, where mutant spermatogonial stem cells have a proliferative advantage over unmutated cells. Combined with data on several other PAE mutations, our results support the idea that the PAE, associated with a number of Mendelian disorders, may be explained primarily by a selective mechanism.