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.

Measurement of FGFR3 signaling at the cell membrane via total internal reflection fluorescence microscopy to compare the activation of FGFR3 mutants.

Fibroblast growth factor receptors (FGFRs) initiate signal transduction via the RAS/mitogen-activated protein kinase pathway by their tyrosine kinase activation known to determine cell growth, tissue differentiation, and apoptosis. Recently, many missense mutations have been reported for FGFR3, but we only know the functional effect for a handful of them. Some mutations result in aberrant FGFR3 signaling and are associated with various genetic disorders and oncogenic conditions. Here, we employed micropatterned surfaces to specifically enrich fluorophore-tagged FGFR3 (monomeric GFP [mGFP]-FGFR3) in certain areas of the plasma membrane of living cells. We quantified receptor activation via total internal reflection fluorescence microscopy of FGFR3 signaling at the cell membrane that captured the recruitment of the downstream signal transducer growth factor receptor-bound 2 (GRB2) tagged with mScarlet (GRB2-mScarlet) to FGFR3 micropatterns. With this system, we tested the activation of FGFR3 upon ligand addition (fgf1 and fgf2) for WT and four FGFR3 mutants associated with congenital disorders (G380R, Y373C, K650Q, and K650E). Our data showed that ligand addition increased GRB2 recruitment to WT FGFR3, with fgf1 having a stronger effect than fgf2. For all mutants, we found an increased basal receptor activity, and only for two of the four mutants (G380R and K650Q), activity was further increased upon ligand addition. Compared with previous reports, two mutant receptors (K650Q and K650E) had either an unexpectedly high or low activation state, respectively. This can be attributed to the different methodology, since micropatterning specifically captures signaling events at the plasma membrane. Collectively, our results provide further insight into the functional effects of mutations to FGFR3.

High prevalence of somatic PIK3CA and TP53 pathogenic variants in the normal mammary gland tissue of sporadic breast cancer patients revealed by duplex sequencing.

The mammary gland undergoes hormonally stimulated cycles of proliferation, lactation, and involution. We hypothesized that these factors increase the mutational burden in glandular tissue and may explain high cancer incidence rate in the general population, and recurrent disease. Hence, we investigated the DNA sequence variants in the normal mammary gland, tumor, and peripheral blood from 52 reportedly sporadic breast cancer patients. Targeted resequencing of 542 cancer-associated genes revealed subclonal somatic pathogenic variants of: PIK3CA, TP53, AKT1, MAP3K1, CDH1, RB1, NCOR1, MED12, CBFB, TBX3, and TSHR in the normal mammary gland at considerable allelic frequencies (9 × 10-2- 5.2 × 10-1), indicating clonal expansion. Further evaluation of the frequently damaged PIK3CA and TP53 genes by ultra-sensitive duplex sequencing demonstrated a diversified picture of multiple low-level subclonal (in 10-2-10-4 alleles) hotspot pathogenic variants. Our results raise a question about the oncogenic potential in non-tumorous mammary gland tissue of breast-conserving surgery patients.

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.

Processing Hundreds of SARS-CoV-2 Samples with an In-House PCR-Based Method without Robotics.

The SARS-CoV-2 pandemic has required the development of multiple testing systems to monitor and control the viral infection. Here, we developed a PCR test to screen COVID-19 infections that can process up to ~180 samples per day without the requirement of robotics. For this purpose, we implemented the use of multichannel pipettes and plate magnetics for the RNA extraction step and combined the reverse transcription with the qPCR within one step. We tested the performance of two RT-qPCR kits as well as different sampling buffers and showed that samples taken in NaCl or PBS are stable and compatible with different COVID-19 testing systems. Finally, we designed a new internal control based on the human RNase P gene that does not require a DNA digestion step. Our protocol is easy to handle and reaches the sensitivity and accuracy of the standardized diagnostic protocols used in the clinic to detect COVID-19 infections.

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

[This corrects the article DOI: 10.1093/nargab/lqab002.].

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.

Pathogenic postzygotic mosaicism in the tyrosine receptor kinase pathway: potential unidentified human disease hidden away in a few cells.

Mutations occurring during embryonic development affect only a subset of cells resulting in two or more distinct cell populations that are present at different levels, also known as postzygotic mosaicism (PZM). Although PZM is a common biological phenomenon, it is often overlooked as a source of disease due to the challenges associated with its detection and characterization, especially for very low-frequency variants. Moreover, PZM can cause a different phenotype compared to constitutional mutations. Especially, lethal mutations in receptor tyrosine kinase (RTK) pathway genes, which exist only in a mosaic state, can have completely new clinical manifestations and can look very different from the associated monogenic disorder. However, some key questions are still not addressed, such as the level of mosaicism resulting in a pathogenic phenotype and how the clinical outcome changes with the development and age. Addressing these questions is not trivial as we require methods with the sensitivity to capture some of these variants hidden away in very few cells. Recent ultra-accurate deep-sequencing approaches can now identify these low-level mosaics and will be central to understand systemic and local effects of mosaicism in the RTK pathway. The main focus of this review is to highlight the importance of low-level mosaics and the need to include their detection in studies of genomic variation associated with disease.

