Your new experience awaits. Try the new design now and help us make it even better

METHODS article

Front. Genet., 26 September 2025

Sec. Cancer Genetics and Oncogenomics

Volume 16 - 2025 | https://doi.org/10.3389/fgene.2025.1674138

CNAdjust: enhancing CNA calling accuracy through systematic baseline adjustment

  • 1Department of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
  • 2Computational Oncogenomics Group, Swiss Institute of Bioinformatics, Zurich, Switzerland

Accurate determination of the genomic copy number baseline is crucial for identifying copy number alterations (CNAs) in cancer, yet it remains a significant challenge in tumors with complex karyotypes. To address this, we present CNAdjust, an integrated method to systematically detect and correct baseline inaccuracies in CNA data. CNAdjust employs a Bayesian framework that integrates cohort-specific CNA frequency priors with a data-driven plausibility score, ensuring that adjusted calls align with both biological cohort patterns and study-specific data. Performance validation using the TCGA pan-cancer dataset demonstrated improved alignment with absolute copy number estimates and enhanced CNA pattern interpretation. Furthermore, we revealed a strong correlation between chromosomal aneuploidy and baseline abnormalities, underscoring the prevalence of this issue in cancer genomics. By systematically improving the precision of CNA calls, CNAdjust serves as a critical tool for constructing harmonized reference datasets and advancing the progress of precision oncology. Its implementation as a standard, portable workflow enables the reproducible and scalable analysis of large, heterogeneous datasets, supporting large-scale genomic research. Source codes are available at: https://github.com/baudisgroup/CNAdjust.

1 Introduction

The term “Copy Number Alterations” (CNAs1) refers to changes in the copy number of genomic segments, representing deviations from the normal karyotype. CNAs are a hallmark of cancer cells and play a crucial role in cancer development (Sansregret et al., 2018). These variations are typically expressed as relative copy number ratios between test and reference samples. Normalizing this relative ratio is a critical step in data processing, as it corrects systematic technical variations and enables meaningful biological comparisons (Quackenbush, 2002).

For cancer cells with extensive copy number imbalances, traditional normalization methods, such as median or lowess normalization, often fail to maintain an accurate baseline. These shifts from the expected log ratio (logR) of 0 complicate relative CNA calling and lead to inaccuracies in downstream analyses. While density-based adjustments offer a solution (Staaf et al., 2007; Marzouka et al., 2016), they rely on the assumption that most genomic markers cluster near the true baseline—an assumption that fails in cancers with pervasive CNAs, such as adrenocortical tumors (Ronchi et al., 2013). Even more sophisticated single-sample calibration tools ultimately rely on similar assumptions, attempting to identify an anchor baseline from the sample’s internal data alone, implicitly relying on diploid, allelic-balanced regions or marker density dominance (Fu et al., 2015; Gao and Baudis, 2020). In highly aneuploid tumors, however, these assumptions are frequently violated, leading to systematic baseline errors.

The resulting unreliability of CNA detection tools is a well-documented problem. A comprehensive benchmarking study (Luo, 2019), for example, revealed that no single method performs consistently well across both near-diploid and highly aneuploid tumors, highlighting the critical need for self-adaptive strategies that can adjust to varying karyotypes. Moreover, the difficulty is not merely technical but also conceptual, especially in aneuploid samples. For instance, in a genome where half the regions have three copies and the other half have 2, setting the baseline at the average ploidy of 2.5 would label the entire genome as altered. While seemingly accurate in isolation, this becomes problematic if prior biological knowledge suggests a diploid background, as this incorrect baseline diminishes the significance of true alterations and obscures meaningful patterns.

