Response to chemotherapy, recurrence and survival in advanced-stage ovarian, fallopian tube and primary peritoneal cancer patients with non-Ashkenazi Jewish BRCA mutations, compared to those without

Objectives. Previous studies have shown a survival advantage in ovarian cancer patients with Ashkenazi-Jewish (AJ) BRCA founder mutations, compared to sporadic ovarian cancer patients. The purpose of this study was to determine if this association exists in ovarian cancer patients with non-Ashkenazi Jewish BRCA mutations. In addition, we sought to account for possible "survival bias" by minimizing any lead time that may exist between diagnosis and genetic testing. ^ Methods. Patients with stage III/IV ovarian, fallopian tube, or primary peritoneal cancer and a non-Ashkenazi Jewish BRCA1 or 2 mutation, seen for genetic testing January 1996-July 2007, were identified from genetics and institutional databases. Medical records were reviewed for clinical factors, including response to initial chemotherapy. Patients with sporadic (non-hereditary) ovarian, fallopian tube, or primary peritoneal cancer, without family history of breast or ovarian cancer, were compared to similar cases, matched by age, stage, year of diagnosis, and vital status at time interval to BRCA testing. When possible, 2 sporadic patients were matched to each BRCA patient. An additional group of unmatched, sporadic ovarian, fallopian tube and primary peritoneal cancer patients was included for a separate analysis. Progression-free (PFS) & overall survival (OS) were calculated by the Kaplan-Meier method. Multivariate Cox proportional hazards models were calculated for variables of interest. Matched pairs were treated as clusters. Stratified log rank test was used to calculate survival data for matched pairs using paired event times. Fisher's exact test, chi-square, and univariate logistic regression were also used for analysis. ^ Results. Forty five advanced-stage ovarian, fallopian tube and primary peritoneal cancer patients with non-Ashkenazi Jewish (non-AJ) BRCA mutations, 86 sporadic-matched and 414 sporadic-unmatched patients were analyzed. Compared to the sporadic-matched and sporadic-unmatched ovarian cancer patients, non-AJ BRCA mutation carriers had longer PFS (17.9 & 13.8 mos. vs. 32.0 mos., HR 1.76 [95% CI 1.13–2.75] & 2.61 [95% CI 1.70–4.00]). In relation to the sporadic- unmatched patients, non-AJ BRCA patients had greater odds of complete response to initial chemotherapy (OR 2.25 [95% CI 1.17–5.41]) and improved OS (37.6 mos. vs. 101.4 mos., HR 2.64 [95% CI 1.49–4.67]). ^ Conclusions. This study demonstrates a significant survival advantage in advanced-stage ovarian cancer patients with non-AJ BRCA mutations, confirming the previous studies in the Jewish population. Our efforts to account for "survival bias," by matching, will continue with collaborative studies. ^

ACCURACY OF THE BRCAPRO RISK ASSESSMENT MODEL IN MALES PRESENTING TO MD ANDERSON FOR BRCA TESTING

ACCURACY OF THE BRCAPRO RISK ASSESSMENT MODEL IN MALES PRESENTING TO MD ANDERSON FOR BRCA TESTING Publication No. _______ Carolyn A. Garby, B.S. Supervisory Professor: Banu Arun, M.D. Hereditary Breast and Ovarian Cancer (HBOC) syndrome is due to mutations in BRCA1 and BRCA2 genes. Women with HBOC have high risks to develop breast and ovarian cancers. Males with HBOC are commonly overlooked because male breast cancer is rare and other male cancer risks such as prostate and pancreatic cancers are relatively low. BRCA genetic testing is indicated for men as it is currently estimated that 4-40% of male breast cancers result from a BRCA1 or BRCA2 mutation (Ottini, 2010) and management recommendations can be made based on genetic test results. Risk assessment models are available to provide the individualized likelihood to have a BRCA mutation. Only one study has been conducted to date to evaluate the accuracy of BRCAPro in males and was based on a cohort of Italian males and utilized an older version of BRCAPro. The objective of this study is to determine if BRCAPro5.1 is a valid risk assessment model for males who present to MD Anderson Cancer Center for BRCA genetic testing. BRCAPro has been previously validated for determining the probability of carrying a BRCA mutation, however has not been further examined particularly in males. The total cohort consisted of 152 males who had undergone BRCA genetic testing. The cohort was stratified by indication for genetic counseling. Indications included having a known familial BRCA mutation, having a personal diagnosis of a BRCA-related cancer, or having a family history suggestive of HBOC. Overall there were 22 (14.47%) BRCA1+ males and 25 (16.45%) BRCA2+ males. Receiver operating characteristic curves were constructed for the cohort overall, for each particular indication, as well as for each cancer subtype. Our findings revealed that the BRCAPro5.1 model had perfect discriminating ability at a threshold of 56.2 for males with breast cancer, however only 2 (4.35%) of 46 were found to have BRCA2 mutations. These results are significantly lower than the high approximation (40%) reported in previous literature. BRCAPro does perform well in certain situations for men. Future investigation of male breast cancer and men at risk for BRCA mutations is necessary to provide a more accurate risk assessment.

