Functional Analysis of BRCA1 3'UTR Variants Predisposing to Breast Cancer.

Purpose: Breast Cancer (BC) is the main female cancer diagnosed worldwide, and it has been described that few genes, such as BRCA1, have a high penetrance for this type of cancer. In this manuscript, we were interested in evaluating the effect of 3'UTR variants on BRCA1 expression. Patients and Methods: To accomplish this objective, Whole Exome Sequencing (WES) data of 400 patients with unselected BC was used to filter variants located in the region of interest of BRCA1 gene, finding two of them (c.*36C>G and c.*369_373del). miRGate and miRanda in silico tools were used to predict microRNA (miRNA) interaction. Results: The two variants (c.*36C>G, c.*369_373del) were predicted to affect miRNA interaction. After cloning of BRCA1 3'UTR into pMIR-Report vector, the construct was transfected into two BC cell lines (MDA-MB-231 and MCF-7), and the variant c.*36C>G evidenced overexpression of reporter gene luciferase, showing that the transcript was not being degraded by the miRNA in MDA-MB-231 cells. Conclusion: The variant seems to protect against Triple Negative BC probably due to the expression level of miRNA in this particular cell line (MDA-MB-231). This is consistent with the clinical history of the patients who harbor BC Hormone Receptors positive (HR+).

Next-generation sequencing of host genetics risk factors associated with COVID-19 severity and long-COVID in Colombian population.

Coronavirus disease 2019 (COVID-19) was considered a major public health burden worldwide. Multiple studies have shown that susceptibility to severe infections and the development of long-term symptoms is significantly influenced by viral and host factors. These findings have highlighted the potential of host genetic markers to identify high-risk individuals and develop target interventions to reduce morbimortality. Despite its importance, genetic host factors remain largely understudied in Latin-American populations. Using a case-control design and a custom next-generation sequencing (NGS) panel encompassing 81 genetic variants and 74 genes previously associated with COVID-19 severity and long-COVID, we analyzed 56 individuals with asymptomatic or mild COVID-19 and 56 severe and critical cases. In agreement with previous studies, our results support the association between several clinical variables, including male sex, obesity and common symptoms like cough and dyspnea, and severe COVID-19. Remarkably, thirteen genetic variants showed an association with COVID-19 severity. Among these variants, rs11385942 (p < 0.01; OR = 10.88; 95% CI = 1.36-86.51) located in the LZTFL1 gene, and rs35775079 (p = 0.02; OR = 8.53; 95% CI = 1.05-69.45) located in CCR3 showed the strongest associations. Various respiratory and systemic symptoms, along with the rs8178521 variant (p < 0.01; OR = 2.51; 95% CI = 1.27-4.94) in the IL10RB gene, were significantly associated with the presence of long-COVID. The results of the predictive model comparison showed that the mixed model, which incorporates genetic and non-genetic variables, outperforms clinical and genetic models. To our knowledge, this is the first study in Colombia and Latin-America proposing a predictive model for COVID-19 severity and long-COVID based on genomic analysis. Our study highlights the usefulness of genomic approaches to studying host genetic risk factors in specific populations. The methodology used allowed us to validate several genetic variants previously associated with COVID-19 severity and long-COVID. Finally, the integrated model illustrates the importance of considering genetic factors in precision medicine of infectious diseases.

Linking genotype to trophoblast phenotype in preeclampsia and HELLP syndrome associated with STOX1 genetic variants.

Preeclampsia is a major hypertensive pregnancy disorder with a 50% heritability. The first identified gene involved in the disease is STOX1, a transcription factor, whose variant Y153H predisposes to the disease. Two rare mutations were also identified in Colombian women affected by the hemolysis, elevated liver enzyme, low platelet syndrome, a complication of preeclampsia (T188N and R364X). Here, we explore the effects of these variants in trophoblast cell models (BeWo) where STOX1 was previously invalidated. We firstly showed that STOX1 knockout alters response to oxidative stress, cell proliferation, and fusion capacity. Then, we showed that mutant versions of STOX1 trigger alterations in gene profiles, growth, fusion, and oxidative stress management. The results also reveal alterations of the STOX interaction with DNA when the mutations affected the DNA-binding domain of STOX1 (Y153H and T188N). We also reveal here that a major contributor of these effects appears to be the E2F3 transcription factor.