To address these challenges, we introduce CNAdjust, a novel method designed to automatically detect and correct baseline abnormalities in tumor CNA profiles. CNAdjust employs a Bayesian framework that systematically evaluates potential baseline corrections. It integrates prior probabilities derived from cohort-level CNA patterns with a data-driven plausibility score calculated from the intra-study samples’ logR value distribution. This synergistic approach ensures that adjusted CNA calls are consistent with both broad biological context and study-specific evidence, leading to more accurate and reliable interpretations. Built on the Nextflow workflow management system, CNAdjust offers robust reproducibility, parallelism, and portability, making it an effective tool for calibrating large-scale and heterogeneous CNA datasets for oncogenomic studies.

2 Methods

2.1 Input data

Our method is built on the observation that samples from the same cohort, such as those sharing a tumor (sub)type, typically display similar CNA patterns (Ni et al., 2013; Komori, 2022). Therefore, to refine the accuracy of CNA calls, the workflow requires external information, including cohort assignments for each sample and prior probabilities of gain and loss specific to each cohort (Figure 1a).

Figure 1
A detailed overview of the CNAdjust workflow for systematic baseline correction. Panel a presents a flowchart illustrating the entire process, starting from the input of segmented CNA data, through the extraction of features to identify problematic profiles, to the core adjustment step. This adjustment is guided by a Bayesian framework that leverages external information (prior probabilities) to decide if a proposed baseline shift improves the CNA calls. Panel b visualizes this Bayesian decision for a single genomic segment. It displays the segment’s logR value and the two critical inputs: the cohort-specific prior probability and the study-specific, data-driven plausibility score.

Figure 1. Overview of the workflow. (a) Execution steps, with segment colors indicating CNA states: blue for deletions, orange for duplications, and green for copy-neutral states. (b) Decision-making process for a segment identified as a duplication within a sample profile. In the middle prior pattern plot, dashed lines represent the regions spanned by the segment. In the bottom weighted logR distribution plot (by segment length), the dashed line indicates the logR value of the segment.

To demonstrate the applicability of our approach, we analyzed 9,846 masked copy number segment profiles from The Cancer Genome Atlas (TCGA), spanning 33 tumor types for which matched absolute copy number data was available. Two distinct sets of CNA calls are generated from this data set: the first derived from a straightforward cutoff approach using a threshold of 0.1, and the second using the same threshold but integrated with a density-based baseline correction using CopyNumber450kCancer (Marzouka et al., 2016). All samples from the same project were assigned to the same cohort group. To establish the prior probabilities essential for our model, we utilized CNA frequency data from Progenetix (Cai et al., 2012; Huang et al., 2021; Baudis Group, 2024), a comprehensive database that aggregates genomic mutation profiles from various sources, including Gene Expression Omnibus (GEO) (Barrett et al., 2012), cBioPortal (Cerami et al., 2012), and TCGA. The Progenetix database provides computed frequencies of CNAs for 1 MB genomic bins across samples from the same cohorts, thus providing the necessary prior information to reflect the CNA landscapes of specific populations. In our analysis, we used the National Cancer Institute Thesaurus (NCIt) terminology (National Cancer Institute, 2024) to match NCIt codes to different TCGA projects. Based on these NCIt codes, we retrieved the corresponding CNA frequency data from Progenetix, ensuring that the frequency values were derived from cohorts without any or containing only a small minority fraction of TCGA samples. Detailed information on the TCGA data and the corresponding NCIt codes can be found in the Supplementary Section 2 and Supplementary Table S2.

2.2 Workflow steps

The initial step is to individually analyze the segment profiles to identify those with potentially abnormal baselines. Aberrant baselines are indicated by pervasive CNA calls throughout the genome or extreme skewness in the CNA ratio. To detect these, this workflow extracts key features including the overall fraction of CNAs, the standard deviation of the segment value (logR) weighted by the number of markers, and the ratio of CNA fractions between duplications and deletions. The criteria for flagging problematic samples are based on an examination of 36 randomly selected GEO series that covered a variety of disease types and measurement platforms. For a detailed account of the test series and the derivation of default parameters, see the Supplementary Section 1. Our pipeline uses these default parameters, but they are designed to be user-adjustable.