Family reunion via error correction: an efficient analysis of duplex sequencing data.

BACKGROUND: Duplex sequencing is the most accurate approach for identification of sequence variants present at very low frequencies. Its power comes from pooling together multiple descendants of both strands of original DNA molecules, which allows distinguishing true nucleotide substitutions from PCR amplification and sequencing artifacts. This strategy comes at a cost-sequencing the same molecule multiple times increases dynamic range but significantly diminishes coverage, making whole genome duplex sequencing prohibitively expensive. Furthermore, every duplex experiment produces a substantial proportion of singleton reads that cannot be used in the analysis and are thrown away. RESULTS: In this paper we demonstrate that a significant fraction of these reads contains PCR or sequencing errors within duplex tags. Correction of such errors allows "reuniting" these reads with their respective families increasing the output of the method and making it more cost effective. CONCLUSIONS: We combine an error correction strategy with a number of algorithmic improvements in a new version of the duplex analysis software, Du Novo 2.0. It is written in Python, C, AWK, and Bash. It is open source and readily available through Galaxy, Bioconda, and Github: https://github.com/galaxyproject/dunovo.

PRDM9 forms a trimer by interactions within the zinc finger array.

PRDM9 is a trans-acting factor directing meiotic recombination to specific DNA-binding sites by its zinc finger (ZnF) array. It was suggested that PRDM9 is a multimer; however, we do not know the stoichiometry or the components inducing PRDM9 multimerization. In this work, we used in vitro binding studies and characterized with electrophoretic mobility shift assays, mass spectrometry, and fluorescence correlation spectroscopy the stoichiometry of the PRDM9 multimer of two different murine PRDM9 alleles carrying different tags and domains produced with different expression systems. Based on the migration distance of the PRDM9-DNA complex, we show that PRDM9 forms a trimer. Moreover, this stoichiometry is adapted already by the free, soluble protein with little exchange between protein monomers. The variable ZnF array of PRDM9 is sufficient for multimerization, and at least five ZnFs form already a functional trimer. Finally, we also show that only one ZnF array within the PRDM9 oligomer binds to the DNA, whereas the remaining two ZnF arrays likely maintain the trimer by ZnF-ZnF interactions.

The impact of poly-A microsatellite heterologies in meiotic recombination.

Meiotic recombination has strong, but poorly understood effects on short tandem repeat (STR) instability. Here, we screened thousands of single recombinant products with sperm typing to characterize the role of polymorphic poly-A repeats at a human recombination hotspot in terms of hotspot activity and STR evolution. We show that the length asymmetry between heterozygous poly-A's strongly influences the recombination outcome: a heterology of 10 A's (9A/19A) reduces the number of crossovers and elevates the frequency of non-crossovers, complex recombination products, and long conversion tracts. Moreover, the length of the heterology also influences the STR transmission during meiotic repair with a strong and significant insertion bias for the short heterology (6A/7A) and a deletion bias for the long heterology (9A/19A). In spite of this opposing insertion-/deletion-biased gene conversion, we find that poly-A's are enriched at human recombination hotspots that could have important consequences in hotspot activation.

LDJump: Estimating variable recombination rates from population genetic data.

As recombination plays an important role in evolution, its estimation and the identification of hotspot positions is of considerable interest. We propose a novel approach for estimating population recombination rates based on genotyping or sequence data that involves a sequential multiscale change point estimator. Our method also permits demography to be taken into account. It uses several summary statistics within a regression model fitted on suitable scenarios. Our proposed method is accurate, computationally fast, and provides a parsimonious solution by ensuring a type I error control against too many changes in the recombination rate. An application to human genome data suggests a good congruence between our estimated and experimentally identified hotspots. Our method is implemented in the R-package LDJump, which is freely available at https://github.com/PhHermann/LDJump.

The consequences of sequence erosion in the evolution of recombination hotspots.

Meiosis is initiated by a double-strand break (DSB) introduced in the DNA by a highly controlled process that is repaired by recombination. In many organisms, recombination occurs at specific and narrow regions of the genome, known as recombination hotspots, which overlap with regions enriched for DSBs. In recent years, it has been demonstrated that conversions and mutations resulting from the repair of DSBs lead to a rapid sequence evolution at recombination hotspots eroding target sites for DSBs. We still do not fully understand the effect of this erosion in the recombination activity, but evidence has shown that the binding of trans-acting factors like PRDM9 is affected. PRDM9 is a meiosis-specific, multi-domain protein that recognizes DNA target motifs by its zinc finger domain and directs DSBs to these target sites. Here we discuss the changes in affinity of PRDM9 to eroded recognition sequences, and explain how these changes in affinity of PRDM9 can affect recombination, leading sometimes to sterility in the context of hybrid crosses. We also present experimental data showing that DNA methylation reduces PRDM9 binding in vitro Finally, we discuss PRDM9-independent hotspots, posing the question how these hotspots evolve and change with sequence erosion.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Water transport through the intestinal epithelial barrier under different osmotic conditions is dependent on LI-cadherin trans-interaction.