Knowledge, Attitudes, and Utilization of BRCA Testing Among Obstetricians and Gynecologists

Hereditary breast and ovarian cancer (HBOC) is an inherited cancer syndrome that is associated with mutations in the BRCA1 and BRCA2 genes. Carriers of BRCA mutations, both men and women, are at an increased risk for developing certain cancers. Carriers are most notably at an increased risk to develop breast and ovarian cancers; however an increased risk for prostate cancer, melanoma, and pancreatic cancers has also been associated with these mutations. In 2009 the American Congress of Obstetricians and Gynecologists (ACOG) released a practice bulletin stating that evaluating a patient’s risk for HBOC should be a routine part of obstetric and gynecologic practice. A survey was created and completed by 83 obstetricians and gynecologists in the greater Houston, TX area. The survey consisted of four sections designed to capture demographic information, attitudes towards HBOC and BRCA testing, utilization of BRCA testing, and the overall knowledge of respondents with regards to HBOC and BRCA testing. This study found that the majority of participants indicated that they felt that obstetricians and gynecologists should have the primary responsibility of identifying patients who may be at increased risk of carrying a BRCA mutation. Moreover, this study found that the majority of participants indicated that they felt comfortable or very comfortable in identifying patients at an increased risk of carrying a BRCA mutation. However, only about a quarter of participants indicated that they order BRCA genetic testing one to two times per month or more. Lastly, this study demonstrates that the overall knowledge of HBOC and BRCA testing among this population of obstetricians and gynecologists is poor. The results of this study stress the need for more education regarding HBOC, genetic testing, and strategies for identifying patients that may be at risk for having a mutation in a BRCA gene. Furthermore, it reiterates the importance of raising awareness to current practice guidelines and recommendations that can assist obstetricians and gynecologist to better identify and manage patients that may be at an increased risk of having HBOC.

Factors Associated with Early Versus Late Development of Breast and Ovarian Cancer in BRCA1 and BRCA2 Positive Women

Hereditary breast and ovarian cancer (HBOC) is caused by a mutation in the BRCA1 or BRCA2 genes. Women with a BRCA1/2 mutation are at increased risks for breast and ovarian cancer and often develop cancer at an earlier age than the general population. However, some women with a BRCA1/2 mutation do not develop breast or ovarian cancer under the age of 50 years. There have been no specific studies on BRCA positive women with no cancer prior to age 50, therefore this study sought to investigate factors within these women with no cancer under age 50 with respect to reproductive risk factors, BMI, tumor pathology, screening history, risk-reducing surgeries, and family history. 241 women were diagnosed with cancer prior to age 50, 92 with cancer at age 50 or older, and 20 women were over age 50 with no cancer. Data were stratified based on BRCA1 and BRCA2 mutation status. Within the cohorts we investigated differences between women who developed cancer prior to age 50 and those who developed cancer at age 50 or older. We also investigated the differences between women who developed cancer at age 50 or older and those who were age 50 or older with no cancer. Of the 92 women with a BRCA1/2 mutation who developed cancer at age 50 or older, 46 developed ovarian cancer first, 45 developed breast cancer, and one had breast and ovarian cancer diagnosed synchronously. BRCA2 carriers diagnosed age 50 or older were more likely to have ER/PR negative breast tumors when compared to BRCA2 carriers who were diagnosed before age 50. This is consistent with one other study that has been performed. Ashkenazi Jewish women with a BRCA1 mutation were more likely to be diagnosed age 50 or older than other ethnicities. Hispanic women with a BRCA2 mutation were more likely to be diagnosed prior to age 50 when compared to other ethnicities. No differences in reproductive factors or BMI were observed. Further characterization of BRCA positive women with no cancer prior to age 50 may aid in finding factors important in the development of breast or ovarian cancer.