For profiles identified with potential baseline issues, CNAdjust systematically corrects CNA calls by iteratively shifting the baseline upward or downward (Figure 1a). This core function is powered by our developed method, labelSeg (Zhao and Baudis, 2024), which uses unsupervised clustering to group segments based on their logR values. This clustering-based design makes the entire CNAdjust framework inherently adaptive to variations in tumor purity. Because purity primarily affects the amplitude of logR signals, the data-driven clustering identifies CNA states based on each sample’s own signal distribution, removing the need for an explicit purity prior. Rather than applying a single, fixed correction, CNAdjust conducts a systematic evaluation. The process begins by using the original calls to identify the initial baseline cluster. Depending on the detected bias, CNAdjust then directs labelSeg to select an adjacent cluster as a candidate for the new baseline, generating a new set of calls. This process can be repeated to test multiple candidate shifts.

Each set of candidate calls is then evaluated to determine which offers the greatest improvement. This decision-making process is automated by a Bayesian framework, as illustrated in Figure 1b and formalized in Equations 1,2. Equation 1 shows the general principle of Maximum a Posteriori (MAP) estimation, where the most probable set of CNA states, θ̂MAP is identified. To make this practical, we define an objective function, J(θ), shown in Equation 2. This function assumes independence between the n genomic segments in a profile and is calculated by summing two components for each segment: a log-prior-probability and a log-plausibility-score. The prior probability, p(θi|ri,c) reflects the expected occurrence of a CNA state θi (duplication, deletion, or neutrality) in genomic region ri for a given cohort c. For demonstration purposes, we use CNA frequency data from Progenetix as our prior; however, the framework supports integrating any belief about CNA occurrence. The plausibility score, S(di|θi,s), substitutes for a formal likelihood and is derived from the tail probability of the observed logR value, di, within the state-specific Gaussian distribution fitted on samples from the same study s. The final CNA calls are determined by finding the set of states θ that maximizes the objective function J(θ), thus ensuring that the chosen calls reflect the most probable states given both the observed data and our prior knowledge. If no adjusted set is deemed superior to the original, the original calls are retained.

θ̂MAP=argmaxθpθ|DargmaxθpD|θpθargmaxθlogpD|θpθ(1)
Jθ=i=1nlogSdi|θi,s+logpθi|ri,c(2)

3 Results

3.1 Performance validation

To validate the performance of our approach and emphasize its utility, we analyzed TCGA pan-cancer segment data using two distinct methods for generating input callings: logR-based calling, which applies a direct 0.1 cutoff value and assumes the theoretical baseline (0) as correct, and density-based calling, which employs a density-based baseline correction using CopyNumber450kCancer (Marzouka et al., 2016) before applying the same cutoff, assuming the true baseline corresponds to the region where the majority of markers are located. These input groups provide insights into diverse baseline normalization challenges, as depicted in Figure 2a.

Figure 2
Panel a conceptually illustrates four challenging CNA profile types and includes a table summarizing the effectiveness of different calling strategies. Panel b presents the quantitative results of applying CNAdjust. The histograms on the left show the distribution of logR values for the problematic samples, often deviating from the expected peak at zero. The central boxplots show that the adjustment significantly reduces the overall fraction of the genome called as altered, particularly for logR-based inputs. The violin plots on the right demonstrate a marked improvement in accuracy, showing that after adjustment (orange), the ploidy-corrected copy number values for deleted, neutral, and duplicated segments become more distinct and better aligned with their expected values compared to the original, ambiguous calls (blue).