Functional analysis of CTLA4 promoter variant and its possible implication in colorectal cancer immunotherapy.

Background: Colorectal cancer (CRC) is a prevalent cancer, ranking as the third most common. Recent advances in our understanding of the molecular causes of this disease have highlighted the crucial role of tumor immune evasion in its initiation and progression. CTLA4, a receptor that acts as a negative regulator of T cell responses, plays a pivotal role in this process, and genetic variations in CTLA4 have been linked to CRC susceptibility, prognosis, and response to therapy. Methods: We conducted a case-control study involving 98 CRC patients and 424 controls. We genotyped the CTLA4 c.-319C > T variant (rs5742909) and performed an association analysis by comparing allele frequencies between the patients and controls. To assess the potential functional impact of this variant, we first performed an In Silico analysis of transcription factor binding sites using Genomatix. Finally, to validate our findings, we conducted a luciferase reporter gene assay using different cell lines and an electrophoretic mobility shift assay (EMSA). Results: The case-control association analysis revealed a significant association between CTLA4 c.-319C > T and CRC susceptibility (p = 0.023; OR 1.89; 95% CI = 1.11-3.23). Genomatix analysis identified LEF1 and TCF7 transcription factors as specific binders to CTLA4 c.-319C. The reporter gene assay demonstrated notable differences in luciferase activity between the c.-319 C and T alleles in COS-7, HCT116, and Jurkat cell lines. EMSA analysis showed differences in TCF7 interaction with the CTLA4 C and T alleles. Conclusion: CTLA4 c.-319C > T is associated with CRC susceptibility. Based on our functional validation results, we proposed that CTLA4 c.-319C > T alters gene expression at the transcriptional level, triggering a stronger negative regulation of T-cells and immune tumoral evasion.

Pharmacogenomic profile of actionable molecular variants related to drugs commonly used in anesthesia: WES analysis reveals new mutations.

Background: Genetic interindividual variability is associated with adverse drug reactions (ADRs) and affects the response to common drugs used in anesthesia. Despite their importance, these variants remain largely underexplored in Latin-American countries. This study describes rare and common variants found in genes related to metabolism of analgesic and anaesthetic drug in the Colombian population. Methods: We conducted a study that included 625 Colombian healthy individuals. We generated a subset of 14 genes implicated in metabolic pathways of common medications used in anesthesia and assessed them by whole-exome sequencing (WES). Variants were filtered using two pipelines: A) novel or rare (minor allele frequency-MAF <1%) variants including missense, loss-of-function (LoF, e.g., frameshift, nonsense), and splice site variants with potential deleterious effect and B) clinically validated variants described in the PharmGKB (categories 1, 2 and 3) and/or ClinVar databases. For rare and novel missense variants, we applied an optimized prediction framework (OPF) to assess the functional impact of pharmacogenetic variants. Allelic, genotypic frequencies and Hardy-Weinberg equilibrium were calculated. We compare our allelic frequencies with these from populations described in the gnomAD database. Results: Our study identified 148 molecular variants potentially related to variability in the therapeutic response to 14 drugs commonly used in anesthesiology. 83.1% of them correspond to rare and novel missense variants classified as pathogenic according to the pharmacogenetic optimized prediction framework, 5.4% were loss-of-function (LoF), 2.7% led to potential splicing alterations and 8.8% were assigned as actionable or informative pharmacogenetic variants. Novel variants were confirmed by Sanger sequencing. Allelic frequency comparison showed that the Colombian population has a unique pharmacogenomic profile for anesthesia drugs with some allele frequencies different from other populations. Conclusion: Our results demonstrated high allelic heterogeneity among the analyzed sampled, enriched by rare (91.2%) variants in pharmacogenes related to common drugs used in anesthesia. The clinical implications of these results highlight the importance of implementation of next-generation sequencing data into pharmacogenomic approaches and personalized medicine.

