HRDetect

BRCA1 and BRCA2 play crucial roles in maintaining genome integrity, mainly through homologous recombination (HR) for DNA double-strand breaks (DSB)repair. The mutations of BRCA1 and BRCA2 can lead to a reduced capacity of HR machinery, increased genomic instability, and elicit a predisposition to malignancies.[2] People with BRCA1 and BRCA2 deficiency have higher risks of developing certain cancers such as breast and ovarian cancers. Germline defects in BRCA1/BRCA2 genes account for up to 5% of breast cancer cases.[1]

Diagram illustrating synthetic lethality[3]

Poly (ADP ribose) polymerase (PARP) inhibitors are designed to treat BRCA1- and BRCA2- defect tumors owing to their homologous recombination deficiency.[4] These drugs have been majorly implemented in breast and ovarian cancers,[5] and their clinical efficacy among patients with other types of cancers, such as pancreatic cancer, is still being investigated.[6] It is vital to identify adequate patients with BRCA1/BRCA2 deficiency to utilize PARP inhibitors optimally. PARP inhibitors operate on the concept of synthetic lethality where by selectively causing cell death in BRCA-mutant cells while sparing normal cells.

HRDetect was implemented to detect tumors with BRCA1/BRCA2 deficiency using the data from whole-genome sequencing. This model quantitatively aggregates six HRD-associated signatures into a single score called HRDetect to accurately classify breast cancers by their BRCA1 and BRCA2 status. The machine learning algorithm assigns weight values to these signatures prior to computing the final score. The six signatures, ranked by decreasing weight, include microhomology-mediated indels, the HRD index, base- substitution signature 3, rearrangement signature 3, rearrangement signature 5, and base- substitution signature 8. Additionally, this weighted approach is able to identify BRCAness, which refers to mutational phenotypes displaying homologous recombination deficiency similar to tumors with BRCA1/BRCA2 germline defects.[7]

HRDetect requires four types of inputs:

1. Counts of mutations associated with each signature of single-base substitutions
2. Indels with microhomology at the indel breakpoint junction, indels at polynucleotide-repeat tracts and other complex indels as proportions
3. Counts of rearrangements associated with each signature
4. HRD index (Arithmetic sum of loss of heterozygosity (LOH), telomeric-allelic imbalance (TAI), and large-scale state transitions (LST) scores)

It is based on a supervised learning method using a lasso logistic regression model to distinguish samples into those with and without BRCA 1/2 deficiency. Optimal coefficients are obtained by minimizing the objective function.

To account for a high substitution count in samples, the genomic data is first log transformed:

${\displaystyle x=\ln(x+1)}$

The transformed data is then standardized to make mutational class values comparable giving each object a mean of 0 and a standard deviation (sd) of 1:

${\displaystyle \mathrm {x} ={\frac {x-\operatorname {mean} \left(x\right)}{\mathrm {s} \mathrm {d} \cdot \left(x\right)}}}$

Weight distributions of input signatures

To be able to distinguish between those affected and not affected by BRCA1/BRCA2 deficiency, a lasso logistic regression model is used:

${\displaystyle \min _{((\beta _{0},\,\beta ))\in \mathbb {R} ^{p+1}}{\left(-\left[{\frac {1}{N}}\sum _{i=1}^{N}y_{i}\cdot \left(\beta _{0}+x_{i}^{T}\beta \right)-\log \left(1+e^{\left(\beta _{0}+x_{i}^{T}\beta \right)}\right)\right]+\lambda \|\beta \|_{1}\right)}}$

where:

${\displaystyle y_{i}}$: BRCA status of a sample || y_i = 1 for BRCA1/BRCA2-null samples || y_i = 0 otherwise
${\displaystyle \beta _{0}}$: Intercept, interpreted as the log of odds of ${\displaystyle y_{i}}$ = 1 when ${\displaystyle x_{i}^{T}}$ = 0
${\displaystyle \beta }$: Vector of weights
${\displaystyle p}$: Number of features characterizing each sample
${\displaystyle N}$: Number of samples
${\displaystyle x_{i}^{T}}$: Vector of features characterizing the i^th sample
${\displaystyle \lambda }$: Penalty promoting the sparseness of the weights
${\displaystyle \|\beta \|}$: L1 norm of the vector of weights