Figure 2. Performance Validation. (a) Different types of segment profiles. The dashed line represents the theoretical baseline. The accompanying table on the right summarizes the performance of various calling strategies across these profile types, emphasizing their challenges and the effectiveness of CNAdjust in addressing them. Checkmarks indicate reliable performance, while crosses denote problematic situations. (b) Samples with baseline issues. The top and bottom panels correspond to the logR-based and density-based calling groups, respectively. The left plot shows weighted histograms of logR values (weighted by segment length). The middle plot illustrates CNA fractions; statistical significance was determined using a paired t-test with Benjamini–Hochberg correction. The right plot depicts ploidy-corrected copy changes across different calling states; statistical significance was determined using a t-test with Benjamini–Hochberg correction.

In cases with a correctly positioned baseline and moderate CNA load, the logR-based and density-based methods are expected to produce largely concordant results, as the core assumptions of both approaches are met. However, complications arise when these assumptions are violated. The accuracy of these calls can, of course, still be affected by other factors such as data quality, but our focus here is on discrepancies arising from baseline estimation. The logR-based calling performs well with correct baselines but struggles with shifted baselines, while the density-based calling is robust to baseline shifts but fails when the assumption of marker clustering around the true baseline is invalid. For such challenging profiles, CNAdjust enhances calling quality by aligning results with higher posterior probabilities, producing baseline-corrected outputs. Substantial discrepancies between adjusted and original calls serve as indicators of baseline abnormalities in the original profiles. In this study, the calling state for segments spanning at least half of a chromosome arm was considered as the arm-level calling state, and samples with at least one arm-level discrepancy between original and adjusted calls were classified as having baseline issues.

3.1.1 Closer alignment with absolute ploidy-corrected copy changes

Our validation, illustrated in Figure 2b, focused on these samples with baseline issues. Specifically, 447 samples failed under the logR-based approach, 2,035 failed under the density-based approach, and 202 failed under both methods. Analysis of segment value distributions revealed deviations from the expected maximum signal peak of 0 in the logR-based group, consistent with the third and fourth segment profile types depicted in Figure 2a. Similarly, the emergence of multiple major peaks in the density-based group aligns with expectations, as shown in the second and fourth profile types in Figure 2a. Additionally, other non-adjusted samples exhibited a distinct normal distribution of segment values (Supplementary Figure S2), further supporting the necessity and efficacy of CNAdjust. These findings underscore the unique challenges of each calling method and demonstrate the robust ability of CNAdjust to resolve baseline issues effectively.

Application of the workflow to the different input calling methods yielded distinct effects on the overall CNA fraction, as illustrated in the middle panel of Figure 2b. For the logR-based group, the adjustment led to a universally pronounced decrease in the CNA fraction. In contrast, the effect on the density-based group was more project-dependent. For example, while adjustments affected over 30% of calls in the TCGA-SARC and TCGA-CHOL datasets for the logR-based group, the impact on the density-based group was more variable: in the TCGA-KICH and TCGA-ACC projects, over 50% of samples experienced an increase in their CNA fraction post-adjustment, whereas in the TCGA-LUSC, TCGA-UCS, and TCGA-ESCA projects, over 30% of samples showed a decrease.

Validation against absolute copy number estimates further confirmed the effectiveness of our adjustment strategy (Figure 2b, right panel). For this comparison, we used absolute copy number and ploidy estimates generated for the same samples by the ABSOLUTE algorithm (Carter et al., 2012). To ensure consistency, these absolute estimates were summarized at the chromosome-arm level by calculating weighted averages based on segment lengths. In problematic samples, the original calling states exhibited ambiguous copy change amplitudes that lacked clear alignment with these absolute estimates. After adjustment, our workflow significantly enhanced the clarity of ploidy-corrected copy change amplitudes, achieving a much closer alignment with the absolute copy number data. In contrast, non-adjusted samples retained distinct ploidy-corrected copy change patterns across different calling states (Supplementary Figure S2), underscoring the value of the baseline correction.

