GIInger predicts homologous recombination deficiency and patient response to PARPi treatment from shallow genomic profiles.

Homologous recombination deficiency (HRD) is a predictive biomarker for poly(ADP-ribose) polymerase 1 inhibitor (PARPi) sensitivity. Routine HRD testing relies on identifying BRCA mutations, but additional HRD-positive patients can be identified by measuring genomic instability (GI), a consequence of HRD. However, the cost and complexity of available solutions hamper GI testing. We introduce a deep learning framework, GIInger, that identifies GI from HRD-induced scarring observed in low-pass whole-genome sequencing data. GIInger seamlessly integrates into standard BRCA testing workflows and yields reproducible results concordant with a reference method in a multisite study of 327 ovarian cancer samples. Applied to a BRCA wild-type enriched subgroup of 195 PAOLA-1 clinical trial patients, GIInger identified HRD-positive patients who experienced significantly extended progression-free survival when treated with PARPi. GIInger is, therefore, a cost-effective and easy-to-implement method for accurately stratifying patients with ovarian cancer for first-line PARPi treatment.

Benefits of applying molecular barcoding systems are not uniform across different genomic applications.

BACKGROUND: Despite the wide variety of Next Generation Sequencing (NGS)-based methods, it remains challenging to detect mutations present at very low frequencies. This problem is particularly relevant in oncology, where the limiting amount of input material, and its low quality, often limit the performance of the assays. Unique Molecular Identifiers (UMIs) are a molecular barcoding system often coupled with computational methods of noise suppression to improve the reliability of detection of rare variants. Although widely adopted, UMI inclusion imposes additional technical complexity and sequencing cost. Currently, there are no guidelines on UMI usage nor a comprehensive evaluation of their advantage across different applications. METHODS: We used DNA sequencing data generated by molecular barcoding and hybridization-based enrichment, from various types and quantities of input material (fresh frozen, formaldehyde-treated and cell-free DNA), to evaluate the performance of variant calling in different clinically relevant contexts. RESULTS: Noise suppression achieved by read grouping based on fragment mapping positions ensures reliable variant calling for many experimental designs even without exogenous UMIs. Exogenous barcodes significantly improve performance only when mapping position collisions occur, which is common in cell-free DNA. CONCLUSIONS: We demonstrate that UMI usage is not universally beneficial across experimental designs and that it is worthwhile to critically consider the comparative advantage of UMI usage for a given NGS application prior to experimental design.

Improved NGS-based detection of microsatellite instability using tumor-only data.

Microsatellite instability (MSI) is a molecular signature of mismatch repair deficiency (dMMR), a predictive marker of immune checkpoint inhibitor therapy response. Despite its recognized pan-cancer value, most methods only support detection of this signature in colorectal cancer. In addition to the tissue-specific differences that impact the sensitivity of MSI detection in other tissues, the performance of most methods is also affected by patient ethnicity, tumor content, and other sample-specific properties. These limitations are particularly important when only tumor samples are available and restrict the performance and adoption of MSI testing. Here we introduce MSIdetect, a novel solution for NGS-based MSI detection. MSIdetect models the impact of indel burden and tumor content on read coverage at a set of homopolymer regions that we found are minimally impacted by sample-specific factors. We validated MSIdetect in 139 Formalin-Fixed Paraffin-Embedded (FFPE) clinical samples from colorectal and endometrial cancer as well as other more challenging tumor types, such as glioma or sebaceous adenoma or carcinoma. Based on analysis of these samples, MSIdetect displays 100% specificity and 96.3% sensitivity. Limit of detection analysis supports that MSIdetect is sensitive even in samples with relatively low tumor content and limited microsatellite instability. Finally, the results obtained using MSIdetect in tumor-only data correlate well (R=0.988) with what is obtained using tumor-normal matched pairs, demonstrating that the solution addresses the challenges posed by MSI detection from tumor-only data. The accuracy of MSI detection by MSIdetect in different cancer types coupled with the flexibility afforded by NGS-based testing will support the adoption of MSI testing in the clinical setting and increase the number of patients identified that are likely to benefit from immune checkpoint inhibitor therapy.

High-Throughput Nucleotide Resolution Predictions of Assay Limitations Increase the Reliability and Concordance of Clinical Tests.