The β weights are constrained to be positive to reflect the presence of mutational actions due to BRCA1/BRCA2 defects. Setting the constraint of nonnegative weights ensures that all samples would be scored on the basis of the presence of relevant mutational signatures associated with BRCA1/BRCA2 deficiency, irrespective of whether these signatures are the dominant mutational process in the cancer.

Lastly, the weights obtained from the lasso regression are used to give a new sample a probabilistic score using the normalized mutational data ${\displaystyle x_{i}^{T}}$and application of the model parameters(${\displaystyle \beta }$, ${\displaystyle \beta _{0}}$):

${\displaystyle P\left(C_{i}=BRCA\right)={\frac {1}{1+e^{-\left(\beta _{0}+x_{i}^{T}\beta \right)}}}}$

where:

${\displaystyle C_{i}}$ : variable encoding the status of the i^th sample
${\displaystyle \beta _{0}}$ : Intercept weight
${\displaystyle x_{i}^{T}}$: Vector encoding features of the i^th sample
${\displaystyle \beta }$: Vector of weights

The probability value quantifies the degree of BRCA1/BRCA2 defectiveness. A cut-off probability value should be chosen while maintaining a high sensitivity. These scores can be utilized to guide therapy.

Mutations in genes responsible for HR are prevalent among human cancers. The BRCA1 and BRCA2 genes are centrally involved in HR, DNAdamage repair, end resection, and checkpoint signaling. Mutational signatures of HRD have been identified in over 20% of breast cancers, as well as pancreatic, ovarian, and gastric cancers. BRCA1/2 mutations confer sensitivity to platinum-based chemotherapies. HRDetect can independently trained to predict BRCA1/2 status, and has the capacity to predict outcomes on platinum-based chemotherapies.[8]

HRDetect was initially developed to detect tumors with BRCA1 and BRCA2 deficiency based on the data from whole-genome sequencing of a cohort of 560 breast cancer samples. Within this cohort, 22 patients were known to carry germline BRCA1/BRCA2 mutations. BRCA1/BRCA2- deficiency mutational signatures were found in more breast cancer patients than previously known. This model was able to identify 124 (22%) breast cancer patients showing BRCA1/2 mutational signatures in this cohort of 560 samples. Apart from the 22 known cases, an additional 33 patients showed deficiency with germline BRCA1/2 mutations, 22 patients displayed somatic mutation of BRCA1/2, and 47 were recognized to show functional defect without detected BRCA1/2 mutation. As a result, with an application of a probabilistic cut-off 0.7, HRDetect was able to demonstrate a 98.7% sensitivity recognizing BRCA1/2- deficient cases.

In contrast, germline mutations of BRCA1/2 are present in only 1~5% of breast cancer cases. Furthermore, these findings suggest that more breast cancer patients, as many as 1 in 5 (20%), may benefit from PARP inhibitors than a small percentage of patients currently given with the treatment. Cohort of 80 Breast cancer patients. 6 out of 7 are above HRDetect score 0.7.

Cohort of 80 Breast Cancer Samples

HRDetect was tested in 80 breast cancer cases with mainly ER positive and HER2 negative. The tool was able to find ones that exceed HRDetect score 0.7, including one germline BRCA1 mutation carrier, four germline BRCA2 mutation carriers and one somatic BRCA2 mutation carrier. The sensitivity of this tool also reached 86%.

HRDetect can be applied to other cancer types and yields adequate sensitivity.[6]

In a cohort of 73 patients with ovarian cancer, 30 patients were known to carry BRCA1/BRCA2 mutations and 46 (63%) patients were assessed by HRDetect to have HRDetect score over 0.7. The sensitivity of detecting BRCA1/2-deficient cancer was almost 100%, with an additional 16 cases identified.

In a cohort of 96 patients with pancreatic cancers, 6 cases were known to have mutation or allele loss and 11 (11.5%) patients were identified by HRDetect to an exceed cutoff of 0.7. The study observed a similar result of sensitivity approaching 100%, with five other cases identified.