Performance Comparison of a Duplex Implementation of the CDC EUA 2019-nCoV Assay with the Seegene Allplex-SARS-CoV-2 Assay for the Detection of SARS-CoV-2 in Nasopharyngeal Swab Samples.

RT-PCR tests have become the gold standard for detecting the SARS-CoV-2 virus in the context of the COVID-19 pandemic. Because of the extreme number of cases in periodic waves of infection, there is a severe financial and logistical strain on diagnostic laboratories. For this reason, alternative implementations and validations of academic protocols that employ the lowest cost and the most widely available equipment and reagents found in different regions are essential. In this study, we report an alternative implementation of the EUA 2019-nCoV CDC assay which uses a previously characterized duplex PCR reaction for the N1 and RNAse P target regions and an additional uniplex reaction for the N2 target region. Taking advantage of the Abbott m2000 Sample Preparation System and NEB Luna Universal Probe One-Step RT-qPCR kit, some of the most widely available and inexpensive nucleic acid extraction and amplification platforms, this modified test shows state-of-the-art analytical and clinical sensitivities and specificities when compared with the Seegene Allplex-SARS-CoV-2 assay. This implementation has the potential to be verified and implemented by diagnostic laboratories around the world to guarantee low-cost RT-PCR tests that can take advantage of widely available equipment and reagents.

Characterization of ADME Gene Variation in Colombian Population by Exome Sequencing.

In genes related to drug pharmacokinetics, molecular variations determine interindividual variability in the therapeutic efficacy and adverse drug reactions. The assessment of single-nucleotide variants (SNVs) is used with growing frequency in pharmacogenetic practice, and recently, high-throughput genomic analyses obtained through next-generation sequencing (NGS) have been recognized as powerful tools to identify common, rare and novel variants. These genetic profiles remain underexplored in Latin-American populations, including Colombia. In this study, we investigated the variability of 35 genes included in the ADME core panel (absorption, distribution, metabolism, and excretion) by whole-exome sequencing (WES) of 509 unrelated Colombian individuals with no previous reports of adverse drug reactions. Rare variants were filtered according to the minor allele frequencies (MAF) <1% and potential deleterious consequences. The functional impact of novel and rare missense variants was assessed using an optimized framework for pharmacogenetic variants. Bioinformatic analyses included the identification of clinically validated variants described in PharmGKB and ClinVar databases. Ancestry from WES data was inferred using the R package EthSEQ v2.1.4. Allelic frequencies were compared to other populations reported in the public gnomAD database. Our analysis revealed that rare missense pharmacogenetic variants were 2.1 times more frequent than common variants with 121 variants predicted as potentially deleterious. Rare loss of function (LoF) variants were identified in 65.7% of evaluated genes. Regarding variants with clinical pharmacogenetic effect, our study revealed 89 sequence variations in 28 genes represented by missense (62%), synonymous (22.5%), splice site (11.2%), and indels (3.4%). In this group, ABCB1, ABCC2, CY2B6, CYP2D6, DPYD, NAT2, SLC22A1, and UGTB2B7, are the most polymorphic genes. NAT2, CYP2B6 and DPYD metabolizer phenotypes demonstrated the highest variability. Ancestry analysis indicated admixture in 73% of the population. Allelic frequencies exhibit significant differences with other Latin-American populations, highlighting the importance of pharmacogenomic studies in populations of different ethnicities. Altogether, our data revealed that rare variants are an important source of variability in pharmacogenes involved in the pharmacokinetics of drugs and likely account for the unexplained interindividual variability in drug response. These findings provide evidence of the utility of WES for pharmacogenomic testing and into clinical practice.

Association Between the LZTFL1 rs11385942 Polymorphism and COVID-19 Severity in Colombian Population.