To confirm that these improvements are driven by the integration of biologically accurate information, we performed a control experiment using a “least-likely prior.” This counter-evidential prior was constructed by identifying, for each genomic region, the CNA state with the lowest probability in the original cohort data and assigning it a very high probability (0.98), while the other two states were assigned low probabilities (0.01). When CNAdjust was run with this actively misleading prior, performance degraded catastrophically, as illustrated in Supplementary Figure S3. Instead of reducing the overall CNA fraction, the adjustment process now significantly increased it across both input methods. Furthermore, the alignment with absolute copy number estimates was largely lost. This result validates that the success of CNAdjust relies on the correct interplay between the cohort-specific prior and the study-specific data, and is not an artifact of the adjustment framework itself.

3.1.2 Enhanced CNA pattern interpretation

To illustrate the improved interpretation achieved through baseline adjustment, we examined data from the TCGA testicular germ cell tumor (TCGA-TGCT) project specifically. Among the 155 samples analyzed, 91 displayed consistent CNA callings across both input groups after adjustment, as these samples did not exhibit baseline issues. Despite originating from distinct calling strategies, these samples showed highly similar CNA patterns (Figure 3b). Notably, key CNA features, such as deletions on chromosomes 4, 5, 11, 13q, and 18, along with duplications on chromosomes 7, 8, 12p, and 21q, were readily identifiable. These features aligned closely with prior CNA patterns from Progenetix used for adjustment (Figure 3a).

Figure 3
A case study on the TCGA Testicular Germ Cell Tumor (TGCT) cohort, demonstrating how CNAdjust enhances biological pattern interpretation. Panel a displays the reference CNA frequency plot for TGCT from the Progenetix database, which serves as the biological prior. Panel b shows frequency plots for two subgroups of TCGA samples: those with difficult baselines and those with easy baselines. For the difficult group, the “Adjusted“ plot successfully restores the characteristic CNA pattern, making it highly similar to both the easy group and the reference prior. Panel c provides a single-sample view, contrasting a difficult sample before and after adjustment with a stable, easy sample, visually confirming the correction of significant baseline errors.

Figure 3. Effects of adjustment in the TCGA-TGCT cohort. (a) The prior CNA frequency plot for testicular germ cell tumors, sourced from the Progenetix database. (b) CNA frequency plots from TCGA-TGCT samples, categorized by input calling method (rows) and baseline quality (columns). The top and bottom rows correspond to the logR-based and density-based methods, respectively. The samples are grouped into 13 with “difficult baselines” (left, corresponding to profile type 4 in Figure 2a) and 91 with “easy baselines” (right, profile type 1). (c) Individual segment profiles for representative samples from the groups in panel (b). The left column shows a difficult sample before and after adjustment, illustrating the correction of baseline errors. The right column shows an easy sample for comparison. CNA states are indicated by color: blue for deletions, orange for duplications, and green for copy-neutral states.

In contrast, 13 samples with baseline issues posed calling challenges, as both methods initially failed to produce reliable results. Their original CNA patterns deviated substantially from those of consistent samples, with the two calling methods yielding conflicting outputs: the logR-based input tended to overestimate total CNAs, while the density-based input introduced a duplication bias. These discrepancies obscured key characteristic patterns, complicating biological interpretation.

Following adjustment, our workflow resolved these inconsistencies by integrating prior cohort-level CNA patterns with study-specific logR distributions (Figures 3b,c). The adjusted CNA patterns aligned more closely with expected biological characteristics, effectively revealing features such as the hallmark 12p duplication—a genetic signature of testicular germ cell tumors (Batool et al., 2019)—that had been previously obscured. These results underscore the capability of CNAdjust in handling complex cases and provide more accurate and biologically meaningful CNA interpretations.

3.2 Aneuploidy’s impact on baseline ambiguity

