INTRODUCTION
One of the most prevalent neoplastic pathologies associated with male sex is prostate cancer. The estimated prevalence is 1.1 million people worldwide1-3 and it is impacted by ethnicity and geographic location.4 Populations of African descent are the most affected, showing an 11% increase in prevalence in recent years.5 The Southwest region of Colombia is inhabited by populations of Latin American and African descent in approximately the same proportion but with different rates of disease incidence.6
Variants of certain genes have been associated with a higher frequency of prostate cancer (BRCA1-2, ATM, NBN, TMPRSS2, among others).7 Serine proteases, such as the TMPRSS2 gene, are recognized through their mechanisms of action in inflammatory and immune processes. That gene is located on chromosome 21q22.3 and is expressed at the apex of the secretory epithelium of the glands. Fusion with members of the ETS family is the most frequent chromosomal re-arrangement found in 50% of prostate cancers, mainly produced by the microdeletion of a portion of the TMPRSS2 gene.8 The TMPRSS2 gene and the fusion gene (TMPRSS2:ERG) have been associated with the severity and prognosis of prostate cancer, although the actual pathophysiologic process or the variant associated with that condition are not very well known.9 The fusion gene has been widely studied and at present has been postulated as one of the most important biomarkers for diagnostic and prognostic purposes in the prostate cancer population.10 There are reports in the literature on the single nucleotide polymorphisms (SNPs) most frequently related to those clinical scenarios.
The present study is important because there are no similar descriptive studies characterizing the presence of the TMPRSS2 gene and its variants in a population from Southwestern Colombia. Our study focuses on describing the frequency of the allelic variants of the TMPRSS2 gene in that population.
MATERIALS AND METHODS
A descriptive, observational study was conducted on persons, regardless of age, from Southwestern Colombia (Nariño, Cauca, Putumayo, and Valle), within the time frame of 2014 to 2016.
Sample size
According to the expected frequency for hereditary prostate cancer (≈ 15%), alpha 5%, and an expected error of 5%, the calculated sample size was 162 people and convenience sampling was carried out.
Complete exome sequencing was performed, which enabled the sequencing of all protein-coding regions (exome) in the genome, thus identifying the variants that could alter the sequence of a protein. It was carried out as follows:
DNA extraction
Blood was drawn from each patient for DNA extraction. All drops of blood were collected and dried on filter paper. The filter paper was then immersed in a phosphate buffer utilizing the DNeasy kit from the QIAGEN® company (Hilden, Germany-Operational). Each extraction was quantified, and its quality was verified, to continue the sequencing processing.
Sequencing protocol
DNA aliquots from each sample underwent a preparation process with the TruSeq Exome Library Prep®. The resulting libraries were then normalized for sequencing using the TruSeq Rapid Exome®. The kits were provided by Illumina® from San Diego, California, USA. The normalized fragments with their corresponding adaptors for sequencing were charged in a HiSeq2500 machine.
We sequenced the full exome and identified the related variants, specifically the SNPs for the TMPRSS2 gene that is associated with prostate cancer (PCa).
The present project was conducted following all ethical international standards. Descriptive statistics were performed in R and the results are shown in frequency tables for each gene and its associated variants. Finally, we looked for the variants in the following public databases: Exome Aggregation Consortium (ExAC), PharmGKB,11 Clinvar,12 Ensemble, and dbSNP,13 searching for a pattern through which we could use the variants we found as markers.
RESULTS
One hundred sixty-two patients were included in the study, providing 7,315,466 sequence data, and the TMPRSS2 gene was found in 414 data (4.3%). Missense variants were identified in 23% of the data, although the most frequent variants were synonymous variants and introns. Only one stop variant was found in those data (Table 1).
Variant | Absolute Frequency | Percentage (%) |
---|---|---|
5UTR | 3 | 0.72 |
Intron | 144 | 34.78 |
Missense | 98 | 23.67 |
Stop | 1 | 0.24 |
Synonymous | 168 | 40.58 |
In addition, the most common variants for the TMPRSS2 gene were: rs140530035 (32.12%), rs17854725 (19.8%), and rs2298659 (13.5%) (Table 2).
Variants | Absolute Frequency | Percentage (%) |
No identifier available | 36 | 8.70 |
rs12329760 | 44 | 10.63 |
rs140530035 | 133 | 32.13 |
rs143049780 | 1 | 0.24 |
rs148125094 | 1 | 0.24 |
rs149527323 | 1 | 0.24 |
rs17854725 | 82 | 19.81 |
rs181414852 | 1 | 0.24 |
rs2298659 | 56 | 13.53 |
rs3787950 | 15 | 3.62 |
rs61735789 | 2 | 0.48 |
rs61735790 | 1 | 0.24 |
rs61735792 | 1 | 0.24 |
rs61735793 | 1 | 0.24 |
rs61735794 | 3 | 0.72 |
rs61735795 | 1 | 0.24 |
rs75603675 | 35 | 8.45 |
DISCUSSION
Transmembrane protease serine 2, also called TMPRSS2, is a protease composed of 492 amino acids expressed on the cell surface of multiple organs and they are theorized to be strategically located to regulate cell-cell interactions. The TMPRSS2 gene has been shown to be positively regulated by androgenic hormones in neoplastic tissue, possibly modulating the inflammatory response of prostate cells through the activation of PAR-2.14,15
Prostate cancer is one of the most frequent cancers in males and the TMPRSS2 gene has historically been associated with that malignant tumor. Numerous authors have conducted studies over the past decades in an attempt to link the presence of the TMPRSS2 gene with the frequency of cancer and its prognosis.16 Although there are studies that have found that the TMPRSS2 gene does not represent a worse prognosis for prostate cancer,17 an important fusion of that gene with the ERG gene was described, with an increasing relation to the diagnosis and aggressiveness of prostate cancer (present in 50% of high-risk prostate cancers).18,19
We found a low frequency of the allelic variant associated with the TMPRSS2:ERG fusion gene in our cancer-free population from Southwestern Colombia. The rs12329760 variant, albeit not the most frequent SNP found in the present study, is reported in the literature to have a non-negligible allele frequency (AF) in populations from East Asia and Northern Europe (0.38 and 0.37, respectively), with a major homozygote ratio (> 7%). Frequency in the Hispanic population is 0.155, with a low number of homozygotes.20 It should be noted that Southwestern Colombia has a large population of African descent, in which a higher frequency of said allelic variant (0.29) has been identified. That is an important fact to keep in mind when identifying new biomarkers for prostate cancer.20
The most sequenced polymorphism in the present study was rs140530035. It is a very common intron in the world population and the allele frequency of that variant reaches 0.9.21 Comparing populations, inhabitants of northern Europe (Finland) have an AF of 0.99, whereas it is only 0.66 in the so-called Latino population, according to Lek et al. 2016.21
The rs17854725 and rs2298659 polymorphisms are synonymous variants that are rare in the Latin American population, according to the literature, with an AF of 0.15 or less, and they have no known pathologic associations. Likewise, the rs75603675 polymorphism is not known to be associated with any pathology.21
Strengths and limitations
The present study is the first to describe the relation of the TMPRSS2 gene and its allelic variants to a Southwestern Colombian cancer-free population. An advantage of the project was the quality of the study’s samples, analyses, and data. Several variants associated with the TMPRSS2 gene were identified. That is very important information for the performance of future longitudinal studies in cancer-free patients to determine the risk for that disease.
A limitation of the present study was the fact that we did not find any information associated with the presence of a pathologic relationship to prostate cancer.