Genetic and non-genetic factors are responsible for the high interindividual variability in the response to SARS-CoV-2. Although numerous genetic polymorphisms have been identified as risk factors for severe COVID-19, these remain understudied in Latin-American populations. This study evaluated the association of non-genetic factors and three polymorphisms: ACE rs4646994, ACE2 rs2285666, and LZTFL1 rs11385942, with COVID severity and long-term symptoms by using a case-control design. The control group was composed of asymptomatic/mild cases (n = 61) recruited from a private laboratory, while the case group was composed of severe/critical patients (n = 63) hospitalized in the Hospital Universitario Mayor-Méderi, both institutions located in Bogotá, Colombia. Clinical follow up and exhaustive revision of medical records allowed us to assess non-genetic factors. Genotypification of the polymorphism of interest was performed by amplicon size analysis and Sanger sequencing. In agreement with previous reports, we found a statistically significant association between age, male sex, and comorbidities, such as hypertension and type 2 diabetes mellitus (T2DM), and worst outcomes. We identified the polymorphism LZTFL1 rs11385942 as an important risk factor for hospitalization (p < 0.01; OR = 5.73; 95% CI = 1.2-26.5, under the allelic test). Furthermore, long-term symptoms were common among the studied population and associated with disease severity. No association between the polymorphisms examined and long-term symptoms was found. Comparison of allelic frequencies with other populations revealed significant differences for the three polymorphisms investigated. Finally, we used the statistically significant genetic and non-genetic variables to develop a predictive logistic regression model, which was implemented in a Shiny web application. Model discrimination was assessed using the area under the receiver operating characteristic curve (AUC = 0.86; 95% confidence interval 0.79-0.93). These results suggest that LZTFL1 rs11385942 may be a potential biomarker for COVID-19 severity in addition to conventional non-genetic risk factors. A better understanding of the impact of these genetic risk factors may be useful to prioritize high-risk individuals and decrease the morbimortality caused by SARS-CoV2 and future pandemics.

Preeclampsia, Natural History, Genes, and miRNAs Associated with the Syndrome.

Preeclampsia (PE) is a hypertensive disease that affects pregnant women after 20 weeks of gestation. This disease is associated with an important risk of maternal and fetal mortality. PE is described as a placental pathology because, after delivery, most women recover normal arterial pressure. Poor invasion of the spiral arteries is a phenomenon well described in PE; this leads to a hypoxic uterine bed and imbalance of antiangiogenic and proangiogenic factors in the uteroplacental region, which in turn triggers the disease phenotype. The causes of the pathology are unclear; nevertheless, numerous approaches, including next-generation sequencing, association, and case control and miRNA studies, have shed light on the genetic/molecular basis of PE. These studies help us better understand the disease to advance new treatment strategies.

Whole-Exome Sequencing in Patients Affected by Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis Reveals New Variants Potentially Contributing to the Phenotype.

BACKGROUND: Adverse drug reactions (ADRs) are frequent occurring events that can essentially be defined as harmful or unpleasant symptoms secondary to the use of a medicinal product. ADRs involve a wide spectrum of clinical manifestations ranging from minor itching and rash to life-threatening reactions. Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are rare ADRs. SJS-TEN may be considered a polygenic pathology due to additive/epistatic effects caused by sequence variants in numerous genes. Next-generation sequencing (NGS) represents a potentially interesting exploration tool in such scenario as it facilitates the simultaneous analysis of large genomic regions and genes at affordable cost. METHODS: The present study has involved using whole-exome sequencing (WES) for the first time on SJS-TEN patients. It involved robust and innovative multistep bioinformatics analysis focusing on 313 candidate genes potentially participating in the disease's aetiology, specific drugs' metabolism and gene regulation. RESULTS: We identified combinations of frequently occurring and rare variants that may contribute to the disease's pathogenesis. Depending on the specific drug being taken, different variants (and alleles) in NAT2, CYP2D8, CYP2B6, ABCC2, UGT2B7 and TCF3 were identified as coherent candidates representing potential future markers for SJS-TEN. CONCLUSION: The present study proposed and has described (for the first time) a large-scale genomic analysis of patients affected by SJS-TEN. The genes and variants identified represent relevant candidates potentially participating in the disease's pathogenesis. Corroborating that proposed by others, we found that complex combinations of frequently occurring and rare variants participating in particular drug metabolism molecular cascades could be associated with the phenotype. TCF3 TF may be considered a coherent candidate for SJS-TEN that should be analysed in new cohorts of patients having ADRs.