To investigate the influence of aneuploidy on baseline normalization challenges, we analyzed its association with the baseline abnormalities identified by our method. Samples with ploidy 2.5 or 1.5 were classified as aneuploid. We then tested for a statistical association between aneuploidy and baseline ambiguity within each TCGA project. As summarized in Figure 4, we observed that samples with challenging baselines (red bars) consistently had a higher proportion of aneuploidy than those with stable baselines (blue bars). This association between aneuploidy and baseline ambiguity was statistically significant in 31 of the 33 projects. These findings underscore the critical impact of aneuploidy on baseline accuracy and highlight the need for adjustment methods, like CNAdjust, that can address complex karyotypic variations.

Figure 4
Bar chart illustrating the strong statistical association between aneuploidy and baseline ambiguity across 33 different TCGA cancer projects. The y-axis quantifies the proportion of samples within a group that are aneuploid. For each project on the x-axis, the chart compares two groups: samples identified by CNAdjust as having “challenging baselines” (red bars) and those considered “baseline-stable” (blue bars). The consistently and significantly higher red bars across nearly all projects provide strong evidence that samples with complex, aneuploid genomes are far more likely to suffer from baseline calling errors, underscoring the critical need for methods like CNAdjust.

Figure 4. Association between aneuploidy and baseline complexity. The Y-axis represents the proportion of aneuploid samples in each group. Samples with arm-level calling differences post-adjustment, regardless of the calling strategy, were categorized as having baseline issues (red), while other samples were considered baseline-stable (blue). P-values were calculated using a Chi-squared test, except for the 13 rightmost projects where Fisher’s exact test was used due to low event counts.

4 Discussion

CNAdjust provides an automated and systematic solution to address baseline inaccuracies in CNA calling, a prevalent and often overlooked challenge in cancer genomics. This method is designed to enhance calling accuracy and effectively handle the complexities of large, heterogeneous datasets. Its robustness was demonstrated through its application to the TCGA pan-cancer dataset, where we analyzed around 10,000 samples across 33 tumor types derived from divergent calling methodologies. This strategy highlighted various instances of ambiguous baselines and validated the fundamental logic of our approach: samples flagged for baseline issues exhibited abnormal logR distributions, and our adjustments led to improved alignment with absolute copy number estimates. A focused analysis of the TCGA-TGCT cohort further underscored CNAdjust’s utility, showing how it can restore characteristic CNA patterns obscured by baseline errors. Finally, the strong association we found between aneuploidy and baseline ambiguity across nearly all TCGA projects emphasizes that addressing karyotypic complexity is critical for accurate cancer CNA studies.

Despite its strengths, CNAdjust has limitations. A key assumption is that segment values for a given CNA state are drawn from a common distribution within a study. This may not hold if a study contains highly aberrant samples that do not appear as statistical outliers, which could complicate the adjustment. To mitigate this, our method incorporates a reference logR distribution when more than 50% of a study’s samples are flagged as abnormal and adjusts the fitted distributions when this fraction exceeds 25% (detailed in Supplementary Section 3).

Another consideration is the reliance on appropriate prior selection, which raises a valid concern about intra-cohort heterogeneity. It is true that applying a single, genome-wide prior to a diverse cancer cohort is a simplification, as molecular subtypes with distinct CNA patterns exist (Yang et al., 2024). However, it is crucial to emphasize that the prior in our Bayesian framework does not rigidly force profiles to conform to a cohort average. Instead, it acts as a probabilistic guide that is balanced against the data-driven plausibility score from the sample itself. For this pan-cancer study, we employed broad NCIt-based cohort definitions with frequency data from Progenetix (accessed via the pgxRpi (Zhao and Baudis, 2025) package) as a pragmatic and scalable approach to demonstrate the method’s utility. Recognizing that the relationship between CNA heterogeneity and tumor classification is not always straightforward, the CNAdjust framework was designed for flexibility. Its performance can be potentially enhanced by incorporating more granular priors, such as those defined by known molecular subtypes or other user-defined biological knowledge. This adaptability is a core strength, allowing for a more tailored and precise adjustment process as cohort definitions become more refined.

5 Conclusion