PURPOSE: The ability of next-generation sequencing (NGS) assays to interrogate thousands of genomic loci has revolutionized genetic testing. However, translation to the clinic is impeded by false-negative results that pose a risk to patients. In response, regulatory bodies are calling for reliability measures to be reported alongside NGS results. Existing methods to estimate reliability do not account for sample- and position-specific variability, which can be significant. Here, we report an approach that computes reliability metrics for every genomic position and sample interrogated by an NGS assay. METHODS: Our approach predicts the limit of detection (LOD), the lowest reliably detectable variant fraction, by taking technical factors into account. We initially explored how LOD is affected by input material amount, library conversion rate, sequencing coverage, and sequencing error rate. This revealed that LOD depends heavily on genomic context and sample properties. Using these insights, we developed a computational approach to predict LOD on the basis of a biophysical model of the NGS workflow. We focused on targeted assays for cell-free DNA, but, in principle, this approach applies to any NGS assay. RESULTS: We validated our approach by showing that it accurately predicts LOD and distinguishes reliable from unreliable results when screening 580 lung cancer samples for actionable mutations. Compared with a standard variant calling workflow, our approach avoided most false negatives and improved interassay concordance from 94% to 99%. CONCLUSION: Our approach, which we name LAVA (LOD-aware variant analysis), reports the LOD for every position and sample interrogated by an NGS assay. This enables reliable results to be identified and improves the transparency and safety of genetic tests.

Recommendations for accurate genotyping of SARS-CoV-2 using amplicon-based sequencing of clinical samples.

OBJECTIVES: Genotyping of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been instrumental in monitoring viral evolution and transmission during the pandemic. The quality of the sequence data obtained from these genotyping efforts depends on several factors, including the quantity/integrity of the input material, the technology, and laboratory-specific implementation. The current lack of guidelines for SARS-CoV-2 genotyping leads to inclusion of error-containing genome sequences in genomic epidemiology studies. We aimed to establish clear and broadly applicable recommendations for reliable virus genotyping. METHODS: We established and used a sequencing data analysis workflow that reliably identifies and removes technical artefacts; such artefacts can result in miscalls when using alternative pipelines to process clinical samples and synthetic viral genomes with an amplicon-based genotyping approach. We evaluated the impact of experimental factors, including viral load and sequencing depth, on correct sequence determination. RESULTS: We found that at least 1000 viral genomes are necessary to confidently detect variants in the SARS-CoV-2 genome at frequencies of ≥10%. The broad applicability of our recommendations was validated in over 200 clinical samples from six independent laboratories. The genotypes we determined for clinical isolates with sufficient quality cluster by sampling location and period. Our analysis also supports the rise in frequencies of 20A.EU1 and 20A.EU2, two recently reported European strains whose dissemination was facilitated by travel during the summer of 2020. CONCLUSIONS: We present much-needed recommendations for the reliable determination of SARS-CoV-2 genome sequences and demonstrate their broad applicability in a large cohort of clinical samples.

The tobacco genome sequence and its comparison with those of tomato and potato.

The allotetraploid plant Nicotiana tabacum (common tobacco) is a major crop species and a model organism, for which only very fragmented genomic sequences are currently available. Here we report high-quality draft genomes for three main tobacco varieties. These genomes show both the low divergence of tobacco from its ancestors and microsynteny with other Solanaceae species. We identify over 90,000 gene models and determine the ancestral origin of tobacco mosaic virus and potyvirus disease resistance in tobacco. We anticipate that the draft genomes will strengthen the use of N. tabacum as a versatile model organism for functional genomics and biotechnology applications.

Protein kinases and phosphatases involved in the acclimation of the photosynthetic apparatus to a changing light environment.

Photosynthetic organisms are subjected to frequent changes in light quality and quantity and need to respond accordingly. These acclimatory processes are mediated to a large extent through thylakoid protein phosphorylation. Recently, two major thylakoid protein kinases have been identified and characterized. The Stt7/STN7 kinase is mainly involved in the phosphorylation of the LHCII antenna proteins and is required for state transitions. It is firmly associated with the cytochrome b(6)f complex, and its activity is regulated by the redox state of the plastoquinone pool. The other kinase, Stl1/STN8, is responsible for the phosphorylation of the PSII core proteins. Using a reverse genetics approach, we have recently identified the chloroplast PPH1/TAP38 and PBPC protein phosphatases, which counteract the activity of STN7 and STN8 kinases, respectively. They belong to the PP2C-type phosphatase family and are conserved in land plants and algae. The picture that emerges from these studies is that of a complex regulatory network of chloroplast protein kinases and phosphatases that is involved in light acclimation, in maintenance of the plastoquinone redox poise under fluctuating light and in the adjustment to metabolic needs.

The phosphorylation status of the chloroplast protein kinase STN7 of Arabidopsis affects its turnover.