A Novel Splice-Site Mutation in the ELN Gene Suggests an Alternative Mechanism for Vascular Elastinopathies.

The ELN gene encodes elastin, a fundamental protein of the extracellular matrix that confers elasticity to different tissues including blood vessels. The formation of elastin fibers is a complex process involving monomer coacervation and subsequent crosslinking. Mutations in exons 1-29 of the ELN gene have been linked to supravalvular aortic stenosis (SVAS) whereas mutations in exons 30-33 are associated with autosomal dominant cutis laxa (ADCL). This striking segregation has led to the hypothesis that distinct molecular mechanisms underlie both diseases. SVAS is believed to arise through haploinsufficiency while ADCL is hypothesized to be caused by a dominant negative effect. Here, we describe a patient with SVAS harboring a novel splice-site mutation in the last exon of ELN. The location of this mutation is not consistent with current knowledge of SVAS, since all mutations reported in the C-terminus have been found in ADCL patients, and a thorough evaluation did not reveal significant skin involvement in this case. RT-PCR analysis of skin tissue showed that C-terminal mutations in the region can lead to the production of aberrant transcripts through intron retention and activation of cryptic splice sites and suggest that disruption of the very last exon can lead to functional haploinsufficiency potentially related to SVAS.

Identifying new potential genetic biomarkers for HELLP syndrome using massive parallel sequencing.

BACKGROUND: Preeclampsia (PE) is a frequently occurring multisystemic disease affecting ~5% of pregnancies. PE patients may develop HELLP syndrome (haemolysis, elevated liver enzymes, and low platelet), a mother and foetus life-threatening condition. Research into HELLP's genetic origin has been relatively unsuccessful, mainly because normal placental function and blood pressure regulation involve the fine-regulation of hundreds of genes. OBJECTIVE: To identify new genes and mutations constituting potential biomarkers for HELLP syndrome. STUDY DESIGN: The present case-control study involved whole-exome sequencing of 79 unrelated HELLP women. Candidate variants were screened in a control population constituted by 176 individuals. Stringent bioinformatics filters were used for selecting potentially etiological sequence variants in a subset of 487 genes. We used robust in silico mutation modelling for predicting the potential effect on protein structure. RESULTS: We identified numerous sequence variants in genes related to angiogenesis/coagulation/blood pressure regulation, cell differentiation/communication/adhesion, cell cycle and transcriptional gene regulation, extracellular matrix biology, lipid metabolism and immunological response. Five sequence variants generated premature stop codons in genes playing an essential role in placental physiology (STOX1, PDGFD, IGF2, MMP1 and DNAH11). Six variants (ERAP1- p.Ile915Thr, ERAP2- p.Leu837Ser, COMT-p.His192Gln, CSAD-p.Pro418Ser, CDH1- p.Ala298Thr and CCR2-p.Met249Lys) led to destabilisation of protein structure as they had significant energy and residue interaction-related changes. We identified at least two mutations in 57% of patients, arguing in favour of a polygenic origin for the HELLP syndrome. CONCLUSION: Our results provide novel evidence regarding PE/HELLP's genetic origin, leading to new biomarkers, having potential clinical usefulness, being proposed.

DNA demethylation agent 5azadC downregulates HPV16 E6 expression in cervical cancer cell lines independently of TBX2 expression.