CNAdjust effectively improves the accuracy of relative CNA calls by systematically identifying and correcting baseline inaccuracies. This study underscores that proper baseline determination is critical for the accurate interpretation of CNA patterns, which in turn is essential for advancing our understanding of tumor biology and informing clinical applications. In the era of big data, the demand for scalable, efficient, and automated solutions like CNAdjust is paramount. By enhancing the quality and reliability of CNA analyses, our method represents a significant advancement towards the construction of harmonized reference datasets and provides a valuable tool to support the progress of large-scale genomic research and precision medicine.

Data availability statement

The datasets analyzed in this study are publicly available. The TCGA pan-cancer datasets used for validation were downloaded from the GDC Data Portal https://portal.gdc.cancer.gov/). The data used for method development, including the GEO series for parameter optimization and the cohort-specific CNA frequency data, were sourced from the Progenetix oncogenomic resource (https://progenetix.org). Further details on the specific TCGA projects and GEO series used are provided in the Supplementary Material.

Ethics statement

Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

Author contributions

HZ: Writing – original draft, Software, Visualization, Conceptualization, Methodology, Writing – review and editing, Validation. MB: Supervision, Writing – review and editing, Funding acquisition, Project administration, Resources.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that Generative AI was used in the creation of this manuscript. Generative AI was used to improve the grammar, clarity, and readability of the text. All scientific content, including the concepts, study design, methodology, results, and conclusions, was conceived and written by the authors.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2025.1674138/full#supplementary-material

Footnotes

1Alternatively the term Copy Number Variations (CNV) can be employed, potentially with added “somatic” specifier to disambiguate from inherited variants.

References

Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2012). Ncbi geo: archive for functional genomics data Sets—Update. Nucleic acids Res. 41 (D1), D991–D995. doi:10.1093/nar/gks1193

PubMed Abstract | CrossRef Full Text | Google Scholar

Batool, A., Karimi, N., Wu, X.-N., Chen, S.-R., and Liu, Y.-X. (2019). Testicular germ cell tumor: a comprehensive review. Cell. Mol. Life Sci. 76, 1713–1727. doi:10.1007/s00018-019-03022-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Baudis Group (2024). Progenetix.

Google Scholar

Cai, H., Kumar, N., and Baudis, M. (2012). Arraymap: a reference resource for genomic copy number imbalances in human malignancies. PLoS One 7 (5), e36944. doi:10.1371/journal.pone.0036944

PubMed Abstract | CrossRef Full Text | Google Scholar

Carter, S. L., Cibulskis, K., Helman, E., McKenna, A., Shen, H., Zack, T., et al. (2012). Absolute quantification of somatic dna alterations in human cancer. Nat. Biotechnol. 30 (5), 413–421. doi:10.1038/nbt.2203

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerami, E., Gao, J., Dogrusoz, U., Gross, B. E., Sumer, S. O., Aksoy, B. A., et al. (2012). The cbio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2 (5), 401–404. doi:10.1158/2159-8290.CD-12-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, Yi, Yu, G., Levine, D. A., Wang, N., Shih, I.-M., Zhang, Z., et al. (2015). Bacom2. 0 facilitates absolute normalization and quantification of somatic copy number alterations in heterogeneous tumor. Sci. Rep. 5 (1), 13955. doi:10.1038/srep13955

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, Bo, and Baudis, M. (2020). Minimum error calibration and normalization for genomic copy number analysis. Genomics 112 (5), 3331–3341. doi:10.1016/j.ygeno.2020.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Q., Carrio-Cordo, P., Gao, Bo, Paloots, R., and Baudis, M. (2021). The progenetix oncogenomic resource in 2021. Database 2021, baab043. doi:10.1093/database/baab043

PubMed Abstract | CrossRef Full Text | Google Scholar