Predictors of Contralateral Breast Cancer in BRCA Negative Women

Breast cancer is the most common cancer diagnosis and second leading cause of death in women. Risk factors associated with breast cancer include: increased age, alcohol consumption, cigarette smoking, white race, physical inactivity, benign breast conditions, reproductive and hormonal factors, dietary factors, and family history. Hereditary breast and ovarian cancer syndrome (HBOC) is caused by mutations in the BRCA1 and BRCA2 genes. Women carrying a mutation in these genes are at an increased risk to develop a second breast cancer. Contralateral breast cancer is the most common second primary cancer in patients treated for a first breast cancer. Other risk factors for developing contralateral breast cancer include a strong family history of breast cancer, age of onset of first primary breast cancer, and if the first primary was a lobular carcinoma, which has an increased risk of being bilateral. A retrospective chart review was performed on a select cohort of women in an IRB approved database at MD Anderson Cancer Center. The final cohort contained 572 women who tested negative for a BRCA1 or BRCA2 mutation, had their primary invasive breast cancer diagnosed under the age of 50, and had a BRCAPro risk assessment number over 10%. Of the 572 women, 97 women developed contralateral breast cancer. A number of predictors of contralateral breast cancer were looked at between the two groups. Using univariable Cox Proportional Hazard model, thirteen statistically interesting risk factors were found, defined as having a p-value under 0.2. Multivariable stepwise Cox Proportional Hazard model found four statistically significant variables out of the thirteen found in the univariable analysis. In our study population, the incidence of contralateral breast cancer was 17%. Four statistically significant variables were identified. Undergoing a prophylactic mastectomy was found to reduce the risk of developing contralateral breast cancer, while not having a prophylactic mastecomy, a young age at primary diagnosis, having a positive estrogen receptor status of the primary tumor, and having a family history of breast cancer increased a woman’s risk to develop contralateral breast cancer.

Estimation of the incidence of a rare genetic disease through a two-tier mutation survey.

Recent attempts to detect mutations involving single base changes or small deletions that are specific to genetic diseases provide an opportunity to develop a two-tier mutation-screening program through which incidence of rare genetic disorders and gene carriers may be precisely estimated. A two-tier survey consists of mutation screening in a sample of patients with specific genetic disorders and in a second sample of newborns from the same population in which mutation frequency is evaluated. We provide the statistical basis for evaluating the incidence of affected and gene carriers in such two-tier mutation-screening surveys, from which the precision of the estimates is derived. Sample-size requirements of such two-tier mutation-screening surveys are evaluated. Considering examples of cystic fibrosis (CF) and medium-chain acyl-CoA dehydrogenase deficiency (MCAD), the two most frequent autosomal recessive disease in Caucasian populations and the two most frequent mutations (delta F508 and G985) that occur on these disease allele-bearing chromosomes, we show that, with 50-100 patients and a 20-fold larger sample of newborns screened for these mutations, the incidence of such diseases and their gene carriers in a population may be quite reliably estimated. The theory developed here is also applicable to rare autosomal dominant diseases for which disease-specific mutations are found.

Mutation and association analysis of the PVR and PVRL2 genes in patients with non-syndromic cleft lip and palate.

Orofacial clefts (OFC; MIM 119530) are among the most common major birth defects. Here, we carried out mutation screening of the PVR and PVRL2 genes, which are both located at an OFC linkage region at 19q13 (OFC3) and are closely related to PVRL1, which has been associated with both syndromic and non-syndromic cleft lip and palate (nsCLP). We screened a total of 73 nsCLP patients and 105 non-cleft controls from the USA for variants in PVR and PVRL2, including all exons and encompassing all isoforms. We identified four variants in PVR and five in PVRL2. One non-synonymous PVR variant, A67T, was more frequent among nsCLP patients than among normal controls, but this difference did not achieve statistical significance.