The chloroplast serine-threonine protein kinase STN7 of Arabidopsis (Arabidopsis thaliana) is required for the phosphorylation of the light-harvesting system of photosystem II and for state transitions, a process that allows the photosynthetic machinery to balance the light excitation energy between photosystem II and photosystem I and thereby to optimize the photosynthetic yield. Because the STN7 protein kinase of Arabidopsis is known to be phosphorylated at four serine-threonine residues, we have changed these residues by site-directed mutagenesis to alanine (STN7-4A) or aspartic acid (STN7-4D) to assess the role of these phosphorylation events. The corresponding mutants were still able to phosphorylate the light-harvesting system of photosystem II and to perform state transitions. Moreover, we noticed a marked decrease in the level of the STN7 kinase in the wild-type strain under prolonged state 1 conditions that no longer occurs in the STN7-4D mutant. The results suggest a possible role of phosphorylation of the STN7 kinase in regulating its turnover.

Comparative phosphoproteome profiling reveals a function of the STN8 kinase in fine-tuning of cyclic electron flow (CEF).

Important aspects of photosynthetic electron transport efficiency in chloroplasts are controlled by protein phosphorylation. Two thylakoid-associated kinases, STN7 and STN8, have distinct roles in short- and long-term photosynthetic acclimation to changes in light quality and quantity. Although some substrates of STN7 and STN8 are known, the complexity of this regulatory kinase system implies that currently unknown substrates connect photosynthetic performance with the regulation of metabolic and regulatory functions. We performed an unbiased phosphoproteome-wide screen with Arabidopsis WT and stn8 mutant plants to identify unique STN8 targets. The phosphorylation status of STN7 was not affected in stn8, indicating that kinases other than STN8 phosphorylate STN7 under standard growth conditions. Among several putative STN8 substrates, PGRL1-A is of particular importance because of its possible role in the modulation of cyclic electron transfer. The STN8 phosphorylation site on PGRL1-A is absent in both monocotyledonous plants and algae. In dicots, spectroscopic measurements with Arabidopsis WT, stn7, stn8, and stn7/stn8 double-mutant plants indicate a STN8-mediated slowing down of the transition from cyclic to linear electron flow at the onset of illumination. This finding suggests a possible link between protein phosphorylation by STN8 and fine-tuning of cyclic electron flow during this critical step of photosynthesis, when the carbon assimilation is not commensurate to the electron flow capacity of the chloroplast.

Phosphorylation of photosystem II controls functional macroscopic folding of photosynthetic membranes in Arabidopsis.

Photosynthetic thylakoid membranes in plants contain highly folded membrane layers enriched in photosystem II, which uses light energy to oxidize water and produce oxygen. The sunlight also causes quantitative phosphorylation of major photosystem II proteins. Analysis of the Arabidopsis thaliana stn7xstn8 double mutant deficient in thylakoid protein kinases STN7 and STN8 revealed light-independent phosphorylation of PsbH protein and greatly reduced N-terminal phosphorylation of D2 protein. The stn7xstn8 and stn8 mutants deficient in light-induced phosphorylation of photosystem II had increased thylakoid membrane folding compared with wild-type and stn7 plants. Significant enhancement in the size of stacked thylakoid membranes in stn7xstn8 and stn8 accelerated gravity-driven sedimentation of isolated thylakoids and was observed directly in plant leaves by transmission electron microscopy. Increased membrane folding, caused by the loss of light-induced protein phosphorylation, obstructed lateral migration of the photosystem II reaction center protein D1 and of processing protease FtsH between the stacked and unstacked membrane domains, suppressing turnover of damaged D1 in the leaves exposed to high light. These findings show that the high level of photosystem II phosphorylation in plants is required for adjustment of macroscopic folding of large photosynthetic membranes modulating lateral mobility of membrane proteins and sustained photosynthetic activity.

Analysis of the chloroplast protein kinase Stt7 during state transitions.

State transitions allow for the balancing of the light excitation energy between photosystem I and photosystem II and for optimal photosynthetic activity when photosynthetic organisms are subjected to changing light conditions. This process is regulated by the redox state of the plastoquinone pool through the Stt7/STN7 protein kinase required for phosphorylation of the light-harvesting complex LHCII and for the reversible displacement of the mobile LHCII between the photosystems. We show that Stt7 is associated with photosynthetic complexes including LHCII, photosystem I, and the cytochrome b6f complex. Our data reveal that Stt7 acts in catalytic amounts. We also provide evidence that Stt7 contains a transmembrane region that separates its catalytic kinase domain on the stromal side from its N-terminal end in the thylakoid lumen with two conserved Cys that are critical for its activity and state transitions. On the basis of these data, we propose that the activity of Stt7 is regulated through its transmembrane domain and that a disulfide bond between the two lumen Cys is essential for its activity. The high-light-induced reduction of this bond may occur through a transthylakoid thiol-reducing pathway driven by the ferredoxin-thioredoxin system which is also required for cytochrome b6f assembly and heme biogenesis.