HPV16 is the most carcinogenic human papillomavirus and causes >50% of cervical cancers, the majority of anal cancers and 30% of oropharyngeal squamous cell carcinomas. HPV carcinogenesis relies on the continuous expression of the two main viral oncoproteins E6 and E7 that target >150 cellular proteins. Among them, epigenetic modifiers, including DNA Methyl Transferases (DNMT), are dysregulated, promoting an aberrant methylation pattern in HPV-positive cancer cells. It has been previously reported that the treatment of HPV-positive cervical cancer cells with DNMT inhibitor 5-aza-2'-deoxycytidine (5azadC) caused the downregulation of E6 expression due to mRNA destabilization that was mediated by miR-375. Recently, the T-box transcription factor 2 (TBX2) has been demonstrated to repress HPV LCR activity. In the current study, the role of TBX2 in E6 repression was investigated in HPV16 cervical cancer cell lines following 5azadC treatment. A decrease of E6 expression was accompanied by p53 and p21 restoration. While TBX2 mRNA was upregulated in 5azadC-treated SiHa and Ca Ski cells, TBX2 protein was not detectable. Furthermore, the overexpression of TBX2 protein in cervical cancer cells did not allow the repression of E6 expression. The TBX2 transcription factor is therefore unlikely to be associated with the repression of E6 following 5azadC treatment of SiHa and Ca Ski cells.

Structural and functional responses of invertebrate communities to climate change and flow regulation in alpine catchments.

Understanding and predicting how biological communities respond to climate change is critical for assessing biodiversity vulnerability and guiding conservation efforts. Glacier- and snow-fed rivers are one of the most sensitive ecosystems to climate change, and can provide early warning of wider-scale changes. These rivers are frequently used for hydropower production but there is minimal understanding of how biological communities are influenced by climate change in a context of flow regulation. This study sheds light on this issue by disentangling structural (water temperature preference, taxonomic composition, alpha, beta and gamma diversities) and functional (functional traits, diversity, richness, evenness, dispersion and redundancy) effects of climate change in interaction with flow regulation in the Alps. For this, we compared environmental and aquatic invertebrate data collected in the 1970s and 2010s in regulated and unregulated alpine catchments. We hypothesized a replacement of cold-adapted species by warming-tolerant ones, high temporal and spatial turnover in taxa and trait composition, along with reduced taxonomic and functional diversities in consequence of climate change. We expected communities in regulated rivers to respond more drastically due to additive or synergistic effects between flow regulation and climate change. We found divergent structural but convergent functional responses between free-flowing and regulated catchments. Although cold-adapted taxa decreased in both of them, greater colonization and spread of thermophilic species was found in the free-flowing one, resulting in higher spatial and temporal turnover. Since the 1970s, taxonomic diversity increased in the free flowing but decreased in the regulated catchment due to biotic homogenization. Colonization by taxa with new functional strategies (i.e. multivoltine taxa with small body size, resistance forms, aerial dispersion and reproduction by clutches) increased functional diversity but decreased functional redundancy through time. These functional changes could jeopardize the ability of aquatic communities facing intensification of ongoing climate change or new anthropogenic disturbances.

Functional evidence implicating NOTCH2 missense mutations in primary ovarian insufficiency etiology.

Primary ovarian insufficiency (POI) is a frequently occurring disease affecting women under 40 years old. Recently, we have analyzed unrelated POI women via whole exome sequencing (WES) and identified NOTCH2 mutations underlying possible functional effects. The present study involved reanalyzing of WES assays. We used in the KGN granulosa-like cell model, a synthetic gene reporter construct driving luciferase gene expression to assess the functional effects of five NOTCH2 mutations identified in POI patients. We found that NOTCH2-p.Ser1804Leu, p.Ala2316Val, and p.Pro2359Ala mutations had a functional impact on the protein's transcriptional activity. The results have demonstrated for the first time that NOTCH2 mutations contribute to POI etiology. We therefore recommend sequencing NOTCH2's open reading frame in large panels of POI patients to establish an accurate genotype-phenotype correlation. We cannot rule out the fact that patients affected by Alagille syndrome carrying NOTCH2 mutations may suffer ovarian dysfunction.

5azadC treatment upregulates miR-375 level and represses HPV16 E6 expression.