Mutation analysis of the PVRL1 gene in caucasians with nonsyndromic cleft lip/palate.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P, MIM 119530) is perhaps the most common major birth defect. Homozygous PVRL1 loss-of-function mutations result in an autosomal recessive CL/P syndrome, CLPED1, and a PVRL1 nonsense mutation is associated with sporadic nsCL/P in Northern Venezuela. To address the more general role of PVRL1 variation in risk of nsCL/P, we carried out mutation analysis of PVRL1 in North American and Australian nsCL/P cases and population-matched controls. We identified a total of 15 variants, 5 of which were seen in both populations and 1 of which, an in-frame insertion at Glu442, was more frequent in patients than in controls in both populations, though the difference was not statistically significant. Another variant, which is specific to the PVRL1 beta (HIgR) isoform, S447L, was marginally associated with nsCL/P in North American Caucasian patients, but not in Australian patients, and overall variants that affect the beta-isoform were significantly more frequent among North American patients. One Australian patient had a splice junction mutation of PVRL1. Our results suggest that PVRL1 may play a minor role in susceptibility to the occurrence of nsCL/P in some Caucasian populations, and that variation involving the beta (HIgR) isoform might have particular importance for risk of orofacial clefts. Nevertheless, these results underscore the need for studies that involve very large numbers when assessing the possible role of rare variants in risk of complex traits such as nsCL/P.

Severe aortic and arterial aneurysms associated with a TGFBR2 mutation.

BACKGROUND: A 24-year-old man presented with previously diagnosed Marfan's syndrome. Since the age of 9 years, he had undergone eight cardiovascular procedures to treat rapidly progressive aneurysms, dissection and tortuous vascular disease involving the aortic root and arch, the thoracoabdominal aorta, and brachiocephalic, vertebral, internal thoracic and superior mesenteric arteries. Throughout this extensive series of cardiovascular surgical repairs, he recovered without stroke, paraplegia or renal impairment. INVESTIGATIONS: CT scans, arteriogram, genetic mutation screening of transforming growth factor beta receptors 1 and 2. DIAGNOSIS: Diffuse and rapidly progressing vascular disease in a patient who met the diagnostic criteria for Marfan's syndrome, but was later rediagnosed with Loeys-Dietz syndrome. Genetic testing also revealed a de novo mutation in transforming growth factor beta receptor 2. MANAGEMENT: Regular cardiovascular surveillance for aneurysms and dissections, and aggressive surgical treatment of vascular disease.

DYSREGULATION OF MEOX2 FOLLOWING WT1 MUTATION IN KIDNEY DEVELOPMENT AND WILMS TUMORIGENESIS

Wilms tumor (WT) is a childhood tumor of the kidney and a productive model for understanding the role of genetic alteration and interactions in tumorigenesis. The Wilms tumor gene 1 (WT1) is a transcriptional factor and one of the few genes known to have genetic alterations in WT and has been shown be inactivated in 20% of WTs. However, the mechanisms of how WT1 mutations lead to Wilms tumorigenesis and its influence on downstream genes are unknown. Since it has been established that WT1 is a transcriptional regulator, it has been hypothesized that the loss of WT1 leads to the dysregulation of downstream genes, in turn result in the formation of WTs. To identify the dysregulated downstream genes following WT1 mutations, an Affymetrix GeneChip Human Genome Array was previously conducted to assess the differentially expressed genes in the WT1-wildtype human and WT1-mutant human WTs. Approximately 700 genes were identified as being significantly dysregulated. These genes were further prioritized based on their statistical significance, fold change, chromosomal region, spatial pattern of gene expression and known or putative cellular functions. Mesenchyme homeobox 2 (MEOX2) was one of the most significantly upregulated genes in WT1-mutant WT. MEOX2 is known to play a role in cell proliferation, apoptosis, and differentiation. In addition to its biological roles, it is expressed during early kidney development in the condensed mesenchyme similar to WT1. Furthermore, the use of the Match® web-based tool from the BIOBASE Biological Data base identified a significant predicted WT1 binding site within the first intron of MEOX2. The similarity in spatial gene expression in the developing kidney and the significant predicted WT1 binding site found in the first intron of MEOX2 lead to the development of my hypothesis that MEOX2 is upregulated via a WT1-dependent manner. Here as a part of my master’s work, I have validated the Affymetrix GeneChip Human Genome Array data using an independent set of Wilms tumors. MEOX2 remained upregulated in the mutant WT1 Wilms tumor by 41-fold. Wt1 and Meox2 gene expression were assessed in murine newborn kidney; both Wt1 and Meox2 were expressed in the condensed, undifferentiated metanephric mesenchyme. I have shown that the in vivo ablation of Wt1 during embryonic development at embryonic day (E) 13.5 resulted in the slight increase of Meox2 gene expression by two fold. In order to functionally demonstrate the effect of the loss of Wt1 on Meox2 gene expression in undifferentiated metanephric mesenchyme, I have generated a kidney mesenchymal cell line to genetically ablate Wt1 in vitro by adenoviral infection. The ablation of Wt1 in the kidney mesenchymal cell line resulted in the upregulation of Meox2 by 61-fold. Moreover, the upregulation of Meox2 resulted in the significant induction of p21 and Itgb5. In addition to the dysregulation of these genes the ablation of Wt1 in the kidney mesenchymal cells resulted in decrease in cell growth and loss of cellular adherence. However, it is uncertain whether the upregulation of Meox2 caused this particular cellular phenotype. Overall, I have demonstrated that the upregulation of Meox2 is Wt1-dependent during early kidney development.