In the intestine water has to be reabsorbed from the chymus across the intestinal epithelium. The osmolarity within the lumen is subjected to high variations meaning that water transport often has to take place against osmotic gradients. It has been hypothesized that LI-cadherin is important in this process by keeping the intercellular cleft narrow facilitating the buildup of an osmotic gradient allowing water reabsorption. LI-cadherin is exceptional among the cadherin superfamily with respect to its localization along the lateral plasma membrane of epithelial cells being excluded from adherens junction. Furthermore it has 7 but not 5 extracellular cadherin repeats (EC1-EC7) and a small cytosolic domain. In this study we identified the peptide VAALD as an inhibitor of LI-cadherin trans-interaction by modeling the structure of LI-cadherin and comparison with the known adhesive interfaces of E-cadherin. This inhibitory peptide was used to measure LI-cadherin dependency of water transport through a monolayer of epithelial CACO2 cells under various osmotic conditions. If LI-cadherin trans-interaction was inhibited by use of the peptide, water transport from the luminal to the basolateral side was impaired and even reversed in the case of hypertonic conditions whereas no effect could be observed at isotonic conditions. These data are in line with a recently published model predicting LI-cadherin to keep the width of the lateral intercellular cleft small. In this narrow cleft a high osmolarity can be achieved due to ion pumps yielding a standing osmotic gradient allowing water absorption from the gut even if the faeces is highly hypertonic.

The long zinc finger domain of PRDM9 forms a highly stable and long-lived complex with its DNA recognition sequence.

PR domain containing protein 9 (PRDM9) is a meiosis-specific, multi-domain protein that regulates the location of recombination hotspots by targeting its DNA recognition sequence for double-strand breaks (DSBs). PRDM9 specifically recognizes DNA via its tandem array of zinc fingers (ZnFs), epigenetically marks the local chromatin by its histone methyltransferase activity, and is an important tether that brings the DNA into contact with the recombination initiation machinery. A strong correlation between PRDM9-ZnF variants and specific DNA motifs at recombination hotspots has been reported; however, the binding specificity and kinetics of the ZnF domain are still obscure. Using two in vitro methods, gel mobility shift assays and switchSENSE, a quantitative biophysical approach that measures binding rates in real time, we determined that the PRDM9-ZnF domain forms a highly stable and long-lived complex with its recognition sequence, with a dissociation halftime of many hours. The ZnF domain exhibits an equilibrium dissociation constant (K D) in the nanomolar (nM) range, with polymorphisms in the recognition sequence directly affecting the binding affinity. We also determined that alternative sequences (15-16 nucleotides in length) can be specifically bound by different subsets of the ZnF domain, explaining the binding plasticity of PRDM9 for different sequences. Finally, longer binding targets are preferred than predicted from the numbers of ZnFs contacting the DNA. Functionally, a long-lived complex translates into an enzymatically active PRDM9 at specific DNA-binding sites throughout meiotic prophase I that might be relevant in stabilizing the components of the recombination machinery to a specific DNA target until DSBs are initiated by Spo11.

Haplotyping of Heterozygous SNPs in Genomic DNA Using Long-Range PCR.

To study meiotic recombination products, cis- or trans-association of disease polymorphisms, or allele-specific expression patterns, it is necessary to phase heterozygous polymorphisms separated by several kilobases. Haplotyping using long-range polymerase chain reaction (PCR) is a powerful, cost-effective method to directly obtain the phase of multiple heterozygous sites with standard laboratory equipment in a handful of loci for many samples. The method is based on the amplification of large genomic DNA regions (up to ~40 kb) with a reaction mixture that combines a proofreading polymerase with allele-specific primer pairs that preferentially amplify matched templates. The analysis of two heterozygous SNPs requires four reactions, each containing one of the four possible allele-specific primer combinations (two forward and two reverse primers), with the mismatches occurring at the 3' ends of the primers. The two correct primer combinations will more efficiently elongate the matching alleles than the alternative alleles, and the difference in amplification efficiency can be monitored with real-time PCR.

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.

High-Throughput Genotyping with TaqMan Allelic Discrimination and Allele-Specific Genotyping Assays.

Real-time PCR-based genotyping methods, such as TaqMan allelic discrimination assays and allele-specific genotyping, are particularly useful when screening a handful of single nucleotide polymorphisms in hundreds of samples; either derived from different individuals, tissues, or pre-amplified DNA. Although real-time PCR-based methods such as TaqMan are well-established, alternative methods, like allele-specific genotyping, are powerful alternatives, especially for genotyping short tandem repeat (STR) length polymorphisms. Here, we describe all relevant aspects when developing an assay for a new SNP or STR using either TaqMan or allele-specific genotyping, respectively, such as primer and probe design, optimization of reaction conditions, the experimental procedure for typing hundreds of samples, and finally the data evaluation. Our goal is to provide a guideline for developing genotyping assays using these two approaches that render reliable and reproducible genotype calls involving minimal optimization.

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-5) 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.