Komori, T. (2022). Grading of adult diffuse gliomas according to the 2021 who classification of tumors of the central nervous system. Lab. Investig. 102 (2), 126–133. doi:10.1038/s41374-021-00667-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, F. (2019). A systematic evaluation of copy number alterations detection methods on real snp array and deep sequencing data. BMC Bioinforma. 20 (Suppl. 25), 692. doi:10.1186/s12859-019-3266-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Marzouka, N. al-dain, Nordlund, J., Bäcklin, C. L., Lönnerholm, G., Syvänen, A.-C., and Almlöf, J. C. (2016). Copynumber450kcancer: baseline correction for accurate copy number calling from the 450k methylation array. Bioinformatics 32 (7), 1080–1082. doi:10.1093/bioinformatics/btv652

PubMed Abstract | CrossRef Full Text | Google Scholar

National Cancer Institute (2024). Nci thesaurus.

Google Scholar

Ni, X., Zhuo, M., Su, Z., Duan, J., Gao, Y., Wang, Z., et al. (2013). Reproducible copy number variation patterns among single circulating tumor cells of lung cancer patients. Proc. Natl. Acad. Sci. 110 (52), 21083–21088. doi:10.1073/pnas.1320659110

PubMed Abstract | CrossRef Full Text | Google Scholar

Quackenbush, J. (2002). Microarray data normalization and transformation. Nat. Genet. 32 (4), 496–501. doi:10.1038/ng1032

PubMed Abstract | CrossRef Full Text | Google Scholar

Ronchi, C. L., Sbiera, S., Leich, E., Henzel, K., Rosenwald, A., Allolio, B., et al. (2013). Single nucleotide polymorphism array profiling of adrenocortical tumors-evidence for an adenoma carcinoma sequence? PloS one 8 (9), e73959. doi:10.1371/journal.pone.0073959

PubMed Abstract | CrossRef Full Text | Google Scholar

Sansregret, L., Vanhaesebroeck, B., and Swanton, C. (2018). Determinants and clinical implications of chromosomal instability in cancer. Nat. Rev. Clin. Oncol. 15 (3), 139–150. doi:10.1038/nrclinonc.2017.198

PubMed Abstract | CrossRef Full Text | Google Scholar

Staaf, J., Jönsson, G., Ringner, M., and Vallon-Christersson, J. (2007). Normalization of array-cgh data: influence of copy number imbalances. BMC genomics 8, 1–18. doi:10.1186/1471-2164-8-382

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Carrio-Cordo, P., and Baudis, M. (2024). Copy number variation heterogeneity reveals biological inconsistency in hierarchical cancer classifications. Mol. Cytogenet. 17 (1), 26. doi:10.1186/s13039-024-00692-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H., and Baudis, M. (2024). Labelseg: segment annotation for tumor copy number alteration profiles. Briefings Bioinforma. 25 (2), bbad541. doi:10.1093/bib/bbad541

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, H., and Baudis, M. (2025). Pgxrpi: an r/bioconductor package for user-friendly access to the beacon v2 Api. Bioinforma. Adv. 5 (1), vbaf172. doi:10.1093/bioadv/vbaf172

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: copy number alterations, baseline correction, bayesian framework, cancer genomics, nextflow workflow

Citation: Zhao H and Baudis M (2025) CNAdjust: enhancing CNA calling accuracy through systematic baseline adjustment. Front. Genet. 16:1674138. doi: 10.3389/fgene.2025.1674138

Received: 27 July 2025; Accepted: 10 September 2025;
Published: 26 September 2025.

Edited by:

Nicoletta Coccaro, University of Bari Aldo Moro, Italy

Reviewed by:

Stefan Kirov, Flare Therapeutics Inc., United States
Seokhyun Yoon, Dankook University, Republic of Korea

Copyright © 2025 Zhao and Baudis. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hangjia Zhao, aGFuZ2ppYS56aGFvQHV6aC5jaA==; Michael Baudis, bWljaGFlbC5iYXVkaXNAbWxzLnV6aC5jaA==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.