High-risk human papillomaviruses are the etiological agents of cervical cancer and HPV16 is the most oncogenic genotype. Immortalization and transformation of infected cells requires the overexpression of the two viral oncoproteins E6 and E7 following HPV DNA integration into the host cell genome. Integration often leads to the loss of the E2 open reading frame and the corresponding protein can no longer act as a transcriptional repressor on p97 promoter. Recently, it has been proposed that long control region methylation also contributes to the regulation of E6/E7 expression.To determine which epigenetic mechanism is involved in HPV16 early gene regulation, 5-aza-2'-deoxycytidine was used to demethylate Ca Ski and SiHa cell DNA. Decreased expression of E6 mRNA and protein levels was observed in both cell lines in an E2-independent manner. E6 repression was accompanied by neither a modification of the main cellular transcription factor expression involved in long control region regulation, nor by a modification of the E6 mRNA splicing pattern. In contrast, a pronounced upregulation of miR-375, known to destabilize HPV16 early viral mRNA, was observed. Finally, the use of miR-375 inhibitor definitively proved the involvement of miR-375 in E6 repression. These results highlight that cellular DNA methylation modulates HPV16 early gene expression and support a role for epigenetic events in high-risk HPV associated-carcinogenesis.

Long-term storage and safe retrieval of human papillomavirus DNA using FTA elute cards.

Biobanking or collection and storage of specimens for future research purposes have become an essential tool in many fields of biomedical research and aims to provide a better understanding of disease mechanisms as well as the identification of disease-specific biomarkers that can navigate in complex diseases. In this study, we assessed the use of Flinders Technology Associates (FTA) cards as a long-term storage device for cervical specimens with suspected human papillomavirus (HPV) infections. HPV detection and genotyping results in liquid-based transport media were compared to HPV results from FTA cards. The overall agreement for the presence of any HPV infection between liquid-based medium and FTA cards stored for 1 year at ambient temperature was 100%. Reproducibility analysis of HPV detection and genotyping from FTA cards demonstrated that FTA cards are a reliable medium to store and preserve viral nucleic acids. Biobanking of cervical cells on FTA cards may provide a key resource for epidemiological and retrospective HPV studies.

Methylation of human papillomavirus Type 16 CpG sites at E2-binding site 1 (E2BS1), E2BS2, and the Sp1-binding site in cervical cancer samples as determined by high-resolution melting analysis-PCR.

High-risk (HR) human papillomavirus (HPV)-associated carcinogenesis is driven mainly by the overexpression of E7 and E6 oncoproteins following viral DNA integration and the concomitant loss of the E2 open reading frame (ORF). However, the integration of HR-HPV DNA is not systematically observed in cervical cancers. The E2 protein acts as a transcription factor that governs viral oncogene expression. The methylation of CpGs in the E2-binding sites (E2BSs) in the viral long control region abrogates E2 binding, thus impairing the E2-mediated regulation of E7/E6 transcription. Here, high-resolution melting (HRM)-PCR was developed to quantitatively analyze the methylation statuses of E2BS1, E2BS2, and the specificity protein 1 (Sp1)-binding site in 119 HPV16-positive cervical smears. This is a rapid assay that is suitable for the analysis of cervical samples. The proportion of cancer samples with methylated E2BS1, E2BS2, and Sp1-binding site CpGs was 47%, whereas the vast majority of samples diagnosed as being within normal limits, low-grade squamous intraepithelial lesions (LSIL), or high-grade squamous intraepithelial lesions (HSIL) harbored unmethylated CpGs. Methylation levels varied widely, since some cancer samples harbored up to 60% of methylated HPV16 genomes. A pyrosequencing approach was used as a confirmation test and highlighted that quantitative measurement of methylation can be achieved by HRM-PCR. Its prognostic value deserves to be investigated alone or in association with other biomarkers. The reliability of this single-tube assay offers great opportunities for the investigation of HPV16 methylation in other HPV-related cancers, such as head and neck cancers, which are a major public health burden.