Detection and sequence analysis of a nickel-induced mutation in a retroviral model system

We have developed a novel way to assess the mutagenicity of environmentally important metal carcinogens, such as nickel, by creating a positive selection system based upon the conditional expression of a retroviral transforming gene. The target gene is the v-mos gene in MuSVts110, a murine retrovirus possessing a growth temperature dependent defect in expression of the transforming gene due to viral RNA splicing. In normal rat kidney cells infected with MuSVts110 (6m2 cells), splicing of the MuSVts110 RNA to form the mRNA from which the transforming protein, p85$\sp{\rm gag-mos}$, is translated is growth-temperature dependent, occurring at 33 C and below but not at 39 C and above. This splicing "defect" is mediated by cis-acting viral sequences. Nickel chloride treatment of 6m2 cells followed by growth at 39 C, allowed the selection of "revertant" cells which constitutively express p85$\sp{\rm gag-mos}$ due to stable changes in the viral RNA splicing phenotype, suggesting that nickel, a carcinogen whose mutagenicity has not been well established, could induce mutations in mammalian genes. We also show by direct sequencing of PCR-amplified integrated MuSVts110 DNA from a 6m2 nickel-revertant cell line that the nickel-induced mutation affecting the splicing phenotype is a cis-acting 70-base duplication of a region of the viral DNA surrounding the 3$\sp\prime$ splice site. These findings provide the first example of the molecular basis for a nickel-induced DNA lesion and establish the mutagenicity of this potent carcinogen. ^

Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

A coalescent analysis for modeling the mutation process in colorectal cancer

Colorectal cancer is the forth most common diagnosed cancer in the United States. Every year about a hundred forty-seven thousand people will be diagnosed with colorectal cancer and fifty-six thousand people lose their lives due to this disease. Most of the hereditary nonpolyposis colorectal cancer (HNPCC) and 12% of the sporadic colorectal cancer show microsatellite instability. Colorectal cancer is a multistep progressive disease. It starts from a mutation in a normal colorectal cell and grows into a clone of cells that further accumulates mutations and finally develops into a malignant tumor. In terms of molecular evolution, the process of colorectal tumor progression represents the acquisition of sequential mutations. ^ Clinical studies use biomarkers such as microsatellite or single nucleotide polymorphisms (SNPs) to study mutation frequencies in colorectal cancer. Microsatellite data obtained from single genome equivalent PCR or small pool PCR can be used to infer tumor progression. Since tumor progression is similar to population evolution, we used an approach known as coalescent, which is well established in population genetics, to analyze this type of data. Coalescent theory has been known to infer the sample's evolutionary path through the analysis of microsatellite data. ^ The simulation results indicate that the constant population size pattern and the rapid tumor growth pattern have different genetic polymorphic patterns. The simulation results were compared with experimental data collected from HNPCC patients. The preliminary result shows the mutation rate in 6 HNPCC patients range from 0.001 to 0.01. The patients' polymorphic patterns are similar to the constant population size pattern which implies the tumor progression is through multilineage persistence instead of clonal sequential evolution. The results should be further verified using a larger dataset. ^